




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、數(shù)值分析實(shí)驗報告姓名學(xué)號 日期2014-3-22實(shí)驗項目簡單迭代法與加速方法指導(dǎo)教師一、上機(jī)實(shí)驗的問題和要求(需求分析):目的與要求:1、掌握求解非線性方程實(shí)根的簡單迭代法的編程運(yùn)算2、會分析迭代步數(shù),設(shè)計容許誤差實(shí)驗內(nèi)容:1、方程可以等價化成以下三種形式:(i) (ii) (iii) 針對三種等價形式給出三種不同的簡單迭代格式并使用每種格式計算方程在區(qū)間上的解,初值選為1.5,容許誤差選為1.0E-5,即;分析每種格式的收斂性;分析收斂格式的迭代步數(shù)與計算時間.2、結(jié)合上述問題中(ii)相應(yīng)的迭代格式,利用Stenffenson迭代法求原方程的解。初值選為1.5,容許誤差選為1.0E-5,分
2、析迭代步數(shù)與計算時間,并與上述簡單迭代法作比較.2、 程序設(shè)計的基本思想,原理和算法描述: 算法的描述:1 . 輸入初值,以及各種迭代形式2 .采用for循環(huán)實(shí)現(xiàn)反復(fù)迭代。3 . 以誤差選為1.0E-5,即為終止循環(huán)的條件。4 . 輸出 :最后一個y即為收斂值;y的輸出次數(shù),即為收斂次數(shù)三、主要程序代碼或命令:1.(1)x=1.5;y=(10/x-4*x)(1/2);while abs(y-x)>=1.0e-5x=y;y=(10/x-4*x)1/2endy (2)x=1.5;y=(10/(4+x)(1/2)while abs(y-x)>=1.0e-5x=y;y=(10/(4+x)(
3、1/2)end (3)x=1.5;y=x-(x3+4*x2-10)/(3*x2+8*x)while abs(y-x)>=1.0e-5x=y;y=x-(x3+4*x2-10)/(3*x2+8*x)end 2.(1)x1=1.5;y=(10/x1-4*x1)(1/2);z=y-(10/y-4*y)(1/2);x2=x1-(y-x1)2/(z-2*y+x1)while abs(x2-x1)>=1.0e-5x1=x2;y=(10/x1-4*x1)(1/2);z=y-(10/y-4*y)(1/2);x2=x1-(y-x1)2/(z-2*y+x1)end(2)x1=1.5;y=(10/(4+x
4、1)(1/2);z=(10/(4+y)(1/2);x2=x1-(y-x1)2/(z-2*y+x1)while abs(x2-x1)>=1.0e-5x1=x2;y=(10/(4+x1)(1/2);z=(10/(4+y)(1/2);x2=x1-(y-x1)2/(z-2*y+x1)end(3)x1=1.5;y=x1-(x13+4*x12-10)/(3*x12+8*x1);z=y-(y3+4*y2-10)/(3*y2+8*y);x2=x1-(y-x1)2/(z-2*y+x1)while abs(x2-x1)>=1.0e-5x1=x2;y=x1-(x13+4*x12-10)/(3*x12+8
5、*x1);z=y-(y3+4*y2-10)/(3*y2+8*y);x2=x1-(y-x1)2/(z-2*y+x1)end四、調(diào)試和運(yùn)行程序過程中產(chǎn)生的問題及采取的措施:1的輸入結(jié)果不在區(qū)間上,出現(xiàn)錯誤,經(jīng)檢查,不是程序出現(xiàn)問題;5、 運(yùn)行輸出結(jié)果及分析:1. (1) 的結(jié)果如圖所示, 得出結(jié)論該迭代不收斂。 (2)的結(jié)果如下圖分析: 迭代次數(shù)即y的輸出次數(shù)n = 6收斂于p = 1.365230575673434。(3)分析:形式( 迭代次數(shù)n = 4并且收斂于p = 1.36523001341410。故第三種形式收斂速度更快,時間代價也更少。2. 利用Stenffenson迭代法求原方程的解時形式 (i) 通過編程迭代發(fā)現(xiàn)其收斂,n = 4,收斂于p = 1.365230013196788。形式(ii) 迭代次數(shù)n 3,p =1.365230013414097。形式(iii) 迭代
溫馨提示
- 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)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 農(nóng)村土地整治工作指南
- 2025年黨校入黨積極分子培訓(xùn)班培訓(xùn)結(jié)業(yè)考試試卷及答案(共六套)
- 鄉(xiāng)村振興背景下高職院校新型職業(yè)農(nóng)民工教育的路徑與機(jī)制研究
- 刮膠合同范例
- 2025年安徽省安全員-A證考試題庫及答案
- 乙方解除運(yùn)輸合同范本
- 2025海南省安全員知識題庫及答案
- 農(nóng)村庫房建房合同范本
- 二年級口算題目總匯100道
- 三年級口算題目匯編1000道
- 新能源船舶動力系統(tǒng)的工程實(shí)踐
- SHAFER氣液聯(lián)動執(zhí)行機(jī)構(gòu)培訓(xùn)
- 小學(xué)生守則、日常行為規(guī)范教育實(shí)施方案
- 湖南省六年級上冊數(shù)學(xué)期末試卷(含答案)
- 部編版小學(xué)六年級道德與法治下冊課堂達(dá)標(biāo)檢測試卷全冊含答案
- 巖土工程中的非線性問題分析
- 他們創(chuàng)造了數(shù)學(xué):50位著名數(shù)學(xué)家的故事
- 《普洱茶的定義》課件
- 2023山東經(jīng)貿(mào)職業(yè)學(xué)院教師招聘考試真題題庫
- 《定向運(yùn)動》教學(xué)大綱(含課程思政要素)
- 學(xué)技能如何打逃生繩結(jié)固定繩結(jié)
評論
0/150
提交評論