東北大學(xué) 數(shù)值分析常微分方程數(shù)值解法

《東北大學(xué) 數(shù)值分析常微分方程數(shù)值解法》由會(huì)員分享，可在線(xiàn)閱讀，更多相關(guān)《東北大學(xué) 數(shù)值分析常微分方程數(shù)值解法（48頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、會(huì)計(jì)學(xué)1東北大學(xué)東北大學(xué) 數(shù)值分析數(shù)值分析 常微分方程數(shù)值解常微分方程數(shù)值解法法 所謂數(shù)值解法,就是設(shè)法將常微分方程離散化,建立差分方程,給出解在一些離散點(diǎn)上的近似值. a=x0 x1x2xnxN=b其中剖分節(jié)點(diǎn)xn=a+nh,n=0,1,N, h稱(chēng)為剖分步長(zhǎng).數(shù)值解法就是求精確解y(x)在剖分節(jié)點(diǎn)xn上的近似值yny(xn), n=1,2,N. 假設(shè)初值問(wèn)題(8.1)的解y=y(x)唯一存在且足夠光滑.對(duì)求解區(qū)域a,b做剖分 我們采用數(shù)值積分方法來(lái)建立差分公式. 1.2 構(gòu)造數(shù)值解法的基本思想構(gòu)造數(shù)值解法的基本思想 在區(qū)間xn,xn+1上對(duì)方程(8.1)做積分,則有第1頁(yè)/共48頁(yè)對(duì)右邊的積

