This paper presents a new interpretation of the Lyapunov and Riccati equations from the perspective of positive system theory. We show it is possible to construct positive systems related to these equations, and then certain conclusions -- such as the existence and uniqueness of solutions -- can be drawn from positive systems theory. Specifically, under standard observability assumptions, a strictly positive linear system can be constructed for Lyapunov equations, leading to exponential convergence in Hilbert metric to the Perron-Frobenius vector -- closely related to the solution of the Lyapunov equation. For algebraic Riccati equations, homogeneous strictly positive systems can be constructed, which exhibit more complex dynamical behaviors. While the existence and uniqueness of the solution can still be proven, only asymptotic convergence can be obtained.
- ID Articolo: 2511.17243
- Titolo: Lyapunov and Riccati Equations from a Positive System Perspective
- Autori: Dongjun Wu (Lund University), Yankai Lin (Wuhan Institute of Technology)
- Classificazione: math.OC (Ottimizzazione e Controllo)
- Data di Presentazione: Sottomesso ad arXiv il 24 novembre 2025
- Link Articolo: https://arxiv.org/abs/2511.17243
Questo articolo propone una nuova interpretazione delle equazioni di Lyapunov e Riccati dalla prospettiva della teoria dei sistemi positivi (positive system theory). La ricerca dimostra che è possibile costruire sistemi positivi associati a queste equazioni e derivare dalla teoria dei sistemi positivi conclusioni riguardanti l'esistenza e l'unicità delle soluzioni. Nello specifico, sotto ipotesi standard di osservabilità, è possibile costruire sistemi lineari strettamente positivi per l'equazione di Lyapunov, che convergono esponenzialmente nella metrica di Hilbert al vettore di Perron-Frobenius, strettamente correlato alla soluzione dell'equazione di Lyapunov. Per l'equazione algebrica di Riccati, è possibile costruire sistemi strettamente positivi omogenei, che sebbene esibiscano dinamiche più complesse, consentono comunque di provare l'esistenza e l'unicità della soluzione, con la convergenza che passa da esponenziale ad asintotica.
Le equazioni di Lyapunov e l'equazione algebrica di Riccati sono fondamentali nella teoria del controllo, ampiamente applicate nell'analisi di sistemi, riduzione di modelli, stabilizzazione e controllo ottimale. Sebbene le proprietà di queste equazioni (come l'esistenza, l'unicità e il calcolo numerico delle soluzioni) siano state ampiamente studiate, i metodi tradizionali presentano limitazioni:
- L'equazione di Lyapunov ha soluzioni esplicite (forma integrale per il tempo continuo, serie infinita per il tempo discreto), ma l'equazione di Riccati raramente ha soluzioni esplicite
- Le due classi di equazioni richiedono tipicamente metodi di trattamento diversi, mancando di un quadro unificato
Queste due classi di equazioni sono fondamentali nella teoria dei sistemi lineari, e le proprietà delle loro soluzioni influenzano direttamente l'analisi di stabilità dei sistemi di controllo e la progettazione del controllo ottimale. Stabilire un quadro teorico unificato non solo aiuta ad approfondire la comprensione teorica, ma può anche fornire nuove prospettive per gli algoritmi numerici.
- Metodo di Pates (2024): Considera l'equazione di Lyapunov come un problema di punto fisso su un insieme convesso compatto, utilizzando il teorema del punto fisso di Brouwer per provare l'esistenza della soluzione. Tuttavia, questa tecnica non può essere direttamente estesa all'equazione di Riccati, e il teorema di Brouwer non può garantire l'unicità
- Metodi tradizionali: L'equazione di Lyapunov e l'equazione di Riccati sono trattate separatamente, mancando di una prospettiva unificata di sistemi dinamici
Questo articolo propone di trattare unificatamente queste due classi di equazioni mediante la costruzione di sistemi dinamici, sfruttando i potenti strumenti della teoria dei sistemi positivi per analizzare le proprietà delle soluzioni. L'innovazione chiave è l'introduzione di variabili ausiliarie (funzioni lineari scalari) per omogeneizzare le equazioni, rendendole adatte al quadro dei sistemi positivi.
- Stabilire Nuovi Legami Teorici: Per la prima volta, collegare la teoria dei sistemi positivi alle equazioni di controllo classiche (equazioni di Lyapunov e Riccati), fornendo una prospettiva teorica completamente nuova
- Tecnica di Omogeneizzazione: Proporre la tecnica chiave di omogeneizzazione delle equazioni introducendo funzioni lineari a valori scalari (come tr(P)), consentendo il trattamento unificato di equazioni non omogenee nel quadro dei sistemi positivi
- Costruzione di Sistemi Dinamici Positivi:
- Costruire sistemi lineari strettamente positivi per l'equazione di Lyapunov, provando la convergenza esponenziale
- Costruire sistemi strettamente positivi omogenei per l'equazione algebrica di Riccati, provando la convergenza asintotica
- Quadro Unificato di Risoluzione: Caratterizzare la soluzione dell'equazione come lo stato stazionario di un sistema dinamico positivo opportunamente costruito, fornendo un metodo di risoluzione dalla prospettiva dei sistemi dinamici
- Analisi di Convergenza: Utilizzare l'analisi della metrica di Hilbert per studiare le proprietà di convergenza del sistema dinamico, fornendo garanzie teoriche per gli algoritmi numerici
Equazione di Lyapunov in Tempo Continuo:
A⊤P+PA+C⊤C=0
Equazione di Lyapunov in Tempo Discreto:
A⊤QA−Q+C⊤C=0
Equazione Algebrica di Riccati in Tempo Discreto:
P=A⊤PA−A⊤PB(R+B⊤PB)−1B⊤PA+C⊤C
Equazione Algebrica di Riccati in Tempo Continuo:
A⊤P+PA−PBR−1B⊤P+C⊤C=0
dove P,Q∈S+n (insieme di matrici semidefinite positive), A,B,C,R sono matrici note.
Definizione di Cono: Questo articolo utilizza principalmente due classi di coni:
- Ortante positivo: R+n={x∈Rn:xi≥0,∀i}
- Insieme di matrici semidefinite positive: S+n={X∈Rn×n:X⪰0}
Metrica di Hilbert: Per x,y=0 nel cono K, si definisce
dH(x,y)=lnm(x/y)M(x/y)
dove M(x/y)=inf{λ:x≤λy}, m(x/y)=sup{μ:μy≤x}
Teorema Chiave (Teorema 1 - Teoria Nonlineare di Perron-Frobenius):
Sia K un cono regolare, F:V→V una mappa continua, omogenea e che preserva l'ordine. Allora esiste x∗∈K∖{0} tale che F(x∗)=ρ(F)x∗. Se F preserva fortemente l'ordine, allora x∗ è unica nel senso di scalari positivi, e dH(Fk(x),x∗)→0; se F è inoltre lineare, la convergenza è esponenziale.
Caso in Tempo Continuo:
Passo 1: Costruzione del Sistema Linearizzato
Considerare la semplice equazione lineare Ap+c=0 (A è una matrice di Metzler, c≥0). Il metodo tradizionale costruisce il sistema dinamico x˙=Ax+c, ma questo non è omogeneo. Mediante la tecnica di "linearizzazione":
x˙=(A+αc1⊤)x
dove α>0 è un parametro di regolazione, 1 è il vettore di tutti gli uni.
Passo 2: Costruzione del Sistema Positivo per l'Equazione di LyapunovP˙=A⊤P+PA+αtr(P)C⊤C(5)
Lemma Chiave (Lemma 1): Se (A,C) è osservabile, allora il sistema (5) è strettamente positivo su S+n.
Idea della Dimostrazione: L'espressione della soluzione è
P(t)=eA⊤tP0eAt+α∫0ttr(P(s))eA⊤(t−s)C⊤CeA(t−s)ds
Se esiste v=0 tale che P(t)v=0, allora CeAsv=0 per tutti gli s∈[0,t], il che contraddice l'osservabilità. Pertanto P(t)≻0.
Passo 3: Regolazione dei Parametri e Costruzione della Soluzione
Dal Corollario 1, esiste un unico P∗≻0 e λ∈R tale che
A⊤P∗+P∗A+αtr(P∗)C⊤C=λP∗
Per un autovettore v di A (autovalore η), si ha
λ=2Re(η)+αtr(P∗)v∗P∗v∥Cv∥2
- Quando α è sufficientemente piccolo, λ<0 (poiché A è stabile)
- Quando α è sufficientemente grande, λ>0 (poiché ∥Cv∥2=0)
- Per continuità, esiste α>0 tale che λ=0
In questo caso P#=αtr(P∗)P∗ è l'unica soluzione definita positiva dell'equazione di Lyapunov.
Convergenza: Il sistema (5) converge esponenzialmente nella metrica di Hilbert:
dH(P(t),P#)≤ke−λtdH(P0,P#)
Caso in Tempo Discreto:
Costruire il sistema iterativo
Xk+1=F(Xk)=A⊤XkA+αtr(Xk)C⊤C
Un'analisi simile mostra che, sotto l'osservabilità e la stabilità di Schur, esiste una soluzione unica e il sistema converge esponenzialmente.
Caso in Tempo Discreto:
Passo 1: Costruzione del Sistema OmogeneoXk+1=F(Xk)=A⊤S(Xk)A+αtr(Xk)C⊤C
dove
S(X):=X−XB(αtr(X)R+B⊤XB)−1B⊤X
Lemma Chiave (Lemma 3): L'operatore F è continuo, omogeneo, che preserva l'ordine e concavo. Se (A,C) è osservabile, allora Fn preserva fortemente l'ordine.
Punti Chiave della Dimostrazione:
- La proprietà di preservare l'ordine deriva dalla monotonia del complemento di Schur (Lemma 6)
- La concavità è una proprietà nota del complemento di Schur
- La proprietà di preservare fortemente l'ordine si ottiene mediante argomenti ricorsivi: se v⊤Xk+1v=0, allora Cv=0 e v⊤A⊤XkAv=0, iterando n volte si ottiene
col(C,CA,…,CAn−1)v=0
Per osservabilità si ha v=0, quindi Fn(X)∈IntS+n
Passo 2: Regolazione dei Parametri
- Quando α è piccolo, se A è stabile di Schur (garantito dall'ipotesi di controllabilità), allora ρ(F)<1
- Quando α→∞, ρ(F)→∞ (per osservabilità)
- Per continuità del raggio spettrale del cono, esiste α>0 tale che ρ(F)=1
In questo caso P=αtr(X∗)X∗ è l'unica soluzione definita positiva dell'equazione di Riccati.
Convergenza: Nella metrica di Hilbert converge asintoticamente (non esponenzialmente):
dH(Xk,P)→0 quando k→∞
Caso in Tempo Continuo:
Costruire il sistema
P˙=αA⊤P+αPA−tr(P)1PBR−1B⊤P+α2tr(P)C⊤C
Equivalentemente,
P˙=Aα(t)⊤P+PAα(t)+α2tr(P)C⊤C
dove Aα(t)=αA−2tr(P)1BR−1B⊤P(t)
Lemma Chiave (Lemma 4): Se (A,C) è osservabile, allora il sistema è strettamente positivo.
Un argomento di regolazione dei parametri simile mostra che esiste α>0 tale che il sistema ha un'unica soluzione stazionaria P, e converge asintoticamente nella metrica di Hilbert.
- Tecnica di Omogeneizzazione: Introdurre tr(P) come fattore di scalatura, trasformando l'equazione originale non omogenea in un sistema dinamico omogeneo, che è la chiave per applicare la teoria dei sistemi positivi
- Quadro Unificato: Mediante la costruzione di sistemi dinamici positivi, unificare l'analisi dell'equazione di Lyapunov e dell'equazione di Riccati
- Metodo di Regolazione dei Parametri: Regolando il parametro α per far raggiungere al raggio spettrale del cono un valore specifico (0 per l'equazione di Lyapunov, 1 per l'equazione di Riccati), costruire abilmente la soluzione dell'equazione originale
- Analisi della Metrica di Hilbert: Utilizzare la metrica di Hilbert piuttosto che la metrica euclidea tradizionale per analizzare la convergenza, più adatta alla struttura geometrica dei sistemi positivi
- Nuovo Ruolo dell'Osservabilità: Rivelare il ruolo chiave dell'osservabilità nel garantire la positività stretta del sistema, fornendo una nuova interpretazione geometrica dell'ipotesi classica
Questo articolo è principalmente un lavoro teorico, verificato mediante esempi numerici:
Esempio 1 (Figura 1): Per la semplice equazione lineare Ap+c=0, mostrare la curva dell'autovalore massimo della matrice A+αc1⊤ al variare di α, verificando che esiste α tale che l'autovalore sia 0.
Esempio 2 (Figura 2): Per l'equazione di Lyapunov in tempo discreto, mostrare la curva del raggio spettrale del cono dell'operatore L(X)=A⊤XA+αtr(X)C⊤C al variare di α. Quando A è stabile di Schur e (A,C) è osservabile, esiste α tale che il raggio spettrale sia 1.
- Continuità dei Parametri: Le Figure 1 e 2 mostrano chiaramente la continuità del raggio spettrale/autovalore rispetto al parametro α, supportando l'applicazione del teorema del valore intermedio nella teoria
- Comportamento di Convergenza: Sebbene l'articolo non fornisca curve di convergenza dettagliate, la teoria garantisce la convergenza esponenziale per l'equazione di Lyapunov e la convergenza asintotica per l'equazione di Riccati
Risultati Teorici:
- Fatto 1 (Equazione di Lyapunov in Tempo Continuo): Se (A,C) è osservabile e A è stabile, allora esiste un'unica soluzione definita positiva P#, e il sistema dinamico costruito converge esponenzialmente a P# nella metrica di Hilbert
- Fatto 2 (Equazione di Lyapunov in Tempo Discreto): Se A è stabile di Schur e (A,C) è osservabile, allora esiste un'unica soluzione definita positiva Q, e il sistema converge esponenzialmente
- Fatto 3 (Equazione Algebrica di Riccati in Tempo Discreto): Se (A,B) è controllabile e (A,C) è osservabile, allora esiste un'unica soluzione definita positiva P, e il sistema converge asintoticamente nella metrica di Hilbert
- Fatto 4 (Equazione Algebrica di Riccati in Tempo Continuo): Se (A,B) è controllabile e (A,C) è osservabile, allora esiste un'unica soluzione definita positiva P, e il sistema converge asintoticamente
- Ruolo Duale dell'Osservabilità e della Controllabilità (Osservazione 5):
- L'osservabilità tende ad aumentare il raggio spettrale del cono
- La controllabilità tende a diminuire il raggio spettrale del cono
- Entrambe sono necessarie, lavorando insieme per garantire l'esistenza e l'unicità della soluzione
- Differenze di Convergenza:
- Equazione di Lyapunov: sistema lineare → convergenza esponenziale
- Equazione di Riccati: sistema omogeneo non lineare → convergenza asintotica (ma non esponenziale)
- Estensibilità: Il metodo può essere esteso a equazioni di Lyapunov più generali, come:
- Sistemi stocastici: A⊤P+PA+∑i=1rAi⊤PAi+C⊤C=0
- Sistemi singolari: A⊤PE+E⊤PA+C⊤C=0
Per l'equazione di Bellman dei sistemi di controllo positivi:
p=s+A⊤p−E⊤∣r+B⊤p∣
Sotto le ipotesi H1 (esistenza di un guadagno stabilizzante) e H2 (condizione simile all'osservabilità), costruendo il sistema
pk+1=α(1⊤pk)s+A⊤pk−E⊤∣α(1⊤pk)r+B⊤pk∣
si provano l'esistenza e l'unicità della soluzione, mostrando l'applicazione del metodo nel controllo ottimale di sistemi positivi.
- Letteratura 9,11,18,1,14: Studiano sistematicamente le proprietà di queste equazioni, gli algoritmi numerici e le applicazioni
- Metodi Tradizionali: Soluzioni esplicite (Lyapunov) vs algoritmi iterativi (Riccati)
- Pates 19: Considera l'equazione di Lyapunov come un problema di punto fisso su un insieme convesso compatto, utilizzando il teorema del punto fisso di Brouwer
- Limitazioni: Non può garantire l'unicità, difficile da estendere all'equazione di Riccati
- Rantzer e Valcher 21: Rassegna della teoria dei sistemi positivi
- Lemmens e Nussbaum 16: Teoria nonlineare di Perron-Frobenius
- Applicazioni: Problemi di consenso 24,22, controllo ottimale 20,6
- Bertsekas 5, Lee 15, Lai e Xiong 12: Convergenza dell'iterazione di valore per l'equazione di Riccati
- Contributo di questo Articolo: Mediante omogeneizzazione, rendere l'iterazione di valore applicabile al quadro della teoria dei sistemi positivi
- Rantzer 20, Li e Rantzer 17: Equazione di Bellman esplicita per sistemi di controllo positivi
- Estensione di questo Articolo: Analizzare unificatamente mediante la teoria dei sistemi positivi
- Unificazione Teorica: Successo nel portare l'equazione di Lyapunov e l'equazione algebrica di Riccati nel quadro della teoria dei sistemi positivi, fornendo un metodo di analisi unificato
- Proprietà della Soluzione: Sotto ipotesi standard (osservabilità, controllabilità, stabilità), provare l'esistenza, l'unicità e la convergenza della soluzione
- Caratteristiche di Convergenza:
- Equazione di Lyapunov: convergenza esponenziale (sistema lineare)
- Equazione di Riccati: convergenza asintotica (sistema omogeneo non lineare)
- Contributo Metodologico: La tecnica di omogeneizzazione fornisce un nuovo strumento per il trattamento di equazioni algebriche
- Velocità di Convergenza: Per l'equazione di Riccati, si può garantire solo la convergenza asintotica piuttosto che esponenziale, il che potrebbe portare a una convergenza più lenta nell'implementazione numerica
- Scelta dei Parametri: Sebbene la teoria garantisca l'esistenza di un α appropriato, l'articolo non fornisce un algoritmo per scegliere efficientemente α
- Esperimenti Numerici Mancanti: L'articolo è principalmente un lavoro teorico, mancano esperimenti numerici su larga scala per verificare l'efficienza pratica del metodo e il confronto con gli algoritmi esistenti
- Estensibilità: L'estensione a sistemi tempo-varianti, sistemi non lineari e altri casi più generali non è ancora stata esplorata
- Limiti di Errore: L'articolo non quantifica i limiti di errore e le costanti specifiche della velocità di convergenza degli algoritmi iterativi
Direzioni di ricerca chiaramente indicate nell'articolo:
- Equazioni Tempo-Varianti: Estendere alle equazioni di Lyapunov e Riccati tempo-varianti
- Quantificazione dei Limiti di Errore: Stabilire limiti di errore precisi per i metodi iterativi proposti
- Implementazione Algoritmica: Sviluppare algoritmi numerici efficienti, incluse strategie di scelta adattiva del parametro α
- Estensione delle Applicazioni: Esplorare le applicazioni nel controllo predittivo del modello, controllo robusto e altri campi
- Forte Innovazione Teorica:
- Per la prima volta, stabilire un profondo collegamento tra la teoria dei sistemi positivi e le equazioni di controllo classiche
- La tecnica di omogeneizzazione ha universalità e potrebbe ispirare la ricerca su altre equazioni algebriche
- Fornire una nuova interpretazione geometrica delle ipotesi classiche (osservabilità, controllabilità)
- Unificazione del Metodo:
- Trattare l'equazione di Lyapunov e l'equazione di Riccati con lo stesso quadro
- Trattamento unificato di sistemi in tempo continuo e discreto
- Estensibile a forme di equazioni più generali
- Rigore Matematico:
- Prove complete e logica chiara
- Utilizzo completo dei risultati profondi della teoria dei sistemi positivi (teoria di Perron-Frobenius, metrica di Hilbert)
- Tutti i lemmi chiave hanno prove dettagliate
- Intuizione Teorica Profonda:
- Rivelare il ruolo duale dell'osservabilità e della controllabilità nella regolazione del raggio spettrale
- Spiegare l'essenza delle differenze di convergenza tra l'equazione di Lyapunov e l'equazione di Riccati
- Fornire un nuovo paradigma per comprendere le equazioni algebriche dalla prospettiva dei sistemi dinamici
- Scrittura Chiara:
- Transizione graduale da esempi semplici (equazioni lineari) a casi complessi
- Struttura chiara, progressione logica
- Sistema di simboli unificato, facile da comprendere
- Praticità da Verificare:
- Mancanza di confronto di efficienza con gli algoritmi numerici esistenti (algoritmo di Bartels-Stewart, iterazione di Newton)
- Nessuna strategia pratica fornita per la scelta del parametro α
- Le costanti di convergenza k,λ non hanno metodi di stima
- Esperimenti Numerici Insufficienti:
- Solo due semplici esempi numerici
- Mancanza di test su sistemi ad alta dimensione, sistemi mal condizionati e altri casi complessi
- Nessuna verifica della coerenza tra la velocità di convergenza teorica e quella pratica
- Convergenza dell'Equazione di Riccati:
- Può garantire solo la convergenza asintotica, senza stima quantitativa della velocità di convergenza
- Per le applicazioni pratiche, la convergenza esponenziale è più ideale
- Condizioni di Ipotesi:
- Dipendere dalle ipotesi standard di osservabilità, controllabilità e stabilità
- Non discutere il caso in cui le ipotesi non sono soddisfatte
- Complessità Computazionale:
- Non analizzare la complessità computazionale dell'algoritmo
- Ogni iterazione richiede il calcolo della traccia della matrice, che potrebbe aumentare il costo computazionale
- Generalizzabilità:
- L'estensione a equazioni non omogenee, problemi di ottimizzazione vincolata, ecc. non è ancora chiara
- L'estensione a casi tempo-varianti, stocastici e non lineari richiede ulteriori ricerche
- Contributo Teorico:
- Stabilire un nuovo ponte tra la teoria del controllo e la teoria dei sistemi positivi
- Potrebbe ispirare la ricerca su metodi di sistemi dinamici per altre equazioni algebriche
- La tecnica di omogeneizzazione ha valore metodologico
- Valore Accademico:
- Fornire una nuova prospettiva teorica per l'equazione di Lyapunov e l'equazione di Riccati
- Approfondire la comprensione dei concetti classici come osservabilità e controllabilità
- Potrebbe stimolare ricerche successive sull'applicazione della teoria dei sistemi positivi
- Valore Pratico:
- Valore pratico limitato nel breve termine (richiede ulteriore sviluppo algoritmico)
- A lungo termine, potrebbe fornire nuove prospettive per sistemi su larga scala e calcolo distribuito
- Ha applicazione diretta nel controllo ottimale di sistemi positivi (Sezione 4.3)
- Riproducibilità:
- La prova teorica è completa e facile da verificare
- La descrizione dell'algoritmo è chiara, ma mancano i dettagli di implementazione
- Si consiglia agli autori di rendere pubblico il codice per promuovere la ricerca
- Ricerca Teorica:
- Studiare le proprietà dei sistemi dinamici delle equazioni algebriche
- Esplorare nuove applicazioni della teoria dei sistemi positivi
- Analizzare la struttura geometrica dei sistemi di controllo
- Applicazioni Specifiche:
- Controllo ottimale di sistemi positivi (come modelli epidemiologici, reti di traffico)
- Sistemi sparsi su larga scala (la struttura positiva potrebbe portare vantaggi computazionali)
- Ottimizzazione distribuita (la metrica di Hilbert è adatta all'analisi della coerenza)
- Uso Didattico:
- Mostrare la potenza della teoria dei sistemi positivi
- Illustrare il collegamento tra equazioni algebriche e sistemi dinamici
- Come materiale supplementare per corsi avanzati di teoria del controllo
- 16 Lemmens & Nussbaum (2012): Nonlinear Perron-Frobenius Theory - La principale fonte degli strumenti teorici di questo articolo
- 19 Pates (2024): Lyapunov equations: a (fixed) point of view - Lavoro direttamente correlato
- 21 Rantzer & Valcher (2018): A tutorial on positive systems and large scale control - Rassegna della teoria dei sistemi positivi
- 20 Rantzer (2022): Explicit solution to bellman equation for positive systems - Controllo ottimale di sistemi positivi
- 8 Bushell (1973): Hilbert's metric and positive contraction mappings - Riferimento classico della metrica di Hilbert
Valutazione Complessiva: Questo è un lavoro con forte innovazione teorica, che ha successo nell'applicare la teoria dei sistemi positivi all'analisi delle equazioni di controllo classiche, fornendo una prospettiva teorica completamente nuova e un quadro unificato. La tecnica di omogeneizzazione ha valore metodologico e potrebbe ispirare la ricerca in altri campi. La principale carenza è che la praticità rimane da verificare, richiedendo lavori successivi per integrare l'implementazione algoritmica e gli esperimenti numerici. Per i ricercatori teorici e gli studiosi interessati alla teoria dei sistemi positivi, questo è un articolo che merita una lettura approfondita.