I'm performing system identification of the lateral closed-loop dynamics of a quadrotor. My model receives a setpoint and should return position and acceleration. I've proposed a second order model of the type, see equations, and solved it with MATLAB's greyest():
$\begin{bmatrix} v_y \\ a_y \end{bmatrix} = \begin{bmatrix}0&1 \\-\omega^2&-2\xi\omega \end{bmatrix} \begin{bmatrix}y \\ v_y \end{bmatrix} + \begin{bmatrix} 0 \\ \omega^2 \end{bmatrix} y_0$
$\begin{bmatrix} y \\ a_y \end{bmatrix} = \begin{bmatrix}1&0 \\-\omega^2&-2\xi\omega \end{bmatrix} \begin{bmatrix}y \\ v_y \end{bmatrix} + \begin{bmatrix} 0 \\ \omega^2 \end{bmatrix} y_0$
The model performs very well for the position, 97% on a cross-validation dataset, but the performance is very low for acceleration 43%.
Note the acceleration dataset comes out from an accelerometer so it's noisy, I've already tried low pass filtering the noise to improve performance with no positive results.
Any ideas on how to improve this?
Thanks!
What is $\xi$? Maybe the following (related to dead-reckoning) is way different than what you are actually asking? but...
References