徳城余市総合研究所トップページへ戻る
高度（Alt,m）と単位質量あたりの全エネルギー（Et,m^2/s^2）からマッハ数（M）を計算
Octaveスクリプトファイル名：AEtoM.m
-----------------------------------------------------------
# 高度（Alt,m）と単位質量あたりの全エネルギー（Et,m^2/s^2）からマッハ数（M）を計算
# 適用範囲：Alt=0m〜47000m
output_precision(5);
# ーーー定数ーーー
g=9.80665;
# ーーー入力ーーー
Alt=input("Input:Alt(m)=");
Et=input("Input:Et(m^2/s^2)=");
# ーーー計算ーーー
if(Alt<11000)
Ts=288.15-0.0065*Alt;
elseif(Alt<20000)
Ts=216.65;
elseif(Alt<32000)
Ts=216.65+0.001*(Alt-20000);
elseif(Alt<=47000)
Ts=228.65+0.0028*(Alt-32000);
endif
#
Ep=g*Alt;
Ek=Et-Ep;
V=sqrt(2*Ek);
a=340.294*sqrt(Ts/288.15);
M=V/a;
# ーーー出力ーーー
disp("Output:M=")
disp(M)
-----------------------------------------------------------
(タイトル訂正：MM20130122)(MM20120810)