2、分應(yīng)用左矩形公式，則有)2 . 8()(,()()(11nnxxnndxxyxfxyxy第2頁(yè)/共48頁(yè)梯形公式oxyab左矩形公式y(tǒng)=(x)babfafabdxxf)()(2)(baafabdxxf)()()(右矩形公式babfabdxxf)()()(中矩形公式babafabdxxf)2()()(第3頁(yè)/共48頁(yè)對(duì)右邊的積分應(yīng)用左矩形公式，則有)2 . 8()(,()()(11nnxxnndxxyxfxyxy因此，建立節(jié)點(diǎn)處近似值yn滿(mǎn)足的差分公式稱(chēng)之為Euler公式公式. 稱(chēng)為梯形公式梯形公式. )(,()()(1nnnnxyxhfxyxy),(1nnnnyxhfyy1,2 , 1 , 0

3、,0Nny 若對(duì)(8.2)式右邊的積分應(yīng)用梯形求積公式，則可導(dǎo)出差分公式1,2 , 1 , 0,0Nny),(),(2111nnnnnnyxfyxfhyy第4頁(yè)/共48頁(yè) 利用Euler方法求初值問(wèn)題 解解 此時(shí)的Euler公式為稱(chēng)為Euler中點(diǎn)公式中點(diǎn)公式或稱(chēng)雙步雙步Euler公式公式. 若在區(qū)間xn-1,xn+1上對(duì)方程(8.1)做積分,則有11)(,()()(11nnxxnndxxyxfxyxy對(duì)右邊的積分應(yīng)用中矩形求積公式，則得差分公式),(211nnnnyxhfyy1,2 , 1 , 0,0Nny20 ,21122xyxy 0)0(y的數(shù)值解.此問(wèn)題的精確解是y(x)=x/(1+x

4、2).第5頁(yè)/共48頁(yè)分別取步長(zhǎng)h=0.2 ,0.1 ,0.05,計(jì)算結(jié)果如下)211(221nnnnyxhyy2 , 1 , 0,00ny第6頁(yè)/共48頁(yè)hxnyny(xn)y(xn)-ynh=0.20.000.400.801.201.602.000.000000.376310.542280.527090.466320.406820.000000.344830.487800.491800.449440.400000.00000-0.03148-0.05448-0.03529-0.01689-0.00682h=0.10.000.400.801.201.602.000.000000.360850

5、.513710.509610.458720.404190.000000.344830.487800.491800.449440.400000.00000-0.01603-0.02590-0.01781-0.00928-0.00419h=0.050.000.400.801.201.602.000.000000.352870.500490.500730.454250.402270.000000.344830.487800.491800.449440.400000.00000-0.00804-0.01268-0.00892-0.00481-0.00227第7頁(yè)/共48頁(yè)Euler中點(diǎn)公式則不然, 計(jì)

6、算yn+1時(shí)需用到前兩步的值yn , yn-1 ,稱(chēng)其為兩步方法兩步方法,兩步以上的方法統(tǒng)稱(chēng)為多步法多步法. 在Euler公式和梯形公式中,為求得yn+1,只需用到前一步的值yn,這種差分方法稱(chēng)為單步法單步法,這是一種自開(kāi)始方法. 隱式公式中,每次計(jì)算yn+1都需解方程,要比顯式公式需要更多的計(jì)算量,但其計(jì)算穩(wěn)定性較好. 在Euler公式和Euler中點(diǎn)公式中,需要計(jì)算的yn+1已被顯式表示出來(lái),稱(chēng)這類(lèi)差分公式為顯式公式顯式公式,而梯形公式中,需要計(jì)算的yn+1隱含在等式兩側(cè),稱(chēng)其為隱式公式隱式公式.第8頁(yè)/共48頁(yè) 從數(shù)值積分的角度來(lái)看,梯形公式計(jì)算數(shù)值解的精度要比Euler公式好,但它屬于

7、隱式公式,不便于計(jì)算. 實(shí)際上,常將Euler公式與梯形公式結(jié)合使用: 2.1 改進(jìn)的改進(jìn)的Euler方法方法),(),(2111nnnnnnyxfyxfhyy1,2 , 1 , 0,0Nny),(01nnnnyxhfyy),(),(21111knnnnnknyxfyxfhyy1,2 , 1 , 0,0Nny第9頁(yè)/共48頁(yè) 由迭代法收斂的角度看，當(dāng) (是給定的精度要求)時(shí), 取 就可以保證迭代公式收斂, 而當(dāng)h很小時(shí), 收斂是很快的. 而且, 只要|111knknyy.111knnyy, 12Lyfh),(1nnnnyxhfyy),(),(2111nnnnnnyxfyxfhyy1,2 , 1

8、 , 0,0Nny 通常采用只迭代一次的算法:稱(chēng)之為改進(jìn)的改進(jìn)的Euler方法方法. 這是一種單步顯式方法. 改進(jìn)的Euler方法也可以寫(xiě)成第10頁(yè)/共48頁(yè))(2211KKhyynn),(1nnyxfK 1,2 , 1 , 0,0Nny y=y-2x/y , 0 x1的數(shù)值解, 取步長(zhǎng)h=0.1 . 精確解為y(x)=(1+2x)1/2.),(12hKyhxfKnn y(0)=1 解解 (1) 利用Euler方法nnnnyxyy/2 . 01 . 119 ,2 , 1 , 0,10ny第11頁(yè)/共48頁(yè))(05. 0211KKyynnnnnyxyK/219 ,2 , 1 , 0,10ny計(jì)算

9、結(jié)果如下:1121 . 0) 1 . 0(21 . 0KyxKyKnnn (2) 利用改進(jìn)Euler方法第12頁(yè)/共48頁(yè)nxnEuler方法yn改進(jìn)Euler法yn精確解y(xn)01234567891000.10.20.30.40.50.60.70.80.9111.11.1918181.2774381.3582131.4351331.5089661.5803381.6497831.7177791.78477011.0959091.1840961.2662011.3433601.4164021.4859561.5525151.6164761.6781681.73786911.0954451.

10、1832161.2649911.3416411.4142141.4832401.5491931.6124521.6733201.732051第13頁(yè)/共48頁(yè) 在節(jié)點(diǎn)xn+1的誤差y(xn+1)-yn+1 ,不僅與yn+1這一步計(jì)算有關(guān),而且與前n步計(jì)算值yn,yn-1,y1都有關(guān). 為了簡(jiǎn)化誤差的分析,著重研究進(jìn)行一步計(jì)算時(shí)產(chǎn)生的誤差.即假設(shè)yn=y(xn),求誤差y(xn+1)-yn+1,這時(shí)的誤差稱(chēng)為局部截?cái)嗾`差局部截?cái)嗾`差,它可以反映出差分公式的精度. 如果單步差分公式的局部截?cái)嗾`差為O(hp+1),則稱(chēng)該公式為p階方法階方法.這里p為非負(fù)整數(shù).顯然,階數(shù)越高,方法的精度越高. 研究差

11、分公式階的重要手段是Taylor展開(kāi)式,一元函數(shù)和二元函數(shù)的Taylor展開(kāi)式為:第14頁(yè)/共48頁(yè)另外,在yn=y(xn)的條件下,考慮到y(tǒng)(x)=(x,y(x),則有 321! 3)(! 2)()()()()(hxyhxyhxyxyhxyxynnnnnn2222222),(),(2),(! 21),(),(),(),(kyyxfhkyxyxfhxyxfkyyxfhxyxfyxfkyhxfnnnnnnnnnnnnnn y(xn)=(xn,y(xn)=(xn,yn)=n y(xn)=nnnnnfyfxfxyxfdxd)(,(nnnnnnnnnfyfyfxffyffyxfxf2222222)(2

12、)(nnnnfyfxfdxdxy 第15頁(yè)/共48頁(yè) yn+1=yn+h(xn,yn) 對(duì)Euler方法，有 21! 2)()()()()(hxyhxyxyhxyxynnnnn =yn+(xn,yn)h+O(h2)從而有: y(xn+1)-yn+1=O(h2)所以Euler方法是一階方法.再看改進(jìn)Euler方法, 因?yàn)?,(12hKyhxfKnn1hKyfhxffnnn)(221321222122222hOKhyfKhyxfhxfnnn可得第16頁(yè)/共48頁(yè)所以, 改進(jìn)的Euler方法是二階方法.而nnnnnnfyxfhhfyy221)(244222223hOfyffyxfxfhnnnnn)(

13、! 3)(2)()()()(4321hOhxyhxyhxyxyxynnnnn nnnnnfyfxfhhfy22)()(26422222223hOfyfyfxffyffyxfxfhnnnnnnnnn從而有: y(xn+1)-yn+1=O(h3) 設(shè)y(x)是初值問(wèn)題(8.1)的精確解, 利用Taylor展開(kāi)式可得第17頁(yè)/共48頁(yè)稱(chēng)之為p階Taylor展開(kāi)方法. 1)1()(21)!1()(!)(! 2)()()()( pppnpnnnnhpyhPxyhxyhxyxyxy)()(,(!)(,(! 2)(,()(1)1()1(2pnnppnnnnnhOxyxfPhxyxfhxyxhfxy因此，可建

14、立節(jié)點(diǎn)處近似值yn滿(mǎn)足的差分公式),(!),(! 2),()1()1(21nnppnnnnnnyxfPhyxfhyxhfyy1,2 , 1 , 0,0Nny),(),(),(),()1(yxfyyxfxyxfyxffyfyfxffyffyxfxfyxf2222222)2()(2),(其中第18頁(yè)/共48頁(yè)所以,此差分公式是p階方法. 由于Taylor展開(kāi)方法涉及很多復(fù)合函數(shù)(x,y(x)的導(dǎo)數(shù)的計(jì)算,比較繁瑣,因而很少直接使用,經(jīng)常用它為多步方法提供初始值.然而, Taylor展開(kāi)方法給出了一種構(gòu)造單步顯式高階方法的途徑. Euler方法可寫(xiě)為 可見(jiàn),公式的局部截?cái)嗾`差為: y(xn+1)-y

15、n+1=O(hp+1). 3.1 Runge-Kutta方法的構(gòu)造方法的構(gòu)造hKyynn1),(nnyxfK第19頁(yè)/共48頁(yè) 構(gòu)造差分公式 改進(jìn)的Euler方法可寫(xiě)為)(2211KKhyynn),(1nnyxfK ),(12hKyhxfKnn)(22111ppnnKKKhyy),(1nnyxfK ),(12122KhyhxfKnn),(11piipinPnPKhyhxfK其中i,i,ij為待定參數(shù). 若此公式的局部截?cái)嗾`差為第20頁(yè)/共48頁(yè)由于 yn+1=yn+h1n+h2(n+hxn+hn yn)+O(h3)O(hp+1),稱(chēng)公式為p階階Runge-kutta方法方法,簡(jiǎn)稱(chēng)p階階R-K方

16、法方法. 對(duì)于p=2的情形, 應(yīng)有)(22111KKhyynn),(1nnyxfK ),() 3 . 8(12hKyhxfKnn =yn+h(1+2)n+h22(xn+n yn)+O(h3)(2)(321hOfffhhfyxyynnxnnnn所以,只要令 1+2=1, 2=1/2, 2=1/2 (8.4)第21頁(yè)/共48頁(yè) 一般地一般地, 參數(shù)由參數(shù)由(8.4)確定的一族差分公式確定的一族差分公式(8.3)統(tǒng)稱(chēng)為二階統(tǒng)稱(chēng)為二階R-K方法方法.稱(chēng)之為中點(diǎn)公式中點(diǎn)公式,或可寫(xiě)為若取=1,則得1=2=1/2,=1,此時(shí)公式(8.3)就是改進(jìn)的Euler公式; 若取1=0,則得2=1,=1/2,公式(

17、8.3)為21hKyynn),(1nnyxfK ),(121212hKyhxfKnn),(,(21211nnnnnnyxhfyhxhfyy 高階R-K公式可類(lèi)似推導(dǎo). 下面列出常用的三階、四階R-K公式.第22頁(yè)/共48頁(yè) 四階標(biāo)準(zhǔn)四階標(biāo)準(zhǔn)R-K公式公式 三階三階R-K公式公式)4(63211KKKhyynn),(1nnyxfK )2,(213hKhKyhxfKnn)22(643211KKKKhyynn),(1nnyxfK ),(121212hKyhxfKnn),(34hkyhxfKnn),(121212hKyhxfKnn),(221213hKyhxfKnn第23頁(yè)/共48頁(yè) 解解 四階標(biāo)準(zhǔn)R

18、-K公式為)22(4321611KKKKhyynnnnnyxyK/21 y=y-2x/y , 0 x1 y(0)=1的數(shù)值解, 取步長(zhǎng)h=0.2 .)/()2(2212213hKyhxhKyKnnn)/()2(1211212hKyhxhKyKnnn)/()(2334hKyhxhKyKnnn計(jì)算結(jié)果如下:第24頁(yè)/共48頁(yè)nxnyny(xn)nxnyny(xn)0120.00.20.41.001.18321.34171.001.18321.34163450.60.81.01.48331.61251.73211.48321.61251.7321 也可以構(gòu)造隱式R-K方法,其一般形式為prrrnnK

19、hyy11prKhyhxfKpiirinrnr, 2 , 1,),(1稱(chēng)之為p級(jí)隱式級(jí)隱式R-K方法方法,同顯式R-K方法一樣確定參數(shù).如)(21211KKhyynn),(1nnyxfK ),(2211212hKhKyhxfKnn第25頁(yè)/共48頁(yè)是二級(jí)二階隱式R-K方法,也就是梯形公式.但是p級(jí)隱式R-K方法的階可以大于p,例如,一級(jí)隱式中點(diǎn)公式為11hKyynn),(121211hKyhxfKnn或?qū)憺?(,(121211nnnnnyyhxhfyy它是二階方法. 以p階R-K方法為例討論.設(shè)從xn以步長(zhǎng)h計(jì)算y(xn+1)的近似值為)(1hny ,局部截?cái)嗾`差為1)(11)(phnnChy

20、xy其中,C是與h無(wú)關(guān)的常數(shù).第26頁(yè)/共48頁(yè) 如果將步長(zhǎng)減半,取h/2為步長(zhǎng), 從xn經(jīng)兩步計(jì)算得到y(tǒng)(xn+1)的近似值記為 ,其局部截?cái)嗾`差為于是有從而,得到事后誤差估計(jì))(12hny11)(1121)2(2)(2pppnnChhCyxyh可見(jiàn),當(dāng)phnnnnyxyyxyh21)()()(11)(112)(121)()(1)(1)(1122hnnpnnyyyxyhh|)(1)(12hnnyyh 成立時(shí),可取)(112)(hnnyxy .否則,應(yīng)將步長(zhǎng)再次減半進(jìn)行計(jì)算.第27頁(yè)/共48頁(yè) 求解初值問(wèn)題的單步顯式方法可一統(tǒng)一寫(xiě)為如下形式 yn+1=yn+h(xn,yn,h) (8.5) 對(duì)

21、于Euler方法,有 4.1 單步方法的收斂性單步方法的收斂性 y=(x,y) ,axb y(a)= 其中(x,y,h)稱(chēng)為增量函數(shù)增量函數(shù). (x,y,h)=(x,y)對(duì)于改進(jìn)的Euler方法,有 (x,y,h)=1/2(x,y)+(x+h,y+h(x,y)第28頁(yè)/共48頁(yè) 設(shè)y(x)是初值問(wèn)題(8.1)的解 ,yn是單步法 (8.5)產(chǎn)生的近似解.如果對(duì)任意固定的點(diǎn)xn,均有y(xn),則稱(chēng)單步法(8.5)是收斂的. 可見(jiàn),若方法(8.5)是收斂的,則當(dāng)h0時(shí),整體截?cái)嗾`差en=y(xn)-yn將趨于零. 定理定理8.1 設(shè)單步方法(8.5)是p1階方法, 增量函數(shù)(x,y,h)在區(qū)域a

22、xb,-yn)的變化均不超過(guò) ,則稱(chēng)此差分方法是絕對(duì)穩(wěn)定絕對(duì)穩(wěn)定的. 討論數(shù)值方法的穩(wěn)定性,通常僅限于典型的試驗(yàn)方程 y=y 其中是復(fù)數(shù)且Re()0. 在復(fù)平面上,當(dāng)方法穩(wěn)定時(shí)要求變量h的取值范圍稱(chēng)為方法的絕對(duì)穩(wěn)定域絕對(duì)穩(wěn)定域,它與實(shí)軸的交集稱(chēng)為絕對(duì)穩(wěn)定區(qū)間絕對(duì)穩(wěn)定區(qū)間. 第32頁(yè)/共48頁(yè) 將Euler方法應(yīng)用于方程y=y, 得到 設(shè)在計(jì)算yn時(shí)產(chǎn)生誤差n,計(jì)算值yn=yn+n,則n將對(duì)以后各節(jié)點(diǎn)值計(jì)算產(chǎn)生影響.記ym=ym+m ,mn,由上式可知誤差m滿(mǎn)足方程 m=(1+h)m-1=(1+h)m-nn , mn 對(duì)隱式單步方法也可類(lèi)似討論.如將梯形公式用于方程y=y,則有 yn+1=yn+

