On the Lyapunov equation with the state matrix in companion form

Resumen Técnico: Sobre la ecuación de Lyapunov con la matriz de estado en forma de compañía

Planteamiento del problema
Este artículo investiga la no negatividad elemento a elemento de la solución única $X$ de la ecuación de Lyapunov en tiempo continuo:
$$XA + A^\top X = -Q$$
donde $A \in \mathbb{R}^{n \times n}$ es una matriz de compañía Hurwitz asociada con un polinomio $a(s)$, y $Q = Q^\top \succeq 0$ es una matriz semidefinida positiva.

Si bien la teoría de Lyapunov estándar garantiza que $X$ es simétrica y semidefinida positiva (en el orden de Loewner), esta propiedad no implica que $X$ sea elemento a elemento no negativa ($X_{ij} \geq 0$ para todo $i,j$). El artículo aborda la Conjetura 1.1, propuesta por el Prof. Luigi Fortuna, la cual postula que si $A$ es una matriz de compañía Hurwitz y $Q \succeq 0$, entonces la solución $X$ es elemento a elemento no negativa. Esta propiedad es distinta de la semidefinición positiva y es de especial interés en la teoría de sistemas positivos, donde la no negatividad componente a componente asegura la preservación de magnitudes físicas como masas, concentraciones o poblaciones.

Metodología
La estrategia de prueba se basa en la reducción del problema general a una estructura algebraica específica que involucra polinomios simétricos elementales y matrices de tipo Cauchy. La metodología procede a través de los siguientes pasos:

1. Simplificación del problema:

   • Reducción de rango 1: Mediante la linealidad de la ecuación de Lyapunov y la descomposición de cualquier matriz semidefinida positiva en una suma de matrices de rango 1 ($Q = \sum q_i q_i^\top$), los autores asumen sin pérdida de generalidad que $Q = qq^\top$.
   • Diagonalizabilidad: Los autores argumentan que la solución $X$ es una función continua de $A$. Por lo tanto, basta con probar el resultado para matrices de compañía diagonalizables (donde los autovalores son distintos). Los casos no diagonalizables pueden aproximarse mediante matrices diagonalizables con coeficientes arbitrariamente cercanos.

2. Descomposición espectral:

   • Asumiendo que $A$ es diagonalizable con autovalores reales distintos $-\lambda_i$ (ordenados $0 < \lambda_1 < \dots < \lambda_n$), los autores utilizan los autovectores de $A^\top$. Estos autovectores, denotados como $v_j$, se construyen utilizando polinomios simétricos elementales $\sigma_m^{(j)}$ de los autovalores.
   • El vector de forzamiento $q$ se expande en este basis de autovectores: $q = \sum a_i v_i$.

3. Representación integral y formas cuadráticas:

   • Utilizando la fórmula integral para la solución de Lyapunov, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, los autores derivan una expresión explícita para las entradas $X_{kl}$.
   • Se demuestra que la entrada $X_{kl}$ es una forma cuadrática $a^\top M a$, donde la matriz $M$ tiene elementos definidos por:
     $$M_{ij} = \frac{\sigma_{n-k}^{(i)}\sigma_{n-l}^{(j)} + \sigma_{n-k}^{(j)}\sigma_{n-l}^{(i)}}{\lambda_i + \lambda_j}$$
   • Esta matriz $M$ se factoriza como $M = \Delta C \Delta$, donde $\Delta$ es una matriz diagonal y $C$ es una matriz de tipo Cauchy con entradas $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$.

4. Verificación de la positividad:

   • La no negatividad de $X_{kl}$ se sigue si la matriz $M$ es semidefinida positiva. Dado que $\Delta$ es no singular, esto es equivalente a probar que $C \succeq 0$.
   • Los autores citan resultados establecidos [8] que confirman que tales matrices de tipo Cauchy son semidefinidas positivas.

Contribuciones clave y resultados

• Teorema 1 (Resultado principal): El artículo demuestra la Conjetura 1.1 bajo la condición específica de que la matriz de compañía $A$ tiene solo autovalores reales. Bajo este supuesto, la solución única $X$ es elemento a elemento no negativa.
• Análisis de positividad total: Los autores investigan una propiedad más fuerte: la positividad total (todos los menores son no negativos).
  • Proporcionan un contraejemplo que muestra que $X$ no es generalmente totalmente no negativa, incluso cuando los autovalores son reales.
  • Teorema 2: Establecen una condición suficiente para la positividad total: si el vector de forzamiento $q$ tiene una expansión no negativa en el basis de autovectores de $A^\top$ (específicamente, $q = \sum \alpha_j v_j$ con $\alpha_j \geq 0$), entonces $X$ es totalmente positiva. Esto se demuestra factorizando $X$ en un producto de matrices totalmente positivas (que involucran matrices de Vandermonde y matrices de Cauchy).

Significancia y pretensiones
El artículo afirma resolver una conjetura de larga data sobre las propiedades estructurales de las soluciones de Lyapunov para matrices de compañía, específicamente dentro del ámbito de los sistemas positivos.

• Motivación: La no negatividad elemento a elemento de $X$ está motivada por aplicaciones en control óptimo para sistemas positivos. Específicamente, para sistemas de entrada única con matrices de entrada no negativas, una solución elemento a elemento no negativa de la ecuación de Riccati algebraica asociada (reescrita como ecuación de Lyapunov) implica que el regulador óptimo puede implementarse utilizando hardware con solo suministros de voltaje positivos.
• Alcance y limitaciones: Los autores son modestos respecto a la generalidad de su prueba. Indican explícitamente que el supuesto de autovalores reales es una restricción. El artículo deja abierta la cuestión de si la conjetura se mantiene para matrices de compañía Hurwitz con autovalores complejos.
• Valor teórico: El trabajo conecta la teoría de Lyapunov con la teoría de matrices de tipo Cauchy y los polinomios simétricos elementales, demostrando que la no negatividad elemento a elemento de $X$ surge de la semidefinición positiva de una clase específica de matrices duales a las matrices de Kwong (o anti-Loewner).

En resumen, el artículo proporciona una prueba rigurosa de que las soluciones de Lyapunov para matrices de compañía Hurwitz con espectros reales son elemento a elemento no negativas, mientras aclara que la positividad total requiere restricciones adicionales sobre el vector de forzamiento.