Discussion :

[Ionosphere]

Bonjour, je sollicite votre aide concernant un problème sur lequel je bute depuis quelques jours.

Ma fonction de départ est une double exponentielle de la forme suivante : y = a * ( exp (-t/tau1) - exp (-t/tau2) )

Je souhaite retrouver les paramètres de cette fonction qui sont a , tau1 et tau2.

Pour a, il me suffit de chercher le max de la courbe, donc je pose a = max (y) .

Pour tau1, si je sépare les exponentielle, j'ai y = a* exp(-t/tau1).
y(t=tau1) = a* exp(-tau1/tau1) = a* exp(-1) = 0.37*a où a = max(y).
Idem pour tau2.

On a donc une ordonnée qui vaut 0.37*max(y) et j'aimerai retrouver les points d'abscisse correspondant à cette ordonnée pour tau1 et tau2.

Si je pose dans Matlab :

Code :

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

1
2
3
x=[0:1:250] ;
y = 0.5* (exp(-x/100)-exp(-x/50));
plot(x,y)

et que je regarde à l'ordonnée 0.37*max(y) je retombe bien sur les abscisses x=50 et x=100, qui correspondent à mes tau2 et tau1.

Pourtant, lorsque je veux faire :

Code :

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

1
2
3
4
5
x=[0:1:250] ;
y = 0.5* (exp(-x/100)-exp(-x/50));
plot(x,y)
a=max(y);
tau1= find (y==0.37*a)

J'ai comme résultat : tau1 = empty matrix 1-by-0


Je coince complétement...

Merci beaucoup pour votre aide.

[Gooby]

Bonjour,


Je ne suis pas sûr que dans ton vecteur y existe la valeur 0.37*a. A ta place je diminuerai le pas de x afin d'améliorer la précision et je remplacerai cette ligne:

Code :

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

tau1= find (y==0.37*a)

par quelque chose du genre

Code :

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

tau1= find(y>0.36*a & y<0.38*a)

Possible que je me trompe, je suis en phase d'apprentissage!

Gooby

[Jerome Briot]

J'ai comme résultat : tau1 = empty matrix 1-by-0

Voir la FAQ : Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?

[Ionosphere]

Gooby, cette solution est bonne, je n'avais pas pensé au fait que la courbe que je voyais ne contenait pas de points pour toutes les valeurs, donc il faut passer par un encadrement !

Merci beaucoup !

Dut, merci pour le lien sur la FAQ, j'avais déjà lu la chose, mais j'ai eu du mal à l'adapter à ma situation.