Extended Kalman Filter Jacobians

  • Context: Engineering 
  • Thread starter Thread starter Master1022
  • Start date Start date
  • Tags Tags
    Filter Kalman filter
Click For Summary
SUMMARY

The discussion focuses on calculating Jacobian matrices for the Extended Kalman Filter (EKF) applied to a nonlinear system defined by state transition and measurement equations. The user seeks clarification on forming the A and C matrices, which represent the Jacobians of the system's state transition function and measurement function, respectively. The correct formulation of these matrices involves partial derivatives of the functions with respect to the state variables. The user’s proposed matrix structures for A_k and C_k are accurate, reflecting the necessary linearization for EKF implementation.

PREREQUISITES
  • Understanding of Extended Kalman Filter (EKF) principles
  • Knowledge of Jacobian matrix calculations
  • Familiarity with state-space representation of dynamic systems
  • Basic proficiency in calculus, particularly partial derivatives
NEXT STEPS
  • Study the derivation of Jacobian matrices in nonlinear systems
  • Learn about state-space modeling techniques
  • Explore numerical methods for implementing EKF
  • Investigate common pitfalls in EKF implementation and how to avoid them
USEFUL FOR

Engineers, data scientists, and researchers working on robotics, navigation systems, or any applications involving nonlinear state estimation using Extended Kalman Filters.

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:
\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*}
where the state x_k comprises of the three variables p_1, p_2, and p_3. The input u_k comprises of the variables q_1 and q_2. 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:
A_k = \frac{\partial f_k}{\partial x_k} |_{\hat x_{k|k} , u_k}
and
C_k = \frac{\partial h_k}{\partial x_k} |_{\hat x_{k|k} , u_k}

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

Similar threads

  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 7 ·
Replies
7
Views
22K
  • · Replies 17 ·
Replies
17
Views
10K
Replies
2
Views
3K
  • · Replies 175 ·
6
Replies
175
Views
27K
  • · Replies 71 ·
3
Replies
71
Views
13K