A piecewise ellipsoidal reachable set estimation method for continuous bimodal piecewise affine systems

Resumen Técnico: Un Método de Estimación de Conjuntos Alcanzables Elipsoidales por Partes para Sistemas Afines por Partes Bimodales Continuos

Planteamiento del Problema
El artículo aborda el problema de estimar los conjuntos alcanzables para sistemas afines por partes (PWA) bimodales continuos sujetos a perturbaciones exógenas. Aunque los sistemas PWA se utilizan ampliamente para modelar fenómenos físicos (p. ej., relés, histéresis, saturación) y aproximar la dinámica no lineal, calcular sus conjuntos alcanzables es, por lo general, un problema indecidible. Los métodos existentes para sistemas híbridos PWA generales suelen sufrir de una alta complejidad computacional o no son adecuados para la estructura dinámica específica de los sistemas bimodales. Además, aunque existen técnicas de acotación elipsoidal para sistemas lineales y conmutados, estas a menudo dependen de funciones de Lyapunov comunes o requieren descomposiciones matriciales específicas que pueden conducir a estimaciones conservadoras o fallar al manejar la conmutación dependiente del estado con perturbaciones no evanescentes de manera efectiva.

Metodología
Los autores proponen un nuevo enfoque dentro del marco de la acotación elipsoidal, utilizando funciones de Lyapunov cuadráticas por partes adaptadas a la naturaleza bimodal del sistema. El sistema se define por dos dinámicas afines que conmutan basándose en un hiperplano $c^T x + f = 0$. Las perturbaciones están contenidas dentro de un elipsoide.

La metodología procede a través de los siguientes pasos técnicos:

1. Condiciones de Positividad de Semiespacio: Los autores derivan primero condiciones (Lema 3.1) para asegurar que una función cuadrática sea positiva en un semiespacio específico. Esto se logra caracterizando el semiespacio mediante una transformación matricial específica que involucra los parámetros de la superficie de conmutación.
2. Restricciones de Continuidad: Para asegurar que la función de Lyapunov por partes sea continua a través del hiperplano de conmutación, los autores establecen condiciones (Lema 3.2) que relacionan los parámetros de las funciones cuadráticas en los dos modos.
3. Homogeneización y Procedimiento S: Los autores emplean el concepto de homogeneización para conjuntos y un procedimiento S generalizado (Lema 3.5) para manejar simultáneamente las restricciones impuestas por el elipsoide de perturbación y la superficie de conmutación.
4. Formulación de LMI: Al combinar estos lemas con una caracterización de Lyapunov del conjunto alcanzable (Lema 2.3), los autores derivan condiciones de Desigualdades Matriciales Lineales (LMI) (Teorema 3.6 y Corolario 3.7). Estas condiciones garantizan que la derivada de la función de Lyapunov por partes sea definida negativa fuera de una región específica, delimitando así el conjunto alcanzable dentro de un conjunto elipsoidal por partes.

Contribuciones Clave

• Nuevas Caracterizaciones LMI: El artículo introduce nuevas condiciones LMI (Lemas 3.1, 3.2, 3.4 y Teorema 3.5) diseñadas específicamente para explotar la dinámica bimodal. A diferencia de los métodos anteriores que asumen una matriz de Lyapunov común o requieren descomposiciones especiales (p. ej., $P_i = F_i^T T F_i$), este enfoque construye el elipsoide por partes sin tales suposiciones restrictivas.
• Manejo de Perturbaciones No Evanescentes: El método acomoda perturbaciones limitadas y no evanescentes, una característica donde la literatura existente (p. ej., [28]) requiere que las perturbaciones se vuelvan nulas en el infinito.
• Conmutación Dependiente del Estado: El enfoque está formulado para la conmutación dependiente del estado, distinguiéndose de los métodos que dependen de la conmutación dependiente del tiempo.
• Eficiencia Computacional: El método propuesto resulta en un problema de optimización convexo que puede resolverse eficientemente, ofreciando una reducción significativa en el tiempo de computación en comparación con los enfoques basados en zonotopos para sistemas de orden superior.

Resultos
El método propuesto se valida mediante dos ejemplos numéricos:

1. Ejemplo 4.1 (Sistema 2D): Se analiza un sistema bimodal con planos de fase "ortogonales". Los resultados demuestran que la función de Lyapunov cuadrática por partes proporciona un límite más ajustado (menos conservador) para el conjunto alcanzable en comparación con una función de Lyapunov cuadrática común. El método por partes captura con éxito la geometría del sistema que el método común pasa por alto.
2. Ejemplo 4.2 (Sistema Mecánico 4D): Se utiliza un sistema mecánico con un resorte de un solo lado para comparar el método propuesto contra la herramienta CORA (que utiliza zonotopos).
   • Precisión: Se muestra que la estimación elipsoidal por partes es más pequeña (más ajustada) que la estimación elipsoidal común.
   • Desempeño: El método propuesto resuelve el sistema de 4º orden en aproximadamente 0.1 segundos. En contraste, la herramienta CORA requirió más de 200 segundos para un zonotopo de orden 40 y más de 215 segundos para un zonotopo de orden 50. Notablemente, aumentar el orden del zonotopo en CORA incrementó el tiempo de computación sin mejorar significativamente la estimación, mientras que el método propuesto logró una alta precisión con un costo computacional mínimo.

Significancia y Reclamaciones
El artículo sostiene que su principal significancia radica en explotar la estructura dinámica específica de los sistemas PWA bimodales continuos para proponer un enfoque más simple y menos conservador que los métodos generales de sistemas híbridos. Los autores enfatizan que sus condiciones son novedosas porque evitan la necesidad de descomposiciones especiales de las matrices de Lyapunov, las cuales suelen ser difíciles de derivar de manera óptima en la literatura existente.

Los autores señalan modestamente que, si bien el método es efectivo para sistemas bimodales continuos, extenderlo a sistemas PWA multimodales generales (trimodales o superiores) no es trivial y requiere investigación futura, involucrando potencialmente representaciones especiales como hiperplanos de articulación (hinging hyperplanes). Este trabajo se posiciona como una herramienta especializada y eficiente para una clase de sistemas que incluye importantes aplicaciones de ingeniería, como los sistemas de relé lineales.