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

《最小二乘拟合多项式实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、本科实验报告课程名称：数值计算方法B实验项目：最小二乘拟合多项式实验地点：ZSA401专业班级：软件10班学号：2学生姓名：指导教师：李志2012年4月9日一．方程求根一、实验内容给定数据点（xi,yi），用最小二乘法拟合数据的多项式，并求平方误差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713.00二、计算公式y=a0+a1x三、结构程序设计#include"stdafx.h"#include"stdio.h"#include"math.h"#definenum10floatneiji(floatb[num],floatc[

2、num])/*内积函数*/{intp;floatnj=0;for(p=1;p

3、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+

4、+){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",index+1);for(i=0;i<=index;i++){for(j=i;j<=index+i;j++)fprintf(f,"%f",sumX[j]

5、);fprintf(f,"%f",sumY[i+1]);}fclose(f);//利用求解线性方程组的方法求解…}二、实验结果与讨论、心得讨论、心得:刚开始调试代码的时候有时候就是很小的错误导致整个程序不能运行，需要我们一步一步慢慢来，经过无数次的检查程序错误的原因，以及在老师的帮助下，完成了这次实验。这段时间的实验课提高了我的分析问题，解决问题的能力，特别提高了对一个程序的整体操作能力，对程序的细微之处有了明显的提高认识的地方，追求最实用的程序，弥补学习上的不足，同时认识到还应深入理解课本上的知识，学过的东西要知道理论与实践相结合，增加动手能力。