sample3_12.cs

来自「C#的数值计算的几乎所有算法」· CS 代码 · 共 46 行

/*
 * 示例程序Sample3_12: Matrix类的一般实矩阵的奇异值分解
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
		[STAThread]
		static void Main(string[] args)
		{
            // 矩阵数据
            double[] mtxData13 = {
                                     1.0,1.0,-1.0,
                                     2.0,1.0,0.0,
                                     1.0,-1.0,0.0,
                                     -1.0,2.0,1.0};

            // 构造矩阵
            Matrix mtx13 = new Matrix(4, 3, mtxData13);

            // 一般实矩阵的奇异值分解
            Matrix mtxU1 = new Matrix();
            Matrix mtxV = new Matrix();
            if (mtx13.SplitUV(mtxU1, mtxV, 0.0001))
            {
                Console.WriteLine("分解后的U矩阵=");
                Console.WriteLine(mtxU1);
                Console.WriteLine("-------------------------------"); 
                Console.WriteLine("分解后的V矩阵=");
                Console.WriteLine(mtxV);
                Console.WriteLine("-------------------------------"); 
                Console.WriteLine("分解后的A矩阵=");
                Console.WriteLine(mtx13);
            }
            else
            {
                Console.WriteLine("失败");
            }
        }
	}
}