- 7
- 0
- 约8.99千字
- 约 25页
- 2022-12-11 发布于广西
- 举报
数值分析
Mr — z±z± 川一早
2?当X=l,—1,2时,f(x)=O, -3,4,求f(x)的二次插值多项式。
解:X _l,x
解:
X _l,x
0 —
1, X _2,
2 —
f(Xo )= 0, f (xj = -3, f (x2)= 4;lo(X
f(Xo )= 0, f (xj = -3, f (x2)= 4;
lo(X)亠二-1 (x 1) (x-2)
(x 2
o -X/Xo J2)
(x -Xo) (x -X2) 1
li(x) 0 (x-l)(X-2)
6
(x i yo)(xi \)
L(x)
(x一X0)(X一X,) 1
0 (x-l) (x 1)
(x2 —x。) (x2 - Xi) 3
则二次拉格朗日插值多项式为
L y 1
2(X)=,- k
kzO
(X)
二-31° (x)
4L (x)
4
4
= (x_l) (x一2) (x-l) (x 1)
2
2 3
xx~
2
6?设 Xj, j =0,b
n为互异节点,求证:
n
(1 ) 7 x: lj (x) =xx
j=0
(k =0, 1 川,n);
n
2
()7 (Xj -x)kl j(x)三
j £
证明
⑴令 f (x) =xk
若插值节点为Xj,
n
j =0, 1, 1(, n,则函数f (x)的n次插值多项式为Ln(x?l xklj(x)。
j=0
f (n 十)(?)插值余项为 Rn(X)二 f (X) Ln(X)
nl (X)
(n +1) !
? f(nl)Rn (X)二0n二瓦
? f(nl)
Rn (X)二0
n
二瓦 xklj(x) =xk (k =0, 1 川,n);
j : o
n
⑵、(Xj ~x)klj (x)
卫
n n
二為(C?Xj (―x) k-L) lj (x)
j =0 i =0
n n
i k i i
二為 Ck( -x) (、Xjl j (x))
i =0 j=0
又70 i n 由上题结论可知
n
.原式二厂 Ck (-x) k_Lx,
=(X -x)k
二 0
-得证。
7 设 f (x) c2 la, b 1 且 f
=f (b)二 0,求
(a)
证:
max f (x)兰一 (b~a) max
f
a $31 小 -症
(x)?
解:令x八a,x~,以此为插值节点,则线性插值多项式为
二二 f (a)
f(b)-
又 T f Q)二 f (b)二 0
Li(x)二 0
1
插值余项为 R(x)二 f (x) - L, (x) f (x) (x - Xq)(X - xj 1
f (x) = 2 f (x) (x -Xg) (X -xj
又;(X —Xo) (X —Xx)
声乞
声
乞 2 ((X -Xo)
-.2
(Xi - X)]
二 g(Xx-Xo )2
4
二](b-a)2
4
1 2
f (x) _ (r) f ” (x
二 max —8 醴 ).
在-4乞x岂4上给出f (x) =ex的等距节点函数表,若用二次插值求截断误 h的近似值,要使 差不超过10:问使用函数表的步长h应取多少? 解:若插值节点为Xj丄Xj和X”则分段二次插值多项式的插值余项为
1 X
R(x) =3 f ” (-) (X —XiJ ) (X—X) (x —Xy) 3!
1
二 R(X)M6(x—XiJ ) (x—x) (x—xy)max f (X)
设步长为h,即Xi」二Xj-h, x i二人h
一eh3.
27
若截断误差不超过10$,则
R2(X)〈10
山 eh^lO 11
27
.仁 0. 0065.
9?若 yn =2:求「yn 公? yn.,
解:根据向前差分算子和中心差分算子的定义进行求解。
n
yn =2
:Yn = (E -1) yn
1
1
1
1
4
4 £
ElJ y
ElJ yn
J
?y
j
?g
j=0
4
汽X-1)J
j=0
4
二s (-i)J
:n
-(2T)yn
二 yn
二 2n
11、. yn = (E2 —E J yn
1
二(E)(ET)Sn?%
二 ynd
二 2no
16. f (x) +X‘+3X+1,求F [芳门川公]及F [芳门川/]。
解:T f (x) = x1 x1 3x 1
若 x = 2」,i 二 0, 1J |,8
/)(亡、
则f 。,为,丨1(凶丄一一
n!
? f (Xo,Xb 1(, X7 丄小丄二 1
7! 7!
f⑻心
f (Xo, Xi, 1 (, xs i =— =0
8!
19 .求一个次数不高于4次的多项式P ( X ),使它满足
P(0)二 P (0) =0,P(l) = P (1) = 0,P(2) =0
解法一:利用埃米尔特插值可得到次数不高于 4的多项式
x 0, x 1
0 — 1
原创力文档

文档评论(0)