# [愚翁专栏]如何在C#去求矩阵的逆矩阵

private double[,] ReverseMatrix( double[,] dMatrix, int Level )

{

double dMatrixValue = MatrixValue( dMatrix, Level );

if( dMatrixValue == 0 ) return null;

double[,] dReverseMatrix = new double[Level,2*Level];

double x, c;

// Init Reverse matrix

for( int i = 0; i < Level; i++ )

{

for( int j = 0; j < 2 * Level; j++ )

{

if( j < Level )

dReverseMatrix[i,j] = dMatrix[i,j];

else

dReverseMatrix[i,j] = 0;

}

dReverseMatrix[i,Level + i ] = 1;

}

for( int i = 0, j = 0; i < Level && j < Level; i++, j++ )

{

if( dReverseMatrix[i,j] == 0 )

{

int m = i;

for( ; dMatrix[m,j] == 0; m++ );

if( m == Level )

return null;

else

{

// Add i-row with m-row

for( int n = j; n < 2 * Level; n++ )

dReverseMatrix[i,n] += dReverseMatrix[m,n];

}

}

// Format the i-row with "1" start

x = dReverseMatrix[i,j];

if( x != 1 )

{

for( int n = j; n < 2 * Level; n++ )

if( dReverseMatrix[i,n] != 0 )

dReverseMatrix[i,n] /= x;

}

// Set 0 to the current column in the rows after current row

for( int s = Level - 1; s > i;s-- )

{

x = dReverseMatrix[s,j];

for( int t = j; t < 2 * Level; t++ )

dReverseMatrix[s,t] -= ( dReverseMatrix[i,t]* x );

}

}

// Format the first matrix into unit-matrix

for( int i = Level - 2; i >= 0; i-- )

{

for( int j = i + 1 ; j < Level; j++ )

if( dReverseMatrix[i,j] != 0 )

{

c = dReverseMatrix[i,j];

for( int n = j; n < 2*Level; n++ )

dReverseMatrix[i,n] -= ( c * dReverseMatrix[j,n] );

}

}

double[,] dReturn = new double[Level, Level];

for( int i = 0; i < Level; i++ )

for( int j = 0; j < Level; j++ )

dReturn[i,j] = dReverseMatrix[i,j+Level];

return dReturn;

}

private double MatrixValue( double[,] MatrixList, int Level )

{

double[,] dMatrix = new double[Level, Level];

for( int i = 0; i < Level; i++ )

for( int j = 0; j < Level; j++ )

dMatrix[i,j] = MatrixList[i,j];

double c, x;

int k = 1;

for( int i = 0, j = 0; i < Level && j < Level; i++, j++ )

{

if( dMatrix[i,j] == 0 )

{

int m = i;

for( ; dMatrix[m,j] == 0; m++ );

if( m == Level )

return 0;

else

{

// Row change between i-row and m-row

for( int n = j; n < Level; n++ )

{

c = dMatrix[i,n];

dMatrix[i,n] = dMatrix[m,n];

dMatrix[m,n] = c;

}

// Change value pre-value

k *= (-1);

}

}

// Set 0 to the current column in the rows after current row

for( int s = Level - 1; s > i;s-- )

{

x = dMatrix[s,j];

for( int t = j; t < Level; t++ )

dMatrix[s,t] -= dMatrix[i,t]* ( x/dMatrix[i,j] );

}

}

double sn = 1;

for( int i = 0; i < Level; i++ )

{

if( dMatrix[i,i] != 0 )

sn *= dMatrix[i,i];

else

return 0;

}

return k*sn;

}

函数的调用如下：

double[,] dMatrix = new double[3,3]{{0,1,2},{1,0,1},{4,2,1}};

double[,] dReturn = ReverseMatrix( dMatrix, 3 );

if( dReturn != null )

{

for( int i=0; i < 3; i++ )

Debug.WriteLine( string.Format( "{0} {1} {2}",

dReturn[i,0], dReturn[i,1],dReturn[i,2] ) );

}

• 本文已收录于以下专栏：

## c# 矩阵求逆（转载）

1.代码思路1）对矩阵进行合法性检查：矩阵必须为方阵2）计算矩阵行列式的值（Determinant函数）3）只有满秩矩阵才有逆矩阵，因此如果行列式的值为0（在代码中以绝对值小于1E-6做判断），则终止...
• caimagic
• 2016年05月14日 18:35
• 1720

## 求逆矩阵的c++实现

• Zhangbei_
• 2015年08月11日 20:03
• 653

## 求解一个矩阵的逆矩阵

c++ #include using namespace std; const int N=20; const double E=1e-10; double calculate_A(double s...
• XerxesSimon
• 2017年04月12日 14:46
• 366

## Java之求逆矩阵

public class MatrixInverse { public static double Det(double [][]Matrix,int N)//计算n阶行列式（N=n-1） { ...
• u014581901
• 2016年03月04日 19:39
• 2331

## 手把手教你用Execel求解一个矩阵的逆矩阵

Excel中有很多执行线性代数运算的函数，只要轻点鼠标，各种计算都手到擒来。这一部分，我们将介绍矩阵求逆的方法。...
• baimafujinji
• 2015年11月08日 21:11
• 5976

## 矩阵伪逆介绍及C++/OpenCV/Eigen的三种实现

• fengbingchun
• 2017年06月06日 09:24
• 3692

## C#矩阵求逆

• wmmthu
• 2012年04月24日 08:55
• 5553

## 求矩阵的逆矩阵 c++

//计算A的逆矩阵保存到C中 void inverse_matrix(double A[],int n,double C[]) { int i,j,k,m=2*n; double mik,temp...
• hjq376247328
• 2015年12月26日 03:11
• 2720

## C++计算逆矩阵

#include #include using namespace std; double Fun(int n,double *ab) { double *aa = new double[n*...
• sinat_36219858
• 2017年01月04日 19:12
• 1127

## 4x4矩阵求逆

• birdflyto206
• 2016年04月07日 20:48
• 1143

举报原因： 您举报文章：[愚翁专栏]如何在C#去求矩阵的逆矩阵 色情 政治 抄袭 广告 招聘 骂人 其他 (最多只允许输入30个字)