standardatmosphere.m

来自「matlab程序」· M 代码 · 共 754 行 · 第 1/3 页

    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]",