Discussion :

[fracx]

petit souci : je voudrais definir une fonction par morceaux comme ceci

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
function y = hum(t)
d=365*24*3600;
if t<9.833,
    y = 7.9516+52.0484*exp(-1.0433*10^(-8)*t*d);
else 
    y = 6;
end

probleme, quand j'applique ma fonction à un vecteur, elle ne renvoie que l'image du premier élément du vecteur. Impossible par exemple d'afficher le graphe de la fonction.
une solution?

[rostomus]

Bonjour,

le probleme est là:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

if t<9.833

cette condition n'est vraie que si tout les element de "t" sont inferieur à 9.833.
je te propose ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
function y = hum(t)
d=365*24*3600;
y = 7.9516+52.0484*exp(-1.0433*10^(-8)*t*d);
y(t>=9.833)=6;

[fracx]

merci, c'est impeccable !
quelle rapidité !

[fracx]

mince, une idée pourquoi ceci ne marche pas ??

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
function y = morc(t)
d=365*24*3600;
y = -11.8635-0.05321*exp(1.1468*10^(-8)*t*d);
y(t<=9.833) = 7.9516+52.0484*exp(-1.0433*10^(-8)*t*d);
y(16.6419<=t)=-34.03225-44.03225*exp(-5.2163*10^(-9)*t*d);

[rostomus]

parce que le nombre d'elements dans "y(t<=9.833)" est different de celui dans "t".
tu dois faire ça:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
function y = morc(t)
d=365*24*3600;
y = -11.8635-0.05321*exp(1.1468*10^(-8)*t*d);
y(t<=9.833) = 7.9516+52.0484*exp(-1.0433*10^(-8)*t(t<=9.833)*d);
y(16.6419<=t) =-34.03225-44.03225*exp(-5.2163*10^(-9)*t(16.6419<=t)*d);

[fracx]

Genial, grand merci !