23、h/2 (yn+yn+1) yn+1=(1+h)yn 可見(jiàn),若要|m|n|,必須且只須|1+h|1 ,因此Euler法的絕對(duì)穩(wěn)定域?yàn)閨1+h|1,絕對(duì)穩(wěn)定區(qū)間是-2Re()h0.解出yn+1得 第33頁(yè)/共48頁(yè)nnyhhy2121111類(lèi)似前面分析,可知絕對(duì)穩(wěn)定區(qū)域?yàn)?112121hh由于Re()0,所以此不等式對(duì)任意步長(zhǎng)h恒成立,這是隱式公式的優(yōu)點(diǎn). 一些常用方法的絕對(duì)穩(wěn)定區(qū)間為方 法方法的階數(shù)穩(wěn) 定 區(qū) 間Euler方法梯形方法改進(jìn)Euler方法二階R-K方法三階R-K方法四階R-K方法122234(-2 , 0)(- , 0)(-2 , 0)(-2 , 0)(-2.51 , 0)(-2

24、.78 , 0)第34頁(yè)/共48頁(yè) 解解 因y0=1,計(jì)算得y10=1024,而y(1)=9.35762310-14. y=-30y , 0 x1 y(0)=1取步長(zhǎng)h=0.1 ,利用Euler方法計(jì)算y10y(1). y(x)=e-30 x 這是因?yàn)閔=-3不屬于Euler方法的絕對(duì)穩(wěn)定區(qū)間. 若取h=0.01,計(jì)算得y100=3.23447710-16. 若取h=0.001,計(jì)算得y1000=5.91199810-14. 若取h=0.0001,計(jì)算得y10000=8.94505710-14. 若取h=0.00001,計(jì)算得y100000=9.315610-14.第35頁(yè)/共48頁(yè) 單步顯式

