atmosprops.m

来自「用MATLAB编写关于飞机设计的matlab程序」· M 代码 · 共 27 行

   function [Rho, Nu, Acousvel]=ATMOSPROPS(Height)
%	*** Solves for atmospheric density and kinematic viscosity at a 
%        given altitude ***
%
%   *** Real Height, D1, D2, Rho, Nu, Acousvel
%   ***Common /logunits/imon,ikey
%
      
      D2 = 518.7;
      D1 = D2 - 3.56*Height;
      Rho = (2.3769e-3)*(D1/D2)^4.2561;
      if Height>36.0
         D1 = 390.0;
         Rho = (7.1028e-4)*exp(4.7997e-2 * (36.0 - Height));
      end
      Nu = (6.549e-8 + D1*5.9055e-10)/Rho;
      Acousvel = 49.0*sqrt(D1);
%	write(imon,*)
      fprintf( ' **** computed atmospheric properties ****\n');
      fprintf( ' speed of sound (ft/sec) =%6f\n',Acousvel);
fprintf(' density of air (slug/ft3) = %6f\n',Rho);
fprintf( ' viscosity (ft2/sec) =%6f\n ', Nu);
%	write(imon,*)
%      Return
   
%	
%