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

垂直衝撃波上流の速度(V1,m/s)とマッハ数(M1)を与えて衝撃波通過後の速度(V2,m/s)を計算


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

(MM20120901)