0

For a standard Kalman filter

\begin{align} x_{t+1}&= Ax_t + w_t\\ y_t &= Cx_t + v_t, \end{align}

I am interested in designing the matrix $C$ such that the error covariance is minimized and $\operatorname{tr}(CC^T)\le p$. Can it be formulated as a semidefinite program (SDP)?


My attempt

The optimization can be written using the Riccati equation as:

$$ \begin{array}{ll} \underset {P \succeq0,C, K} {\text{minimize}} & \operatorname{tr} (P) \\ \text{subject to} & (A - K C) P (A - K C)^T + K R K^T + Q = P, \operatorname{tr}(CC^T)\le p \end{array} $$

Using my previous question in Kalman filter solution via LMI: summary of methods, we can write:

\begin{align} \begin{array}{ll} \underset {X,Y,Z,C} {\text{minimize}} & \operatorname{tr} (Z) \\ \text{such that} & \begin{bmatrix} X & XA - YC & Y & X\\ (XA - YC)^T & X & 0& 0\\ Y^T & 0 & R^{-1}&0\\ X& 0 & 0 & Q^{-1} \end{bmatrix} \succeq 0, \\ & \begin{bmatrix} Z & I\\ I & X \end{bmatrix} \succeq 0, \\ & \operatorname{tr}(CC^T)\le p \end{array} \end{align}

but there are still bilinear terms.


Related

Morad
  • 690
  • @RodrigodeAzevedo thanks for your edit. This question is also related to the one you solved a week ago :) – Morad Mar 11 '25 at 11:06
  • 1
    @RodrigodeAzevedo The question before that on the completion of the square for simple polynomials to find roots :) they are all around the same challenges but I try to simplify it as much as I can at each step – Morad Mar 11 '25 at 11:16
  • Can't one make $\text{tr}(C,C^\top)$ arbitrarily small using a state space model coordinate transformation? – Kwin van der Veen Mar 11 '25 at 11:20
  • @KwinvanderVeen. I am familiar with coordinate transf. but not sure what you mean by making its norm small (A is not necessarily stable). Can you please elaborate? – Morad Mar 11 '25 at 11:23
  • 1
    Are the dimensions of $C$ fixed? – KBS Mar 11 '25 at 16:04
  • I think we can assume it is square w.l.o.g. since it subsumes zero entries in the measurements vector – Morad Mar 11 '25 at 16:20
  • Normally, when solving physical problems, a measurement matrix is dictated by the structure of the system that you are building an observer for. Making any substantive change in that matrix involves making substantive, physical, changes to the system itself. This may involve some rather vigorous conversations with mechanical engineers, and whoever is holding the purse strings when you go to buy sensors. – TimWescott Mar 11 '25 at 20:07
  • @TimWescott. I understand your practical point of view, but it is not directly related to my mathematical question. Practically, you can think of a system where all sensors exist but their sampling is costly in some sense. – Morad Mar 12 '25 at 08:47

0 Answers0