资源描述:
《龙格现象实验报告.docx》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、精品文档数值计算实验报告(四)实验名称:观察龙格(Rung*现象实验姓名班级:学号:实验日期:周次:实验地点:实验目的.观察拉格朗日插值的龙格(Runge)现象.实验内容_.5对于函数f(x)——2进行拉格朗日插值,取不同的节点数n,在区间[-5,5]上取等距间隔的a+x节点为插值点,把f(x)和插值多项式的曲线画在同一张图上进行比较。(a可以取任意值)具体步骤:1、a=1时,1)取n=4,作出f(x)和插值多项式的曲线图;2)取n=10,作出f(x)和插值多项式的曲线图;2、a=0.25时,1)取n=4,作出f(x)和插值多项式的曲线图;2)取n=10
2、,作出f(x)和插值多项式的曲线图;3、分析上述曲线图,你可以得出什么结论?实验步骤拉格朗日插值的调用函数functiony=lagrange(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);L=0.0;forj=1:nT=1.0;fork=1:nifk~=jT=T*(z-x0(k))/(x0(j)-x0(k));endendL=T*y0(j)+L;endy(i)=L;enda=1时,1)取n=4x0=[-5:2:5];y0=5./(1+x0.A2);x=[-5:0.1:5];y=lagrange(x0,y
3、0,x);y1=5./(1+x.A2);.。1欢迎下载精品文档plot(x,y,'--r')holdonplot(x,y1,'-b')holdoff2)取r=10x0=[-5:1:5];y0=5./(1+x0.A2);x=[-5:0.1:5];y=lagrange(x0,y0,x);y1=5./(1+x.A2);plot(x,y,'--r')holdorplot(x,y1,'-b')holdoff2、a=0.25时,1)取r=4x0=[-5:2:5];y0=5./(0.25*0.25+x0.A2);x=[-5:0.1:5];y=lagrange(x0,y
4、0,x);y1=5./(1+x.A2);plot(x,y,'--r')holdorplot(x,y1,'-b')holdoff2)取r=10x0=[-5:1:5];y0=5./(0.25*0.25+x0.A2);x=[-5:0.1:5];y=lagrange(x0,y0,x);y1=5./(1+x.A2);plot(x,y,'--r')holdorplot(x,y1,'-b')holdoffa=1时,_02欢迎下载精品文档1)取n=4实验结果分析-5-4-3-2-1012345a=1时,2)取n=10-5-4-3-2-1012345a=0.25时,3欢迎
5、下载精品文档1)取n=45004003002001000a=0.25时,2)取n=10-5-4-3-2-100上述现象和定理告诉我们,并不是插值多项式的次数越高(即插值节点越多)精度越高,_。4欢迎下载精品文档从数值计算上可解释为高次插值多项式的计算会带来舍入误差的增大,从而弓1起计算失真。因此,实际应用做插值时一般只用一次、二次最多用三次插值多项式。那么如何提高插值精度呢?采用分段插值是一种办法。教师评语6欢迎下载精品文档欢迎您的下载,资料仅供套考!致力为企业和个人提供合同协议,策划案计划书,学习资料等等打造全网一站式需求6欢迎下载