您的位置首页百科问答

C#矩阵求逆的运算方法

C#矩阵求逆的运算方法

的有关信息介绍如下:

C#矩阵求逆的运算方法

在编程的处理数据过程中通常需要定义数组来实现矩阵运算,矩阵的求逆也经常会使用到,那么我们在使用C#时怎么求矩阵的逆呢。

双击运行visual studio 2015,文件-新建一个项目,可以是窗体也可以是控制台应用程序,看个人的编程题目要求而定。

14怎么创建一个C#窗体应用程序

定义一个矩阵(数组)a,c#中都是用数组来存储矩阵的,下面就求矩阵a的逆:

double[,] a = new double[3, 3] { { 1, 2, 3 }, { 4, 5, 6 }, { 7, 8, 9 } };

下面我们写一个类来实现矩阵的求逆(也可以写函数求逆),因为类只要调用就行了,而且不在主程序的代码中,可以简化我们的程序。创建方法如下:解决方案资源管理器-项目名称-添加-类。

类的名称改为“qiuni”。

在“求逆”类中添加矩阵求逆的代码(图片只截取了部分,全部代码如下):

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;

}

最后,回到主程序,定义矩阵(数组)b为a的逆矩阵,并调用我们的“求逆”类,就可以实现对矩阵的求逆运算了: