計算機支擋建筑設計方程求解

1、計算機支擋建筑設計方程求解1 1、 通用函數(shù)的程序編寫通用函數(shù)的程序編寫v問題：如何求一個函數(shù)的解，而函數(shù)是可以隨時問題：如何求一個函數(shù)的解，而函數(shù)是可以隨時替換的。替換的。v例如：例如：vF(x)=xF(x)=x3 3+4x+4x2 2+5+5vF(x)=eF(x)=ex x-5-5實現(xiàn)步驟：實現(xiàn)步驟：v1 1在窗體上加一個在窗體上加一個ScriptcontrolScriptcontrol控件控件v2 2 在工程中添加一個類模塊在工程中添加一個類模塊class1.class1.v 類模塊中的代碼：類模塊中的代碼：v Public K1 As Double Public K1 As Doubl

2、e1 1、 通用函數(shù)的程序編寫通用函數(shù)的程序編寫v Private Function Hanshu(X As Double, Str1 As String) As Private Function Hanshu(X As Double, Str1 As String) As DoubleDoublev On Error Resume Next On Error Resume Next v Dim shar As New Class1 Dim shar As New Class1v ScriptControl1.AddObject sh, shar ScriptControl1.AddObjec

3、t sh, sharv ScriptControl1.AddCode Function F(x) & ScriptControl1.AddCode Function F(x) & vbCrLf & kk= & Str1 & vbCrLf & sh.k1 = kk & vbCrLf & kk= & Str1 & vbCrLf & sh.k1 = kk & vbCrLf & End FunctionvbCrLf & End Functionv v ScriptControl1.Run ScriptControl1.Run “F, XF, Xv Hanshu = shar.K1 Hanshu = s

4、har.K1v ScriptControl1.Reset ScriptControl1.Reset v End FunctionEnd Function1 1、 通用函數(shù)的程序編寫通用函數(shù)的程序編寫作業(yè)作業(yè)1 1：試設計一個工程，可以計算函數(shù)的值。試設計一個工程，可以計算函數(shù)的值。2 2 分別輸出分別輸出x x1,2,3101,2,310函數(shù)的值。函數(shù)的值。3 3 測試替換了函數(shù)表達式后，你編寫的程測試替換了函數(shù)表達式后，你編寫的程序是否有效。序是否有效。1 1、 通用函數(shù)的程序編寫通用函數(shù)的程序編寫課題課題2 2 ：高次方程求解：高次方程求解：v 在解決科學研究和工程領域中的實際問題時，往在

5、解決科學研究和工程領域中的實際問題時，往往涉及到高次方程。高次方程解的個數(shù)和方程的往涉及到高次方程。高次方程解的個數(shù)和方程的次數(shù)相對應，即使只在實數(shù)范圍內(nèi)考慮，大多也次數(shù)相對應，即使只在實數(shù)范圍內(nèi)考慮，大多也有幾個解。但是，對于從解決土木工程問題歸納有幾個解。但是，對于從解決土木工程問題歸納得到的高次方程，其有意義的解大多僅有一個，得到的高次方程，其有意義的解大多僅有一個，并且是正數(shù)解。并且是正數(shù)解。2.1 2.1 高次方程求解原理高次方程求解原理v計算機解高次方程的基本過程大致如計算機解高次方程的基本過程大致如下下: :v以以x x2 2-2=0-2=0為例，令為例，令y=xy=x2 2-2

6、-2v其函數(shù)圖像如圖。其函數(shù)圖像如圖。y=0y=0時的時的x x值值(x=(x=土土2 20.50.5) )即為原方程的解。在即為原方程的解。在x =2x =20.50.5附近附近y y值由負變正。由于函數(shù)是連續(xù)的，所值由負變正。由于函數(shù)是連續(xù)的，所以必定存在以必定存在y=0y=0對應的點，該點即為原對應的點，該點即為原方程的一個解。其他方程也是如此，方程的一個解。其他方程也是如此，如果知道了一個解的區(qū)間，就可以通如果知道了一個解的區(qū)間，就可以通過代人不同的過代人不同的x x值觀察值觀察y y的變化，從而的變化，從而逐次逼近該解。逐次逼近該解。高次方程求解常用方法高次方程求解常用方法v高次方程

