Synthesis of Limit Cycles and Reference Tracking via Switching Affine Systems

Technische Samenvatting: Synthese van Limietcycli en Referentievolgeling via Schakelende Affiene Systemen

Probleemstelling
Periodiek gedrag is alomtegenwoordig in gebieden variërend van biologie tot techniek, doch de wiskundige analyse van niet-lineaire oscillatoren (bijv. Kuramoto, Van der Pol) blijft beperkt. Het karakteriseren van limietcycli, met name wat betreft uniciteit en stabiliteit, is vaak alleen haalbaar in speciale gevallen. Bestaande datagedreven technieken, zoals machine learning, kunnen periodieke dynamica benaderen maar missen doorgaans de structurele eigenschappen die vereist zijn voor formele verificatie. Bovendien waren eerdere benaderingen voor het benaderen van limietcycli met behulp van schakelende of stuksgewijs-affiene systemen (PAS) grotendeels beperkt tot planaire ($\mathbb{R}^2$) omgevingen met eenvoudige partities in twee gebieden, of ze leverden geen garanties voor globale stabiliteit in hogere dimensies.

Het primaire doel van dit werk is tweeledig:

1. Een synthesemethode ontwikkelen die de limietcycli van niet-lineaire gewone differentiaalvergelijkingen (ODE's) benadert met behulp van schakelende affiene systemen in willekeurige dimensies ($\mathbb{R}^{n_x}$), waarbij wordt gegarandeerd dat de resulterende limietcyclus globaal stabiel is.
2. Een referentievolgeling ontwerpen voor dergelijke systemen die asymptotische convergentie naar een periodieke trajectorie bereikt, zelfs wanneer het referentiesysteem zelf geen garanties voor globale stabiliteit biedt.

Methodologie

Partitie van de Toestandsruimte en Synthese van Dynamica
De voorgestelde methode begint met een set van bemonsterde toestandspunten $\mathcal{F}$ die zijn verzameld langs een stabiele limietcyclus van een niet-lineair systeem. De toestandsruimte wordt gepartitioneerd in een eindig aantal polytopen $P_i$. In tegenstelling tot eerdere planaire methoden, behandelt deze aanpak willekeurige dimensies ($n_x = 2, 3, >3$) door een gemeenschappelijk middelpunt $x_s$ en een set van grensvlakken (of stralen in 2D) te definiëren die samenkomen in $x_s$.

• Voor $n_x=2$ worden lijnen getrokken vanuit $x_s$ door geselecteerde steekproefpunten om gebieden te definiëren.
• Voor $n_x \geq 3$ wordt eerst een hypervlak bepaald om de afstand tot de steekproefpunten te minimaliseren, waardoor een snijlijn $\Gamma$ wordt gedefinieerd die door $x_s$ loopt. Gebieden worden vervolgens gevormd door vlakken die $\Gamma$ en specifieke steekproefpunten bevatten.

De dynamica binnen elk gebied $P_i$ wordt gemodelleerd als $\dot{x} = A_i x + b_i + Bu(t)$. De synthese wordt geformuleerd als een beperkt numeriek optimalisatieprobleem. Het doel is om de fout tussen de trajectorie van het model en de bemonsterde datapunten te minimaliseren. Cruciaal is dat de optimalisatie contractiviteit afdwingt om globale stabiliteit te garanderen:

1. Continuïteit: De vectorvelden moeten continu zijn over de schakelgrenzen.
2. Gemeenschappelijke Lyapunov-functie: Er moet een positief-definiete matrix $Q$ bestaan zodat $A_i^T Q + Q A_i \prec 0$ voor alle gebieden. Dit zorgt ervoor dat het systeem contractief is, wat garandeert dat elke trajectorie convergeert naar een unieke limietcyclus met dezelfde periode als het ingangssignaal $u(t)$.

Om de niet-lineariteit van de matrixexponentiële in de kostenfunctie te hanteren, maakt de methode gebruik van Taylorreeksbenaderingen en tegenmaatregelen (zoals het waarborgen van korte overgangstijden tussen gebieden en het begrenzen van de singuliere waarden van $A_i$) om nauwkeurigheid te behouden zonder excessive computationele complexiteit.

Referentievolgeling
Voor het volgprobleem adresseert het artikel scenario's waarin de geïdentificeerde periodieke referentietrajectorie mogelijk niet globaal stabiel is (bijv. als stabiliteitsbeperkingen tijdens de identificatie zijn weggelaten om de benaderingsnauwkeurigheid te verbeteren). De regellaw wordt ontworpen als $u_c(t) = u_r(t) + K_i x_c(t) + w_i - (K_j x_r(t) + w_j)$.
Om asymptotische volgeling te bereiken ($\lim_{t\to\infty} ||x_c(t) - x_r(t)|| = 0$), stellen de auteurs een syntheseschema voor met meerdere Lyapunov-functies in plaats van één gemeenschappelijke.

• Een stuksgewijs kwadratische Lyapunov-functie $V_i$ wordt gedefinieerd voor elk gebied dat de referentietrajectorie bevat.
• Stabiliteitsvoorwaarden worden afgeleid om te waarborgen dat de Lyapunov-functie strikt daalt binnen gebieden en begrensd blijft over schakelgrenzen, mits de verblijftijd in elk gebied een minimumdrempel $T_{min}$ overschrijdt.
• Deze aanpak verlicht de vereiste voor een gemeenschappelijke matrix $Q$, waardoor minder conservatieve convergentieresultaten mogelijk zijn in vergelijking met methoden die vertrouwen op een enkele Lyapunov-functie.

Belangrijkste Bijdragen

1. Synthese in Hoge Dimensies: Het artikel breidt de benadering van limietcycli via schakelende affiene systemen uit van planaire, partities in twee gebieden naar toestandsruimten in hogere dimensies ($\mathbb{R}^{n_x}$) met algemene partities.
2. Garantie voor Globale Stabiliteit: Door gebruik te maken van het contractieprincipe en het afdwingen van een gemeenschappelijke Lyapunov-functie tijdens de synthese-optimalisatie, garandeert de methode dat de gegenereerde limietcyclus globaal stabiel en uniek is.
3. Referentievolgeling met Meerdere Lyapunov-functies: Het werk introduceert een volgelingregelaar voor periodieke schakelende affiene systemen die meerdere Lyapunov-functies gebruikt. Dit wordt gepresenteerd als de eerste toepassing van een dergelijke techniek voor periodieke referentievolgeling in deze context, en biedt een minder conservatief alternatief voor benaderingen met een gemeenschappelijke Lyapunov-functie.
4. Constructieve Partitieregels: Het artikel biedt specifieke geometrische regels voor het partitioneren van toestandsruimten in 2D, 3D en hogere dimensies op basis van bemonsterde data, waarbij wordt gewaarborgd dat de resulterende polytopen compatibel zijn met de synthesebeperkingen.

Resultaten
De effectiviteit van de voorgestelde methoden wordt aangetoond door middel van numerieke voorbeelden:

• 2D en 3D Benadering: De synthesemethode reconstrueerde met succes limietcycli uit bemonsterde data in zowel 2D- als 3D-ruimten. Simulaties toonden aan dat trajecten die starten vanuit verschillende beginvoorwaarden (zowel binnen als buiten de limietcyclus) convergeerden naar de unieke, globaal stabiele limietcyclus. Het verhogen van de orde van de Taylor-benadering voor de matrixexponentiële verminderde de benaderingsfout aanzienlijk.
• Referentievolgeling: Een numeriek voorbeeld in 2D demonstreerde het vermogen van de volgelingregelaar om asymptotisch een periodieke referentietrajectorie te volgen. Het referentiesysteem werd expliciet geconstrueerd zonder garanties voor globale stabiliteit (in strijd met de voorwaarde voor een gemeenschappelijke Lyapunov-functie), toch slaagde het geregelde systeem erin om de referentie te volgen vanaf een ver verwijderde beginstaat, waarbij de volgelingfout convergeerde naar nul.

Betekenis en Beweringen
De auteurs beweren dat dit werk een kritieke kloof aanpakt in de systematische onderzoeking van oscillatoire fenomenen, met name in toepassingen zoals biologische circadiane ritmen waar het begrijpen van stabiliteit en synchronisatie onmisbaar is. Door een hanteerbare modelklasse (schakelende affiene systemen) te bieden die essentiële dynamische eigenschappen (stabiliteit, uniciteit) behoudt terwijl het rigoureuze verificatie toelaat, faciliteert de methode de analyse van complexe niet-lineaire oscillatoren.

Het artikel noteert bescheiden dat hoewel het optimalisatieprobleem niet-convex is en globale oplossingen niet a-priori gegarandeerd zijn, elke haalbare oplossing een convergente trajectorie oplevert. Bovendien erkennen de auteurs dat de vereiste voor een gemeenschappelijke Lyapunov-functie beperkend kan zijn voor de benaderingsnauwkeurigheid; daarom stellen zij een tweestapsstrategie voor: identificeer het model zonder strikte stabiliteitsbeperkingen om de datafit te maximaliseren, en synthetiseer vervolgens een volgelingregelaar met behulp van de aanpak met meerdere Lyapunov-functies om convergentie te waarborgen. Toekomstig werk wordt geïdentificeerd als het verkennen van partitieschema's zonder een gemeenschappelijk middelpunt en het koppelen van meerdere oscillatoren.