Synthesis of Limit Cycles and Reference Tracking via Switching Affine Systems

Resumo Técnico: Síntese de Ciclos Limites e Rastreamento de Referência via Sistemas Afins comutados

Declaração do Problema
O comportamento periódico é ubíquo em áreas que vão da biologia à engenharia, contudo a análise matemática de osciladores não lineares (por exemplo, Kuramoto, Van der Pol) permanece restrita. A caracterização de ciclos limites, particularmente no que tange à unicidade e estabilidade, é frequentemente viável apenas em casos especiais. Técnicas existentes baseadas em dados, como aprendizado de máquina, podem aproximar dinâmicas periódicas, mas geralmente carecem das propriedades estruturais necessárias para verificação formal. Além disso, abordagens anteriores para aproximar ciclos limites usando sistemas comutados ou por partes afins (PAS) foram amplamente confinadas a configurações planares ($\mathbb{R}^2$) com partições simples de duas regiões, ou falham em fornecer garantias de estabilidade global em dimensões superiores.

O objetivo principal deste trabalho é duplo:

1. Desenvolver um método de síntese que aproxime os ciclos limites de equações diferenciais ordinárias (EDOs) não lineares usando sistemas afins comutados em dimensões arbitrárias ($\mathbb{R}^{n_x}$), garantindo que o ciclo limite resultante seja globalmente estável.
2. Projetar um controlador de rastreamento de referência para tais sistemas que alcance convergência assintótica a uma trajetória periódica, mesmo quando o próprio sistema de referência carece de garantias de estabilidade global.

Metodologia

Particionamento do Espaço de Estados e Síntese de Dinâmica
O método proposto inicia-se com um conjunto de pontos de estado amostrados $\mathcal{F}$ coletados ao longo de um ciclo limite estável de um sistema não linear. O espaço de estados é particionado em um número finito de politopos $P_i$. Diferentemente de métodos planares anteriores, esta abordagem lida com dimensões arbitrárias ($n_x = 2, 3, >3$) definindo um ponto central comum $x_s$ e um conjunto de planos de fronteira (ou raios em 2D) que se intersectam em $x_s$.

• Para $n_x=2$, linhas são traçadas de $x_s$ através de pontos de amostra selecionados para definir regiões.
• Para $n_x \geq 3$, um hiperplano é primeiro determinado para minimizar a distância aos pontos de amostra, definindo uma linha de interseção $\Gamma$ passando por $x_s$. As regiões são então formadas por planos contendo $\Gamma$ e pontos de amostra específicos.

A dinâmica dentro de cada região $P_i$ é modelada como $\dot{x} = A_i x + b_i + Bu(t)$. A síntese é formulada como um problema de otimização numérica com restrições. O objetivo é minimizar o erro entre a trajetória do modelo e os pontos de dados amostrados. Crucialmente, a otimização impõe contratividade para garantir estabilidade global:

1. Continuidade: Os campos vetoriais devem ser contínuos através das fronteiras de comutação.
2. Função de Lyapunov Comum: Deve existir uma matriz definida positiva $Q$ tal que $A_i^T Q + Q A_i \prec 0$ para todas as regiões. Isso garante que o sistema seja contrativo, assegurando que qualquer trajetória convirja para um ciclo limite único com o mesmo período do sinal de entrada $u(t)$.

Para lidar com a não linearidade da exponencial de matriz na função de custo, o método emprega aproximações de séries de Taylor e contramedidas (como garantir tempos de transição pequenos entre regiões e limitar os valores singulares de $A_i$) para manter a precisão sem complexidade computacional excessiva.

Rastreamento de Referência
Para o problema de rastreamento, o artigo aborda cenários onde a trajetória de referência periódica identificada pode não ser globalmente estável (por exemplo, se restrições de estabilidade foram omitidas durante a identificação para melhorar a precisão da aproximação). A lei de controle é projetada como $u_c(t) = u_r(t) + K_i x_c(t) + w_i - (K_j x_r(t) + w_j)$.
Para alcançar o rastreamento assintótico ($\lim_{t\to\infty} ||x_c(t) - x_r(t)|| = 0$), os autores propõem um esquema de síntese usando múltiplas funções de Lyapunov em vez de uma única comum.

