逆矩阵（伴随矩阵法）C++

算法过程：
-计算 判断|A|是否为0
-利用原矩阵生成A*（伴随）矩阵，具体：A*二维数组中第[i][j]个元素，除去该行该列，其他元素进入临时数组，计算临时数组行列式值，即为A*[i][j]

最后矩阵A*/|A| 即为该矩阵的逆矩阵

源码：

#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
//By Vove.

float Cal_Det(float s[9][9],int n);//计算行列式|A|

class Matrix{//矩阵 方阵
public:
    friend void Cal_Astar(Matrix &A);//计算A*
    friend float Cal_Aij(Matrix &A,int i,int j);//计算Aij->A*
    Matrix();
    void Input();
    void Cheak_IsSolvable();
    float CalDet_A();
    void Inverse_of_A();
    void ToAstar(int j,int i,float s){
        Astar[j][i]=s;
    }
    float getAkl(int k,int l){
        return A[k][l];
    }
    void Dis_A(float s[9][9]);
    void Display(int t);
private:
    int n;              //行，列r=c
    float A[9][9];      //矩阵A
    float Astar[9][9];  // A*
    float Det_A;        // |A|值
    float IOA[9][9];    //inverse of A  A的逆矩阵
};


int main(){
    Matrix A;
    A.Input();
    A.Display(1);       //输出矩阵A
    A.Cheak_IsSolvable();//检查是否可解
    Cal_Astar(A);       //计算A*
    A.Inverse_of_A();   //计算A*/|A|
    A.Display(3);       //输出逆矩阵
    return 1;
}

void Cal_Astar(Matrix &A){
    for(int i=0;i<A.n;i++)
        for(int j=0;j<A.n;j++)
            A.ToAstar(j,i,pow(-1,i+j+2)*Cal_Aij(A,i,j));
        A.Display(2);
}

float Cal_Aij(Matrix &A,int i,int j){
    float Aij[9][9]={0};
    int x=0,y=0;                    //Aij的下标
    for(int k=0;k<A.n;k++){         //k l -> A的下标
        for(int l=0;l<A.n;l++){
            if(k!=i&&l!=j){
                    Aij[x][y]=A.getAkl(k,l);
                    y++;
                    if(!(y%((A.n)-1))){//进位
                        x++;
                        y=0;
                    }
            }
        }
    }
    return Cal_Det(Aij,(A.n)-1);
}

float Cal_Det(float s[9][9],int n){
    float sum=0; //存|A|值
    for(int i=0;i<n-1;i++){
        if(s[i][i]){//对角线元素不为0
            for(int m=i+1;m<n;m++){
                if(s[m][i]){
                    float temp=-(s[m][i])/s[i][i];
                    for(int p=0;p<n;p++)//R(m)+tempRi
                        s[m][p]=s[m][p]+temp*s[i][p];
                }
                else continue;
            }
        }
        else {//若对角线元素为0
            int m;
            for(m=i+1;m<n;m++){//
                if(s[m][i]){//使对角线元素非0
                    for(int p=0;p<n;p++)//Ri+Rm
                        s[i][p]=s[m][p]+s[i][p];
                    break;
                }
                else continue;
            }
            if(m==n) return sum;   //sum=0;
            i--;
        }
    }
    sum=s[0][0];//求对角积
    for(i=1;i<n;i++)
        sum=sum*s[i][i];
    return sum;
}

void Matrix::Dis_A(float s[9][9]){
    for(int i=0;i<n;i++){
        cout<<"|";
        for(int j=0;j<n;j++)
            cout<<setw(8)<<s[i][j];
        cout<<"    |"<<endl;
    }
    cout<<endl;
}
Matrix::Matrix(){
    n=9;
    for(int i=0;i<9;i++)
        for(int j=0;j<9;j++)
            A[i][j]=Astar[i][j]=IOA[i][j]=0;
}
void Matrix::Input(){
    cout<<"\t\t矩阵>>>>>>>>>>>逆矩阵"<<endl<<endl;
    cout<<"方阵A的行列数n:";
    while(1){
        cin>>n;
        if(n>1&&n<10) break;
        else cout<<"n值错误，重新输入：";
    }
    cout<<"矩阵各值:"<<endl;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            cin>>A[i][j];
}
void Matrix::Cheak_IsSolvable(){
    while(1){
        if(Det_A=CalDet_A()) break;
        else{
            cout<<"此矩阵的行列式为0，无解"<<endl;
            Input();
        } 
    }
}
float Matrix::CalDet_A(){
    float temp[9][9]={0};
    for(int i=0;i<n;i++)    //临时赋值
        for(int j=0;j<n;j++)
            temp[i][j]=A[i][j];

        float s=Cal_Det(temp,n);
        cout<<"|A|="<<s<<endl<<endl;
        return s;
}
void Matrix::Inverse_of_A(){
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            IOA[i][j]=Astar[i][j]/Det_A;
}
void Matrix::Display(int t){
    switch (t){
    case 1: cout<<"---A :"<<endl;Dis_A(A);break;
    case 2: cout<<"---A* :"<<endl;Dis_A(Astar);break;
    case 3: cout<<"---inverse of A :"<<endl;Dis_A(IOA);break;
    }
}

原博客：
http://blog.sina.com.cn/s/blog_155d4a9470102wji3.html