D方程近似解PPT學(xué)習(xí)教案

1、會計學(xué)1D方程近似解方程近似解，內(nèi)只有一個根在若方程,0)(baxf內(nèi)嚴(yán)格單調(diào)）（在且baxf,)(為則稱,ba.其隔根區(qū)間0)()(, ,)(bfafbaCxf為隔根區(qū)間,ba(1) 作圖法 1. 求隔根區(qū)間的一般方法求隔根區(qū)間的一般方法 ;)(估計隔根區(qū)間的草圖由xfy 轉(zhuǎn)化為等價方程將0)(xfOxy)(xfy .)(, )(的草圖估計隔根區(qū)間由xyxyab)()(xx)(xy)(xyOxyab第1頁/共16頁01,3 xx方程例如13 xx由圖可見只有一個實(shí)根, )5 . 1, 1 (可轉(zhuǎn)化為.)5 . 1, 1 (即為其隔根區(qū)間,的左端點(diǎn)出發(fā)從區(qū)間ba以定步長 h 一步步向右搜索,

2、若0) 1()(hjafjhaf) 1(;, 1,0(bhjaj.) 1(內(nèi)必有根，則區(qū)間hjajha搜索過程也可從 b 開始 , 取步長 h 0 .xy213xy 1 xyO第2頁/共16頁只有且方程0)(xf1a1b，設(shè),)(baCxf，0)()(bfaf，一個根),(ba取中點(diǎn)，21ba1，若0)(1f.1即為所求根則，若0)()(1faf, ),(1a則根;,111baa令, ),(1b否則對新的隔根區(qū)間,11ba重復(fù)以上步驟,反復(fù)進(jìn)行,得 ,111bba令,11nnbababa的中點(diǎn)若取,nnba則誤差滿足)(211nnnab )(121abnab)(211nnnba ,的近似根作為

3、0 n1a1b第3頁/共16頁04 . 19 . 01 . 123xxx的近似實(shí)根時,要使誤差不超過,103至少應(yīng)對分區(qū)間多少次 ?解解: 設(shè) ,4 . 19 . 01 . 1)(23xxxxf),()(Cxf則9 . 02 . 23)(2xxxf)067. 5(0,),()(單調(diào)遞增在xf又,04 . 1)0(f06 . 1) 1 (f故該方程只有一個實(shí)根 , 1,0為其一個隔根區(qū)間欲使)01 (1211nn310必需,100021n即11000log2n96. 8可見只要對分區(qū)間9次 ,即可得滿足要求的實(shí)根近似值10第4頁/共16頁:)(滿足xf0)()(,) 1bfafba上連續(xù)在不變號

4、及上在)()(,)2xfxfba .),(0)(內(nèi)有唯一的實(shí)根在方程baxf有如下四種情況:xbayOxbayOxbayO00 ff00 ff00 ff00 ffxbayO第5頁/共16頁)()(0001xfxfxx程的近似根 .記縱坐標(biāo)與)(xf 同號的端點(diǎn)為,)(,(00 xfx用切線近似代替曲線弧求方1x在此點(diǎn)作切線 ,其方程為)()(000 xxxfxfy令 y = 0 得它與 x 軸的交點(diǎn), )0,(1x其中再在點(diǎn))(,(11xfx作切線 ,可得近似根.2x如此繼續(xù)下去, 可得求近似根的迭代公式 :)()(111nnnnxfxfxx),2, 1(n2x稱為牛頓迭代公式牛頓迭代公式 y

5、xabO0 x第6頁/共16頁1x2xyxabO0 x)()(111nnnnxfxfxx由微分中值定理得)()()(nnxffxf)(之間與在nx,0)(f)()(fxfxnn,0則得mxfxnn)(說明說明: 用牛頓法時,若過縱坐標(biāo)與)(xf 異號的端點(diǎn)作切線 ,則切線與 x 軸焦點(diǎn)的橫坐標(biāo)未必在.,內(nèi)ba)(min,xfmba記第7頁/共16頁(1) 簡化牛頓法簡化牛頓法若用一常數(shù)代替, )(1nxf即用平行, )()(10nxfxf代替例如用則得簡化牛頓迭代公式. 線代替切線,得)()(011xfxfxxnnn),2, 1(n優(yōu)點(diǎn):，避免每次計算)(1nxf因而節(jié)省計算量.缺點(diǎn): 逼近根

6、的速度慢一些. yxaO0 x第8頁/共16頁yx0 x1x為避免求導(dǎo)運(yùn)算 , )(1nxf用割線代替切線,2121)()(nnnnxxxfxf即用差商代替從而得迭代公式:)()()()(212111nnnnnnnxxxfxfxfxx2x3x(雙點(diǎn)割線法),3,2(n特點(diǎn)特點(diǎn): 逼近根的速度快于簡化牛頓法, 但慢于牛頓法.說明說明: 若將上式中,02xxn換為則為單點(diǎn)割線法,逼近根的速度與簡化牛頓法相當(dāng).O第9頁/共16頁074223xxx的近似解, 使誤差不超過 0.01 .解解:.742)(23xxxxf設(shè)由草圖可見方程有唯一的正實(shí)根 ,且9)4(,10)3(ff.43為一隔根區(qū)間，因此上

7、，由于在43443)(2xxxf)2)(23(xx046)( xxf)23(2x0)(min4, 3xfm11)3( fyx3 4O第10頁/共16頁,40 x故取得)4()4(41ffx289468. 3而mxfx)(111103. 109. 0，精度不夠故1x再求)68. 3()68. 3(68. 32ffx9 .2103. 168. 363. 3mxfx)(2211042. 001. 0004. 0因此得滿足精度要求的近似解63. 3yx3 4O第11頁/共16頁(補(bǔ)充補(bǔ)充) , )(0)(xxxf轉(zhuǎn)化為等價方程將方程在隔根區(qū),0 x間內(nèi)任取一點(diǎn)按遞推公式),2, 1()(1nxxnn,

8、nx生成數(shù)列,limnnx若則 即為原方程的根 .式稱為迭代格式 ,)(稱為迭代函數(shù)x稱為迭代0 x,lim存在稱迭代收斂若nnx初值 .否則稱為發(fā)散 .第12頁/共16頁.2, 1 013內(nèi)的實(shí)根在 xx解法解法1 將方程變形為, 13 xx迭代格式為, 131nnxx5 . 10 x取123nnx05 . 1375. 2396.12779.1903發(fā)散 !解法解法2 將方程變形為,13xx迭代格式為, 131nnxx5 . 10 x取12nnx05 . 135721. 133086. 17832472. 132472. 1迭代收斂 ,1.32472 為計算精度范圍內(nèi)的所求根 .第13頁/共16頁:,)(上滿足在區(qū)間方程baxxbxaxx)()(1且，連續(xù)）1)()(2Lxx且，存在），上有唯一解在方程）,)(1baxx nnnxxbax)(,210）(證明略)迭代法的斂散性與迭代函數(shù)的特性