徳城余市総合研究所トップページへ戻る
マッハ数（M）と衝撃波角（Beta,deg）を与えて偏角（Theta,deg）を計算
Octaveスクリプトファイル名：BetoTh.m
-----------------------------------------------------------
# マッハ数（M）と衝撃波角（Beta,deg）を与えて偏角（Theta,deg）を計算
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
M=input("Input:M=");
Beta=input("Input:Beta(deg)=");
# ーーー計算ーーー
j=Beta*pi/180;
Th=atan((2*cot(j)*((M*sin(j))^2-1))/(M^2*(k+cos(2*j))+2));
Theta=Th*180/pi
# ーーー出力ーーー
disp("Theta(deg)=")
disp(Theta)
-----------------------------------------------------------
(MM20120818)