• Uma função de Lyapunov quadrática por partes $V_i$ é definida para cada região contendo a trajetória de referência.
• Condições de estabilidade são derivadas para garantir que a função de Lyapunov diminua estritamente dentro das regiões e permaneça limitada através das fronteiras de comutação, desde que o tempo de permanência em cada região exceda um limiar mínimo $T_{min}$.
• Esta abordagem relaxa o requisito de uma matriz comum $Q$, permitindo resultados de convergência menos conservadores em comparação com métodos que dependem de uma única função de Lyapunov.

Contribuições Principais

1. Síntese em Alta Dimensão: O artigo estende a aproximação de ciclos limites via sistemas afins comutados de partições planares de duas regiões para espaços de estados de dimensões superiores ($\mathbb{R}^{n_x}$) com partições gerais.
2. Garantia de Estabilidade Global: Ao alavancar o princípio da contração e impor uma função de Lyapunov comum durante a otimização de síntese, o método garante que o ciclo limite gerado seja globalmente estável e único.
3. Rastreamento de Referência com Múltiplas Funções de Lyapunov: O trabalho introduz um controlador de rastreamento para sistemas afins comutados periódicos que utiliza múltiplas funções de Lyapunov. Isso é apresentado como a primeira aplicação de tal técnica ao rastreamento de referência periódica neste contexto, oferecendo uma alternativa menos conservadora às abordagens de Lyapunov comum.
4. Regras de Particionamento Construtivas: O artigo fornece regras geométricas específicas para particionar espaços de estados em 2D, 3D e dimensões superiores com base em dados amostrados, garantindo que os politopos resultantes sejam compatíveis com as restrições de síntese.

Resultados
A eficácia dos métodos propostos é demonstrada através de exemplos numéricos:

• Aproximação em 2D e 3D: O método de síntese reconstruiu com sucesso ciclos limites a partir de dados amostrados em espaços tanto 2D quanto 3D. Simulações mostraram que trajetórias iniciadas a partir de várias condições iniciais (tanto dentro quanto fora do ciclo limite) convergiram para o ciclo limite único e globalmente estável. O aumento da ordem da aproximação de Taylor para a exponencial de matriz reduziu significativamente o erro de aproximação.
• Rastreamento de Referência: Um exemplo numérico em 2D demonstrou a capacidade do controlador de rastreamento de seguir assintoticamente uma trajetória de referência periódica. O sistema de referência foi construído explicitamente para carecer de garantias de estabilidade global (violando a condição de Lyapunov comum), contudo o sistema controlado rastreou com sucesso a referência a partir de um estado inicial distante, com o erro de rastreamento convergindo para zero.

Significado e Alegações
Os autores afirmam que este trabalho aborda uma lacuna crítica na investigação sistemática de fenômenos oscilatórios, particularmente em aplicações como ritmos circadianos biológicos, onde a compreensão da estabilidade e sincronização é indispensável. Ao fornecer uma classe de modelos tratável (sistemas afins comutados) que preserva propriedades dinâmicas essenciais (estabilidade, unicidade) enquanto permite verificação rigorosa, o método facilita a análise de osciladores não lineares complexos.

O artigo nota modestamente que, embora o problema de otimização seja não convexo e soluções globais não sejam garantidas a priori, qualquer solução viável produz uma trajetória convergente. Além disso, os autores reconhecem que o requisito de uma função de Lyapunov comum pode ser restritivo para a precisão da aproximação; assim, propõem uma estratégia de dois passos: identificar o modelo sem restrições estritas de estabilidade para maximizar o ajuste aos dados e, em seguida, sintetizar um controlador de rastreamento usando a abordagem de múltiplas funções de Lyapunov para garantir a convergência. Trabalhos futuros são identificados como a exploração de esquemas de particionamento sem um ponto central comum e o acoplamento de múltiplos osciladores.