徳城余市総合研究所トップページへ戻る

垂直衝撃波上流のマッハ数(M1)を与えて衝撃波通過後のマッハ数(M2)を計算


Octaveスクリプトファイル名:M1toM2.m
-----------------------------------------------------------
# 垂直衝撃波上流のマッハ数(M1)を与えて衝撃波通過後のマッハ数(M2)を計算
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
M1=input("Input:M1=");
#  ーーー計算ーーー
M2=sqrt((1+(k-1)/2*M1^2)/(k*M1^2-(k-1)/2));
# ーーー出力ーーー
disp("M2=")
disp(M2)
-----------------------------------------------------------

(MM20120826)