Synthesis of Limit Cycles and Reference Tracking via Switching Affine Systems

技術的概要：スイッチングアフィン系によるリミットサイクルの合成と参照追従

問題定義
周期的振る舞いは生物学から工学に至るまで広範な分野に普遍的に見られるが、非線形振動子（例えば、Kuramoto モデルや Van der Pol 振動子）の数学的解析は依然として制約されている。特に、一意性と安定性に関するリミットサイクルの特性評価は、特殊な場合を除いてしばしば困難である。機械学習などの既存のデータ駆動型手法は周期的ダイナミクスを近似できるが、形式的検証に必要な構造的性質を欠くことが一般的である。さらに、スイッチングまたは区画アフィン系（PAS）を用いたリミットサイクルの近似に関する先行研究は、主に単純な 2 領域分割を有する平面（$\mathbb{R}^2$）設定に限定されるか、あるいは高次元における大域安定性に関する保証を提供できていない。

本研究の主要な目的は以下の 2 点である：

1. 任意次元（$\mathbb{R}^{n_x}$）のスイッチングアフィン系を用いて、非線形常微分方程式（ODE）のリミットサイクルを近似する合成手法を開発し、生成されるリミットサイクルが大域安定であることを保証する。
2. 参照システム自体が大域安定性の保証を持たない場合であっても、周期的軌道への漸近収束を達成する、そのようなシステム向けの参照追従コントローラを設計する。

手法

状態空間の分割とダイナミクスの合成
提案手法は、非線形システムの安定なリミットサイクルに沿って収集されたサンプル状態点の集合 $\mathcal{F}$ から開始する。状態空間は有限個の多面体 $P_i$ に分割される。先行する平面手法とは異なり、本アプローチは共通の中心点 $x_s$ と、$x_s$ で交差する境界平面（2 次元では半直線）の集合を定義することで、任意次元（$n_x = 2, 3, >3$）を扱う。

• $n_x=2$ の場合、$x_s$ から選択されたサンプル点を通る直線を引き、領域を定義する。
• $n_x \geq 3$ の場合、まずサンプル点との距離を最小化する超平面を決定し、$x_s$ を通る交線 $\Gamma$ を定義する。その後、$\Gamma$ と特定のサンプル点を含む平面によって領域を形成する。

各領域 $P_i$ 内のダイナミクスは $\dot{x} = A_i x + b_i + Bu(t)$ としてモデル化される。合成は、制約付き数値最適化問題として定式化される。目的は、モデルの軌道とサンプルデータ点との間の誤差を最小化することである。重要なのは、最適化が大域安定性を保証するために収縮性を強制することである：

1. 連続性: ベクトル場はスイッチング境界を跨いで連続でなければならない。
2. 共通リヤプノフ関数: すべての領域に対して $A_i^T Q + Q A_i \prec 0$ となる正定行列 $Q$ が存在しなければならない。これによりシステムが収縮的となり、入力信号 $u(t)$ と同じ周期を持つ一意のリミットサイクルへの任意の軌道の収束が保証される。

コスト関数における行列指数関数の非線形性に対処するため、本手法はテイラー級数近似と（領域間の遷移時間を小さく保つ、$A_i$ の特異値を有界にするなどの）対策を併用し、過度な計算複雑性を伴わずに精度を維持する。

参照追従
追従問題において、同定された周期的参照軌道が（近似精度を向上させるために安定性制約を同定段階で省略した場合など）大域安定性を有しないシナリオに対処する。制御則は $u_c(t) = u_r(t) + K_i x_c(t) + w_i - (K_j x_r(t) + w_j)$ として設計される。
漸近追従（$\lim_{t\to\infty} ||x_c(t) - x_r(t)|| = 0$）を達成するため、著者は共通の単一リヤプノフ関数ではなく、複数のリヤプノフ関数を用いる合成手法を提案する。

• 参照軌道を含む各領域に対して、区画二次リヤプノフ関数 $V_i$ を定義する。
• 各領域での滞在時間が最小閾値 $T_{min}$ を超えることを前提として、リヤプノフ関数が領域内で厳密に減少し、スイッチング境界を跨いで有界に保たれることを保証する安定性条件を導出する。
• このアプローチは共通行列 $Q$ の要件を緩和し、単一リヤプノフ関数に依存する手法と比較して、より保守的ではない収束結果をもたらす。

主要な貢献

1. 高次元合成: 本論文は、リミットサイクル近似を、平面の 2 領域分割から、一般の分割を有する高次元状態空間（$\mathbb{R}^{n_x}$）へ拡張するスイッチングアフィン系による手法を提示する。
2. 大域安定性の保証: 収縮原理を活用し、合成最適化中に共通リヤプノフ関数を強制することで、生成されるリミットサイクルが大域安定かつ一意であることを保証する。
3. 複数のリヤプノフ関数による参照追従: 周期的スイッチングアフィン系向けの追従コントローラを導入し、複数のリヤプノフ関数を利用する。これは、この文脈における周期的参照追従へのそのような手法の初適用として提示され、共通リヤプノフ関数アプローチに対するより保守的ではない代替案を提供する。
4. 構築的分割規則: サンプリングデータに基づき、2 次元、3 次元、および高次元における状態空間の分割のための具体的な幾何学的規則を提供し、生成される多面体が合成制約と互換性があることを保証する。

結果
提案手法の有効性は、数値例を通じて実証された：

• 2 次元および 3 次元近似: 合成手法は、2 次元および 3 次元空間においてサンプルデータからリミットサイクルの再構成に成功した。シミュレーションにより、リミットサイクルの内外を含む様々な初期条件から始まる軌道が、一意の大域安定リミットサイクルに収束することが示された。行列指数関数のテイラー近似の次数を上げることで、近似誤差が大幅に減少した。
• 参照追従: 2 次元の数値例は、追従コントローラが周期的参照軌道を漸近的に追従する能力を実証した。参照システムは明示的に大域安定性の保証を持たないよう（共通リヤプノフ条件に違反するように）構築されたが、制御されたシステムは遠隔の初期状態から参照を成功裏に追従し、追従誤差はゼロに収束した。

意義と主張
著者は、この研究が、特に安定性と同期の理解が不可欠である生物学的概日リズムなどの応用において、振動現象の体系的な調査における重要なギャップを埋めると主張する。本手法は、本質的な動的性質（安定性、一意性）を保持しつつ厳密な検証を可能にする扱いやすいモデルクラス（スイッチングアフィン系）を提供することで、複雑な非線形振動子の解析を容易にする。

本論文は控えめに、最適化問題は非凸であり、大域解が事前に保証されるわけではないものの、実行可能解はすべて収束する軌道をもたらすと指摘している。さらに、著者は共通リヤプノフ関数の要件が近似精度に対して制限的となり得ることを認めており、したがって、データ適合を最大化するために厳密な安定性制約なしにモデルを同定し、その後、複数のリヤプノフ関数アプローチを用いて追従コントローラを合成して収束を確保する、という 2 段階戦略を提案している。今後の課題として、共通中心点を有さない分割スキームの探求と、複数の振動子の結合が挙げられている。