C#实现将一个矩阵分解为对称矩阵与反称矩阵之和

一颗松果

发布时间：2026-10-13 | 1人浏览过

原创

本文将详细介绍如何在C#中将任意方阵分解为对称矩阵和反对称矩阵（也称斜对称矩阵），提供核心理论、完整代码实现及验证步骤。

理论依据

对于任意 n 阶方阵 A，都可以唯一地分解为一个对称矩阵和一个反对称矩阵的和：

A = (A + A^T) / 2 + (A - A^T) / 2

(A + A^T) / 2 是对称部分，满足 M = M^T。
(A - A^T) / 2 是反对称部分，满足 S = -S^T。
其中 A^T 表示矩阵 A 的转置。

C#代码实现

以下代码使用双重数组（double[][]）表示矩阵，并提供了从校验到计算再到打印的完整功能模块。

第1步：定义辅助函数与数据校验

目的：确保输入的是合法的二维矩阵，并能方便获取矩阵维度，防止运行时错误。

<summary>
/// 判断是否为有效矩阵
</summary>
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>
private static int[] MatrixCR(double[][] matrix)
{
 if (!isMatrix(matrix)) throw new Exception("接收到的参数不是矩阵");
 if (matrix.Length == 0) return new int[] { 0, 0 };
 return new int[] { matrix.Length, matrix[0].Length };
}

<summary>
/// 格式化打印矩阵
</summary>
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.WriteLine();
 }
}

第2步：实现对称矩阵的计算逻辑

目的：根据公式 (A + A^T) / 2 提取对称部分。注意：必须先校验矩阵是否为方阵。

private static double[][] SymmetricPart(double[][] matrix)
{
 //合法性校验：矩阵必须为方阵
 if (MatrixCR(matrix)[0] != MatrixCR(matrix)[1])
 throw new Exception("matrix 不是一个方阵");
 
 if (matrix.Length == 0)
 return new double[][] { };

 //生成同型空矩阵
 double[][] result = new double[matrix.Length][];
 for (int i = 0; i < result.Length; i++)
 result[i] = new double[matrix[i].Length];

 //计算对称部分：(A + Transpose(A)) / 2
 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步：实现反对称矩阵的计算逻辑

目的：根据公式 (A - A^T) / 2 提取反对称部分。

private static double[][] SkewSymmetricPart(double[][] matrix)
{
 //合法性校验：矩阵必须为方阵
 if (MatrixCR(matrix)[0] != MatrixCR(matrix)[1])
 throw new Exception("matrix 不是一个方阵");
 
 if (matrix.Length == 0)
 return new double[][] { };

 //生成同型空矩阵
 double[][] result = new double[matrix.Length][];
 for (int i = 0; i < result.Length; i++)
 result[i] = new double[matrix[i].Length];

 //计算反对称部分：(A - Transpose(A)) / 2
 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步：主函数验证与运行示例

目的：通过具体数值验证分解的正确性。我们将创建一个 3x3 矩阵，分别输出原矩阵、对称部分和反对称部分。

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("\n对称矩阵:");
 PrintMatrix(SymmetricPart(matrix));
 
 Console.WriteLine("\n反对称矩阵:");
 PrintMatrix(SkewSymmetricPart(matrix));

 Console.ReadLine();
}