기술 요약: 분수 차수 시스템의 카오스 어트랙터와 리아푸노프 지수

1. 문제 제기

비정수 차수 $\nu \in (0, 1)$의 카푸토 미분(Caputo derivative)으로 정의되는 분수 차수 동역학 시스템은 현재의 변화율이 과거 전체 이력에 멱법칙 가중치를 두어 의존하는 과정을 모델링한다. 이러한 시스템(예: 분수 로렌츠, 로슬러, 첸 시스템)은 카오스를 나타내는 것으로 알려져 있지만, 리아푸노프 지수를 통해 이 카오스를 측정하는 것은 여

동안 문제가 되어 왔다. 표준 리아푸노프 알고리즘은 미래의 진화를 결정하는 "상태(state)"라는 개념에 의존한다. 그러나 분수 차수 시스템에는 그러한 상태가 존재하지 않는다. 이들의 미래는 전체 이력에 의존하기 때문이다. 결과적으로 이 시스템은 전통적인 의미의 흐름(flow)을 갖지 못한다. 기존 방법들은 메모리를 절단(truncation)하여 상당한 오차를 유발하거나, 절단된 스킴에 표준 알고리즘을 적용하여 완전히 다른 시스템의 지수를 계산한다. 또한, 메모리가 노화됨에 따라 메모리의 노름(norm)이 유계된 어트랙터 위에서도 무한히 성장할 수 있기 때문에, 분수 차수 시스템에 대한 리아푸노프 지수의 수학적 정의는 매우 까다롭다.

2. 방법론: 증강 시스템 알고리즘 (Augmented-System Algorithm)

본 논문은 "확산 표현(diffusive representation)"에 기반한 실용적인 해결책을 제안하며, 이는 분수 차수 시스템을 무한 차원 상미분 방정식(ODE) 시스템으로 변환한다. 핵심 방법론은 다음 세 가지 기둥에 기초한다.

A. 확산 표현

멱법칙 메모리 커널 $k_\nu(\tau) = \tau^{\nu-1}/\Gamma(\nu)$는 감쇠하는 지수 함수들의 중첩으로 표현된다:
$$k_\nu(\tau) = \int_0^\infty \mu_\nu(\omega) e^{-\omega\tau} d\omega$$
여기서 $\mu_\nu(\omega) \propto \omega^{-\nu}$이다. 이를 통해 분수 차수 시스템을 "증강 시스템"으로 재작성할 수 있으며, 여기서 상태는 메모리 모드를 나타내는 함수 $\omega \mapsto z(\omega, t)$가 된다. 각 모드는 다음과 같은 1차 ODE를 만족한다:
$$\partial_t z(\omega, t) = -\omega z(\omega, t) + f(x(t))$$
원래의 변수 $x(t)$는 이러한 모드들의 가중치 합을 통해 복구된다. 이 정식화는 분수 차수 시스템을 무한 차원 상태를 가진 자율 ODE 시스템으로 변환하여, 표준 동역학 도구를 적용할 수 있게 한다.

B. 이산화 및 타임 스텝

이를 수치적으로 구현하기 위해, 모드의 연속체를 적분을 근사하는 쿼드러처 규칙(quadrature rule)을 사용하여 이산화한다. 본 논문은 $\omega=0$에서의 $\mu_\nu(\omega)$의 특이성을 처리하기 위해 특정 전략을 사용한다:

• 특이 패널(Singular Panel): 장기적 거동에 필수적인 느리고 노화되는 모드들을 포착하기 위해 $[0, \omega_{min}]$ 구간에 Gauss-Jacobi 규칙을 사용한다.
• 기하학적 패널(Geometric Panels): $\omega_{max}$까지의 이후 십진수 구간들에 대해 $\ln(\omega)$ 변수에서 Gauss-Legendre 규칙을 사용한다.
• 지수 시간 차분(Exponential Time Differencing, ETD): 시스템이 매우 뻣뻣하기(stiff) 때문에(변수 $\omega$의 범위가 $10^{-7}$에서 $10^7$에 달함), 모드 방정식의 선형 부분은 해석적으로 적분하고, 비선형 결합은 2차 룬게-쿠타(Runge-Kutta) 스킴을 사용하여 수치적으로 처리한다.

C. 접선 시스템과 노름

리아푸노프 지수를 계산하기 위해, 증강 시스템의 선형화된 역학을 궤적과 함께 적분한다. 직교화(QR 분해를 통한)를 위한 노름의 선택에는 결정적인 이론적 통찰이 필요하다:

