数值分析实验报告2——Runge现象

1、精品文档数值分析课程实验报告插值逼近题目一 .Runge 函数的插值1. Runge函数Runge函数的表达式为：R( x)1125 x2其在 -1,1 区间上的函数图像如图 1.1 。在课程学习中我们知道，对 Runge函数进行高次插值时有可能在两端出现不收敛的情况， 即 Runge现象。下面将分别用四种不同的插值方法在 -1,1 区间上对 Runge 函数进行插值，并分析是否产生Runge现象，比较插值效果。图函数在 -1,1区间的函数图像2.Newton 插值首先根据课本上的 Newton 插值算法进行编程 （代码略）。核心思想就是用符号变量进行中间运算， 以便将最终的插值函数用符号表达

2、式表示出来， 并进一步生成图像。此处插值节点选择为等距插值节点，即：。1欢迎下载精品文档xi0.1ih (i0,1,2,20)其中 h=0.1 。插值曲线与原曲线的对比如图 1.2 （蓝色为原曲线，红色为插值曲线）。从图中看出，在区间中部，二者吻合较好；但在区间两端二者则产生了明显偏差，甚至可以达到一个非常大的数值（ e20 量级）。因此，在等距节点的 20 次 Newton 插值下，产生了明显的 Runge现象。图插值曲线与原曲线对比3. Lagrange 插值此处同样是根据 Lagrange 插值的具体算法进行编程。但插值节点不再是等距分布，而是如下形式：2i1xicos()(i0,1,2

3、,20)插值曲线与原曲线的对比如图1.3 （蓝色为原曲线，红色为插值曲线） 。从图中看出，插值曲线与原曲线吻合的很好， 没有产生明显的 Runge现象。对比产生了明显 Runge现象的 20 次 Newton 插值， Lagrange 插值的最高次数虽然也是 20，但由于此处的插值节点不是等距分布的（事实上，此处采用的插值节点正是Chebyshev 多项式的零点），而是中间疏两边密，因此两侧较密的节点很好地抑制了 Runge现象。2欢迎下载精品文档图 1.3. Lagrange 插值曲线与原曲线对比4. 分段线性插值分段线性插值是这几种插值方法中最容易处理的一个， 只需要将每个节点对应的函数值

4、求出再将相邻的数据点两两用直线相连即可。此处采用了等距节点，所得插值曲线与原曲线对比如图 1.4 （蓝色为原曲线，红色为插值曲线） 。从图中。3欢迎下载Runge现象的倾向。精品文档图 1.4.分段线性插值曲线与原曲线对比看出，此处分段线性插值的效果也还是不错的，二者只在区间中部略微存在一些偏差，而在其他区域整体上吻合的很好， 并且不存在 Runge现象。这是由于分段线性插值通过对插值区间分段的方法将插值函数的次数有效降低， 因而即使是等距节点分布，也很好地避免了出现5. 三次样条插值三次样条插值是这四种插值方法中编程最麻烦的， 但并不是说存在多大的技术难度，只是因为插值过程中的步骤比较繁琐，

5、 因而代码也显得较为冗长。 此处依然采用等距节点，所得插值曲线与原曲线对比如图1.5 （蓝色为原曲线，红色为插值曲线）。从图中看出，三次样条插值的效果比分段线性插值更胜一筹，三次样条插值曲线和原曲线在整个插值区间都基本处于重合状态，几乎没有肉眼可见的偏差。同样，由于三次样条插值的插值函数最高次数只有3，在等距节点下也没有产生 Runge现象。4欢迎下载精品文档图 1.5. 三次样条插值曲线与原曲线对比。5欢迎下载精品文档题目二 . 分段函数的插值1. 分段函数定义在 -1,1区间的分段函数的函数表达式为：sinx, 1x0f ( x)cos x,01x20, 1x 12其函数图像如图 2.1

6、。分段函数最大的特点就是在个别点上函数值或导数值存在突变，因此可以预计， 除了可能出现的 Runge现象外，在那些突变点附近的插值结果也可能会出现较大的偏差。 下面将分别采用之前的四种插值方法在该函数的-1,1定义域内对其进行插值。图 2.1. 分段函数图像2.Newton 插值首先根据课本上的 Newton 插值算法进行编程。此处插值节点选择为等距插值节点，即：。6欢迎下载精品文档xi0.1 ih (i0,1,2,20)其中 h=0.1 。插值曲线与原曲线的对比如图2.2（蓝色为原曲线，红色为插值曲线）。从图中看出，与 Newton 法对 Runge函数的插值结果相比， Newton 法对于

