darkside wrote:
>
>
> Hello everyone:
> My question is not only about matplotlib, but I wonder if I can do what
> I want with some matplotlib subroutine, instead of the one for scipy.
> I'm new using scipy, so I'm sorry if any of my questions are silly.
>
> I'm trying to find the maxima, absolut and local, of a function, in order to
> fit an exponencial curve and get de exponecial argument.
>
> My function if the soluction of a couple equations system:
> 
> def derivs3(x,t,gamma,omega,dl):
> d1 = omega*x[2]  gamma *x[0]
> d2 = dl*x[2]  (gamma/2.)* x[1]
> d3 = omega *x[0]  dl*x[1]  (gamma/2.)* x[2] + (omega/2.)
> return d1,d2,d3
> def solucion(a,t,gamma, omega, dl):
> sol=odeint(derivs3,a,t,(gamma,omega,dl))
> return sol
> 
> The case I'm interesting in, the soluction have the form of a sin*exp, so I
> want to find the function that envolves it, a exponencial function.
> To do this, I can find the maximas, and fit them, so I use:
Why don't you fit the solution to sin*exp? You could estimate frequency/phase
via fft and use it as an initial guess.
Christian
ps: It's considered to be rude to post the same message three times within one
day. Btw. many people here read both scipy and matplotlib mailing lists.
