Простая функция определения плотности воздуха
Простая функция аппроксимирующая плотность воздуха в диапазоне высот от 0 до 100 км
% Модель атмосферы в диапазоне от 0 до 100 км с погрешностью 1,5%
% В. А. ЯРОШЕВСКИЙ "АППРОКСИМАЦИЯ МОДЕЛИ СТАНДАРТНОЙ АТМОСФЕРЫ"
% Ученые записки ЦАГИ, т. XL, №3, 2009.
%
% density = rho(h)
% Atmospheric density
% h [meters]
% density [kg/m^3]
function density = rho(h)
a = [-6.3759,-7.3012,-1.1817];
b = [-0.4754,-0.0096,-0.0068,-0.0120,0.0042];
c = [ 0.1803, 0.0872,-0.0153, 0.0145,0 ];
x = (h*0.001)/50.0-1;
sa = a(1) + a(2)*x + a(3)*x*x;
%rho = sum(arrayfun(@(i) b(i)*cos(i*pi*x) + c(i)*sin(i*pi*x), 1:5));
rho = 0;
for i=1:5
rho = rho + b(i)*cos(i*pi*x) + c(i)*sin(i*pi*x);
end
density = exp(sa + rho);
end