




下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、迭代函數(shù)對收斂性的影響一、實(shí)驗(yàn)?zāi)康模撼醪搅私夥蔷€性方程的簡單迭代法及其收斂性,體會迭代函數(shù)對收斂性的影響,知道當(dāng)?shù)瘮?shù)滿足什么條件時,迭代法收斂。二、實(shí)驗(yàn)內(nèi)容:用簡單迭代法求方程 的根。方案一: 化為等價方程 方案二: 化為等價方程 三、實(shí)驗(yàn)要求:(1)分別對方案一、方案二取初值,迭代10次,觀察其計算值,并加以分析。(2)用MATLAB內(nèi)部函數(shù)solve直接求出方程的所有根,并與(1)的結(jié)果進(jìn)行比較。四、迭代法程序functionk,piancha,xdpiancha,xk=diedai(x0,k)x(1)=x0;for i=1:kx(i+1)=fun1(x(i);piancha=abs(
2、x(i+1)-x(i);xdpiancha=piancha/(abs(x(i+1)+eps);i=i+1;xk=x(i);(i-1) piancha xdpiancha xkendif (piancha1)&(xdpiancha0.5)&(k3) disp(此迭代序列發(fā)散,請重新輸入新的迭代公式)return;endif (piancha0.001)&(xdpiancha3)disp(此迭代序列收斂,且收斂速度較快)return;endp=(i-1) piancha xdpiancha xk五、實(shí)驗(yàn)結(jié)果:方案一: 化為等價方程 建立M文件fun1.m的文件function y1=fun1(x)
3、y1=(x+1)./2).(1/3)在MATLAB窗口輸入程序 k,piancha,xdpiancha,xk=diedai(0,10)運(yùn)行后輸出結(jié)果y1 =0.7937ans =1.0000 0.7937 1.0000 0.7937y1 =0.9644ans =2.0000 0.1707 0.1770 0.9644y1 = 0.9940ans =3.0000 0.0297 0.0298 0.9940y1 =0.9990ans =4.0000 0.0050 0.0050 0.9990y1 = 0.9998ans =5.0000 0.0008 0.0008 0.9998y1 =1.0000ans
4、=6.0000 0.0001 0.0001 1.0000y1 = 1.0000ans =7.0000 0.0000 0.0000 1.0000y1 = 1.0000ans =8.0000 0.0000 0.0000 1.0000y1 =1.0000ans =9.0000 0.0000 0.0000 1.0000y1 =1.0000ans =10.0000 0.0000 0.0000 1.0000此迭代序列收斂,且收斂速度較快k =10 piancha = 1.0685e-07 xdpiancha =1.0685e-07 xk =1.00002、方案二: 化為等價方程 (1)建立M文件fun,m
5、的文件function y1=fun1(x)y1=2.*(x.3)-1(2)在MATLAB窗口輸入程序 k,piancha,xdpiancha,xk=diedai(0,10)(3)運(yùn)行后輸出結(jié)果y1 =-1ans =1.0000 1.0000 1.0000 -1.0000y1 =-3ans =2.0000 2.0000 0.6667 -3.0000y1 =-55ans =3.0000 52.0000 0.9455 -55.0000y1 =-ans =1.0e+05 * 0.0000 3.3270 0.0000 -3.3275y1 =-7.3687e+16ans =1.0e+16 * 0.000
6、0 7.3687 0.0000 -7.3687y1 =-8.0019e+50ans =1.0e+50 * 0.0000 8.0019 0.0000 -8.0019y1 =-1.0247e+153ans =1.0e+153 * 0.0000 1.0247 0.0000 -1.0247y1 =-Infans =8 Inf NaN -Infy1 =-Infans =9 NaN NaN -Infy1 =-Infans =10 NaN NaN -Infp = 10 NaN NaN -Infk = 10 piancha =NaN xdpiancha = NaN xk =-Inf3、用MATLAB內(nèi)部函數(shù)s
7、olve直接求(1)輸入程序 x=solve(2.*(x.3)-x-1=0)運(yùn)行后輸出的結(jié)果x = 1. -.+.*i -.-.*i表1 2個迭代公式的計算結(jié)果迭代第10次的偏差迭代第10次的相對偏差方案一1.0685e-071.0685e-070.79370.96440.99400.99900.99981.00001.00001.00001.00001.0000方案二NaNNaN-1-3-55-7.3687e+16-8.0019e+50-1.0247e+153-Inf-Inf-Inf6、 實(shí)驗(yàn)結(jié)果分析:從表1可以看出,方案一收斂很快,偏差和偏差的相對誤差幾乎為零;方案二根本不收斂,它的偏差piancha已經(jīng)NaN且相對誤差xdwucha的知也已經(jīng)NaN。由此可見,迭代序列的斂散性與迭代公式有關(guān),也與相鄰兩次迭代的偏差和偏差的相對誤差有關(guān)
溫馨提示
- 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)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 【正版授權(quán)】 ISO/IEC 15045-4-2:2024 EN Information technology - Home Electronic System (HES) gateway - Part 4-2: Structure - Simple gateway
- 【正版授權(quán)】 ISO 10993-4:2017/Amd 1:2025 EN Biological evaluation of medical devices - Part 4: Selection of tests for interactions with blood - Amendment 1
- 2025年度兒童接送服務(wù)與社區(qū)共建合作協(xié)議
- 2025年度出租車大包合同范本與合同法解讀
- 2025年新型不銹鋼罐體設(shè)計、制造與安裝集成合同
- 財務(wù)報表披露規(guī)定計劃
- 主管如何激勵高表現(xiàn)者計劃
- 倉庫精益管理的實(shí)施效果計劃
- 班級外聯(lián)活動的經(jīng)驗(yàn)分享計劃
- 課程創(chuàng)新與教學(xué)實(shí)驗(yàn)計劃
- 留置胃管課件
- 《紡織服裝材料》課件-項(xiàng)目6 紡織材料的水分及檢測
- (正式版)SH∕T 3548-2024 石油化工涂料防腐蝕工程施工及驗(yàn)收規(guī)范
- AQ/T 2059-2016 磷石膏庫安全技術(shù)規(guī)程(正式版)
- 四川省宜賓市中學(xué)2025屆九上數(shù)學(xué)期末統(tǒng)考模擬試題含解析
- 貴州人民版五年級勞動下冊教案
- 2024年包頭市水務(wù)(集團(tuán))有限公司招聘筆試沖刺題(帶答案解析)
- 知識庫管理規(guī)范大全
- 華能分布式光伏項(xiàng)目EPC總承包工程投標(biāo)文件-技
- 蜜雪冰城營銷策略研究開題報告
- 《游戲界面設(shè)計》課件-01單元-游戲界面設(shè)計概論
評論
0/150
提交評論