資源描述：
《常微分方程 數(shù)值解法》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。
1、第8章 常微分方程實際中，很多問題的數(shù)學模型都是微分方程。我們可以研究它們的一些性質(zhì)。但是，只有極少數(shù)特殊的方程有解析解。對于絕大部分的微分方程是沒有解析解的。常微分方程作為微分方程的基本類型之一，在自然界與工程界有很廣泛的應用。很多問題的數(shù)學表述都可以歸結為常微分方程的定解問題。很多偏微分方程問題，也可以化為常微分方程問題來近似求解。本章討論常微分方程的數(shù)值解法對于一個常微分方程：通常會有無窮個解。如：因此，我們要加入一個限定條件。通常會在端點出給出，如下面的初值問題：為了使解存在唯一，一般，要加限制條件在f上，要求f對y滿足Lipschitz條件：常
2、微分方程的解是一個函數(shù)，但是，計算機沒有辦法對函數(shù)進行運算。因此，常微分方程的數(shù)值解并不是求函數(shù)的近似，而是求解函數(shù)在某些節(jié)點的近似值。例：我們對區(qū)間做等距分割：設解函數(shù)在節(jié)點的近似為由數(shù)值微分公式，我們有，則：向前差商公式可以看到，給出初值，就可以用上式求出所有的基本步驟如下：③解差分方程，求出格點函數(shù)①對區(qū)間作分割：求在上的近似值。稱為分割上的格點函數(shù)②由微分方程出發(fā)，建立求格點函數(shù)的差分方程。這個方程應該滿足：A、解存在唯一；B、穩(wěn)定，收斂；C、相容數(shù)值方法，主要研究步驟②，即如何建立差分方程，并研究差分方程的性質(zhì)。這種方法，稱為數(shù)值離散方法。求的
3、是在一系列離散點列上，求未知函數(shù)y在這些點上的值的近似。我們的目的，就是求這個格點函數(shù)為了考察數(shù)值方法提供的數(shù)值解，是否有實用價值，需要知道如下幾個結論：①步長充分小時，所得到的數(shù)值解能否逼近問題得真解；即收斂性問題②誤差估計③產(chǎn)生得舍入誤差，在以后得各步計算中，是否會無限制擴大；穩(wěn)定性問題8.1Euler公式做等距分割，利用數(shù)值微分代替導數(shù)項，建立差分方程。1、向前差商公式所以，可以構造差分方程稱為局部截斷誤差。顯然，這個誤差在逐步計算過程中會傳播，積累。因此還要估計這種積累定義在假設yi=y(xi)，即第i步計算是精確的前提下，考慮的截斷誤差Ri=y
4、(xi+1)?yi+1稱為局部截斷誤差/*localtruncationerror*/。定義若某算法的局部截斷誤差為O(hp+1)，則稱該算法有p階精度。記為2、收斂性考察局部誤差的傳播和積累稱為整體截斷誤差是1階方法3、穩(wěn)定性－誤差在以后各步的計算中不會無限制擴大。是格式對舍入誤差的抑止作用我們考慮一種簡單情況，即僅初值有誤差，而其他計算步驟無誤差。設是初值有誤差后的計算值，則所以，我們有：可以看出，向前差商公式關于初值是穩(wěn)定的。當初始誤差充分小，以后各步的誤差也充分小4、向后差商公式是隱格式，要迭代求解可以由向前差商公式求出5、中心差商公式是多步，2
5、階格式，該格式不穩(wěn)定6、梯形法－基于數(shù)值積分的公式對微分方程做積分，則：類似，可以算出其誤差估計式：2階的方法所以，有格式為：是個隱式的方法，要用迭代法求解局部截斷誤差8.2Runge－Kutta法由Taylor展開記為所以，可以構造格式這種格式使用到了各階偏導數(shù)，使用不便。從另一個角度看，取(x,y)及其附近的點做線性組合，表示F，問題就好辦了。當然，要求此時的展開精度相同。這種方法稱為Runge－Kutta法在(x,y)處展開，比較以2階為例，設有：1、改進的Euler公式2、Heun公式一般的Runge－Kutta法構造常見的為3階，4階公式8.3
6、線性多步法用若干節(jié)點處的y及y’值的線性組合來近似y(xn+1)。)...(...110111101knknnnknknnnffffhyyyy--+---+++++++++=bbbbaaa其通式可寫為：當??1?0時，為隱式公式;??1=0則為顯式公式。?基于數(shù)值積分的構造法將在上積分，得到只要近似地算出右邊的積分，則可通過近似y(xn+1)。而選用不同近似式Ik，可得到不同的計算公式。若積分用節(jié)點作為積分點，則有積分系數(shù)這是顯格式，q+1階r+1步格式。r=max{p,q}為積分節(jié)點，可以構造r+1步q+1階隱格式局部截斷誤差同樣，若以例：建立p=1,
7、q=2的顯格式p=1，q=2，顯格式，積分區(qū)間為積分節(jié)點為所以例：建立p=2,q=2的隱格式p=2，q=2，隱格式，積分區(qū)間為積分節(jié)點為所以它的截斷誤差較顯格式小，通常也具有更好的穩(wěn)定性。?Adams公式－－p=0時候的多步法參見書§8.4方程組和高階方程的數(shù)值解法寫成向量的形式：各種方法都可以直接運用過來。Euler公式以兩個方程的方程組為例Runge-Kutta公式1、2、確定方法，然后求解（0.202760.0881157）（0.2130070.0934037）（0.2237630.0988499）（0.2350520.104437）（0.2469
8、020.110146）4階Runge-Kutta法，h=1高階方程則有：令例：考
