Characterization of stability radii for robustly asymptotically stable dissipative Hamiltonian differential-algebraic systems

Resumo Técnico: Caracterização dos Raios de Estabilidade para Sistemas de Equações Diferenciais-Algebraicas Dissipativas Hamiltonianas Robustamente Assintoticamente Estáveis

Enunciação do Problema
Este artigo aborda a estabilidade assintótica robusta de sistemas de equações diferenciais-algebraicas (EDAs) dissipativas hamiltonianas lineares invariantes no tempo da forma $E \dot{x} = (J - R)Qx$. Embora seja bem estabelecido que os autovalores finitos de tais sistemas residem no semiplano esquerdo fechado e que os autovalores puramente imaginários (excluindo o zero) são semissimples, esses sistemas não são necessariamente robustamente assintoticamente estáveis. Pequenas perturbações nas matrizes de coeficiente podem destruir a regularidade, aumentar o índice do sistema ou deslocar autovalores para o eixo imaginário, levando à instabilidade.

O problema central é quantificar a "distância" de um sistema dHDAE robustamente assintoticamente estável dado até a fronteira do conjunto de tais sistemas. Especificamente, os autores buscam caracterizar as menores perturbações estruturadas necessárias para causar:

1. Singularidade (perda de regularidade).
2. Um aumento no índice diferencial-algebraico para maior que um.
3. O aparecimento de autovalores puramente imaginários (perda da estabilidade assintótica).

A análise considera dois tipos de perturbações estruturadas:

• $S_d$ (Dissipativa): Perturbações que preservam a definição de sinal de $E$ e $R$ (isto é, $\Delta E \le 0, \Delta R \le 0$) e a estrutura anti-hermitiana de $J$.
• $S_i$ (Indefinida): Perturbações que preservam as estruturas hermitiana/anti-hermitiana, mas não necessariamente preservam as restrições de definição de sinal nas próprias perturbações, apenas que as matrizes resultantes permaneçam semidefinidas positivas.

Metodologia
Os autores empregam uma combinação de teoria matricial estrutural, análise espectral e técnicas de otimização:

1. Caracterização Estrutural: O artigo utiliza a forma canônica de Kronecker e uma "forma em escada" específica (Lema 8) derivada via transformações unitárias para analisar a regularidade e o índice do feixe dHDAE. A estabilidade assintótica robusta é caracterizada por três condições: regularidade, índice $\le 1$ e uma abscesso espectral negativa (Teorema 11).
2. Redução a Formas Padrão: A análise frequentemente assume que $Q$ é invertível, permitindo que o sistema seja transformado na forma $E \dot{x} = (J - R)x$ com $E=E^* \ge 0$, $R=R^* \ge 0$, e $J=-J^*$.
3. Problemas de Mapeamento: Para calcular a distância até a instabilidade, os autores formulam o problema como encontrar a perturbação de norma mínima $(\Delta E, \Delta J, \Delta R)$ que satisfaça restrições algébricas específicas (por exemplo, tornar uma matriz singular ou criar um autovalor no eixo imaginário). Eles aproveitam soluções de norma mínima para problemas de mapeamento hermitiano e anti-hermitiano (Teoremas 15 e 16).
4. Formulações de Otimização: As distâncias são expressas como problemas de otimização envolvendo quocientes de Rayleigh generalizados. Para casos em que $E$ é fixo (perturbações parciais), caracterizações exatas são derivadas. Para perturbações completas (onde $E$ varia), os autores derivam limites.
   • A distância à singularidade e a um índice mais alto envolve a minimização de somas de normas quadradas de perturbações sujeitas a condições de posto.
   • A distância a autovalores imaginários envolve a minimização sobre um parâmetro de frequência $\omega$ e um vetor $x$.
5. Implementação Numérica: Os problemas de otimização derivados são resolvidos usando uma estratégia de dois níveis: um loop interno resolvendo um problema de autovalor não linear (NEPv) via iterações de campo autoconsistente com deslocamento de nível (SCF), e um loop externo minimizando sobre o parâmetro de frequência $\omega$ usando rotinas padrão de otimização (por exemplo, fminsearch, fmincon).

Principais Contribuições e Resultados

• Caracterizações Exatas para Perturbações Parciais: Quando a matriz $E$ é fixa e apenas $J$ e $R$ são perturbados, o artigo fornece fórmulas exatas para as distâncias estruturadas à singularidade, a um índice mais alto e a autovalores imaginários (Teoremas 24–26).
  • Para perturbações $S_i$ com $R > 0$, a distância à singularidade é exatamente $\sqrt{\lambda_{\min}(N^*(R^2 - J^2)N)}$, onde $N$ abrange o espaço nulo de $E$.
  • A distância a autovalores imaginários para $S_i$ é dada por $\inf_{\omega} \sigma_{\min}([R, i\omega E - J])$.
• Limites para Perturbações Completas: Quando perturbações são permitidas em todas as matrizes ($E, J, R$), os autores derivam limites inferiores e superiores explícitos:
  • Distância à Singularidade: Fórmulas exatas são fornecidas para $S_d$ (Teorema 20) e limites para $S_i$ (Teorema 21).
  • Distância a um Índice Mais Alto: Limites superiores são derivados (Teorema 22) empregando uma abordagem de perturbação em duas etapas (primeiro perturbando $E$ para aumentar a dimensão do espaço nulo, depois perturbando $J$ e $R$). Os autores observam que, devido a termos de perturbação de segunda ordem, estes são atualmente limites superiores e não caracterizações exatas.
  • Distância a Autovalores Imaginários: Limites inferiores são derivados para ambos $S_d$ e $S_i$ (Teorema 18). Estes envolvem otimização complexa sobre quocientes de Rayleigh generalizados e autovalores de feixes matriciais específicos.
• Estrutura Algorítmica: Um algoritmo passo a passo (Seção 2.7) é apresentado para verificar a estabilidade assintótica robusta e calcular os limites de distância associados usando decomposições de álgebra linear padrão (SVD, decomposição espectral).

Significado e Alegações
O artigo alega fornecer a primeira caracterização abrangente dos raios de estabilidade para sistemas dHDAE sob perturbações que preservam a estrutura. Embora o cálculo dos raios de estabilidade para EDAs gerais não estruturadas seja conhecido por ser difícil, os autores demonstram que as propriedades geométricas e algébricas específicas dos sistemas hamiltonianos dissipativos permitem a derivação de fórmulas explícitas e limites computáveis.

O significado reside na capacidade de quantificar a tolerância de modelos físicos (como aqueles em mecânica estrutural, dinâmica de fluidos e fluxo em meios porosos) a erros de modelagem e discretização. Ao determinar a distância até o sistema singular, de alto índice ou instável mais próximo, engenheiros e cientistas podem avaliar a confiabilidade de seus modelos dHDAE.

Os autores são modestos quanto à precisão dos limites para perturbações completas (especificamente para a distância a um índice mais alto), reconhecendo que os resultados atuais são limites superiores devido à complexidade dos termos de perturbação de segunda ordem. Eles identificam a determinação de fórmulas exatas para esses casos e o desenvolvimento de métodos de otimização numérica mais eficientes como direções necessárias para pesquisas futuras.