基于数值方法的龙格函数分析.pdf

基于数值方法的龙格函数分析 口蔡杰燕云 (东北大学秦皇岛分校河北·秦皇岛066004) 摘要：通过使用不同的数值插值算法和函数逼近方法对提出的龙格函数进行分析总结，得出在对应不同类型 的算法时，龙格函数是否显出病态性，插值方法是否对于高阶的问题仍然适用的相关结论。对拉格朗日插值、牛 顿插值、切比雪夫多项式零点插值、三次样条插值、曲线拟合等数值算法通过具体的Matlab程序，做出他们的相 应的函数图形并进行结果分析，以得出每个数值方法的性质，适用范围和优缺点。 关键词：龙格函数 三次样条插值 多项式插值切比雪夫多项武零点 中图分类号：0175．8 文献标识码：A 文章编号：100r7．3973(2013)002-112．03 l问题提出 fork=l：N 用不同的数值方法对龙格函数f∞气乏亏舻插值和函数 la=y1∞： 逼近： forj=l：k-1 (I)在区间【．1，1】上分别取n=lO，20用俩组等距离节点对 la=la*(x-xl(j))／(xl∞·xl(j))； end． 龙格函数作拉格朗日插值； forj=k+l：n (2)在区间【．1，l】上分别取n=10，20用俩组等距离节点对 la=la*(x-xlO))／(xl(k)-xl(j))； 龙格函数作牛顿插值； end (3)在区间【-l，l】上利用11次切比雪夫多项式零点作插值 L=L+Ia； 点，构造龙格函数的10次拉格朗日插值多项式，最后和L。。∞ end 做比较； 2．3作出两种情况下的图形 (4)在区间【．1，1】上分别取n=ll，25用俩组等距离节点对 龙格函数作三次样条插值。 2对龙格函数作拉格朗日插值以及相应的Matlab实现 2．1拉格朗日插值多项式基本知识 下面讨论构造通过n+1个勋x，…五节点的n次插 值多项式乙∞，假设它满足条件： L．∞嘞，j=o，l。2，…，n 若11．次多项式／j(x)(iO，1，2，…，n)在n+1个节点勘而 …矗上满足条件： 五扣ol’2…n 懈)_{甚Z 就称这n+1个n次多项式lo(x)，，IQ)，…，，-Q)为节点知 知…五上的n次插值基函数。 经推导得n次插值基函数为： ‘(J)：j竺盐些_氲竖立立生霉4， … (‘-xo)…(薯一以一IX黾-x,+。)…吒一毛)’ k=O，1，2，…。n 于是满足条件的插值多项式可表示为： k(工)=D让∽． 2．2拉格朗日插值的实现算法 xl=linspaee(-1，l，n)； yl=1．／(I+25xi．^2)： symsx； N=length(x1)； L--0； 图2 20次Lagrange差值与原函数图像 ——斟协论瑶·2013年第2期(下)—— 万方数据 使用函数：fplot(fun，limits)在指定的范围limits内画出函P 数名为fun的图像。其中limits是一个指定X轴范围的向量3．3作出两种情况下的图形 [xminxmax](见图2)。 2．4