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

Resumen Técnico: Caracterización de los Radios de Estabilidad para Sistemas de Ecuaciones Diferenciales-Algebraicas Hamiltonianas Disipativas Asintóticamente Estables de Manera Robusta

Enunciado del Problema
Este artículo aborda la estabilidad asintótica robusta de sistemas de ecuaciones diferenciales-algebraicas (EDAs) Hamiltonianas disipativas lineales e invariantes en el tiempo de la forma $E \dot{x} = (J - R)Qx$. Si bien está bien establecido que los eigenvalores finitos de tales sistemas se encuentran en el semiplano izquierdo cerrado y que los eigenvalores puramente imaginarios (excluyendo el cero) son semisimples, estos sistemas no son necesariamente robustamente asintóticamente estables. Pequeñas perturbaciones en las matrices de coeficientes pueden destruir la regularidad, aumentar el índice del sistema o desplazar los eigenvalores al eje imaginario, lo que conduce a la inestabilidad.

El problema central consiste en cuantificar la "distancia" desde un sistema dHDAE dado que es asintóticamente estable de manera robusta hasta el borde del conjunto de dichos sistemas. Específicamente, los autores buscan caracterizar las perturbaciones estructuradas más pequeñas requeridas para causar:

1. Singularidad (pérdida de regularidad).
2. Un aumento en el índice diferencial-algebraico a un valor mayor que uno.
3. La aparición de eigenvalores puramente imaginarios (pérdida de estabilidad asintótica).

El análisis considera dos tipos de perturbaciones estructuradas:

• $S_d$ (Disipativa): Perturbaciones que preservan la definición de signo de $E$ y $R$ (es decir, $\Delta E \le 0, \Delta R \le 0$) y la estructura anti-Hermítica de $J$.
• $S_i$ (Indefinida): Perturbaciones que preservan las estructuras Hermíticas/anti-Hermíticas pero que no necesariamente preservan las restricciones de definición de signo en las propias perturbaciones, solo que las matrices resultantes permanezcan semidefinidas positivas.

Metodología
Los autores emplean una combinación de teoría matricial estructural, análisis espectral y técnicas de optimización:

1. Caracterización Estructural: El artículo utiliza la forma canónica de Kronecker y una "forma escalonada" específica (Lema 8) derivada mediante transformaciones unitarias para analizar la regularidad y el índice del lápiz dHDAE. La estabilidad asintótica robusta se caracteriza mediante tres condiciones: regularidad, índice $\le 1$ y un abscisa espectral negativa (Teorema 11).
2. Reducción a Formas Estándar: El análisis a menudo asume que $Q$ es invertible, permitiendo transformar el sistema a la forma $E \dot{x} = (J - R)x$ con $E=E^* \ge 0$, $R=R^* \ge 0$, y $J=-J^*$.
3. Problemas de Mapeo: Para calcular la distancia a la inestabilidad, los autores formulan el problema como la búsqueda de la perturbación de norma mínima $(\Delta E, \Delta J, \Delta R)$ que satisfaga restricciones algebraicas específicas (por ejemplo, hacer que una matriz sea singular o crear un eigenvalor en el eje imaginario). Aprovechan soluciones de norma mínima a problemas de mapeo Hermítico y anti-Hermítico (Teoremas 15 y 16).
4. Formulaciones de Optimización: Las distancias se expresan como problemas de optimización que involucran cocientes de Rayleigh generalizados. Para casos en los que $E$ es fijo (perturbaciones parciales), se derivan caracterizaciones exactas. Para perturbaciones completas (donde $E$ varía), los autores derivan cotas.
   • La distancia a la singularidad y a un índice mayor implica minimizar sumas de normas al cuadrado de las perturbaciones sujetas a condiciones de rango.
   • La distancia a eigenvalores imaginarios implica minimizar sobre un parámetro de frecuencia $\omega$ y un vector $x$.
5. Implementación Numérica: Los problemas de optimización derivados se resuelven utilizando una estrategia de dos niveles: un bucle interno que resuelve un problema de eigenvalor no lineal (NEPv) mediante iteraciones de campo autoconsistente con desplazamiento de nivel (SCF), y un bucle externo que minimiza sobre el parámetro de frecuencia $\omega$ utilizando rutinas estándar de optimización (por ejemplo, fminsearch, fmincon).

Contribuciones y Resultados Clave

• Caracterizaciones Exactas para Perturbaciones Parciales: Cuando la matriz $E$ es fija y solo $J$ y $R$ son perturbados, el artículo proporciona fórmulas exactas para las distancias estructuradas a la singularidad, índice mayor y eigenvalores imaginarios (Teoremas 24–26).
  • Para perturbaciones $S_i$ con $R > 0$, la distancia a la singularidad es exactamente $\sqrt{\lambda_{\min}(N^*(R^2 - J^2)N)}$, donde $N$ abarca el espacio nulo de $E$.
  • La distancia a eigenvalores imaginarios para $S_i$ viene dada por $\inf_{\omega} \sigma_{\min}([R, i\omega E - J])$.
• Cotas para Perturbaciones Completas: Cuando se permiten perturbaciones en todas las matrices ($E, J, R$), los autores derivan cotas inferiores y superiores explícitas:
  • Distancia a la Singularidad: Se proporcionan fórmulas exactas para $S_d$ (Teorema 20) y cotas para $S_i$ (Teorema 21).
  • Distancia a Índice Mayor: Se derivan cotas superiores (Teorema 22) empleando un enfoque de perturbación en dos etapas (primero perturbando $E$ para aumentar la dimensión del espacio nulo, luego perturbando $J$ y $R$). Los autores señalan que, debido a los términos de perturbación de segundo orden, estas son actualmente cotas superiores en lugar de caracterizaciones exactas.
  • Distancia a Eigenvalores Imaginarios: Se derivan cotas inferiores tanto para $S_d$ como para $S_i$ (Teorema 18). Estas implican una optimización compleja sobre cocientes de Rayleigh generalizados y eigenvalores de lápices matriciales específicos.
• Marco Algorítmico: Se presenta un algoritmo paso a paso (Sección 2.7) para verificar la estabilidad asintótica robusta y calcular las cotas de distancia asociadas utilizando descomposiciones de álgebra lineal estándar (SVD, descomposición espectral).

Significado y Afirmaciones
El artículo afirma proporcionar la primera caracterización integral de los radios de estabilidad para sistemas dHDAE bajo perturbaciones que preservan la estructura. Si bien se sabe que el cálculo de los radios de estabilidad para DAEs generales no estructuradas es difícil, los autores demuestran que las propiedades geométricas y algebraicas específicas de los sistemas Hamiltonianos disipativos permiten la derivación de fórmulas explícitas y cotas computables.

El significado radica en la capacidad de cuantificar la tolerancia de los modelos físicos (como los de la mecánica estructural, la dinámica de fluidos y el flujo en medios porosos) a los errores de modelado y discretización. Al determinar la distancia al sistema singular, de alto índice o inestable más cercano, los ingenieros y científicos pueden evaluar la fiabilidad de sus modelos dHDAE.

Los autores son modestos respecto a la precisión de las cotas para perturbaciones completas (específicamente para la distancia a un índice mayor), reconociendo que los resultados actuales son cotas superiores debido a la complejidad de los términos de perturbación de segundo orden. Identifican la determinación de fórmulas exactas para estos casos y el desarrollo de métodos de optimización numérica más eficientes como direcciones necesarias para la investigación futura.