7、的常用解法有掃描法、對分法、優(yōu)選法、高次方程的常用解法有掃描法、對分法、優(yōu)選法、迭代法、牛頓法等。迭代法、牛頓法等。2.2 2.2 掃描法求解方程的根掃描法求解方程的根v 掃描法的主要過程是掃描法的主要過程是: :v 1. 1.尋找方程解所在的區(qū)域?qū)ふ曳匠探馑诘膮^(qū)域a,ba,b。由初始點。由初始點a a出發(fā)，根據(jù)步長出發(fā)，根據(jù)步長h h逐次迭代找到逐次迭代找到b b點，使點，使f (a ) X f (a ) X (b) = 0 ,(b) = 0 ,前一點即為前一點即為a,f(a)a,f(a)與與f(b)f(b)的值為一正一的值為一正一負。由于從實際土木問題中提煉出的函數(shù)負。由于從實際土木問題

8、中提煉出的函數(shù)f(x)f(x)在在a,ba,b連續(xù)，所以在連續(xù)，所以在a,ba,b內(nèi)必有一解內(nèi)必有一解x0 x0，使，使f(x0)=0f(x0)=0。該過程也可固定。該過程也可固定a a點不變，單純擴大區(qū)點不變，單純擴大區(qū)域找域找b b點點.2.2 2.2 掃描法求解方程的根掃描法求解方程的根 2. 2.縮小解所在的區(qū)域縮小解所在的區(qū)域a,ba,b。十等分。十等分a,ba,b，逐點求，逐點求f(xi)f(xi)值，直至相鄰兩點的值符號相反，這兩點即值，直至相鄰兩點的值符號相反，這兩點即為新的區(qū)域為新的區(qū)域a,ba,b。3.3.重復上述過程直至收斂，達到一定精度為止。重復上述過程直至收斂，達到一

9、定精度為止。2.2 2.2 掃描法求解方程的根掃描法求解方程的根v 注意，求實際土木工程問題的解時，起始點注意，求實際土木工程問題的解時，起始點a a應應取足夠小取足夠小( (小于解小于解) )或干脆取或干脆取a =0a =0。由于其有意義。由于其有意義的解是一個正數(shù)解，搜索方向只要向正向進行，的解是一個正數(shù)解，搜索方向只要向正向進行，甚至根據(jù)實際情況，直接選取適當區(qū)域甚至根據(jù)實際情況，直接選取適當區(qū)域a,ba,b，然，然后從步驟后從步驟2 2開始求解。而一般的高次方程則要復雜開始求解。而一般的高次方程則要復雜些，解不一定比些，解不一定比a a值大，因此解的區(qū)域也可能在初值大，因此解的區(qū)域也可

