作业：van der Waals equation体积求解.DOC

計算機程式 數學方法介紹講義 () 目錄 壹、 Least-squares Method最小二平方法 2 一. 數據 2 二. 作圖結果 2 三. 數學模式之係數的計算之推導 2 四. 設計計算程式 3 貳、 van der Waals 立方型狀態方程式體積求解 4 一. 問題解析 4 二. 體積求解之數學表示式 4 三. 可利用下列之試算，找出根的坐落範圍 5 參、 The Bisection method；Halving the interval method求解 6 一. 計算程序解析 6 二. 解題的邏輯觀念示意圖 6 三. 流程圖範例一 7 四. 流程圖範例二 8 肆、 Linear interpolation methods; The Secant method求解 9 一. 計算程序解析 9 二. 解題的邏輯觀念示意圖 9 三. 流程圖範例 10 伍、 Linear interpolation methods; The False position method求解 11 一. 計算程序解析 11 二. 解題的邏輯觀念示意圖 11 三. 流程圖範例 12 陸、 Newton’s methods求解 13 一. 計算程序解析 13 二. 解題的邏輯觀念示意圖 13 三. 流程圖範例 14 柒、 Fixed-point iteration求解 15 一. 計算程序解析 15 二. 流程圖範例 15 Least-squares Method最小二平方法 數據 0.05 0.11 0.15 0.31 0.46 0.52 0.70 0.74 0.82 0.98 1.17 0.956 0.890 0.832 0.717 0.571 0.539 0.378 0.370 0.306 0.242 0.104 作圖結果 數學模式之係數的計算之推導 數據變數有兩個：自變數()與因變數() 函數關係 或 實驗數據與模式理論數值間有誤差取誤差的平方和為 當期望誤差的平方和為最小時，在微積分中，函數一次微分等於0為極值。數學模式的變數此時為a與b。 整理上兩個一次微分式 再整理i (1) (2) 各項計算數值 項目 欄總計值 6.01 5.905 4.6545 2.1839 36.1201 35.48905 第（1）式乘以N 減去 第（2）式乘以 可解得 a 的算式（或者數值） a的算式帶入第(2)式可解得b 的算式（或者數值） 設計計算程式 回目錄 van der Waals 立方型狀態方程式體積求解 問題解析 van der Waals equation體積立方型方程式之各種類型： 立方型狀態方程式之物質參數a、b的取得： 在臨界點時 ，所以 在臨界點時，van der Waals equation 比較係數： 、 、 整理，用臨界壓力與臨界溫度表示： 、 體積求解之數學表示式 嘗試錯誤法（Try and error method）：，須採用嘗試錯誤法求解；在數值方法中，通常有下列數種嘗試錯誤之迫近法尋求適當解。 The Bisection method Linear interpolation methods; the Secant method Linear interpolation methods; Linear Interpolation, (False position) Newton’s method Muller’s method 疊代法；Fixed-point iteration: x = g(x) method ：，須採用疊代法求解。 ：，須採用疊代法求解。 ：，所以 ，須採用疊代法求解。 ：，所以 ，須採用疊代法求解。 可利用下列之試算，找出根的坐落範圍 V f(P,V,T) V f(P,V,T) 30000 1,369,259,433,207 1000 (11,739,485,707) 29000 456,303,568,846 900 (8,339,498,556) 28000 (340,842,637,515) 800 (5,467,414,826) 27000 (1,028,179,185,878) 700 (3,129,234,515) 26000 (1,611,706,076,242) 600 (1,330,957,624) 25000 (2,097,423,308,607) 500 (78,584,154) 24000 (2,491,330,882,973) 400 621,885,897 23000 (2,799,428,799,341)