C语言解线性方程组.doc

实验一、用C语言解线性方程组 一、实验名称：用C语言解线性方程组。 二、实验原理：用初等行变化的方法将线性方程组的增广矩阵化简成最简型；从最简型中可以看出线性方程组解的情况。 三、实验工具：计算机和Microsoft Visual C++ 6.0软件。 四、实验原程序代码如下： #includestdio.h #includemath.h #includestdlib.h #define MAXR 20 #define MAXC 40 float M[MAXC]; float M_B[MAXR][MAXC]; void f1(int m,int n) {int i,j,r=0,c=0; float x0; printf(\n以下过程是将矩阵化为行阶梯型\n); printf(\n输出初始增广矩阵!\n); for(i=0;im;i++){ for(j=0;jn;j++)printf(%8.2f,M_B[i][j]); /*原始增广矩阵*/ printf(\n); } printf(按任一健继续...\n); getchar(); while(cn) {for(i=r;im;i++)if(fabs(M_B[i][c])=0.0001)break; if(im) { if(i!=r)for(j=0;jn;j++) { x0=M_B[i][j]; M_B[i][j]=M_B[r][j]; M_B[r][j]=x0; } x0=M_B[r][c]; for(j=0;jn;j++)M_B[r][j]/=x0; for(i=r+1;im;i++) { x0=M_B[i][c]; for(j=0;jn;j++)M_B[i][j]-=x0*M_B[r][j]; } r++; for(i=0;im;i++) { for(j=0;jn;j++)printf(%8.2f,M_B[i][j]); printf(\n); } printf(按任一健继续...\n); getchar(); } c++; if(r==m)break; } printf(\n矩阵的行阶梯型如上!\n); printf(\n以下过程是将行阶梯型化为行最简型\n); while(r) { r--; for(j=0;jn-1;j++)if(fabs(M_B[r][j])=0.0001)break; c=j; for(i=0;ir;i++) { x0=M_B[i][c]; for(j=0;jn;j++)M_B[i][j]-=x0*M_B[r][j]; } for(i=0;im;i++) { for(j=0;jn;j++)printf(%8.2f,M_B[i][j]); printf(\n); } printf(按任一健继续...\n); getchar(); } float M[MAXC]={0}; int l1,l2; for (l1=0;l1m;l1++) { for(l2=0;l2n;l2++)if(M_B[l1][l2]!=0){M[l2]=1;break;}} if(M[n-1]=1){printf(该方程无解 !);return;} for(l1=0;l1m;l1++) {for(l2=0;l2n;l2++){if(M_B[l1][l2]!=0){M[l2]=M_B[l1][n-1];}}} printf(输出线性方程组的解:); for(l2=0;l2n-1;l2++) printf(x%d(l2+1)=\n,M[l2]); printf(行最简型即为所求!按任一健退出...\n); getchar(); } void main (void) { int p,q,m,n; printf(\n 请输入线性方程组的行数m和列数n:); scanf(%d %d,m,n); printf(\n现在输入系数矩阵A和向量b:); for(p=1;p=m;p++) { printf(\n请输入a%d1--a%d%d系数和向量b%d: