計算機在材料科學(xué)與工程中的應(yīng)用ch04 非線性方程的解法.ppt_第1頁
計算機在材料科學(xué)與工程中的應(yīng)用ch04 非線性方程的解法.ppt_第2頁
計算機在材料科學(xué)與工程中的應(yīng)用ch04 非線性方程的解法.ppt_第3頁
計算機在材料科學(xué)與工程中的應(yīng)用ch04 非線性方程的解法.ppt_第4頁
計算機在材料科學(xué)與工程中的應(yīng)用ch04 非線性方程的解法.ppt_第5頁
已閱讀5頁,還剩14頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、第四章非線性方程的解法,主要內(nèi)容,1、引言,2、二分法,3、迭代法,(1),問題的提出 (2),求解步驟,計算機處理非線性方程的基本方法,常用的計算方法,4、牛頓迭代法和弦截法,內(nèi)容3的兩個例子,引言,哪些情況下,非線性方程通常用計算機來求解?,不易用解析的方法求解,通??捎脠D解法獲得近似解,在精確的場合則需要用數(shù)值方法,一般情況下,上述方程可寫成f(x)=0的形式,求根步驟,1.確定根的初始近似值。(粗略判斷有根范圍) 2.根的精確化(幾種方法:二分法、迭代法、牛頓迭代法、弦截法等等),第一步驟粗略判斷有根范圍,(1) 判斷f(x)在a,b連續(xù),且f(a)*f(b)0,且f(x0)f(xk+

2、1)0 ,則判斷根x* 在xk,xk+1內(nèi)。,例: 方程 f(x)=x3-x-1=0 ,找出有根范圍,要求范圍大小為h=0.5,(b)從x=0出發(fā),h=0.5,判斷f(x)的符號,(c)則粗略判斷根 x* 在 1,1.5 內(nèi),(d)再利用其他方法進行精確化,二分法,基本思想,a,平分含根的區(qū)間,判斷根的位置在哪個區(qū)間 b,舍去無根的區(qū)間,再進行a的判斷 c,重復(fù)a,b過程,直到達到事先要求的精度為止,計算過程,例: 方程 f(x)=x3-2x2-4x-7=0在3,4內(nèi)的根,精確到10-4,float a,b,eps,x,f0,f; a=3;b=4;eps=0.0001; f0=f(a); wh

3、ile(1) x=(a+b)/2; f=f(x); if(f0*f0) b=x; else a=x; if(fabs(b-a)eps) printf(x=%f,a); exit(0); ,說明,1,二分法是比較常用的求解非線性方程的方法,2,優(yōu)點1:原理簡單,容易操作,3,優(yōu)點2:不存在收斂性的問題,4,缺點:速度較慢,5,特殊情況,f(x)的單調(diào)性很重要,迭代法,基本思想 1,把一般方程f(x)=0化成等價形式x=g(x),構(gòu)造迭代式xk+1=g(xk),k=0,1,2, 2,由近似初始值x0,得到一系列迭代值:x0,x1,x2,xk,。記為xk,xk迭代序列 g(x)迭代函數(shù),如xk收斂于

4、x*,則x*=g(x*),即f(x*)=0??芍獂*是方程f(x)=0的解。,這種方法稱為迭代法,又叫逐次逼近法。,例:求方程f(x)=x-10 x+2=0的一個根,解:因為f(0)0,f(x)函數(shù)連續(xù),則0,1中必然有一個根。,將方程寫成等價形式x=lg(x+2),構(gòu)造迭代式:xk+1=lg(xk+2),取初始迭代值x0=1,進行迭代,則所求根為x*=0.3758,說明:這種方法有時無效。 例:本題如果選迭代式為xk+1=10 xk-2, 則:x0=1,x1=8,x2=108-2,溢出。xk不收斂,計算不出結(jié)果。,迭代法的收斂條件,計算步驟和流程圖,例:求方程x=e-x在x=0.5附近的一個根,計算精度10-3。,解:,牛頓迭代法和弦截法,牛頓迭代法 基本思想,例:求方程x=e-x在x=0.5附近的一個根,計算精度10-3。 (參考前面迭代法的例子),迭代過程,只需3次迭代,比普通迭代法(10次)快!,有關(guān)說明,弦截法 基本思想,與牛頓法的區(qū)別: 牛頓法只用到前一步迭代的數(shù)值單點迭代法。 弦截法用到前兩步迭代的結(jié)果多點迭代法,本章小結(jié),1/4、引言,2/4、

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論