[SciPy-user] odeint with digital data

Ryan Krauss ryanlists@gmail....
Mon Jun 4 12:04:38 CDT 2007


Thanks Rob.  Perhaps the time has finally come.  Yes, please help me
set this up in PyDSTool.  I will send a schematic description of the
system in a few minutes.

On 6/4/07, Rob Clewley <rhc28@cornell.edu> wrote:
> Hi Ryan,
>
> To my knowledge you cannot do this with odeint unless you change a
> constant value on the RHS (i.e. technically change your system) to
> reflect the changing input value after every time-step, but it is very
> easy to do in PyDSTool. There you can also force integration to be
> only at the discrete mesh points of your input signal (and would
> linearly interpolate in-between otherwise). I am certainly happy to
> help you set up your script to do this if you wish to try it in
> PyDSTool.
>
> HTH!
> Rob
>
> On 04/06/07, Ryan Krauss <ryanlists@gmail.com> wrote:
> > I have a vector of experimental data that I need to use as part of a
> > system of ode's.  I would like to solve this system using
> > integrate.odeint.  Can odeint be forced to solve only at the discrete
> > points in time where the experimentally measured signal is available?
> > Or do I need to set up some interpolation function to find the
> > measured signal at any time?  If that signal is the only thing that
> > explictly depends on time and I set up a digital look-up table that
> > returns the same constant value for the range from t to t+dt, would I
> > effectively force odeint to do what I want?  Am I making any sense?
> > Is there a better way?
> >
> > Thanks,
> >
> > Ryan
> > _______________________________________________
> > SciPy-user mailing list
> > SciPy-user@scipy.org
> > http://projects.scipy.org/mailman/listinfo/scipy-user
> >
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fbd.jpg
Type: image/jpeg
Size: 63348 bytes
Desc: not available
Url : http://projects.scipy.org/pipermail/scipy-user/attachments/20070604/a3ed8b42/attachment-0001.jpg 


More information about the SciPy-user mailing list