数值计算课后答案
习 题 四 解 答
1、设010,1x x ==,写出()x f x e -=的一次插值多项式1()L x ,并估计插值误差。 解:根据已知条件,有
设插值函数为1()L x ax b =+,由插值条件,建立线性方程组为
1
01
1a b a b e -?+=???+=? 解之得11
1a e b -?=-?=?
则11()(1)1L x e x -=-+ 因为(),()x x y x e y x e --'''=-= 所以,插值余项为
(1)(2)
(2)011
()()()()()
(1)!
1()()2!1
()()()2!1
(0)(1)((0,1))2n r x f x p x f x n f x f x x x x e x x ξξπξπξξ+-=-=+=
=--=--∈
所以
01
0101
()max max (1)
2111248x r x e x x e ξξ-≤≤≤≤-≤-=??=。
2、给定函数表
i x
()i f x
选用合适的三次插值多项式来近似计算f 和f 。
解:设三次插值多项式为230123()f x a a x a x a x =+++,由插值条件,建立方程组为
23012323
012323
01232301
23(0.1)(0.1)(0.1)0.9950.30.30.30.995
0.70.70.70.7651.1 1.1 1.10.454
a a a a a a a a a a a a a a a a ?+?-+?-+?-=?+?+?+?=??+?+?+?=??+?+?+?=? 即
012301230123
123012312301230.10.010.0010.9950.10.010.0010.9950.30.090.0270.9950.40.080.02800.70.490.3430.7650.80.480.344 1.761.1 1.21 1.3310.454a a a a a a a a a a a a a a a a a a a a a a a a a a -+-=-+-=??+++=++=???
+++=++=??+++=?12301231232330.40.720.9880.3110.10.010.0010.9950.40.080.02800.320.288 1.760.384 3.831a a a a a a a a a a a a a ??????++=-?
-+-=??++=???
+=?
?-=-?
解之得 01
230.416.293.489.98
a a a a =??=-??
=-??=? 则所求的三次多项式为23()0.41 6.29 3.489.98f x x x x =--+。 所以
2323
(0.2)0.41 6.290.2 3.480.29.980.20.91
(0.8)0.41 6.290.8 3.480.89.980.8 1.74f f =-?-?+?=-=-?-?+?=-
3、设(0,1,2,,)i x i n =是 n+1个互异节点,证明: (1)0()(0,1,2,
,)n
k k i i i x l x x k n ===∑;
(2)0
()()0(0,1,2,
,)n k i i i x x l x k n =-==∑。
证明: (1)由拉格朗日插值定理,以x 0,x 1,x 2,…x n 为插值节点,对y=f(x)=x k 作n 次插值,插值多项式为 0()()n
n i i i p x l x y ==∑,
而y i =x i k ,
所以0
()()()n
n
k n i i i i i i p x l x y l x x ====∑∑
同时,插值余项
(1)(1)11
()()()()()()0(1)!(1)!
n k n k n r x x p x f x x x n n ξξππ++=-=
==++
所以0
()n
k k i i i l x x x ==∑
结论得证。
(2)取函数()(),0,1,2,,k f x x t k n =-=
对此函数取节点(0,1,2,,)i x i n =,则对应的插值多项式为
0()()()n
k n i i i p x x t l x ==-∑,
由余项公式,得
(1)
(1)011
()()()()()()()()0(1)!(1)!
n
n k
k n k
i i i r x x t x t l x f x x t x n n ξ
ξππ++==---=
=-=++∑所以
()()()n
k
k i i i x t x t l x =-=-∑
令t=x ,
()()0n
k
i
i
i x x l x =-=∑
4
、给定数据(()f x =
x
f(x)
(1)试用线性插值计算f 的近似值,并估计误差;
(2)试用二次Newton 插值多项式计算f 的近似值,并估计误差。 解:用线性插值计算f ,取插值节点为和,则相应的线性插值多项式是
1.54919 1.48320
() 1.48320( 2.2)
2.4 2.2
1.483200.32995(
2.2)
p x x x -=+--=+- 用x=代入,得
(2.3) 1.483200.32995(2.3 2.2) 1.450205f ≈+?-= (2)
根据定理2,
f(x)=f(x 0)+f[x 0,x 1](x-x 0)+f[x 0,x 1,x 2](x-x 0)(x-x 1)+…
+f[x 0,x 1,…,x n ](x-x 0)(x-x 1)…(x-x n -1)
+f[x 0,x 1,…,x n ,x]π(x) 。 以表中的上方一斜行中的数为系数,得 f =+ × × × = 指出: 误差未讨论。 5、给定函数表
x 0
1
2
4
5 y
16
46
88
试求各阶差商,并写出牛顿插值多项式和插值余项。
57
()0167(1)(1)(2)(1)(2)(4)26
p x x x x x x x x x x x =++--------。
指出: 余项未讨论。 5*、给定函数表
x 0
1
2
3
4 y
16
46
88
试求各阶差分,并求等距节点插值。
解:由已知条件,显然,x 0=0,h=1,x=t 。 0(1)(1)(2)(1)(2)(3)
()()01614(2)(140)2!3!4!(1)(2)35
167(1)(1)(2)(3)
36
n n t t t t t t t t t p x th p t t t t t t t t t t t t ------+==+?+
?+?-+?---=+------指出:
在本题这种情况下,实际上()()n n p t p x =,也就是说,在这样的条件下,t 的多项式就是x 的多项式,可以直接转换。
一般情况下,把t 的关系转换为x 的关系需要根据x=x 0+th ,将t 用x 表示,即将
x x t h
-=代入得到的多项式。 6
解:所给节点是等距结点:
000.125,0.125,,0,1,2,3,4,5i x h x x ih i ===+=。
令0
0()x x th t h
=+=
,根据等距结点插值公式,得 0(1)
()()0.79618(0.02284)(0.00679)
2!
(1)(2)(1)(2)(3)(1)(2)(3)(4)(0.00316)0.00488(0.00460)
3!4!5!
n n t t p x th p t t t t t t t t t t t t t t -+==+?-+?----------+?-+?+?-则
(0.1581)(0.1581)(0.1250.2648)0.790294822,
(0.636)(0.6363)(0.125 4.088)0.651804826n n n n f p p h f p p h ≈=+=≈=+=。
7、设f(x)在[-4,4]有连续的4阶导数,且 (1)1,(0)2,(0)0,(3)1,(3)1f f f f f ''-=====
(1)试构造一个次数最低的插值多项式p(x),使其满足
(1)(1)1,(0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f p f ''''-=-=-======== (2)给出并证明余项f(x)-p(x)的表达式。 解:
(1)由7*可以求出满足
(0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''======== 的三次埃尔米特插值多项式
32
52()2273
H x x x =-+。
设2232
2252()()(3)2(3)273
p x H x a x x x x a x x =+-=
-++-,则p(x)满足 (0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''========, 由(1)1f -=得 3222521
(1)(1)2(13)(1)1273108a a ?--?-++---=?=-
, 所以
223222
432
521
()()(3)2(3)273108
11332108544p x H x a x x x x x x x x x =+-=-+--=-++-+。
(2)余项具有如下结构
22()()()()(1)(3)r x f x p x k x x x x =-=+- 作辅助函数
22()()()()(1)(3)t f t p t k x t t t ?=--+-
则显然()t ?在点,1,0,3x -处有6个零点(其中0,3是二重零点),即 ()0,(1)0,(0)0,(0)0,(3)0,(3)0x ??????''=-=====, 不妨假设(1,0)x ∈-。
由罗尔定理,存在123(1,),(,0),(0,3)x x ξξξ∈-∈∈, 使得123()0,()0,()0?ξ?ξ?ξ'''===,
再注意到(0)0,(3)0??''==,即()t ?'有5个互异的零点12303ξξξ<<<< 再次由罗尔定理得,存在111223343(,),(,0),(0,),(,3)ηξξηξηξηξ∈∈∈∈, 使得1234()0,()0,()0,()0?η?η?η?η''''''''====
第三次应用罗尔定理得,存在112223334(,),(,),(,)ξηηξηηξηη∈∈∈ 使得123()0,()0,()0?ξ?ξ?ξ'''''''''===,
第四次应用罗尔定理得,存在112223(,),(,)μξξμξξ∈∈ 使得(4)(4)12()0,()0?μ?μ==,
第五次应用罗尔定理得,存在12(,)τμμ∈
使得(5)()0?τ= 注意到
(5)(5)(5)()()5!()()5!()t r t k x f t k x ?=-=-
(()()()r t f t p t =-中p(t)是4次函数,其5次导数为0)。 所以
(5)(5)(5)
()()()5!()=0()=5!
f f k x k x ξ?ττ=-?,
代入余项表达式,有
(5)22()
()()()(1)(3)5!
f r x f x p x x x x ξ=-=+-。
指出:
本题是非标准插值问题,比较简单的求解方法有:
①求插值问题的基本方法是待定系数法。以本题来说,有5个条件,可以确定一个4次的插值多项式,设为23301233y a a x a x a x a x =++++,将条件代入,建立一个5元的
线性方程组,求出各参数,就可以求出插值多项式。
②求插值问题的第二种方法是基函数法,即根据给定条件设定插值多项式的结构和各基函数的结构,根据条件确定基函数即可。具体方法与拉格朗日插值基函数构造和埃尔米特插值基函数构造相似。
③以标准插值为基础的方法是一种更简单的方法,本题中,首先利用4个条件构造一个埃尔米特插值,在此基础上设定所求插值多项式的一般形式,保证其满足埃尔米特插值条件,代入未利用条件解方程(组),求出其中的未知参数,即可求出插值多项式。
本题也可以先利用(1)(1)1,(0)(0)2,(3)(3)1p f p f p f -=-=-====构造一个2次插值多项式2()p x ,以此为基础构造4次插值多项式4()p x ,4()p x 的结构是 42()()()(1)(3)p x p x ax b x x x =+++-,
满足
(1)(1)1,(0)(0)2,(3)(3)1p f p f p f -=-=-====
再根据(0)(0)0,(3)(3)1p f p f ''''====列出两个线性方程组成的方程组,求出a 、b 两个参数,即可求出所求的插值多项式。 求插值函数余项()r x 的常用方法是:
()()()r x f x p x =-应具有如下形式(以本题为例) 22()()()()(1)(3)r x f x p x k x x x x =-=+- 作辅助函数
22()()()()(1)(3)t f t p t k x t t t ?=--+-
则()t ?在点,1,0,3x -处有6个零点(其中0,3是二重零点)。反复应用罗尔定理,直到至少有一个(4,4)τ∈-,使得(5)()0?τ=。此时即有
(5)(5)(5)
()()()5!()=0()=5!
f f k x k x ξ?ττ=-?
代入余项表达式即可求出。
7*、设f(x)在[-4,4]有连续的4阶导数,且 (0)2,(0)0,(3)1,(3)1f f f f ''====
试用两种方法构造三次埃尔米特插值多项式H(x),使其满足 (0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''========。 解一(待定系数法):
解:设230123()H x a a x a x a x =+++,则
2123()23H x a a x a x '=++,
由插值条件得 01
0123123
2(0)0(0)1(1)1(1)23H a H a H a a a a H a a a ==??'==??
==+++??'==++? 解之得012325
2,0,,327
a a a a ===-=,
所以32
52()2273
H x x x =-+。 解二(基函数法):
解:设300110011()()()()()()()()()H x f x x f x x f x x f x x ααββ''=+++,
因为线性拉格朗日插值基函数为100133()033
x x x x
l x x x ---=
==--,01100()303
x x x x
l x x x --=
==--, 由④得
200001
2
1001012
231
()[12()
]()1[12()]13[12(0)]033279227
x x x l x x x x x x x x x x x x x x x α=---??
-=-- ?
--??
-??
=-- ?
-??
-+=
同理
2
23
01101102()92()[1]27
x x x x x x x x x x x α??---=+=
?--?? 由⑤得
2
2
100013()()3x x x x x x x x x β??--??
=-= ? ?-????
2
32
011103()()9x x x x x x x x x β??--=-=
?-??
则
32
52()2273
H x x x =
-+。 8、设()(01)x f x e x =≤≤,试作一个二次多项式p(x),使其满足 (0)(0),(0)(0),(1)(1)p f p f p f ''===,并导出余项估计式。 解:设此二次式为2()p x a bx cx =++, 因为(),()x x f x e f x e '==, 所以,由已知条件
(0)(0)1,(0)(0)1,(1)(1)p f p f p f e ''====== 将其代入2(),()2p x a bx cx p x b cx '=++=+,得 11112a a b b a b c e c e ==????=?=????++==-??
所以,要求的二次多项式为 2()1(2)p x x e x =++-。
因为0是2重零点,1是1重零点,因此可以设余项具有如下形式: 2()()()()(0)(1)r x f x p x K x x x =-=--, 其中K(x)为待定函数。 固定x ,作辅助函数
2()()()(0)(1)t r t K x t t ?=--- 显然
(0)0,(0)0,()0,(1)0x ????'====, 不妨假设(0,1)x ∈。
由罗尔定理,存在12(0,),(,1)x x ξξ∈∈, 使得12()0,()0?ξ?ξ''==, 再注意到(0)0?'=
再次由罗尔定理得,存在11212(0,)(0,1),(,)(0,1)ηξηξξ∈?∈?, 使得12()0,()0?η?η''''==
再次应用罗尔定理,存在12(,)(0,1)ξηη∈? 使得
()0?ξ'''=。 注意到
()()3!()()3!()t r t K x f t K x ?'''''''''=-=-
(()()()r t f t p t =-中p(t)是2次函数,其3次导数为0)。 所以
()
()()3!()=0()=3!
f f K x K x ξ?ξξ'''''''''=-?, 代入余项表达式,有
2
2()()()()(0)(1)=(1)3!3!
f e r x f x p x x x x x ξξ'''=-=---。
指出:
石瑞民《数值计算》关于余项讨论很清楚。
9、给出sinx 在[0,π]上的等距结点函数表,用线性插值计算sinx 的近似值,使其
截断误差为41
102
-?,问该函数表的步长h 取多少才能满足要求?
解:设(0,1,)k x k =为等距结点,步长为h ,则1k k x x h +=+ 当1[,]k k x x x +∈时,作f(x)的线性插值
11111()()()k k
k k k k k k
x x x x L x f x f x x x x x ++++--=+--
则有
11()
()()()()2
k k f f x L x x x x x ξ+''-=--,
由此易知
12
11111()()max ()()(),[,]224
k k k k k k x x x h f x L x f x x x x x x x x +++≤≤''-≤--≤?∈
因此
2
1()()8
h f x L x -≤
由241
1082
h -≤?,得0.02h ≤。 指出:关于最大值的计算与12题相同。
10、求4()f x x =在区间[a,b]上的分段埃尔米特插值,并估计误差。
解:由分段三次埃尔米特插值多项式
30()[()()()()]n
i i i i i H x f x x f x x αβ='=+∑
则4()f x x =的分段埃尔米特插值为
30
430
()[()()()()]
[()4()]
n
i i i i i n
i i i i i H x f x x f x x x x x x αβαβ=='=+=+∑∑
其中
2
11112
1
1112
1112
112()[1],,0
2()()[1],,0,(),,0
()(),i i i i i i i i i i i i i i i i i i i i i i i i i i i i i x x x x x x x i x x x x x x x x x x x x i n
x x x x x x x x x x x i x x x x x x x x x x αβ----++++---++???--?+≤≤≠ ?--??????--?=+≤≤≠? ?--???
????
??--≤≤≠ ?-??
??-=- ?-??其他1,0,i x x i n
+?????
≤≤≠??
????
其他
其余项估计式为
444(4)
()max ()4!38438416
a x
b h h h r x f x ≤≤≤=?=。
解:这是第一类边界条件,要求解方程组 001
111222102012M g M g M g μλ?????? ??? ?= ??? ? ??? ??????? 其中
0101210101117.5 2.55107.5 2.52
0.66667
3
10.33333
h x x h x x h h h μλμ=-=-==-=-====+=-= 10
011
102111221
2122
226()0.5646
() 1.120026
() 1.608y y g y h h y y y y g h h h h y y g y h h ?-'=-=???--=
-=-?+??-'=-=??
将以上数据代入方程组 001111222102012M g M g M g μλ?????? ??? ?= ??? ? ??? ??????? 解之得
012
0.8073391.0506781.329334M M M =??
=-??=? 将获得的数据代入到
3322111
11()()()()()
6666i i i i i i i i i i i i i i i i i
x x x x M x x M x x S x M M y h y h h h h h ---------=++-+- 中,得
33
33
0.026911(7.5)0.035023( 2.5)0.127218(7.5) 2.275565( 2.5)()0.070045(10)0.088622(7.5) 3.237783(5.0) 1.446111(7.5)x x x x s x x x x x ??--?-+?-+?-?=?-?--?-+?-+?-??
12、设2()[,]f x C a b ∈(具有二阶连续导数),且f(a)=f(b)=0,证明:
21
max ()()max ()8
a x
b a x b f x b a f x ≤≤≤≤''≤- 证明:以a 、b 为节点进行插值,得
1()()()
1
()()()()()2!1
()()()()2!f x p x r x x b x a f a f b f x a x b a b b a f x a x b a b ξξξ=+--''=
++----''=--<< 因为()()x a x b --在1
()2
x a b =+处取得最大值,故
21
max ()max ()max ()()
21
()max ()
8
a x b
a x b
a x
b a x b f x f x x a x b b a f x ≤≤≤≤≤≤≤≤''≤
--''=-
13.给定数据表
x -2
-1
1
2 y
-0.1
0.1
0.4
0.9
1.6
用两种方法求其二次拟合曲线。 解一:
设所求的拟合函数为2y a bx cx =++, 则5
221[()]i i i i L a bx cx y ==++-∑。
对a 、b 、c 分别求偏导,并令偏导数等于0,得 5
21
5
21555
21
1
1
2[()]0[()]0
50
i i i i i i i i i i
i i i i L a bx cx y a a bx cx y a b x c x y =====?=++-=??++-=?++-=∑∑∑∑∑
5
21
5
215555
2
3
1
1
1
1
2[()]0[()]0
i i i i i i i i i i i i i i i i i i i L a bx cx y x b a bx cx y x a x b x c x x y ======?=++-=??++-=?++-=∑∑∑∑∑∑
5
221
5
2215555
23
4
21
1
1
1
2[()]0[()]0
i i i i i i i i i i i
i i i i i i i i L a bx cx y x c a bx cx y x a x b x c x x y ======?=++-=??++-=?++-=∑∑∑∑∑∑
将各数据点的数值代入,得方程组为
510 2.910 4.210347a c b a c +=??
=??+=?
解之得a=,b=0。42,c=,
所以数据点所反映的函数的近似关系为 20.40860.420.0857y x x =++
解二:设所求的拟合函数为2y a bx cx =++, 将数据代入方程得
240.10.10.4
0.924 1.6
a b c a b c a a b c a b c -+=-??-+=??
=??++=?++=?? 方程组的系数矩阵和右端向量为
1240.11110.1,1000.41110.9124 1.6A B --???? ? ?- ? ? ? ?== ? ? ? ? ? ?????
因为
124111115010111210120100,
10041014100341111240.111111 2.90.121012 4.20.44101470.91.6T T A A A B -?? ?
-???? ? ? ? ?=--= ? ? ? ? ????
? ? ???
-?? ?
???? ? ? ? ?=--= ? ?
? ? ?????
? ???
所以
5010 2.90100 4.2100347a b c ?????? ??? ?= ??? ? ??? ???????
解之得a=,b=0。42,c=,
所以数据点所反映的函数的近似关系为 20.40860.420.0857y x x =++ 14、已知试验数据
x 19
25
31
38
44 y
19.0
32.3
49.0
73.3
97.8
用最小二乘法求形如2y a bx =+的经验公式,并计算均方误差。 解:设2y a bx =+ 则
5
221[()]i i i L a bx y ==+-∑
对a 、b 分别求偏导,并令偏导数等于0,得
5
21
5
2155
21
1
2[()]0[()]0
50
i i i i i i i i i i L a bx y a a bx y a b x y ====?=+-=??+-=?+-=∑∑∑∑
5
221
5
221555
24
21
1
1
2[()]0[()]0
i i i i i i i i i
i i i i i i L a bx y x b a bx y x a x b x x y =====?=+-=??+-=?+-=∑∑∑∑∑
将数据代入得
22222222224444422
2225(1925313844)(19.032.349.073.397.8)0(1925313844)(1925313844)(1919.02532.33149.03873.34497.8)0a b a b ?+?++++-++++=??+++++?++++-?+???+?+?+?=?
化简得
55327271.40
53277277699369321.50
a b a b +-=??
+-=? 第二个方程减去第一个方程乘以1065进一步化简得 55327271.40
2160444480279.50a b a b +-=??
+-=? 解之得 1.01
0.05
a b =??
=? 则x 与y 的函数关系是 y=+。
此时,平方逼近误差为
5
221[()]0.017i i i L a bx y ==+-=∑
所以,
0.13=。 指出:
均方误差实际上就是按最小二乘法则确定的残差。 15
求运动方程。
解:设运动方程为s =a+bt 则
6
666
21
1
1
1
14.7,53.63,280,1078i i
i i i i i i i t
t s t s ========∑∑∑∑
将上述数据代入方程组
66
11
666
21
116i i i i i i i i i i i a b t s a t b t t s =====?
+=????+=??∑∑∑∑∑ 得方程组
614.7280
14.753.631078a b a b +=??
+=?
解之得
7.8550478
22.25376
a b =-??
=? 所以,7.855047822.25376s t =-+。 指出:
利用统计型计算器,有关中间数据可以简单求出。
用最小二乘法求y=f(t)。
解:描草图,观察草图可以发现,该组数据分布近似于指数函数曲线,而且随着t 的增大,y 的增速放缓,故设
b t
y ae =。
两边取对数,得
1
ln ln y a b t =+,
令1
ln ,,ln y z s a c t
===,
则拟合函数转化为线性拟合关系z c bs =+。
11
11
2110.6039755,0.06232136i
i i i s
s ====∑∑
1111
1
1
13.639649,0.5303303i
i i i i z
s z ====∑∑。
将上述数据代入
1111
11
111111
211111i i i i i i i i i i i c b s z c s b s s z =====?
+=????+=??∑∑∑∑∑ 得
110.603975513.639649
0.60397550.062321360.5303303c b c b +=??
+=? 解之得
7.4961692, 1.6515592 5.2151048b c a =-=?= 所以
7.4961692
5.2151048t
y e -
=。 指出:
(1)T=0,该拟合函数不适用。
(2)专业的变化规律(经验函数)应当由专业人员给出。仅仅从有限数据的草图得出的规律可能不具普遍性。 17
用最小二乘法求形如y ae =的经验公式。
解:对bx y ae =两边取对数,得 ln ln y a bx =+,
令01ln ,ln ,y Y a a b a ===, 则
01Y a a x =+,
代入数据,建立方程组为 0101517.319.97968815
17.364.2368.55117703a a a a +=??
+=? 解之得
00114.45380
85.95290.1323290.132329a a a e a b a =??==??
?=-==-?? 所以
0.13232985.9529x y e -=。
18、用最小二乘法求方程组
241135326214
x y x y x y x y +=??-=?
?
+=??+=? 的近似解。
分析:这是方程个数多于未知数个数的超定方程组,是矛盾方程组,用最小二乘法求解。
解:设方程组中各个方程的一般形式为i i i a x b y c +=,则
4
21[()]i i i i L a x b y c ==+-∑
对x 、y 分别求偏导,并令偏导数等于0,得
4
1
4
1444
21
1
1
2[()]0[()]0
i i i i i i i i i i i i i i i i i i L a x b y c a x a x b y c a x a y a b a c =====?=+-=??+-=?+-=∑∑∑∑∑
4
1
4
1444
21
1
1
2[()]0[()]0
i i i i i i i i i i i i i
i i i i i L a x b y c b y a x b y c b x a b y b b c =====?=+-=??+-=?+-=∑∑∑∑∑
将数据代入得 153510
349690x y x y --=??
-+-=? 解之得 3.727
1.636x y =??
=?
它有形如()p x a bx =+的拟合函数,试求本问题的最小二乘解。
解:令1
y z
=,则拟合函数变形为
z a bx =+,原拟合问题转化为线性拟合问题。 则8
21[()]i i i L a bx y ==+-∑。
对a 、b 分别求偏导,并令偏导数等于0,得 8
1
8
188
1
1
2[()]0[()]0
80
i i i i i i i i i i L a bx y a a bx y a b x y ====?=+-=??+-=?+-=∑∑∑∑
8
1
8
1888
2
1
1
1
2[()]0[()]0
i i i i i i i i i i i i i i i L a bx y x b a bx y x a x b x x y =====?=+-=??+-=?+-=∑∑∑∑∑
将数据代入,得
88
11
888
211180836419.90836419.90362042479.4848737.900i i i i i i i i i i i a b x y a b a b a b b a x b x x y =====?
+-=?+-=+-=????????+--+=?
??+-=??∑∑∑∑∑ 解之得 520.58
104.02
a b =??
=-? 所以,所求的拟合函数为
1
()520.58104.02p x x
=
-。 20、在平面上给出三个点,它们的坐标是123(1,1),(2,0),(1.5,3)T T T x x x ===,每个点
对应一个函数值1231.8, 2.6, 3.1z z z ===,找出一个通过这三个点的平面。 解:这实际上是求过三个点(1,1,1.8),(2,0,2.6),(1.5,3,3.1)的平面方程。 由解析几何知识可知,平面的三点式方程为 1112223
2
31
1011
x y z x y z x y z x y z =
将三点坐标代入,解此方程就可求出所求平面方程。 (以下从略)
补充题(一)
1、求次数不超过2和3的多项式p 2(x)和p 3(x)。使得
p 2(0)=p 3(0)=0,p 2(1)=p 3(1)=1,p 2(2)=p 3(2)=8,p 3(3)=27。 解一:设二次多项式为p 2(x)=a 0+a 1x+a 2x 2 ,则有 20122
0122012000111228a a a a a a a a a ?+?+?=?+?+?=??+?+?=?
解之得,0120,2,3a a a ==-=。所以
22()23p x x x =-+。
设三次多项式为p 3(x)=a 0+a 1x+a 2x 2+a 3x 3 ,则有 23012323
012323
01232301
2300001111
222833327
a a a a a a a a a a a a a a a a ?+?+?+?=?+?+?+?=??+?+?+?=??+?+?+?=? 解之得,01230,0,0,1a a a a ====。所以
32()p x x =。
解二:由题6,可以直接利用插值多项式公式求出所要求的多项式来。 解三:在学习了差商和差分后,也可以利用牛顿插值公式或等距节点插值
公式求出所求多项式。
对f(x)在0,1,2,3处求差商得
所以,p 2(x)=p 2(0)+1×(x -0)+3×(x -0)(x-1)=3x 2-2x ,
p 3(x)=p 3(0)+1×(x -0)+3×(x-0)(x-1)+1×(x -0)(x-1)(x-2)=x 3。
2、已知函数f(x)在节点-1,0,1处的值分别是,,,用待定系数法和插值基函数法两种方法求出拉格朗日插值。
解1:设所求的多项式为
22012()p x a a x a x =++,把已知条件代入得
20122
0122012(1)(1)0.3679(0)(0) 1.000
(1)(1) 2.7182a a a a a a a a a ?+?-+?-=?+?+?=??+?+?=?
解之得
0121, 1.751,0.5431a a a === 所以
22()1 1.17510.5431p x x x =++。 解2:由插值基函数公式
00()
()()
n
k
k k i i n
i
k
k k i
x x l x x x =≠=≠-=
-∏∏
0(0)(1)(1)
()(10)(11)2x x x x l x ---=
=----
1[(1)](1)()(1)(1)[0(1)](01)
x x l x x x ---==-+----
2[(1)](0)(1)()[1(1)](10)2
x x x x l x ---+==---
代入插值公式得
2012()0.3679() 1.000() 2.7182()p x l x l x l x =++ 即
22()1 1.17510.5431p x x x =++。
3、设f(x)=x 4,试利用拉格朗日插值余项定理写出以-1,0,1,2为插值节点的三次插值多项式。
解:记三次插值多项式为p(x),由插值余项定理
(1)(4)1
()()()()
(1)!
1
()(1)(0)(1)(2)4!
(1)(0)(1)(2)
n f x p x f x n f x x x x x x x x ξπξ+-=+=+---=+---
所以,
4232()()(1)(0)(1)(2)(1)(2)
22p x f x x x x x x x x x x x x =-+---=---=+-
思考:
用插值多项式公式直接求插值多项式与本题求出的多项式比较一下。 4、已知= 567,= 487,= 274,用抛物线插值计算。 解:= 374。
5、设l k (x)(k=0,1,2,…,n)是n+1个互异节点x 0,x 1,x 3,…x n 上的n 次基本插值多项式,证明下面的恒等式成立
()(0,1,2,
,)n
m
m k k
k x
l x x m n ===∑
证明: 由拉格朗日插值定理,以x 0,x 1,x 2,…x n 为插值节点,对y=f(x)=x m 作n 次插值,插值多项式为
0()()n
n i i i p x l x y ==∑,
而 y i =x i m , 所以
()()()n
n
m n i i i i i i p x l x y l x x ====∑∑
同时,插值余项
(1)(1)11
()()()()()()0(1)!(1)!
n m n m n r x x p x f x x x n n ξξππ++=-=
==++
所以
()n
m
m i i
i l x x
x ==∑