多项式插值的振荡现象

1、ZHEHANG NORMAL UNIVERSITY数值分析实验报告多项式插值的振荡现象姓 名：学 院：数理与信息工程学院班 级：学 号：数值分析实验报告实验名称多项式插值的振荡现象实验时间2013年10月23日姓名班级学号成绩实验目的1 理解多项式插值，懂得它的振荡现象2. 研究样条插值，并分析它的收敛性。3. 学会在实际生活中使用二维插值。实验内容1. 设区间-1,1上函数f(x)11 25x2考虑区间-1,1的一个等距划分，分点为Xi2i .1 -i 0,1,2, ,n n则拉格朗日插值多项式为Ln(X)125x2li(x)其中的li（x）,i0,1,2, ,n是n次拉格朗日插值基函数2.

2、 请按一定的规则分别选择等距或者非等距的插值节点，并不断增加插值节点的个数。考虑实验1中的函数或选择其他你有兴趣的函数，可以用MATL的函数“ spline ”作此函数的三次样条插值。3. 在一丘陵地带测量高程，x和y方向每隔100米测一个点，得高程数据如下。试用 MATLAB的二维插值函数“ in terp2 ”进行插值，并由此找出最高点 和该点的高程。三、算法描述（1） 编写好拉格朗日插值函数，保存在M文件中；（2） 考虑到：1、一幅图中太多的曲线会相互覆盖；2、n取奇偶数可能结果不同；3、不同的节点选取方法可能导致不同的结果。故而n的选择分为 n=2:2:8、n=3:2:9 或者 n=2

3、:4:10、n=3:4:11 与 n=40 三种情况；（3）节点的选取分为均匀节点、切比雪夫节点两种四、程序流程图由于实验方案明显、简单，实现步骤及流程图省略五、实验结果具体结果在实验分析里：整理的结果如下1实验一的结果:1.f(x)11 25x2当节点为均匀节点时：插值点数目为奇数、偶数、40时，图像对称，但是不收敛,但是节点数越多，0附近的拟合效果越好，但是两端误差较大。当节点为切比雪夫点时：插值点数目为奇数、偶数、40时，图像对称，但是可以收敛，节点数越多，拟合效果越好。2.h(x)x1 x4当节点为均匀节点时：插值点数目为奇数、偶数、40时，图像对称，也是不收敛,但是节点数越多，0附近

4、的拟合效果越好，同时两端的误差较大。当节点为切比雪夫点时：插值点数目为奇数、偶数、40时，图像对称，但是可以收敛，节点数越多，拟合效果越好。3. g(x) arctan x当节点为均匀节点时：插值点数目为奇数、偶数、40时，图像对称，也是不收敛,但是节点数越多，0附近的拟合效果越好，同时两端的误差较大。当节点为切比雪夫点时：插值点数目为奇数、偶数、40时，图像对称，但是可以收敛，节点数越多，拟合效果越好。2实验二的结果通过作图可以发现：插值点数目增加时，三次样条插值光滑度依然很好，而且精度比以前更高，收敛性很好；但是发现lagra nge 插值却出现偏离，即存在误差，而且随着节点的增加，偏离越

5、明显。由此，可以发现,三次样条插值的收敛性比lagra nge插值好。3思考题结果通过分析计算可知，最咼点为:166178该点的高程为：721.098六、实验结果分析1实验一结果分析首先尝试了一些 n值，发现振荡明显，而且还有覆盖现象，由下图可见:1.5n=101n=20.5n=400.2-0.5-1-0.8-0.6-0.4-0.20.40.60.81y=1/(1+25*xn=3故针对上述现象，我们可以采用分开讨论测试的方法;°)f(X)1 25x21.节点为均匀节点时：Xi1 ?,i 0,1,2,n,na)当节点为奇数时，即n=2 : 2: 8，可以得到如下图像1.2i.LLLiL

6、LI.hI>n=7卡in=5=L2y=1/(1+25*x2)/ t fn=3/J7n=1J-rrriTrritr10.80.60.40.200.20.40.60.81-0.2-1-0.8-0.6-0.4-0.2从图中可以看到：节点数为基数个并且对称时，插值函数也是对称的；节点数越多，附近的区域拟合越好；节点数越多，两端误差越大；b)当节点为偶数时，即 n=3 : 2: 9，可以得到如下图像从图中可以看到：节点数为偶数个并且对称时，插值函数也是对称的；节点数越多；附近的区域拟合越好；节点数越多，两端误差越大；和奇数结果大致相同。c)当 n=40 时：16000140001200010000

