sample3_19.cs

用于矩阵的各种运算

/*
 * 示例程序Sample3_19: Matrix类的求实对称矩阵特征值与特征向量的雅可比过关法
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
		[STAThread]
		static void Main(string[] args)
		{
            // 矩阵数据
            double[] mtxData20 = {
                                     10.0,1.0,2.0,3.0,4.0,
                                     1.0,9.0,-1.0,2.0,-3.0,
                                     2.0,-1.0,7.0,3.0,-5.0,
                                     3.0,2.0,3.0,12.0,-1.0,
                                     4.0,-3.0,-5.0,-1.0,15.0};

            // 构造矩阵
            Matrix mtx20 = new Matrix(5, 5, mtxData20);

            // 求实对称矩阵特征值与特征向量的雅可比过关法
            double[] valArray1 = new double[mtx20.GetNumColumns()];
            Matrix mtxVV1 = new Matrix();
            if (mtx20.ComputeEvJacobi(valArray1, mtxVV1, 0.000001))
            {
                Console.WriteLine("特征值=");
                String s = "";
                for (int i=0; i<mtx20.GetNumColumns(); ++i)
                {
                    s += valArray1[i] + ", ";
                }
                Console.WriteLine(s);
                Console.WriteLine("-------------------------------"); 
                Console.WriteLine("对应的特征向量=");
                for (int i=0; i<mtx20.GetNumColumns(); ++i)
                {
                    Console.WriteLine(mtxVV1.ToStringCol(i, ", "));
                }
            }
            else
            {
                Console.WriteLine("失败");
            }
        }
	}
}