二阶微分方程式的电脑数值解.PDF

Numerical Methods for Solving the Second-Order Differential Equations (二階微分方程式的電腦數值解) 南台電機 趙春棠 整理 重點1 ：任何一個 N 階 微分方程式，都可以表示為 N 個 1 階 微分方程式 的 聯立 。如 此一來 ，電腦只要會解 1 階 微分方程式 ，電腦 就可以 幫我們解 任何一個 N 階 微分方程 式 了。 重點2 ：利用電腦解 1 階 微分方程式，首先必須了解 EulersMethod(尤拉 法) 。尤拉 法 非 常簡單 ，例如已知初始值，先計算初始值該處的函數微分值 ，就可以估測出下一點 解函數 的 位置了，再以下一點當做起點 ，則可估測下下一點解函數的位置 。由於只利用到一點的函 數微分值，故可視為 Runge-Kutta-1方法。哈 ！由於只利用到一點的函數微分值，故誤差頗大喔 ！ 重點3 ：休恩法(HuensMethod) 改進了尤拉 法。休恩法 除了 採用 目前這一點的函數微分值， 還採用了原本預估點位置的函數微分值喔 ，他最後真正是把兩個函數微分值取平均作為預測斜 率 ，來做估測的喔。由於它利用到兩點的函數微分值，故可視為 Runge-Kutta-2方法。 重點4 ：哈 ！如此說來 ，真正的 Runge-Kutta-4法 ，可知是利用到四個點位置 (簡單來說，起點 1 個，中點取2 個，終點取 1個) 的函數微分值，給它們不同的權重 ，最後以算出的結果作為預 測斜率 ，來做估測的喔 。效果非常準喔 ！  試利用 Euler ’s Method ，撰寫Matlab 程式解聯立方程組 （假定 step size h=0.25 ）  y 1 3y 1 5 y 1 (0) 1  y 2 3y 5 y (0) 2  2 2 【看吧 ！如以上重點1所說 ，一個二階系統變成兩個一階系統的聯立式了 ！ 】 eulerp_2A.m % 歐拉 法 : 解聯立方程組 clear; % initial condition x0 = 0; y10 = 1; yd_1A.m y20 = 2; xmax = 2; h = 0.25; function fxy=yd_1(x,y1,y2); % calculate fxy = 5-3*y1; n = (xmax-x0)/h; for i=1:1:n+1 if i==1 x(i) = x0; y1(i) = y10; y2(i) = y20; yd_2A.m else x(i) = x0+(i-1)*h; % function yd=f(x,y) y1(i) = y1(i-1)+h*yd_1A(x(i-1),y1(i-1),y2(i-1)); function fxy = yd_2(x,y1,y2); y2(i) = y2(i-1)+h*y