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

Resumen Técnico: Estabilidad Exponencial de Sistemas Positivos Homogéneos de Grado Uno con Retardos Variantes en el Tiempo

Planteamiento del Problema
Si bien la estabilidad asintótica de los sistemas lineales positivos en presencia de retardos acotados está bien establecida —específicamente la propiedad de que tales sistemas son estables de forma independiente al retardo si y solo si sus contrapartes sin retardo son estables—, el marco teórico para los sistemas positivos no lineales está menos desarrollado. La literatura existente suele depender de supuestos de retardo constante o técnicas como funcionales de Lyapunov-Krasovskii lineales, que pueden conducir a un conservadurismo excesivo o fallar cuando los retardos son variantes en el tiempo. Además, aunque la estabilidad asintótica es conocida, existe una falta de caracterizaciones precisas sobre cómo los límites máximos de retardo afectan la tasa de decaimiento exponencial de estos sistemas. Este artículo aborda la brecha en el establecimiento de la estabilidad exponencial independiente del retardo y en la cuantificación de las tasas de decaimiento para una clase significativa de sistemas positivos no lineales con retardos variantes en el tiempo y acotados.

Metodología
El artículo analiza dos clases primarias de sistemas:

1. Sistemas en tiempo continuo: Definidos por $\dot{x}(t) = f(x(t)) + g(x(t - \tau(t)))$, donde $f$ es un campo vectorial cooperativo y $g$ es preservador de orden. Ambos se asumen homogéneos de grado uno.
2. Sistemas en tiempo discreto: Definidos por $x(k+1) = f(x(k)) + g(x(k - d(k)))$, donde $f$ y $g$ son homogéneos de grado uno y preservadores de orden.

Los autores utilizan las propiedades de los campos vectoriales cooperativos (Jacobiano de Metzler) y preservadores de orden (monótonos). Un paso metodológico clave consiste en el uso de una norma $l_\infty$ ponderada, $\|x\|_v^\infty = \max_i |x_i|/v_i$, definida por un vector positivo $v$. El análisis de estabilidad se basa en encontrar un vector $v > 0$ que satisfaga desigualdades algebraicas específicas involucrando los campos vectoriales evaluados en $v$.

Para los sistemas lineales, los autores transforman las condiciones de estabilidad en problemas de optimización. Al aplicar un cambio de variables logarítmico ($z_i = \ln v_i$), demuestran que la búsqueda de la tasa de decaimiento óptima y el vector $v$ asociado puede formularse como un problema de optimización convexa. Esto permite la computación eficiente de los límites de tasa de decaimiento más ajustados proporcionados por su marco teórico.

Contribuciones Clave y Resultados

1. Condiciones Necesarias y Suficientes para la Estabilidad Exponencial:

   • Tiempo continuo: El artículo establece que un sistema cooperativo homogéneo con retardos variantes en el tiempo y acotados es globalmente exponencialmente estable si y solo si existe un vector $v > 0$ tal que $f(v) + g(v) < 0$. Esta condición es idéntica a la del caso sin retardo, confirmando que la estabilidad exponencial es independiente del retardo para esta clase de sistemas.
   • Tiempo discreto: De manera similar, para sistemas monótonos homogéneos en tiempo discreto, la estabilidad exponencial global se cumple si y solo si existe $v > 0$ tal que $f(v) + g(v) < v$.

2. Límites Explícitos de la Tasa de Decaimiento:

   • Los autores derivan expresiones explícitas para cuantificar el impacto de los retardos en la tasa de decaimiento.
   • Para sistemas en tiempo continuo, la tasa de decaimiento $\eta$ está acotada por la solución positiva única de una ecuación trascendente que involucra el límite de retardo $\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 en tiempo discreto, el factor de decaimiento $\gamma$ se deriva de:
     $$\left(\frac{f_i(v)}{v_i}\right) + \left(\frac{g_i(v)}{v_i}\right)\gamma_i^{-d_{max}} = \gamma_i$$
   • Los resultados muestran que a medida que el límite de retardo aumenta, la tasa de decaimiento garantizada se deteriora (acercándose a cero para tiempo continuo y a uno para tiempo discreto).

3. Optimización para Sistemas Lineales:

   • Para sistemas lineales positivos (donde $f(x)=Ax, g(x)=Bx$), el artículo demuestra que la mejor tasa de decaimiento garantizada puede encontrarse mediante optimización convexa. Las restricciones son convexas en las variables transformadas, lo que permite una solución numérica eficiente.
   • El marco se extiende a sistemas lineales generales (no necesariamente positivos) utilizando los valores absolutos de las matrices del sistema ($A_M$ y $|B|$), asegurando la estabilidad si $(A_M + |B|)v < 0$.

4. Validación:

   • Los resultados teóricos se validan mediante ejemplos ilustrativos que involucran tanto sistemas no lineales como lineales (continuos y discretos). Las simulaciones confirman que los límites superiores derivados de la tasa de decaimiento son ajustados y que el enfoque de optimización convexa produce límites significativamente mejores que el uso de vectores factibles arbitrarios.

Significancia y Reivindicaciones
El artículo afirma extender la propiedad fundamental de la estabilidad independiente del retardo de los sistemas lineales positivos a una clase más amplia de sistemas positivos no lineales (homogéneos cooperativos y monótonos) con retardos variantes en el tiempo. A diferencia de trabajos previos que se centraban únicamente en la estabilidad asintótica o en retardos constantes, este artículo proporciona:

• Una prueba rigurosa de que la estabilidad exponencial para esta clase de sistemas es independiente de la magnitud de los retardos variantes en el tiempo y acotados.
• Límites cuantitativos explícitos sobre la tasa de decaimiento como función de la magnitud del retardo.
• Un método constructivo (optimización convexa) para encontrar la tasa de decaimiento óptima para sistemas lineales positivos, abordando la falta de caracterizaciones precisas en la literatura existente.

Los autores señalan que, si bien los resultados para sistemas cooperativos homogéneos extienden el trabajo previo sobre retardos constantes, la extensión a retardos variantes en el tiempo y la provisión de límites explícitos de la tasa de decaimiento representan un avance significativo, particularmente porque las técnicas estándar como los funcionales de Lyapunov-Krasovskii lineales suelen ser insuficientes o excesivamente conservadoras en este contexto.