7、80006000400020000-2000-1由图可知：插值函数也是左右对称，而且0附近几乎和被插值函数重合,但是两端误差很大，所以结论可以算是准确的。n=39-0.8-0.6-0.4-0.20.20.40.60.82.当节点为切比雪夫节点时 ：Xk 山 口cos （2k ° ,k222（n 1）即升=匚皿（二：；），£ = 1A- f十1,节点是对称的a）当节点为奇数个时，即n=2:2:8时，可以得到：1,2, ,n 110.80.60.40.20n=7n=5n=3 y=1/(1+25*x 2)-0.2 - - n=1-0.4 =:=:=:-1-0.8-0.6-0.4-

8、0.200.20.40.60.81从图中可以看出：节点数为基数个并且对称时，插值函数也是对称的;所有区域拟合都越好；b）当节点为偶数个时，即n=3:2:9时，可以得到：10.90.80.70.60.50.40.30.20.10-1-0.8-0.6-0.4-0.200.20.40.60.81n=7f / V/n=52 y=1/(1+25*x 2)n=3n=1节点数越多,此时，节点的选取也是对称的，同样我们也看到插值函数的图像是对称的；观察结论与 节点数为奇数时几乎一样：节点数越多，所有区域拟合都越好；c)当n=40时，得到：1.41.20.60.8n=70.40.20 IIttL-1-0.8-0

9、.6-0.4-0.200.20.40.60.8从图中我们看到，插值函数左右对称，插值函数几乎和被插值函数重合。故而，上面的观察结论是正确的。(2)h(x)x1 x41.节点为均匀节点时:Xi5 g ,i0,1,2,nna)当节点为奇数时，即n=2 : 4: 10，可以得到如下图像从图中可以看到:节点数为基数个并且对称时，插值函数也是对称的；节点数越多，0附近的区域拟合越好；节点数越多，两端误差越大；b)n=10n=2n=6y=x/(1+x )n=40当节点为偶数时，即n=3 : 4: 11，可以得到如下图像10.80.60.40.20-0.2-0.4-0.6-0.8-1-5-4-3-2-101

10、2345从图中可以看到：节点数为偶数个并且对称时，插值函数也是对称的；节点数越多；附近的区域拟合越好；节点数越多，两端误差越大；和奇数结果大致相同。c）当 n=40 时：5x 1043210-1-2-3-4-5-4-3-2-1012345由图可知：插值函数也是左右对称，而且0附近几乎和被插值函数重合,但是两端误差很大，所以结论可以算是准确的。2.当节点为切比雪夫节点时:Xk晋 穿cos轨，k 1,2, ,n 1即 ：:-：,节点是对称的a）当节点为奇数个时，即n=2:4:10时，可以得到:从图中可以看出，插值函数过两端和原点，并且也是奇函数；出现误差增大的现象；b）当节点为偶数个时，即n=3:

11、4:11时，可以得到:从图中可以看出，插值函数不经过两端，但也是奇函数；节点数越多，拟合度也越好c）当n=40时，得到：n越大拟合度越好，没有N取得很大的时候，插值函数和被插值函数几乎重合(3) g(x)arctan x10i1.节点为均匀节点时： Xi 5,i0,1,2,nna)当节点为奇数时，即n=2 : 4: 10，可以得到如下图像1.50.5n=9n=5y=ata n xn=1-0.5-1-1.5-5-4-3-2-101234从图中可以看出：节点数为基数个并且对称时，插值函数也是对称的；节点数越多, 近的区域拟合越好；节点数越多，两端误差越大；b)当节点为偶数时，即n=3 : 4: 1

12、1，可以得到如下图像从图中可以看到：节点数为偶数个并且对称时，插值函数也是对称的；节点数越多；附近的区域拟合越好；节点数越多，两端误差越大；和奇数结果大致相同。c）当 n=40 时：25002000 1500 1000500n=40-500-1000 -1500-2000-2500 1L-5-4-3-2-1012345由图可知：插值函数也是左右对称，而且0附近几乎和被插值函数重合,但是两端误差很大，所以结论可以算是准确的。2.当节点为切比雪夫节点时 ：Xk山 口cos （2k ° ,k 1,2, ,n 1222（n 1）即升=匚皿（：二；），£ = 1A- w十1,节点是对

13、称的a）当节点为奇数个时，即n=2:4:10时，可以得到：从图中可以看出，插值函数过两端和原点，并且也是奇函数；出现误差增大的现象；b）当节点为偶数个时，即n=3:4:11时，可以得到:n越大拟合度越好，没有从图中可以看出，插值函数不经过两端，但也是奇函数；节点数越多，拟合度也越好c）当n=40时，得到：N取得很大的时候，插值函数和被插值函数几乎重合综合分析上面的图像和数据可以发现：节点数目的奇偶对实验没有什么影响，而且节点不是越多拟合越好，可能会发生发 散现象，对称的节点选取， 得到的插值函数的对称性与被插值函数相同节点的位置不对 称，贝U得到的插值函数也不对称节点位置的选取会影响插值函数的

