Linear Lyapunov Functions for Nonlinear Compartmental Systems

Resumen Técnico: Funciones de Lyapunov Lineales para Sistemas Compartimentales No Lineales

Planteamiento del Problema
Esta nota técnica investiga la estabilidad exponencial de la solución nula (el origen) para una amplia clase de sistemas compartimentales no lineales y no autónomos gobernados por ecuaciones diferenciales ordinarias (EDO). Los sistemas compartimentales —ampliamente utilizados en el flujo de tráfico, la epidemiología y el modelado biológico— son conocidos por describir ecuaciones de balance de masa para cantidades no negativas; sin embargo, la literatura existente se ha centrado principalmente en la estabilidad asintótica o en la existencia de puntos de equilibrio. Existe una notable escasez de trabajos que aborden la tasa de convergencia (específicamente la estabilidad exponencial) para sistemas que son simultáneamente no lineales y no autónomos. El artículo tiene como objetivo llenar este vacío estableciendo condiciones suficientes para la estabilidad exponencial y, en casos específicos, derivando una equivalencia entre la atractividad y la existencia de una función de Lyapunov lineal.

Metodología
Los autores emplean funciones de Lyapunov lineales como la herramienta principal para el análisis de estabilidad. El enfoque central consiste en:

1. Formulación del Sistema: El sistema se define como $\dot{x}(t) = f(t, x(t))x(t)$, donde $f$ es una función de matriz que satisface propiedades compartimentales (elementos fuera de la diagonal no negativos y sumas de columnas no positivas).
2. Condición de Lyapunov: La estabilidad se evalúa mediante la Proposición 1, la cual establece que si existen un vector $v \in (0, +\infty)^n$ y una tasa $\lambda > 0$ tales que $v^\top f(t, q) \leq -\lambda v^\top$, la solución nula es exponencialmente estable.
3. Clasificación de Matrices: El artículo introduce el concepto de matrices compartimentales "conectadas aguas abajo" ($D_n$) y define una clase específica de matrices $A_n(a, b)$ basada en un problema de valores propios restringido (EVP).
4. Procedimientos Constructivos:
   • Caso General: Los autores proponen una prueba (Teorema 1) para verificar si un sistema admite una función de Lyapunov lineal. Esto implica comprobar si la matriz del sistema puede permutarse a una forma conectada aguas abajo que satisfaga límites específicos en las sumas de la diagonal y fuera de la diagonal.
   • Caso Especial (Sistemas Perturbados): Para sistemas donde la matriz $f$ puede descomponerse en una matriz compartimental constante invertible $F$ y una perturbación acotada $p$, los autores proporcionan un procedimiento algorítmico (Proposición 3). Este método construye el vector de Lyapunov $v$ ordenando los elementos de un vector derivado de $F^{-1}$, evitando la necesidad de búsquedas exhaustivas de permutación.

Contribuciones Clave y Resultados

• Condiciones Suficientes para la Estabilidad Exponencial: El artículo proporciona condiciones suficientes explícitas (Teorema 1) bajo las cuales un sistema compartimental no lineal admite una función de Lyapunov lineal. Estas condiciones se verifican comprobando los límites en funcionales de matriz específicos ($d_n$ y $u_n$) y, si es necesario, encontrando una matriz de permutación que haga al sistema "conectado aguas abajo".
• Tasas de Decaimiento Basadas en Valores Propios: Los coeficientes de la función de Lyapunov y la tasa de decaimiento exponencial resultante $\lambda$ se obtienen resolviendo un problema de valores propios (EVP) asociado con una matriz canonical conectada aguas abajo $g_n(a, b)$.
• Teorema Inverso y Equivalencia: Para una clase especial de sistemas que involucran una matriz constante invertible $F$ y una perturbación estructurada $p$, los autores establecen un resultado inverso (Teorema 2). Demuestran que los siguientes puntos son equivalentes:
  1. La solución nula es atractiva.
  2. La matriz $F$ es no singular.
  3. Existe una función de Lyapunov lineal que satisface la condición de decaimiento exponencial.
  4. La solución nula es exponencialmente estable.
• Robustez: El análisis demuestra que las condiciones de estabilidad derivadas son robustas ante perturbaciones acotadas en los elementos fuera de la diagonal de la matriz del sistema.

Demostración Numérica
Los resultados teóricos se ilustran mediante un sistema no lineal de 3 dimensiones donde los elementos de la matriz dependen de la norma $L_1$ del estado. Al aplicar el Teorema 1, los autores determinan los parámetros $a$ y $b$ para definir el conjunto $A_3(a, b)$. La resolución del EVP asociado produce una tasa de decaimiento $\lambda = 0.04$ y un vector de Lyapunov $v$. Las simulaciones numéricas confirman que las soluciones que parten de un conjunto acotado satisfacen el límite exponencial derivado, aunque los autores señalan que la tasa calculada es conservadora respecto al intervalo de tiempo observado.

Significancia y Pretensiones
El artículo pretende proporcionar un marco constructivo para analizar la estabilidad exponencial de sistemas compartimentales no lineales y no autónomos, un tema previamente poco explorado.

• En el caso general, el método ofrece una forma sistemática (aunque potencialmente costosa computacionalmente debido a las búsquedas de permutación) de verificar la estabilidad y calcular las tasas de decaimiento.
• En el caso especial de matrices constantes perturbadas, el artículo afirma una simplificación significativa: la matriz de permutación requerida para la construcción de la función de Lyapunov puede encontrarse simplemente ordenando los elementos de un vector, en lugar de mediante una búsqueda exhaustiva.
• Crucialmente, para este caso especial, el artículo establece que la estabilidad exponencial es equivalente a la invertibilidad de la matriz constante subyacente $F$, una condición que es numéricamente fácil de verificar.

Los autores mantienen un alcance modesto, centrándose estrictamente en la derivación teórica de estas condiciones y su aplicación a la clase definida de EDOs, sin extender sus afirmaciones a clases más amplias de sistemas dinámicos ni proponer nuevas aplicaciones experimentales más allá del ejemplo numérico proporcionado.