• 선형화된 역학은 순수 완화(노름 1)와 컴팩트한 결합으로 구성된다.
• 이론은 분수 차수 시스템이 음의 리아푸노프 지수를 갖지 않는다고 예측한다. 섭동이 감쇠할 때는 지수적으로가 아니라 대수적으로(시간 $t$의 멱함수로서) 감쇠한다.
• 따라서 접벡터들은 특정 **가중 노름(weighted norm)**에서 직교화되어야 한다: $\|( \delta a, \delta z )\|^2 = |\delta a|^2 + \sum w_i(1+\omega_i)|\delta z_i|^2$. 표준 유클리드 노름을 사용하면 잘못된 결과를 얻게 된다.

3. 주요 기여

본 논문은 이 알고리즘을 구현하는 완전하고 자기 완결적인 Python 프로그램(약 150행)을 제공한다. 구체적인 기여는 다음과 같다:

1. 작동하는 구현체: 분수 카오스 시스템을 시뮬레이션하고 리아푸노프 지수를 계산하는 스크립트를 제공하며, 알려진 정확한 해를 가진 선형 시스템을 통해 검증되었다.
2. 코드 설계에 대한 이론적 정당성: 왜 특정 설계 선택(특이 패널 포함 및 특정 가중 노름 등)이 필요한지를 설명하며, 이를 확산 표현의 수학적 특성과 직접 연결한다.
3. 출력 해석의 명확화: 분수 차수 계산 결과를 읽는 방법을 확립한다:
   • 양의 지수는 잘 정의되어 있으며 이산적이다.
   • 음이 아닌 지수는 음의 숫자로 수렴하지 않고 0을 향해 천천히 표류(drift)한다.
   • 소멸하는 지수($\lambda \approx 0$)의 경우, 평형점, 주기적 창(periodic window), 또는 과도기(transient)를 구분하기 위해 궤적 검사가 필요하다.

4. 결과

프로그램은 다음과 같은 시스템들에 적용되어 검증되었다:

• 선형 검증: $\nu=0.8$인 선형 시스템 $CD^\nu x = \text{diag}(0.05, -1, -2)x$에 대해, 프로그램은 최대 지수를 $0.023636$으로 계산하였으며, 이는 이론값인$0.05^{1/0.8} \approx 0.023644$와$0.04%$이내의 오차로 일치하였다. 다른 지수들은$-1$과$-2$로 수렴하는 대신 0을 향해 올바르게 표류하였다.
• 분수 로렌츠 시스템 ($\nu=0.995$): 프로그램은 최대 지수 $\lambda_1 \approx 0.856$을 계산하였다(정수 차수 시스템의 $0.906$과 비교). 이는 카오스가 시작되기 전 궤적이 평형점 근처에서 머무는 현상을 정확히 포착하였다.
• 분수 첸 시스템 ($\nu=0.9$): 프로그램은 $\lambda_1 \approx 2.33$을 찾아냈다.
• 주기적 창(Periodic Windows): 알고리즘은 첸 시스템의 $\nu=0.99$ 지점에서 주기적 창을 성공적으로 식별하였으며, 여기서 지수는 0으로 표류하고 궤적은 긴 과도기 후에 주기적이 되었다.

5. 의의 및 주장

본 논문은 최근의 이론적 논문의 실용적인 동반자 역할을 자처하며, 분수 차수 시스템의 카오스 측정에 대한 모호함을 해결하는 것을 목표로 한다. 그 의의는 다음과 같다:

• 오해 교정: 표준 리아푸노프 알고리즘을 절단된 메모리 스킴에 적용하는 것은 잘못된 시스템을 측정하는 것임을 입증한다.
• 상태의 정의: "분수 차수 시스템의 상태는 그 메모리이다"라는 개념을 구체화하여, 표준 동역학 이론의 적용을 가능하게 한다.
• 절단 오차 방지: 메모리 절단(느린 모드 제거)이 시스템을 근본적으로 변화시키고 지수에 큰 오차를 유발함을 커널 오차 분석을 통해 명시적으로 경고한다.
• 제한된 범위: 저자는 모든 분수 차수 동역학 문제를 해결하려는 것이 아니라, $CD^\nu x = f(x)$ 형태의 시스템에 대한 카오스 어트랙터와 그 지수를 분석하기 위한 견고한 도구를 제공하는 데 집중한다. 논문은 독자가 출력값을 올바르게 해석하기 위해 (분수 차수 지수의 특징인) 0 근처의 띠(band)를 이해해야 함을 강조한다.