




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、目錄第一章 非線性方程求根····················31.1迭代法··························
2、3;···········31.2牛頓法·····································
3、83;41.3弦截法······································51.4二分法·········&
4、#183;····························6第二章 插值····················
5、···········72.1線性插值····································72.2二次插值
6、183;···································82.3拉格朗日插值············
7、3;···················92.4分段線性插值·····························
8、··102.5分段二次插值·······························11第三章 數值積分·············
9、83;············133.1復化矩形積分法·····························133.2復化梯形積分法····
10、83;························143.3辛普森積分法·······················
11、3;·······153.4變步長梯形積分法···························16第四章 線性方程組數值法···········&
12、#183;······174.1約當消去法·································17 4.2高斯消去法······
13、3;··························184.3三角分解法······················
14、···········204.4雅可比迭代法································214.5高斯賽德爾迭代法···
15、·······················23第五章 常積分方程數值法··················255.1顯示歐拉公式法····
16、3;·························255.2歐拉公式預測校正法······················
17、183;···26 5.3改進歐拉公式法······························275.4四階龍格庫塔法············&
18、#183;···············28數值計算方法第一章 非線性方程求根1.1迭代法程序代碼:Private Sub Command1_Click()x0 = Val(InputBox("請輸入初始值x0")ep = Val(InputBox(請輸入誤差限ep)f = 0While f = 0X1 = (Exp(2 * x0) - x0) / 5If Abs(X1 - x0) < ep ThenPrint X1f
19、= 1Elsex0 = X1End IfWendEnd Sub例:求f(x)=e2x-6x=0在x=0.5附近的根(ep=10-10)1.2牛頓法程序代碼:Private Sub Command1_Click()b = Val(InputBox("請輸入被開方數x0")ep = Val(InputBox(請輸入誤差限ep)f = 0While f = 0X1 = x0 - (x0 2 - b) / (2 * b)If Abs(X1 - x0) < ep ThenPrint X1f = 1Elsex0 = X1End IfWendEnd Sub例:求的值。(ep=10-
20、10)1.3弦截法程序代碼:Private Sub Command1_Click()x0 = Val(InputBox("請輸入第一個初始值x0")X1 = Val(InputBox("請輸入第二個初始值x1")ep = Val(InputBox("請輸入誤差限ep")f = 0While f = 0X2 = X1 - (X1 8 - 13) * (X1 - x0) / (X1 8 - 13) - (x0 8 - 13)If Abs(X2 - X1) < ep ThenPrint X2f = 1Elsex0 = X1X1 = X
21、2End IfWendEnd Sub例:求f(x)=x8-13的正根(初始值x1=1,x2=10,ep=10-10)1.4二分法程序代碼:Private Sub Command1_Click()a = Val(InputBox("請輸入區間端點a")b = Val(InputBox("請輸入區間端點b")ep = Val(InputBox("請輸入誤差限ep")f = 0While f = 0x = (a + b) / 2fx = Exp(-x / 7) * (9 - 2 * x) - 8fa = Exp(-a / 7) * (9 -
22、 2 * a) - 8If fx = 0 Thenf = 1Print "方程的根是", xElseIf fa * fx > 0 Thena = xElseb = xEnd IfIf Abs(b - a) < ep Thenx = (b + a) / 2f = 1Print "方程的根是", xEnd IfEnd IfWendEnd Sub例:求方程f(x)=e-7/x(9-2x)-8在區間0,1內的實根。(ep=10-10)第2章 插 值2.1線性插值程序代碼:Private Sub Command1_Click()X0 = Val(Inp
23、utBox("請輸入第一個結點X:")Y0 = Val(InputBox("請輸入第一個結點Y:")X1 = Val(InputBox("請輸入第二個結點X:")Y1 = Val(InputBox("請輸入第二個結點Y:")f = 0While f = 0x = Val(InputBox("請輸入未知點的自變量值X:")L0 = (x - X1) / (X0 - X1)L1 = (x - X0) / (X1 - X0)y = L0 * Y0 + L1 * Y1Print "x=&quo
24、t; x, "y=" yf = Val(InputBox("是否繼續(0/1):")WendEnd Sub例:已知兩點(13 , 1)、(49 , 8),求30處的值。2.2二次插值程序代碼:Private Sub Command1_Click()X0 = Val(InputBox("請輸入第一個結點X:")Y0 = Val(InputBox("請輸入第一個結點Y:")X1 = Val(InputBox("請輸入第二個結點X:")Y1 = Val(InputBox("請輸入第二個結點Y
25、:")X2 = Val(InputBox("請輸入第三個結點X:")Y2 = Val(InputBox("請輸入第三個結點Y:")f = 0While f = 0x = Val(InputBox("請輸入未知點的自變量值X:")L0 = (x - X1) * (x - X2) / (X0 - X1) / (X0 - X2)L1 = (x - X0) * (x - X2) / (X1 - X0) / (X1 - X2)L2 = (x - X0) * (x - X1) / (X2 - X0) / (X2 - X1)y = L0
26、* Y0 + L1 * Y1 + L2 * Y2Print "x=" x, "y=" yf = Val(InputBox("是否繼續(0/1):")WendEnd Sub例:已知三點(81 ,9)、(100 ,10)、(121 ,10),求98處的值。2.3拉格朗日插值程序代碼:Private Sub Command1_Click()Dim x(), y()n = Val(InputBox("請輸入插值節點數N")ReDim x(n), y(n)For i = 0 To nx(i) = Val(InputBox(&
27、quot;請輸入插值節點x(" + Str(i) + ")")y(i) = Val(InputBox("請輸入插值節點y(" + Str(i) + ")")Next if = 0While f = 0xx = Val(InputBox("請輸入未知點的自變量x:")Sum = 0For i = 0 To nt = 1For j = 0 To nIf j <> i Thent = t * (xx - x(j) / (x(i) - x(j)End IfNext jSum = Sum + t * y
28、(i)Next iPrint "x=" xx, "y=" Sumf = Val(InputBox("是否繼續(0/1)")WendEnd Sub例:已知四點(100 ,10)、(81 ,9)、(64 ,8)、(49 ,7),求87 處的值。2.4分段線性插值程序代碼:Private Sub Command1_Click()Dim x(), y()n = Val(InputBox("請輸入插值節點數N")ReDim x(n), y(n)For i = 0 To nx(i) = Val(InputBox("請
29、輸入插值節點x(" + Str(i) + ")")y(i) = Val(InputBox("請輸入插值節點y(" + Str(i) + ")")Next if = 0While f = 0xx = Val(InputBox("請輸入未知點的自變量x:")L = 0j = 1While L = 0If xx < x(j) Thenk = j + 1L = 1Elsej = j + 1If j > n - 1 Thenk = n - 1L = 1End IfEnd IfWendl0 = (xx -
30、 x(k) / (x(k - 1) - x(k)l1 = (xx - x(k - 1) / (x(k) - x(k - 1)yy = l0 * y(k - 1) + l1 * y(k)Print "x=" xx, "y=" yyf = Val(InputBox("是否繼續(0/1)")WendEnd Sub例:已知三點(361 , 19)、(324 ,18)、(289 ,17),N=2,求300處的值。2.5分段二次插值程序代碼:Private Sub Command1_Click()Dim x(), y()n = Val(Input
31、Box("請輸入插值節點數N")ReDim x(n), y(n)For i = 0 To nx(i) = Val(InputBox("請輸入插值節點x(" + Str(i) + ")")y(i) = Val(InputBox("請輸入插值節點y(" + Str(i) + ")")Next if = 0While f = 0xx = Val(InputBox("請輸入未知點的自變量x:")If x0 < x(1) Thenk = 1f = 1End Ifi = 2Do W
32、hile f = 0 And i >= n - 1If x0 < x(i) ThenIf x0 - x(i - 1) < x(i) - x0 Thenk = i - 1f = 1Elsek = if = 1End IfElsei = i + 1End IfLoopIf f = 0 Thenk = n - 1End Ifl1 = (xx - x(k + 1) * (xx - x(k) / (x(k - 1) - x(k + 1) * (x(k - 1) - x(k)l2 = (xx - x(k + 1) * (xx - x(k - 1) / (x(k) - x(k + 1) *
33、 (x(k) - x(k - 1)l3 = (xx - x(k) * (xx - x(k - 1) / (x(k + 1) - x(k) * (x(k + 1) - x(k - 1)yy = l1 * y(k - 1) + l2 * y(k) + l3 * y(k + 1)Print "x=" xx, "y=" yyf = Val(InputBox("是否繼續(0/1)")WendEnd Sub例:已知三點(225 , 15)、(196 ,14)、(169 ,13),求180處的值。第3章 數值積分3.1 復化矩形積分法程序代碼:Pr
34、ivate Sub Command1_Click()a = Val(InputBox("請輸入積分下限a")b = Val(InputBox("請輸入積分上限b")n = Val(InputBox("請輸入積分區間等分數N")h = (b - a) / nSum = 0For i = 1 To nSum = Sum + (a + (i - 0.5) * h) 2Next ir = h * SumPrint "復化矩形積分法計算結果:" rEnd Sub 例:求X2在區間3,5的積分值,等分區間為100。3.2復化梯
35、形積分法程序代碼:Private Sub Command1_Click()a = Val(InputBox("請輸入積分下限a")b = Val(InputBox("請輸入積分上限b")n = Val(InputBox("請輸入積分區間等分數N")h = (b - a) / nSum = 0For i = 1 To n - 1Sum = Sum + (a + i * h) 2Next it = h * (a 2 + b 2) / 2 + h * SumPrint "復化梯形積分法計算結果:" tEnd Sub例:求
36、X2在區間3,5的積分值,等分區間為100。3.3辛普生積分法程序代碼:Private Sub Command1_Click()a = Val(InputBox("請輸入積分下限A:")b = Val(InputBox("請輸入積分上限B:")n = Val(InputBox("請輸入積分區間等分數N:")h = (b - a) / nSum = 0w = a + h / 2While w < bSum = Sum + (w - h / 2) 2 + 4 * w 2 + (w + h / 2) 2w = w + hWends =
37、 Sum * h / 6Print "辛普生積分法計算結果:" sEnd Sub例:求X2在區間3,5的積分值,等分區間為100。3.4變步長梯形積分法程序代碼:Private Sub Command1_Click()a = Val(InputBox("請輸入積分下限a:")b = Val(InputBox("請輸入積分上限b:")ep = Val(InputBox("請輸入誤差限ep:")n = 1h = b - at1 = h * (a 2 + b 2) / 2f = 0While f = 0Sum = 0Fo
38、r i = 1 To n - 1Sum = Sum + (a + (i - 1 / 2) * h) 2Next iT2 = t1 / 2 + h * Sum / 2If Abs(T2 - t1) < ep ThenPrint T2f = 1Elset1 = T2h = h / 2n = 2 * nEnd IfWendEnd Sub例:求X2在區間3,5的積分值,誤差限ep=0.0001。第4章 線性方程組數值解法4.1 約當消去法程序代碼:Private Sub Command1_Click()Dim a()n = Val(InputBox("請輸入方程的個數N:")
39、ReDim a(n, n + 1)For i = 1 To nFor j = 1 To n + 1a(i, j) = Val(InputBox("請輸入增廣矩陣A(" + Str(i) + "," + Str(j) + ")=")Next jNext iFor k = 1 To nm = a(k, k)For j = k To n + 1a(k, j) = a(k, j) / mNext jFor i = 1 To nIf i <> k Thenm = a(i, k)For j = k To n + 1a(i, j) =
40、a(i, j) - a(k, j) * mNext jEnd IfNext iNext kFor i = 1 To nPrint "x(" + Str(i) + ")=" a(i, n + 1)Next iEnd Sub例:已知方程組,求X1 X2 X3的值。 X1+ X2+ X3=20 2X1+3X2+ X3=38 X1+4X2+ X3=35 4.2 高斯消去法程序代碼:Private Sub Command1_Click()Dim a(), x()n = Val(InputBox("請輸入方程個數N:")ReDim a(n, n
41、+ 1), x(n)For i = 1 To nFor j = 1 To n + 1a(i, j) = Val(InputBox("請輸入增廣矩陣A(" + Str(i) + "," + Str(j) + "):")Next jNext iFor k = 1 To n - 1m = a(k, k)For j = k To n + 1a(k, j) = a(k, j) / mNext jFor i = k + 1 To nm = a(i, k)For j = k To n + 1a(i, j) = a(i, j) - a(k, j) *
42、 mNext jNext iNext kx(n) = a(n, n + 1) / a(n, n)For i = n - 1 To 1 Step -1Sum = 0For j = i + 1 To nSum = Sum + a(i, j) * x(j)Next jx(i) = a(i, n + 1) - SumNext iFor i = 1 To nPrint "x(" + Str(i) + ")=" x(i)Next iEnd Sub例:已知方程組,求X1 X2 X3的值。2X1+5X2+ X3=21 4X1+3X2+3X3=31 X1+3X2+2X3=
43、16 4.3 三角分解法程序代碼:Private Sub Command1_Click()Dim a(), b(), l(), u(), x(), y()n = Val(InputBox("請輸入方程個數N:")ReDim a(n, n), b(n), l(n, n), u(n, n), x(n), y(n)For i = 1 To nFor j = 1 To na(i, j) = Val(InputBox("請輸入系數矩陣A(" + Str(i) + "," + Str(j) + "):")Next jb(i)
44、= Val(InputBox("請輸入右端常數項B(" + Str(i) + "):")Next iFor i = 1 To nFor j = 1 To nIf i > j ThenSum = 0For k = 1 To j - 1Sum = Sum + l(i, k) * u(k, j)Next kl(i, j) = (a(i, j) - Sum) / u(j, j)ElseSum = 0For k = 1 To i - 1Sum = Sum + l(i, k) * u(k, j)Next ku(i, j) = a(i, j) - SumEnd
45、IfNext jNext iFor i = 1 To nSum = 0For j = 1 To i - 1Sum = Sum + l(i, j) * y(j)Next jy(i) = b(i) - SumNext iFor i = n To 1 Step -1Sum = 0For j = i + 1 To nSum = Sum + u(i, j) * x(j)Next jx(i) = (y(i) - Sum) / u(i, i)Next iFor i = 1 To nPrint "x(" + Str(i) + ")=" x(i)Next iEnd Sub
46、例:已知方程組,求X1 X2 X3的值。2X1+3X2+4X3=38 X1+7X2+3X3=65 5X1+2X2+ X3=33 4.4 雅可比迭代法程序代碼:Dim a(), b(), x0(), x()n = Val(InputBox("請輸入方程個數N:")ReDim a(n, n), b(n), x0(n), x(n)ep = Val(InputBox("請輸入誤差限EP:")nmax = Val(InputBox("請輸入最大迭代次數Nmax:")For i = 1 To nFor j = 1 To na(i, j) = Va
47、l(InputBox("請輸入系數矩陣A(" + Str(i) + "," + Str(j) + "):")Next jb(i) = Val(InputBox("請輸入右端常數項B(" + Str(i) + "):")Next iFor i = 1 To nx0(i) = Val(InputBox("請輸入初始解X0(" + Str(i) + "):")Next ik = 0f = 0While f = 0Max = 0For i = 1 To nSum
48、= 0For j = 1 To nSum = Sum + a(i, j) * x0(j)Next jd = (b(i) - Sum) / a(i, i)x(i) = x0(i) + dIf Max < Abs(d) ThenMax = Abs(d)End IfNext ik = k + 1If Max < ep Then For i = 1 To n Print "x(" + Str(i) + ")=" x(i) Next i f = 1Else If k > nmax Then Print "迭代失敗" f = 1
49、 Else For i = 1 To n x0(i) = x(i) Next i End IfEnd IfWendEnd Sub例:已知方程組,求X1 X2 X3 X4的值。(ep=0.0001,最大迭代次數1000,初始值均為0)7X1+2X2+3X3+ X4=5 X1+5X2+2X3+ X4=42X1+ X2+8X3+3X4=2 X1+4X2+ X3+7X4=74.5 高斯-賽德爾迭代法程序代碼:Private Sub Command1_Click()Dim a(), b(), x()n = Val(InputBox("請輸入方程個數N:")ep = 0.0001nma
50、x = 10000ReDim a(n, n), b(n), x(n)For i = 1 To nFor j = 1 To na(i, j) = Val(InputBox("請輸入系數矩陣A(" + Str(i) + "," + Str(j) + "):")Next jb(i) = Val(InputBox("請輸入右端常數項b(:" + Str(i) + "):")Next iFor i = 1 To nx(i) = Val(InputBox("請輸入初始解X(:" + St
51、r(i) + "):")Next ik = 0f = 0While f = 0Max = 0For i = 1 To nSum = 0For j = 1 To nSum = Sum + a(i, j) * x(j)Next jd = (b(i) - Sum) / a(i, i)x(i) = x(i) + dIf Max < Abs(d) ThenMax = Abs(d)End IfNext ik = k + 1If Max < ep ThenFor i = 1 To nPrint "x(" + Str(i) + ")="
52、x(i)Next if = 1ElseIf k > nmax ThenPrint "迭代失敗"f = 1End IfEnd IfWendEnd Sub例:已知方程組,求X1 X2 X3的值。(初始值均為0)5X1+ X2+2X3=2 X1+3X2+ X3=8 3X1+ X2+7X3=6第5章 常微分方程的數值解5.1顯示歐拉法程序代碼:Private Sub Command1_Click()x0 = Val(InputBox("請輸入自變量初始值x0")y0 = Val(InputBox("請輸入自變量初始值y0")xn = Val(InputBox("請輸入自變量值終值xn")n = Val(InputBox("請輸入區間等分數")h = (xn - x0) / nX1 = x0 + hIf
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 農業種植確權管理辦法
- 高層建筑火災模擬與人員安全疏散策略研究
- 教育的進階之路:學校改進策略與實踐探索
- 普通設備租賃管理辦法
- 液氯企業安全風險隱患排查表
- 景區物業收費管理辦法
- 電氣工程雙創人才培養模式探討與實踐
- 重點高校自學考試課程體系優化研究
- 幕墻工程工作總結
- 高校數字化資源服務系統用戶體驗優化
- GB/T 8312-2002茶咖啡堿測定
- 2023年蘇州國發創業投資控股有限公司招聘筆試題庫及答案解析
- 通信線路工程施工組織設計方案【實用文檔】doc
- 護士注冊健康體檢表下載【可直接打印版本】
- 預計財務報表編制及分析課件
- 學生集體外出活動備案表
- Q∕SY 1347-2010 石油化工蒸汽透平式壓縮機組節能監測方法
- 西門子順序功能圖語言S7-Graph的應用
- 中醫治療室工作制度管理辦法
- 提花裝造工藝技術培訓課程
- 直播傳媒公司簡介PPT課件(參考)
評論
0/150
提交評論