作业vanderWaalsequation体积求解.doc

數值微分解析之資訊與數據 - - 原始方程式： 數據：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 N 26 27 28 29 30 5.2 5.4 5.6 5.8 6.0 74.68 79.52 84.52 89.68 95.00 在x = 2.4時，一次導數的預期計算結果 = (26.32 - 21.68) / (2.6 - 2.2) = 5.04 / 0.4 = 12.6 = 4 * x + 3 = 4 * 2.4 + 3 = 12.6 第一種選擇，以求取x=2.4的一階導數 第二種選擇，以求取x=2.4的一階導數 0 0 5；(1,1) 0 0 5；(1,1) 1 1 10；(2,1) 1 2 19；(2,1) 2 2 19；(3,1) 2 3 32；(3,1) 3 3 32；(4,1) 3 4 49；(4,1) 4 4 49；(5,1) 4 5 70；(5,1) 5 5 70；(6,1) 5 6 95；(6,1) 利用差分法求取某一點的一次導數值 有一個非線性函數可表示成n階的多項式表示利用泰勒級數展開的觀念，但採用分離差分法處裡，其表示式為：級數展開適用於等間距，分離差分法則可適用於不等間距數據的系統。 若可用合理的表示，所以一階導數只要對作一次微分即可。 其中皆為分離差分所得的數值函數的數值：一次微分的數值：；下標(0)表示其起算點為二次微分的數值：三次微分的數值： 各多項式為多階函數，為差分項相乘之之運算，由，所以高階項微分之結果為： 計算一次導數之使用通式一 計算一次導數之使用通式二 計算使用式三，應用於本範例的計算 使用第 1、2、3 三個點 使用第 2、3、4 三個點 插分表的建立 選擇一 0 0 5 (10-5)/(1-0) = 5 (9-5)/(2-0) = 2 0 0 0 1 1 10 (19-10)/(2-1) = 9 (13-9)/(3-1) = 2 0 0 2 2 19 (32-19)/(3-2) = 13 (17-13)/(4-2) = 2 0 3 3 32 (49-32)/(4-3) = 17 (21-17)/(5-3) = 2 4 4 49 (70-19)/(5-4) = 21 5 5 70 2.4 使用第 0、1、2 點的計算結果： 12.6 2.4 使用第 1、2、3 點的計算結果： 12.6 2.4 使用第 2、3、4 點的計算結果： 12.6 2.4 使用第 3、4、5 點的計算結果： 12.6 2.4 使用第 4、5點的計算結果： 21 選擇二 0 0 5 (19-5)/(2-0) = 7 (13-7)/(3-0) = 2 0 0 0 1 2 19 (32-19)/(3-2) = 13 (17-13)/(4-2) = 2 0 0 2 3 32 (49-32)/(4-3) = 17 (21-17)/(5-3) = 2 0 3 4 49 (70-49)/(5-4) = 21 (25-21)/(6-4) = 2 4 5 70 (95-70)/(6-5) = 25 5 6 95 2.4 使用第 0、1、2 點的計算結果： 12.6 2.4 使用第 1、2、3 點的計算結果： 12.6 2.4 使用第 2、3、4 點的計算結果： 12.6 2.4 使用第 3、4、5 點的計算結果： 12.6 2.4 使用第 4、5點的計算結果： 25 計算程序採用『偽碼法』範例 宣告矩陣：x(10)、f(10,