LMIs in Control/pages/H2 Optimal Filter

Optimal filtering is a means of adaptive extraction of a weak desired signal in the presence of noise and interfering signals. Optimal filters normally are free from stability problems. There are simple operational checks on an optimal filter when it is being used that indicate whether it is operating correctly. Optimal filters are probably easier to make adaptive to parameter changes than suboptimal filters.The goal of optimal filtering is to design a filter that acts on the output $z$ of the generalized plant and optimizes the transfer matrix from w to the filtered output.

The System:

Consider the continuous-time generalized LTI plant with minimal states-space realization

$\begin{array}{cc}\dot{x}& =Ax+B_{1}w\\ z& =C_{1}x+D_{11}w,\\ y& =C_{2}x+D_{21}w,\end{array}$

where it is assumed that $A$ is Hurwitz.

The Data

The matrices needed as inputs are $A,B_1,C_2,D_{21}$.

The Optimization Problem:

An $H_2$-optimal filter is designed to minimize the $H_2$ norm of $\tilde{P}(s)$ in following equation.

$\begin{array}{c}\tilde{P}(s)={\tilde{C}}_1(sI-\tilde{A}{)}^{-}1{\tilde{B}}_1+{\tilde{D}}_{11},\\ \text{where}\\ \tilde{A}=\left[\begin{array}{ccc}A& & 0\\ B_f{C}_2& & A_f\end{array}\right]& <0\\ {\tilde{B}}_1=\left[\begin{array}{c}{B}_1\\ B_f{D}_{21}\end{array}\right]& <0\\ {\tilde{C}}_1=\left[\begin{array}{c}{C}_1-D_f{C}_2-C_f\end{array}\right]& <0\\ {\tilde{D}}_{11}=D_{11}-D_f{D}_{21}\\ \end{array}$

To ensure that $\tilde{P}(s)$ has a finite $H_2$ norm, it is required that $D_f=D_{11}$, which results in ${\tilde{D}}_{11}=D_{11}-D_f=0$

The LMI: $H_2$- Optimal filter

Solve for $A_n\in {ℝ}^{n_x\times n_x},B_n\in {ℝ}^{n_x\times n_y},C_f\in {ℝ}^{n_x\times n_x}$, $X,Y\in {𝕊}^{n_x},Z\in {𝕊}^{n_z}$ and $\nu \in {ℝ}_{>0}$ that minimize $\zeta (\nu )=\nu$ subject to $X>0,Y>0,Z>0$.

$\begin{array}{cc}\left[\begin{array}{ccccc}YA+A^{T}Y+B_n{C}_2& & A_n+C_2^T{B}_n^T+A^{T}X& & Y{B}_1+B_n{D}_{21}\\ \star & & A_n+A_n^T& & X{B}_1+B_n{D}_{21}\\ \star & & \star & & -1\end{array}\right]& <0\\ \left[\begin{array}{ccccc}-Z& & C_1-D_f{C}_2& & -C_f\\ \star & & -Y& & -X\\ \star & & \star & & -X\end{array}\right]& <0\\ Y-X>0\\ trZ</nu\end{array}$

Conclusion:

The filter is recovered by $A_f=X^{-1}{A}_n$ and $B_f=X^{-1}{B}_n$.

Implementation

MATLAB code of $H_2$ Optimal filter

External links

• [1]- Optimal Filtering
• LMI Properties and Applications in Systems, Stability, and Control Theory - A List of LMIs by Ryan Caverly and James Forbes.
• [http://users.cecs.anu.edu.au/~john/papers/BOOK/B02.PDF}- Optimal Filtering by Brian D.O. Anderson