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

MATLAB應用在常微分方程式 尋找華航Ci611失事班機 大綱 建立華航Ci611班機的數學模型 介紹MATLAB如何處理常微分方程式 飛機水平速度250m/s 度 高度10000公尺 空氣阻力: Ψ = (?·ρair·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 A = 292.55 (m2) m = 374850 (kg) B = 0.0004 ( kg/m) MATLAB 華航Ci611範例 運算成果 另一推演 END * * 林筱珊 Boeing 747-200 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 : 空氣密度 V : 飛機速度 資料來源 : 高宏明 CD : 阻力係數 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.8 a fy = d Vy / dt = 9.8 – 0.0004 (Vx 2+Vy2) 1/2 Vy a fy = d Vy / dt = g - B V(V sinθ) a fy = d Vy / dt = g - B V Vy a fy = d Vy / dt = g - B (Vx 2+Vy2) 1/2 Vy q sin × × = 2 y V B - g af q cos × × - = 2 x V B af 垂直飛機加速度 水平飛機加速度 a fy = d Vy / dt = 9.8 – 0.0004 (Vx 2+Vy2) 1/2 Vy a fx =d Vx / dt = -0.0004 (Vx 2+Vy2) 1/2 Vx t=0 : Vx = 飛機水平速度250m/s 度 t=0 : Vy = 飛機垂直速度?0m/s 度 資料來源：陳達毅 函式 function dx=ci611(t,v) ax=-0.0004.*(((v(1).^2+v(2).^2).^(1/2)).*v(1)); ay=9.8-0.0004.*((v(1).^2+v(2).^2).^(1/2)).*v(2); dx=[ax;ay]; 註解： % main.m % function dx=ci611(t,v) % ax=-0.0004.*(((v(1).^2+v(2).^2).