25、方法的穩(wěn)定性與步長(zhǎng)密切相關(guān), 在一種步長(zhǎng)下是穩(wěn)定的差分公式,取大一點(diǎn)步長(zhǎng)就可能是不穩(wěn)定的. 收斂性是反映差分公式本身的截?cái)嗾`差對(duì)數(shù)值解的影響;穩(wěn)定性是反映計(jì)算過(guò)程中舍入誤差對(duì)數(shù)值解的影響.只有即收斂又穩(wěn)定的差分公式才有實(shí)用價(jià)值. 由于在計(jì)算yn+1時(shí) ,已經(jīng)知道yn ,yn-1 ,及(xn,yn), (xn-1,yn-1),利用這些值構(gòu)造出精度高、計(jì)算量小的差分公式就是線(xiàn)性多步法.5.1 利用待定參數(shù)法構(gòu)造線(xiàn)性多步方法利用待定參數(shù)法構(gòu)造線(xiàn)性多步方法 r+1步線(xiàn)性多步方法的一般形式為第36頁(yè)/共48頁(yè)當(dāng)-10時(shí),公式為隱式公式,反之為顯式公式.參數(shù)i,i的選擇原則是使方法的局部截?cái)嗾`差為 y(x

26、n+1)-yn+1=O(h)r+2 選取參數(shù),0,1,2,使三步方法 yn+1=yn+h(0n+1n-1+2n-2) 這里,局部截?cái)嗾`差是指 ,在yn-i=y(xn-i),i=0,1,r的前提下,誤差y(xn+1)-yn+1.為三階方法. ririiniininfhyy011 解解 設(shè)yn=y(xn),yn-1=y(xn-1),yn-2=y(xn-2),則有 第37頁(yè)/共48頁(yè) n=(xn,y(xn)=y(xn) y(xn+1)=y(xn)+hy(xn)+1/2h2y(xn)+1/6h3y(xn) 于是有若使: y(xn+1)-yn+1=O(h4) ,只要,0,1,2滿(mǎn)足: n-1=(xn-1

27、,y(xn-1)=y(xn-1)=y(xn-h) =y(xn)-hy(xn)+1/2h2y(xn)-1/6h3y(4)(xn)+O(h4) n-2=y(xn)-2hy(xn)+2h2y(xn)-4/3h3y(4)(xn)+O(h4) yn+1=y(xn)+h(0+1+2)y(xn)-h2(1+22)y(xn) +h3(1/21+22)y(xn)-h4/6(1+82)y(4)(xn)+O(h5) +1/24h4y(4)(xn)+O(h5) 第38頁(yè)/共48頁(yè) =1, 0+1+2=1, 1+22=-1/2 , 1+42=1/3于是有三步三階顯式差分公式設(shè)pr(x)是函數(shù)(x,y(x)的某個(gè)r次插值

28、多項(xiàng)式,則有解之得: yn+1=yn+h/12(23n-16n-1+5n-2) 因?yàn)?125,34,1223, 12101)(,()()(1nnxxnndxxyxfxyxy1)()()(1nnxxnrnnRdxxpxyxy其中 第39頁(yè)/共48頁(yè) 選取不同的插值多項(xiàng)式pr(x),就可導(dǎo)出不同的差分公式.下面介紹常用的Adams公式公式. 設(shè)已求得精確解y(x)在步長(zhǎng)為h的等距節(jié)點(diǎn)xn-r,xn上的近似值yn-r ,yn , 記k=(xk,yk) ,利用r+1個(gè)數(shù)據(jù)(xn-r,n-r),(xn,n)構(gòu)造r次Lagrange插值多項(xiàng)式由此,可建立差分公式 1.Adams顯式公式顯式公式 1)(1n

29、nxxrnndxxpyy1)()(,(nnxxrndxxpxyxfR其中 rjjnjnrfxlxp0)()(第40頁(yè)/共48頁(yè)由此,可建立差分公式 由于 rjxxxxxlrjkkknjnknjn, 1 , 0)()()(0 hrj jnxxjnrjnnfdxxlyynn1)(01)()()()(110thxxdxxxxxdxxlnxxrjkkknjnknxxjnnnnn令100, 1 , 0,)()!( !) 1(rjkkjrjdtktjrjh則有 rjjnrjnnfhyy01稱(chēng)之為r+1步步Adams顯式公式顯式公式. 第41頁(yè)/共48頁(yè)下面列出幾個(gè)帶有局部截?cái)嗾`差主項(xiàng)的Adams顯式公式

