高斯-约当消去法解线性方程组.docVIP

计算机科学与技术 蔡梦卿 学号高斯—约当消去法解线性方程组 源代码： #includeiostream #includecmath #includeiomanip #define kk 50 //定义最大方程元数 int n; using namespace std; void changer(double *p,double *q); main() { int i,j,k,s; double A[kk][kk],r,l[kk][kk],b[kk],u,temp; aa: cout输入的方程元数endl; cinn; cout请输入方程系数矩阵:\n; for(i=0;in;i++) for(j=0;jn;j++) cinA[i][j]; cout请输入右边向量:\n; for(i=0;in;i++) cinb[i]; for(k=0;kn-1;k++) { s=k; r=fabs(A[k][k]); for(i=k;in;i++) if(rfabs(A[i][k])) { r=fabs(A[i][k]); s=i; }//获取主元 if(s!=k) { changer(*(A+k),*(A+s));//交换行元素 temp=b[k];b[k]=b[s];b[s]=temp; //右边向量的交换 } cout第k+1次选取主元并交换行的顺序endl; for (i=0;in;i++) { for (j=0;jn;j++) coutsetw(4)setprecision(4)A[i][j] ; coutb[i]endl; }// coutendl; for(i=k+1;in;i++) { l[i][k]=A[i][k]/A[k][k];//计算消去系数 for(j=k;jn;j++) A[i][j]=A[i][j]-l[i][k]*A[k][j]; b[i]=b[i]-l[i][k]*b[k]; }//消元过程 cout第k+1次消元endl; for (int e=0;en;e++) { for (j=0;jn;j++) coutsetw(4)setprecision(4)A[e][j] ; coutb[e]endl; } coutendl; } if(A[k][k]0.0001) { cout系数矩阵是非奇异矩阵!\n; cout所以：方程无唯一解!\n; cout请问是否？endl; cout退出请输入1endl; cout要继续请输入任意非1的整数endl; int g; cing; if(g!=1)goto aa; else exit(1); }//判断系数矩阵的奇异性 for(i=n-1;i=0;i--) { u=0; for(j=i+1;jn;j++) u=u+A[i][j]*b[j]; b[i]=(b[i]-u)/A[i][i]; } cout最后得方程的根为:\n; for(i=0;in;i++) coutxi+1=b[i]\n; system(pause); return 0; } void changer(double *p,double *q) { int i; double temp; for(i=0;in;i++) { temp=*(p+i); *(p+i)=*(q+i); *(q+i)=temp; } } //交换行元素的函数 例题：方程组： 3X1 ＋ 4X2 － X3 ＋ X4 ≒ 38 ﹣X1 ＋ 2X2 ＋ 5X3 － X4 ≒ 38 X1 － 4X2 ＋ 3X3 ＋2X4 ≒﹣3 2X1 － X2 ＋ 4X3 ＋5X4 ≒ 32 X1≒5, X2≒7, X3≒6, X4≒1