비가우시안 1/f 잡음 모델링을 위한 확률 미분 방정식

초록

본 논문은 완화율이 확산적으로 변동하는 선형 확률 미분 방정식을 도입하여, 비가우시안 파워‑러프 분포와 1/f^β 형태의 잡음을 동시에 구현하는 단일 프랙탈 신호 모델을 제시한다. 모델의 수학적 유도와 수치 시뮬레이션을 통해 확률적 완화율 변동이 신호의 꼬리 분포와 스펙트럼 특성에 미치는 영향을 분석하고, 기존 모델과의 차별성을 강조한다.

상세 분석

논문은 먼저 전통적인 선형 SDE dx = −γ x dt + σ dW_t 를 회복한다. 여기서 γ는 고정된 완화율이며, 이 경우 x(t)는 가우시안 정규분포를 따르고 파워 스펙트럼은 1/f^2 형태의 로렌츠 곡선을 보인다. 저자들은 γ 자체를 확률 과정 γ(t)로 확장한다는 점에서 혁신적이다. 구체적으로 γ(t)는 평균 γ₀와 확산 계수 D_γ를 갖는 오스코르츠(Ornstein‑Uhlenbeck) 과정으로 모델링되며, dγ = −λ(γ − γ₀)dt + √(2D_γ) dW’_t 로 기술된다. 이때 λ는 완화율 변동의 회복 속도, W’_t는 독립적인 백색 잡음이다.

이중 확률 과정은 x(t)의 확률 밀도 함수(PDF)를 조건부 가우시안 형태로 유지하면서, γ의 변동에 의해 전체 PDF가 혼합 가우시안으로 전이한다. 즉, P(x) = ∫ P(x|γ) P(γ) dγ 로 표현되며, γ가 넓은 분포를 가질 경우 P(x)는 꼬리가 두꺼운 레비‑스틸러형 파워‑러프를 나타낸다. 저자들은 이론적 계산을 통해 P(x) ∝ |x|^{−(1+α)} (α > 0) 형태의 비가우시안 분포가 도출됨을 보였다.

스펙트럼 분석에서는 γ(t)의 저주파 변동이 x(t)의 장기 상관성을 강화한다는 점을 강조한다. γ의 자기상관함수는 exp(−λ|τ|) 형태이므로, λ가 작을수록 γ는 느리게 변동하고, x(t)의 자기상관함수는 ∝ |τ|^{−β} (0 < β < 1) 로 근사된다. 따라서 파워 스펙트럼 S(f) ∝ 1/f^β 가 나타나며, β는 λ와 D_γ의 비율에 의해 조절된다. 이는 기존의 단일 γ 모델이 제공할 수 없는 1/f 잡음 현상을 자연스럽게 재현한다.

수치 시뮬레이션에서는 Euler‑Maruyama 방법을 이용해 10⁶ 단계 이상을 수행했으며, γ₀ = 1, λ = 10⁻³, D_γ = 10⁻⁴, σ = 0.5 등 다양한 파라미터 조합을 실험했다. 결과는 이론적 예측과 일치했으며, 특히 γ의 변동 폭이 커질수록 PDF의 꼬리가 두꺼워지고 β가 0.5 ~ 0.9 사이로 변동함을 확인했다. 또한, γ의 초기값 분포가 균등이든 정규이든 최종 통계적 특성에는 큰 차이가 없음을 보고했다.

논문의 주요 기여는 다음과 같다. 첫째, 완화율을 확률 과정으로 취급함으로써 비가우시안 분포와 1/f 잡음을 동시에 생성하는 최소 모델을 제시했다. 둘째, γ(t)의 동역학 파라미터(λ, D_γ)가 신호의 꼬리 지수 α와 스펙트럼 지수 β를 직접 제어한다는 정량적 관계를 도출했다. 셋째, 기존의 비선형 SDE 기반 모델과 달리 선형 구조를 유지하면서도 복잡한 통계적 현상을 설명할 수 있음을 증명했다. 이러한 접근은 전자소자, 생물학적 신호, 금융 시계열 등 다양한 분야에서 관측되는 1/f 잡음과 비가우시안 통계의 근본 메커니즘을 탐구하는 데 유용한 이론적 틀을 제공한다.