10、能在初始點的另一邊，搜索方向需先判別。始點的另一邊，搜索方向需先判別。2.2 2.2 掃描法求解方程的根掃描法求解方程的根v掃描法程序流程圖見圖掃描法程序流程圖見圖2-22-2，圖中，圖中h h為步長，為步長，e e為要為要求的精度。求的精度。v掃描法程序中，因為要判斷的是兩函數(shù)的符號，掃描法程序中，因為要判斷的是兩函數(shù)的符號，而不是它們的大小，利用符號函數(shù)可使計算量減而不是它們的大小，利用符號函數(shù)可使計算量減少，即用少，即用m msgn(f(a),n=sgn(f(b)sgn(f(a),n=sgn(f(b)替代原來的替代原來的計算更合理。計算更合理。掃描法求方程根的流程圖掃描法求方程根的流程圖

11、2.2 2.2 掃描法求解方程的根掃描法求解方程的根用掃描法求方程用掃描法求方程x2-2=0 x2-2=0的正數(shù)解，的正數(shù)解，其計算結果見表其計算結果見表2-12-12.2 2.2 掃描法求解方程的根掃描法求解方程的根v由計算結果知，方程由計算結果知，方程x x2 2-2=0-2=0的一個解在的一個解在1. 4141. 414至至1. 1. 415415之間，只要進一步縮小步長之間，只要進一步縮小步長. .重復上述過程，重復上述過程，可獲得更高精度?？色@得更高精度。2.2 2.2 掃描法求解方程的根掃描法求解方程的根v Dim M%, N%, i As Double Dim M%, N%, i

12、 As Doublev Dim A#, B#, H#, E#, X#, Str2$ Dim A#, B#, H#, E#, X#, Str2$v Dim Str1$ Dim Str1$v i = 1 i = 1v Str1 = Str1 = “計算次數(shù)計算次數(shù) a= b= m a= b= mn n & vbCrLf & vbCrLfv A = Val(Text2) A = Val(Text2)初始值初始值v H = Val(Text3) H = Val(Text3)求解步長求解步長v E = Val(Text4) E = Val(Text4)求解精度求解精度v Str2 = Text1.Tex

13、t Str2 = Text1.Text函數(shù)表達式函數(shù)表達式v M = Sgn(Hanshu(A, Str2) M = Sgn(Hanshu(A, Str2)2.2 2.2 掃描法求解方程的根掃描法求解方程的根v Do Dov B = A + H B = A + Hv N = Sgn(Hanshu(B, Str2) N = Sgn(Hanshu(B, Str2)v Str1 = Str1 & i & Space(3) & Str1 = Str1 & i & Space(3) & A & Space(3) & B & Space(3) & M A & Space(3) & B & Space(3)

14、 & M * * N & vbCrLf N & vbCrLfv If M If M * * N 0 Then N 0 Thenv A = B A = Bv M = N M = Nv i = i + 1 i = i + 1v ElseIf M = 0 Then ElseIf M = 0 Thenv X = A X = Av Exit Do Exit Dov ElseIf N = 0 Then ElseIf N = 0 Thenv X = B X = Bv Exit Do Exit Dov v ElseIf H = E Then ElseIf H 5000 Then If i 5000 Thenv

15、 GoTo err1: GoTo err1:v End If End Ifv v Loop Until H = E Loop Until H 0 Then Fb 0 Thenv Text4 = Text4 = “方程根不在你方程根不在你規(guī)定的區(qū)間規(guī)定的區(qū)間 v Exit Sub Exit Subv End If End Ifv If (Fa = 0 Or Fb = 0) Then If (Fa = 0 Or Fb = 0) Then 如果如果FaFa，F(xiàn)bFb0,0,則就是方程的則就是方程的根根v If (Fa = 0) Then If (Fa = 0) Thenv X = A X = Av

16、GoTo 20 GoTo 20v Else Elsev X = B X = Bv GoTo 20 GoTo 20結束程序結束程序v End If End Ifv End IfEnd If2.3 2.3 對分法優(yōu)選法求解對分法優(yōu)選法求解v Strtmp = “計算次數(shù)I= A= B= 中點X= 函數(shù)值 & vbCrLfv i = 1vDo Until Abs(A - B) = Epsv v X = (A + B) / 2v v If Abs(Hanshu(X, Str1) = 0 Then GoTo 20 如果Hanshu(X, Str1)=0，說明此時的x已經(jīng)是方程的解v If Hanshu(

17、X, Str1) * Hanshu(A, Str1) 0 Thenv B = (A + B) / 2v ElseIf Hanshu(X, Str1) * Hanshu(B, Str1) 5000 Then If i 5000 Thenv Text4 = Text4 =“達不到你達不到你所要求的精度，已經(jīng)累加所要求的精度，已經(jīng)累加計算了計算了50005000次次 v Exit Do Exit Dov End If End Ifv LoopLoop2.3.4 2.3.4 對分法關鍵程序?qū)Ψ址P鍵程序2.3 2.3 對分法優(yōu)選法求解對分法優(yōu)選法求解2.3.5 2.3.5 程序界面程序界面2.3 2.

18、3 對分法優(yōu)選法求解對分法優(yōu)選法求解2.3.6 2.3.6 作業(yè)作業(yè)3 3：v本節(jié)已經(jīng)給出了對分的調(diào)用程序，請編寫程序界本節(jié)已經(jīng)給出了對分的調(diào)用程序，請編寫程序界面，完成輸入和輸出。面，完成輸入和輸出。v要求求解時，列出每次求解時的中間過程。要求求解時，列出每次求解時的中間過程。2.3 2.3 對分法優(yōu)選法求解對分法優(yōu)選法求解2.4 2.4 優(yōu)選法優(yōu)選法v 優(yōu)選法又叫優(yōu)選法又叫 0. 618 0. 618法，與對分法流程基本相同，法，與對分法流程基本相同，但每次插人點。不是區(qū)域的中點，而是在區(qū)域的但每次插人點。不是區(qū)域的中點，而是在區(qū)域的0. 0. 618618處，即處，即c=a+0.618(

19、b-a)c=a+0.618(b-a)或或c=b-0.618(b-a),c=b-0.618(b-a),一一般收斂速率更快。般收斂速率更快。v用對分法求解方程用對分法求解方程x x2 2-2=0-2=0時的計算結果見表時的計算結果見表2 22 22.4.1 2.4.1 優(yōu)選法計算表優(yōu)選法計算表2.4 2.4 優(yōu)選法優(yōu)選法2.4.2 2.4.2 作業(yè)作業(yè)v按照對分法求解方法，編寫優(yōu)化法求解方程根的按照對分法求解方法，編寫優(yōu)化法求解方程根的程序，程序，v要求：每一個循環(huán)中間步驟要輸出。要求：每一個循環(huán)中間步驟要輸出。2.4 2.4 優(yōu)選法優(yōu)選法2.5 2.5 迭代法迭代法v 高次方程也可利用迭代方法求

20、解。把原方程高次方程也可利用迭代方法求解。把原方程f(x)f(x)0 0。進行適當變換，建立迭代方程進行適當變換，建立迭代方程xixi1=G(x)1=G(x)依次迭代，依次迭代，直至收斂。直至收斂。v 例例2-12-1求方程求方程x x3 3-x-1 =0-x-1 =0在在x=1.5x=1.5附近的根。附近的根。v 方程由方程由f(x)=0f(x)=0形式改寫成迭代方程形式改寫成迭代方程: :311iixxv 以以1.51.5作作x x的初值的初值x x0 0, ,代入上式求出代入上式求出x x1 11.3572,1.3572,再以再以x x1 1代入，進一步求出代入，進一步求出x x2 2.

21、，直到，直到v 計算結果見下表計算結果見下表(0.00001)0.00001)iiixxx12.5 2.5 迭代法迭代法( (原理）原理）2.5 2.5 迭代法迭代法( (原理）原理）2.5.1 2.5.1 注意事項注意事項v迭代法的關鍵在于選定合適的迭代方程，要求迭迭代法的關鍵在于選定合適的迭代方程，要求迭代方程收斂，否則，因發(fā)散而得不到方程的解。代方程收斂，否則，因發(fā)散而得不到方程的解。如上例，若用如上例，若用x=xx=x3 3-1-1為迭代方程，同樣取為迭代方程，同樣取x0=1.5x0=1.5，則則x1= 2. 375x1= 2. 375，x2= 12.3965.x2= 12.3965.

22、顯然顯然xkxk的值越來的值越來越大，即迭代方程發(fā)散，正因為迭代法有此缺陷越大，即迭代方程發(fā)散，正因為迭代法有此缺陷而限制了它的使用。確定迭代方程是否收斂，可而限制了它的使用。確定迭代方程是否收斂，可以用以用|G(x)|1|G(x)|1進行判別，當進行判別，當x x滿足此式時一般收滿足此式時一般收斂，否則發(fā)散。斂，否則發(fā)散。2.6 2.6 牛頓法牛頓法v 牛頓法又稱弦切法，也是迭代法的一種。牛頓法又稱弦切法，也是迭代法的一種。v 求求f(x) = 0f(x) = 0的解，即求曲線的解，即求曲線f(x)f(x)與橫坐標的交與橫坐標的交點點a(a(如圖如圖2-52-5所示所示) )，在，在a a點

23、附近點附近x x0 0處的斜率處的斜率0100)(0)( xxxfxydxdyxf所以所以)( )(0001xfxfxx)( )(111iiiixfxfxx所以迭代方程為：所以迭代方程為：2.62.6牛頓法牛頓法v結束條件：結束條件：v直到直到f(xf(xi i)=0,)=0,即即x xi i+1=x+1=xi i或或|x|xi+1i+1-x-xi i|為止。為止。2.6.1 2.6.1 牛頓法舉例牛頓法舉例v思考題思考題1 1：v用牛頓迭代法求解用牛頓迭代法求解x x2 2-25=0-25=0的正數(shù)解。的正數(shù)解。v自己用筆先計算。自己用筆先計算。2.6.1 2.6.1 牛頓法舉例牛頓法舉例v

24、 牛頓迭代收斂快，尤其是土木方程，總有一實數(shù)牛頓迭代收斂快，尤其是土木方程，總有一實數(shù)解，且知道解的基本范圍，故使用很方便。解，且知道解的基本范圍，故使用很方便。2.6.3 2.6.3 牛頓迭代法的近似解法牛頓迭代法的近似解法v有時，有時，f(x)f(x)的表達式很復雜或無法求時的表達式很復雜或無法求時v怎么辦呢？怎么辦呢？2.6.3 2.6.3 牛頓迭代法的近似解法牛頓迭代法的近似解法可用數(shù)值近似法替代可用數(shù)值近似法替代: :式中，式中，為遠小于為遠小于xixi的一個小數(shù)。的一個小數(shù)。)()()( iiixfxfxf本節(jié)本節(jié) 上機作業(yè)上機作業(yè)v 4. 4.寫出用寫出用0.6180.618法求

25、解一元三次方程的程序框圖。法求解一元三次方程的程序框圖。按流程求出方程按流程求出方程2 2x3 3+3+3x2 2-17-17x-30=0-30=0在在1,101,10區(qū)域的區(qū)域的解，精度為解，精度為0.01.0.01.v5 5，用迭代法解方程，用迭代法解方程f(x)x-sinx-0.5=0。精度。精度0. 0. 001001v6.6.用牛頓近似迭代法求用牛頓近似迭代法求x2 2+10cos+10cosx=0=0的根的根要求能夠輸出每一步計算的結果。要求能夠輸出每一步計算的結果。本節(jié)本節(jié) 上機作業(yè)上機作業(yè)v牛頓迭代法求解時應注意：作業(yè)牛頓迭代法求解時應注意：作業(yè)7 7v給的初值不同時，求解的解

26、并不同，這時因為如給的初值不同時，求解的解并不同，這時因為如果方程有多個解，它是趨于最近的解。例如作業(yè)果方程有多個解，它是趨于最近的解。例如作業(yè)中中2 2x3 3+3+3x2 2-17-17x-30=0 -30=0 ，分別初值代入，分別初值代入-4,-3,-4,-3,-2,0.1,22,0.1,2分別代入，找到方程的解。分別代入，找到方程的解。思考題思考題v本節(jié)講的掃描法，對分法，優(yōu)選法及迭代法，當本節(jié)講的掃描法，對分法，優(yōu)選法及迭代法，當求到一個合適的解后，程序就自動完成了，求到一個合適的解后，程序就自動完成了，v問題提出，如何求出方程所有的解：問題提出，如何求出方程所有的解：數(shù)值積分數(shù)值積

27、分1 1 定積分的求解定積分的求解v在高等數(shù)學中，對一個定積分在高等數(shù)學中，對一個定積分v的求解，可以利用牛頓一萊布尼茲的求解，可以利用牛頓一萊布尼茲(Newton-(Newton-Leibniz)Leibniz)公式，即公式，即v這里這里F(x)F(x)是是f(x)f(x)的一個原函數(shù)。的一個原函數(shù)。 badxxfI)()()()()(aFbFxFdxxfIbabav 但是，在科學研究和工程技術中常會遇到下述幾種情況，但是，在科學研究和工程技術中常會遇到下述幾種情況，v 如：如：v 1.1.被積函數(shù)的結構復雜，求原函數(shù)困難被積函數(shù)的結構復雜，求原函數(shù)困難v 2.2.原函數(shù)不能用初等函數(shù)表示原

28、函數(shù)不能用初等函數(shù)表示; ;v 3.3.被積函數(shù)不存在原函數(shù)。被積函數(shù)不存在原函數(shù)。v 這時用牛頓一萊布尼茲公式精確計算積分值相當困難，需要這時用牛頓一萊布尼茲公式精確計算積分值相當困難，需要建立積分的近似計算方法。數(shù)值積分就是一種常用的近似計建立積分的近似計算方法。數(shù)值積分就是一種常用的近似計算方法。數(shù)值積分不受被積函數(shù)有無解析解的限制，土木工算方法。數(shù)值積分不受被積函數(shù)有無解析解的限制，土木工程中都可用數(shù)值積分解決。程中都可用數(shù)值積分解決。v 本章介紹最基本的數(shù)值積分法，它們是矩形法、梯形法、辛普森法。本章介紹最基本的數(shù)值積分法，它們是矩形法、梯形法、辛普森法。2 2 矩形積分法矩形積分法

29、v4-14-1矩形矩形v在滿足一定精度的范圍內(nèi)在滿足一定精度的范圍內(nèi). .數(shù)值積分可以只利用被數(shù)值積分可以只利用被積函數(shù)來求得積分值，為什么可以這樣進行呢積函數(shù)來求得積分值，為什么可以這樣進行呢? ?v我們知道，函數(shù)我們知道，函數(shù)f(x)f(x)在區(qū)間在區(qū)間a,ba,b之間求積，就是之間求積，就是求圖求圖4-14-1中的中的f(x)f(x)曲線與曲線與x x軸以及兩直線軸以及兩直線x=a , x x=a , x = b= b所圍成的幾何圖形的面積。所圍成的幾何圖形的面積。2.1 2.1 基本原理基本原理v如圖如圖4-24-2所示，我們把區(qū)間所示，我們把區(qū)間a,ba,b分為分為n n個等分，每個

30、等分，每等分長為等分長為x x，那么其面積，那么其面積10)(niixxfSv作業(yè)作業(yè)1 1：v試編程實現(xiàn)上述矩形積分試編程實現(xiàn)上述矩形積分v要求是個通用程序要求是個通用程序v通過傳遞參數(shù)計算通過傳遞參數(shù)計算2.2 2.2 關鍵程序代碼關鍵程序代碼v Public Function Public Function JuxingJifen(Str1 As String, A JuxingJifen(Str1 As String, A As Double, B As Double, As Double, B As Double, DengFen As Integer) As DoubleDengF

31、en As Integer) As Doublev Dim Sum1 As Double, i As Dim Sum1 As Double, i As Long, H As Double Long, H As Double v H = Abs(B - A) / DengFen H = Abs(B - A) / DengFenv Sum1 = 0 Sum1 = 0v For i = 1 To DengFen For i = 1 To DengFenv Sum1 = Sum1 + Sum1 = Sum1 + Hanshu(A + i Hanshu(A + i * * H, Str1) H, Str

32、1) * * H H v Next Nextv JuxingJifen = Sum1 JuxingJifen = Sum1 v End FunctionEnd FunctionvStr1 Str1 函數(shù)表達式，但必函數(shù)表達式，但必須符合須符合VBVB的表示方法。的表示方法。v A A，B B 積分的上下限積分的上下限vDengfenDengfen積分區(qū)間的等分積分區(qū)間的等分份數(shù)份數(shù)2.3 2.3 精度控制精度控制v在計算時，曲線以下、矩形以上形如三角形的面在計算時，曲線以下、矩形以上形如三角形的面積都未被考慮，故有一定的誤差。當然，若將積都未被考慮，故有一定的誤差。當然，若將x x取得小一些，

33、誤差就可減小，不過計算量卻要大取得小一些，誤差就可減小，不過計算量卻要大大增加，從而也增大了誤差的積累。大增加，從而也增大了誤差的積累。3 3 梯形積分法梯形積分法v減少誤差的一種方法是，將每個等分用梯形代替減少誤差的一種方法是，將每個等分用梯形代替矩形，如圖矩形，如圖4 4 3 3所示。計算方法相應地改為：所示。計算方法相應地改為：xxxfxfSniii02)()(3.1 3.1 梯形積分法舉例梯形積分法舉例v 例例4-14-1利用梯形法求積分利用梯形法求積分的近似解。102dxxv 解解: :該積分的解析解是該積分的解析解是0.333.0.333.。按梯形法，將積。按梯形法，將積分區(qū)間二等

34、分（分區(qū)間二等分（n = 2 ),n = 2 ),5 . 0201x375. 02) 1 ()21(212)21()0(21ffffS3.1 3.1 梯形積分法舉例梯形積分法舉例v隨隨n n增大，增大，x x減小，近似結果趨近于解析解，見減小，近似結果趨近于解析解，見表表xSxS40.250.344160.06250.33480.1260.334640.01560.3333.2 3.2 程序關鍵代碼程序關鍵代碼v但是梯形法以直線但是梯形法以直線來代替曲線仍有誤來代替曲線仍有誤差，如果用某種曲差，如果用某種曲線來代替，就能進線來代替，就能進一步減少誤差，為一步減少誤差，為此而發(fā)展了辛普森此而發(fā)展

35、了辛普森求積方法。求積方法。v Public Function Public Function TixingJifen(Str1 As String, TixingJifen(Str1 As String, A As Double, B As Double, A As Double, B As Double, DengFen As Integer) As DengFen As Integer) As DoubleDoublev Dim Sum1 As Double, i As Dim Sum1 As Double, i As Long, H As Double Long, H As Doubl

36、e v H = Abs(B - A) / DengFen H = Abs(B - A) / DengFenv Sum1 = 0 Sum1 = 0v For i = 1 To DengFen For i = 1 To DengFen v Sum1 = Sum1 + (Hanshu(A Sum1 = Sum1 + (Hanshu(A + (i - 1) + (i - 1) * * H, Str1) + Hanshu(A H, Str1) + Hanshu(A + i + i * * H, Str1) / 2 H, Str1) / 2 * * H Hv Next Nextv TixingJifen = Sum1 TixingJifen = Sum1 v End FunctionEnd Function4 4辛普森法求積辛普森法求積 v辛普森求積公式又叫做三點求積公式，或叫做拋辛普森求積公式又叫做三點求積公式，或叫做拋物線求積公式。如圖所示，將曲線物線求積公式。如圖所示，將曲線f(x)f(x)視為拋物視為拋物線，表示為線，表示為rqxpxxf2)()()22(.)4()2( 2) 12(.)5()3()( 4)(3bfhnafhafhafhnafhafha