14、收敛性和误差切比雪 肤插值节点确实比以上用到的均匀节点、不均匀节点要好，对于不同的被插值函数，同 样的插值节点选取往往能得到类似的结果。2实验二结果分析11.令f（X）2则下面就以这个函数做具体分析:1 25x节点为5时，图如下节点为10时，图如下节点为20时，图如下:由上面三幅图可以发现，随着节点数的增加，三次样条插值函数的拟合越来越好，与原 函数越来越逼近；相反 I agrange插值函数的拟合却越来越差，误差越来越明显，因此，我 们可以得出结论，随着节点数的增加，三次样条插值函数的收敛性比 lagra nge插值函数好。故而我们可以证明了样条插值的收敛性。2.针对二中的一些数据我们可以进

15、行拟合一下，在图上面观察的更加明显。3.53 _-2.5 _2 1.5 _1 -0.5-_0rrrrrrrrr，012345678910以上就是在给定数据的条件下，做出的三次样条插值设计车门的曲线。3思考题结果分析可以利用 matlab进行编程，分别利用linear ， cubic ， spline 计算此题,计算的结果：利用lin ear 可以计算得到x =166 y =178zmax = 7.210983601652320e+02利用cubic可以计算得到：x=166 y=178zmax=7.210983601652320e+02利用spline可以计算得到：x=166 y=178zmax

16、=7.210983601652320e+02721.098所以我们可以得到，在不同的条件下，最高点应该是在 166 178附近，且最高程为总之，通过以上实验我们可以知道插值节点的数目不一定是越多拟合得越好，很多 时候会出现发散现象；对称的节点选取，得到的插值函数的对称性与被插值函数相同； 节点的位置不对称，贝U得到的插值函数也不对称；节点位置的选取会影响插值函数的收 敛性和误差；切比雪夫插值节点确实比以上用到的；均匀节点效果好；对于不同的被插 值函数，同样的插值节点选取往往能得到类似的结果。而且通过实验二我们可以知道体会到有些三次样条插值的收敛性很难证明时，有时 候利用图片，观察图像，也同样可

