MATLAB應用在常微分方程式.pptVIP

MATLAB應用在常微分方程式 尋找華航Ci611失事班機 大綱 建立華航Ci611班機的數學模型 飛機水平速度250m/s 度 高度10000公尺 空氣阻力: Ψ = (?·ρair·CD·A)·V2 假設飛機質量為m Ψ = m · a f {F =M*A} 阻力加速度: a f = Ψ / m a f = [(?·ρair·CD·A)/m]· V2 令 (?·ρair·CD·A)/m = B a f = B · V2 ρair = 1.21 ( kg/m3) CD = 0.85 (1/kg) A = 292.55 (m2) m = 374850 (kg) B = 0.0004 ( kg/m) MATLAB 介紹M檔案 ODE範例 華航範例介紹 M 檔案 以MATLAB程式碼所撰寫的檔案通常以 “ m ” 為副檔名,所以這些檔案又稱為M檔案 M檔案又可分為兩類—”主程式”和”函式” 主程式: 是一個副檔名為m的檔案,其內容包含MATLAB的各種指令,我們只要輸入主檔名,即可逐一執行在此檔案內的所有指令 函式: 也是M檔案的一種,可以自行定義一個函式,定義好之後經過存檔,不必執行,當對此函式輸入變數時,程式會自動讀取此已定義的函式 函數對應解法 ode45：4階/5階Rung-Kutta(4,5)法。採單步法，不需附加初始值，速度快，作為解問題的初試法 ode23：2階/3階Rung-Kutta(2,3)法。對於誤差容許範圍較寬且稍帶剛性問題比ode45效果好 ode113：採可變階的Adams-Bashforth-Moulton PECE法，屬於多步解法，比ode45更適合於誤差容許範圍要求較嚴格時 ode23t：採自由內插法實現的梯形規則，適用剛性並要求無數值衰減的問題 ode23s：採用改進的2階Rosenbrock法，採單步法，比ode15s更適於誤差範圍容許較寬的情況。 ode15s：可變階的數值微分公式法(NDFs)，屬於多步解法，如果用在ode45效果很差或失敗，可考慮用此方式。可解剛性問題。 ode23tb：採用TR-BDF2方法，和ode23s一樣，對於誤差容許範圍較寬的情況，比ode15s效果好，可解剛性問題 ODE範例 ODE範例 華航Ci611範例 運算成果 習題 END * * 蔡宏睿 介紹MATLAB如何處理常微分方程式 Boeing 747-200 225 passenger78 combination73 freighter13 convertible4 military393 Total Deliveries 搭載項目 195 ft 8 in (59.6 m) 231 ft 10.2 in (70.6 m)63 ft 5 in (19.3 m) 20 ft (6.1 m) Basic Dimensions??Wing Span (翼寬)??Overall Length (機身長度) ??Tail Height(尾翼高度)??Interior Cabin Width(內部寬度) Mach 0.84555 mph (895 km/h) Typical Cruise Speed (at 35,000 feet) 正規航速 833,000 lb (374,850 kg) Maximum Takeoff Weight 最大載重 747-200 資料來源 : 高宏明 Ψ = (?·ρair·CD·A)·V2 L D A = π·(D/2)2 空氣阻力 圓面積 ρair : 空氣密度 CD : 阻力係數 V : 飛機速度 資料來源 : 高宏明 a f = 0.0004· V2 資料來源 : 高宏明 ρair：空氣密度、CD：阻力係數、A：飛機機頭面積、m：質量 a ax q sin × × = 2 y V B a q cos × × = 2 x V B a 2 V B a f × = 資料來源 : 高宏明 ay 阻力加速度a1=a ax ay q sin × × = 2 y V B - g af q cos × × - = 2 x V B af 垂直飛機加速度 水平飛機加速度 a fx =d Vx / dt = -B V(V cosθ) a fx =d Vx / dt = -B V Vx a fx =d Vx / dt = -B (Vx 2+Vy2) 1/2 Vx 已知 B =0.0004 , g=9.8 a fx =d Vx / dt = -0.0004 (Vx 2+Vy2) 1/2 Vx 已知 B =0.0004 , g=9.