📄 standardatmosphere.m
字号:
8903.4, 8809.6, 8715.5, 8622.1, 8528.3, 8434.5, 8340.7, 8246.9, 8153.0,
8059.2, 7965.3, 7871.4, 7777.5, 7683.6, 7589.7, 7495.7, 7401.8, 7307.8,
7213.8, 7119.8, 7025.8, 6931.7, 6837.7, 6743.6, 6649.5, 6555.4, 6461.3,
6367.2, 6364.6, 6365.6, 6366.6, 6367.6, 6368.6, 6369.6, 6370.6, 6371.6,
6373.6, 6375.6, 6377.6, 6379.6, 6381.6, 6411.0, 6442.3, 6473.6, 6504.9,
6536.2, 6567.5, 6598.9, 6630.2, 6661.6, 6692.0, 6724.3, 6755.7, 6831.2,
6915.4, 6999.5, 7083.7, 7252.1, 7420.6, 7589.2, 7757.9, 7926.7, 8042.4,
8047.4, 8004.3, 7845.3, 7686.2, 7527.0, 7367.8, 6969.1, 6569.9, 6244.9,
5961.7, 5678.0, 5649.6, 5621, 5636, 5727, 6009, 7723, 12091, 16288,
20025, 23380, 26414, 29175, 31703, 34030, 36183, 38113, 40043, 43405,
46346, 48925, 51193, 53199, 54996, 56637, 58178, 59678, 63644, 68785,
76427, 88244, 105992, 130630, 161074, 193862, 224737, 250895, 271754,
288203};
AtmosphericPlot[PressureScaleHeight, opts___?OptionQ]:=
ListPlot[
Transpose[{$PressureScaleHeight, $Alts/1000.}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"PressureScaleHeight, Km",
" Altitude, Km"},
PlotRange->All]
$NumberDensity= (* Meter^(-3) *) {4.0151 10^25, 3.8445 10^25, 3.6795 10^25,
3.5201 10^25, 3.3660 10^25, 3.2171 10^25, 3.1017 10^25, 2.9346 10^25,
2.8007 10^25, 2.6715 10^25, 2.5470 10^25, 2.4269 10^25, 2.3113 10^25,
2.2000 10^25, 2.0928 10^25,1.9897 10^25, 1.8905 10^25, 1.7952 10^25,
1.7036 10^25, 1.6156 10^25, 1.5312 10^25, 1.4502 10^25, 1.3725 10^25,
1.2980 10^25, 1.2267 10^25, 1.1585 10^25, 1.0932 10^25, 1.0308 10^25,
9.7110 10^24, 9.1413 10^24, 8.5976 10^24, 8.0790 10^24, 7.5854 10^24,
7.0157 10^24, 6.4857 10^24, 5.9958 10^24, 5.5430 10^24, 5.1244 10^24,
4.7375 10^24, 4.3799 10^24, 4.0493 10^24, 3.4612 10^24, 2.9587 10^24,
2.5292 10^24, 2.1622 10^24,1.8486 10^24, 1.5742 10^24, 1.3413 10^24,
1.1437 10^24, 9.7591 10^23, 8.3341 10^23, 7.1225 10^23, 6.0916 10^23,
5.2138 10^23, 4.4657 10^23, 3.8278 10^23, 3.2833 10^23, 2.8133 10^23,
2.4062 10^23, 2.0558 10^23, 1.7597 10^23, 1.5090 10^23, 1.1158 10^23,
8.3077 10^22, 6.2266 10^22, 4.6965 10^22, 3.5640 10^22, 2.7376 10^22,
2.1351 10^22, 1.6750 10^22, 1.3286 10^22, 1.0488 10^22, 8.2390 10^21,
6.4387 10^21, 3.3934 10^21, 1.7222 10^21, 8.3003 10^20, 3.8378 10^20,
1.7090 10^20, 1.5727 10^20, 1.447 10^20, 7.116 10^19, 2.920 10^19,
1.189 10^19, 2.144 10^18, 5.107 10^17, 1.930 10^17, 9.322 10^16,
5.186 10^16, 3.162 10^16, 2.055 10^16, 1.400 10^16, 9.887 10^15,
7.182 10^15, 5.611 10^15, 4.040 10^15, 2.420 10^15, 1.515 10^15,
9.807 10^14, 6.509 10^14, 4.405 10^14, 3.029 10^14, 2.109 10^14,
1.485 10^14, 1.056 10^14, 4.678 10^13, 2.192 10^13, 1.097 10^13,
5.950 10^12, 3.540 10^12, 2.311 10^12, 1.637 10^12, 1.234 10^12,
9.717 10^11, 7.876 10^11, 6.505 10^11, 5.442 10^11} ;
AtmosphericPlot[NumberDensity, opts___?OptionQ]:=
ListPlot[
Transpose[{Log[10,$NumberDensity], $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"Log[10, NumberDensity]",
"Altitude, Km"},
PlotRange->All]
AtmosphericPlot[GravityAcceleration, opts___?OptionQ]:=
ListPlot[
Transpose[{$GravityAcceleration, $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"GravityAcceleration, m/s^2",
"Altitude, Km"},
Axes->False,
PlotRange->All]
$MeanParticleSpeed=(* (Meter/Second) *) {484.15,481.69, 479.22, 476.73,
474.23, 471.71, 469.69, 466.65, 464.09, 461.53, 458.94, 456.35, 453.74,
451.12, 448.48, 445.82, 443.15, 440.47, 437.76, 435.05, 432.31, 429.56,
426.79, 424.00, 421.20, 418.37, 415.53, 412.67, 409.79, 406.89, 403.97,
401.03, 398.07, 397.95, 397.95,397.95, 397.95, 397.95 , 397.95 , 397.95 ,
397.95 , 397.95 , 397.95 , 397.95 , 397.95 , 397.95, 398.81, 399.72,
400.62, 401.53, 402.43, 403.33, 404.23, 405.12, 406.01, 406.91, 407.79,
408.68, 410.90, 413.35, 415.79, 418.22, 423.03, 427.78, 432.48, 437.13,
441.72, 444.79, 444.79, 443.46, 438.90, 434.29, 429.63, 424.93, 412.95,
400.64, 390.30, 381.05, 371.59, 370.63, 369.7, 369.9, 372.6, 381.4,
431.7, 539.3, 625.0, 691.9, 746.5, 792.2, 831.3, 865.3, 895.1, 921.6,
944.05, 966.5, 1003.2, 1033.5, 1058.7, 1079.7, 1097.4, 1112.4, 1125.5,
1137.4, 1148.5, 1177.4, 1215.0, 1271.5, 1356.4, 1476.0, 1627.0, 1793.9,
1954.3, 2089.6, 2192.6, 2266.4, 2318.1};
AtmosphericPlot[MeanParticleSpeed, opts___?OptionQ]:=
ListPlot[
Transpose[{$MeanParticleSpeed, $Alts/1000}],
opts,
Axes->False,
PlotJoined->True,
Frame->True,
FrameLabel->{"MeanParticleSpeed, m/s",
"Altitude, Km"},
PlotRange->All]
$CollisionFrequency =(* (Second^-1) *) {1.506 10^10, 1.0961 10^10,
1.0437 10^10, 9.9328 10^9, 9.4481 10^9, 8.9824 10^9, 8.6231 10^9,
8.1056 10^9, 7.6934 10^9, 7.2980 10^9, 6.9189 10^9, 6.5555 10^9,
6.2075 10^9, 5.8743 10^9, 5.5554 10^9, 5.2504 10^9, 4.9588 10^9,
4.6802 10^9, 4.4141 10^9, 4.1602 10^9, 3.9180 10^9, 3.6871 10^9,
3.4671 10^9, 3.2577 10^9, 3.0584 10^9, 2.8689 10^9, 2.6888 10^9,
2.5178 10^9, 2.3555 10^9, 2.2016 10^9, 2.0558 10^9, 1.9177 10^9,
1.7871 10^9, 1.6525 10^9, 1.5277 10^9, 1.4123 10^9, 1.3056 10^9,
1.2070 10^9, 1.1159 10^9, 1.0317 10^9, 9.5380 10^8, 8.1528 10^8,
6.9691 10^8, 5.9576 10^8, 5.0931 10^8, 4.3543 10^8, 3.7161 10^8,
3.1733 10^8, 2.7119 10^8, 2.3194 10^8, 1.9852 10^8, 1.7004 10^8,
1.4575 10^8, 1.2502 10^8, 1.0732 10^8, 9.2192 10^7, 7.9251 10^7,
6.8175 10^7, 5.8522 10^7, 5.0297 10^7, 4.3307 10^7, 3.7356 10^7,
2.7939 10^7, 2.1036 10^7, 1.5939 10^7, 1.2152 10^7, 9.3182 10^6,
7.2075 10^6, 5.6201 10^6, 4.3966 10^6, 3.4515 10^6, 2.6961 10^6,
2.0952 10^6, 1.6195 10^6, 8.2945 10^5, 4.0839 10^5, 1.9175 10^5,
8.6559 10^4, 9.8858 10^3, 3.4501 10^3, 3.17 10^3, 1.56 10^3, 6.44 10^2,
2.68 10^2, 5.48 10^1, 1.63 10^1, 7.1 10^0, 3.8 10^0, 2.3 10^0, 1.5 10^0,
1. 10^0, 7.2 10^-1, 5.2 10^-1, 3.9 10^-1, 3.1 10^-1, 2.3 10^-1,
1.4 10^-1, 9.3 10^-2, 6.1 10^-2, 4.2 10^-2, 2.9 10^-2, 2. 10^-2,
1.4 10^-2, 1 10^-2, 7.2 10^-3, 3.3 10^-3, 1.6 10^-3, 8.3 10^-4,
4.8 10^-4, 3.1 10^-4, 2.2 10^-4, 1.7 10^-4, 1.4 10^-4, 1.2 10^-4,
1. 10^-4, 8.7 10^-5, 7.5 10^-6} ;
AtmosphericPlot[CollisionFrequency, opts___?OptionQ]:=
ListPlot[
Transpose[{Log[10, $CollisionFrequency], $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"Log[10, CollisionFrequency]",
"Altitude, Km"},
PlotRange->All,
Axes->False]
$MeanFreePath=(* Meter *) {4.2078 10^-8, 4.395 10^-8, 4.5915 10^-8,
4.7995 10^-8, 5.0193 10^-8, 5.2515 10^-8, 5.4469 10^-8, 5.7571 10^-8,
6.0324 10^-8, 6.3240 10^-8, 6.6332 10^-8, 6.9613 10^-8, 7.3095 10^-8,
7.6795 10^-8, 8.0728 10^-8, 8.4912 10^-8, 8.9367 10^-8, 9.4113 10^-8,
9.9173 10^-8, 1.0457 10^-7, 1.1034 10^-7, 1.1650 10^-7, 1.2310 10^-7,
1.3016 10^-7, 1.3772 10^-7, 1.4583 10^-7, 1.5454 10^-7, 1.6390 10^-7,
1.7397 10^-7, 1.8482 10^-7, 1.9651 10^-7, 2.0912 10^-7, 2.2274 10^-7,
2.4081 10^-7, 2.6049 10^-7, 2.8178 10^-7, 3.0479 10^-7, 3.2969 10^-7,
3.5662 10^-7, 3.8574 10^-7, 4.1723 10^-7, 4.8812 10^-7, 5.7102 10^-7,
6.6797 10^-7, 7.8135 10^-7, 9.1393 10^-7, 1.0732 10^-6, 1.2596 10^-6,
1.4772 10^-6, 1.7312 10^-6, 2.0272 10^-6, 2.3720 10^-6, 2.7734 10^-6,
3.2404 10^-6, 3.7832 10^-6, 4.4137 10^-6, 5.1456 10^-6, 5.9945 10^-6,
7.0212 10^-6, 8.2182 10^-6, 9.6010 10^-6, 1.1196 10^-5, 1.5141 10^-5,
2.0336 10^-5, 2.7133 10^-5, 3.5973 10^-5, 4.7404 10^-5, 6.1713 10^-5,
7.9130 10^-5, 1.0086 10^-4, 1.2716 10^-4, 1.6108 10^-4, 2.0506 10^-4,
2.6239 10^-4, 4.9787 10^-4, 9.8104 10^-4, 2.0354 10^-3, 4.4022 10^-3,
9.8858 10^-3, 1.0743 10^-2, 1.17 10^-2, 2.37 10^-2, 5.79 10^-2,
1.42 10^-1, 7.88 10^-1, 3.31, 8.8, 18, 33, 53, 82, 1.2 10^2, 1.7 10^2,
2.4 10^2,3.3 10^2, 4.2 10^2, 7 10^2, 1.1 10^3, 1.7 10^3, 2.6 10^3,
3.8 10^3, 5.6 10^3, 8. 10^3, 1.1 10^4, 1.6 10^4, 3.6 10^4, 7.7 10^4,
1.5 10^5, 2.8 10^5, 4.8 10^5, 7.3 10^5, 1. 10^6, 1.4 10^6, 1.7 10^6,
2.1 10^6, 2.6 10^6, 3.1 10^6} ;
AtmosphericPlot[MeanFreePath, opts___?OptionQ]:=
ListPlot[
Transpose[{Log[10, $MeanFreePath], $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"Log[10, MeanFreePath]",
"Altitude, Km"},
PlotRange->All,
Axes->False]
$MolecularWeight=(* (Kilo Gram/(Kilo Mole))*) {28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964, 28.964,
28.964, 28.964, 28.964, 28.964, 28.964, 28.95, 28.91, 28.73, 28.40,
27.27, 26.20, 25.44, 24.75, 24.10, 23.49, 22.90, 22.34, 21.81, 21.30,
20.835, 20.37, 19.56, 18.85, 18.24, 17.73, 17.29, 16.91, 16.57, 16.27,
15.98, 15.25, 14.33, 13.09, 11.51, 9.72, 8.00, 6.58, 5.54, 4.85, 4.4,
4.12, 3.94};
AtmosphericPlot[MeanMolecularWeight, opts___?OptionQ]:=
ListPlot[ Transpose[{$MolecularWeight, $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"MolecularWeight, Kg/kmol",
"Altitude, Km"},
PlotRange -> All,
Axes->False]
$Temperature =(* Kelvin *) {320.676, 317.421, 314.166, 310.913, 307.659,
304.406, 301.154, 297.902, 294.651, 291.400, 288.150, 284.900, 281.651,
278.402, 275.154, 271.906, 268.659, 265.413, 262.166, 258.921, 255.676,
252.431, 249.187, 245.943, 242.700, 239.457, 236.215, 232.974, 229.733,
226.492, 223.252, 220.013, 216.774, 216.650, 216.650, 216.650, 216.650,
216.650, 216.650, 216.650, 216.650, 216.650, 216.650, 216.650, 216.650,
216.650, 217.581, 218.574, 219.567, 220.560, 221.552, 222.544, 223.536,
224.527, 225.518, 226.509, 227.500, 228.490, 230.973, 233.473, 236.513,
239.282, 244.818, 250.350, 255.878, 261.403, 266.925, 270.650, 270.650,
269.031, 263.524, 258.019, 252.518, 247.021, 233.292, 219.585, 208.399,
198.639, 188.893, 187.920, 186.87, 186.87, 188.42, 195.08, 240.00,
360.00, 469.27, 559.63, 634.39, 696.29, 747.57, 790.07, 825.16, 854.56,
878.84, 899.01, 929.73, 950.99, 965.75, 976.01, 983.16, 988.15, 991.65,
994.10, 995.83, 998.22, 999.24, 999.67, 999.85, 999.93, 999.97, 999.00,
999.99, 1000, 1000, 1000, 1000};
AtmosphericPlot[KineticTemperature, opts___?OptionQ]:=ListPlot[
Transpose[{$Temperature, $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
Axes->False,
FrameLabel->{"KineticTemperature, K",
"Altitude, Km"},
PlotRange->All]
$Pressure=(* Milli*Bar *) {1.7776 10^3, 1.6848 10^3, 1.5959 10^3, 1.5109 10^3,
1.4297 10^3, 1.3520 10^3, 1.2778 10^3, 1.2069 10^3, 1.1393 10^3,
1.0747 10^3, 1.01325 10^3, 9.5461 10^2, 8.9876 10^2, 8.4559 10^2,
7.9501 10^2, 7.4691 10^2, 7.0121 10^2, 6.5780 10^2, 6.1660 10^2,
5.7752 10^2, 5.4048 10^2, 5.0539 10^2, 4.7217 10^2, 4.4075 10^2,
4.1105 10^2, 3.8299 10^2, 3.5651 10^2, 3.3154 10^2, 3.0800 10^2,
2.8584 10^2, 2.6499 10^2, 2.4540 10^2, 2.2699 10^2, 2.0984 10^2,
1.9299 10^2, 1.7934 10^2, 1.6579 10^2, 1.5327 10^2, 1.4170 10^2,
1.3100 10^2, 1.2111 10^2, 1.0352 10^2, 8.8497 10^1, 7.5652 10^1,
6.4674 10^1, 5.5293 10^1, 4.7289 10^1, 4.0475 10^1, 3.4668 10^1,
2.9717 10^1, 2.5492 10^1, 2.1883 10^1, 1.8799 10^1, 1.6161 10^1,
1.3904 10^1, 1.1970 10^1, 1.0312 10^1, 8.8906 , 7.6730 , 6.6341 ,
5.7459 , 4.9852 ,3.7713, 2.8713, 2.1996, 1.6949, 1.3134, 1.0229,
7.9779 10^-1, 6.2214 10^-1, 4.8337 10^-1, 3.7362 10^-1, 2.8723 10^-1,
2.1958 10^-1, 1.0929 10^-1, 5.2205 10^-2, 2.3881 10^-2, 1.0524 10^-2,
4.4568 10^-3, 4.0802 10^-3, 3.7338 10^-3, 1.8359 10^-3, 7.5966 10^-4,
3.2011 10^-4, 7.1042 10^-5, 2.5382 10^-5, 1.2505 10^-5,7.2028 10^-6,
4.5422 10^-6, 3.0395 10^-6, 2.1210 10^-6, 1.5279 10^-6, 1.1266 10^-6,
8.4736 10^-7, 6.4756 10^-7, 5.0149 10^-7, 3.1059 10^-7, 1.9894 10^-7,
1.3076 10^-7, 8.7704 10^-8, 5.9796 10^-8, 4.1320 10^-8, 2.8878 10^-8,
2.0384 10^-8, 1.4518 10^-8, 6.4468 10^-9, 3.0236 10^-9, 1.4137 10^-9,
8.2130 10^-10, 4.8865 10^-10, 3.1908 10^-10, 2.2599 10^-10,
1.7036 10^-10, 1.3415 10^-10, 1.0873 10^-10, 8.9816 10^-11,
7.5138 10^-11} ;
AtmosphericPlot[Pressure, opts___?OptionQ]:=
ListPlot[
Transpose[{Log[10, $Pressure], $Alts/1000}],
opts,
PlotJoined->True,
Frame->True,
FrameLabel->{"Log[10, Pressure]",
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -