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

说明：深 圳 大 学 实 验 报 告实验课程名称： 计算机数值方法 实验项目名称： Guass列主元消去法及其应用 学院： 计算机软件学院 专业： 计算机科学与技术 报告人： 学号： 班级：04 指导教师： 实验时间： 2010.4.20 实验报告提交时间： 2010.5.3 一、 实验目的与要求了解Guass列主元消去法及其应用，掌握利用Guass列主元消去法解方程，解决实际问题二、 模型求解0.x+0.4y+0.3z+60000=x0.3x+0.1y+0.2z+30000=y0.2x+0.1y+0.2z+50000=z三、 源代码：#include #include #include #define N 20#define EPS 1.0e-4double aNN,bN,cNN,dNN;void swap(double &a,double &b)/换行double temp;temp=a;a=b;a=temp;void gauss(double aN,double bN,int n)/gauss求解int i,j,k,l;double temp,m;for(k=0;k=n-2;k+)temp=fabs(akk); /选主元l=k;for(i=k+1;itemp)/fabs为取绝对值函数。包含在头文件math.h中 temp=fabs(aik); l=i;if(fabs(alk)EPS) printf(方程有无数组解!);exit(0);if(l!=k) /换行for(j=k;j=n-1;j+)swap(alj,akj);swap(bl,bk);for(i=k+1;in;i+) /按行进行消元m=aik/akk;/行乘数for(j=k;j=0;i-) for(j=n-1;ji;j-)bi-=aij*bj; bi/=aii;void main()int i,j,n;double temp=0.0;printf(请输入要求解的方程数:); /矩阵维数为nscanf(%d,&n);printf(请输入矩阵a%d%d:n,n,n); for(i=0;in;i+) /系数矩阵a for(j=0;jn;j+)scanf(%lf,&aij);printf(请输入矩阵b%d:n,n); for(i=0;in;i+) /矩阵b的输入 scanf(%lf,&bi);gauss(a,b,n); /高斯求解法for(i=0;in;i+) /输出结果 printf(X%d= %+.4fn,i+1,bi);printf(请输入矩阵c%d%d:n,n,n); for(i=0;in;i+) /系数矩阵c for(j=0;jn;j+)scanf(%lf,&cij);for(i=0;in;i+) for(j=0;jn;j+)dji=cji*bi;for(i=0;in;i+) for(j=0;jn;j+) printf(%.4f,dij); printf(n);四、 实验小结：通过本次实验，对如何利用Guass列主元消去法解决实际问题有了一定的经验。指导教师批阅意见：成绩评定： 指导教师签字： 年 月