数值分析第一次上机练习实验报告
——Lagrange插值与三次样条插值
1、 问题的描述
设
,
,取
,
.试求出10次Lagrange插值多项式
和三次样条插值函数
(采用自然边界条件),并用图画出
,
,
.
2、
方法
快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载
描述——Lagrange插值与三次样条插值
我们取
,
,通过在
点的函数值
来对原函数进行插值,我们记插值函数为
,要求它满足如下条件:
(1)
我们在此处要分别通过Lagrange插值(即多项式插值)与三次样条插值的方法对原函数
进行插值,看两种方法的插值结果,并进行结果的比较。
10次的Lagrange插值多项式为:
(2)
其中:
以及
我们根据(2)进行程序的编写,我们可以通过几个循环很容易实现函数的Lagrange插值。
理论上我们根据区间
上给出的节点做出的插值多项式
近似于
,而多项式
的次数
越高逼近
的精度就越好。但实际上并非如此,而是对任意的插值节点,当
的时候
不一定收敛到
;而是有时会在插值区间的两端点附近会出现严重的
偏离
的现象,即所谓的Runge现象。因此用高次插值多项式
近似
的效果并不总是好的,因而人们通常在选择插值方式的时候不用高次多项式插值,而用分段低次插值,而这样的插值效果往往是非常好的,能够克服高次多项式插值的弱点,达到令人满意的效果。
分段低次插值包括分段线性插值、分段三次Hermite插值、三次样条插值等。前两种插值函数都具有一致收敛性,但是光滑性较差,而在实际问题中我们往往要求函数具有二阶光滑度,即有二阶连续导数。而对第三种插值方式,我们得到的是一个样条曲线,它是由分段三次曲线拼接而成,在连接点(即样点)上二阶导数连续。
我们记三次样条插值函数为
,它在每个小区间
上是三次函数,因此在每个区间上需要确定4个参数,总共有10个小区间,因此共需确定40个未知参数。首先我们有插值条件:
(3)
其次在每个节点
上满足连续性条件:
(4)
此外在端点处满足自然边界条件:
(5)
我们假设
。则在每个小区间
上:
(6)
其中:
及
我们利用边界条件(3)(4)(5)可以得到:
(7)
其中:
以及
两端点处的边界条件为:
(8)
将边界条件写成矩阵形式为:
(9)
其中根据自然边界条件(8)有:
我们解方程(9)就可以得到
,将他们代入(6)就可以得到各段区间上的
的值。
3、
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
我们通过编写Matlab程序来进行10次Lagrange插值与三次样条插值的工作。在我们的程序文件中interplotion.m文件是主程序文件;L10.m文件是计算10次Lagrange插值多项式
的子程序文件,给它任一个
,此程序将返回
的值;Mspline.m是根据(9)计算各节点二阶导数值
的子程序文件,它将会返回在自然边界条件下的各节点的二阶导数值
;然后spline.m是根据
以及(6)计算三次样条插值函数
的子程序文件。然后运行主程序将给出三幅曲线图,分别是
与
曲线,
与
曲线,以及
、
与
三条曲线共同画在一幅图上得到的图象。
解决这个问题的思路很简单,按部就班的来就可以。首先我们计算各节点
上的函数值
以备后用,然后调用Mspline.m计算
。随后我们给出一系列
的值,计算
,并分别调用L10.m与spline.m分别计算
与
。然后根据我们得到的数据绘图观察插值结果。具体程序的实现可参见所给程序的相关
注释
圣经注释小学小古文100篇及注释小古文100篇及注释简短小古文100篇及注释译文小古文100篇及注释
。
4、 计算结果及其分析
下面是我们根据程序计算结果得到的数据,其中分别给出了在各典型
处的的原函数的值
、Lagrange插值结果
与样条插值结果
;以及绝对误差
和
,相对误差
,
。由于在两端点处进行Lagrange插值插值的时候可能出现Runge现象,因此我们在两端点附近多给了几个点的数据。
x
f(x)
L10(x)
S(x)
误差
L10(x)-f(x)
相对误差
(L10(x)-f(x))/ f(x)
误差
S(x)-f(x)
相对误差
(S(x)-f(x))/ f(x)
-1.00
0.10000
0.10000
0.10000
0.00000
0.00000
0.00000
0.00000
-0.98
0.10370
0.41597
0.10409
0.31228
3.01147
0.00040
0.00383
-0.96
0.10759
0.56881
0.10823
0.46122
4.28677
0.00064
0.00592
-0.94
0.11170
0.61150
0.11245
0.49980
4.47440
0.00075
0.00669
-0.92
0.11604
0.58445
0.11680
0.46841
4.03656
0.00076
0.00651
-0.90
0.12063
0.51779
0.12131
0.39716
3.29244
0.00069
0.00569
-0.88
0.12548
0.43329
0.12604
0.30781
2.45313
0.00057
0.00452
-0.86
0.13061
0.34608
0.13103
0.21547
1.64974
0.00042
0.00320
-0.84
0.13605
0.26605
0.13631
0.13001
0.95560
0.00026
0.00192
-0.82
0.14181
0.19906
0.14193
0.05725
0.40368
0.00012
0.00082
-0.80
0.14793
0.14793
0.14793
0.00000
0.00000
0.00000
0.00000
-0.70
0.18484
0.11010
0.18484
-0.07474
-0.40433
0.00000
0.00000
-0.60
0.23585
0.23585
0.23585
0.00000
0.00000
0.00000
0.00000
-0.50
0.30769
0.33383
0.30640
0.02614
0.08494
-0.00129
-0.00419
-0.40
0.40984
0.40984
0.40984
0.00000
0.00000
0.00000
0.00000
-0.30
0.55249
0.53935
0.55629
-0.01314
-0.02378
0.00380
0.00689
-0.20
0.73529
0.73529
0.73529
0.00000
0.00000
0.00000
0.00000
-0.10
0.91743
0.92247
0.91516
0.00504
0.00549
-0.00227
-0.00247
0.00
1.00000
1.00000
1.00000
0.00000
0.00000
0.00000
0.00000
0.10
0.91743
0.92247
0.91516
0.00504
0.00549
-0.00227
-0.00247
0.20
0.73529
0.73529
0.73529
0.00000
0.00000
0.00000
0.00000
0.30
0.55249
0.53935
0.55629
-0.01314
-0.02378
0.00380
0.00689
0.40
0.40984
0.40984
0.40984
0.00000
0.00000
0.00000
0.00000
0.50
0.30769
0.33383
0.30640
0.02614
0.08494
-0.00129
-0.00419
0.60
0.23585
0.23585
0.23585
0.00000
0.00000
0.00000
0.00000
0.70
0.18484
0.11010
0.18484
-0.07474
-0.40433
0.00000
0.00000
0.80
0.14793
0.14793
0.14793
0.00000
0.00000
0.00000
0.00000
0.82
0.14181
0.19906
0.14193
0.05725
0.40368
0.00012
0.00082
0.84
0.13605
0.26605
0.13631
0.13001
0.95560
0.00026
0.00192
0.86
0.13061
0.34608
0.13103
0.21547
1.64974
0.00042
0.00320
0.88
0.12548
0.43329
0.12604
0.30781
2.45313
0.00057
0.00452
0.90
0.12063
0.51779
0.12131
0.39716
3.29244
0.00069
0.00569
0.92
0.11604
0.58445
0.11680
0.46841
4.03656
0.00076
0.00651
0.94
0.11170
0.61150
0.11245
0.49980
4.47440
0.00075
0.00669
0.96
0.10759
0.56881
0.10823
0.46122
4.28677
0.00064
0.00592
0.98
0.10370
0.41597
0.10409
0.31228
3.01147
0.00040
0.00383
1.00
0.10000
0.10000
0.10000
0.00000
0.00000
0.00000
0.00000
尽管从数据可以看出一些端倪,但是通过图象我们更能清楚地看到最终插值结果的定性情况。首先我们给出
与
曲线:
其中蓝色的曲线代
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
曲线,绿色的曲线代表
曲线。可见此时两者之间具有很大的差别,尤其在端点附近会出现严重的
偏离
的现象,即出现了所谓的Runge现象。而此时
曲线与我们用样条插值得到的
的曲线为:
其中蓝色的曲线代表
曲线,绿色的曲线代表
曲线,可见两条曲线几乎完全重合,
与
符合的很好。
上面我们由曲线定性看到的结论也可以通过表中的数据定量的看出。
5、 结论
插值方法中最基本的是多项式插值,而我们可以通过Lagrange多项式来方便的实现这种插值方式。理论上我们根据给定区间上的给定的节点做出的插值多项式
近似于
,而多项式
的次数
越高逼近
的精度就越好。但实际上对任意的插值节点,当
的时候
不一定收敛到
;而是有时会在插值区间的两端点附近会出现严重的
偏离
的现象,即所谓的Runge现象。因此用高次插值多项式
近似
的效果并不总是好的,而我们通过本次试验中的实际计算发现对本次试验中的函数确实出现了Runge现象,插值结果很不令人满意;我们转而采用分段的三次样条插值,得到了非常好的插值效果。
_1285829558.unknown
_1285835391.unknown
_1285837061.unknown
_1285852019.unknown
_1285852583.unknown
_1285855858.unknown
_1285856858.unknown
_1285856952.unknown
_1285856963.unknown
_1285856792.unknown
_1285855814.unknown
_1285855834.unknown
_1285855787.unknown
_1285852210.unknown
_1285852557.unknown
_1285852420.unknown
_1285852511.unknown
_1285852087.unknown
_1285837451.unknown
_1285852014.unknown
_1285851910.unknown
_1285851943.unknown
_1285837185.unknown
_1285837352.unknown
_1285837125.unknown
_1285836005.unknown
_1285836109.unknown
_1285836319.unknown
_1285836015.unknown
_1285835981.unknown
_1285835989.unknown
_1285835830.unknown
_1285835947.unknown
_1285832641.unknown
_1285833215.unknown
_1285833291.unknown
_1285833402.unknown
_1285833304.unknown
_1285833248.unknown
_1285833283.unknown
_1285832827.unknown
_1285833191.unknown
_1285832706.unknown
_1285832624.unknown
_1285829712.unknown
_1285829805.unknown
_1285829595.unknown
_1285829682.unknown
_1285685134.unknown
_1285829471.unknown
_1285829520.unknown
_1285685825.unknown
_1285680642.unknown
_1285684351.unknown
_1285684362.unknown
_1285681182.unknown
_1285680344.unknown