最小二乘拟合多项式实验报告.docVIP

本科实验报告 课程名称： 数值计算方法B 实验项目： 最小二乘拟合多项式 实验地点： ZSA401 专业班级： 软件10班 学号： 2 学生姓名： 指导教师： 李志 2012年 4月 9日 一．方程求根 一、实验内容 给定数据点（xi ,yi），用最小二乘法拟合数据的多项式，并求平方误差。 xi 0 0.5 0.6 0.7 0.8 0.9 1.0 yi 1 1.75 1.96 2.19 2.44 2.71 3.00 计算公式 y=a0+a1x 结构程序设计 #include stdafx.h #include stdio.h #include math.h #define num 10 float neiji(float b[num],float c[num]) /*内积函数*/ { int p; float nj=0; for (p=1;pnum;p++) nj+=c[p]*b[p]; return nj; } float power(float a,int n) { float b=1; for(int i=0;in;i++) { b*=a; } return b; } float x[num],y[num],sumX[num],sumY[num]; void main() { int i,j,n,index; char conti; FILE *f; conti= ; printf(请输入已知点的个数n=\n); scanf(%d,n); printf(请输入x和y:); sumX[1]=sumY[1]=0; for(i=1;i=n;i++) { printf(x[%d]=,i); scanf(%f,x[i]); sumX[1]+=x[i]; printf(y[%d]=,i); scanf(%f,y[i]); sumY[1]+=y[i]; } printf(sumX[1]=%f,sumX[1]); printf(sumY[1]=%f,sumY[1]); printf(请输入拟和次数=); scanf(%d,index); i = n; sumX[0]=i; for(i=2;i=2*index;i++) { for(j=1;j=n;j++) { sumX[i]+=power(x[j],i); } printf(sumX[%d]=%f,i,sumX[i]); } for(i=2;i=index+1;i++) { for(j=1;j=n;j++) { sumY[i]+=power(x[j],i-1)*y[j]; } printf(sumY[%d]=%f,i,sumY[i]); } f=fopen(Gaussdata.txt,w); fprintf(f,%i\n,index+1); for(i=0;i=index;i++) { for(j=i;j=index+i;j++) fprintf(f,%f ,sumX[j]); fprintf(f,%f\n,sumY[i+1]); } fclose(f); //利用求解线性方程组的方法求解… } 实验结果与讨论、心得 讨论、心得: 刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。 这段时间的实验课提高了我的分析问题，解决问题的能力，特别提高了对一个程序的整体操作能力，对程序的细微之处有了明显的提高认识的地方，追求最实用的程序，弥补学习上的不足，同时认识到还应深入理解课本上的知识，学过的东西要知道理论与实践相结合，增加动手能力。