In control theory, we may need to find out whether or not a system such as
\begin{array}{c} \boldsymbol{x |
is observable, where
\boldsymbol{A}
\boldsymbol{B}
\boldsymbol{C}
\boldsymbol{D}
n x n
n x p
q x n
q x p
One of the many ways one can achieve such goal is by the use of the Observability Gramian.
Linear Time Invariant (LTI) Systems are those systems in which the parameters
\boldsymbol{A}
\boldsymbol{B}
\boldsymbol{C}
\boldsymbol{D}
One can determine if the LTI system is or is not observable simply by looking at the pair
(\boldsymbol{A},\boldsymbol{C})
1. The pair
(\boldsymbol{A},\boldsymbol{C})
2. The
n x n
\boldsymbol{Wo
is nonsingular for any
t>0
3. The
nq x n
\left[\begin{array}{c} \boldsymbol{C}\\ \boldsymbol{CA}\\ \boldsymbol{CA}2\\ \vdots\\ \boldsymbol{CA}n-1\end{array}\right]
has rank n.
4. The
(n+q) x n
\left[\begin{array}{c} \boldsymbol{A}\boldsymbol{-λ}\boldsymbol{I}\\ \boldsymbol{C} \end{array}\right]
has full column rank at every eigenvalue
λ
\boldsymbol{A}
If, in addition, all eigenvalues of
\boldsymbol{A}
\boldsymbol{A}
\boldsymbol{AT
is positive definite, then the system is observable. The solution is called the Observability Gramian and can be expressed as
\boldsymbol{Wo
In the following section we are going to take a closer look at the Observability Gramian.
The Observability Gramian can be found as the solution of the Lyapunov equation given by
\boldsymbol{AT
In fact, we can see that if we take
\boldsymbol{Wo
as a solution, we are going to find that:
\begin{array}{ccccc} \boldsymbol{AT
Where we used the fact that
e\boldsymbol{At}=0
t=infty
\boldsymbol{A}
\boldsymbol{W}o
We can see that
\boldsymbol{CTC}
\boldsymbol{W}o
We can use again the fact that, if
\boldsymbol{A}
\boldsymbol{W}o
\boldsymbol{AT
and they are given by
\boldsymbol{W}o1
\boldsymbol{W}o2
\boldsymbol{AT
Multiplying by
\boldsymbol{AT | |
e |
t}
e\boldsymbol{At}
\boldsymbol{AT | |
e |
t}[\boldsymbol{AT
Integrating from
0
infty
\boldsymbol{AT | |
[e |
t}[(\boldsymbol{W}o1-\boldsymbol{W}o2)e\boldsymbol{A
infty | |
t}]| | |
t=0 |
=\boldsymbol{0}
using the fact that
e\boldsymbol{At} → 0
t → infty
\boldsymbol{0}-(\boldsymbol{W}o1-\boldsymbol{W}o2)=\boldsymbol{0}
In other words,
\boldsymbol{W}o
Also, we can see that
\boldsymbol{xTWo
infty | |
x}=\int | |
0 |
\boldsymbol{x}T
\boldsymbol{AT | |
e |
t}\boldsymbol{CTC}e\boldsymbol{A
infty | |
t}\boldsymbol{x}dt=\int | |
0 |
\left\Vert\boldsymbol{Ce\boldsymbol{At}\boldsymbol{x}}\right\Vert
2 | |
2 |
dt
is positive for any
\boldsymbol{x}
{\displaystyle{\boldsymbol{Ce{\boldsymbol
\boldsymbol{W}o
More properties of observable systems can be found in,[1] as well as the proof for the other equivalent statements of "The pair
(\boldsymbol{A},\boldsymbol{C})
For discrete time systems as
\begin{array}{c} \boldsymbol{x}[k+1]\boldsymbol{=Ax}[k]+\boldsymbol{Bu}[k]\\ \boldsymbol{y}[k]=\boldsymbol{Cx}[k]+\boldsymbol{Du}[k] \end{array}
One can check that there are equivalences for the statement "The pair
(\boldsymbol{A},\boldsymbol{C})
We are interested in the equivalence that claims that, if "The pair
(\boldsymbol{A},\boldsymbol{C})
\boldsymbol{A}
1
\boldsymbol{A}
\boldsymbol{AT
is positive definite and given by
\boldsymbol{W}do
infty | |
=\sum | |
m=0 |
(\boldsymbol{A}T)m\boldsymbol{C}T\boldsymbol{C}\boldsymbol{A}m
That is called the discrete Observability Gramian. We can easily see the correspondence between discrete time and the continuous time case, that is, if we can check that
\boldsymbol{W}dc
\boldsymbol{A}
1
(\boldsymbol{A},\boldsymbol{B})
Linear time variant (LTV) systems are those in the form:
\begin{array}{c} \boldsymbol{x |
That is, the matrices
\boldsymbol{A}
\boldsymbol{B}
\boldsymbol{C}
(\boldsymbol{A}(t),\boldsymbol{C}(t))
The system
(\boldsymbol{A}(t),\boldsymbol{C}(t))
t0
t1>t0
n x n
\boldsymbol{W}o(t0,t1
| |||||
)=\int | |||||
t0 |
\boldsymbol{\Phi}T(\tau,t0)\boldsymbol{C}T(\tau)\boldsymbol{C}(\tau)\boldsymbol{\Phi}(\tau,t0)d\tau
where
\boldsymbol{\Phi}(t,\tau)
\boldsymbol{x |
Again, we have a similar method to determine if a system is or not an observable system.
\boldsymbol{W}o(t0,t1)
We have that the Observability Gramian
\boldsymbol{W}o(t0,t1)
\boldsymbol{W}o(t0,t1)=\boldsymbol{W}o(t0,t)+\boldsymbol{\Phi}T(t,t0)\boldsymbol{W}o(t,t0)\boldsymbol{\Phi}(t,t0)
that can easily be seen by the definition of
\boldsymbol{W}o(t0,t1)
\boldsymbol{\Phi}(t0,t1)=\boldsymbol{\Phi}(t1,\tau)\boldsymbol{\Phi}(\tau,t0)
More about the Observability Gramian can be found in.[3]