[논문 리뷰] A New Test for Chaos

Georg A. Gottwald, Ian Melbourne|arXiv (Cornell University)|2002. 08. 26.

Chaos control and synchronization참고 문헌 16인용 수 309

한 줄 요약

이 논문은 단일 시간 시리즈 데이터에 직접 적용되며 위상 공간 재구성이나 시스템의 기본 방정식에 대한 지식이 필요 없는 결정론적 동역학계에서의 혼돈을 진단하기 위한 새로운 0–1 테스트를 제안한다. 이 방법은 관측 가능량과 단위 각도 변수로부터 유도된 시간 평균 함수를 사용하여 평균 제곱 이동 성장률 K를 계산하며, K ≈ 0은 비혼돈 동역학, K ≈ 1은 혼돈 동역학을 나타낸다. 이는 ODE, PDE, 맵에 모두 적용 가능한 보편적이고 저비용의 진단 방법을 제공한다.

ABSTRACT

We describe a new test for determining whether a given deterministic dynamical system is chaotic or nonchaotic. (This is an alternative to the usual approach of computing the largest Lyapunov exponent.) Our method is a 0-1 test for chaos (the output is a 0 signifying nonchaotic or a 1 signifying chaotic) and is independent of the dimension of the dynamical system. Moreover, the underlying equations need not be known. The test works equally well for continuous and discrete time. We give examples for an ordinary differential equation, a partial differential equation and for a map.

연구 동기 및 목표

• 시스템의 방정식이나 위상 공간 재구성에 대한 지식이 필요 없는 보편적이고 계산 비용이 낮은 혼돈 진단법을 개발하기 위해.
• ODE, PDE, 맵 포함 고차원 또는 복잡한 시스템에서 시간 시리즈 데이터만으로도 혼돈 탐지를 가능하게 하기 위해.
• 비혼돈 또는 혼돈 행동을 나타내는 이진 진단(0 또는 1)을 제공하며, 성장률 K를 통해 명확한 수치적 구분을 제공하기 위해.
• 고차원 또는 실험 데이터에서 라플라스 지수 계산의 한계를 극복하기 위해 선형화 및 임bedding을 피하기 위해.

제안 방법

• 이 방법은 $ p(t) = \int_0^t \phi(\mathbf{x}(s)) \cos(\theta(s)) \, ds $ 를 정의하며, 여기서 $ \theta(s) = cs + \int_0^s \phi(\mathbf{x}(u)) \, du $ 이고, $ c > 0 $ 는 임의로 선택된다.
• 관측 가능량 $ \phi(\mathbf{x}) $ 가 선택되며, 일반적으로 $ x_1 $ 과 같은 단순한 성분이 사용된다. 이는 $ p(t) $ 의 동역학을 이끌어낸다.
• 평균 제곱 이동(MSD) $ \mathbf{M}(t) $ 는 $ p(t) $ 의 제곱 증분의 시간 평균으로 계산되며, 일시적인 영향에 대한 강건성을 확보한다.
• 성장률 $ K = \lim_{t \to \infty} \log(\mathbf{M}(t) + 1) / \log t $ 는 $ \log \mathbf{M}(t) $ 와 $ \log t $ 의 선형 회귀를 통해 추정된다.
• 진단값 $ K \approx 0 $ 은 $ p(t) $ 가 유계임을 나타내며, 비혼돈 동역학을 의미한다. $ K \approx 1 $ 은 확산 성장임을 나타내며, 혼돈 동역학을 의미한다.
• 선형 이동을 제거하기 위해 $ \mathbb{R}^2 $ 확장(즉, $ \theta(t) $ 를 통한)을 사용하여 정규 동역학과 혼돈 동역학 간의 구분이 흐려지는 것을 방지한다.

실험 결과

연구 질문

• RQ1 위상 공간 재구성이나 시스템 방정식에 대한 지식이 없이도 시간 시리즈 데이터에 직접 적용 가능한 이진 혼돈 테스트를 구성할 수 있는가?
• RQ2 p(t) 의 평균 제곱 이동 성장률 K 는 혼돈 대비 비혼돈 동역학을 신뢰성 있고 보편적으로 나타내는가?
• RQ3 이 테스트는 ODE, PDE, 맵 포함 고차원 시스템에서 낮은 계산 오버헤드로 혼돈을 구분할 수 있는가?
• RQ4 p(t) 정의에 단위 각도 변수 $ \theta(t) $ 를 포함함으로써, 기존에 존재하던 선형 이동이 어떻게 정규 운동과 혼돈 운동 간의 구분을 가리키지 않도록 제거되는가?

주요 결과

• 비혼돈 동역학에서는 $ p(t) $ 가 유계를 유지하므로 $ \mathbf{M}(t) $ 도 유계이며 결과적으로 $ K \approx 0 $ 이다.
• 혼돈 동역학에서는 $ p(t) $ 가 渐近적으로 브라운 운동과 유사하게 행동하여 $ \mathbf{M}(t) \sim t $ 가 되며, 따라서 $ K \approx 1 $ 이다.
• 이 방법은 차원이나 방정식 형태에 관계없이 어떤 결정론적 시스템에도 보편적으로 적용 가능하다. ODE, PDE, 맵 모두 포함된다.
• 이 테스트는 $ p(t) $ 와 $ \theta(t) $ 를 위한 추가로 두 개의 미분방정식만 필요로 하여, n차원 시스템에 대해 $ n^2 $ 개의 방정식이 필요한 라플라스 지수 계산보다 훨씬 효율적이다.
• 관측 가능량 $ \phi $ 의 선택에 대해 매우 강건하며, 거의 모든 비퇴화된 $ \phi $ 가 정확한 분류를 제공한다.
• 강제 진동자인 반데르폴 진동자에 대한 수치 결과는 $ \omega $ 가 변화함에 따라 $ K \approx 0 $ (비혼돈)과 $ K \approx 1 $ (혼돈) 영역 간 명확한 분리가 있음을 보여준다.