数值分析实验五new

《数值分析实验五new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、牛顿插值法一、实验目的：学会牛顿插值法，并应用算法于实际问题。二、实验内容：给定函数，已知：三、实验要求：（1）用牛顿插值法求4次Newton插值多项式在2.15处的值，以此作为函数的近似值。在MATLAB中用内部函数ezplot绘制出4次Newton插值多项式的函数图形。（2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次进行比较。四、算法描述插值与逼近都是用某个简单函数在满足一定条件下，在某个范围内近似代替真实函数的一部分，以此来达到简化计算的的目的。牛顿多项式插值法就是其中一种利用这种近似逼近思想的计算方法。五、实验步骤（1）用课本

2、上的先列差商表，再根据差商表得出牛顿差值多项式，再代入x=2.15解得f(2.15)。然后再用MATLAB得出牛顿差值多项式。用MATLAB中内部函数ezplot绘制出4次Newton插值多项式的函数图形。（2）在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形，并与作出的4次Newton插值多项式的图形进行比较。六、实验结果（1）(①)解：差商表xF(x)一阶差商二阶差商三阶差商四阶差商X0=21.414214X1=2.11.4491380.34924X2=2.21.4832400.34102-0.0411X3=2.31.5165750.33335-0.0

3、38350.0084X4=2.41.5491930.32618-0.035850.00830.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=[22.12.22.32.4];Y=[1.4142141.4491381.4832401.5165751.549193];[A,C,L,wcgs,Cw]=newploy(X,Y)A=1.414200001.44910.34920001.48320.3

4、410-0.0411001.51660.3334-0.03830.009201.54920.3262-0.03590.0083-0.0021C=-0.00210.0271-0.15660.72150.4141L=-4803839603609061/2305843009213693952*x^4+7806239355294329/288230376151711744*x^3-176292469178709/1125899906842624*x^2+1624739243112817/2251799813685248*x+1865116246031207/45035996273

5、70496Wcgs=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.09170.4029-0.88460.9700-0.4250即L=-0.0021+0.0271-0.1566+0.7215x+0.4141(②)则N(2.15)=1.46628861用MATLAB中内部函数ezplot绘制出4次Newton插值多项式的函数图形（a）和用plot作出的4次Newton插值多项式的图形(b)的程序为>>y=-0.2169.*x.

6、^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)输入n个节点七、程序流程图按一定的原则取两个点得到i阶差商I=n??NY构造n阶差商表由选定的基和得到的差商表得到牛顿插值函数，化简八、实验结果分析（1）=1.46628861，在x=2.1和x=2.2之间，误差也很小，可见牛顿插值函数的现实可行性。（2）MATL

7、AB中用内部函数ezplot绘制出4次Newton插值多项式的函数图形和在MATLAB中用内部函数ezplot可直接绘制出以上函数的图形只有很小的差别，二者在[1,3]这个区间上整体相似，仅仅很小的地方不同，如函数的最顶端部分和个别的地方的斜率等。因为是同一个函数，所以其图形差别不大，但是隐函数是真实的图形，而牛顿插值法得到的函数是去真实函数上的一些点取的相似函数，与真实函数之间还有一些误差的。