C#矩阵求逆

来源：http://zhidao.baidu.com/link?url=DiqAbq9YUYn3z7QjxGGoF0PLZwN-Y9ecqKB7Gy38JWRD1riMIYukVKXKq88pxtWLwIl6-iOUf2p3liE51phEe_

private double[,] ReverseMatrix( double[,] dMatrix, int Level )
{
    double dMatrixValue = MatrixValue( dMatrix, Level );
    if( dMatrixValue ==  ) return null;
 
    double[,] dReverseMatrix = new double[Level,*Level];
    double x, c;
    // Init Reverse matrix
    for( int i = ; i < Level; i++ )
    {
        for( int j = ; j <  * Level; j++ )
        {
            if( j < Level )
                dReverseMatrix[i,j] = dMatrix[i,j];
            else
                dReverseMatrix[i,j] = ;
        }
 
        dReverseMatrix[i,Level + i ] = ;
    }
 
    for( int i = , j = ; i < Level && j < Level; i++, j++ )
    {
        if( dReverseMatrix[i,j] ==  )
        {
            int m = i;
            for( ; dMatrix[m,j] == ; m++ );
            if( m == Level )
                return null;
            else
            {
                // Add i-row with m-row
                for( int n = j; n <  * Level; n++ )
                    dReverseMatrix[i,n] += dReverseMatrix[m,n];
            }
        }
 
        // Format the i-row with "1" start
        x = dReverseMatrix[i,j];
        if( x !=  )
        {
            for( int n = j; n <  * Level; n++ )
                if( dReverseMatrix[i,n] !=  )
                    dReverseMatrix[i,n] /= x;
        }
 
        // Set 0 to the current column in the rows after current row
        for( int s = Level - ; s > i;s-- )
        {
            x = dReverseMatrix[s,j];
            for( int t = j; t <  * Level; t++ )
                dReverseMatrix[s,t] -= ( dReverseMatrix[i,t]* x );
        }
    }
 
    // Format the first matrix into unit-matrix
    for( int i = Level - ; i >= ; i-- )
    {
        for( int j = i +  ; j < Level; j++ )
            if( dReverseMatrix[i,j] !=  )
            {
                c = dReverseMatrix[i,j];
                for( int n = j; n < *Level; n++ )
                    dReverseMatrix[i,n] -= ( c * dReverseMatrix[j,n] );
            }
    }
 
    double[,] dReturn = new double[Level, Level];
    for( int i = ; i < Level; i++ )
        for( int j = ; 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 = ; i < Level; i++ )
        for( int j = ; j < Level; j++ )
            dMatrix[i,j] = MatrixList[i,j];
    double c, x;
    int k = ;
    for( int i = , j = ; i < Level && j < Level; i++, j++ )
    {
        if( dMatrix[i,j] ==  )
        {
            int m = i;
            for( ; dMatrix[m,j] == ; m++ );
            if( m == Level )
                return ;
            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 *= (-);
            }
        }
 
        // Set 0 to the current column in the rows after current row
        for( int s = Level - ; 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 = ;
    for( int i = ; i < Level; i++ )
    {
        if( dMatrix[i,i] !=  )
            sn *= dMatrix[i,i];
        else
            return ;
    }
    return k*sn;
}

调用如下：

double[,] dMatrix = new double[,]{{,,},{,,},{,,}};
    double[,] dReturn = ReverseMatrix( dMatrix,  );
    if( dReturn != null )
    {
        for( int i=; i < ; i++ )
            Debug.WriteLine( string.Format( "{0} {1} {2}",
                dReturn[i,], dReturn[i,],dReturn[i,] ) );
    }