On the Lyapunov equation with the state matrix in companion form

Resumo Técnico: Sobre a Equação de Lyapunov com a Matriz de Estado em Forma de Companheiro

Enunciado do Problema
Este artigo investiga a não negatividade elemento a elemento da solução única $X$ da equação de Lyapunov em tempo contínuo:
$$XA + A^\top X = -Q$$
onde $A \in \mathbb{R}^{n \times n}$ é uma matriz de companheiro de Hurwitz associada a um polinômio $a(s)$, e $Q = Q^\top \succeq 0$ é uma matriz semidefinida positiva.

Embora a teoria de Lyapunov padrão garanta que $X$ seja simétrica e semidefinida positiva (na ordem de Loewner), essa propriedade não implica que $X$ seja elemento a elemento não negativa ($X_{ij} \geq 0$ para todo $i,j$). O artigo aborda a Conjectura 1.1, proposta pelo Prof. Luigi Fortuna, que postula que, se $A$ é uma matriz de companheiro de Hurwitz e $Q \succeq 0$, então a solução $X$ é elemento a elemento não negativa. Esta propriedade é distinta da semidefinitude positiva e é de particular interesse na teoria de sistemas positivos, onde a não negatividade componente a componente garante a preservação de grandezas físicas como massas, concentrações ou populações.

Metodologia
A estratégia de prova baseia-se na redução do problema geral a uma estrutura algébrica específica envolvendo polinômios simétricos elementares e matrizes do tipo Cauchy. A metodologia procede através das seguintes etapas:

1. Simplificação do Problema:

   • Redução de Posto-1: Pela linearidade da equação de Lyapunov e pela decomposição de qualquer matriz semidefinida positiva em uma soma de matrizes de posto-1 ($Q = \sum q_i q_i^\top$), os autores assumem, sem perda de generalidade, que $Q = qq^\top$.
   • Diagonalizabilidade: Os autores argumentam que a solução $X$ é uma função contínua de $A$. Assim, basta provar o resultado para matrizes de companheiro diagonalizáveis (onde os autovalores são distintos). Casos não diagonalizáveis podem ser aproximados por diagonalizáveis com coeficientes arbitrariamente próximos.

2. Decomposição Espectral:

   • Assumindo que $A$ é diagonalizável com autovalores reais distintos $-\lambda_i$ (ordenados $0 < \lambda_1 < \dots < \lambda_n$), os autores utilizam os autovetores de $A^\top$. Esses autovetores, denotados por $v_j$, são construídos usando polinômios simétricos elementares $\sigma_m^{(j)}$ dos autovalores.
   • O vetor de forçamento $q$ é expandido nesta base própria: $q = \sum a_i v_i$.

3. Representação Integral e Formas Quadráticas:

   • Utilizando a fórmula integral para a solução de Lyapunov, $X = \int_0^\infty e^{A^\top t} Q e^{At} dt$, os autores derivam uma expressão explícita para as entradas $X_{kl}$.
   • A entrada $X_{kl}$ é mostrada como uma forma quadrática $a^\top M a$, onde a matriz $M$ possui 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$ é fatorada como $M = \Delta C \Delta$, onde $\Delta$ é uma matriz diagonal e $C$ é uma matriz do tipo Cauchy com entradas $C_{ij} = \frac{r_i + r_j}{\lambda_i + \lambda_j}$.

4. Verificação de Positividade:

   • A não negatividade de $X_{kl}$ segue se a matriz $M$ for semidefinida positiva. Como $\Delta$ é não singular, isso é equivalente a provar que $C \succeq 0$.
   • Os autores citam resultados estabelecidos [8] confirmando que tais matrizes do tipo Cauchy são semidefinidas positivas.

Principais Contribuições e Resultados

• Teorema 1 (Resultado Principal): O artigo prova a Conjectura 1.1 sob a condição específica de que a matriz de companheiro $A$ possui apenas autovalores reais. Sob esta suposição, a solução única $X$ é elemento a elemento não negativa.
• Análise de Positividade Total: Os autores investigam uma propriedade mais forte: a total não negatividade (todos os menores são não negativos).
  • Eles fornecem um contraexemplo mostrando que $X$ não é geralmente totalmente não negativa, mesmo quando os autovalores são reais.
  • Teorema 2: Eles estabelecem uma condição suficiente para a positividade total: se o vetor de forçamento $q$ possui uma expansão não negativa na base própria de $A^\top$ (especificamente, $q = \sum \alpha_j v_j$ com $\alpha_j \geq 0$), então $X$ é totalmente positiva. Isso é provado fatorando $X$ em um produto de matrizes totalmente positivas (envolvendo matrizes de Vandermonde e matrizes de Cauchy).

Significância e Alegações
O artigo afirma resolver uma conjectura de longa data sobre as propriedades estruturais das soluções de Lyapunov para matrizes de companheiro, especificamente dentro do âmbito dos sistemas positivos.

• Motivação: A não negatividade elemento a elemento de $X$ é motivada por aplicações em controle ótimo para sistemas positivos. Especificamente, para sistemas de entrada única com matrizes de entrada não negativas, uma solução elemento a elemento não negativa da equação algébrica de Riccati associada (reescrita como uma equação de Lyapunov) implica que o regulador ótimo pode ser implementado utilizando hardware com apenas tensões de alimentação positivas.
• Escopo e Limitações: Os autores são modestos quanto à generalidade de sua prova. Eles declaram explicitamente que a suposição de autovalores reais é uma restrição. O artigo deixa em aberto a questão de se a conjectura se mantém para matrizes de companheiro de Hurwitz com autovalores complexos.
• Valor Teórico: O trabalho conecta a teoria de Lyapunov com a teoria de matrizes do tipo Cauchy e polinômios simétricos elementares, demonstrando que a não negatividade elemento a elemento de $X$ surge da semidefinitude positiva de uma classe específica de matrizes duais a matrizes de Kwong (ou anti-Loewner).

Em resumo, o artigo fornece uma prova rigorosa de que as soluções de Lyapunov para matrizes de companheiro de Hurwitz com espectros reais são elemento a elemento não negativas, enquanto esclarece que a positividade total requer restrições adicionais no vetor de forçamento.