常微分方程的MATLAB求解實用教案

1、14.1 微分方程(wi fn fn chn)的基本概念微分方程：一般的，凡表示未知函數、未知函數的導數與自變量之間的關系的方程，叫做微分方程，有時也簡稱( jinchng)方程。微分方程的階：微分方程中所出現的未知函數的最高階導數的階數，叫做微分方程的階微分方程的解：找出這樣的函數，把這函數代入微分方程能使該方程成為恒等式。這個函數就叫做微分方程的解。微分方程的通解：如果微分方程的解中含有任意常數，且任意常數的個數與微分方程的階數相同，這樣的解叫做微分方程的通解。初始條件：設微分方程中的未知函數為 ，如果微分方程是一階的，通常用來確定任意常數的條件是 時， 或寫成 其中 都是給定的值；如果微

2、分方程是二階的，通常用來確定任意常數的條件是其中 和 都是給定的值，上述這種條件叫做初始條件。確定了通解中的任意常數以后，就得到微分方程的特解。求微分方程 滿足初始條件 的特解是這樣一個問題，叫做一階微分方程的初值問題，記作微分方程的解的圖形是一條曲線，叫做微分方程的積分曲線。第1頁/共12頁第一頁，共13頁。14.2 幾種(j zhn)常用微分方程類型1.可分離變量的微分方程 一般的，如果一個一階微分方程能寫成的形式，就是說，能把微分方程寫成一端只含 的函數和 ，另一端只含 的函數和 ，那么原方程就稱為可分離變量的微分方程。2.齊次方程 如果一階微分方程可化成 的形式，那么就稱這方程為齊次方

3、程。3.一階線性微分方程 線性方程：方程 叫做一階線性微分方程因為它對于未知函數y 及其導數(do sh)是一次方程。如果 ， 則上述方程稱為齊次的；如果 ， 則上述方程稱為非齊次的。為了求出非齊次線性方程的解，我們先把 換成零而寫出方程 該方程叫做對應于非齊次線性方程的齊次線性方程。齊次線性方程的通解為 非齊次線性方程的通解為伯努利方程：方程 叫做伯努利（Bernoulli）方程。當 時，該方程是線性微分方程，當 時，該方程不是線性的，但是通過變量的替換，便可把它化為線性的第2頁/共12頁第二頁，共13頁。4.可降階的高階微分方程 型的微分方程：微分方程 的右端僅含有自變量 x ，容易看出，

4、只要把 作為新的未知函數，那么微分方程 即化為新未知函數 的一階微分方程，兩邊積分，就得到一個 階的微分方程同理可得依此法繼續(xù)進行，接連積分 n次，便得到方程 的含有 n 個任意常數的通解。 型的微分方程:方程的右端不顯含未知函數 y。如果我們設 ，那么因此，方程 就成為 ，這是一個關于變量 的一階微分方程，設其通解為 ，又 因此又得到一個一階微分方程對它進行積分，便得到方程 的通解為 型的微分方程：方程中不顯含自變量x ，為了求出它的解，我們令 ，并利用復合函數求導法則把 化為對 的導數，即這樣(zhyng)，方程 就成為 這是一個關于變量 的一階微分方程，設它的通解為 分離變量并積分，便得

5、方程 的通解為 第3頁/共12頁第三頁，共13頁。14.3 高階線性微分方程高階線性微分方程(wi fn fn chn)1.線性微分方程解的結構 在 n 階微分方程 中， 若 是 的一次有理整式，則稱此方程為 n 階線性微分方程。一般形式可寫成：線性微分方程解的結構定理： 如果 是方程的n個線性無關的解，則該方程的通解為 其中 是任意常數。 設 是方程 的一個特解， 是對應的齊次線性方程的通解，則是上述方程的通解。 若 和 分別是方程與 的特解，則是方程 的特解2.常系數線性微分方程的MATLAB符號求解 MATLAB中提供了dsolve函數求解微分方程（組）。該函數允許用字符串的形式描述微分

6、方程及初值、邊值條件，最終(zu zhn)將給出微分方程的解析解。第4頁/共12頁第四頁，共13頁。14.4 一階微分方程一階微分方程(wi fn fn chn)初值問題的數值解初值問題的數值解1.歐拉法及其MATLAB實現(shxin) 對于一階微分方程的初值問題 ，若要求其數值解，我們可以采用離散化方法。在求解區(qū)間 上取一組節(jié)點：稱 為步長。為簡單起見，僅考慮等距步長 ，即 將方程 的兩端在區(qū)間 上積分，得到 即 應用左矩形公式 ： ，則有略去上式中的 ，得 考慮到 ，設已求得 ， 的1個近似值 ，則由上式可得 由可依次求出 。稱上式即為求解初值問題的Euler公式。 第5頁/共12頁第五

