Synthesis of Limit Cycles and Reference Tracking via Switching Affine Systems

기술 요약: 스위칭 아핀 시스템을 통한 리미트 사이클 합성 및 참조 추종

문제 제기
주기적 행동은 생물학부터 공학에 이르기까지 다양한 분야에서 보편적으로 나타나지만, 비선형 발진기 (예: 쿠라모토, 반 더 폴) 의 수학적 분석은 여전히 제한적입니다. 특히 유일성과 안정성과 관련된 리미트 사이클의 특성 규명은 종종 특수한 경우에만 가능합니다. 머신러닝과 같은 기존 데이터 기반 기법은 주기적 동역학을 근사할 수 있지만, 형식적 검증을 위해 필요한 구조적 속성을 일반적으로 갖추지 못합니다. 또한, 스위칭 또는 조각별 아핀 시스템 (PAS) 을 사용하여 리미트 사이클을 근사하는 이전 접근법들은 주로 단순한 두 영역 분할을 가진 평면 ($\mathbb{R}^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)$로 모델링됩니다. 합성은 제약 조건이 있는 수치 최적화 문제로 공식화됩니다. 목적 함수는 모델의 궤적과 샘플 데이터 점들 사이의 오차를 최소화하는 것입니다. 핵심적으로, 이 최적화는 전역 안정성을 보장하기 위해 **수축성 (contractivity)**을 강제합니다:

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$) 을 달성하기 위해, 저자들은 단일 공통 리아푸노프 함수 대신 다중 리아푸노프 함수를 사용하는 합성 계획을 제안합니다.

• 참조 궤적을 포함하는 각 영역에 대해 조각별 2 차 리아푸노프 함수 $V_i$가 정의됩니다.
• 각 영역 내에서는 리아푸노프 함수가 엄격히 감소하고 스위칭 경계를 가로질러는 유계로 유지되도록 하는 안정성 조건이 유도되며, 이는 각 영역에서의 체류 시간이 최소 임계값 $T_{min}$을 초과할 때 보장됩니다.
• 이 접근법은 공통 행렬 $Q$에 대한 요구 사항을 완화하여, 단일 리아푸노프 함수에 의존하는 방법들에 비해 덜 보수적인 수렴 결과를 제공합니다.

주요 기여

1. 고차원 합성: 이 논문은 평면, 두 영역 분할을 통한 리미트 사이클 근사를 일반 분할을 가진 고차원 상태 공간 ($\mathbb{R}^{n_x}$) 으로 확장합니다.
2. 전역 안정성 보장: 수축 원리를 활용하고 합성 최적화 과정에서 공통 리아푸노프 함수를 강제함으로써, 생성된 리미트 사이클이 전역적으로 안정적이고 유일함을 보장합니다.
3. 다중 리아푸노프 함수를 통한 참조 추종: 이 연구는 다중 리아푸노프 함수를 활용하는 주기적 스위칭 아핀 시스템을 위한 추종 제어기를 도입합니다. 이는 이러한 맥락에서 주기적 참조 추종에 이러한 기법을 적용한 첫 번째 사례로 제시되며, 공통 리아푸노프 접근법에 비해 덜 보수적인 대안을 제공합니다.
4. 구성적 분할 규칙: 이 논문은 샘플 데이터에 기반하여 2 차원, 3 차원 및 고차원에서 상태 공간을 분할하기 위한 구체적인 기하학적 규칙을 제공하여, 결과적으로 생성된 다면체가 합성 제약 조건과 호환되도록 보장합니다.

결과
제안된 방법의 효과는 수치 예시를 통해 입증되었습니다:

• 2 차원 및 3 차원 근사: 합성 방법은 2 차원 및 3 차원 공간에서 샘플 데이터로부터 리미트 사이클을 성공적으로 재구성했습니다. 시뮬레이션은 리미트 사이클 내부와 외부 모두의 다양한 초기 조건에서 시작하는 궤적들이 유일하고 전역적으로 안정적인 리미트 사이클로 수렴함을 보여주었습니다. 행렬 지수함수에 대한 테일러 근사의 차수를 높이면 근사 오차가 크게 감소했습니다.
• 참조 추종: 2 차원의 수치 예시는 추종 제어기가 주기적 참조 궤적을 점근적으로 추종할 수 있음을 보여주었습니다. 참조 시스템은 전역 안정성 보장 (공통 리아푸노프 조건 위반) 이 결여되도록 명시적으로 구성되었음에도 불구하고, 제어된 시스템은 먼 초기 상태에서도 참조를 성공적으로 추종하여 추종 오차가 0 으로 수렴했습니다.

의의 및 주장
저자들은 이 연구가 생물학적 일주기 리듬과 같은 안정성과 동기화를 이해하는 것이 필수적인 응용 분야에서 진동 현상에 대한 체계적 조사에서 중요한 격차를 해소했다고 주장합니다. 본질적인 동적 속성 (안정성, 유일성) 을 보존하면서도 엄격한 검증을 허용하는 처리 가능한 모델 클래스 (스위칭 아핀 시스템) 를 제공함으로써, 이 방법은 복잡한 비선형 발진기의 분석을 용이하게 합니다.

이 논문은 최적화 문제가 비볼록하며 전역 해가 사전에 보장되지 않지만, 어떤 실현 가능한 해라도 수렴하는 궤적을 산출한다는 점을 겸손하게 지적합니다. 또한, 저자들은 공통 리아푸노프 함수에 대한 요구 사항이 근사 정확도에 제한적일 수 있음을 인정하며, 따라서 데이터 적합을 극대화하기 위해 엄격한 안정성 제약 없이 모델을 식별한 후, 수렴을 보장하기 위해 다중 리아푸노프 함수 접근법을 사용하여 추종 제어기를 합성하는 2 단계 전략을 제안합니다. 향후 작업은 공통 중심점 없이 분할 체계를 탐구하고 여러 발진기를 결합하는 것으로 식별되었습니다.