30、r=0 yn+1=yn+hn+(1/2)h2y(xn) 2.Adams隱式公式隱式公式 r=1 yn+1=yn+(h/2)(3n-n-1)+(5/12)h3y(xn) r=2 yn+1=yn+(h/12)(23n-16n-1+5n-2) +(3/8)h4y(4)(xn) r=3 yn+1=yn+(h/24)(55n-59n-1+37n-2-9n-3) +(251/720)h5y(5)(xn) 如果利用r+1個(gè)數(shù)據(jù)(xn-r+1,n-r+1),(xn+1,n+1)構(gòu)造r次Lagrange插值多項(xiàng)式pr(x),則可導(dǎo)出數(shù)值穩(wěn)定性好的隱式公式,稱(chēng)為Adams隱式公式隱式公式,其一般形式為第42頁(yè)/共

31、48頁(yè)其中系數(shù)為 010*, 1 , 0,)()!( !) 1(rjkkjrjrjdtktjrjrjjnrjnnfhyy01*1下面列出幾個(gè)帶有局部截?cái)嗾`差主項(xiàng)的Adams隱式公式 r=0 yn+1=yn+hn+1-(1/2)h2y(xn) r=1 yn+1=yn+(h/2)(n+n+1)-(1/12)h3y(xn) r=2 yn+1=yn+(h/12)(5n+1+8n-n-1) -(1/24)h4y(4)(xn) r=3 yn+1=yn+(h/24)(9n+1+19n-5n-1+n-2) -(19/720)h5y(5)(xn)第43頁(yè)/共48頁(yè) 3.Adams預(yù)估預(yù)估-校正公式校正公式 由顯

32、式公式提供一個(gè)預(yù)估值，再用隱式公式校正一次，求得數(shù)值解，稱(chēng)為預(yù)估校正方法預(yù)估校正方法。 校正 yn+1=yn+(h/24)(9n+1+19n-5n-1+n-2) 一般預(yù)估公式和校正公式都采用同階公式。例如： 預(yù)估 yn+1=yn+(h/24)(55n-59n-1+37n-2-9n-3) n+1=(xn+1,yn+1) , n=3,4,稱(chēng)為四階Adams預(yù)估校正公式.實(shí)際計(jì)算時(shí)通常用四階單步方法(如四階R-K公式)為它提供起始值y1,y2,y3 . 用四階Adams預(yù)估校正公式求解初值問(wèn)題 第44頁(yè)/共48頁(yè) y=y-2x/y , 0 x1 y(0)=1取步長(zhǎng)h=0.1. 解 用四階R-K公式提

33、供起始值,計(jì)算結(jié)果如下xnR-k法yn預(yù)估值yn校正值yn精確值y(xn)00.10.20.30.40.50.60.70.80.9111.0954461.1832171.2649121.3415511.4140451.4830171.5489171.6121141.6729141.7315661.3416411.4142131.4832391.5491921.6124501.6733181.73204811.0954451.1832161.2649911.3416411.4142141.4832401.5491931.6124521.6733201.732051第45頁(yè)/共48頁(yè)第第250頁(yè)頁(yè) 習(xí)題習(xí)題88-5, 8-7, 8-8, 8-11, 8-12, 8-13, 8-15 第46頁(yè)/共48頁(yè)課間休息課間休息第47頁(yè)/共48頁(yè)