7、頁，共13頁。2. Runge-Kutta法及其MATLAB實現 考慮微分方程 ，由Lagrange微分中值定理，存在 ，使得(sh de)于是，由 得 記 ，則稱 為區(qū)間 上的平均斜率。這樣，只要給出了 的一種算法，就可以得到求解微分方程初值問題的一種計算公式。顯然，顯式Euler公式就是以 作為平均斜率 的近似。經典四階Runge-Kutta方法的迭代公式：第6頁/共12頁第六頁，共13頁。14.5 一階微分方程一階微分方程(wi fn fn chn)組和高階微組和高階微分方程分方程(wi fn fn chn)的數值解的數值解1. 一階微分方程組 前面研究的是求解單微分方程 的數值解法，對

8、于微分方程組，只需將y 理解成向量， 理解成向量函數，那么對前面研究過的各種計算公式即可用到一階微分方程組上來。2. 高階微分方程 對于高階微分方程組的數值求解，首先應將其變換成一階顯式常微分方程組。其具體轉換方法如下：（1）將微分方程的最高階變量移到等式的左邊，其他移到右邊，并按階次從低到高排列，（這里以兩個高階微分方程的轉換為例）假設兩個高階微分方程最后能夠顯式的表達成下述形式：（2）為每一階微分式選擇狀態(tài)變量，最高階除外（3）根據（2）中選用的狀態(tài)變量，寫出所有狀態(tài)變量的一階微分的表達式最后，對初值進行(jnxng)相應的變換，就可以得到所期盼的一階微分方程組了。3.微分方程組的MATL

9、AB求解函數 MATLAB提供了一系列的函數來求解微分方程組，包括ode系列函數，另外還提供了幾類特殊的微分方程的求解函數，例如ode15s，ode15i等。第7頁/共12頁第七頁，共13頁。14.6 邊值問題的數值邊值問題的數值(shz)解解1.打靶法 打靶法也稱為試射法，其基本思想是把邊值問題作初值問題來求解，從滿足左端邊界條件的解曲線中尋找也滿足右端邊界條件的解。線性方程邊值問題的打靶法：考慮如下給出的二階線性邊值問題該邊值問題的打靶法求解過程可以由如下步驟完成：（1）計算( j sun)下面齊次微分方程在區(qū)間 上的數值解 ， ，初值條件： ；（2）計算( j sun)下面齊次微分方程在

10、區(qū)間 上的數值解 ， ，初值條件： ；（3）計算( j sun)下面初值問題在區(qū)間 上的數值解 ， ，初值條件 ： ；（4）若 ，則計算( j sun) ；（5）計算( j sun)下面初值問題的數值解，則 即為原邊值問題的數值解 ，初值條件：第8頁/共12頁第八頁，共13頁。非線性方程邊值問題的打靶法：考慮二階常微分方程 的邊值問題，邊界條件為 。假定該問題可以轉換為下面的初值問題則問題轉化為求解 ，這是一個復雜(fz)的超越方程，可以考慮引入牛頓迭代法求解參數 m。具體的迭代公式為：式中 通過這些關系可以建立方程具體計算中可以指定一個m 值，然后求解上面的初值問題，將結果代入上面的迭代公式

11、中迭代一步，并將結果代入上式中重新計算，直至兩次計算出來的 m值的誤差在允許的范圍內為止，最后將 m值代入初值問題 即可求解原始問題。第9頁/共12頁第九頁，共13頁。2.邊值問題的MATLAB函數求解 MATLAB能求解的邊值問題的一般形式如下 其中y 為狀態(tài)變量向量 為方程中其他未知參數向量。該方程已知的邊界值為 MATLAB提供了專門(zhunmn)求解邊值問題的bvp解算器bvpslover。要想求解一個常微分方程的邊值問題，一般應該遵循以下幾個步驟：（1）參數初始化（2）微分方程和邊值問題的MATLAB函數描述（3）邊值問題的求解第10頁/共12頁第十頁，共13頁。謝謝(xi xie)大家！第11頁/共12頁第十一頁，共13頁。感謝您的觀看(gunkn)！第12頁/共12頁第十二頁，共13頁。NoImage內容(nirng)總結14.1 微分方