后小霞2012214023.docVIP

《数值分析》上机实验报告 ———— 函数插值方法 现代农业工程学院 农业水土工程 2012214023 后小霞 任课教师：罗志强 课题名称 课题五 函数插值法 一、 问题提出 对于给定的一元函数 的n+1个节点值 。试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。 数据如下： （1） 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式L，和分段三次插值多项式，计算 的值。 （2） 1 2 3 4 5 6 7 0.368 0.135 0.050 0.018 0.007 0.002 0.001 试构造Lagrange多项式L，计算的值。 结果0.165299 0 要求 1、利用Lagrange插值公式 编写出插值多项式程序； 2、给出插值多项式或分段三次插值多项式的表达式； 3、根据节点选取原则，对问题（2）用三点插值或二点插值，其结果如何； 4、对此插值问题用Newton插值多项式其结果如何。 三、班级、姓名、学号 现代农业工程学院 农业水土工程 后小霞 2012214023 四、目的与意义 1、学会常用的插值方法，求函数的近似表达式，以解决其它实际问题； 2、明确插值多项式和分段插值多项式各自的优缺点； 3、熟悉插值方法的程序编制； 五、 n次的拉格朗日插值多项式为，其中是插值基函数，且 插值基函数计算公式为： 2. Newton插值法 Newton均差插值多项式为： 六、结构程序设计function [lk,L]=lagrange(X,Y) a=length(X); L=ones(a,a); for n=1: a; V=1; for i=1:a if n~=i V=conv(V,poly(X(i)))/(X(k)-X(i)); end end L1(k,:)=V; l(k,:)=poly2sym (V); end lk=Y*L1;L=Y*l; 其中L为拉格朗日插值多项式, lk为基函数。. (2)求表1中五次Lagrange多项式L 在窗口中输入： X=[0.4 0.55 0.65 0.80 0.95 1.05]; Y=[0.41075 0.57815 0.69675 0.90 1.00 1.25382]; [lk,L]=lagrange(X,Y) 运行结果C和L： lk = 121.6264 -422.7503 572.5667 -377.2549 121.9718 -15.0845 L = (3913328*x^5)/32175-(116018793466243525133*x^4)/274438102292889600 + (654725549453611782859*x^3)/1143492092887040000- (142358028293743156047463*x^2)/377352390652723200000+ (276158189076188508336341*x)/2264114343916339200000 -3750000 取L＝121.6264-422.7503+572.5667-377.2549+121.9718-15.0845 计算表1 lagranger插值多项式对应插值点的值 x=0.596; M=1; X=[0.4,0.55,0.65,0.80,0.95,1.05]; Y=[0.41075,0.57815,0.69675,0.90,1.00,1.25382]; [y,R]=lagranger(X,Y,x,M) 运行结果： y = 0.6257 R = 2.2170e-008 在命令窗口输入程序： x=0.99; M=1; X=[0.4,0.55,0.65,0.80,0.95,1.05]; Y=[0.41075,0.57815,0.69675,0.90,1.00,1.25382]; [y,R]=lagranger(X,Y,x,M) 运行结果： y = 1.0542 R = 5.5901e-008 （）求表二六次lagrange多项式 M=[1:1:7]; Y=[0.368 0.135 0.050 0.018 0.007 0.002 0.001]; [lk, L]=lagrang