2019級碩士研究生數值分析上機實習報告

2021-10-05 01:42:57 字數 1741 閱讀 6552

2015級碩士研究生數值分析上機實習 (第一次)姓名:李友龍學號:15s030015 學院:

船舶與海洋工程學院實習題目:分別用二分法、newton迭代法和割線法求方程的實根.

實習目的:掌握3種解法,體會3種解法的收斂速度.

實習要求:用c程式語言程式設計上機進行計算,精確到8位有效數字.

報告內容:

1. 確定實根的個數以及所在區間.

令,則,由可得從而有在(,)和(0,)上單調增加,在()上單調遞減。,,而;。因此,在(0,)上存在唯一個零點,又=-5, =14.故在(1,2)內有且僅有乙個根。

2. 將最後兩次計算結果填入下表:

3. 實習過程中遇到哪些問題?如何解決?有何心得體會?

在程式編寫過程中對有效數字和輸出結果的形式比較模糊,導致精度出現過小,但是在輸出形式上卻沒有反應出來,經過和同學的討論,以及不斷地改變精度和輸出結果的形式最終得到解決,並且體會到了計算機結果中精度的重要性,從這三種解法來看,收斂速度不一樣,newton法收斂最快,割線法次之,二分法最慢,從編寫的程式來看,newton法和割線法程式相對來說要簡單,綜之,newton法和割線法要優於二分法。

4. 在下面分別寫出3種解法的計算程式(此頁寫不下時可以加頁):

二分法程式:

program erfenfa

real(8)::a=1,b=2,d=5e-8real(8) c,m,n

integer i

do i=1,50

c=(a+b)/2

m=c*c*c+4*c*c-10

if(m==0) then

write(*,*)i

write(*,1)c

stop

else

if(b-c<=d) then

write(*,*)i

write(*,1)c

stop

else

n=a*a*a+4*a*a-10

if(m*n<0) then

b=celse

a=cend if

end if

end if

write(*,*)i

write(*,1)c

end do

1 format("c=",f10.8)

end program

newton:

program newton

implicit none

real(8)::x=1.5,y

integer i

do i=1,50

y=x-(x*x*x+4*x*x-10)/(3*x*x+8*x)do while(abs(y-x)>5e-8)x=y

end do

write(*,*)i

write(*,1)y

end do

1 format("y=",f10.8)

end program

gexianfa:

program gexianfa

real(8)::a=1,b=2,c

integer i

do i=1,50

c=b-(b-a)*(b*b*b+4*b*b-10)/(b*b*b+4*b*b-a*a*a-4*a*a)

do while(abs(c-b)>5e-8)a=b

b=cend do

write(*,*)i

write(*,1)c

end do

1 format("c=",f10.8)

end program

2019級碩士研究生數值分析上機實習報告

姓名 李友龍學號 15s030015 學院 船舶與海洋工程 實習題目 分別用復化梯形公式和復化simpson公式計算積分 的近似值和.實習目的 體會兩種復化求積公式的收斂性與收斂速度.實習要求 用c程式語言程式設計上機進行計算,結果要有八位有效數字.報告內容 1.寫出求和的復化求積公式 2.寫出具有...

2019級碩士研究生數值分析上機實習報告 答案

姓名學號學院 實習題目 分別用二分法和newton迭代法求方程的根.實習目的 掌握兩種解法,體會兩種解法的收斂速度.實習要求 用c程式語言程式設計上機進行計算,精確到8位有效數字.報告內容 1.確定實根的個數以及所在區間.2.將最後兩次計算結果填入下表 保留8位數字 3.實習過程中遇到哪些問題?如何...

2019級化工類碩士研究生現代數值分析試題全集

一 選擇題 簡述題 1 具有4位有效數字的近似值為 d a 3.1415 b.3.142 c 3.141 d.3.1416.2.若,求值的比較準確的算式為 c a.b.cd.3.使用gauss消去法求解乙個n元線性方程組 ax b所需乘 除法 運算次數約為 b a.b.c.d.4.若實方陣a滿足 d...