C#计算矩阵的秩实例分析

 更新时间:2015年08月13日 09:44:33   作者:北风其凉  
这篇文章主要介绍了C#计算矩阵的秩实现方法,以实例形式较为详细的分析了C#计算矩阵秩的原理与实现技巧,具有一定参考借鉴价值,需要的朋友可以参考下

本文实例讲述了C#计算矩阵的秩的方法。分享给大家供大家参考。具体如下:

1.代码思路

计算矩阵的秩,即把矩阵进行行初等变换,得出的行最简矩阵的非零行数。过程如下
1)将矩阵各行按第一个非零元素出现的位置升序排列(Operation1函数)
2)查看矩阵是否为行最简矩阵(isFinished函数),是则到第6步,不是则到第3步
3)如果有两行第一个非零元素出现的位置相同,则做消法变换,让下面行的第一个非零元素位置后移(Operation2函数)
4)将矩阵各行按第一个非零元素出现的位置升序排列(Operation1函数)
5)返回第2步
6)判断误差,对趋近与0的元素(如1E-5)按0处理,以免在第7步误判(Operation3函数)
7)统计非零行的数目(Operation4函数),即为矩阵的秩

2.函数代码

(注:本段代码只实现了一个思路,可能并不是该问题的最优解)

/// <summary>
/// 计算矩阵的秩
/// </summary>
/// <param name="matrix">矩阵</param>
/// <returns></returns>
private static int Rank(double[][] matrix)
{
  //matrix为空则直接默认已经是最简形式
  if (matrix == null || matrix.Length == 0) return 0;
  //复制一个matrix到copy,之后因计算需要改动矩阵时并不改动matrix本身
  double[][] copy = new double[matrix.Length][];
  for (int i = 0; i < copy.Length; i++)
  {
    copy[i] = new double[matrix[i].Length];
  }
  for (int i = 0; i < matrix.Length; i++)
  {
    for (int j = 0; j < matrix[0].Length; j++)
    {
      copy[i][j] = matrix[i][j];
    }
  }
  //先以最左侧非零项的位置进行行排序
  Operation1(copy);
  //循环化简矩阵
  while (!isFinished(copy))
  {
    Operation2(copy);
    Operation1(copy);
  }
  //过于趋近0的项,视作0,减小误差
  Operation3(copy);
  //行最简矩阵的秩即为所求
  return Operation4(matrix);
}
/// <summary>
/// 判断矩阵是否变换到最简形式(非零行数达到最少)
/// </summary>
/// <param name="matrix"></param>
/// <returns>true:</returns>
private static bool isFinished(double[][] matrix)
{
  //统计每行第一个非零元素的出现位置
  int[] counter = new int[matrix.Length];
  for (int i = 0; i < matrix.Length; i++)
  {
    for (int j = 0; j < matrix[i].Length; j++)
    {
      if (matrix[i][j] == 0)
      {
        counter[i]++;
      }
      else break;
    }
  }
  //后面行的非零元素出现位置必须在前面行的后面,全零行除外
  for (int i = 1; i < counter.Length; i++)
  {
    if (counter[i] <= counter[i - 1] && counter[i] != matrix[0].Length)
    {
      return false;
    }
  }
  return true;
}
/// <summary>
/// 排序(按左侧最前非零位位置自上而下升序排列)
/// </summary>
/// <param name="matrix">矩阵</param>
private static void Operation1(double[][] matrix)
{
  //统计每行第一个非零元素的出现位置
  int[] counter = new int[matrix.Length];
  for (int i = 0; i < matrix.Length; i++)
  {
    for (int j = 0; j < matrix[i].Length; j++)
    {
      if (matrix[i][j] == 0)
      {
        counter[i]++;
      }
      else break; 
    }
  }
  //按每行非零元素的出现位置升序排列
  for (int i = 0; i < counter.Length; i++)
  {
    for (int j = i; j < counter.Length; j++)
    {
      if(counter[i]>counter[j])
      {
        double[] dTemp = matrix[i];
        matrix[i] = matrix[j];
        matrix[j] = dTemp;
      }
    }
  }
}
/// <summary>
/// 行初等变换(左侧最前非零位位置最靠前的行,只保留一个)
/// </summary>
/// <param name="matrix">矩阵</param>
private static void Operation2(double[][] matrix)
{
  //统计每行第一个非零元素的出现位置
  int[] counter = new int[matrix.Length];
  for (int i = 0; i < matrix.Length; i++)
  {
    for (int j = 0; j < matrix[i].Length; j++)
    {
      if (matrix[i][j] == 0)
      {
        counter[i]++;
      }
      else break;
    }
  }
  for (int i = 1; i < counter.Length; i++)
  {
    if (counter[i] == counter[i - 1] && counter[i] != matrix[0].Length)
    {
      double a = matrix[i - 1][counter[i - 1]];
      double b = matrix[i][counter[i]]; //counter[i]==counter[i-1]
      matrix[i][counter[i]] = 0;
      for (int j = counter[i] + 1; j < matrix[i].Length; j++)
      {
        double c = matrix[i - 1][j];
        matrix[i][j] -= (c * b / a);
      }
      break;
    }
  }
}
/// <summary>
/// 将和0非常接近的数字视为0
/// </summary>
/// <param name="matrix"></param>
private static void Operation3(double[][] matrix)
{
  for (int i = 0; i < matrix.Length; i++)
  {
    for (int j = 0; j < matrix[0].Length; j++)
    {
      if (Math.Abs(matrix[i][j]) <= 0.00001)
      {
        matrix[i][j] = 0;
      }
    }
  }
}
/// <summary>
/// 计算行最简矩阵的秩
/// </summary>
/// <param name="matrix"></param>
/// <returns></returns>
private static int Operation4(double[][] matrix)
{
  int rank = -1;
  bool isAllZero = true;
  for (int i = 0; i < matrix.Length; i++)
  {
    isAllZero = true;
    //查看当前行有没有0
    for (int j = 0; j < matrix[0].Length; j++)
    {
      if (matrix[i][j] != 0)
      {
        isAllZero = false;
        break;
      }
    }
    //若第i行全为0,则矩阵的秩为i
    if (isAllZero)
    {
      rank = i;
      break;
    }
  }
  //满秩矩阵的情况
  if (rank == -1)
  {
    rank = matrix.Length;
  }
  return rank;
}

