山東大學(xué)數(shù)值計算試驗報告9

1、山東大學(xué)軟件學(xué)院數(shù)值計算課程實驗報告學(xué)號：姓名：班級：5班實驗題目：1、對教材例5.8的方程，用不動點迭代法，列出迭代過程，畫圖顯示迭代過程（類似圖5.5,不必畫箭頭）。迭代最大次數(shù)N=30,誤差限E=10-2、用庫函數(shù)求解（選做）實驗學(xué)時：2實驗日期：2017年12月5日實驗?zāi)康模悍蔷€性方程求解硬件環(huán)境:PCW軟件環(huán)境:MATLAB實驗步驟與內(nèi)容:等價不動點函數(shù)為g(x)=xA2-2代碼如下：x=linspace(0,4);y=x;y1=x.A2-2;plot(x,y)%ffl出y=xholdonplot(x,y1)%®出g(x)holdonx1=2.1;P=；q=;fori=1:

2、3履代三次a=x1.A2-2;P=Px1;q=qa;P=Pa;q=qa;x1=a;endplot（p,q,'-bo'）holdon運行結(jié)果如下：4kFigure1X文件輯輯國直看國）mx（I）工具田桌面©國口他）帝助（H）Qdd、入踮0居0口日1412108602468101214等價不動點函數(shù)為g(x)=(x+2)A(1/2)代碼如下：x=linspace(0,4);y=x;y1=(x+2).A(1/2);plot(x,y)%®出y=xholdonplot(x,y1)%®出g(x)holdonx1=1;p=；q=;fori=1:3履代三次a=(x

3、1+2).A(1/2);P=Pxl；q=qa；P=Pa；q=qa;x1=a;endplot(p,q,'-bo')holdon運行結(jié)果：，F(xiàn)igure1X文件由案第直看MKA(I)IRffi桌面©國口兇)帚助(H)3日UiR二，-J®0乂-息口明等價不動點函數(shù)為g(x)=1+2/x代碼如下：x=linspace(0.1,4);y=x;y1=1+2./x;Plot(x,y)holdon%®出y=xplot(x,y1)%®出g(x)holdonx1=1;P=；q=;fori=1:10%S代十次a=1+2/x1;P=Px1;q=qa;P=Pa;q

4、=qa;x1=a;endplot(p,q,'-bo')holdon運行結(jié)果：KFigure1X文件(日編信查看M瓶入0)工具(D點面窗口也)幫助但)1.4等價不動點函數(shù)為g(x)=(xA2+2)/(2x-1)代碼如下：x=linspace(0.9,4);y=x;y1=(x.A2+2)./(2*x-1);plot(x,y)%ffl出y=xholdonplot(x,y1)%®出g(x)holdonx1=1;P=；q=;fori=1:6%迭代六次a=(x1A2+2)/(2*x1-1);P=Px1;q=qa;P=Pa;q=qa;x1=a;endPlot(P,q,'-bo')holdon運行結(jié)果:AFigure1文件由強國堂盲N)插入5工具CD桌面電)幫助(H),0)運行結(jié)果:心口Rir.gcriteriadetails)ansfsolvecompletedbecausethevectoroffunctLonvaluesIsnearzeroasmeasuredbythedefaultvalueofthefunctisntdlEran七丘,andthe口reblrmappa:rsregularasmeasuredb