高斯约旦求矩阵的逆

 

#include<iostream>
using namespace std;

#define n 3
int main()
{
    int i,j,k;
    double m;
    double a[n][n],E[n][n];

    //单位矩阵E[n][n]
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            if(i==j)
                E[i][j]=1;
            else
                E[i][j]=0;
        }
    }

    //输入矩阵a[n][n]
    cout<<"输入矩阵："<<endl;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            cin>>a[i][j];
        }
    }
    //上三角变换
    for(k=0;k<n-1;k++)
    {
        for(i=k+1;i<n;i++)
        {
            m=a[i][k]/a[k][k];
            for(j=0;j<n;j++)
            {
                a[i][j]=a[i][j]-m*a[k][j];
                E[i][j]=E[i][j]-m*E[k][j];
            }
        }
    }
    //下三角变换
    for(k=n-1;k>0;k--)
    {
        for(i=k-1;i>=0;i--)
        {
            m=a[i][k]/a[k][k];
            for(j=0;j<n;j++)
            {
                a[i][j]=a[i][j]-m*a[k][j];
                E[i][j]=E[i][j]-m*E[k][j];
            }
        }
    }
    //单位矩阵……
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            E[i][j]=E[i][j]/a[i][i];
        }
    }
    //输出矩阵的逆
    cout<<"矩阵的逆是："<<endl;
    for(i=0;i<n;i++)
    {
        for(j=0,k=0;j<n;j++)
        {
            cout<<E[i][j]<<"\t";k++;
            if(k==n)
                cout<<endl;
        }
    }
}