3.Main函数调用

static void Main(string[] args)
{
  //示例矩阵1:秩为3
  double[][] matrix1 = new double[][] 
  {
    new double[] { 1, 1, 1 },
    new double[] { 1, 1, 0 },
    new double[] { 0, 1, 1 } 
  };
  Console.WriteLine(Rank(matrix1));
  //示例矩阵2:秩为3
  double[][] matrix2 = new double[][] 
  {
    new double[] { 3, 2, 0, 5, 0 }, 
    new double[] { 3, -2, 3, 6, -1 },
    new double[] { 2, 0, 1, 5, -3 },
    new double[] { 1, 6, -4, -1, 4 } 
  };
  Console.WriteLine(Rank(matrix2));
  //示例矩阵3:秩为3
  double[][] matrix3 = new double[][] 
  {
    new double[] { 2, 3, 1, -3, -7 }, 
    new double[] { 1, 2, 0, -2, -4 },
    new double[] { 3, -2, 8, 3, 0 },
    new double[] { 2, -3, 7, 4, 3 }
  };
  Console.WriteLine(Rank(matrix3));
  Console.ReadLine();
}

4.执行结果

希望本文所述对大家的C#程序设计有所帮助。

相关文章

  • C# cmd中修改显示(显示进度变化效果)的方法

    C# cmd中修改显示(显示进度变化效果)的方法

    好多人想在运行或者调试含有大量数据或者比较慢C#程序的时候能够显示自己的程序完成的程度,这里有一个方法能发不断地修改cmd的同一行,以达到显示完成百分比的目的
    2013-04-04
  • 简单学习C#中的泛型方法使用

    简单学习C#中的泛型方法使用

    这篇文章主要介绍了C#中的泛型方法使用,需要的朋友可以参考下
    2016-02-02
  • 轻松学习C#的装箱与拆箱

    轻松学习C#的装箱与拆箱

    轻松学习C#的装箱与拆箱,在之前的文章简单的提到了轻松学习C#的装箱与拆箱,本文带着大家更加详细的介绍轻松学习C#的装箱与拆箱,感兴趣的小伙伴们可以参考一下
    2015-11-11
  • Unity制作图片字体的方法

    Unity制作图片字体的方法

    这篇文章主要为大家详细介绍了Unity制作图片字体的方法,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2020-12-12
  • C#程序员统计自己的代码行数

    C#程序员统计自己的代码行数

    这篇文章给大家讲解了下作为程序员如何统计自己写过的代码的行数,这个也是证明自己程序员能力的一个表现,一起来看下。
    2017-12-12
  • 深入DropDownList用法的一些学习总结分析

    深入DropDownList用法的一些学习总结分析

    本篇文章是对DropDownList的用法进行了详细的分析介绍,需要的朋友参考下
    2013-06-06
  • C#自定义日志记录

    C#自定义日志记录

    这篇文章主要为大家详细介绍了C#自定义日志记录的相关资料,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2016-09-09
  • C#正则匹配RegexOptions选项的组合使用方法

    C#正则匹配RegexOptions选项的组合使用方法

    本文主要简单介绍RegexOptions各种选项的作用,并介绍如何组合使用,为初学者解除一些疑惑。
    2016-04-04
  • C# Winform 禁止用户调整ListView的列宽

    C# Winform 禁止用户调整ListView的列宽

    在使用 ListView 的时候, 有时我们不想让别人随意调整列宽, 或者某几列的列宽, 以便达到美观, 或者隐藏数据的作用. 那么可以用一下代码来实现
    2011-05-05
  • C#中AS和IS关键字的用法

    C#中AS和IS关键字的用法

    这篇文章主要介绍了C#中AS和IS关键字的用法的相关资料,需要的朋友可以参考下
    2016-03-03

最新评论