数值分析实验五.docVIP

牛顿插值法 一、实验目的： 学会牛顿插值法，并应用算法于实际问题。 实验内容： 给定函数 ，已知： 三、实验要求： （1）用牛顿插值法求4次Newton插值多项式在2.15处的值，以此作为函数的近似值。在MATLAB中用内部函数ezplot绘制出4次Newton插值多项式的函数图形。 （2） 在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次 进行比较。 四、算法描述 插值与逼近都是用某个简单函数在满足一定条件下，在某个范围内近似代替真实函数的一部分，以此来达到简化计算的的目的。牛顿多项式插值法就是其中一种利用这种近似逼近思想的计算方法。 五、实验步骤 （1）用课本上的先列差商表，再根据差商表得出牛顿差值多项式，再代入x=2.15解得f(2.15)。然后再用MATLAB得出牛顿差值多项式。用MATLAB中内部函数ezplot绘制出4次Newton插值多项式的函数图形。 （2） 在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。 六、实验结果 （1）(①)解：差商表 x F(x) 一阶差商 二阶差商 三阶差商 四阶差商 X0=2 1.414214 X1=2.1 1.449138 0.34924 X2=2.2 1.483240 0.34102 -0.0411 X3=2.3 1.516575 0.33335 -0.03835 0.0084 X4=2.4 1.549193 0.32618 -0.03585 0.0083 0.00025 所以，=f(2)+0.34924(x-2)-0.0411(x-2)(x-2.1)+0.0084(x-2)(x-2.1)(x-2.2)+0.00025(x-2)(x-2.1)(x-2.2)(x-2.3) MATLAB算法为： X=[2 2.1 2.2 2.3 2.4];Y=[1.414214 1.449138 1.483240 1.516575 1.549193];[A,C,L,wcgs,Cw]= newploy (X,Y) A =1.4142 0 0 0 0 1.4491 0.3492 0 0 0 1.4832 0.3410 -0.0411 0 0 1.5166 0.3334 -0.0383 0.0092 0 1.5492 0.3262 -0.0359 0.0083 -0.0021 C =-0.0021 0.0271 -0.1566 0.7215 0.4141 L=-4803839603609061/2305843009213693952*x^4+7806239355294329/288230376151711744*x^3-176292469178709/1125899906842624*x^2+1624739243112817/2251799813685248*x+1865116246031207/4503599627370496 Wcgs=1/120*M*(x^5-11*x^4+967/20*x^3-2123/20*x^2+2047772676816529/17592186044416*x-31878/625) Cw =0.0083 -0.0917 0.4029 -0.8846 0.9700 -0.4250即L=-0.0021+0.0271-0.1566+0.7215x+0.4141 (②)则N(2.15)=1用MATLAB中内部函数ezplot绘制出4次Newton插值多项式的函数图形（a）和用plot作出的4次Newton插值多项式的图形(b)的程序为 y=-0.2169.*x.^4+0.0648.*x.^3+2.1076.*x.^2+3.3960.*x-4.5745; plot(x,y); figure(3); ezplot(@(x)[-0.2169.*x.^4+0.0648.*x.^3+2.1076.*x.^2+3.3960.*x-4.5745],[-2,2]); 两个图形为 图 (a) 图(b) 七、程序流程图