Équation de Liapounov

En théorie du contrôle, l'équation discrète de Liapounov (également connue sous le nom d'équation de Stein ) est une équation de la forme

${\displaystyle AXA^{H}-X+Q=0}$,

où ${\displaystyle Q}$ est une matrice hermitienne et ${\displaystyle A^{H}}$ est la matrice adjointe de ${\displaystyle A}$.

L'équation continue de Liapounov est de la forme

${\displaystyle AX+XA^{H}+Q=0}$.

L'équation de Liapounov apparaît dans de nombreuses branches de la théorie du contrôle, telles que la stabilité de Liapounov et la commande optimale. Cette équation et des équations associées portent le nom du mathématicien russe Alexandre Liapounov^[1],[2].

Dans les énoncés suivants ${\displaystyle A,P,Q\in \mathbb {R} ^{n\times n}}$, et ${\displaystyle P}$ et ${\displaystyle Q}$ sont des matrices symétriques. La notation ${\displaystyle P>0}$ signifie que la matrice ${\displaystyle P}$ est définie positive.

Théorème (version temps continu) — Étant donné ${\displaystyle Q>0}$, il existe un unique ${\displaystyle P>0}$ satisfaisant

${\displaystyle A^{T}P+PA+Q=0}$

si et seulement si le système linéaire ${\displaystyle {\dot {x}}=Ax}$ est globalement asymptotiquement stable.

La fonction quadratique ${\displaystyle V(x)=x^{T}Px}$ est une fonction de Liapounov qui peut être utilisée pour vérifier la stabilité.

Théorème (version en temps discret) — Étant donné ${\displaystyle Q>0}$, il existe un unique ${\displaystyle P>0}$ satisfaisant

${\displaystyle A^{T}PA-P+Q=0}$

si et seulement si le système linéaire ${\displaystyle x_{t+1}=Ax_{t}}$ est globalement asymptotiquement stable.

Comme ci-dessus, ${\displaystyle x^{T}Px}$ est une fonction de Liapounov.

L'équation de Liapounov est linéaire, et donc si ${\displaystyle X}$ est de taille ${\displaystyle n}$, il peut être calculé en temps ${\displaystyle {\mathcal {O}}(n^{3})}$ en utilisant les méthodes standard de factorisation matricielle.

Cependant, la structure spécifique de l'équation de Liapounov permet l'usage d'algorithmes beaucoup plus rapides. Dans le cas discret, la méthode de Schur de Kitagawa est souvent utilisée^[3]. Dans le cas de l'équation de Liapounov continue, l'algorithme de Bartels-Stewart peut être utilisé^[4].

On considère l'opérateur de vectorisation ${\displaystyle \operatorname {vec} (A)}$ qui empile les colonnes d'une matrice ${\displaystyle A}$, et on note ${\displaystyle A\otimes B}$ le produit de Kronecker de ${\displaystyle A}$ et de ${\displaystyle B}$. Les équations de Liapounov en temps continu et en temps discret peuvent être exprimées comme des solutions d'une équation matricielle. De plus, si la matrice ${\displaystyle A}$ est stable, la solution peut également être exprimée sous la forme d'une intégrale (cas du temps continu) ou d'une somme infinie (cas du temps discret).

En utilisant le fait que ${\displaystyle \operatorname {vec} (ABC)=(C^{T}\otimes A)\operatorname {vec} (B)}$, on a

${\displaystyle (I_{n^{2}}-{\bar {A}}\otimes A)\operatorname {vec} (X)=\operatorname {vec} (Q)}$,

où ${\displaystyle I_{n^{2}}}$ est la matrice identité de taille ${\displaystyle n^{2}}$ et ${\displaystyle {\bar {A}}}$ est la matrice adjointe de ${\displaystyle A}$^[5]. On peut alors résoudre ${\displaystyle \operatorname {vec} (X)}$ par inversion ou en résolvant les équations linéaires. Pour obtenir ${\displaystyle X}$, il suffit de recomposer la matrice depuis ${\displaystyle \operatorname {vec} (X)}$.

De plus, si ${\displaystyle A}$ est stable, la solution ${\displaystyle X}$ peut aussi s'écrire sous la forme

