习题5.3 Tay lor公式和插值多项式 1.由 Lagrange中值定理知 00,得到 8.+(x) 1) 再由(m(x)≠0,两边消去r(x),即得到lmb=-1 3.设f(x)=Vx,取结点为x=1、1728、2744,求f(x)的二次插值多项 式P2(x)及其余项的表达式,并计算p(2)(v2=12599210)
习 题 5.3 Taylor 公式和插值多项式 1．由 Lagrange 中值定理知 x x x x 1 ( ) ln(1 ) +θ + = ，0 < θ (x) < 1， 证明：lim ( ) 1/ 2 0 = → x x θ 。 证 由 ln(1 ) ln(1 ) ( ) x x x x x + − + θ = ，取极限即得到 2 2 0 0 0 0 ln(1 ) ln(1 ) lim ( ) lim lim lim ln(1 ) ln(1 ) x x x x x x x x x x x x x x θ → → → → x − + − + = ⋅ = ⋅ + + 0 0 0 1 1 1 1 1 lim lim (lim ) 1 2 2 1 (1 1 x x x x x x x → → → − + = ⋅ = ⋅ + + 1 ) 2 = 。 2．设 1 2 ( 1 ) ( ) ( ) '( ) "( ) ( ) 2! ! n n f x h f x f x h f x h f x h h n + = + + +"+ +θ ，(0 < θ < 1)， 且 f (n+1) (x) ≠ 0，证明： 1 1 lim 0 + = h→ n θ 。 证 n n f x h h n f x h f x f x h f x h ( ) ! 1 "( ) 2! 1 ( ) ( ) '( ) 2 ( ) + = + + +"+ +θ ( ) ( ) ( 1)! 1 ( ) ! 1 "( ) 2! 1 ( ) '( ) 2 ( ) ( +1) +1 +1 + + = + + + + + n n n n n f x h h n f x h n f x f x h f x h " D ， 于是 ( ) (1) 1 ( ) ( ) 1 ( 1) ( ) ( ) + D + = + − ⋅ + f x h n f x h f x n n n θ θ θ 。 令h → 0，得到 ( 1) ( 1) 0 1 lim ( ) ( ) 1 n n h f x f n θ + + → ⋅ = x + ， 再由 f (n+1) ( x) ≠ 0 ，两边消去 ( 1) ( ) n f x + ，即得到 1 1 lim 0 + = h→ n θ 。 3．设 f x( ) = x 3 ，取结点为 ，求 的二次插值多项 式 及其余项的表达式，并计算 （ x = 1、 、 1.728 2.744 f x( ) p x 2 ( ) p2 (2) 2 12599210 3 = . "）。 116
解f()=1,f(1.728)=1.2,(2.744)=14,由 Lagrange插值公式 f(x)≈P2(x) (x-1.728(x-2744) +1.2 (x-1)(x-2.744) +14 (x-1)(x-1.728) (1-1.728)(1-2.744) (1.728-1)(1.728-2.744)(2.744-1)(2.744-1.728) 0.7876(x-1.728)(x-2.744)-1.6224(x-1)(x-2.744)+0.7901(x-1)(x-1.728) =-004465x2+0.3965x+0.6481。 f()=0x÷,余项()=-5(-1x-18x-214) P2(2)≈1.2626。 4.设f(x)=2,取结点为x=-1、0、1,求f(x)的二次插值多项式P2(x) 及其余项的表达式,并计算P(3/。请与上题的计算结果相比较并分 析产生差异的原因 解f(-1)=0.5,f(0)=1,f()=2,由 Lagrange插值公式 f(x)≈P2(x) =05.(x-0(x-1)1,(x+1)(x-1),(x+1)(x-0) (-1-0)(-1-1)(0+1)(0-1)(1+1)(1-0) =0.25x(x-1)-(x-1)(x+1)+(x+1)x 0.25x2+0.75x+1。 f(x)=1n2,余项r(x)=2(x+)x(x-1 6 P2()≈1.2778。 与上题相比,本题误差较大的原因是2不在所取的三点 x=-1、0、1之间,而上题2在所取的三点x=1、1728、2744之间,因 而误差较小 5.设f(x)在若干个测量点处的函数值如下 14|1.72.3|31 117
解 f f (1) = = 1, (1.728) 1.2, f (2.744) =1.4，由 Lagrange 插值公式 2 ( ) ( ) ( 1.728)( 2.744) ( 1)( 2.744) ( 1)( 1.728) 1 1.2 1.4 (1 1.728)(1 2.744) (1.728 1)(1.728 2.744) (2.744 1)(2.744 1.728) 0.7876( 1.728)( 2.744) 1.6224( 1)( 2.744) 0.7901( 1)( 1.728) f x p x x x x x x x x x x x x x ≈ − − − − − − = ⋅ + ⋅ + ⋅ − − − − − − ≈ − − − − − + − − = − 2 0.04465x +0.3965x+0.6481。 8 3 10 ( ) 27 f x x − ′′′ = , 余项 2 8 3 5 ( ) ( 1)( 1.728)( 2.744) 81 r x x x x ξ = − − − 。 2 p (2) ≈1.2626。 4．设 f x( ) = 2x ，取结点为 x = −1 0 、 、1，求 f (x)的二次插值多项式 p x 2 ( ) 及其余项的表达式，并计算 ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ 3 1 p2 。请与上题的计算结果相比较并分 析产生差异的原因。 解 f f ( 1− =) 0.5, (0) =1, f (1) = 2 ，由 Lagrange 插值公式 2 2 ( ) ( ) ( 0)( 1) ( 1)( 1) ( 1)( 0) 0.5 1 2 ( 1 0)( 1 1) (0 1)(0 1) (1 1)(1 0) 0.25 ( 1) ( 1)( 1) ( 1) 0.25 +0.75 +1 f x p x x x x x x x x x x x x x x x ≈ − − + − + − = ⋅ + ⋅ + ⋅ − − − − + − + − = − − − + + + = 。 3 ( ) ln 2 2x f x ′′′ = ⋅ , 余项 3 ln 2 ( ) 2 ( 1) ( 1) 6 nr x x x x ξ = + − 。 2 1 ( ) 1.2778 3 p ≈ 。 与上题相比，本题误差较大的原因是 2 不在所取的三点 之间，而上题 2 在所取的三点 之间，因 而误差较小。 x = −1 0 、 、1 x = 1、 、 1.728 2.744 5．设 f (x)在若干个测量点处的函数值如下： x 1.4 1.7 2.3 3.1 117
0)|65384 试求f(28)的近似值。 解由 Lagrange插值公式 f(x)≈p3(x) (x-1.7)(x-2.3)(x-3.1) (x-1.4)(x-2.3)x-3.1) (14-1.7)14-2.3)(1.4-3.1)(1.7-141.7-23)(1.7-3.1) (x-1.4)(x-1.7)(x-3.1) +36.(x-14)x-1.7)(x-23) (23-14)(2.3-1.7)(23-3.1)(3.1-14)3.1-1.7(31-2.3) f(2.8)≈P3(28)≈36647。 6.若h是小量,问如何选取常数a、b、c,才能使得 (x+h)+b(x)+c/(x-h)与f"(x)近似的阶最高? 解刂(x+h)+b(x)+(x-h) =d(x)+f(x)h+1r(x)2]+b(x)+c(x)-f(x)h+"xh2]+o(H2) (a+b+cf(x)+(a-cf'(xh+(a+c)f"(x)h+o(h), b+C=0 得到方程组{a-c=0,解之得到a=c=1.b= 7.将插值条件取为n+1个结点上的函数值和一阶导数值,即 pn(x)满足 P,(x)=f(x) i=0,1,2,…,n p(x2)=f'(x) 的插值多项式称为 Hermite插值多项式,在微分方程数值求解等研究 领域中具有重要作用。它可以取为 p2(x)=∑[xk(x)+(x)(x)] 这里,{qg(x)q"(x)是满足条件
f (x) 65 58 44 36 试求 f ( . 2 8)的近似值。 解 由 Lagrange 插值公式 3 ( ) ( ) ( 1.7)( 2.3)( 3.1) ( 1.4)( 2.3)( 3.1) 65 58 (1.4 1.7)(1.4 2.3)(1.4 3.1) (1.7 1.4)(1.7 2.3)(1.7 3.1) ( 1.4)( 1.7)( 3.1) ( 1.4)( 1.7)( 2.3) 44 3.6 (2.3 1.4)(2.3 1.7)(2.3 3.1) (3.1 1.4 f x p x x x x x x x x x x x x x ≈ − − − − − − = ⋅ + ⋅ − − − − − − − − − − − − + ⋅ + ⋅ − − − − , )(3.1− − 1.7)(3.1 2.3) 3 f p (2.8) ≈ (2.8) ≈ 36.647。 ) 6．若h是小量，问如何选取常数a、 、b c，才能使得 af ( ) x + + h bf (x) + cf (x − h 与 f ′′(x)近似的阶最高？ 解 2 2 2 2 ( ) ( ) ( ) 1 1 [ ( ) '( ) ''( ) ] ( ) [ ( ) '( ) ''( ) ] ( ) 2 2 1 ( ) ( ) ( ) '( ) ( ) ''( ) ( ) 2 af x h bf x cf x h a f x f x h f x h bf x c f x f x h f x h o h a b c f x a c f x h a c f x h o h + + + − = + + + + − + + = + + + − + + + ， 2 得到方程组 0 0 2 abc a c a c ⎧ + + = ⎪ ⎨ − = ⎪ ⎩ + = ，解之得到a c = =1,b = −2。 n 7．将插值条件取为n + 1个结点上的函数值和一阶导数值，即 p (x)满足 ， p x f x p x f x n i i n i i ( ) ( ) ( ) ( ) = ′ = ′ ⎧ ⎨ ⎩ i n = 0,1, 2,", 的插值多项式称为 Hermite 插值多项式，在微分方程数值求解等研究 领域中具有重要作用。它可以取为 (0) (1) 0 ( ) ( ) ( ) ( ) ( ) n n k k k k k p x f x q x f x q = = + ⎡ ′ x ⎤ ∑⎣ ⎦, 这里，{ ( q x k k (0) ),q(1) (x)} n k=0是满足条件 118
q0(x)=6k,[q](x)=0,1,k=0,1,2…,n 和 qe(x)=0,[q]x)=k, i,k=0,1,2…,n 的基函数。试仿照 Lagrange插值多项式的情况构造{qg(x)q"(x)=0 解显然当i≠k时,q(x)=[q(x)=0,q(x)=1,[q(x)=0,设 "(xp-a(-x),由Yx)=2.-=0解出,得 到 )(x-x 同理可得到 q(x)
(0) (0) ( ) , [ ]'( ) 0, k i ik k i q x = = δ q x i k , 0, = 1, 2,", n 和 (1) (1) ( ) 0, [ ]'( ) , k i k i ik q x = = q x δ i k , 0, = 1, 2,", n 的基函数。试仿照 Lagrange 插值多项式的情况构造{ ( q x k k (0) ),q(1) (x)} n k=0。 解 显然当i ≠ k 时， (0) (0) ( ) [ ] ( ) 0, k i k i q x = q ′ x = (0) (0) ( ) 1,[ ]'( ) 0 k k k k q x = q x = ，设 (0) 2 0 ( ) ( ) [1 ( )] n i k k i k i i k x x q x c x x = x x ≠ ⎡ ⎤ ⎢ ⎥ − = − ⎢ ⎥ − ⎢ ⎥ ⎣ ⎦ ∏ − ，由 (0) 0 2 [ ]'( ) n k k i k i i k q x c = x x ≠ = − = 0 − ∑ 解出 ，得 到 c (0) 2 0 0 2 ( ) ( ) [1 ( )( )], 0,1,2, , n n i k k i k i i k i i k i k x x q x x x k n = x x = x x ≠ ≠ ⎡ ⎤ ⎢ ⎥ − = − − = ⎢ ⎥ − − ⎢ ⎥ ⎣ ⎦ ∏ ∑ " ； 同理可得到 (1) 2 0 ( ) ( ) ( ) 0,1,2, , n i k k i k i i k x x q x x x k n = x x ≠ ⎡ ⎤ ⎢ ⎥ − = − = ⎢ ⎥ − ⎢ ⎥ ⎣ ⎦ ∏ " 。 119
