迭代法隨想-文庫(kù)吧資料
2024-10-25 10:12本頁(yè)面
【正文】 0 %格式 : x=nanewton(fname,dfname,x0,e,N) fname 和 dfname分別表示 f(x)及其導(dǎo)函數(shù)的 M % 函數(shù)句柄或內(nèi)嵌函數(shù) ,x0為迭代初值 ,e為精度要求 (默認(rèn)值為10^{4}),x為返回解 , % 并顯示計(jì)算過(guò)程設(shè)置迭代次數(shù)上限 N以防發(fā)散 (默認(rèn) 500) if nargin5,N=500。( )nnnnfxxxfx? ??擬 Newton迭代法： ( 0 , 1 , 2 , )n ? …10()39。( )x? 則迭代法 1 ()nnxx?? ? ( 0 , 1 , 2 , )n ? …*| 39。1n n nx x x xL ?? ? ??*10| | | | .1nnLx x x xL? ? ?? 簡(jiǎn)單迭代法 定理 2（收斂充分條件） 若存在區(qū)間（ c,d） ,使 （ 1）方程 x= (x)在 (c,d)內(nèi)有實(shí)根 （ 2） 在 (c,d)內(nèi)連續(xù)且 *。( )x?（ 1） [a,b]上 存在 , 且 | 39。end x=(a+b)/2 end 二分法 Matlab程序 先將方程 簡(jiǎn)單迭代法 ()xx??轉(zhuǎn)化為等價(jià)方程 ( ) 0fx ?然后從某個(gè)數(shù) 出發(fā) ,通過(guò)計(jì)算 0x1 ()nnxx?? ? ( 0 , 1 , 2 , )n ? …構(gòu)造序列 。else a=x。 if fa*fx0,b=x。)。 if fa*fb0,error(39。 fa=feval(fname,a)。 function x=nabisect(fname,a,b,e) %用途 : 二分法解非線性方程 f(x)=0
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1