[논문 리뷰] Stable limit cycles perturbed by noise

Stilianos Louca|arXiv (Cornell University)|2015. 06. 02.

Ecosystem dynamics and resilience참고 문헌 12인용 수 4

한 줄 요약

이 논문은 추가적인 가우시안 화이트 노이즈에 의해 흐트러진 안정한 리미트 사이클을 가진 연속시간 동역학계에 대해 형식적인 확률적 근사법을 유도한다. 이는 동역학을 단계(정적 속도에 의한 구동과 브라운 운동 유사한 확산을 포함)와 편차(다음 차원의 옴스타인-울렌벡 과정으로 모델링됨)로 분해한다. 주요 기여는 초임계 호프 분기 시스템에서 정(stationary) 자기공분산 및 파wer 스펙트럼 밀도에 대한 명시적 공식을 제공함으로써, 노이즈가 있는 순환 모델과 시계열 데이터 간의 정량적 비교를 가능하게 한다.

ABSTRACT

Many physical and biological systems exhibit intrinsic cyclic dynamics that are altered by random external perturbations. We examine continuous-time autonomous dynamical systems exhibiting a stable limit cycle, perturbed by additive Gaussian white noise. We derive a formal approximation for the dynamics of sample paths that stay close to the limit cycle, in terms of a phase coordinate and a deviation perpendicular to the limit cycle. To leading order in the deviation, the phase advances at the deterministic speed superimposed by a Brownian-motion-like drift. The deviation itself takes the form of an (n-1)-dimensional Ornstein-Uhlenbeck process. We apply these results to the case of limit cycles emerging through a supercritical Hopf bifurcation, which is widespread in ecological and epidemiological models. We derive approximation formulas for the system's stationary autocovariance and power spectral density. The latter two reflect the effects of perturbations on the temporal coherence and spectral bandwidth of perturbed limit cycles. We verify our results using numerical simulations and exemplify their application to the El Niño Southern Oscillation.

연구 동기 및 목표

• 연속시간 동역학계에서 안정한 리미트 사이클에 작용하는 추가적 가우시안 화이트 노이즈의 영향을 모델링하기 위해.
• 표본 경로를 리미트 사이클 근처에서 위상 및 수직 편차 좌표로 분해하여 분석적 취급 가능성을 확보하기 위해.
• 이전의 쿼asi-사이클 연구를 확장하여 위상 및 편차 동역학에 대한 근사적인 확률적 미분 방정식을 도출하기 위해.
• 초임계 호프 분기 시스템에서 정 자기공분산 및 파워 스펙트럼 밀도에 대한 명시적 공식을 제공하기 위해.
• 생물학적, 전염병학적, 기후학적 시스템과 같은 실제 시계열 데이터와 노이즈가 있는 리미트 사이클 모델 간 직접 비교를 가능하게 하기 위해.

제안 방법

• 리미트 사이클을 따라 움직이는 프레임(프레네티 유사)을 사용하여, 시간에 따라 변하는 직교 행렬 U(t)로 매개변수화함으로써 위상 및 편차 좌표를 정의한다.
• 이 프레임을 사용하여 SDE를 종방향(위상) 및 횡방향(편차) 성분으로 분해하며, 위상은 정적 속도에 의한 구동과 브라운 운동 유사한 확산을 포함한다.
• 편차를 (n−1)차원의 옴스타인-울렌벡 과정으로 모델링하여 사이클에 수직인 방향의 노이즈 유도된 진동을 포착한다.
• 정상형 초임계 호프 분기의 정형에 이 형식을 적용하여 노이즈가 있는 순환 동역학의 분석적 처리를 가능하게 한다.
• 위상-편차 분해를 사용하여 시스템의 정 자기공분산 및 파워 스펙트럼 밀도에 대한 폐쇄형 근사식을 도출한다.
• 2차 룬게-쿠타 스킴을 사용한 수치 시뮬레이션과 표본 경로로부터의 자기공분산 및 스펙트럼 밀도 추정을 통한 근사식의 검증을 수행한다.

실험 결과

연구 질문

• RQ1 추가적 가우시안 화이트 노이즈는 결정론적 시스템에서 안정한 리미트 사이클의 위상 동역학에 어떻게 영향을 미치는가?
• RQ2 리미트 사이클에서의 편차의 확률적 행동은 무엇이며, 다변량 옴스타인-울렌벡 과정으로 모델링할 수 있는가?
• RQ3 노이즈에 의한 영향, 예를 들어 주파수 이동 및 스펙트럼 넓어짐은 노이즈가 있는 리미트 사이클의 파워 스펙트럼 밀도에 어떻게 나타나는가?
• RQ4 초임계 호프 분기 하에서 노이즈가 작용하는 시스템에 대해 정 자기공분산 및 파워 스펙트럼 밀도에 대한 분석적 공식을 도출할 수 있는가?
• RQ5 유도된 근사식이 시뮬레이션된 노이즈가 있는 리미트 사이클 궤적의 통계적 성질을 얼마나 정확하게 재현하는가?

주요 결과

• 위상 좌표는 정적 속도에 의한 구동과 브라운 운동 유사한 확산을 포함하여 진화하며, 노이즈에 의한 위상 확산을 포착한다.
• 리미트 사이클에서의 편차는 (n−1)차원의 옴스타인-울렌벡 과정으로 잘 근사되며, 평균 회귀 동역학을 갖는다.
• 초임계 호프 분기 시스템에서는 정 자기공분산 함수가 진동 주파수에 의해 조절되는 지수적 감쇠를 나타낸다.
• 파워 스펙트럼 밀도는 결정론적 주파수에서 특징적인 피크를 보이며, 이 밴드폭은 노이즈 강도에 비례하고 시스템의 안정성 마진에 반비례한다.
• 수치 시뮬레이션을 통해 도출된 분석적 근사식이 다양한 노이즈 대 신호 비율에서 노이즈가 있는 리미트 사이클의 스펙트럼 및 상관 성질을 정확하게 예측함을 확인하였다.
• 모델은 엘니뇨 남부 진동에서의 분해효과를 성공적으로 포착하여 실제 기후 시계열에의 적용 가능성을 입증하였다.