Как я могу получить dy/dt от odeint от SciPy

Question

Я хочу получить dy/dt от odeint от SciPy, а также от самого себя.

Я понимаю, что odeint получает решение от dy/dt=func(y0,t). Однако мне также нужно получить dy/dt.

Я имею дело со структурной динамической системы, в которой:

y : displacement 
dy/dt : velocity 
d(dy/dt)/dt : acceleration (force) 

И я поставил систему в две системы:

dy/dt = x 
dx/dt = func(y0,x0,t) 

Затем odeint(func, [y0's, x0's]) работает отлично. Однако конечная цель - получить силы в системе.

Большое спасибо заранее ^^

Answer 1

Просто чтобы быть уверенным, я понимаю вашу проблему правильно:

(d2x)/(dt2) = F(x,t) // your second order system that you wish to solve, to simplify the problem I assume that dim(x) = 1 

/* your system transformed to first order system of 2 equations 
* x[0] is your dx/dt (velocity), x[1] is your x (displacement) 
*/ 

dxdt[1] = x[0]; 
dxdt[0] = F(x[1],t); 

Результаты вы получите являются:

• х [0] - скорость
• x [1] - смещение
• t - время

Для вычисления ускорения вы просто должны сделать:

acc = F(x[1],t); 

Я надеюсь, что это то, что вы просили. Удачи.