插值法与最小二乘法.docVIP

插值法与最小二乘法 实 验 报 告 实验名称 （教师填写） 插值法与最小二乘法 实验目的 （教师填写） 掌握插值法与最小二乘法。 实验题目 （教师填写） 完成以下两题： 试利用在100，121，144处的函数值求的近似值并估计误差. 给定数据点如下： -2 -1 0 1 2 3 1 2 2 3 3 5 分别以一次、二次、三次多项式曲线拟合以上数据. 实验报告要求 （教师填写） 完成编程与计算后写出实验报告，主要内容包括：所用算法语言，算法步骤描述，重要变量说明，程序清单，输出计算结果，结果分析。 实验内容与步骤 （学生填写） 所用算法语言：Visual C++ 算法步骤描述：求得的二次Lagrange插值多项式 Step1：Step2：解出Lagrange插值多项式； Step3：输出结果。 重要变量说明：数组x[3]为函数自变量x；y[3]为函数值；l[3]为Lagrange插值多项式的基函数组；Sum为最终由Lagrange插值法求得的近似值；R为近似误差；i、j为控制变量。 程序清单：#includestdio.h void main() { double x[3]={100,121,144}; double y[3]={10,11,12}; double l[3]={1,1,1}; double Sum=0; double R; int i,j; double X=115; for(i=0;i3;i++) { for(j=0;j3;j++) if(i!=j) { l[i]*=(X-x[j])/(x[i]-x[j]); continue; } } for(i=0;i3;i++) Sum+=y[i]*l[i]; printf(Sum=%lf\n,Sum); R=1; for(i=0;i3;i++) R*=X-x[i]; R*=1e-5/16; printf(R=%le\n,R); } 实验内容与步骤 （学生填写） 算法步骤描述： Step1：将所给数据进行预处理，分别列出列出可获得最小二乘解的正则方程组： Step2：用Gauss消去法分别解三个方程，得出相应参数； Step3：输出结果。 重要变量说明：数组x[6]、y[6]即题目所给的原始数据；a[2][3]、b[3][4]、c[4][5]分别为一次、二次、三次多项式的x矩阵，用来求解参数；m[6]中存放的为x的i次方和；n[4]中存放的是y与x的i次方乘积和；a1[2]、a2[3]、a3[4]分别为一次、二次、三次多项式的系数；i、j为控制变量；k为数据的个数6。 程序清单：#includestdio.h #includemath.h void main() { double x[6]={-2,-1,0,1,2,3}; double y[6]={1,2,2,3,3,5}; double a[2][3],b[3][4],c[4][5]; double m[6]={0},n[4]={0}; double a1[2],a2[3],a3[4]; int i,j,k=6; for(j=0;j6;j++) for(i=0;i6;i++) m[j]+=pow(x[i],(j+1)); for(j=0;j4;j++) for(i=0;i6;i++) n[j]+=y[i]*pow(x[i],j); a[0][0]=k; a[0][1]=a[1][0]=m[0]; a[1][1]=m[1]; a[0][2]=n[0]; a[1][2]=n[1]; for(i=0;i3;i++) a[1][i]+=-a[0][i]/k*m[0]; a1[1]=a[1][2]/a[1][1]; a1[0]=(a[0][2]-a1[1]*a[0][1])/a[0][0]; printf(The first order:\n); printf(%f\t%f\n,a1[0],a1[1]); 实验内容与步骤 （学生填写） b[0][0]=k; b[0][1]=b[1][0]=m[0]; b[0][2]=b[1][1]=b[2][0]=m[1]; b[1][2]=b[2][1]=m[2]; b[2][2]=m[3]; b[0][3]=n[0]; b[1][3]=n[1]; b[2][3]=n[2]; for(i=0;i4;i++) { b[1][i]+=-b[0][i]/k*m[0]; b[2][i]+=-b[0][i]/k*m[1]