Re bonjour,

J'ai un gros problème sur la génération de vecteurs dans mon programme.
Au tout début de mon programme, j'ai 2 vecteurs X et Y qui possèdent 27 valeurs chacun.

Après avoir déterminer mon time2 (discussion précédente) de dimension k, je cherche à déterminer un vecteur Y1 à partir d'autres paramètre. Y et Y1 sont de même unités.

Le problème survient lorsque Y1(i)>Ymax (valeur = 12500). Dans ce cas, au terme Y1(i+1), j'obtiens un NaN.
Dans ma boucle for i=1:length(time2), j'ai inclus un "if" pour tenter de supprimer les NaN
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
for i=1:length(time2) 
 Y1(i)=round(psi(i)/R/K/(h2(i)) %psi(i) est calculé avant Y1(i);h2(i), valeur de h2 à chaque pas;R et K sont des constantes
if  (Y1(i)>12500) && (h2(i)~=6)
     h2(i)=h2(i)+1;
     Y1(i)=12490;
  end
Y1(i)=round(psi(i)/R/K/(h2(i));
...
end
Y1 dépend entre autre de h2. Si h2 augmente, Y1 diminue (h2 = 4,5 ou 6).

Ce que j'aimerais faire, c'est que si Y1(i)>12500, j'incrémente de 1 h2 et je recalcule Y1(i) pour ne plus avoir de NaN, car j'obtiens également des NaN pour toutes les variables suivantes.

Merci