计算方法考查论文重点.docVIP

盐城师范学院数学统计学院 《数值分析》考查论文 第一组 班级 师范132班 学号 1座机电话号码 姓名 卞秀丽 内容：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,DisplayFunction Identity,PlotSt yle- PointSize[0.02]]; fd Plot[f[x], x,10,16 ,DisplayFunction Indentity]; 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[,,] -0.0035 F[,] 0.0742 F[,,] -0.0029 F[,] 0.0689 F[,,] -0.00265 F[,] 0.0646 F[,,] -0.00215 F[,,] 0.000座机电话号码 F[,,] 0.0002 F[,,]