數(shù)值分析非線性方程求根實驗

《數(shù)值分析非線性方程求根實驗》由會員上傳分享，免費在線閱讀，更多相關內容在學術論文-天天文庫。

1、實驗報告一、實驗目的1.迭代閑數(shù)對收斂性的影響。2.初伉的選擇對收斂性的影響。二、實驗題目1.用簡單迭代法求方程/Cr)=x3-x—1=0的根。最大迭代次數(shù)為500,觀察其計算結果并加以分析。分別化方程為如下等價方程：%=+1；%=X3-1;X=取初值x()=1.5,精度為KT4，2.①用牛頓法求方程/(x)=x3+X-1=0在0.5附近的根，分別取初值乂=-1000-100-2-1-0.5,0.54,2,100,1000觀察并比較計算結果，并加以分析。②用牛頓法求方程/(%)=%3-x=0所有根。三、實驗原

2、理簡單迭代法程序，牛頓迭代法程序。實驗內容及結果1.用簡單選代法求方程=X—1=0的根：1)先輸入簡單迭代法的M程序，在MATLAB的M文件窗口中輸入以下內容并保存:functionx=maiter(phi,x0,ep,N)ifnargin<4N=500;endifnargin<3ep=le-4;endk=0;whilek

3、，num2str(k)])1)再輸入以下程序，在M文件窗口輸入以下內容并保存:clc;clear;formatlongxl=maiter(inline('(x+1)A(1/3)'),1.5,le-5)x2=maiter(inline(*xA3-l*),1.5,le-5)x3=maiter(inline(’(1+1/x)A(1/2)’)，1.5,le-5)x4=maiter(inline('(xA3+x-l)/2')z1.5,le-5)2)主窗口品示的結果為：k=6x=1.324719474534364k=50

4、0x2=Infk=7x=1.32471868894279lk=500x4=Inf2.用牛頓迭代法求方程=x3+%—1=0在0.5附近的根：1)先輸入牛頓迭代法的M程序，在MATLAB的M文件窗口屮輸入以下內容并保存:functionx=manewton(fun,dfun,xO,ep,N)ifnargin<5,N=500;endifnargin<4,ep=le-4;endk=0;whilek

5、=x;k=k+l;endifk==N,warning(，已達迭代次數(shù)上限’)；enddisp(['k='znum2str(k)])2)再輸入以下程序，在M文件窗口輸入以下內容并保存：clc;clear;formatlongfun=inline(■xA3-x-l1);dfun=inline('3*x^2-l');-1000,xl=manewton(fun,dfun,-1000,le-4)-100,x2=manewton(fun,dfun,-100,le-4)一2,x3=manewtor?(fun,dfun,-2

6、,le-4)-1,x4=manewton(fun,dfun,-1,le-4)-0.5,x5=manewton(fun,dfun,-0?5,le-4)0.5,x6=manewton(fun,dfun,0.5,le-4)1,x7=manewton(fun,dfun,1,le-4)1,x8=manewton(fun,dfun,2,le-4)100,x9=manewton(fun,dfun,100,le-4)1000,xlO=manewton(fun,dfun,1000,le-4)1)在MATLAB的主程序窗U輸出以

7、下結果：ans=-1000k=21xl=0.682327804075895ans=-10()k=16x2=0.682327803903413ans=-2k=6x3=0.682327803828020ans=-1k=5x4=0.682327803828020ans=-0.500000000000000k=4x5=0.682327803903932ans=().500000000000000k=3x6=0.682327803828347ans=1k=3x7=0.682327803946513ans=2k=5x8=

8、0.682327803828368ans=100k=14x9=0.682327807263074ans=1000k=20x!0=0.6823278038536191.用牛頓法求方程=x=0所有根：1)由于之前已經(jīng)輸過牛頓迭代法的M程序，所以可以直接調用。輸入以下內容并保存:clc;clear;formatlongfun=inline(1xA3-x1);dfun=inline(*3*xz<2-l*);