数值逼近蒋尔雄3样条插值和曲线拟合答案

第三章 样条插值和曲线拟合
1.x y =
有如下的函数表
试用一次、二次、三次、四次多项式插值函数求8,看哪一个最接近8。 解  先作差商表
4
167
1210
13
9
3
42015
11008
16012
4
60
13
1611
1
10
0-⨯-
--
故:8.2)48(5
1
2)8(1=-+=p
819047619
.2)
98)(48(210
1
)48(512)8(2=----+=p
844444.2)98)(48)(18(3
4201)
48)(18(601
)
18(311)8(3=---⨯+----+=p
6222.2)1(4781008
1478601)
18(86
1
)08(10)8(4=-⨯⨯⨯-⨯⨯+---⨯+=p  已知
828427.28=,因此选定
)8(,16,9,42321p x x x ===最接近8。
利用Neville 方法得:
xi
8-xi
f(xi)      2.8284271
8
8
1
7
1
-1.33333333
3.3333333    2.4  4
4
2
2.866666667    2.6222222
2.8    2.8444444
9
-1    3
2.819047619
2.8571429    16 -8    4
f(8)=    2.828427125  xi
8-xi
f(xi)
8
8
1
7
1
-1  1/3
3  1/3      2  2/5    4
4
2
2  13/15      2  28/45
2  4/5      2  38/45  9
-1    3
2  86/105
书拉密女小站2  6/7      16 -8    4
已知 828427
.28=,故选定)8(,16,9,42321p x x x ====2.819047619最接近8.
1
11
01
20
1hanp
01
1
12
12
1
343
43
明溪蓝宝石
42
121------ 所以:)
())(1())(1()1(1)(21
213421344-++-++++-=x x x x x x x x x p , 故:0232.0)1.0(4=p 与f (0.1)=0.1不相等。
(2)若采用分段插值,则在],0[21
上,x x f x f x L =--+--=00)(0)0()(2
1
212121,所以: )1.0(1.0)1.0(f L ==,结果一样。
3.解  (1)若记L (x )为tgx 在[]4
4,π
π-上的分段线性插值函数,则 ],[,8
)()(122+∈≤-=i i x x x h M
x L tgx x R
其中[]4max 4
4,2=''=-∈x g t M x ππ,欲使4
22210218)()(-<≤≤-=h h M x L tgx x R ,故  014142.01022=⨯<-h
(2)如果采用分段二次插值,若)(2x L 为tgx 在[]44,ππ-上的分段二次插值函数,则
]
,
[,3121
6))()((!3)()()(13312
新加坡一航班折返12+++∈≤---'''≤
-=i i i i i x x x h M x x x x x x f x L tgx x R ξ其中[]
16max 4
4,3='''=-∈x g t M x π
π,欲使432
103
121
616)()(-<≤-=h x L tgx x R ,应使:
解:若对在44π
π中用等距分段Hermite 3次插值,则在每个小区间1+i i 上,由第
二章定理8知:
1212)4(3,)()(!
4)
()(++≤≤--=
-i i i i x x x x x x f x H tgx ξξ
由于x
x x x x tg 5
33)4(cos sin 24cos sin 16)(+=,所以在[]44,π
π-上, 802
11
2421116)()4(=+=≤x tg
4
4212324
5161!480)()(max !480)(1h h x x x x x H tgx i i x x x i i =⨯≤--≤-+≤≤+
注意右端与i 无关,故在[]
4
4,π
π-上,有:4
3
5)(h x H tgx ≤-。 解:因为0
10sin cos sin )(x tx
txdt x x x f =
==⎰,所以 ⎰⎪⎭⎫ ⎝
⎛+==10)
(2cos sin )(dt k tx t x x dx d x f k k k k π,因此11)(10)(+=≤⎰k dt t x f k
k 。 (1) 若在区间[]1,0上用等距线性插值,则误差为:
]1,0[,24
18)()(2
22∈≤=-x h h M x L x f 对任
欲使4
10)()(-≤-x L x f ,只须 048988979.010244=⨯<
-h 。
(2) 若在区间[]1,0上用等距Hermite 3次插值,则误差为:
]1,0[,16
5!41)()(!4)()()(42
12)4(3∈⋅⨯≤--≤-+x h x x x x f x H x f i i 对任ξ
欲使4构造区
310)()(-≤-x H x f ,只须 66195.0)10165!4(4
1
4
=⨯⨯⨯<-h  (3) 若在区间[]1,0上用等距样条插值,则由定理5,有:
4
)
4(4
5
1161)
(16
1)()(h x f h x s x f ⋅⋅≤
-∞
欲使4
10)()(-≤-x s x f ,只须 299069756.0)10165(4
14
=⨯⨯<-h 。 6. 对sin x /x , 在[0,1]上取5个节点,求D 13次样条。 解  因2sin cos )(x x
x x x f -=',所以:
02
sin )0(0
=-='=x x f ,-0.30117)1(='f , 先作
差商表
1 0 -0.16615 0    1
-0.04154 -0.16304 0.25 0.989616 -0.12306
-0.15388 0.5 0.958851 -0.2 -0.13905 0.75
0.908852 -0.26952 -0.12658 1 0.841471 -0.30117  1 0.841471
从而形成三弯距方程组的增广矩阵:
2    1 0 0 0 -0.99688 0.5    2 0.5 0 0 -0.97827 0 0.5    2 0.5 0 -0.92326 0 0 0.5    2 0.5 -0.834
3 0
1
2 -0.7595
解得:M 1=-0.33438,M 2=-0.32813,M 3=-0.30966,M 4=-0.27977,M 5=-0.23987。样条曲线为:
7. 对f(x)=|x|,在[-1,1]上取5个等距节点,求3次自然样条插值。
解先作差商表
-1    1 -1 0
-0.5 0.5 -1    2
0 0    1 0
0.5 0.5    1
1    1
从而形成三弯距方程组的增广矩阵:
2 0.5 0 0
0.5    2 0.5 12
0 0.5    2 0
解得:M1=-0(已知),M2= -1.71429=-12/7,M3= 6.857143=48/7,M4= -1.71429=-12/7,M5=-0(已知)。样条曲线为:
8. 对f(x)=|x|,在[-1,1]上取6个等距节点,求3次周期样条插值函数。
解先作差商表
-1    1 -1 0
-
0.6 0.6 -1    1.25 -0.2 0.2 0    1.25 0.2 0.2    1 0 0.6 0.6    1 -2.5 1    1 -1  1.4
0.6
注意:)()(21x f x f n =
+,4.14.01121=+=-+=+x x x x n n ,从而形成三弯距方程组的增广矩阵:
2 0.5 0 0 0.5 0 0.5    2 0.5 0 0 7.5 0 0.5    2 0.5 0 7.5 0 0 0.5    2 0.5 0 0.5
0.5
2
-15
解得:M 1=M 6= -8.18182,M 2= 1.363636,M 3= 2.727273,M 4= 2.727273,M 5= 1.363636。样
条曲线为:
解:取节点521
432121,作差商表:
1
11
无谓空间01
20
1
01
1
12
12
12
121---- 对于自然样条,051==M M ,按公式(10)形成方程组:

本文发布于:2024-09-21 18:54:49,感谢您对本站的认可!

本文链接:https://www.17tex.com/xueshu/44290.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:插值   样条   等距   商表   方程组   弯距   分段
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议