【正文】 x? ? ? ?115 1 5 8 5 1 5( ) d ( 0 )9 5 9 9 5f x x f f f?? ? ? ?? ? ? ?? ? ? ?? ? ? ?? ? ? ??311P ( ) ( 5 3 )2n x x x? ??07:49:44 Numerical Analysis 62 更多 GL 公式 當(dāng) n 3 時(shí)，可用數(shù)值方法計(jì)算 Pn+1(x) 的零點(diǎn) n 節(jié)點(diǎn)個(gè)數(shù) Gauss點(diǎn) Gauss系數(shù) 0 1 1 2 ? 2 3 ? 3 4 ? ? 4 5 ? ? 5 6 ? ? ? 07:49:44 Numerical Analysis 63 余項(xiàng) ? 余項(xiàng) : ? ?? ?422( 2 2 )32 ( 1 ) ! ( )( 2 3 ) ( 2 2 ) !nnn fnn??????? ? ? (1, 1) ( 2 2 )21()[ ] ( ) ( ) d( 2 2 ) !n bnafR f x x xn? ????? ? ?07:49:44 Numerical Analysis 64 GL 公式 一般區(qū)間上的 GL 求積公式 ? 變量代換 22b a b axt????11 0( ) d ( ) d ( )2nbiiaibaf x x g t t A g t? ???? ???()22b a b ag t f t??????????07:49:44 Numerical Analysis 65 舉例 例： 用四點(diǎn) GL公式 (n=3) 計(jì)算定積分 20 co s ( ) dx x x??解： 令 44xt????? ?2 2( ) 1 co s ( 1 )1 6 4g t t t??? ? ?? ?20 . 5 1 2201co s ( ) d 1 co s ( 1 ) d4 1 6 4x x x t t t? ? ? ??? ? ???[ 0. 34 79 ( 0. 86 11 ) 0. 65 21 ( 0. 34 00 )4 0. 65 21 ( 0. 34 00 ) 0. 34 79 ( 0. 86 11 ) ]gggg?? ? ? ????[ ] = 0. 46 74 01 1 00 27 23 4If07:49:44 Numerical Analysis 66 GC 公式 ? 積分區(qū)間 : [1, 1] ? 權(quán)函數(shù) : GaussChebyshev 求積公式 Gauss 點(diǎn) = Chebyshev 多項(xiàng)式 Tn+1(x) 的零點(diǎn) 211 0((11) d )niiif x x A f xx? ??? ??? GC 求積公式 : 21()1xx? ??07:49:44 Numerical Analysis 67 GC 公式 ? Tn+1(x) 的零點(diǎn) 21co s22iixn ????? ????? (i = 0, 1, … , n) ? Gauss 系數(shù) 1iA n??? (i = 0, 1, … , n) 211 0( ) 11d ( )1niif x x f xnx?? ???? ??? GC 求積公式 : ? 余項(xiàng) : [] ( 2 2 )22 2 ( )2 ( 2 2 ) ! nnR f fn? ???? ?? ? (1, 1) 07:49:44 Numerical Analysis 68 簡(jiǎn)單 GC 公式 ? n = 0 1 2 1 / 21 ( 1 ) ( ) d ( 0 )x f x x f??? ???? n = 1 ? ? ? ?1 2 1 / 21 ( 1 ) ( ) d 2 2 2 22x f x x f f??? ??? ? ? ????? n = 2 兩點(diǎn) GC 公式 三點(diǎn) GC 公式 ? ? ? ? ? ?1 2 1 / 21 ( 1 ) ( ) d 3 2 0 3 23x f x x f f f??? ??? ? ? ? ????一般區(qū)間上的 GC 求積公式 ? 變量代換 22b a b axt????07:49:44 Numerical Analysis 69 舉例 例： 用五點(diǎn) GC公式計(jì)算定積分 211 d1xxex? ??解： 直接代公式可得 ?41 021 21 d co1s5 2 2ix ixf nex? ?? ??? ???????? ?????? ??誤差估計(jì) [] ?? ?? ? ( 2 2 )22 2 ( )2 ( 2 2 ) ! nnR f fn? ??? ? ?? 910 2 4 . 6 1 02 1 0 ! e?07:49:44 Numerical Analysis 70 GaussLaguerre 公式 (自學(xué) ) ? 積分區(qū)間 : [0, ?] ? 權(quán)函數(shù) : GaussLaguerre 求積公式 Gauss 點(diǎn) = Laguerre 多項(xiàng)式 Ln+1(x) 的零點(diǎn) 0 0( ) d ( )xniiif x x A f xe??? ? ??? GaussLaguerre 求積公式 : () xxe? ?? d( ) ( )dnx n xn nL x e x ex??07:49:44 Numerical Analysis 71 幾點(diǎn)注記 ? Gauss 型求積公式的優(yōu)點(diǎn) ? 計(jì)算精度高 ? 可計(jì)算無窮區(qū)間上的積分和奇異積分 ? Gauss 型求積公式的缺點(diǎn) ? 需計(jì)算 Gauss 點(diǎn)和 Gauss 系數(shù) ? 增加節(jié)點(diǎn)時(shí)需重新計(jì)算 ? 復(fù)合 Gauss求積公式 ? 將積分區(qū)間分隔成若干小區(qū)間 ? 在每個(gè)小區(qū)間上使用 Gauss 求積公式 07:49:44 Numerical Analysis 72 本講內(nèi)容 ? 基本思想 ? 計(jì)算方法 ? 多重積分 ? 數(shù)值微分 ? 問題描述 ? 計(jì)算方法 07:49:44 Numerical Analysis 73 多重積分 ( , ) df x y v???? 計(jì)算二重積分 ? 基本思想： 化成累次積分 ( , ) d ( , ) d dbdacf x y v f x y y x???? ? ?[ , ] [ , ]a b c d? ? ?2R??07:49:44 Numerical Analysis 74 舉例 例： 用兩點(diǎn) Gauss求積公式計(jì)算二重積分 222dx y v???? [ 1 , 1 ] [ 1 , 1 ]? ? ? ? ?令 ，可得 令 ，可得 解： 112 2 2 2112 d 2 d dx y v x y y x??? ? ? ??? ? ?1 2 2 213 3 42 d , , 23 3 3x y y f x f x x?? ? ? ?? ? ? ? ? ?? ? ? ?? ? ? ?? ? ? ??22( , ) 2f x y x y??2 4( ) 23g x x??1 1 1221 1 1332 d d ( ) d 433x y y x g x x g g? ? ?? ? ? ?? ? ? ? ? ?? ? ? ?? ? ? ?? ? ? ?? ? ?07:49:44 Numerical Analysis 75 復(fù)合公式 ? 復(fù)合公式 ? 對(duì)變量分別積分時(shí)，采用復(fù)合求積公式 07:49:44 Numerical Analysis 76 數(shù)值微分（自學(xué) ） ? 數(shù)值微分 已知 f(x) 在節(jié)點(diǎn) a ? x0 x1 xn ? b 上的函數(shù)值， 對(duì)于 [a, b] 中的任意一點(diǎn) ， 如何計(jì)算其導(dǎo)數(shù) ? 插值型求導(dǎo)公式 ? 構(gòu)造出 f (x) 的插值多項(xiàng)式 pn(x) ? 用 pn(x) 的導(dǎo)數(shù)來近似 f (x) 的導(dǎo)數(shù) ? 外推算法 07:49:44 Numerical Analysis 77 插值型求導(dǎo)公式 ? 插值型求導(dǎo)公式的余項(xiàng) ? ?( 1)( 1)00() 1( ) ( ) ( ) ( ) ( )( 1 ) ! ( 1 ) !n nnnxn j j xjj39。 39。ff 39。 x P 39。 x x x x x fnn? ?? ?????? ? ? ? ?????????39。( ) 39。 ( )nf x P x?( 1)0()( ) ( ) ( )( 1 ) !n nxi n i i jjjiff 39。 x P 39。 x x xn????? ? ?? ?? 在節(jié)點(diǎn)處的余項(xiàng) ? 插值型求導(dǎo)公式 07:49:44 Numerical Analysis 78 兩點(diǎn)公式 ? 兩點(diǎn)公式 ? 節(jié)點(diǎn) x0 , x1 ， 步長(zhǎng) h = x1 x0 011 0 10 1 1 0( ) ( ) ( )xxxxP x f x f xx x x x??????1 0 0 1( ) ( ) ( ) ( )x x f x x x f xh? ? ? ??? ?00 011( ) ( ()) 2)(f39。 x f x f h f39。 39。xh ????? ?10 111( ) ( ()) 2)(f39。 x f x f h f39。 39。xh ????07:49:44 Numerical Analysis 79 三點(diǎn)公式 ? 三點(diǎn)等距公式 ? 步長(zhǎng) h ， 節(jié)點(diǎn) xi = x0 + ih ， i = 0, 1, 2 12200 1 0 20 2 0 1121 0 1 2 2 0 2 1( )( )( ) ( )( )( )( )( ) ( )( ) ( ) ( )( )( ) ( )( )x x x xP x f xx x x xx x x x x x x xf x f xx x x x x x x x?????? ? ? ???? ? ? ?? ?0 0 1 2 ( 3 ) 021( ) 3 ( ) 4 ( ) ()( 3) 2f39。 x f x f x f hx fh ?? ? ? ? ?? ? 2 (2 )1 3 101( ) ( ) ( ) ()2 6 hf39。 x f x f xh f ??? ??? ? 202 321 () 21( ) ( ) 4 ( ) 3 ( ) () 32 hf39。 x f x f x f fxh ?? ? ?