C#计算矩阵的秩实例分析
更新时间:2020年6月25日 11:27 点击:1764
本文实例讲述了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#实现xml文件的读取与写入简单实例
下一篇: C#求数组中元素全排列的方法
相关文章
- 我们在使用C#做项目的时候,基本上都需要制作登录界面,那么今天我们就来一步步看看,如果简单的实现登录界面呢,本文给出2个例子,由简入难,希望大家能够喜欢。...2020-06-25
- 这篇文章主要介绍了C# 字段和属性的的相关资料,文中示例代码非常详细,供大家参考和学习,感兴趣的朋友可以了解下...2020-11-03
- 这篇文章主要介绍了C#中截取字符串的的基本方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-11-03
- 本文给大家分享C#连接SQL数据库和查询数据功能的操作技巧,本文通过图文并茂的形式给大家介绍的非常详细,需要的朋友参考下吧...2021-05-17
- 这篇文章主要介绍了C#实现简单的Http请求的方法,以实例形式较为详细的分析了C#实现Http请求的具体方法,需要的朋友可以参考下...2020-06-25
- 本文主要介绍了C#中new的几种用法,具有很好的参考价值,下面跟着小编一起来看下吧...2020-06-25
使用Visual Studio2019创建C#项目(窗体应用程序、控制台应用程序、Web应用程序)
这篇文章主要介绍了使用Visual Studio2019创建C#项目(窗体应用程序、控制台应用程序、Web应用程序),小编觉得挺不错的,现在分享给大家,也给大家做个参考。一起跟随小编过来看看吧...2020-06-25- 这篇文章主要介绍了C#开发Windows窗体应用程序的简单操作步骤,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-04-12
- 这篇文章主要介绍了C#从数据库读取图片并保存的方法,帮助大家更好的理解和使用c#,感兴趣的朋友可以了解下...2021-01-16
- 最近做一个小项目不可避免的需要前端脚本与后台进行交互。由于是在asp.net中实现,故问题演化成asp.net中jiavascript与后台c#如何进行交互。...2020-06-25
- 这篇文章主要用实例讲解C#递归算法的概念以及用法,文中代码非常详细,帮助大家更好的参考和学习,感兴趣的朋友可以了解下...2020-06-25
- 本文通过例子,讲述了C++调用C#的DLL程序的方法,作出了以下总结,下面就让我们一起来学习吧。...2020-06-25
- 轻松学习C#的基础入门,了解C#最基本的知识点,C#是一种简洁的,类型安全的一种完全面向对象的开发语言,是Microsoft专门基于.NET Framework平台开发的而量身定做的高级程序设计语言,需要的朋友可以参考下...2020-06-25
- 本文主要介绍了C#变量命名规则小结,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2021-09-09
- 这篇文章主要介绍了c#中(&&,||)与(&,|)的区别详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-06-25
- 本文主要介绍了C# 中取绝对值的函数。具有很好的参考价值。下面跟着小编一起来看下吧...2020-06-25
- 这篇文章主要介绍了C#绘制曲线图的方法,以完整实例形式较为详细的分析了C#进行曲线绘制的具体步骤与相关技巧,具有一定参考借鉴价值,需要的朋友可以参考下...2020-06-25
- 这篇文章主要介绍了c#自带缓存使用方法,包括获取数据缓存、设置数据缓存、移除指定数据缓存等方法,需要的朋友可以参考下...2020-06-25
- 下面小编就为大家带来一篇C#学习笔记- 随机函数Random()的用法详解。小编觉得挺不错的,现在就分享给大家,也给大家做个参考。一起跟随小编过来看看吧...2020-06-25
- 这篇文章主要介绍了C#中list用法,结合实例形式分析了C#中list排序、运算、转换等常见操作技巧,具有一定参考借鉴价值,需要的朋友可以参考下...2020-06-25