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

本科实验报告课程名称： 数值计算方法B 实验项目： 最小二乘拟合多项式 实验地点： ZSA401 专业班级： 软件10班 学号： 2 学生姓名： 指导教师： 李志 2012年 4月 9日一方程求根 一、实验内容给定数据点（xi ,yi），用最小二乘法拟合数据的多项式，并求平方误差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713.002、 计算公式 y=a0+a1x3、 结构程序设计#include stdafx.h#include stdio.h#include math.h#define num 10float neiji(float bnum,float cnum) /*内积函数*/ int p;float nj=0;for (p=1;pnum;p+)nj+=cp*bp;return nj;float power(float &a,int n)float b=1;for(int i=0;in;i+)b*=a;return b;float xnum,ynum,sumXnum,sumYnum;void main() int i,j,n,index;char conti;FILE *f;conti= ;printf(请输入已知点的个数n=n);scanf(%d,&n);printf(请输入x和y:);sumX1=sumY1=0;for(i=1;i=n;i+) printf(x%d=,i);scanf(%f,&xi);sumX1+=xi;printf(y%d=,i);scanf(%f,&yi);sumY1+=yi;printf(sumX1=%f,sumX1);printf(sumY1=%f,sumY1);printf(请输入拟和次数=);scanf(%d,&index);i = n;sumX0=i;for(i=2;i=2*index;i+)for(j=1;j=n;j+)sumXi+=power(xj,i);printf(sumX%d=%f,i,sumXi);for(i=2;i=index+1;i+)for(j=1;j=n;j+)sumYi+=power(xj,i-1)*yj;printf(sumY%d=%f,i,sumYi);f=fopen(Gaussdata.txt,w);fprintf(f,%in,index+1);for(i=0;i=index;i+)for(j=i;j=index+i;j+)fprintf(f,%f ,sumXj);fprintf(f,%fn,sumYi+1);fclose(f);/利用求解线性方程组的方法求解4、 实验结果与讨论、心得 讨论、心得:刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。这段时间的实验课提高了我的分析问题，解决问题的能力，特别提高了对一个程序的整体操作能力，对程序的细微之处有了明显的提高