C#实现Gauss 消去法

一、目的

1、掌握高斯消去法基础原理
2、掌握高斯消去法解方程组的步骤
3、能用程序语言对Gauss消去法进行编程实现

二、代码

1、Gauss.cs

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace 高斯消元
{
    class Gauss
    {
        int n = 0;
        public int N
        {
            get { return n; }
            set { n = value; }
        }
        double[,] a = null;
        public double[,] A 
        { get { return a; } set { a = value; } }
        double[] x = null; 
        public double[] X { get { return x; } set { x = value; } }
        
        //输入
    public void Input()
    {
        Console.WriteLine("请输入方程的阶数");
        n = int.Parse(Console.ReadLine());
        a = new double[n + 1, n + 2];
        //  x=new double[n+1];
        Console.WriteLine("请输入方程的系数");
        for (int i = 1; i <= n; i++)
        {
            string r = Console.ReadLine();
            string[] rs = r.Split(' ');
            for (int j = 1; j <= n + 1; j++)
            {
                a[i, j] = double.Parse(rs[j - 1]);
            }
        }
    }
        //输出A
        public void OutputA()
        {
            Console.WriteLine("A[i,j]如下");
            for (int i = 1; i <= n; i++)
            {
                for (int j = 1; j <= n + 1; j++)
                {
                    Console.Write("{0,8:f2}", a[i, j]);
                }
                Console.WriteLine();
            }
        }

        //消元
        public void CalcuA()
        {
            for (int k = 1; k <= n - 1; k++)
            {
                for (int i = k + 1; i <= n; i++)
                {
                    double lik = a[i, k] / a[k, k];
                    for (int j = k + 1; j <= n + 1; j++)
                    {
                        a[i, j] = a[i, j] - lik * a[k, j];
                    }
                    a[i, k] = 0;
                }
            }
            Console.WriteLine("消元过程");
        }

       public void CalcuX( )
        {
            x = new double[n + 1];
            //正在计算回代过程
            for (int i = n; i >= 1; i--)
            {
                double sum = 0;
                for (int j = i + 1; j <= n; j++)
                {
                    sum += a[i, j] * x[j];
                }
                x[i] = (a[i, n + 1] - sum) / a[i, i];
            }

        }


        public void OutputX( )
        {
            int n = x.Length - 1;
            Console.WriteLine("方程解如下");
            for (int i = 1; i <= n; i++)
            {
                Console.WriteLine("X[{0}]={1,8:f2}", i, x[i]);
            }
        }

        public void Calcu()
        {
          Input();
          OutputA();
            //消元
           CalcuA();
           OutputA();
           CalcuX();
           OutputX();
            Console.Read();
        }
    }

}

2、Program.cs

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;


namespace 高斯消元
{
    class Program
    {
        //输入
        static void Input(ref int n,ref double[,] a)
        {
            Console.WriteLine("请输入方程的阶数");
            n = int.Parse(Console.ReadLine());
            a = new double[n + 1, n + 2];
            //  x=new double[n+1];
            Console.WriteLine("请输入方程的系数");
            for (int i=1;i<=n;i++)
            {
                string r = Console.ReadLine();
                string[] rs = r.Split(' ');
                for (int j=1; j<=n+1; j++)
                {
                    a[i, j] = double.Parse(rs[j - 1]);
                }
            }
        }

        //输出A
        static void OutputA( int n,  double[,] a)
        {
            Console.WriteLine("A[i,j]如下");
            for (int i=1;i<=n;i++)
            {
            for (int j = 1; j <= n+1; j++) 
            {
                Console.Write("{0,8:f2}", a[i, j]);
            }
                Console.WriteLine();
            }
        }
   

          
        //消元
        static void CalcuA(int n,ref double[,]a)
        {
            for(int k=1;k<=n-1;k++)
            {
                for (int i=k+1;i<=n;i++)
                {
                    double lik = a[i, k] / a[k, k];
                    for (int j=k+1;j<=n+1;j++)
                    {
                        a[i, j] = a[i, j] - lik * a[k, j];
                    }
                    a[i, k] = 0;
                }
            }
            Console.WriteLine("消元过程");
        }

        static void CalcuX(int n, double[,] a, ref double[] x)
        {
            x = new double[n + 1];
           //正在计算回代过程
            for (int i = n; i >= 1; i--)
            {
                double sum = 0;
                for (int j = i + 1; j <= n; j++)
                {
                    sum += a[i, j] * x[j];
                }
                x[i] = (a[i, n + 1] - sum) / a[i, i];
            }
            
        }



        static void OutputX(double[] x)
        {
            int n = x.Length - 1;
            Console.WriteLine("方程解如下");
            for(int i=1;i<=n; i++)
            {
                Console.WriteLine("X[{0}]={1,8:f2}", i, x[i]);
            }
        }
        
        static void Main(string[] args)
        {

            Gauss obj = new Gauss();
            obj.Calcu ();

        }
    }
}

三、程序运行截图