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. ------------------------------------------------------------------------- This SF.net email is sponsored by DB2 Express Download DB2 Express C - the FREE version of DB2 express and take control of your XML. No limits. Just data. Click to get it now. http://sourceforge.net/powerbar/db2/ _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
