Master1022
- 590
- 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.
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: