% 本程序用来计算高度为H的大气压力P、密度rou、温度T，音速a
% 所用字母含义及单位——————————————————————————
% P——压力（Pa）
% rou——密度（kg/m^3）
% T——温度（K）
% a——音速（m/s）
% H——海拔高度（m）
function [P,rou,T,a]=aero(H)
% 原始参数——————————————————————————————
T0=288.15;A=6.5e-3;rou0=1.225;g0=9.80665;R0=6356766;R=287.1;P0=101325;k=1.4;
g=g0*R0^2/(R0+H)^2;tao_on=289.1;G=6.328e-3;A1=230.2;B=6.328e-3;C=1.172e-6;
P1=2.2755e+004;rou1= 0.3659;e=2.71828;
% 相关参量的计算————————————————————————————
if H<=11000
T=T0-A*H;
rou=rou0*(1-A*H/T0)^(-1+g/A/R);
P=P0*(1-A*H/T0)^(g/A/R);
elseif H>11000&&H<=20000
T=216.65;
P=P1*e^(g*(-H+11000)/T/R);
rou=rou1*e^(g*(-H+11000)/T/R);
end
if H<=9300
tao=tao_on-G*H;
elseif H>9300&&H<=12000
tao=A1-B*(H-9300)+C*(H-9300)^2;
elseif H>12000&&H<=30000
tao=211.7;
end
a=sqrt(k*R*tao);
评论0