양의 피드백이 결합된 일반화 CEV 과정의 동역학

초록

본 논문은 일반화 CEV 과정 (dX_t = aX_t^{,n}dt + bX_t^{,m}dW_t) 에서 드리프트와 확산에 대한 양의 상태‑피드백((n,m>1))이 동시에 작용할 때, 확률적 정상이 존재하는 조건과 그 특성을 분석한다. (n<2m-1)이면 정상분포가 존재하며 꼬리는 거듭 제곱법칙 (P(x)\sim x^{-2m}) 로 감소한다. 파워 스펙트럼은 (S(f)\sim f^{-\beta}) 형태이며 (\beta = 2-\frac{1+\epsilon}{2(m-1)}) 로 정의된다. 수치 실험을 통해 버스트 현상을 확인하고, 버스트 강도와 지속시간이 (S\propto T^{2}) 관계임을 제시한다.

상세 분석

일반화 CEV(gCEV) 과정은 상태 변수 (X_t)가 시간에 따라 비선형적인 드리프트 항 (aX_t^{,n})와 비선형 확산 항 (bX_t^{,m})에 의해 진화한다는 점에서 기존의 선형 확산 모델을 확장한다. 여기서 (n)과 (m)은 각각 상태‑to‑drift와 상태‑to‑diffusion 피드백의 차수를 나타내며, 두 차수가 1보다 클 때 양의 피드백이 존재한다. 양의 피드백은 작은 값에서는 억제 효과가 약해지지만, 큰 값으로 갈수록 드리프트와 확산이 급격히 강화돼 폭발적 성장(burst) 가능성을 만든다.

정상분포 존재 조건을 분석하기 위해 Fokker‑Planck 방정식에 대한 정적 해를 구하면, 확률 흐름이 무한대에서 수렴하려면 (n<2m-1)이어야 함을 알 수 있다. 이 조건 하에서 정상분포는 꼬리 부분이 (P(x)\sim x^{-\mu}) 형태이며, (\mu=2m)가 된다. 즉, 확산 차수 (m)이 클수록 꼬리가 더 가파르게 감소해 큰 값의 발생 확률이 급격히 억제된다.

주파수 영역에서는 파워 스펙트럼 밀도 (S(f))가 저주파에서 (f^{-\beta}) 형태를 보이며, (\beta)는 확산 차수와 작은 양의 파라미터 (\epsilon)에 의해 결정된다. (\beta=2-\frac{1+\epsilon}{2(m-1)})는 (m)이 증가할수록 2에 가까워져 1/f 노이즈에 근접하고, (m)이 1에 가까우면 (\beta)가 크게 감소해 백색 잡음에 가까워진다. 이는 시스템이 얼마나 장기 기억을 유지하는지를 피드백 차수와 직접 연결시킨다.

수치 시뮬레이션에서는 시간 연속적인 버스트 현상이 관찰된다. 버스트는 (X_t)가 급격히 상승한 뒤 다시 낮은 수준으로 복귀하는 일련의 이벤트이며, 버스트 강도 (S)와 지속시간 (T) 사이에 (S\propto T^{2}) 관계가 성립한다. 이는 드리프트와 확산이 동시에 강화되는 구간에서 에너지(또는 면적) 축적이 시간에 대해 이차적으로 증가한다는 물리적 직관과 일치한다. 이러한 스케일링 법칙은 금융 시장의 급등락, 신경 활동의 스파이크, 혹은 기후 변동 등 다양한 복합 시스템에 적용 가능성을 시사한다.

전체적으로 본 연구는 양의 비선형 피드백이 결합된 확률 과정이 어떻게 정상분포와 파워‑law 스펙트럼을 동시에 구현하면서도 불규칙한 버스트를 생성하는지를 이론적·수치적으로 명확히 밝히며, 비선형 확산 모델링에 새로운 통찰을 제공한다.