C#實(shí)現(xiàn)將一個(gè)矩陣分解為對(duì)稱矩陣與反稱矩陣之和的方法

更新時(shí)間：2015年08月12日 16:27:18 作者：北風(fēng)其涼

這篇文章主要介紹了C#實(shí)現(xiàn)將一個(gè)矩陣分解為對(duì)稱矩陣與反稱矩陣之和的方法,較為詳細(xì)的分析了矩陣分解運(yùn)算的原理與C#實(shí)現(xiàn)技巧,需要的朋友可以參考下

本文實(shí)例講述了C#實(shí)現(xiàn)將一個(gè)矩陣分解為對(duì)稱矩陣與反稱矩陣之和的方法。分享給大家供大家參考。具體如下：

1.理論依據(jù)

對(duì)任意n階方陣A，有 A=(A+T(A))/2+(A-T(A))/2，其中T(A)是A的轉(zhuǎn)置，(A+T(A))/2是一個(gè)對(duì)稱矩陣，(A-T(A))/2是一個(gè)反稱矩陣。

2.求出對(duì)稱矩陣部分的函數(shù)

/// <summary>
/// 把矩陣分解為對(duì)稱矩陣與反稱矩陣之和：對(duì)稱矩陣
/// </summary>
/// <param name="matrix">矩陣</param>
/// <returns></returns>
private static double[][] SymmetricPart(double[][] matrix)
{
 //合法性校驗(yàn)：矩陣必須為方陣
 if ( MatrixCR(matrix)[0] != MatrixCR(matrix)[1])
 {
  throw new Exception("matrix 不是一個(gè)方陣");
 }
 //矩陣中沒(méi)有元素的情況
 if (matrix.Length == 0)
 {
  return new double[][] { };
 }
 //生成一個(gè)與matrix同型的空矩陣
 double[][] result = new double[matrix.Length][];
 for (int i = 0; i < result.Length; i++)
 {
  result[i] = new double[matrix[i].Length];
 }
 //對(duì)稱矩陣為 (A+T(A))/2 其中A為原矩陣，T(A)為A的轉(zhuǎn)置矩陣
 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.求出反稱矩陣部分的函數(shù)

/// <summary>
/// 把矩陣分解為對(duì)稱矩陣與反稱矩陣之和：反稱矩陣
/// </summary>
/// <param name="matrix">矩陣</param>
/// <returns></returns>
private static double[][] SkewSymmetricPart(double[][] matrix)
{
 //合法性校驗(yàn)：矩陣必須為方陣
 if (MatrixCR(matrix)[0] != MatrixCR(matrix)[1])
 {
  throw new Exception("matrix 不是一個(gè)方陣");
 }
 //矩陣中沒(méi)有元素的情況
 if (matrix.Length == 0)
 {
  return new double[][] { };
 }
 //生成一個(gè)與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的轉(zhuǎn)置矩陣
 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.其他函數(shù)

/// <summary>
/// 判斷一個(gè)二維數(shù)組是否為矩陣
/// </summary>
/// <param name="matrix">二維數(shù)組</param>
/// <returns>true:是矩陣 false:不是矩陣</returns>
private static bool isMatrix(double[][] matrix)
{
 //空矩陣是矩陣
 if (matrix.Length < 1) return true;
 //不同行列數(shù)如果不相等，則不是矩陣
 int count = matrix[0].Length;
 for (int i = 1; i < matrix.Length; i++)
 {
  if (matrix[i].Length != count)
  {
   return false;
  }
 }
 //各行列數(shù)相等，則是矩陣
 return true;
}
/// <summary>
/// 計(jì)算一個(gè)矩陣的行數(shù)和列數(shù)
/// </summary>
/// <param name="matrix">矩陣</param>
/// <returns>數(shù)組：行數(shù)、列數(shù)</returns>
private static int[] MatrixCR(double[][] matrix)
{
 //接收到的參數(shù)不是矩陣則報(bào)異常
 if (!isMatrix(matrix))
 {
  throw new Exception("接收到的參數(shù)不是矩陣");
 }
 //空矩陣行數(shù)列數(shù)都為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函數(shù)代碼及程序運(yùn)行示例

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("對(duì)稱矩陣");
 PrintMatrix(SymmetricPart(matrix));
 Console.WriteLine("反稱矩陣");
 PrintMatrix(SkewSymmetricPart(matrix));
 Console.ReadLine();
}

運(yùn)行效果如下圖所示：