徳城余市総合研究所トップページへ戻る
垂直衝撃波上流の密度(rho1,kg/m^3)とマッハ数(M1)を与えて衝撃波通過後の密度(rho2,kg/m^3)を計算
Octaveスクリプトファイル名：rho1torho2.m
-----------------------------------------------------------
# 垂直衝撃波上流の密度(rho1,kg/m^3)とマッハ数(M1)を与えて衝撃波通過後の密度(rho2,kg/m^3)を計算
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
rho1=input("Input:rho1(kg/m^3)=");
M1=input("Input:M1=");
# ーーー計算ーーー
rho2=rho1*((k+1)/2*M1^2)/(1+(k-1)/2*M1^2);
# ーーー出力ーーー
disp("rho2(kg/m^3)=")
disp(rho2)
-----------------------------------------------------------
(MM20120903)