数值分析上机实验报告 插值

发布 2019-08-02 05:31:17 阅读 8068

数值分析第一次上机练习实验报告。

—lagrange插值与三次样条插值。

一、 问题的描述。

设, ,取,.试求出10次lagrange插值多项式和三次样条插值函数(采用自然边界条件),并用图画出, ,

二、 方法描述——lagrange插值与三次样条插值。

我们取,,通过在点的函数值来对原函数进行插值,我们记插值函数为,要求它满足如下条件:

我们在此处要分别通过lagrange插值(即多项式插值)与三次样条插值的方法对原函数进行插值,看两种方法的插值结果,并进行结果的比较。

10次的lagrange插值多项式为:

其中:以及。

我们根据(2)进行程序的编写,我们可以通过几个循环很容易实现函数的lagrange插值。

理论上我们根据区间上给出的节点做出的插值多项式近似于,而多项式的次数越高逼近的精度就越好。但实际上并非如此,而是对任意的插值节点,当的时候不一定收敛到;而是有时会在插值区间的两端点附近会出现严重的偏离的现象,即所谓的runge现象。因此用高次插值多项式近似的效果并不总是好的,因而人们通常在选择插值方式的时候不用高次多项式插值,而用分段低次插值,而这样的插值效果往往是非常好的,能够克服高次多项式插值的弱点,达到令人满意的效果。

分段低次插值包括分段线性插值、分段三次hermite插值、三次样条插值等。前两种插值函数都具有一致收敛性,但是光滑性较差,而在实际问题中我们往往要求函数具有二阶光滑度,即有二阶连续导数。而对第三种插值方式,我们得到的是一个样条曲线,它是由分段三次曲线拼接而成,在连接点(即样点)上二阶导数连续。

我们记三次样条插值函数为,它在每个小区间上是三次函数,因此在每个区间上需要确定4个参数,总共有10个小区间,因此共需确定40个未知参数。首先我们有插值条件:

其次在每个节点上满足连续性条件:

此外在端点处满足自然边界条件:

我们假设。则在每个小区间上:

其中:及。

我们利用边界条件(3)(4)(5)可以得到:

其中:以及。

两端点处的边界条件为:

将边界条件写成矩阵形式为:

其中根据自然边界条件(8)有:

我们解方程(9)就可以得到,将他们代入(6)就可以得到各段区间上的的值。

三、 方案设计。

我们通过编写matlab程序来进行10次lagrange插值与三次样条插值的工作。在我们的程序文件中文件是主程序文件;文件是计算10次lagrange插值多项式的子程序文件,给它任一个,此程序将返回的值;是根据(9)计算各节点二阶导数值的子程序文件,它将会返回在自然边界条件下的各节点的二阶导数值;然后是根据以及(6)计算三次样条插值函数的子程序文件。然后运行主程序将给出三幅曲线图,分别是与曲线,与曲线,以及、与三条曲线共同画在一幅图上得到的图象。

解决这个问题的思路很简单,按部就班的来就可以。首先我们计算各节点上的函数值以备后用,然后调用计算。随后我们给出一系列的值,计算,并分别调用与分别计算与。

然后根据我们得到的数据绘图观察插值结果。具体程序的实现可参见所给程序的相关注释。

四、 计算结果及其分析。

下面是我们根据程序计算结果得到的数据,其中分别给出了在各典型处的的原函数的值、lagrange插值结果与样条插值结果;以及绝对误差和,相对误差,。由于在两端点处进行lagrange插值插值的时候可能出现runge现象,因此我们在两端点附近多给了几个点的数据。

尽管从数据可以看出一些端倪,但是通过图象我们更能清楚地看到最终插值结果的定性情况。首先我们给出与曲线:

其中蓝色的曲线代表曲线,绿色的曲线代表曲线。可见此时两者之间具有很大的差别,尤其在端点附近会出现严重的偏离的现象,即出现了所谓的runge现象。而此时曲线与我们用样条插值得到的的曲线为:

其中蓝色的曲线代表曲线,绿色的曲线代表曲线,可见两条曲线几乎完全重合,与符合的很好。

上面我们由曲线定性看到的结论也可以通过表中的数据定量的看出。

五、 结论。

插值方法中最基本的是多项式插值,而我们可以通过lagrange多项式来方便的实现这种插值方式。理论上我们根据给定区间上的给定的节点做出的插值多项式近似于,而多项式的次数越高逼近的精度就越好。但实际上对任意的插值节点,当的时候不一定收敛到;而是有时会在插值区间的两端点附近会出现严重的偏离的现象,即所谓的runge现象。

因此用高次插值多项式近似的效果并不总是好的,而我们通过本次试验中的实际计算发现对本次试验中的函数确实出现了runge现象,插值结果很不令人满意;我们转而采用分段的三次样条插值,得到了非常好的插值效果。

102440数值实验报告

一 实验内容 用两种不同的顺序计算,分析其误差的变化。二 实验目的 1.通过上机编程,复习巩固以前所学程序设计语言 2.通过上机计算,了解舍入误差所引起的数值不稳定性。3.通过上机计算,了解运算次序对计算结果的影响,从而尽量避免大数吃小数的现象。三 算法描述。取float类型数据,从两个方向即 1大...

东北大学数值分析实验报告

课题一迭代格式的比较。一 问题提出。设方程f x x 3x 1 0 有三个实根 x 1.8793 x 0.34727 x 1.53209现采用下面三种不同计算格式,求 f x 0的根 x 或x 1 x 2 x 3 x 二 要求。1 编制一个程序进行运算,最后打印出每种迭代格式的敛散情况 2 用事后误...

判别分析实验报告SPSS

一 实验目的及要求 1 目的。用spss软件实现判别分析及其应用。2 内容及要求。用spss对实验数据利用fisher判别法和贝叶斯判别法,建立判别函数并判定宿州 广安等13个地级市分别属于哪个管理水平类型。二 仪器用具 三 实验方法与步骤 准备工作 把实验所用数据从word文档复制到excel,并...