How to compensate drift in solving numerical differential equation

2 次查看（过去 30 天）

Vincent Gigot 2012-7-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43112-how-to-compensate-drift-in-solving-numerical-differential-equation

I work with experimental data from a force platform for human experiment. This platform records forces in the 3 dimensions, Fx, Fy, Fz. Sampling frequency is 50 Hz or 500 Hz.

Thanks to theses forces I would like to calculate velocities of the center of mass of the body and then the mechanical work of theses forces. To do this I must solve two differentials equations successively.

1st to get velocities V : dV/dt=F/M (M is the body mass)

2nd to get work W : dW/dt=F*V

I used ode113 or a personal function using Adams numerical method.

What's the best way to solve such equations, with F being vectors of experimental data and not mathematical function ?

Drift appears in the solutions, what kind of method should be used to avoid it ?

Thanks