n维线性方程组的求解程序.doc

C语言程序设计课程设计总结报告设计题目：n维线性方程组的求解程序。姓 名：朱战威学 号：111411249同组人员：谢忠建 闫俊宇 卿帅军 尚亚军指导教师：张芳芳、王康平、薛冰 2012年 6 月 21 日一.设计题目生物工程系C语言程序设计 n维线性方程组的求解二.设计时间及地点设计时间：6月21日 设计地点：3#A506三、 题目内容要求能正确的求解n维线性方程的答案 四、 设计部分：1.掌握C语言程序设计的基本思想.概念.方法和技能；2.初步学会利用 c 语言编制一些程序以解决一般实际应用问题；3.培养学生分析问题和解决问题的能力； 五.设计部分设计程序的主要步骤:1. 了解n维线性方程的求解步骤,可根据高斯消元法进行运算。2. 按消元规则进行运算后，再进行元素回代3. 矩阵在c语言中可用二维数组表示4. 循环使用for语句表示此程序运用高斯消元法：（1）按消元规则进行运算后，对角线以下元素为0，故运算中对角线以下的元素不作计算。（2）对角线下的元素回代无影响，将乘数放在该处。既 aik=aik/akk，i=k+1,k+2,.,n;（3）对角线以上的元素和常数项采用“原地”工作方式，既经变换后的元素仍放在原来的位置上 aij=aij-aik*akj;bi=bi-aik*bk;i,j=k+1,k+2,.,n; 目的是为了节省存储单元（4）回代后的值仍放在常数项存储单元bn=bn/ann;bi-=aij*bj;bi=(1/aii)*bi; 这样，数组b中单元存放的就是未知数的解。*#includeusing namespace std;void main()int i=0,j=0,k,n,sum=0;coutn; /输入线性方程未知数的个数double *a=new double*n+1;/分配存储系数矩阵的二维数组（动态分配）for(i=0;i=n;i+)ai=new doublen+1;double *b=new doublen+1; /分配存储方程组右侧常数项的一维数组cout输入线性方程的系数矩的每一行和等号右侧结果:endl;for(i=1;i=n;i+)for(j=1;jaij; /输入系数矩阵的元素，没有的项输入0cinbi;/输入方程组每个等式的等号右边的常数项k=1;while(1)for(i=k+1;i=n;i+)aik=aik/akk;/对角线下的元素回代无影响，将乘数放在该处。for(i=k+1;i=n;i+)for(j=k+1;j=1;i-)for(j=i+1;j=n;j+)bi-=aij*bj;bi=(1/aii)*bi; /回代后的值仍放在常数项存储单元sum=0;for(i=1;i=n;i+)coutxi = biendl;for(i=0;i=n;i+)delete ai;delete a;delete b;六.程序代码#includeusing namespace std;void main()int i=0,j=0,k,n,sum=0;coutn;double *a=new double*n+1;for(i=0;i=n;i+)ai=new doublen+1;double *b=new doublen+1;cout输入线性方程的系数矩的每一行和等号右侧结果:endl;for(i=1;i=n;i+)for(j=1;jaij;cinbi;k=1;while(1)for(i=k+1;i=n;i+)aik=aik/akk;for(i=k+1;i=n;i+)for(j=k+1;j=1;i-)for(j=i+1;j=n;j+)bi-=aij*bj;bi=(1/aii)*bi;sum=0;for(i=1;i=n;i+)coutxi = biendl;for(i=0;i=n;i+)delete ai;delete a;delete b;如果你要计算方程组：2x1-x2+3x3=14x1+2x2+5x3=4x1+2x2+0x3=7七心得体会 通过本次课程设计我对书本上的理论知识加以运用，了解C语言课程设计的真正用意培养自学能力，养成程序编辑的好习惯。从拿到题目到完成整个编程，从理论到实践学到很多很多的的东西，同时不仅可以巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务。本次设计也遇到了各种各样的问题，同时在设计的过程中发现了自己的不足之处，对以前所学过的知识理解得不够深刻，掌握得不够牢固，这次课程设计之后，一定把以