最小二乘法(第4次试验).docVIP

PAGE  PAGE 8 实 验 报 告 实验课程名称 计算机数值方法 实验项目名称 曲线拟合的最小二乘法 年 级 专 业 学生姓名 学 号 理 学 院 实验时间： 2012 年 月 日 学生所在学院： 专业： 班级： 姓 名学 号实验组1实验时间指导教师成 绩实验项目名称曲线拟合的最小二乘法 实验目的及要求： 掌握数据拟合的思想，清楚数据拟合与插值法的区别及最小二乘原理在数据拟合中的重要作用及最小二乘解的求法。根据一组数据学会编制程序验证并求解拟合多项式。 实验（或算法）原理： 对给定数据 ，在取定的函数类中,求,使误差的平方和最小，即 从几何意义上讲，就是寻求与给定点 的距离平方和为最小的曲线函数称为拟合函数或最小二乘解，求拟合函数的方法称为曲线拟合的最小二乘法。 在这里主要讨论多项式拟合假设给定数据点 ，为所有次数不超过的多项式构成的函数类，现求一多项式使得 (1) 当拟合函数为多项式时，称为多项式拟合，满足式（1）的称为最小二乘拟合多项式。特别地，当n=1时，称为线性拟合或直线拟合。 显然 为的多元函数，因此上述问题即为求的极值问题。由多元函数求极值的必要条件，得 (2) 即  (3) （3）是关于的线性方程组，称为正规方程组或法方程组；方程组（3）存在唯一解。从式（3）中解出 (k=0,1,…，n)，从而可得多项式 （4）可以证明，式（4）中的满足式（1），即为所求的拟合多项式。 实验硬件及软件平台： PC机，vc++6.0，Internet网 实验步骤： 根据算法事先写出相应程序； 2．启动PC机，进入vc集成环境，输入代码； 3．编译调试； 4．调试通过，计算出正确结果。实验内容（包括实验具体内容、算法分析、源代码等等）： 算法设计。 (1)输入数据 （2）生成法方程，即计算系数矩阵G的元素 以及右端分量b （3）解法方程组 （4）输出 编写相应的程序上机调试。 #includemath.h #includestdio.h #define N 2 #define M 7 void gauss(double a[N][N],double b[N],double x[N]) {int i,j,k,l; double max,p,q,s; for(k=0;kN-1;k++) {max=fabs(a[k][k]); l=k; for(i=k+1;iN;i++) if(fabs(a[i][k])max) {max=fabs(a[i][k]); l=i; } if(l!=k) {for(j=k;jN;j++) {p=a[k][j]; a[k][j]=a[l][j]; a[l][j]=p; } p=b[k]; b[k]=b[l]; b[l]=p; } for(i=k+1;iN;i++) {q=a[i][k]/a[k][k]; for(j=k+1;jN;j++) a[i][j]=a[i][j]-q*a[k][j]; b[i]=b[i]-q*b[k]; } }  x[N-1]=b[N-1]/a[N-1][N-1]; for(i=N-2;i=0;i--) {s=0; for(j=i+1;jN;j++) s=s+a[i][j]*x[j]; x[i]=(b[i]-s)/a[i][i]; } } void polysimu(double x[M],double y[M],double w[M],double t[N]) {int i,j,k; double a[N][N],b[N]; for(i=0;iN;i++) for(j=0;jN;j++) a[i][j]=0; for(i=0;iN;i++) for(j=i;jN;j++) {for(k=0;kM;k++)