- #1
Master1022
- 611
- 117
- 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 \\
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 =
\begin{pmatrix}
\frac{\partial f_1}{\partial p_1} & \frac{\partial f_1}{\partial p_2} & \frac{\partial f_1}{\partial p_3} \\
\frac{\partial f_2}{\partial p_1} & \frac{\partial f_2}{\partial p_2} & \frac{\partial f_2}{\partial p_3} \\
\frac{\partial f_3}{\partial p_1} & \frac{\partial f_3}{\partial p_2} & \frac{\partial f_3}{\partial p_3}
\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 =
\begin{pmatrix}
\frac{\partial h_1}{\partial p_1} & \frac{\partial h_1}{\partial p_2} & \frac{\partial h_1}{\partial p_3} \\
\frac{\partial h_2}{\partial p_1} & \frac{\partial h_2}{\partial p_2} & \frac{\partial h_2}{\partial p_3} \\
\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 \\
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 =
\begin{pmatrix}
\frac{\partial f_1}{\partial p_1} & \frac{\partial f_1}{\partial p_2} & \frac{\partial f_1}{\partial p_3} \\
\frac{\partial f_2}{\partial p_1} & \frac{\partial f_2}{\partial p_2} & \frac{\partial f_2}{\partial p_3} \\
\frac{\partial f_3}{\partial p_1} & \frac{\partial f_3}{\partial p_2} & \frac{\partial f_3}{\partial p_3}
\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 =
\begin{pmatrix}
\frac{\partial h_1}{\partial p_1} & \frac{\partial h_1}{\partial p_2} & \frac{\partial h_1}{\partial p_3} \\
\frac{\partial h_2}{\partial p_1} & \frac{\partial h_2}{\partial p_2} & \frac{\partial h_2}{\partial p_3} \\
\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: