反对称矩阵合同.docVIP

反对称矩阵合同 实现 矩阵对称化和反对称化 C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法 2016/1/3 转自：http:///article/70996.htm 作者：北风其凉 字体：[增加减小] 类型：转载 时间：2015-08-12 我要评论 这篇文章主要介绍了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法,较为详细的分析了矩阵分解运算的原理与C#实现技巧,需要的朋友可以参考下 本文实例讲述了C#实现将一个矩阵分解为对称矩阵与反称矩阵之和的方法。分享给大家供大 /// /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 (inti = 0; iresult.Length; i++) { result[i] = new double[matrix[i].Length]; } //对称矩阵为 (A+T(A))/2 其中A为原矩阵，T(A)为A的转置矩阵 for (inti = 0; iresult.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 (inti = 0; iresult.Length; i++) { result[i] = new double[matrix[i].Length]; } //反称矩阵为 (A-T(A))/2 其中A为原矩阵，T(A)为A的转置矩阵 for (inti = 0; iresult.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 /// returnstrue:是矩阵 false:不是矩阵/returns private static boolisMatrix(double[][] matrix) { //空矩阵是矩阵 if (matrix.Length 1) return true; //不同行列数如果不相等，则不是矩阵 int count = matrix[0].Length; for (inti = 1; imatrix.Length; i++) { if (matrix[i].Length != count) { return fal