《线性方程组的直接解法实验报告》由会员分享，可在线阅读，更多相关《线性方程组的直接解法实验报告（8页珍藏版）》请在金锄头文库上搜索。
1、 本科实验报告课程名称：数值计算方法B实验项目：线性方程组的直接解法最小二乘拟合多项式实验地点：ZSA401专业班级：学号：201000学生某：指导教师：李志 2012年4月13日线性方程组的直接解法一、实验目的和要求实验目的：合理利用Gauss消元法、LU分解法或追赶法求解方程组。实验要求：利用高斯消元法，LU分解法或追赶法进行编程，求解题中所给的方程组。二、实验内容和原理实验内容：合理利用Gauss消元法、LU分解法或追赶法求解下列方程组：（n=5,10,100,）实验原理：这个实验我选用的是高斯消元法。高斯消元法：先按照Lik=aik(k-1)/akk(k-1) ,aij(k)=aij(
2、k-1)-likakj(k-1)其中k=1,2,n-1;i=k+1,k+2,n;j=k+1,k+2,n+1 将方程组变为上三角矩阵，再经过回代，即可求解出方程组的解。三.计算公式 通过消元、再回代的求解方法称为高斯消元法。特点是始终消去主对角线 下方的元素。四、操作方法与实验步骤#include Stdio.h#define N 3main() double aNN+1,bN; int i,j,k,x=0; for(i=0;iN;i+) for(j=0;j=N;j+) scanf(%lf,&aij); for(k=0;kN;k+) for(i=x;i=0;j-) ai+1j=ai+1j-axj
3、*ai+1x/axx; x+; for(i=0;iN;i+) for(j=0;j=N;j+) printf(%lf ,aij); printf(n); b2=a23/a22; b1=(a13-a12*b2)/a11; b0=(a03-a02*b2-a01*b1)/a00; for(i=0;iN;i+) printf(x(%d)=%lfn,i+1,bi); getch(); 五、实验数据记录和处理实验结果：六、实验结果与分析大体来说，实现了课程设计的算法要求及功能，有很多还不能很好的处理的问题，需要我们在改进中不断完善。做本次实验中有些地方有困难，就是获得数组中各元素的值，要用到for循环来输入
4、各个元素的值。还有就是将方程组化为上三角矩阵，也要用到好几个for循环，比较容易出错。七、讨论、心得做实验要求我们把基础学扎实，上机实验让我又重新巩固了C语言知识。我的实验是线性方程组的直接解法，刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。最小二乘拟合多项式一、实验内容给定数据点（xi ,yi），用最小二乘法拟合数据的多项式，并求平方误差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713.00二、 计算公式y=a0+a1x三、 结构程序设计/最
5、小二乘法拟合/本程序包含两个函数，主函数与Gauss消元函数#include#include#define N 10#define M 20double ANN+1;void Gauss_eliminate(int n,double* answer) int k,i,j,sum;/*消元重构矩阵使之成为上三角矩阵*/ for(k=1;kn;k+) /k为高斯消元法上标，指示消元代数 for(i=k+1;i=n;i+) Aik=Aik/Akk; for(i=k+1;i=n;i+) for(j=k+1;j=n+1;j+) Aij=Aij-Akj*Aik; for(i=2;i=n;i+) for(j
6、=1;j=1;i-) sum=0; for(j=i+1;j=n;j+) sum+=Aij*answerj; answeri=(Ain+1-sum)/Aii; void main() int k; printf(请输入拟合次数：n); scanf(%d,&k); int n=k+1; int m,i,j; double sumM; double xM,yM; printf(请输入散点的个数:); scanf(%d,&m); printf(请创建数组x:n); for(i=0;im;i+) scanf(%lf,&xi); printf(请创建数组y:n); for(i=0;im;i+) scanf
7、(%lf,&yi); /求系数矩阵 for(i=0;iM;i+) sumi=0; /初始化 for(i=0;i2*k+1;i+) for(j=0;jm;j+) sumi+=pow(xj,i); for(i=0;in;i+) for(j=0;jm;j+) sum2*k+1+i+=pow(xj,i)*yj; /创建系数矩阵 int count; for(i=1;i=n;i+) /先建上三角矩阵 Aii=sum2*(i-1); for(j=i+1,count=1;j=n;j+) Aij=sum2*(i-1)+count; count+; for(i=2;i=n;i+) /通过对称建立下三角 for(
8、j=1;ji;j+) Aij=Aji; for(i=1;i=n;i+) /增广矩阵的一列 Ain+1=sum2*k+i; /输出测试 printf(增广矩阵为:n); for(i=1;i=n;i+) for(j=1;j=n+1;j+) printf(%lft,Aij); printf(n); double answerN; Gauss_eliminate(n,answer); /*输出解*/ for(i=1;i=n;i+) printf(a%d=%lft,i-1,answeri); getchar(); getchar();四、 实验结果与讨论、心得讨论、心得:刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。这段时间的实验课提高了我的分析问题，解决问题的能力，特别提高了对一个程序的整体操作能力，对程序的细微之处有了明显的提高认识的地方，追求最实用的程序，弥补学习上的不足，同时认识到还应深入理解课本上的知识，学过的东西要知道理论与实践相结合，增加动手能力。 /
