计算方法考查论文可研报告.docVIP

盐城师范学院数学统计学院 《数值分析》考查论文(第一组) 班级 师范132班 学号姓名 卞秀丽 内容：1. 已知函数表如下： (40’) x 10 11 12 13 14 15 16 f(x) 2.3026 2.3979 2.4849 2.5649 2.6391 2.7080 2.7726 用下列方法计算f(10.05), f(15.04), f(12.94)，并讨论误差分析. 分段线性插值； 牛顿插值多项式； 全区间上的拉格朗日插值。 实验目的 本实验的目的是熟练《数值分析》第二章“差值逼近”的相关内容，掌握分段插值法、拉格朗日差值法和牛顿插值法。 二、相关背景知识介绍 （1）分段线性插值指的是：在n+1个差值节点a=…=b处，给定函数的值，寻求一函数P（x），使得在每个子区间[,]上，P（x）∈。P（x）=，i=0,1，…n.我们称P（x）为f（x）关于节点，，…,的分段线性插值多项式。 （2）（x）)+f[,](x-)+ f[,,](x-)(x-)+…+ f[,,…,](x-)(x-)…(x-)为f（x）关于节点，，…,的n次牛顿插值多项式。 （3）（x）=l（x），，…,的n次拉格朗日插值多项式。 三、程序 (1)分段插值法 In[1]:=Data={{10,2.3026},{11,2.3979},{12,2.4849},{13, 2.5649},{14,2.6391},{15,2.7080},{16,2.7726}}; F=Interpolation[date] Out[2]=InterpolatingFunction[{{10,16}},] In[3]:=pd=ListPlot[data,DisplayFunctionIdentity,PlotSt yle-PointSize[0.02]]; fd=Plot[f[x],{x,10,16},DisplayFunctionIndentity]; Show[{pd,fd},DispayFunction-$DisplayFunction] Out[5]=-Graphics- In[6]=f[10.05] Out[6]=2.307365 In[7]=f[15.04] Out[7]=2.710584 In[8]=f[12.94] Out[8]=2.5601 (2)牛顿插值法 clear?all clc? syms?x? =[10,11,12,13,14,15,16];? =[2.3026,2.3979,24849,2.5649,2.6391,2.7080,2.7726];? for?k=1:7?? ?for?i=1:k??? ?a=1;???? ?b=0;???? ?for?j=1:k?????? ?if?j~=i? ?a=a*((i)，(j));?????? ?end????? end? b=b+(i)/a;??? end??? A(k)=b;? end? B=[1,(x-(1)),(x-(1))*(x-(2)),(x-(1))*(x- (2)) *(x-(3)),(x-(1))*(x-(2))*(x-(3))*(x- (4))];? L1=A.*B;? l=0;? for?m=1:7?? ?l=l+L1(m);? end? L=expand(l) （3）拉格朗日插值法 syms?x? =[10,11,12,13,14,15,16];? =[2.3026,2.3979,24849,2.5649,2.6391,2.7080,2.7726];? for?i=1:7?? ?a=1;? ??for?j=1:5?? ???if?j~=i? a=expand(a*(x- (j))); ?end?? ?end??? b=1;??? for?k=1:7????? if?k~=i? b=b*(x0(i)-x0(k));????? end??? end? A(i)=expand(a/b);? end? L=0;? for?p=1:7? L=L+y0(p)*A(p);? end? L 四、数值结果 （1）分段插值法 x∈[10,11]时，=0.0953x+1.3496 则f[10.05]=L（10.05）=2.307365 x∈[15,16]时，=0.0646x+1.739 则f[15.04]=2.710584 x∈[12,13]时，=0.08x+1.5249 则f[12.94]=2.5601 (2)牛顿插值法 =10，=11，=12，=13，=14，=15，=16 F[,]=0.0953 F[,]=0.087 F[,,]=-0.00415 F[,]=0.08 F[,,]=-