Extended Kalman Filter Jacobians

  • Context: Engineering 
  • Thread starter Thread starter Master1022
  • Start date Start date
  • Tags Tags
    Filter Kalman filter
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 2K views
Master1022
Messages
590
Reaction score
116
Homework Statement
Calculate the Jacobean matrices for the extended Kalman filter
Relevant Equations
Partial derivatives
Hi,

I have a question about calculating the Jacobian matrices for the Extended Kalman filter.

Question: If we have a system of the form:
[tex]\begin{align*} x_{k+1} =f_k (x_k , u_k) + w_k \\<br /> y_k = h_k (x_k , u_k ) +v_k \end{align*}[/tex]
where the state [itex]x_k[/itex] comprises of the three variables [itex]p_1[/itex], [itex]p_2[/itex], and [itex]p_3[/itex]. The input [itex]u_k[/itex] comprises of the variables [itex]q_1[/itex] and [itex]q_2[/itex]. Form an extended Kalman filter for this system.

Method:
From what I understand, I need to linearize the system to form the A and C matrices for the state space model. I am, however, confused on how to form these Jacobian matrices. The source I am learning from presents the following formulae:
[tex]A_k = \frac{\partial f_k}{\partial x_k} |_{\hat x_{k|k} , u_k}[/tex]
and
[tex]C_k = \frac{\partial h_k}{\partial x_k} |_{\hat x_{k|k} , u_k}[/tex]

I don't understand how I ought to interpret these formulae. Should my matrices have the following form, or have I misunderstood something?
[tex]A_k =<br /> \begin{pmatrix}<br /> \frac{\partial f_1}{\partial p_1} & \frac{\partial f_1}{\partial p_2} & \frac{\partial f_1}{\partial p_3} \\<br /> \frac{\partial f_2}{\partial p_1} & \frac{\partial f_2}{\partial p_2} & \frac{\partial f_2}{\partial p_3} \\<br /> \frac{\partial f_3}{\partial p_1} & \frac{\partial f_3}{\partial p_2} & \frac{\partial f_3}{\partial p_3}<br /> \end{pmatrix}[/tex]
where [itex]f_i[/itex] refers to the function in the [itex]i^{th}[/itex] row of the vector [itex]f[/itex] (which has three rows) and
[tex]C_k =<br /> \begin{pmatrix}<br /> \frac{\partial h_1}{\partial p_1} & \frac{\partial h_1}{\partial p_2} & \frac{\partial h_1}{\partial p_3} \\<br /> \frac{\partial h_2}{\partial p_1} & \frac{\partial h_2}{\partial p_2} & \frac{\partial h_2}{\partial p_3} \\<br /> \end{pmatrix}[/tex]
where [itex]h_i[/itex] refers to the function in the [itex]i^{th}[/itex] row of the vector [itex]h[/itex] (which has two rows)

Thank you in advance for any help and guidance.
 
Last edited:
Physics news on Phys.org
spoto2 said:
Previously, we used a Kalman Filter which could only model linear ... It can be very difficult to implement an EKF correctly: you will most likely get the Jacobian.
Thanks for your reply. Do my expressions above seem correct though?