# what to do when array indicies exceed for , FOR loop comparision for below , trying to plot density of atmosphere, how to do the comparison ?

2 vues (au cours des 30 derniers jours)
SL le 24 Sep 2021
Commenté : Walter Roberson le 25 Sep 2021
h(1) = 0;
T(1) = 273.15;
p(1) = 101.2;
for t = 1:328084
if h(t)<[11000]
T(t+1) = 15.04 - 0.00649.*h(t);
p(t+1) = 101.29 .*((T(t) + 273.15) /288.08).^5.256;
elseif [11000]> h(t) && h(t)<=[25000]
T(t) = -56.46;
p(t+1) = 22.65.*e^(1.73-0.000157.*h(t));
elseif h(t)>[25000]
T(t) = -131.21+0.00299.*h(t);
p(t+1) = 2.488 .* ((T(t) + 273.15)/216.6).^-11.388;
end
rho_h(t) = p(t) ./(0.2869 * (T(t) + 273.15))
end
figure(1)
plot(T,h)
grid on
legend on
figure(2)
plot(rho_h, h)
grid on
legend on
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponse acceptée

Walter Roberson le 24 Sep 2021
h(1) = 0;
That is the only place you assign to h
if h(t)<[11000]
but you need h to be 328084 elements by the end of the loop.
##### 8 commentairesAfficher 6 commentaires plus anciensMasquer 6 commentaires plus anciens
SL le 25 Sep 2021
Modifié(e) : SL le 25 Sep 2021
perfect thank you very much highly appreciate it
Walter Roberson le 25 Sep 2021
please recheck rho_h(1) which is not assigned a value.

Connectez-vous pour commenter.

### Catégories

En savoir plus sur MATLAB dans Help Center et File Exchange

R2021b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by