1、双击运行visual studio 2015,文件-新建一个项目,可以是窗体也可以是控制台应用程序,看个人的编程题目要求而定。
2、定义一个矩阵(数组)a,c#中都是用数组来存储矩阵的,下面就求矩阵a的逆:double[,] a = new doubl髫潋啜缅e[3, 3] { { 1, 2, 3 }, { 4, 5, 6 }, { 7, 8, 9 } };
3、下面我们写一个类来实现矩阵的求逆(也可以写函数求逆),因为类只要调用就行了,而且不在主程序的代码中,可以简化我们的程序。创建方法如下:解决方案资源管理器-项目名称-添加-类。
4、类的名称改为“qiuni”。
5、在“求逆”类中添加矩阵求逆的代码(图片只截取了部分,全部代码如下):public 衡痕贤伎static double[,] 求逆(double[荑樊综鲶,] A) { int n; n = A.GetLength(0); if (n != A.GetLength(1)) return null; double[,] C = new double[n, 2 * n]; double[,] D = new double[n, n]; for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) C[i, j] = A[i, j]; for (int i = 0; i < n; i++) C[i, i + n] = 1; for (int k = 0; k < n; k++) { double max = Math.Abs(C[k, k]); int ii = k; for (int m = k + 1; m < n; m++) if (max < Math.Abs(C[m, k])) { max = Math.Abs(C[m, k]); ii = m; } for (int m = k; m < 2 * n; m++) { if (ii == k) break; double c; c = C[k, m]; C[k, m] = C[ii, m]; C[ii, m] = c; } if (C[k, k] != 1) { double bs = C[k, k]; if (bs == 0) { Console.WriteLine("求逆错误!结果可能不正确!"); break; //return null; } C[k, k] = 1; for (int p = k + 1; p < n * 2; p++) { C[k, p] /= bs; } } for (int q = k + 1; q < n; q++) { double bs = C[q, k]; for (int p = k; p < n * 2; p++) { C[q, p] -= bs * C[k, p]; } } } for (int q = n - 1; q > 0; q--) { for (int k = q - 1; k > -1; k--) { double bs = C[k, q]; for (int m = k + 1; m < 2 * n; m++) { C[k, m] -= bs * C[q, m]; } } } for (int i = 0; i < n; i++) for (int j = 0; j < n; j++) D[i, j] = C[i, j + n]; return D; }
6、最后,回到主程序,定义矩阵(数组)b为a的逆矩阵,并调用我们的“求逆”类,就可以实现对矩阵的求逆运算了: