Exponential Stability of Homogeneous Positive Systems of Degree One With Time-Varying Delays

Sintesi Tecnica: Stabilità Esponenziale di Sistemi Positivi Omogenei di Grado Uno con Ritardi Variabili nel Tempo

Enunciato del Problema
Sebbene la stabilità asintotica dei sistemi lineari positivi in presenza di ritardi limitati sia ben consolidata — nello specifico, la proprietà per cui tali sistemi sono stabili indipendentemente dal ritardo se e solo se i loro controparti privi di ritardo sono stabili — il quadro teorico per i sistemi positivi non lineari è meno sviluppato. La letteratura esistente spesso si basa su ipotesi di ritardo costante o tecniche come i funzionali di Lyapunov-Krasovskii lineari, che possono portare a un eccessivo conservativismo o fallire nell'applicarsi quando i ritardi sono variabili nel tempo. Inoltre, sebbene la stabilità asintotica sia nota, manca una caratterizzazione netta di come i limiti massimi di ritardo influenzino il tasso di decadimento esponenziale di questi sistemi. Questo articolo affronta la lacuna nel determinare la stabilità esponenziale indipendente dal ritardo e nel quantificare i tassi di decadimento per una classe significativa di sistemi positivi non lineari con ritardi limitati e variabili nel tempo.

Metodologia
L'articolo analizza due classi primarie di sistemi:

1. Sistemi a tempo continuo: definiti da $\dot{x}(t) = f(x(t)) + g(x(t - \tau(t)))$, dove $f$ è un campo vettoriale cooperativo e $g$ è preservante l'ordine. Entrambi sono assunti essere omogenei di grado uno.
2. Sistemi a tempo discreto: definiti da $x(k+1) = f(x(k)) + g(x(k - d(k)))$, dove $f$ e $g$ sono omogenei di grado uno e preservanti l'ordine.

Gli autori utilizzano le proprietà dei campi vettoriali cooperativi (Jacobiano Metzler) e preservanti l'ordine (monotoni). Un passaggio metodologico chiave consiste nell'uso di una norma $l_\infty$ pesata, $\|x\|_v^\infty = \max_i |x_i|/v_i$, definita da un vettore positivo $v$. L'analisi della stabilità si basa sul trovare un vettore $v > 0$ che soddisfi specifiche disuguaglianze algebriche coinvolgenti i campi vettoriali valutati in $v$.

Per i sistemi lineari, gli autori trasformano le condizioni di stabilità in problemi di ottimizzazione. Applicando un cambiamento di variabili logaritmico ($z_i = \ln v_i$), dimostrano che la ricerca del tasso di decadimento ottimale e del vettore $v$ associato può essere formulata come un problema di ottimizzazione convessa. Ciò consente il calcolo efficiente dei limiti più stretti del tasso di decadimento forniti dal loro quadro teorico.

Contributi Chiave e Risultati

1. Condizioni Necessarie e Sufficienti per la Stabilità Esponenziale:

   • Tempo continuo: L'articolo stabilisce che un sistema cooperativo omogeneo con ritardi variabili nel tempo limitati è globalmente esponenzialmente stabile se e solo se esiste un vettore $v > 0$ tale che $f(v) + g(v) < 0$. Questa condizione è identica al caso privo di ritardo, confermando che la stabilità esponenziale è indipendente dal ritardo per questa classe di sistemi.
   • Tempo discreto: Similmente, per i sistemi monotoni omogenei a tempo discreto, la stabilità esponenziale globale si ha se e solo se esiste $v > 0$ tale che $f(v) + g(v) < v$.

2. Limiti Espliciti del Tasso di Decadimento:

   • Gli autori derivano espressioni esplicite per quantificare l'impatto dei ritardi sul tasso di decadimento.
   • Per i sistemi a tempo continuo, il tasso di decadimento $\eta$ è limitato dalla soluzione positiva unica di un'equazione trascendente che coinvolge il limite di ritardo $\tau_{max}$:
     $$\left(\frac{f_i(v)}{v_i}\right) + \left(\frac{g_i(v)}{v_i}\right)e^{\eta_i \tau_{max}} + \eta_i = 0$$
   • Per i sistemi a tempo discreto, il fattore di decadimento $\gamma$ è derivato da:
     $$\left(\frac{f_i(v)}{v_i}\right) + \left(\frac{g_i(v)}{v_i}\right)\gamma_i^{-d_{max}} = \gamma_i$$
   • I risultati mostrano che all'aumentare del limite di ritardo, il tasso di decadimento garantito peggiora (avvicinandosi a zero per il tempo continuo e a uno per il tempo discreto).

3. Ottimizzazione per Sistemi Lineari:

   • Per i sistemi lineari positivi (dove $f(x)=Ax, g(x)=Bx$), l'articolo dimostra che il miglior tasso di decadimento garantito può essere trovato tramite ottimizzazione convessa. I vincoli sono convessi nelle variabili trasformate, permettendo una soluzione numerica efficiente.
   • Il quadro è esteso a sistemi lineari generali (non necessariamente positivi) utilizzando i valori assoluti delle matrici del sistema ($A_M$ e $|B|$), garantendo la stabilità se $(A_M + |B|)v < 0$.

4. Validazione:

   • I risultati teorici sono validati attraverso esempi illustrativi che coinvolgono sia sistemi non lineari che lineari (continui e discreti). Le simulazioni confermano che i limiti superiori derivati sul tasso di decadimento sono stretti e che l'approccio di ottimizzazione convessa fornisce limiti significativamente migliori rispetto all'uso di vettori ammissibili arbitrari.

Significato e Rivendicazioni
L'articolo sostiene di estendere la proprietà fondamentale della stabilità indipendente dal ritardo dai sistemi lineari positivi a una classe più ampia di sistemi positivi non lineari (omogenei cooperativi e monotoni) con ritardi variabili nel tempo. A differenza dei lavori precedenti che si concentravano esclusivamente sulla stabilità asintotica o sui ritardi costanti, questo articolo fornisce:

• Una prova rigorosa che la stabilità esponenziale per questa classe di sistemi è indipendente dall'entità dei ritardi limitati variabili nel tempo.
• Limiti quantitativi espliciti sul tasso di decadimento in funzione dell'entità del ritardo.
• Un metodo costruttivo (ottimizzazione convessa) per trovare il tasso di decadimento ottimale per i sistemi lineari positivi, affrontando la mancanza di caratterizzazioni nette nella letteratura esistente.

Gli autori osservano che, sebbene i risultati per i sistemi cooperativi omogenei estendano il lavoro precedente sui ritardi costanti, l'estensione ai ritardi variabili nel tempo e la fornitura di espliciti limiti del tasso di decadimento rappresentano un avanzamento significativo, poiché le tecniche standard come i funzionali di Lyapunov-Krasovik lineari sono spesso insufficienti o eccessivamente conservative in questo contesto.