输入:高度(m)
输出:密度、当地音速,重力加速度(也可以输出气压、温度)
参考资料:USSA1976
function [row,av,g]=OnAir(YH)
H=6356.766*YH/(6356766+YH);if H<0
H=0;
elseif H<=11
Tb=288.15;
Lb=-6.5;
Pb=101325.3;
Hb=0;
elseif H<=20
Tb=216.65;
Lb=0;
Pb=22631.99;
Hb=11.0;
elseif H<=32
Tb=216.65;
Lb=1;
Pb=5474.798;
Hb=20.0;
elseif H<=47
Tb=228.65;
Lb=2.8;
Pb=868.0101;
Hb=32;
elseif H<=51
Tb=270.65;
Lb=0;
Pb=110.9005;
Hb=47;
elseif H<=71
Tb=270.65;
Lb=-2.8;
Pb=66.93803;
Hb=51;
elseif H<=84
Tb=214.65;
Lb=-2.0;
Pb=4.220234;
Hb=71;
else
row =0;
av=340;
end
T=Tb+Lb*(H-Hb);
if Lb==0
P=Pb*exp(-34.1631947*(H-Hb)/Tb);
else
P=Pb*(Tb/T)^(34.1631947/Lb);
end
row=P/T/287.05;
av=20.0468*sqrt(T);
R=6371004;
g=(6371004/(6371004+YH))^2*9.80665;