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

Resumo Técnico: Estabilidade Exponencial de Sistemas Positivos Homogêneos de Grau Um com Atrasos Variáveis no Tempo

Definição do Problema
Embora a estabilidade assintótica de sistemas lineares positivos na presença de atrasos limitados no tempo seja bem estabelecida — especificamente a propriedade de que tais sistemas são estáveis independentemente do atraso se, e somente se, seus equivalentes sem atraso forem estáveis — o arcabouço teórico para sistemas positivos não lineares é menos desenvolvido. A literatura existente frequentemente baseia-se em suposições de atraso constante ou técnicas como funcionais de Lyapunov-Krasovskii lineares, que podem levar a um conservadorismo excessivo ou falhar quando os atrasos são variáveis no tempo. Além disso, embora a estabilidade assintótica seja conhecida, há uma carência de caracterizações precisas sobre como os limites máximos de atraso afetam a taxa de decaimento exponencial desses sistemas. Este artigo aborda a lacuna no estabelecimento da estabilidade exponencial independente de atraso e na quantificação das taxas de decaimento para uma classe significativa de sistemas positivos não lineares com atrasos variáveis no tempo e limitados.

Metodologia
O artigo analisa duas classes primárias de sistemas:

1. Sistemas de tempo contínuo: Definidos por $\dot{x}(t) = f(x(t)) + g(x(t - \tau(t)))$, onde $f$ é um campo vetorial cooperativo e $g$ é preservador de ordem. Ambos são assumidos como homogêneos de grau um.
2. Sistemas de tempo discreto: Definidos por $x(k+1) = f(x(k)) + g(x(k - d(k)))$, onde $f$ e $g$ são homogêneos de grau um e preservadores de ordem.

Os autores utilizam as propriedades de campos vetoriais cooperativos (Jacobiano de Metzler) e preservadores de ordem (monótonos). Uma etapa metodológica fundamental envolve o uso de uma norma $l_\infty$ ponderada, $\|x\|_v^\infty = \max_i |x_i|/v_i$, definida por um vetor positivo $v$. A análise de estabilidade baseia-se na busca de um vetor $v > 0$ que satisfaça desigualdades algébricas específicas envolvendo os campos vetoriais avaliados em $v$.

Para sistemas lineares, os autores transformam as condições de estabilidade em problemas de otimização. Ao aplicar uma mudança de variáveis logarítmica ($z_i = \ln v_i$), eles demonstram que a busca pela taxa de decaimento ótima e pelo vetor $v$ associado pode ser formulada como um problema de otimização convexa. Isso permite o cálculo eficiente dos limites mais estreitos da taxa de decaimento fornecidos pelo seu arcabouço teórico.

Principais Contribuições e Resultados

1. Condições Necessárias e Suficientes para Estabilidade Exponencial:

   • Tempo contínuo: O artigo estabelece que um sistema cooperativo homogêneo com atrasos variáveis no tempo limitados é globalmente exponencialmente estável se, e somente se, existir um vetor $v > 0$ tal que $f(v) + g(v) < 0$. Esta condição é idêntica ao caso sem atraso, confirmando que a estabilidade exponencial é independente do atraso para esta classe de sistemas.
   • Tempo discreto: Da mesma forma, para sistemas monótonos homogêneos de tempo discreto, a estabilidade exponencial global ocorre se, e somente se, existir $v > 0$ tal que $f(v) + g(v) < v$.

2. Limites Explícitos de Taxa de Decaimento:

   • Os autores derivam expressões explícitas para quantificar o impacto dos atrasos na taxa de decaimento.
   • Para sistemas de tempo contínuo, a taxa de decaimento $\eta$ é limitada pela solução positiva única de uma equação transcendente envolvendo o limite de atraso $\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$$
   • Para sistemas de tempo discreto, o fator de decaimento $\gamma$ é derivado de:
     $$\left(\frac{f_i(v)}{v_i}\right) + \left(\frac{g_i(v)}{v_i}\right)\gamma_i^{-d_{max}} = \gamma_i$$
   • Os resultados mostram que, à medida que o limite de atraso aumenta, a taxa de decaimento garantida deteriora-se (aproximando-se de zero para tempo contínuo e de um para tempo discreto).

3. Otimização para Sistemas Lineares:

   • Para sistemas lineares positivos (onde $f(x)=Ax, g(x)=Bx$), o artigo demonstra que a melhor taxa de decaimento garantida pode ser encontrada via otimização convexa. As restrições são convexas nas variáveis transformadas, permitindo a solução numérica eficiente.
   • O arcabouço é estendido a sistemas lineares gerais (não necessariamente positivos) utilizando os valores absolutos das matrizes do sistema ($A_M$ e $|B|$), garantindo a estabilidade se $(A_M + |B|)v < 0$.

4. Validação:

   • Os resultados teóricos são validados através de exemplos ilustrativos envolvendo tanto sistemas não lineares quanto lineares (contínuos e discretos). Simulações confirmam que os limites superiores derivados da taxa de decaimento são estreitos e que a abordagem de otimização convexa produz limites significativamente melhores do que o uso de vetores viáveis arbitrários.

Significância e Alegações
O artigo alega estender a propriedade fundamental de estabilidade independente de atraso de sistemas lineares positivos para uma classe mais ampla de sistemas positivos não lineares (homogêneos cooperativos e monótonos) com atrasos variáveis no tempo. Diferente de trabalhos anteriores que focavam apenas em estabilidade assintótica ou atrasos constantes, este artigo fornece:

• Uma prova rigorosa de que a estabilidade exponencial para esta classe de sistemas é independente da magnitude de atrasos variáveis no tempo limitados.
• Limites quantitativos explícitos sobre a taxa de decaimento como uma função da magnitude do atraso.
• Um método construtivo (otimização convexa) para encontrar a taxa de decaimento ótima para sistemas lineares positivos, abordando a falta de caracterizações precisas na literatura existente.

Os autores observam que, embora os resultados para sistemas cooperativos homogêneos estendam o trabalho anterior sobre atrasos constantes, a extensão para atrasos variáveis no tempo e o fornecimento de limites explícitos de taxa de decaimento representam um avanço significativo, particularmente porque as técnicas padrão, como funcionais de Lyapunov-Krasovskii lineares, são frequentemente insuficientes ou excessivamente conservadoras neste contexto.