matlab列主元消去法求逆矩阵,Gauss-Jordan列主元消元法求逆矩阵

#include

#include

#include

#define N 3

void main()

{

double a[N][N]={{1,1,-1},{2,1,0},{1,-1,0}};//线性方程组的系数矩阵

double c[N][N];//求逆矩阵

double m[N];//辅助乘数

double T;//存储换行临时变量

double temp=0;

double t;//最大列主元

int tap1=0,tap2=0;//最大列主元下标

for(int i=0;i

{

for(int j=0;j

{

if(i==j)

{

c[i][j]=1;

}

else

{

c[i][j]=0;

}

}

}

//for(i=0;i

cout<

for(i=0;i

{

for(int j=0;j

{

cout<

}

cout<

}

cout<

for(int s=0;s

{

t=a[s][s];//赋初值

for(int p=s;p

{

if(fabs(a[p][s])>t)

{

t=a[p][s];

tap1=p;

tap2=s;

}

}

//cout<

if(t==0)

{cout<

if(tap1!=tap2)//换行

{

for(int i=0;i

{

T=a[tap1][i];

a[tap1][i]=a[s][i];

a[s][i]=T;

T=c[tap1][i];//逆矩阵换行

c[tap1][i]=c[s][i];

c[s][i]=T;

}

}

tap1=tap2=0;//置零

/*cout<

cout<

for(int j=0;j

{

a[s][j]=a[s][j]/t;

c[s][j]=c[s][j]/t;

}

for(int i=0;i

{

if(i!=s)

{

m[i]=a[i][s]/a[s][s];

for(int j=0;j

{

a[i][j]=a[i][j]-m[i]*a[s][j];

c[i][j]=c[i][j]-m[i]*c[s][j];

}

}

}

/*cout<

}

cout<

for(i=0;i

{

for(int j=0;j

{

cout<

}

cout<

}

cout<

cout<

for(i=0;i

{

for(int j=0;j

{

cout<

}

cout<

}

}