线性方程组的直接解法实验报告.doc

本科实验报告课程名称： 数值计算方法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(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*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循环来输入各个元素的值。还有就是将方程组化为上三角矩阵，也要用到好几个for循环，比较容易出错。七、讨论、心得做实验要求我们把基础学扎实，上机实验让我又重新巩固了C语言知识。我的实验是线性方程组的直接解法，刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。最小二乘拟合多项式 一、实验内容给定数据点（xi ,yi），用最小二乘法拟合数据的多项式，并求平方误差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713.002、 计算公式 y=a0+a1x3、 结构程序设计/最小二乘法拟合/本程序包含两个函数，主函数与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=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(%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(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();4、 实验结果与讨论、心得讨论、心得: 刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这