My plot doesnt exist
显示 更早的评论
my function:
function result=climbjet(H)
T0=288.15; %Standard atmospheric temperature at MSL
p0=101325; %Standard atmospheric pressure at MSL
ro0=1.225; %Standard atmospheric density at MSL
a0=340.294; %Speed of sound
K=1.4; %Adiabatic index of air
R=287.05287; %Real gas constant for air
g0=9.80665; %Gravitational acceleration
beta=-0.0065; %ISA temperature gradient with altitude below the tropopause
vstall=142; %from aircraft.OPF
vcl1=172; %335 knots
vcl2=172; %335 knots
machcl=0.86;
cvmin=1.3;
vdcl1=5;
vdcl2=10;
vdcl3=30;
vdcl4=60;
vdcl5=80;
%Mach transition altitude=29000 ft
if (H<=1499)
result=(cvmin*vstall)+vdcl1;
elseif (H<2999)
result=(cvmin*vstall)+vdcl2;
elseif (H<3999)
result=(cvmin*vstall)+vdcl3;
elseif (H<4999)
result=(cvmin*vstall)+vdcl4;
elseif (H<5999)
result=(cvmin*vstall)+vdcl5;
elseif (H<9999)
result=min([vcl1,129]);
elseif (H<28999)
result=vcl2;
elseif (H<=35000)
result=machcl*a0;
else
disp('error');
end
end
I am trying to plot with this code :
H=[0:100:35000];
plot(H,climbjet(H))
But plot does not appear. Please help
回答（1 个）
result = zeros(size(H));
for i=1:numel(H)
if H(i)<=1499
result(i) = (cvmin*vstall)+vdcl1;
elseif H(i)>1499 && H(i)<2999
result(i)=(cvmin*vstall)+vdcl2;
...
end % if statements
end % for loop
类别
在 帮助中心 和 File Exchange 中查找有关 Geodesy and Mapping 的更多信息
2021-12-10
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
