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...