C#计算矩阵的逆矩阵方法实例分析
更新时间:2020年6月25日 11:27 点击:2501
本文实例讲述了C#计算矩阵的逆矩阵方法。分享给大家供大家参考。具体如下:
1.代码思路
1)对矩阵进行合法性检查:矩阵必须为方阵
2)计算矩阵行列式的值(Determinant函数)
3)只有满秩矩阵才有逆矩阵,因此如果行列式的值为0(在代码中以绝对值小于1E-6做判断),则终止函数,报出异常
4)求出伴随矩阵(AdjointMatrix函数)
5)逆矩阵各元素即其伴随矩阵各元素除以矩阵行列式的商
2.函数代码
(注:本段代码只实现了一个思路,可能并不是该问题的最优解)
/// <summary> /// 求矩阵的逆矩阵 /// </summary> /// <param name="matrix"></param> /// <returns></returns> public static double[][] InverseMatrix(double[][] matrix) { //matrix必须为非空 if (matrix == null || matrix.Length == 0) { return new double[][] { }; } //matrix 必须为方阵 int len = matrix.Length; for (int counter = 0; counter < matrix.Length; counter++) { if (matrix[counter].Length != len) { throw new Exception("matrix 必须为方阵"); } } //计算矩阵行列式的值 double dDeterminant = Determinant(matrix); if (Math.Abs(dDeterminant) <= 1E-6) { throw new Exception("矩阵不可逆"); } //制作一个伴随矩阵大小的矩阵 double[][] result = AdjointMatrix(matrix); //矩阵的每项除以矩阵行列式的值,即为所求 for (int i = 0; i < matrix.Length; i++) { for (int j = 0; j < matrix.Length; j++) { result[i][j] = result[i][j] / dDeterminant; } } return result; } /// <summary> /// 递归计算行列式的值 /// </summary> /// <param name="matrix">矩阵</param> /// <returns></returns> public static double Determinant(double[][] matrix) { //二阶及以下行列式直接计算 if (matrix.Length == 0) return 0; else if (matrix.Length == 1) return matrix[0][0]; else if (matrix.Length == 2) { return matrix[0][0] * matrix[1][1] - matrix[0][1] * matrix[1][0]; } //对第一行使用“加边法”递归计算行列式的值 double dSum = 0, dSign = 1; for (int i = 0; i < matrix.Length; i++) { double[][] matrixTemp = new double[matrix.Length - 1][]; for (int count = 0; count < matrix.Length - 1; count++) { matrixTemp[count] = new double[matrix.Length - 1]; } for (int j = 0; j < matrixTemp.Length; j++) { for (int k = 0; k < matrixTemp.Length; k++) { matrixTemp[j][k] = matrix[j + 1][k >= i ? k + 1 : k]; } } dSum += (matrix[0][i] * dSign * Determinant(matrixTemp)); dSign = dSign * -1; } return dSum; } /// <summary> /// 计算方阵的伴随矩阵 /// </summary> /// <param name="matrix">方阵</param> /// <returns></returns> public static double[][] AdjointMatrix(double [][] matrix) { //制作一个伴随矩阵大小的矩阵 double[][] result = new double[matrix.Length][]; for (int i = 0; i < result.Length; i++) { result[i] = new double[matrix[i].Length]; } //生成伴随矩阵 for (int i = 0; i < result.Length; i++) { for (int j = 0; j < result.Length; j++) { //存储代数余子式的矩阵(行、列数都比原矩阵少1) double[][] temp = new double[result.Length - 1][]; for (int k = 0; k < result.Length - 1; k++) { temp[k] = new double[result[k].Length - 1]; } //生成代数余子式 for (int x = 0; x < temp.Length; x++) { for (int y = 0; y < temp.Length; y++) { temp[x][y] = matrix[x < i ? x : x + 1][y < j ? y : y + 1]; } } //Console.WriteLine("代数余子式:"); //PrintMatrix(temp); result[j][i] = ((i + j) % 2 == 0 ? 1 : -1) * Determinant(temp); } } //Console.WriteLine("伴随矩阵:"); //PrintMatrix(result); return result; } /// <summary> /// 打印矩阵 /// </summary> /// <param name="matrix">待打印矩阵</param> private static void PrintMatrix(double[][] matrix, string title = "") { //1.标题值为空则不显示标题 if (!String.IsNullOrWhiteSpace(title)) { Console.WriteLine(title); } //2.打印矩阵 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(); } //3.空行 Console.WriteLine(); }
3.Main函数调用
static void Main(string[] args) { double[][] matrix = new double[][] { new double[] { 1, 2, 3 }, new double[] { 2, 2, 1 }, new double[] { 3, 4, 3 } }; PrintMatrix(matrix, "原矩阵"); PrintMatrix(AdjointMatrix(matrix), "伴随矩阵"); Console.WriteLine("行列式的值为:" + Determinant(matrix) + '\n'); PrintMatrix(InverseMatrix(matrix), "逆矩阵"); Console.ReadLine(); }
4.执行结果
希望本文所述对大家的C#程序设计有所帮助。
上一篇: C#实现Ping的方法小结
下一篇: 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#绘制曲线图的方法,以完整实例形式较为详细的分析了C#进行曲线绘制的具体步骤与相关技巧,具有一定参考借鉴价值,需要的朋友可以参考下...2020-06-25
- 本文主要介绍了C# 中取绝对值的函数。具有很好的参考价值。下面跟着小编一起来看下吧...2020-06-25
- 这篇文章主要介绍了c#自带缓存使用方法,包括获取数据缓存、设置数据缓存、移除指定数据缓存等方法,需要的朋友可以参考下...2020-06-25
- 下面小编就为大家带来一篇C#学习笔记- 随机函数Random()的用法详解。小编觉得挺不错的,现在就分享给大家,也给大家做个参考。一起跟随小编过来看看吧...2020-06-25
- 这篇文章主要介绍了C#中list用法,结合实例形式分析了C#中list排序、运算、转换等常见操作技巧,具有一定参考借鉴价值,需要的朋友可以参考下...2020-06-25