美国1976年标准大气模型（matlab）

根据1976年的标准大气甲气压函数，返回密度，声音，温度，压力和速度，给定高度的输入高达86公里。 此功能的目的是为那些设计和分析飞机用的，并具有以下的功能，迄今还没有被合并在一个单一标准环境函数 (1)输入可以是标量，矢量，矩阵或高度的n维数组。该功能是向量化和快速计算大气条件在有大量点同时进行。 (2)温度偏移选项为非标准环境中，如分析大热天飞机性能。 (3)输入和输出可以独立地为SI或英制单位。 (4)单位的一致性可以通过使用DimensionedVariable类，减少错误，使代码更清晰强制执行。 (5)返回所需的一切很容易地确定重要参数，如动压，马赫数，雷诺数，滞止温度等

function [Z Z_L Z_U T P rho c g mu nu k n n_sum] = atmo（altdivisionunits）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Program:    1976 Standard Atmosphere Calculator[0-1000 km]
%   Author:     Brent Lewis（RocketLion@gmail.com）
%               University of Colorado-Boulder
%   History:    Original-1/10/2007
%               Revision-1/12/2007-Corrected for changes in Matlab versions
%               for backward compatability-Many thanks to Rich
%               Rieber（rrieber@gmail.com）
%   Input:      alt:        Final Geometric Altitude[km]
%               division:   Reporting points for output arrays[km]
%                           （.01 km & Divisible by .01 km）
%               units:      1-[Metric]
%                           2-{English}
%   Default:    Values used if no input
%               alt:        1000 km
%               division:   1 km
%               units:      Metric
%   Output:     Each value has a specific region that it is valid in with this model
%               and is only printed out in that region
%               Z:          Total Reporting Altitudes[0<=alt<=1000 km][km]{ft}
%               Z_L:        Lower Atmosphere Reporting Altitudes[0<=alt<=86 km][km]{ft}
%               Z_U:        Upper Atmosphere Reporting Altitudes[86<=alt<=1000 km][km]{ft}
%               T:          Temperature array[0<=alt<=1000 km][K]{R}
%               P:          Pressure array[0<=alt<=1000 km][Pa]{in_Hg}
%               rho:        Density array[0<=alt<=1000 km][kg/m^3]{lb/ft^3}
%               c:          Speed of sound array[0<=alt<=86 km][m/s]{ft/s}
%               g:          Gravity array[0<=alt<=1000 km][m/s^2]{ft/s^2}
%               mu:         Dynamic Viscosity array[0<=alt<=86 km][N*s/m^2]{lb/（ft*s）}
%               nu:         Kinematic Viscosity array[0<=alt<=86 km][m^2/s]{ft^2/s}
%               k:          Coefficient of Thermal Conductivity
%                           array[0<=alt<=86 km][W/（m*K）]{BTU/（ft*s*R）}
%               n:          Number Density of individual gases
%                           （N2 O O2 Ar He H）[86km<=alt<=1000km][1/m^3]{1/ft^3}
%               n_sum:      Number Density of total gases
%                           [86km<=alt<=1000km][1/m^3]{1/ft^3}
%   Acknowledgements:       1976 U.S. Standard Atmosphere
%                           Prof. Adam Norris-Numerical Analysis Class
%                           Steven S. Pietrobon USSA1976 Program
%   Notes:                  Program uses a 5-point Simpson‘s Rule in 10
%                           meter increments.  Results DO vary by less 1%
%                           compared to tabulated values and is probably
%                           caused by different integration techniques
%   Examples:               atmo（） will compute the full atmosphere in 1 km
%                           increments and output in Metric Units
%                           atmo（10） will compute 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        6246  2007-01-12 19:49  atmo.m
     文件        3646  2007-01-12 19:35  atmo_compo.m
     文件        2312  2007-01-12 19:36  atmo_p.m
     文件        1736  2007-01-11 04:47  atmo_temp.m
     文件        2893  2007-01-11 04:40  f_n.m
     文件        1108  2007-01-11 04:43  int_tau.m
     文件         645  2007-01-12 19:43  tester.m