KLAP: Aproximación de bajo rango basada en lema KYP para pasivación óptima en $\mathcal{H}_2$

Información Básica

• ID del Artículo: 2501.05178
• Título: KLAP: KYP lemma based low-rank approximation for $\mathcal{H}_2$-optimal passivation
• Autores: Jonas Nicodemus, Matthias Voigt, Serkan Gugercin, Benjamin Unger
• Clasificación: math.OC (Optimización y Control Matemático)
• Fecha de Publicación: 14 de octubre de 2025 (preimpresión arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2501.05178

Resumen

Este artículo propone un nuevo método de imposición de pasividad denominado KLAP, para la pasivación de sistemas lineales invariantes en el tiempo basado en el lema de Kalman-Yakubovich-Popov (KYP) y las ecuaciones de Lur'e asociadas. El problema de pasivación en este marco corresponde a encontrar una perturbación de un sistema no pasivo dado que lo haga pasivo, minimizando simultáneamente la distancia en norma $\mathcal{H}_2$ o $\mathcal{H}_2$ ponderada en frecuencia entre el sistema no pasivo original y el sistema pasivo resultante. Se demuestra que el problema puede formularse como un problema de optimización sin restricciones cuya función objetivo es eficientemente diferenciable incluso en configuraciones de gran escala. Se prueba que cualquier minimizador del problema sin restricciones produce el mismo sistema pasivo, y que en ausencia de término de realimentación directa, cada minimizador local es también un minimizador global.

Antecedentes y Motivación de la Investigación

Importancia del Problema

1. Necesidades de Modelado de Sistemas Físicos: Numerosos campos físicos como sistemas de circuitos, sistemas de potencia, sistemas mecánicos y elasticidad porosa requieren garantías de pasividad para obtener resultados de simulación físicamente significativos
2. Requisitos de Interconexión en Red: Los sistemas pasivos como bloques constructivos de modelos de red de gran escala, cuando están interconectados preservando potencia, producen un sistema pasivo global
3. Desafíos de Modelado Práctico: Aunque los procesos físicos subyacentes son pasivos, los modelos obtenidos mediante métodos no estructurados de reducción de orden o técnicas de identificación de sistemas basadas en datos frecuentemente resultan no pasivos

Limitaciones de Métodos Existentes

Los métodos de pasivación existentes se dividen principalmente en tres categorías:

1. Métodos LMI Basados en Lema KYP: El costo computacional crece rápidamente con la escala del sistema debido a la necesidad de existencia de matrices de Lyapunov
2. Métodos Basados en Características Espectrales de Matrices Hamiltonianas: Carecen de garantías de convergencia y pueden requerir múltiples iteraciones
3. Métodos Basados en Frecuencias Discretas: Solo pueden garantizar pasividad en rangos de frecuencia específicos

Motivación de la Investigación

Este artículo tiene como objetivo desarrollar un método de pasivación eficiente que pueda:

• Procesar sistemas de gran escala
• Proporcionar garantías de convergencia
• Encontrar soluciones óptimas en el sentido de la norma $\mathcal{H}_2$

Contribuciones Principales

1. Parametrización Explícita: Utilizando la existencia de soluciones de rango mínimo de desigualdades KYP, se obtiene una parametrización explícita de cualquier sistema pasivo con $nm$ variables de decisión
2. Reformulación de Optimización Sin Restricciones: Se reformula el problema de optimización convexa con restricciones como un problema de optimización no convexa sin restricciones, estableciendo métodos de resolubilidad, unicidad y cálculo de gradientes
3. Teoría de Optimalidad Global: Se prueba que en el caso de términos de realimentación directa sesgosimétricos ($D + D^T = 0$), cualquier minimizador local es también un minimizador global
4. Detección de Óptimos Locales: Se proporciona un nuevo criterio utilizando soluciones extremales de desigualdades KYP para verificar si un minimizador local es un minimizador global
5. Estrategias Algorítmicas Prácticas: Se proponen estrategias de inicialización basadas en modificación del término de realimentación directa y estrategias de reinicio

Explicación Detallada del Método

Definición de la Tarea

Dado un sistema dinámico lineal invariante en el tiempo: $\Sigma : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = Cx(t) + Du(t) \end{cases}$

El objetivo es encontrar un sistema modificado: $\hat{\Sigma}(\hat{C}) : \begin{cases} \dot{x}(t) = Ax(t) + Bu(t) \\ y(t) = \hat{C}x(t) + Du(t) \end{cases}$

tal que $\hat{\Sigma}(\hat{C})$ sea pasivo y se minimice la distancia en norma $\mathcal{H}_2$ con respecto al sistema original.

Fundamentos Teóricos Principales

Lema KYP y Parametrización de Pasividad

Basado en el lema KYP, un sistema es pasivo si y solo si existen matrices $L \in \mathbb{R}^{n \times m}$ y $M \in \mathbb{R}^{m \times m}$ tales que: $C = B^T\mathcal{L}^{-1}(-LL^T) + ML^T$$D + D^T = MM^T$

donde $\mathcal{L}$ es el operador de Lyapunov: $\mathcal{L}(X) = A^TX + XA$.

Función Objetivo y Gradiente

La función objetivo puede expresarse como: $J(L) = \text{tr}((C - \hat{C}(L))P(C^T - \hat{C}(L)^T))$

donde $P$ es la Gramiana de controlabilidad. El gradiente es: $\nabla J(L) = 2XL - 2P(C^T - \hat{C}(L)^T)M$

Arquitectura del Algoritmo

Flujo del Algoritmo KLAP

1. Inicialización: Obtener $L_0$ inicial utilizando el Algoritmo 1
2. Optimización: Resolver el problema sin restricciones utilizando L-BFGS
3. Detección de Globalidad: Verificar los valores propios de $Y^* = A - B(D+D^T)^{-1}M(L^*)^T$
4. Estrategia de Reinicio: Si se detecta un óptimo local, ejecutar un paso de gradiente y reiniciar

Estrategia de Inicialización

Perturbar el término de realimentación directa $D$ para hacer el sistema pasivo:

• Calcular $\lambda_{\min} = \min_\omega \lambda_{\min}(\Phi(i\omega))$
• Establecer $D_{\text{pert}} = D - (\lambda_{\min}/2 - \epsilon)I_m$
• Resolver la ecuación algebraica de Riccati correspondiente para obtener la inicialización

Configuración Experimental

Sistemas de Prueba

1. Problema de Referencia ACC: Sistema de pequeña escala ($n=4, m=1$)
2. Brazo de Reproductor de CD: Sistema de escala media ($n=120, m=2$)
3. Enlace de Interconexión de Teléfono Inteligente de Alta Velocidad: Sistema de gran escala ($n=800, m=4$)

Métodos de Comparación

• LMI: Método LMI estándar basado en lema KYP
• LMI-TP: Método LMI con parametrización de traza
• Método Hamiltoniano: Método basado en perturbación de valores propios Hamiltonianos

Indicadores de Evaluación

• Error en norma $\mathcal{H}_2$: \|G - \hat{G}(\cdot; \hat{C})\|_{\mathcal{H}_2
• Tiempo computacional e iteraciones
• Tasa de éxito de convergencia al óptimo global

Resultados Experimentales

Resultados Principales

Comparación de Eficiencia Computacional

Hallazgos Clave

1. Eficiencia Computacional: KLAP es 1-2 órdenes de magnitud más rápido que métodos LMI tradicionales
2. Convergencia Global: En ausencia de término de realimentación directa, todos los óptimos locales son óptimos globales
3. Efectividad de la Estrategia de Reinicio: La estrategia de reinicio puede recuperarse exitosamente de óptimos locales no globales
4. Aplicabilidad a Gran Escala: Funciona efectivamente incluso en sistemas de dimensión 800

Análisis de Casos

Problema de Referencia ACC

• Sin término de realimentación directa: Todas las inicializaciones convergen al óptimo global
• Con término de realimentación directa: El 40% de inicializaciones aleatorias convergen a óptimos locales no globales
• Utilizando estrategia de reinicio: Todas las inicializaciones convergen al óptimo global

Enlace de Interconexión de Teléfono Inteligente de Alta Velocidad

• Mejora del error en norma $\mathcal{H}_2$ de aproximadamente 31% en comparación con métodos de referencia
• Mediante transformación diagonalizante, el tiempo de resolución de una única ecuación de Lyapunov se reduce de 550 ms a 4 ms

Trabajo Relacionado

Clasificación de Métodos de Pasivación

1. Métodos Basados en Lema KYP: Producen problemas de optimización convexa pero con alto costo computacional
2. Métodos Basados en Espectro Hamiltoniano: Carecen de garantías de convergencia
3. Métodos Basados en Muestreo de Frecuencia: Solo son efectivos en rangos de frecuencia específicos

Ventajas del Presente Trabajo

• Evita la resolución de LMI a gran escala
• Proporciona garantías teóricas de convergencia
• Aplicable a sistemas de gran escala
• Posee criterios explícitos de optimalidad global

Conclusiones y Discusión

Conclusiones Principales

1. El método KLAP transforma exitosamente un problema de optimización con restricciones en un problema sin restricciones
2. Garantiza optimalidad global en el caso de términos de realimentación directa sesgosimétricos
3. Proporciona mecanismos efectivos de detección de óptimos locales y reinicio
4. Demuestra eficiencia computacional superior en múltiples pruebas de referencia

Limitaciones

1. Para términos de realimentación directa no triviales, pueden existir múltiples óptimos locales
2. Requiere que el sistema satisfaga la suposición de estabilidad asintótica
3. Actualmente se enfoca principalmente en optimización en norma $\mathcal{H}_2$

Direcciones Futuras

1. Extensión al lema de realidad acotada para encontrar el sistema contractivo más cercano
2. Aplicación a sistemas parametrizados y ecuaciones diferenciales algebraicas
3. Investigación del problema de pasivación óptima en norma $\mathcal{H}_\infty$

Evaluación Profunda

Fortalezas

1. Contribuciones Teóricas Sólidas: Proporciona análisis teórico completo incluyendo existencia, unicidad y optimalidad global
2. Innovación Metodológica Fuerte: Utiliza ingeniosamente la descomposición de bajo rango del lema KYP, evitando cuellos de botella computacionales de métodos tradicionales
3. Practicidad Destacada: El algoritmo es fácil de implementar y aplicable a sistemas de gran escala
4. Experimentación Suficiente: Verifica la efectividad del método en sistemas de referencia de diferentes escalas

Insuficiencias

1. Problema de Óptimos Locales: Para términos de realimentación directa generales, aún puede atrapase en óptimos locales
2. Dependencia de Inicialización: El desempeño del método depende en cierta medida de la calidad de la inicialización
3. Análisis Teórico Limitado: El análisis para el caso $D + D^T \not\succ 0$ no es suficientemente completo

Impacto

1. Valor Académico: Proporciona una nueva perspectiva teórica y método de solución para el problema de pasivación
2. Valor Práctico: Particularmente adecuado para necesidades de pasivación de sistemas de ingeniería de gran escala
3. Reproducibilidad: El código y datos se han publicado, facilitando la verificación y aplicación

Escenarios de Aplicación

1. Pasivación de sistemas lineales de gran escala
2. Recuperación de pasividad después de reducción de orden de modelos
3. Postprocesamiento después de identificación de sistemas basada en datos
4. Diseño de sistemas de interconexión en red

Referencias Bibliográficas

El artículo cita 58 referencias relacionadas, que abarcan principalmente:

• Fundamentos de teoría de sistemas disipativos Willems, 1972
• Teoría del lema KYP y realidad positiva Anderson & Vongpanitlerd, 1973
• Revisión de métodos de pasivación Grivet-Talocia & Gustavsen, 2016
• Métodos de optimización numérica Boyd et al., 1989