高斯列主元素消去法.docVIP

2013-2014（1）专业课程实践论文 题目：高斯列主元素消去法 队员1：马 越，0918180216，R数学09-2 队员2：徐永刚，0918180214，R数学09-2 专业课程实践论文成绩考核表 学生姓名1 马越 专业班级学号 R数学09-2 16 学生姓名2 徐永刚 专业班级学号 R数学09-2 14 学生姓名3 专业班级学号 题 目 高斯列主元素消去法 评 审 考 核 项 目 评分（每项 满分20分） 学生1 学生 2 学生3 1 对平时出勤、态度、道德素质进行成绩评定 2 根据学生对算法编程与调试熟练程度进行成绩评定（三人相同） 3 思路清晰、表达准确、能客观的反映自己的成果优劣进行成绩评定（三人相同） 4 对所抽取的试题能选择合理的程序计算，并根据完成程度进行成绩评定 总评成绩 总评成绩等级（优、良、中、及格、不及格） 指导教师签字： 算法理论 消元过程：消元结果冲掉，行乘数冲掉，det存放行列式。 （1）。 对于做(~⑤： (按列选主元，即确定使 。 (若，则输出错误信息；停机。 (若转④，否则执行换行 ④消元计算： ⑤ 。 若，则判断是否等于零： (，则输出“无解”信息；停机。 (，则输出“非唯一解”，停机。 若，则。 回带过程：解冲掉. 输出及。值的大小可以作为方程组是否病态的参考。 二、算法框图 三、算法程序 #include stdio.h #include math.h #define N 3 //未知函数个数 int I; float max_value(float a[N][N+1],int n,int k) { float max; int i; max=a[k][k]; for(i=k+1;in;i++) if(maxa[i][k]) { max=a[i][k]; I=i; } return(max); } void change(float *p,float *q) { float temp; temp=*p; *p=*q; *q=temp; return; } main() {float a[N][N+1]={{0.101,2.304,3.555,1.183},{-1.347,3.712,4.623,2.137},{-2.835,1.072,5.643,3.035}}; //方程组系数 float x[N],max,m,*c,*d; int k=0,q=0,n=N,i,j; for(k=0;kn q==0;k++) { max=max_value(a,n,k); if(max==0) q=1; else { if(I!=k) { for(j=k;jn+1;j++) { c=a[I][j]; d=a[k][j]; change(c,d); } } for(i=k+1;in;i++) { m=a[i][k]/a[k][k]; for(j=0;jn+1;j++) a[i][j]=a[i][j]-a[k][j]*m; } } } for(i=n-1;i=0;i--) { for(j=i+1;jn;j++) a[i][3]=a[i][3]-a[i][j]*x[j]; //3为未知函数个数 x[i]=a[i][3]/a[i][i]; //3为未知函数个数 } printf(得到的结果如下：\n); for(i=0;in;i++) printf( x[%d]=% 6.4f\n,i+1,x[i]); } 四、算法实现 例1．用高斯列主元素消去法解方程组的解。 解：把上页代码编译运行后，输入方程组系数即可求得如下结果 例2．用高斯列主元素消去法解方程组的解 解：在