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