矩陣求逆算法及程序實現（C++）

　　在做課題時，遇到了求多項式問題，利用了求逆方法。矩陣求逆一般使用簡單的算法，還有快速算法 如全選主元高斯-約旦消元法，但本文程序主要寫了簡單的矩陣求逆算法定義法之伴隨矩陣求逆公式如下，其中A可逆：

　　　　　　　　　　　　　　　　，其中是的伴隨矩陣。。

　　1.給定一個方陣，非奇異（不是也可，程序有考慮）；

　　2.由矩陣得到其行列式，求其值如|A|；

　　3.求其伴隨矩陣；

　　4.得到其逆矩陣。

主要函數如下：

//得到給定矩陣src的逆矩陣保存到des中。
bool GetMatrixInverse(double src[N][N],int n,double des[N][N])
{
    double flag=getA(src,n);
    double t[N][N];
    if(flag==0)
    {
        return false;
    }
    else
    {
        getAStart(src,n,t);
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                des[i][j]=t[i][j]/flag;
            }

        }
    }


    return true;

}

計算|A|:

//按第一行展開計算|A|
double getA(double arcs[N][N],int n)
{
    if(n==1)
    {
        return arcs[0][0];
    }
    double ans = 0;
    double temp[N][N]={0.0};
    int i,j,k;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n-1;j++)
        {
            for(k=0;k<n-1;k++)
            {
                temp[j][k] = arcs[j+1][(k>=i)?k+1:k];

            }
        }
        double t = getA(temp,n-1);
        if(i%2==0)
        {
            ans += arcs[0][i]*t;
        }
        else
        {
            ans -=  arcs[0][i]*t;
        }
    }
    return ans;
}

計算伴隨矩陣：

//計算每一行每一列的每個元素所對應的余子式，組成A*
void  getAStart(double arcs[N][N],int n,double ans[N][N])
{
    if(n==1)
    {
        ans[0][0] = 1;
        return;
    }
    int i,j,k,t;
    double temp[N][N];
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            for(k=0;k<n-1;k++)
            {
                for(t=0;t<n-1;t++)
                {
                    temp[k][t] = arcs[k>=i?k+1:k][t>=j?t+1:t];
                }
            }


            ans[j][i]  =  getA(temp,n-1);
            if((i+j)%2 == 1)
            {
                ans[j][i] = - ans[j][i];
            }
        }
    }
}

　　這三個函數組成了求逆矩陣的步驟，筆者調用並成功運用，無錯，在此共勉。如有問題，請留言。謝謝

　　另有其他矩陣求逆好算法及代碼請大家共享。