徳城余市総合研究所トップページへ戻る
高度（Alt,m）とマッハ数（M）から飛行動圧（DP,kPa）を計算
スクリプトファイル名：AMtoDP.m
-----------------------------------------------------------
# 高度（Alt,m）とマッハ数（M）から飛行動圧（DP,kPa）を計算
# 適用範囲：Alt=0m〜47000m
output_precision(5);
# ーーー定数ーーー
R=287.05287;
g=9.80665;
# ーーー入力ーーー
Alt=input("Input:Alt(m)=")
M=input("Input:M=")
# ーーー計算ーーー
if(Alt<11000)
Ts=288.15-0.0065*Alt;
Ps=101.325*(Ts/288.15)^(-g/(R*(-0.0065)));
elseif(Alt<20000)
Ts=216.65;
Ps=22.632*exp(-g*(Alt-11000)/(R*216.65));
elseif(Alt<32000)
Ts=216.65+0.001*(Alt-20000);
Ps=5.4749*(Ts/216.65)^(-g/(R*0.001));
elseif(Alt<=47000)
Ts=228.65+0.0028*(Alt-32000);
Ps=0.86802*(Ts/228.65)^(-9.80665/(R*0.0028));
endif
#
a=340.294*sqrt(Ts/288.15);
rho=Ps*1000/(R*Ts);
DP=0.5*rho*(a*M)^2/1000;
# ーーー出力ーーー
disp("Output:DP(kPa)=")
disp(DP)
-----------------------------------------------------------
(MM20120722)