${\displaystyle X=\sum _{k=0}^{\infty }A^{k}Q(A^{H})^{k}}$ .

À titre d'exemple, considérons le cas unidimensionnel, où la formule dit simplement que la solution de ${\displaystyle (1-a^{2})x=q}$ est la fraction

${\displaystyle x={\frac {q}{1-a^{2}}}=\sum _{k=0}^{\infty }qa^{2k}}$ .

En utilisant à nouveau le produit de Kronecker et l'opérateur de vectorisation, on a l'équation matricielle

${\displaystyle (I_{n}\otimes A+{\bar {A}}\otimes I_{n})\operatorname {vec} X=-\operatorname {vec} Q,}$

où ${\displaystyle {\bar {A}}}$ est la matrice adjointe de ${\displaystyle A}$.

Comme dans le cas discret, si ${\displaystyle A}$ est stable, la solution ${\displaystyle X}$ peut aussi s'écrire comme

${\displaystyle X=\int _{0}^{\infty }{e}^{A\tau }Q\mathrm {e} ^{A^{H}\tau }d\tau }$ .

À titre d'exemple, considérons le cas unidimensionnel ; l'expression dit simplement que la solution de ${\displaystyle 2ax=-q}$ est

${\displaystyle x={\frac {-q}{2a}}=\int _{0}^{\infty }q{e}^{2a\tau }d\tau }$ .

On considère la dynamique linéaire en temps continu :

${\displaystyle {\dot {\mathbf {x} }}=\mathbf {A} \mathbf {x} }$ .

On la discrétise en :

${\displaystyle {\dot {\mathbf {x} }}\approx {\frac {\mathbf {x} _{t+1}-\mathbf {x} _{t}}{\delta }}}$,

où ${\displaystyle \delta >0}$ indique un petit déplacement dans le temps. En explicitant l'équation et en mélangeant les termes, on obtient une équation en temps discret pour ${\displaystyle \mathbf {x} _{t+1}}$ :

${\displaystyle \mathbf {x} _{t+1}=\mathbf {x} _{t}+\delta \mathbf {A} \mathbf {x} _{t}=(\mathbf {I} +\delta \mathbf {A} )\mathbf {x} _{t}=\mathbf {B} \mathbf {x} _{t}}$

où ${\displaystyle \mathbf {B} \equiv \mathbf {I} +\delta \mathbf {A} }$ . Maintenant, on peut utiliser l'équation de Liapounov en temps discret pour ${\displaystyle \mathbf {B} }$ :

${\displaystyle \mathbf {B} ^{T}\mathbf {M} \mathbf {B} -\mathbf {M} =-\delta \mathbf {Q} }$.

Avec la définition de ${\displaystyle \mathbf {B} }$, on a :

${\displaystyle (\mathbf {I} +\delta \mathbf {A} )^{T}\mathbf {M} (\mathbf {I} +\delta \mathbf {A} )-\mathbf {M} =-\delta \mathbf {Q} }$.

En développant cette expression, on obtient :

${\displaystyle (\mathbf {M} +\delta \mathbf {A} ^{T}\mathbf {M} )(\mathbf {I} +\delta \mathbf {A} )-\mathbf {M} =\delta (\mathbf {A} ^{T}\mathbf {M} +\mathbf {M} \mathbf {A} )+\delta ^{2}\mathbf {A} ^{T}\mathbf {M} \mathbf {A} =-\delta \mathbf {Q} }$.

Comme ${\displaystyle \delta }$ est petit, si ${\displaystyle \delta }$ tend zéro, on s'approche de plus en plus d'une dynamique continue, et c'est ce qu'on obtient à la limite. On peut également récupérer les équations de Lyapounov en temps continu à la limite. Pour cela, on divise par ${\displaystyle \delta }$ des deux côtés, puis quand ${\displaystyle \delta \to 0}$, on trouve que :

${\displaystyle \mathbf {A} ^{T}\mathbf {M} +\mathbf {M} \mathbf {A} =-\mathbf {Q} }$

qui est l'équation de Liapounov en temps continu, comme indiqué.

• Équation de Sylvester
• Équation algébrique de Riccati
• Filtre de Kalman

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Lyapunov equation » (voir la liste des auteurs).

• Portail des mathématiques
• Portail de l'informatique théorique