sample4_4.cs

使用C#编辑的一些矩阵运算的源程序

/*
 * 示例程序Sample4_4: LEquations类的复系数方程组的全选主元高斯消去法
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
		[STAThread]
		static void Main(string[] args)
		{
            // 系数矩阵数据
            // 实部
            double[] mtxDataCoef3Real = {
                                            1.0,3.0,2.0,13.0,
                                            7.0,2.0,1.0,-2.0,
                                            9.0,15.0,3.0,-2.0,
                                            -2.0,-2.0,11.0,5.0};
            // 虚部
            double[] mtxDataCoef3Imag = {
                                            3.0,-2.0,1.0,6.0,
                                            -2.0,7.0,5.0,8.0,
                                            9.0,-3.0,15.0,1.0,
                                            -2.0,-2.0,7.0,6.0};
            // 常数矩阵数据
            // 实部
            double[] mtxDataConst3Real = {
                                             2.0,
                                             7.0,
                                             3.0,
                                             9.0};
            // 虚部
            double[] mtxDataConst3Imag = {
                                             1.0,
                                             2.0,
                                             -2.0,
                                             3.0};
		
            // 构造系数矩阵
            Matrix mtxCoef3Real = new Matrix(4, mtxDataCoef3Real);
            Matrix mtxCoef3Imag = new Matrix(4, mtxDataCoef3Imag);
            // 构造常数矩阵
            Matrix mtxConst3Real = new Matrix(4, 1, mtxDataConst3Real);
            Matrix mtxConst3Imag = new Matrix(4, 1, mtxDataConst3Imag);

            // 构造线性方程组
            LEquations leqs3 = new LEquations(mtxCoef3Real, mtxConst3Real);
		
            // 复系数方程组的全选主元高斯消去法
            Matrix mtxResult3Real = new Matrix();
            Matrix mtxResult3Imag = new Matrix();
            if (leqs3.GetRootsetGauss(mtxCoef3Imag, mtxConst3Imag, mtxResult3Real, mtxResult3Imag))
            {
                for (int i=0; i<mtxConst3Real.GetNumRows(); ++i)
                {
                    Complex cp = new Complex(mtxResult3Real.GetElement(i,0), mtxResult3Imag.GetElement(i,0));
                    Console.WriteLine(cp);
                }
            }
            else
            {
                Console.WriteLine("失败");
            }
        }
	}
}