7、该分段函数的插值效果显得更加糟糕： 不仅在区间两端产生了极强烈的震荡 （即Runge 现象），就连区间中部也存在较小的上下震荡。因此，从整体来看，几乎所有距插值节点稍远的点都存在较大的偏差， 这表明该分段函数在等距节点下的 20 次 Newton插值效果非常不理想。图插值曲线与原曲线对比3. Lagrange 插值此处同样是根据 Lagrange 插值的具体算法进行编程。但插值节点不再是等距分布，而是如下形式：xicos( 2i1)(i0,1,2,20)42插值曲线与原曲线的对比如图 2.3 （蓝色为原曲线，红色为插值曲线） 。从图中看出，与同样次数的 Newton 法相比， Lagrange

8、 法所得的插值曲线虽然在区间中部的分布与其相似，但在区间两端较好地收敛到了原曲线上，即较好地消除了Runge现象。这同样是因为此处的插值节点不是等距分布的（事实上，此处采用的插值节点正是 Chebyshev 多项式的零点），而是中间疏两边密，因此两侧较密的节点很好地抑制了 Runge现象。7欢迎下载精品文档图 2.3. Lagrange 插值曲线与原曲线对比4. 分段线性插值分段线性插值是这几种插值方法中最容易处理的一个， 只需要将每个节点对应的函数值求出再将相邻的数据点两两用直线相连即可。此处采用了等距节点，图 2.4.分段线性插值曲线与原曲线对比。8欢迎下载精品文档所得插值曲线与原曲线对比

9、如图 2.4 （蓝色为原曲线，红色为插值曲线） 。从图中看出，此处分段线性插值的效果较好， 二者只在区间中部函数值的突变点附近存在一些偏差， 而在其他区域整体上吻合的很好， 不存在 Runge现象。这是由于分段线性插值通过对插值区间分段的方法将插值函数的次数有效降低， 因而即使是等距节点分布，也很好地避免了出现 Runge现象的倾向。5. 三次样条插值三次样条插值是这四种插值方法中编程最麻烦的， 但并不是说存在多大的技术难度，只是因为插值过程中的步骤比较繁琐， 因而代码也显得较为冗长。 此处依然采用等距节点，所得插值曲线与原曲线对比如图 2.5 （蓝色为原曲线，红色为插值曲线）。从图中看出，三

10、次样条插值的效果与分段线性插值相近，也是在区间中部的函数值突变处有一定的偏差，而其他区域都吻合较好，也没有产生Runge现象。同样，这也是由于三次样条插值的插值函数最高次数只有 3，因此在等距节点下进行插值也没有产生 Runge现象。图 2.5. 三次样条插值曲线与原曲线对比。9欢迎下载精品文档三. 总结本文通过 matlab 编程分别采用Newton 插值法、 Lagrange 插值法、分段线性插值法以及三次样条插值法对 Runge函数和一个分段函数进行了插值逼近， 插值区间 -1,1 ，插值节点 21 个，并通过分析计算结果主要得到了以下结论：1. 插值多项式次数过高时会产生严重的 Run

11、ge现象。本实验中，无论是 Runge函数还是分段函数的 20 次 Newton 插值多项式都产生了严重的 Runge现象，区间两端处的插值出现剧烈震荡，严重失真。2. 同样是高次插值多项式，若适当选取插值节点能够在一定程度上抑制Runge现象。本实验中， Runge函数和分段函数的 20 次 Lagrange 插值多项式由于采用了中间疏两边密的非等距结点，而不是 Newton 插值多项式所用的等距节点，有效地抑制或消除了本应出现的 Runge现象（这里说“本应出现”是经过计算验证的，若取等距节点，则 20 次的 Lagrange 插值多项式也会出现严重 Runge 现象）。3. 降低插值多项

12、式的次数能有效避免 Runge现象。本实验中，分段线性插值法（各区间上均为 1 次）和三次样条插值法（最高次数为 3）都取得了较为理想的差值逼近效果， 没有出现 Runge现象，且在整个插值区间都与原函数的图像吻合的很好。4. 与连续函数相比，存在不连续点的分段函数的插值逼近误差更大，且更加不稳定。本实验中， 对连续的 Runge函数进行插值逼近时， 除了等距节点的高次Newton 多项式出现严重 Runge 现象，其余三种方法基本都收敛到了原曲线上，取得了不错的插值逼近效果； 而对分段函数进行插值逼近时， 除了等距节点的高次 Newton 多项式的逼近效果非常糟糕外 （巨大偏差， 严重震荡），其余三种方法虽然没有出现 Runge现象，但在不连续点（x=0）的附近区