5runge-kuttamethods工程计算常用的方法.doc

計算機程式 微積分與微分方程介紹講義 (三) 目錄 壹、 數值微分解析之資訊與數據 - - 等間距數據 2 一、 原始範例方程式 2 二、 原始方程式轉換之相當數據 2 三、 所摘列4種一階導數之數值方法計算式與計算結果 2 四、 計算程序採用『偽碼法』範例 2 五、 計算程序採用流程圖範例 3 貳、 數值微分解析之資訊與數據 - - 非等間距數據 4 一、 原始範例方程式 4 二、 原始方程式轉換之數據與數據選取方式： 4 三、 利用差分法求取某一點的一次導數值 4 四、 插分表的建立 6 五、 計算程序採用『偽碼法』範例 6 六、 計算程序採用流程圖範例 7 參、 數值積分解析之資訊與數據 8 一、 原始方程式 8 二、 原始方程式轉換之數據 8 三、 數據之函數與自變數關係圖； 8 四、 計算程序採用『偽碼法』範例 9 五、 計算程序採用流程圖範例 10 六、 建議 10 肆、 微分方程解析之資訊與數據（範例一） 12 一、 一階微分方程式： 12 二、 數據： 12 三、 The Taylor-Series Method 12 四、 Euler and Modified Euler Methods 12 五、 Runge-Kutta Methods 12 伍、 微分方程解析之資訊與數據（範例二） 13 一、 一階微分方程式 13 二、 微分方程數值解參考結果： 13 三、 The Taylor-Series Method 13 四、 Euler and Modified Euler Methods 13 五、 Runge-Kutta Methods，工程計算常用的方法。 13 六、 計算程序採用『偽碼法』，以泰勒法為範例 13 七、 計算程序採用流程圖範例 14 數值微分解析一 - - 小等間距數據 原始範例方程式 原始方程式轉換之相當數據 原始等間距數據 N 0 1 2 3 4 5 6 7 8 9 10 11 12 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 5.00 5.68 6.52 7.52 8.68 10.00 11.48 13.12 14.92 16.88 19.00 21.28 23.72 N 13 14 15 16 17 18 19 20 21 22 23 24 25 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 4.6 4.8 5.0 26.32 29.08 32.00 35.08 38.32 41.72 45.28 49.00 52.88 56.92 61.12 65.48 70.00 在 x = 2.4 時，一次導數的預期計算結果 = 4 * x + 3 = 4 * 2.4 + 3 = 12.6 所摘列4種一階導數之數值方法計算式與計算結果 = (26.32 - 23.72) / 0.2 = 13.0 = (26.32 - 21.28) / (2 * 0.2) = 12.6 = (-29.08 + 4 * 26.32 - 3 * 23.72) / (2 * 0.2) = 12.6 = (-29.08 + 8 * 26.32 - 8 * 21.28 + 19.00) / (12 * 0.2) = 12.6 計算程序採用『偽碼法』範例 宣告矩陣：dimension x(50)、f(50)，兩變數都以一維向量宣告。矩陣位址不可使用0，亦即在數據點n = 0的x = 0.0和f = 5.0點，在矩陣是放在x(1)與f(1)的變數中。但在數值方法中，則由 n = 0 為起始做算式推衍。 宣告變數：依需要宣告 open 開啟輸入與輸出指定路徑之檔案 讀入數據： N=26；數據點 Nn=13；計算點 Xx=2.4；計算點的值 x(i) , f(i)； i=1,n 選擇計算點Xx，確定i=Nn值 選擇計算式；依計算準確度選擇 Fd1i = (f(i+1) - f(i))/h； Fd1i = (f(i+1) - f(i-1))/(2*h) ； Fd1i = (-f(i+2) + 4. * f(i+1) – 3. * f(i)) / (2.*h) ； Fd1i = (-f(i+2) + 8. * f(i+1) – 8. * f(i-1) + f(i-2)) / (12. * h) ； 輸出計算結果 計算程序採用流程圖範例 數值微分解析二 - - 寬間距數據 原始範例方程式 原始方程式轉換之數據與數據選取方式：