matlab实现的标准大气模型

matlab实现的标准大气模型

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%        2014-4-22                %%%%%%
%%%%%%%%%         标准大气模型              %%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % % %温度模型
T0=15;
h1=0:10:11000;
T1=T0-6.5*h1/1000;
h2=11000:100:20000;
T2=-56.5+h2*0;
figure（1）;
axis （[-70 20 0 30000]）;
xlabel（‘T（℃）‘）;
ylabel（‘H（m）‘）;
title（‘图2.1 温度模型‘）;
hold on
plot（T1h1）;
grid ;
plot（T2h2）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % % % 气压模型
P0=101325  ;% 单位 pa
R=287 ; %单位 m^2/（k*s^2） 
g=9.80665; % M/S^2
Tm=216.65; %11000m的标准温度 单位 k
tn1=g/（0.0065*R）;
h1=0:10:11000;
P1=P0*（（1-0.0065*h1/288.15）.^tn1）;
P11=P0*（（1-0.0065*11000/288.15）.^tn1）;
P2=P11*exp（-g*（h2-11000）/（R*Tm））;
figure（2）;
axis （[0 120000 0 30000]）;
xlabel（‘P（pa）‘）;
ylabel（‘H（m）‘）;
title（‘图2.2 气压模型‘）;
hold on
plot（P1h1）;
grid;
plot（P2h2）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % % % 密度模型
h1=0:10:11000;
dst0=1.225 %标准密度 单位kg/m^3
tn2=tn1+1
dst1=dst0*（（1-0.0065*h1/288.15）.^tn2）;
dst11=dst0*（（1-0.0065*11000/288.15）.^tn2）;
dst2=dst11*exp（-g*（h2-11000）/（R*Tm））
figure（3）;
axis （[0 1.4 0 30000]）;
xlabel（‘ρ（kg/m^3）‘）;
ylabel（‘H（m）‘）;
title（‘图2.3 密度模型‘）;
hold on
plot（dst1h1）;
grid;
plot（dst2h2）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1355  2014-05-12 20:13  atmosphere_model.m

----------- ---------  ---------- -----  ----

                 1355                    1