C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法
更新时间:2015年08月12日 16:27:18 作者:北风其凉
这篇文章主要介绍了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法,较为详细的分析了矩阵分解运算的原理与C#实现技巧,需要的朋友可以参考下
本文实例讲述了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法。分享给大家供大家参考。具体如下:
1.理论依据
对任意n阶方阵A,有 A=(A+T(A))/2+(A-T(A))/2,其中T(A)是A的转置,(A+T(A))/2是一个对称矩阵,(A-T(A))/2是一个反称矩阵。
2.求出对称矩阵部分的函数
/// <summary> /// 把矩阵分解为对称矩阵与反称矩阵之和:对称矩阵 /// </summary> /// <param name="matrix">矩阵</param> /// <returns></returns> private static double[][] SymmetricPart(double[][] matrix) { //合法性校验:矩阵必须为方阵 if ( MatrixCR(matrix)[0] != MatrixCR(matrix)[1]) { throw new Exception("matrix 不是一个方阵"); } //矩阵中没有元素的情况 if (matrix.Length == 0) { return new double[][] { }; } //生成一个与matrix同型的空矩阵 double[][] result = new double[matrix.Length][]; for (int i = 0; i < result.Length; i++) { result[i] = new double[matrix[i].Length]; } //对称矩阵为 (A+T(A))/2 其中A为原矩阵,T(A)为A的转置矩阵 for (int i = 0; i < result.Length; i++) { for (int j = 0; j < result.Length; j++) { result[i][j] = (matrix[i][j] + matrix[j][i]) / 2.0; } } return result; }
3.求出反称矩阵部分的函数
/// <summary> /// 把矩阵分解为对称矩阵与反称矩阵之和:反称矩阵 /// </summary> /// <param name="matrix">矩阵</param> /// <returns></returns> private static double[][] SkewSymmetricPart(double[][] matrix) { //合法性校验:矩阵必须为方阵 if (MatrixCR(matrix)[0] != MatrixCR(matrix)[1]) { throw new Exception("matrix 不是一个方阵"); } //矩阵中没有元素的情况 if (matrix.Length == 0) { return new double[][] { }; } //生成一个与matrix同型的空矩阵 double[][] result = new double[matrix.Length][]; for (int i = 0; i < result.Length; i++) { result[i] = new double[matrix[i].Length]; } //反称矩阵为 (A-T(A))/2 其中A为原矩阵,T(A)为A的转置矩阵 for (int i = 0; i < result.Length; i++) { for (int j = 0; j < result.Length; j++) { result[i][j] = (matrix[i][j] - matrix[j][i]) / 2.0; } } return result; }
4.其他函数
/// <summary> /// 判断一个二维数组是否为矩阵 /// </summary> /// <param name="matrix">二维数组</param> /// <returns>true:是矩阵 false:不是矩阵</returns> private static bool isMatrix(double[][] matrix) { //空矩阵是矩阵 if (matrix.Length < 1) return true; //不同行列数如果不相等,则不是矩阵 int count = matrix[0].Length; for (int i = 1; i < matrix.Length; i++) { if (matrix[i].Length != count) { return false; } } //各行列数相等,则是矩阵 return true; } /// <summary> /// 计算一个矩阵的行数和列数 /// </summary> /// <param name="matrix">矩阵</param> /// <returns>数组:行数、列数</returns> private static int[] MatrixCR(double[][] matrix) { //接收到的参数不是矩阵则报异常 if (!isMatrix(matrix)) { throw new Exception("接收到的参数不是矩阵"); } //空矩阵行数列数都为0 if (!isMatrix(matrix) || matrix.Length == 0) { return new int[2] { 0, 0 }; } return new int[2] { matrix.Length, matrix[0].Length }; } /// <summary> /// 打印矩阵 /// </summary> /// <param name="matrix">待打印矩阵</param> private static void PrintMatrix(double[][] matrix) { for (int i = 0; i < matrix.Length; i++) { for (int j = 0; j < matrix[i].Length; j++) { Console.Write(matrix[i][j] + "\t"); //注意不能写为:Console.Write(matrix[i][j] + '\t'); } Console.WriteLine(); } }
5.Main函数代码及程序运行示例
static void Main(string[] args) { double[][] matrix = new double[][] { new double[] { 1, 2, 3 }, new double[] { 4, 5, 6 }, new double[] { 7, 8, 9 } }; Console.WriteLine("原矩阵"); PrintMatrix(matrix); Console.WriteLine("对称矩阵"); PrintMatrix(SymmetricPart(matrix)); Console.WriteLine("反称矩阵"); PrintMatrix(SkewSymmetricPart(matrix)); Console.ReadLine(); }
运行效果如下图所示:
希望本文所述对大家的C#程序设计有所帮助。
最新评论