Boucle if

[turbulent]

Bonjour tout le monde;
Je cherche à avoir des valeurs de la variable S(h>=1)=1. Or, dans mon code Octave, même pour des valeurs de h>=1, j'obtiens des valeurs de S < 1. Je n'arrive pas à détecter l'erreur. Merci d'avance.

clc
clear all
alpha=1.15;
nn=2;
m=1-1/nn;
h=-5:.1:1;
S=(1+(alpha*abs(h).^nn)).^(-m);
  if (h >= 0)
     S=1;
  endif
plot(h,S)

-----

[pm42]

J'ai l'impression que ta condition (le if) est juste une instruction exécutée indépendamment du calcul en 1 fois de S.

[Paraboloide_Hyperbolique]

Bonjour,

La condition dans le "if" n'est pas correcte sur deux points:

1. h est est un vecteur et donc

h >= 0

aussi. Or un "if" demande un scalaire.
2. Il y a une contradiction entre la condition

h >= 0

du code et la demande mentionnant h >= 1...

Il faut soit utiliser une boucle sur les valeurs de h, soit (et c'est plus performant pour Octave/Matlab) tirer partit des capacités de vectorisation d'Octave.

[turbulent]

Même en utilisant une boucle sur les valeurs de h, je n'obtiens sur le graphe qu'un seul point (le dernier point de la liste h)

clear all
alpha=1.15;
nn=2;
m=1-1/nn;
%h=-5:.1:1;
deltah=.1;
h=-5:deltah:1;
N=length(h);
h(1)=-5;
for n=2:N;
   h(n)=h(n-1)+deltah; 
if h(n) < 0
S(n)=(1+(alpha*abs(h(n)).^nn)).^(-m);
else
S(n)=1;
endif
endfor
plot(h(n),S(n))

[A voir en vidéo sur Futura]

[turbulent]

C'est bon, je viens à trouver la bonne réponse en utilisant "plot(h(,S()". Merci.

[Ernum]

Salut,

C'est bon, je viens à trouver la bonne réponse en utilisant "plot(h(,S()". Merci.

il fallait lire ça je supposes:

plot(h(:),S(:))

[umfred]

je pense qu'on peut se passer de la parenthèse => plot(h,S) vu que ce sont des tableaux de même taille