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 :
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.
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)
Pourtant, lorsque je veux faire :
J'ai comme résultat : tau1 = empty matrix 1-by-0
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)
Je coince complétement...
Merci beaucoup pour votre aide.
Partager