[논문 리뷰] On the fractal dimension of the Duffing attractor

Mariusz Tarnopolski|arXiv (Cornell University)|2013. 05. 29.

Mathematical Dynamics and Fractals참고 문헌 9인용 수 5

한 줄 요약

이 연구는 상자 세기 및 상관 차원 방법을 사용하여 수치적으로 Duffing 혼돈 궤도의 분수차원을 추정하며, 샘플링된 점의 수가 증가함에 따라 안정된 값으로 수렴하는 것을 보여준다. 결과적으로 상자 세기 차원은 dC ≈ 1.43, 상관 차원은 dG ≈ 1.38로 도출되었으며, 이는 위상 공간 내 비정수 분수적 구조를 가진 이상한 자기유사 혼돈 궤도임을 시사한다.

ABSTRACT

The box counting dimension $d_C$ and the correlation dimension $d_G$ change with the number of numerically generated points forming the attractor. At a sufficiently large number of points the fractal dimension tends to a finite value. The obtained values are $d_C\approx 1.43$ and $d_G\approx 1.38$.

연구 동기 및 목표

• 수치적으로 생성된 점의 수가 증가함에 따라 Duffing 혼돈 궤도의 분수차원 추정치의 수렴 행동을 조사하기 위해.
• 혼돈 영역에서 Duffing 혼돈 궤도의 상자 세기 차원(dC)과 상관 차원(dG)의 안정된 값을 도출하기 위해.
• 비선형 역학에서 널리 사용되는 두 방법을 이용하여 분수차원을 추정함으로써 Duffing 혼돈 궤도의 분수적 성질을 검증하기 위해.
• 유한한 점 집합에서 상관 차원 알고리즘의 수치적 복잡성과 신뢰성에 대해 상자 세기 방법과 비교하여 평가하기 위해.
• 혼돈 궤도의 위상 공간 내 기하학적 복잡성에 대한 정량적 측도를 제공하여 이상한 궤도로의 분류를 뒷받침하기 위해.

제안 방법

• Mathematica의 기본 방법을 사용하여 Duffing 방정식 (¨x + δẋ + βx + αx³ = γ cos(ωDt))의 수치적 적분을 수행하여 궤도를 생성하기 위해.
• 위상 변수로 위치, 속도, 외부력의 위상 변수를 사용하여 이阶 미분방정식을 3차원 자율 시스템으로 변환하기 위해.
• 외부력 주기마다 샘플링하는 Poincaré 표면을 이용하여 궤도의 구조를 추출하기 위해.
• Mathematica의 BoxCount 함수를 사용하여 상자 세기 차원(dC)을 계산하고, 작은 ε에 대해 ln N(ε) 대 ln(1/ε)의 선형 회귀를 수행하기 위해.
• 작은 값의 r 범위에서 C(r) = (1/N²)Σᵢ<ⱼ H(r − ||xi − xj||)를 계산하는 파이썬 프로그램을 병렬로 구현하여 상관 차원(dG)을 계산하기 위해.
• log C(r) 대 log r 및 log N(ε) 대 log(1/ε)의 선형 회귀를 수행하여 기울기를 추정하고, 이를 분수차원으로 해석하며 오차 최소화 기법을 적용하기 위해.

실험 결과

연구 질문

• RQ1 Duffing 혼돈 궤도의 상자 세기 차원(dC)과 상관 차원(dG)은 샘플링된 점의 수가 증가함에 따라 어떻게 변화하는가?
• RQ2 충분히 큰 점 집합에서 분수차원 추정치가 유한한 값으로 수렴하는가? 이는 진정한 분수적 구조를 시사하는가?
• RQ3 주어진 혼돈 매개변수 조건 하에서 Duffing 혼돈 궤도의 수렴된 상자 세기 차원 dC의 값은 얼마인가?
• RQ4 Duffing 혼돈 궤도의 수렴된 상관 차원 dG의 값은 얼마이며, dC와 비교해보면 어떠한가?
• RQ5 상관 차원 알고리즘의 수치적 복잡성은 점의 수에 따라 어떻게 변화하는가? 대규모 데이터셋에 대해 계산 가능성이 있는가?

주요 결과

• 점의 수가 증가함에 따라 상자 세기 차원 dC는 약 1.43으로 수렴하며, 이는 1과 2 사이의 분수적 구조를 가진다는 것을 시사한다.
• 상관 차원 dG는 약 1.38로 수렴하며, 이는 비정수적이고 자기유사적인 기하학을 가진 이상한 궤도의 존재를 뒷받침한다.
• N ≥ 1592점일 경우, log C(r) 대 log r 및 log N(ε) 대 log(1/ε)의 선형 회귀에서 dG ≈ 1.357 및 dC ≈ 1.330로 안정적이고 일관된 추정치를 도출하였다.
• 추정된 상관 차원의 표준편차는 한 점이 추가될 때 21% 증가하였으며, 이는 N의 미세한 변화에 민감함을 시사하지만, 차원 값 자체는 오직 4% 정도 변화하였다.
• 상관 차원 알고리즘은 O(n²)의 계산 복잡도를 보이며, 실행 시간에 이차 다항식을 피팅하여 확인하였고, 31,838점에 대해 예측 시간은 12,943초였으며, 실제 시간인 12,610초와 매우 유사하였다.
• 두 방법 모두 유한하고 비정수적인 값으로 수렴함을 확인하였으며, 이는 Duffing 혼돈 궤도의 분수적 성질을 확인하고, 소산적 혼돈 시스템 내 이상한 궤도로의 분류를 뒷받침한다.