search

mattja / sdeint

Numerical integration of Ito or Stratonovich SDEs
GNU General Public License v3.0
154 stars 25 forks source link

A question: Mechanical system with white noise external forces #17

Open Peter230655 opened 2 years ago

Peter230655 commented 2 years ago

I have a mechanical system described as d/dt (q) = f(q, t, F), where q …generalized coordinates, F…external forces to be modeled as white noise (Where I used sympy.physics.mechanics to generate these equations of motion)

What I did in order to apply SDEINT, I linearized about the forces F: B = f(q, t, F).diff(F).subs({F: 0.}) F1 = f(q, t, F).subs({F: 0})

Then I calculated as Result = itoint(F1, B, x0, tspan)

Is this the right approach, or what would be the correct one?

Thanks for any help!

Peter230655 commented 1 year ago

Just a reminder, is my idea correct? Thanks a lot for your help!

Peter230655 commented 1 year ago

Has sdeint been 'abandoned'? Would be a shame, it is so nice and useful!