matlab求若尔当标准型,Gauss-Jordan消元法求解线性方程组

#include

#include

#include

#define N 3

void main()

{

double a[N][N]={{1,2,3},{2,4,5},{3,5,6}};//线性方程组的系数矩阵

double b[N]={1,0,0};

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

double T,BT=0;//存储换行临时变量

double temp=0;

double t;//最大列主元

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

cout<

for(int i=0;i

{

for(int j=0;j

{

cout<

}

cout<

}

cout<

for( i=0;i

{

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=s;i

{

T=a[tap1][i];

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

a[s][i]=T;

}

BT=b[tap1];//给B换行

b[tap1]=b[s];

b[s]=BT;

}

tap1=tap2=0;//置零

/*cout<

for(int j=s;j

{

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

}

b[s]=b[s]/t;

for(int i=0;i

{

if(i!=s){

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

for(int j=s;j

{

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

}

b[i]=b[i]-m[i]*b[s];}

}

/*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<

}

cout<

for( i=0;i

{

cout<

}

cout<

}