Guass列主元消去法及其应用.docVIP

说明： 深 圳 大 学 实 验 报 告 实验课程名称： 计算机数值方法 实验项目名称： Guass列主元消去法及其应用 学院： 计算机软件学院 专业： 计算机科学与技术 报告人： 学号： 班级：04 指导教师： 实验时间： 2010.4.20 实验报告提交时间： 2010.5.3 实验目的与要求 了解Guass列主元消去法及其应用，掌握利用Guass列主元消去法解方程，解决实际问题 模型求解 0.x+0.4y+0.3z+60000=x 0.3x+0.1y+0.2z+30000=y 0.2x+0.1y+0.2z+50000=z 源代码： #include stdio.h #include stdlib.h #include math.h #define N 20 #define EPS 1.0e-4 double a[N][N],b[N],c[N][N],d[N][N]; void swap(double a,double b)//换行 { double temp; temp=a; a=b; a=temp; } void gauss(double a[][N],double b[N],int n)//gauss求解 { int i,j,k,l; double temp,m; for(k=0;k=n-2;k++) { temp=fabs(a[k][k]); //选主元 l=k; for(i=k+1;in;i++) { if(fabs(a[i][k]temp))//fabs为取绝对值函数。包含在头文件math.h中 { temp=fabs(a[i][k]); l=i; } } if(fabs(a[l][k])EPS) { printf(方程有无数组解!); exit(0); } if(l!=k) //换行 { for(j=k;j=n-1;j++) swap(a[l][j],a[k][j]); swap(b[l],b[k]); } for(i=k+1;in;i++) //按行进行消元 { m=a[i][k]/a[k][k];//行乘数 for(j=k;jn;j++) a[i][j]-=m*a[k][j]; b[i]-=m*b[k]; } } if(a[n-1][n-1]==0) //矩阵行列式值为0 { printf(矩阵的行列式为0!\n); exit(-1); } b[n-1]/=a[n-1][n-1]; for(i=n-2;i=0;i--) { for(j=n-1;ji;j--) b[i]-=a[i][j]*b[j]; b[i]/=a[i][i]; } } void main() { int i,j,n; double temp=0.0; printf(请输入要求解的方程数:); //矩阵维数为n scanf(%d,n); printf(请输入矩阵a[%d][%d]:\n,n,n); for(i=0;in;i++) //系数矩阵a for(j=0;jn;j++) scanf(%lf,a[i][j]); printf(请输入矩阵b[%d]:\n,n); for(i=0;in;i++) //矩阵b的输入 scanf(%lf,b[i]); gauss(a,b,n); //高斯求解法 for(i=0;in;i++) //输出结果 printf(X%d= %+.4f\n,i+1,b[i]); printf(请输入矩阵c[%d][%d]:\n,n,n); for(i=0;in;i++) //系数矩阵c for(j=0;jn;j++)