拉格朗日实验报告.docVIP

09信息 郭凯 39号 拉格朗日插值实验报告 #includestdio.h int main() { int i,n,j; float a[50],b[50],x,l,y; printf(enter n:); scanf(%d,n); for(i=0;i=n;i++) { printf(x[%d] y[%d]:,i,i); scanf(%f%f,a[i],b[i]); } printf(enter x:); scanf(%f,x); y=0; for(i=0;i=n;i++) { l=1; for(j=0;j=n;j++) { if(j!=i) l=l*(x-a[j])/(a[i]-a[j]); } y=y+l*b[i]; } printf(y=%f\n,y); return 0; } 例：给出函数f(x)=sinx的3组对应的函数值如下表： x 0.32 0.34 0.36 f(x) 0.314567 0.333487 0.352274 利用线性插值及抛物插值计算sin0.3367的值并估计截断误差。 解：（1.1）利用以上程序，由于是利用线性插值，则首先输入n的值为1，由于0.3367介于0.32与0.34之间，故依次输入2对x与相对应的f(x)的值，再输入x的值即0.3367，即可得如下测试结果： 即计算所得的sin0.3367的值为0.330365 误差分析： 余项R1(x)=f(x)-L1(x)=f(2)(q)(x-x0)(x-x1) /2! |R2(x)|≦M2|(0.3367-0.32)(0.3367-0.34) |/2, 其中 M2=max0.32≦x≦0.34| f(2)（x）|=sin0.340.3335 于是 |R2(0.3367)|=|sin0.3367- L1(0.3367)| ≦0.3335*0.0167*0.033/20.92*10-5 （1.2）由于题目中给出了3个x的值，故可选用比0.3367大的两个x的值来计算，测试结果如下： 误差分析同上类似，计算得|R2(0.3367)| 0.1266*10-4 （2）利用以上程序，由于是利用抛物插值，则首先输入n的值为2，然后依次输入3对x与相对应的f(x)的值，再输入x的值即0.3367，即可得如下测试结果： 即计算所得的sin0.3367的值为0.330374 误差分析： 余项R2(x)=f(x)-L2(x)=f(3)(q)(x-x0)(x-x1)(x-x2)/3! | R2(x)|≦M3|(0.3367-0.32)(0.3367-0.34)(0.3367-0.36)|/6, 其中 M3=max0.32≦x≦0.36| f(3)（x）|=cos0.320.9493, 于是 | R2(0.3367)|=|sin0.3367- L2(0.3367)| ≦0.9493*0.0167*0.033*0.0233/62.0132*10-6 实验心得： 由上述例子可以看出，(1.1)与（1.2）比较时，（1.1）的结果优于（1.2），即一般情况下，内插优于外推。而（1.1）与（2）比较时，发现利用2次插值比利用1次插值精度高，误差小，即一般情况下，次数高的插值优于次数低的插值。 牛顿插值实验报告 #includestdio.h int main() { int i,n,j,k; float a[100],b[100],x,l,p; printf(enter n:); scanf(%d,n); for(i=0;i=n;i++) { printf(x[%d] y[%d]:,i,i); scanf(%f%f,a[i],b[i]); } printf(enter x:); scanf(%f,x); k=1; for(i=n;i=1;i--) { for(j=n;jn-i;j--) { b[j]=(b[j]-b[j-1])/(a[j]-a[j-k]); } k++; } for(i=0;i=n;i++) { printf(y[%d]=%f\n,i,b[i]); } p=b[n]; for(j=n;j=1;j--) { p=p*(x-a[j-1])+b[j-1]; } printf(p=%f\n,p); return 0; } 例：给出f（x）的函数表如下，利用4次牛顿插值，计算f（0.596）的近似值。 x 0.40 0.55 0.65 0.80 0.90