实现 矩阵对称化和反对称化.docxVIP

C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法 2016/1/3 转自：/article/70996.htm 作者：北风其凉 字体：[ HYPERLINK javascript:turnbig() 增加 HYPERLINK javascript:turnsmall() 减小] 类型：转载 时间：2015-08-12 HYPERLINK /article/70996.htm \l comments 我要评论 这篇文章主要介绍了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法,较为详细的分析了矩阵分解运算的原理与C#实现技巧,需要的朋友可以参考下 本文实例讲述了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法。分享给大家供大家参考。具体如下： 理论依据 因为反对称矩阵满足 A^T = -A 设A = (aij) 则有 aii = -aii 所以 aii = 0. 即主对角线上元素全为0. 对任意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 /// 判断一个二维数组是否为矩阵 /// /su