Zonotope-Based Elastic Tube Model Predictive Control

Resumen Técnico: Control Predictivo Basado en Modelos de Tubo Elástico con Zonótopos

Planteamiento del Problema
El Control Predictivo Basado en Modelos (MPC) de tipo tubo es un marco robusto para sistemas lineales con restricciones sujetos a perturbaciones aditivas, garantizando que todas las trayectorias futuras permanezcan dentro de un "tubo" de conjuntos. Aunque es efectivo, la tratabilidad numérica del MPC basado en tubos suele verse limitada por la complejidad de las parametrizaciones de conjuntos. Los enfoques poliedrales tradicionales sufren un crecimiento cuadrático en las variables de decisión en relación con la longitud del horizonte o requieren descripciones complejas que aumentan el tamaño del problema. Los tubos homotéticos (escalados) y elásticos (controlados por desplazamiento) ofrecen dominios de atracción ampliados, pero a menudo introducen bilinearidades o requieren la preselección de variables auxiliares para manejar las inclusiones de conjuntos, lo que conduce a una selección de parámetros subóptima o a cargas computacionales elevadas, particularmente en sistemas de dimensiones superiores.

Metodología
El artículo propone un marco novedoso de MPC de Tubo Elástico que utiliza zonótopos escalados para parametrizar las secciones transversales del tubo. La metodología central incluye:

1. Representación de Zonótopos Escalados: El tubo se define como una secuencia de zonótopos $\langle c_k, G\Delta_k \rangle$, donde $c_k$ es el centro, $G$ es una matriz generadora fija y $\Delta_k = \text{diag}(\delta_k)$ representa factores de escala variables en el tiempo. Esto permite que el tubo se estire y contraiga elásticamente a lo largo de direcciones específicas.
2. Condiciones de Inclusión Lineales: Se deriva una nueva clase de condiciones de inclusión para zonótopos escalados. A diferencia de métodos anteriores que dependen de términos bilineales o requieren la preselección de variables auxiliares (como $\Gamma$) para manejar la escala, el enfoque propuesto utiliza una matriz $\Phi$ para codificar linealmente la contención de sumas de Minkowski de zonótopos.
   • Teorema 1 formula la optimización del MPC como un Programa Cuadrático (QP) donde las restricciones de alcanzabilidad en un paso y de admisibilidad de estado/entrada se expresan como igualdades y desigualdades lineales que involucran los factores de escala $\delta_k$ y la matriz auxiliar $\Phi_k$.
   • Crucialmente, los centros de los zonótopos $c_k$ se tratan como variables de decisión independientes, lo que mejora la flexibilidad sin introducir no convexidades.
3. Alternativas de Precomputación: Para reducir aún más la complejidad en tiempo de ejecución, el artículo introduce dos variantes:
   • Basada en $\Phi_0$: Una matriz constante $\Phi_0$ se precomputa fuera de línea, reemplazando la optimización de $\Phi_k$ en cada paso con una restricción lineal sobre los factores de escala.
   • Basada en $\Gamma$: Se utiliza un par precomputado $(\Gamma, \gamma)$, similar a los enfoques poliedrales anteriores pero adaptado para zonótopos.
4. Conjuntos Invariantes Positivamente Robustos (RPI): La matriz generadora $G$ se deriva de una aproximación de conjunto RPI para la dinámica del error. El artículo mejora esto tratando el centro del conjunto RPI como una variable de decisión en el problema de minimización, reduciendo el número de variables de decisión en comparación con trabajos anteriores.
5. Estabilidad y Factibilidad: La factibilidad recursiva y la estabilidad asintótica se garantizan construyendo un conjunto terminal mediante recurrencia de conjuntos y haciendo cumplir condiciones de contracción tipo Lyapunov sobre el costo terminal y los factores de escala.

Contribuciones Clave

• Codificación Lineal Novel: El artículo introduce una codificación lineal para la contención de zonótopos (Corolario 1) que generaliza las condiciones estándar. Esto elimina la necesidad de preseleccionar variables auxiliares para evitar la optimización no lineal, un problema común en la escala poliedral.
• Eficiencia Computacional: La formulación propuesta permanece lineal tanto en variables auxiliares como en factores de escala, facilitando su integración en marcos de optimización a gran escala. Evita las bilinearidades típicamente asociadas con los factores de escala en otros métodos de MPC de tubo.
• Reducción de Complejidad: Al aprovechar los zonótopos, el método reduce significativamente el número de restricciones y factores de escala requeridos en comparación con los tubos elásticos poliedrales, especialmente a medida que aumenta la dimensión del sistema.
• Análisis Explícito de Compensaciones: El artículo proporciona un análisis exhaustivo de la complejidad que compara tubos rígidos, homotéticos y elásticos en configuraciones poliedrales y zonotópicas, ilustrando la compensación entre el dominio de atracción ampliado y el esfuerzo computacional.

Resultados
El enfoque propuesto fue validado mediante experimentos numéricos en dos puntos de referencia:

1. Doble Integrador: Las comparaciones entre tubos elásticos poliedrales, tubos elásticos zonotópicos con $\Phi_0$ precomputado y tubos elásticos zonotópicos con $\Phi$ optimizado mostraron que las variantes zonotópicas producen dominios de atracción más grandes. La variante con $\Phi$ totalmente optimizado ofreció una ligera ventaja sobre la variante con $\Phi_0$ precomputado, pero a un costo computacional mayor.
2. Experimento de Resorte Acoplado (CSE): Este sistema de alta dimensión (hasta 20 estados) demostró la escalabilidad del método. Las implementaciones poliedrales fallaron al encontrar soluciones factibles para sistemas con 3 o más resortes debido a límites computacionales. En contraste, todas las variantes zonotópicas resolvieron exitosamente el problema en todo el rango de parámetros.
   • Tiempo de Ejecución: Los métodos zonotópicos exhibieron complejidad computacional lineal (o cuadrática para la variante $\Phi$), mientras que los métodos poliedrales se volvieron intratables.
   • Centros: Tratar los centros de los zonótopos como variables de decisión resultó en una sobrecarga computacional solo marginal mientras proporcionaba garantías de robustez cruciales.

Significado
El artículo afirma que integrar zonótopos escalados elásticamente en el MPC de tubo reduce significativamente el tiempo de ejecución y mejora la escalabilidad para sistemas de alta dimensión en comparación con los métodos poliedrales tradicionales. Al aliviar la necesidad de especificación a priori de restricciones de contención de conjuntos y evitar bilinearidades, el enfoque logra un equilibrio más eficiente entre el dominio de atracción y el esfuerzo computacional. El trabajo extiende la investigación preliminar sobre MPC de tubo elástico proporcionando una alternativa basada en precomputación computacionalmente eficiente y estableciendo condiciones rigurosas para la factibilidad recursiva y la estabilidad. Los autores señalan que, aunque el método es robusto, el trabajo futuro podría abordar la sensibilidad al conjunto RPI inicial y posibles extensiones al MPC basado en datos.