LMIs in Control/pages/Discrete Time QSR

The Concept

In systems theory the concept of dissipativity was first introduced by Willems which describes dynamical systems by input-output properties. Considering a dynamical system described by its state $x(t)$, its input $u(t)$ and its output $y(t)$, the input-output correlation is given a supply rate $w(u(t),y(t))$. A system is said to be dissipative with respect to a supply rate if there exists a continuously differentiable storage function $V(x(t))$ such that $V(0)=0$, $V(x(t))\ge 0$ and

$\dot{V}(x(t))\le w(u(t),y(t))$

As a special case of dissipativity, a system is said to be passive if the above dissipativity inequality holds with respect to the passivity supply rate $w(u(t),y(t))=u(t{)}^{T}y(t)$.

The physical interpretation is that $V(x)$ is the energy stored in the system, whereas $w(u(t),y(t))$ is the energy that is supplied to the system.

This notion has a strong connection with Lyapunov stability, where the storage functions may play, under certain conditions of controllability and observability of the dynamical system, the role of Lyapunov functions.

Roughly speaking, dissipativity theory is useful for the design of feedback control laws for linear and nonlinear systems. Dissipative systems theory has been discussed by Vasile M. Popov, Jan Camiel Willems, D.J. Hill, and P. Moylan. In the case of linear invariant systems, this is known as positive real transfer functions, and a fundamental tool is the so-called Kalman–Yakubovich–Popov lemma which relates the state space and the frequency domain properties of positive real systems.Dissipative systems are still an active field of research in systems and control, due to their important applications.

The System

Consider a discrete-time LTI system, $𝒢:{𝓁}_{2e}\to {𝓁}_{2e}$, with minimal state-space relization $({𝒜}_d,{ℬ}_d,{𝒞}_d,{𝒟}_d)$, where ${𝒜}_d\in {ℛ}^{n\times n},{ℬ}_d\in {ℛ}^{n\times m},{𝒞}_d\in {ℛ}^{p\times n},$ and ${𝒟}_d\in {ℛ}^{p\times m}$.

$x(k+1)={𝒜}_{d}x(k)+{ℬ}_{d}u(k)$

$y(k)={𝒞}_{d}x(k)+{𝒟}_{d}u(k),k=0,1...$

The Data

The matrices ${𝒜}_d,{ℬ}_d,{𝒞}_d$ and ${𝒟}_d$

The Optimization Problem

The system $𝒢$ is QSR disipative if

$\sum _{i=0}^K({𝓎}_i^{T}Q{𝓎}_i+2{𝓎}_i^{T}S{𝓊}_i+{𝓊}_i^{T}R{𝓊}_i)dt\ge 0,\mathrm{\forall }u\in {𝓁}_{2e},\mathrm{\forall }k\in 𝒵\ge 0$

where ${𝓊}_k$ is the input to $𝒢,{𝓎}_k$ is the output of $𝒢,Q\in {𝒮}^p,S\in {ℛ}^{p\times m},$ and $ℛ\in {𝒮}^m$.

LMI : Discrete-Time KYP Lemma for QSR Dissipative Systems

The system $𝒢$ is also QSR dissipative if and only if there exists $P\in {𝒮}^n,$ where $P>0,$ such that

$\left[\begin{array}{cc}{A}_d^{T}P{A}_d-P-C_d^{T}Q{C}_d& A_d^{T}P{B}_d-C_d^{T}S-C_d^{T}Q{D}_d\\ (A_d^{T}P{B}_d-C_d^{T}S-C_d^{T}Q{D}_d{)}^T& B_d^{T}P{B}_d-D_d^{T}Q{D}_d-(D_d^{T}S+S^T{D}_d)-R\end{array}\right]\le 0.$

Conclusion:

If there exist a positive definite $P$ for the the selected Q,S and R matrices then the system $𝒢$ is QSR dissipative.

Implementation

Code for implementation of this LMI using MATLAB. https://github.com/VJanand25/LMI

Related LMIs

KYP Lemma
KYP Lemma for continous Time QSR Dissipative system

References

1. J. C. Willems, “Dissipative dynamical systems - part I: General theory,” Archive Rational Mechanics and Analysis, vol. 45, no. 5, pp. 321–351, 1972.
2. D. J. Hill and P. J. Moylan, “The stability of nonlinear dissipative systems,” IEEE Transac- tions on Automatic Control, vol. 21, no. 5, pp. 708–711, 1976.
3. LMI Properties and Applications in Systems, Stability, and Control Theory, by Ryan James Caverly1 and James Richard Forbes2