二分法 牛頓迭代法

二分法 牛頓迭代法

ID:12058103

大?。?5.00 KB

頁數(shù):10頁

時(shí)間:2018-07-15

二分法  牛頓迭代法_第1頁
二分法  牛頓迭代法_第2頁
二分法  牛頓迭代法_第3頁
二分法  牛頓迭代法_第4頁
二分法  牛頓迭代法_第5頁
資源描述:

《二分法 牛頓迭代法》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙2014級碩士研究生數(shù)值分析上機(jī)實(shí)習(xí)(第一次)姓名:喬永亮學(xué)號:14S030125學(xué)院:船舶與海洋工程學(xué)院實(shí)習(xí)題目:分別用二分法和Newton迭代法求方程的根.實(shí)習(xí)目的:掌握兩種解法,體會兩種解法的收斂速度.實(shí)習(xí)要求:用C程序語言編程上機(jī)進(jìn)行計(jì)算,精確到8位有效數(shù)字.報(bào)告內(nèi)容:1.確定實(shí)根的個(gè)數(shù)以及所在區(qū)間.解:對函數(shù)求導(dǎo),得。易知恒成立,所以函數(shù)沒有極值,只有一個(gè)實(shí)根。又可以知道,方程在區(qū)間有一個(gè)實(shí)根,且為奇數(shù)重根,可以二分法和Newton求解2.將最后兩次計(jì)算結(jié)果填入下表(保

2、留8位數(shù)字):方法二分法131.3687744141.3688354Newton法21.368814831.36880823.實(shí)習(xí)過程中遇到哪些問題?如何解決?有何心得體會?在編程的過程中由于對基本計(jì)算原理的理解有一定不足,同時(shí)對編程語言的不熟悉,導(dǎo)致在編程過程中錯(cuò)誤百出,耗費(fèi)了大量時(shí)間。但是通過課本以及網(wǎng)絡(luò)對所需知識的不斷學(xué)習(xí),通過嘗試不同的方法,最終還是得到了幾種不同的思路與方法。通過這次編程,深深的感受到自己的不足,同時(shí)也明白了數(shù)學(xué)與計(jì)算機(jī)編程的緊密結(jié)合,不努力提高自己在當(dāng)今社會就要被淘汰。-10-哈爾濱

3、工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙4.兩種解法的計(jì)算程序(此頁寫不下時(shí)可以加頁):二分法(Fortran語言)programAnalysis1real::a,b,c,mreal::fa,fca=1.b=2.m=0.0001!--------------------dowhile(abs(b-a)>=m)c=(a+b)/2fa=a**3+2.*a*a+10.*a-20fc=c**3+2.*c*c+10.*c-20if(fa*fc<0)thenb=celsea=cendifwrite(*,"(f10.7)")cenddop

4、auseendprogramAnslysis1牛頓迭代法(Fortran語言)programAnalysis2implicitnone!定義變量---------------------------------------------------------------externalf,dfrealm,x0,x1,f,dfintegeri!初始化變量-------------------------------------------------------------m=0.0001x0=1.5!牛頓迭代法

5、-------------------------------------------------------------dowhile(abs(f(x0))>=m)x1=x0-f(x0)/df(x0)x0=x1i=i+1write(*,"(i4,f10.7)")i,x0enddo-10-哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙pauseendprogramAnalysis2!定義原函數(shù)-------------------------------------------------------------functi

6、onf(x)implicitnonerealf,xf=x**3+2.*x*x+10.*x-20.returnendfunction!定義原函數(shù)的導(dǎo)函數(shù)-------------------------------------------------------functiondf(x)implicitnonerealdf,xdf=3.*x*x+4.*x+10.returnendfunction-10-哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙2014級碩士研究生數(shù)值分析上機(jī)實(shí)習(xí)(第二次)姓名:學(xué)號:學(xué)院:實(shí)習(xí)題目:計(jì)算8

7、階三對角矩陣的行列式.實(shí)習(xí)目的:掌握計(jì)算行列式的方法.實(shí)習(xí)要求:首先選擇一種算法,然后用C程序語言編程上機(jī)進(jìn)行計(jì)算.報(bào)告內(nèi)容:1.簡單描述所采用的算法:2.計(jì)算結(jié)果:3.實(shí)習(xí)過程中遇到哪些問題?如何解決?有何心得體會?-10-哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙4.寫出C語言計(jì)算程序(此頁寫不下時(shí)可以加頁):-10-哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙2014級碩士研究生數(shù)值分析上機(jī)實(shí)習(xí)(第三次)姓名:學(xué)號:學(xué)院:實(shí)習(xí)題目:分別用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組實(shí)習(xí)目的:感受兩種迭代法的

8、收斂速度.實(shí)習(xí)要求:首先構(gòu)造收斂的Jacobi迭代法和Gauss-Seidel迭代法,然后用C程序語言編程上機(jī)進(jìn)行求解,初始值均取為0,精確到4位小數(shù).報(bào)告內(nèi)容:1.寫出收斂的Jacobi迭代法和Gauss-Seidel迭代法:-10-哈爾濱工業(yè)大學(xué)(威海)實(shí)驗(yàn)報(bào)告紙2.將最后一次迭代次數(shù)與相應(yīng)的迭代結(jié)果填入下表迭代法JacobiGauss-Seidel3.實(shí)習(xí)過程中遇到哪些問題?如何

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會顯示錯(cuò)亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負(fù)責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時(shí)聯(lián)系客服。
3. 下載前請仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請聯(lián)系客服處理。