17、以获得证明的效果。实验二和思考题，重在将实验和理论相结合，学会在实际中运用所学知识，解决实 际问题。教师 评 语指导教师：年月日数值分析上机实验原始记录实验时间： 2013 年实验名称： 多项式插值的振荡现象班级：10 月 23 日 姓名： 学号：实验一的关键程序：事先估计程序m=150; x=-1:2/(m-1):1; y=1./(1+25*x92); z=0*x;plot(x,z,'r',x,y,'k-'), gtext('y=1/(1+25*xA2)'),pausen=3; x0=-1:2/( n-1):1; y0=1./(1+25*x0.

18、A2); y1=lagr1(x0,y0,x);hold on, plot(x,y1,'g'),gtext('n=2'),pause,hold offn=4; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y2=lagr1(x0,y0,x);hold on, plot(x,y2,'b:'),gtext('n=3'),pause,hold offn=5; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y3=lagr1(x0,y0,x);hold on, plot(x,y3,

19、9;r'),gtext('n=4'),pause,hold offn=11; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y4=lagr1(x0,y0,x);hold on, plot(x,y4,'r:'),gtext('n=10'),pause,hold offn=51; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y5=lagr1(x0,y0,x);hold on, plot(x,y5,'m'),gtext('n=50'),pause,hold

20、 offn=81; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y6=lagr1(x0,y0,x);hold on, plot(x,y6,'m'),gtext('n=80'),pause,hold off%1.1.a y=1./(1+25*x.A2) 的程序均匀节点m=150; x=-1:2/(m-1):1; y=1./(1+25*x.A2); z=0*x;plot(x,z,'r',x,y,'k-'), gtext('y=1/(1+25*xA2)'),pausen=2; x0=-1:2

21、/(n-1):1; y0=1./(1+25*x0.A2); y1=lagr1(x0,y0,x);hold on, plot(x,y1,'g'),gtext('n=1'),pause,hold offn=4; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y2=lagr1(x0,y0,x);hold on, plot(x,y2,'b:'),gtext('n=3'),pause,hold offn=6; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y3=lagr1(x0,y0,

22、x);hold on, plot(x,y3,'r'),gtext('n=5'),pause,hold offn=8; x0=-1:2/(n-1):1; y0=1./(1+25*x0.A2); y4=lagr1(x0,y0,x);hold on, plot(x,y4,'r:'),gtext('n=7'),pause,hold off%1.1b 切比雪夫节点m=150; x=-1:2/(m-1):1; y=1./(1+25*x.A2); z=0*x;plot(x,z,'r',x,y,'k-'),gtex

23、t('y=1/(1+25*xA2)'),pausex=q(2); y=1./(1+25*x.A2); n=2; x0=-1:2/( n-1):1;y1=lagr1(x,y,xO);hold on, plot(x0,y1,'g'),gtext('n=1'),pause,hold offx=q(4); y=1./(1+25*x.A2);n=4; x0=-1:2/(n-1):1; y2=lagr1(x,y,x0);hold on, plot(x0,y2,'b:'),gtext('n=3'),pause,hold offx

24、=q(6); y=1./(1+25*x.A2); n=6;x0=-1:2/(n-1):1;y3=lagr1(x,y,x0);hold on, plot(x0,y3,'r'),gtext('n=5'),pause,hold offx=q(8); y=1./(1+25*x.A2);n=8;x0=-1:2/(n-1):1;y4=lagr1(x,y,x0);hold on, plot(x0,y4,'r:'),gtext('n=7'),pause,hold off%2.1 a y=x./(1+x.A4) 的程序均匀节点m=150; x=-5

25、:10/(m-1):5; y=x./(1+x.A4); z=0*x;plot(x,z,'r',x,y,'k-'), gtext('y=x/(1+xA4)'),pausen=2; x0=-5:10/(n-1):5; y0=x0./(1+x0.A4); y1=lagr1(x0,y0,x);hold on, plot(x,y1,'g'),gtext('n=1'),pause,hold offn=6;x0=-5:10/(n-1):5; y0=x0./(1+x0.A4); y3=lagr1(x0,y0,x);hold on,

26、 plot(x,y3,'r'),gtext('n=5'),pause,hold offn=10; x0=-5:10/(n-1):5; y0=x0./(1+x0.A4); y4=lagr1(x0,y0,x);hold on, plot(x,y4,'r:'),gtext('n=9'),pause,hold off%2 1 b切比雪夫节点m=150; x=-5:10/(m-1):5; y=x./(1+x.A4); z=0*x;plot(x,z,'r',x,y,'k-'), gtext('y=x/(

27、1+xA4)'),pausex=q(2); y=x./(1+x.A4);n=3; x0=-5:10/(n-1):5;y1=lagr1(x,y,x0);hold on, plot(x0,y1,'g'),gtext('n=1'),pause,hold offx=q(6); y=x./(1+x.A4);n=5; x0=-5:10/(n-1):5; y2=lagr1(x,y,x0);hold on, plot(x0,y2,'b:'),gtext('n=5'),pause,hold offx=q(10); y=x./(1+x.A4)

28、; n=10;x0=-5:10/(n-1):5;y3=lagr1(x,y,x0);hold on, plot(x0,y3,'r'),gtext('n=9'),pause,hold off%3.1 a y=atan(x) 的程序均匀节点m=150; x=-5:10/(m-1):5; y=atan(x); z=0*x;plot(x,z,'r',x,y,'k-'),gtext('y=atan x'),pausen=2; x0=-5:10/(n-1):5; y0=atan(x0); y1=lagr1(x0,y0,x);ho

29、ld on, plot(x,y1,'g'),gtext('n=1'),pause,hold offn=6;x0=-5:10/(n-1):5; y0=atan(x0); y3=lagr1(x0,y0,x);hold on, plot(x,y3,'r'),gtext('n=5'),pause,hold offn=10; x0=-5:10/(n-1):5; y0=atan(x0); y4=lagr1(x0,y0,x);hold on, plot(x,y4,'y:'),gtext('n=9'),pause,

30、hold off%3.1 b切比雪夫节点m=150;x=-5:10/(m-1):5; y=atan(x); z=0*x;plot(x,z,'r',x,y,'k-'), gtext('y=1/(1+25*xA2)'),pausex=q(1); y=atan(x);n=2; x0=-5:10/(n-1):5;y1=lagr1(x,y,x0);hold on,plot(x0,y1,'g'),gtext('n=1'),pause,hold offx=q(5); y=atan(x); n=6;x0=-5:10/(n-1):5

31、;y3=lagr1(x,y,x0);hold on,plot(x0,y3,'r'),gtext('n=5'),pause,hold offx=q(9); y=atan(x);n=10;x0=-5:10/(n-1):5;y4=lagr1(x,y,x0);hold on,plot(x0,y4,'r:'),gtext('n=9'),pause,hold off实验二的关键程序问题一的程序x0=linspace(-1,1,11);y0=1./(1+9*x0.A2);x=-1:0.02:1;y=lagr1(x0,y0,x);yi=splin

32、e(x0,y0,x)% 求三次样条插值z=1./(1+9*x.A2);Ri=abs(z-yi)./z);% 三次样条插值相对误差R=abs(z-y)./z);%lagrange插值相对误差x,y,z,yi,R,Ri=x',y',z',yi',R',Ri'n=size(x0)plot(x,z,'o',x,y,'-',x,yi,'r*')legend(' 原始图像 ','lagrange 插值 ',' 三次样条插值 ')问题二的程序m=101;x=-1:2/(m-1):1;y=1./(1+25*x.A2);