%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%????????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