如何在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++ )

    {