Reducción de orden de modelo simpléctico de sistemas puerto-Hamiltoniano

Información Básica

• ID del Artículo: 2203.07751
• Título: Symplectic model order reduction of port-Hamiltonian systems
• Autores: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
• Clasificación: math.OC (Optimización y Control), math.CA (Análisis Clásico)
• Fecha de Publicación: Marzo de 2022 (preimpresión arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2203.07751

Resumen

Este artículo propone un novedoso método de reducción de orden de modelo (ROM) que preserva la estructura para sistemas puerto-Hamiltoniano (pH) lineales e invariantes en el tiempo. El objetivo de la investigación es construir un sistema pH de orden reducido que siga siendo interpretable en el dominio físico del modelo de orden completo original. Específicamente, si el sistema pH de alta dimensión inicial es un circuito, deseamos que el modelo reducido siga siendo interpretable como un circuito electrónico. Para sistemas masa-resorte-amortiguador (MSD), existen métodos ROM que garantizan la preservación de esta estructura específica. Este artículo demuestra que el nuevo método ROM que preserva la estructura basado en el método ROM simpléctico recupera el método Arnoldi de segundo orden conocido en el caso de sistemas MSD. Sin embargo, para modelos pH de circuitos y otros modelos con estructura de bloques similar, el método produce técnicas novedosas de reducción de orden de modelo.

Antecedentes y Motivación de la Investigación

Descripción del Problema

Los sistemas puerto-Hamiltoniano son una generalización del marco Hamiltoniano clásico, extremadamente útiles en problemas de ingeniería, ya que combinan el método Hamiltoniano clásico con un marco poderoso para modelar sistemas físicos abiertos. Sin embargo, en el diseño de controladores, el controlador típicamente tiene el mismo orden que el sistema original, lo que motiva la necesidad de modelos de orden reducido que capturen el comportamiento fundamental del sistema original.

Limitaciones de los Métodos Existentes

1. Problema de Pérdida de Estructura: La aplicación directa de métodos ROM tradicionales pierde la estructura pH
2. Falta de Interpretabilidad Física: Los métodos existentes no pueden garantizar la interpretabilidad del modelo reducido en el dominio físico
3. Destrucción de Estructura de Bloques: Para sistemas pH con estructura de bloques especial, los métodos existentes no pueden preservar esta estructura

Motivación de la Investigación

Cuando se reformulan modelos físicos en forma pH, las matrices J, R y H típicamente poseen una estructura de bloques especial. El objetivo de este artículo es mantener esta estructura de bloques durante el proceso de reducción, asegurando que el modelo reducido siga siendo interpretable como un sistema físico del mismo tipo que el modelo de orden completo.

Contribuciones Principales

1. Introducción del Método Simpléctico-pH: Se propone un novedoso método ROM para sistemas pH basado en ROM simpléctico, aplicable a la forma de estructura de bloques (2), de modo que el modelo reducido preserve la forma de sistema pH del mismo tipo
2. Conexión Teórica: Se demuestra que en el caso especial de sistemas MSD, el método propuesto recupera el método Arnoldi de segundo orden (SO-Arnoldi)
3. Verificación Numérica: Se presentan resultados numéricos en circuitos electrónicos (circuito LRCR), comparándose con otros métodos conocidos en la literatura
4. Ventajas de Estabilidad: Se descubre que el método SO-Arnoldi produce modelos de orden reducido inestables en ciertos casos, mientras que el método propuesto mantiene la estabilidad

Explicación Detallada del Método

Definición de la Tarea

Dado un sistema pH lineal e invariante en el tiempo:

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

donde x(t) ∈ C²ⁿ es el estado, u(t) ∈ Cᵖ es la entrada, y(t) ∈ Cᵖ es la salida. El objetivo es encontrar un sistema de la misma forma pero con dimensión de espacio de estados más pequeña, de modo que las dos funciones de transferencia sean iguales en puntos de frecuencia predeterminados.

Marco Teórico Principal

Fundamentos de Geometría Simpléctica

Definición 1 (Forma Simpléctica): Sea V un espacio vectorial de dimensión par sobre C, una forma simpléctica Ω: V × V → C satisface:

1. Antisimetría: Ω(v₁, v₂) = -Ω(v₂, v₁)
2. No degeneración: Para todo v₁ ∈ V, Ω(v₁, v₂) = 0 si y solo si v₂ = 0

Definición 2 (Aplicación Simpléctica): Una matriz Q se denomina (J₂ₙ, J₂ₖ)-simpléctica si:

Q*J₂ₙQ = J₂ₖ

ROM Simpléctico para el Caso sin Disipación

Teorema 4: Para un sistema pH sin disipación (R = 0), si Q ∈ C²ⁿˣ²ᵏ es una matriz (J₂ₙ⁻¹, J₂ₖ⁻¹)-simpléctica, y satisface:

• ran(B) ⊆ ran(Q)
• (sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

entonces el sistema pH de orden reducido es igual al sistema original en las frecuencias dadas.

ROM Simpléctico para el Caso con Disipación

Teorema 5: Para un sistema pH con disipación, se construye la matriz extendida:

QE := [Q₁  0  ]
      [0   Q₂]

donde Q₂ es (J₂ₙ, J₂ₖ)-simpléctica, y Q₁* es la inversa izquierda de Q₂. Bajo condiciones apropiadas, el sistema de orden reducido satisface las condiciones de interpolación.

Algoritmo de Preservación de Estructura de Bloques

Algoritmo 1 (ROM Simpléctico-pH):

1. Construir V₁ basado en datos de interpolación
2. Calcular V₂ := HV₁
3. Dividir V₁ y V₂ según la estructura de bloques
4. Construir Ψ₁, Ψ₂, Φ₁, Φ₂ mediante la fórmula (27)
5. Construir Q₁ y Q₂
6. Aplicar proyección Petrov-Galerkin

Fórmulas de construcción:

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

Puntos de Innovación Técnica

1. Preservación de Estructura Simpléctica: Mediante aplicaciones simplécticas se preserva la estructura geométrica fundamental del sistema Hamiltoniano
2. Mantenimiento de Estructura Diagonal de Bloques: La construcción especial de Q₁ y Q₂ asegura que se preserve la estructura de bloques con significado físico
3. Condiciones de Interpolación: Se garantiza que el sistema de orden reducido coincida exactamente con el sistema original en puntos de frecuencia predeterminados
4. Marco Unificado: Incluye el método SO-Arnoldi para sistemas MSD como caso especial

Configuración Experimental

Ejemplo Numérico: Circuito LRCR

Se considera el circuito LRCR mostrado en la Figura 1, que contiene n=50 unidades, con parámetros establecidos como:

• Capacitancia: Cᵢ = 1 pF
• Inductancia: Lᵢ = 1 nH
• Resistencia: Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

Métricas de Evaluación

• Error H₂:

  H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)
  

• Error H∞:

  H∞error = max‖G(iω) - Gred(iω)‖₂
  

Métodos de Comparación

1. pH-Arnoldi: Método que preserva la estructura basado en el Teorema 2
2. SO-Arnoldi: Método Arnoldi de segundo orden
3. Simpléctico-pH: Método propuesto en este artículo

Selección de Puntos de Interpolación

Se seleccionan puntos de interpolación uniformemente distribuidos en escala logarítmica en el intervalo 10⁻², 10², por ejemplo, 6 puntos de interpolación: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

Resultados Experimentales

Resultados Principales

De la Figura 3 se puede observar que:

1. Simpléctico-pH vs pH-Arnoldi: Simpléctico-pH proporciona errores H₂ y H∞ más bajos en la mayoría de los casos
2. Simpléctico-pH vs SO-Arnoldi:
   • Rendimiento similar cuando la dimensión de orden reducido <24
   • Rendimiento de SO-Arnoldi se degrada drásticamente cuando la dimensión de orden reducido =32
   • SO-Arnoldi produce sistemas inestables cuando la dimensión de orden reducido es 28 y 36

Análisis de Estabilidad

La Figura 4 muestra que el método SO-Arnoldi produce polos inestables en r=28 y r=36, mientras que el método Simpléctico-pH mantiene la estabilidad en todas las dimensiones de prueba.

Verificación Teórica: Equivalencia en Sistemas MSD

El Teorema 7 demuestra que para sistemas MSD, el método Simpléctico-pH es equivalente al método SO-Arnoldi mediante una transformación unitaria: existe una matriz unitaria Θ tal que zred(t) = Θẑred(t)

Trabajo Relacionado

Desarrollo de ROM para Sistemas pH

1. Trabajos Tempranos: Métodos de sistemas pasivos basados en funciones de transferencia positivas reales
2. Coincidencia de Momentos: Trabajo de Polyuga y van der Schaft
3. Aproximación de Estructura de Dirac: Diferentes enfoques para preservar la estructura
4. Métodos de Descomposición Espectral: Métodos basados en la función de Popov
5. Métodos de Optimización: Métodos basados en optimización de parámetros y optimización Riemanniana

Desarrollo de ROM Simpléctico

1. Sistemas Hamiltonianos Autónomos: Trabajo pionero de Peng y Mohseni
2. Métodos No Lineales: Reducción simpléctica basada en variedades
3. Métodos Impulsados por Datos: Inferencia de operadores e identificación de sistemas

Conclusiones y Discusión

Conclusiones Principales

1. Se desarrolló exitosamente un método ROM que preserva la estructura pH y la estructura de bloques
2. Se unificó teóricamente diferentes métodos ROM para sistemas MSD
3. Los experimentos numéricos verificaron la efectividad y las ventajas de estabilidad del método

Limitaciones

1. Método de Interpolación: Actualmente solo coincide con valores de funciones de transferencia, sin considerar coincidencia de derivadas
2. Restricción de Estructura de Bloques: El método está diseñado específicamente para estructura de bloques particular
3. Complejidad Computacional: El costo computacional de construir matrices simplécticas requiere análisis adicional

Direcciones Futuras

1. Coincidencia de Derivadas: Extender a coincidencia de derivadas de funciones de transferencia para mejorar la calidad del ajuste
2. Selección Adaptativa: Desarrollar estrategias de selección adaptativa de puntos de interpolación
3. Extensión No Lineal: Extender a sistemas pH no lineales
4. Aplicaciones a Gran Escala: Optimizar el algoritmo para manejar sistemas de mayor escala

Evaluación Profunda

Ventajas

1. Rigor Teórico: Fundamentos matemáticos sólidos basados en geometría simpléctica
2. Preservación de Estructura: Preserva simultáneamente la estructura pH y la estructura de bloques con significado físico
3. Universalidad: Unifica métodos existentes como casos especiales dentro de un marco único
4. Practicidad: El método posee implementación algorítmica clara y verificación numérica

Insuficiencias

1. Rango de Aplicabilidad: El método está diseñado específicamente para estructura de bloques particular, con universalidad limitada
2. Eficiencia Computacional: El análisis de complejidad computacional de la construcción de matrices simplécticas es insuficiente
3. Verificación a Gran Escala: Falta verificación numérica de sistemas de mayor escala

Impacto

1. Contribución Teórica: Proporciona una nueva perspectiva teórica para ROM de sistemas pH
2. Valor Práctico: Posee valor directo en aplicaciones de ingeniería como análisis de circuitos
3. Significado Metodológico: Demuestra la importancia de métodos geométricos en ROM

Escenarios de Aplicación

1. Análisis de Circuitos: Modelado de orden reducido de circuitos electrónicos a gran escala
2. Sistemas Mecánicos: Diseño de control de sistemas masa-resorte-amortiguador
3. Sistemas Multifísicos: Sistemas físicos complejos con características de preservación de energía

Referencias

El artículo incluye 31 referencias bibliográficas que abarcan múltiples campos incluyendo teoría de ROM, modelado de sistemas pH, geometría simpléctica y teoría de control, proporcionando una base teórica sólida para la investigación.