[논문 리뷰] Optimal stabilization of a cycle in nonlinear discrete systems

Dmitriy Dmitrishin, Anna Khamitova외 2명|arXiv (Cornell University)|2013. 07. 28.

Chaos control and synchronization참고 문헌 5인용 수 4

한 줄 요약

이 논문은 최소한의 이전 상태 깊이를 요구하면서도 비선형 이산 시스템에서 두 주기 사이클을 안정화시키기 위한 최적의 지연 피드백 제어(Delayed Feedback Control, DFC)를 제안한다. 제어를 과거 상태 차분의 선형 조합으로 설정하고, Fejér 다항식을 활용함으로써 필요한 메모리 깊이를 최소화하는 유일한 계수를 도출한다. 이는 (−μ*, −1) 범위 내에서 음수인 승수를 가진 사이클에 대해 점근적 안정성을 보장한다.

ABSTRACT

A method of stabilizing 2-cycles in discrete dynamic systems by Delayed Feedback Control is developed by using classic Harmonic Analysis.

연구 동기 및 목표

• 소규모이고 유계인 제어 조건 하에서 불안정한 두 주기 사이클을 안정화시키는 데 도전하는 것.
• 지연 피드백 제어(DFC)에 필요한 이전 상태 깊이를 최소화하면서도 사이클의 점근적 안정성을 보장하는 것.
• 이전의 평형점 안정화 결과를 주기적 사이클, 특히 혼돈 행동을 나타내는 음수 승수를 가진 사이클로 확장하는 것.
• 주어진 사이클 승수의 상한 μ*에 대해 제어 복잡성(즉, 최소 N*)을 최소화하는 명시적이고 최적의 제어 계수를 도출하는 것.
• 제어 설계와 극값 삼각함수 다항식 및 Fejér 커널 간의 이론적 프레임워크를 구축하는 것.

제안 방법

• 비선형 이산 맵 $x_{n+1} = f_h(x_n)$로 시스템을 모델링하고, 비안정한 두 주기 사이클 $(\eta_1, \eta_2)$과 승수 $\mu = f_h'(\eta_1)f_h'(\eta_2) \in (-\mu^*, -1)$를 고려한다.
• 지연 피드백 제어를 $u = -\sum_{j=1}^{N-1} \varepsilon_j (f_h(x_{n-2j+2}) - f_h(x_{n-2j}))$ 형태로 설계하여, 계수 $a_j$를 가진 선형 조합 $x_{n+1} = \sum_{j=1}^N a_j f_h(x_{n-2j+2})$로 축소한다.
• 사이클 주변에서 선형화를 적용하여 $\lambda$에 대한 특성 방정식을 유도하고, $\lambda - \mu p(\lambda)^2 = 0$을 도출한다. 여기서 $p(\lambda) = \sum_{j=1}^N a_j \lambda^{-(j-1)}$이다.
• 특성 방정식을 $\nu = \lambda^{1/2}$를 사용하여 삼각함수 형태로 변환하여, $\mu$와 관련된 계수를 가진 $\nu$에 대한 다항식을 도출한다.
• Fejér 다항식의 극값 성질을 활용하여 제어 스펙트럼의 최대 편차를 최소화함으로써 최적의 $a_j^0 = \frac{2(N-j)+1}{N^2}$을 도출한다.
• 최적의 $\varepsilon_j$는 $a_j^0$를 바탕으로 $\varepsilon_j = \sum_{k=j+1}^N a_k^0$를 통해 유도되며, 이는 최소한의 이전 상태 깊이 $N^* = 2(N_0 - 1)$를 보장한다.

실험 결과

연구 질문

• RQ1 비선형 이산 시스템에서 승수가 $(-\mu^*, -1)$ 범위 내에 있는 음수인 두 주기 사이클을 최적으로 안정화시키기 위해 필요한 최소한의 이전 상태 깊이는 얼마인가?
• RQ2 메모리 요구량을 최소화하면서도 사이클의 점근적 안정성을 보장하기 위해 제어 계수는 어떻게 선택해야 하는가?
• RQ3 사이클 승수가 0에서 멀리 떨어져 있고 음수일 경우, 지연 피드백 제어의 최적의 구조는 무엇인가?
• RQ4 극값 삼각함수 다항식을 사용하여 최적의 제어를 도출할 수 있으며, 이는 Fejér 커널과 어떻게 관련이 있는가?
• RQ5 제어가 사이클을 안정화하지 못하는 조건는 무엇이며, 로지스틱 맵 예제에서 $h=4$일 때 실패하는 이유는 무엇인가?

주요 결과

• 최적의 제어 계수는 $a_j^0 = \frac{2(N-j)+1}{N^2}$로 유일하게 결정되며, 평균이 $1/N$인 극값 Fejér 다항식에서 유도된다.
• 최소한의 이전 상태 깊이는 $N^* = 2(N_0 - 1)$이며, 여기서 $N_0$는 $N_0 > \sqrt{\mu^*}$를 만족하는 최소 정수로, 모든 $\mu \in (-\mu^*, -1)$에 대해 안정성을 보장한다.
• 로지스틱 맵에서 $\mu^* = 4$일 경우 $N_0 = 3$이므로 $N^* = 4$이며, 최적의 계수는 $\varepsilon_1 = \frac{4}{9}, \varepsilon_2 = \frac{1}{9}$이다.
• 제어가 $h=4$에서 사이클을 안정화하지 못하는 이유는 승수 $\mu = -4$가 안정성 영역의 경계에 위치해 있고, 제어 스펙트럼이 단위 원에 접촉하기 때문이다.
• $\mu^* = 1$일 경우 $N_0 = 2$이므로 $N^* = 2$이며, 최적의 제어는 $u = -\frac{1}{4}(f_h(x_n) - f_h(x_{n-2}))$이다. 이 제어는 $h < 4$에서는 사이클을 안정화하지만 $h=4$에서는 그렇지 않다.
• 수치 시뮬레이션은 제어가 $h=3.95$에서 사이클을 안정화함을 확인하였으며, 해가 주기 궤도의 작은 이웃 영역으로 수렴함을 보여 이론적 결과를 검증한다.