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

垂直衝撃波上流の静圧(Ps1,kPa)とマッハ数(M1)を与えて衝撃波通過後の静圧(Ps2,kPa)を計算


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

(MM20120828)