徳城余市総合研究所トップページへ戻る
高度とマッハ数から単位質量あたりの運動エネルギー、位置エネルギー、全エネルギーを計算
Octaveスクリプトファイル名：AMtoKPT.m
-----------------------------------------------------------
# 高度（Alt,m）とマッハ数（M）から単位質量あたりの
# 運動エネルギー（Ek,m^2/s^2）、位置エネルギー（Ep,m^2/s^2）、全エネルギー（Et,m^2/s^2）を計算
# 適用範囲：Alt=0m〜47000m
output_precision(5);
# ーーー定数ーーー
g=9.80665;
# ーーー入力ーーー
Alt=input("Input:Alt(m)=");
M=input("Input:M=");
# ーーー計算ーーー
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
#
a=340.294*sqrt(Ts/288.15);
V=M*a;
Ek=0.5*V^2;
Ep=g*Alt;
Et=Ek+Ep;
# ーーー出力ーーー
disp("Output:Ek(m^2/s^2)=")
disp(Ek)
disp("Output:Ep(m^2/s^2)=")
disp(Ep)
disp("Output:Et(m^2/s^2)=")
disp(Et)
-----------------------------------------------------------
(MM20120808)