2025-11-19T03:34:12.934208

On Approximate Representation of Fractional Brownian Motion

Rybakov

This paper considers the orthogonal expansion of the fractional Brownian motion relative to the Legendre polynomials. Such an expansion has not only theoretical but also practical interest, since it can be applied to approximate and simulate the fractional Brownian motion in continuous time. The relations for the mean square approximation error are presented, and a comparison with the previously obtained result is carried out.

academic

분수 브라운 운동의 근사 표현에 관하여

기본 정보

논문ID: 2503.04575
제목: On Approximate Representation of Fractional Brownian Motion
저자: Konstantin A. Rybakov
분류: math.PR (확률론)
발표 시간: 2025년 3월 (arXiv 사전인쇄본, 버전 v2는 2025년 10월 14일)
논문 링크: https://arxiv.org/abs/2503.04575

초록

본 논문은 르장드르 다항식에 대한 분수 브라운 운동의 직교 전개를 연구한다. 이러한 전개는 이론적 의의뿐만 아니라 연속 시간에서 분수 브라운 운동의 근사 및 모의실험에 사용될 수 있다는 점에서 실용적 가치를 갖는다. 본 논문에서는 평균제곱 근사 오차의 관계식을 제시하고 이전에 얻은 결과와 비교한다.

연구 배경 및 동기

핵심 문제: 연속 시간에서 분수 브라운 운동(fractional Brownian motion, FBM)을 효과적으로 근사하고 모의실험하는 방법
중요성: 분수 브라운 운동은 금융 모델링, 물리학, 공학 등 여러 분야에서 광범위하게 적용되며, 이의 모의실험은 이론적 결과 검증 및 실제 문제 해결에 필수적
기존 방법의 한계:
- 전통적인 이산 시간 격자 방법은 격자 노드에서만 정확성을 보장
- 기존의 일반적인 가우스 과정 모의실험 방법은 계산 복잡도가 높음
- 근사 오차를 정확히 계산할 수 있는 연속 시간 방법 부재
연구 동기: 저자의 이전 연구2를 기반으로 더 간단한 전개 계수 계산 방법을 제시하여 계산 복잡도를 낮추고 정확도를 향상

핵심 기여

이론적 기여: 분수 브라운 운동의 적분 표현과 관련된 핵 함수 전개 계수의 명시적 관계식 도출 (르장드르 다항식에 대해)
방법 개선: 이전 연구2보다 간단한 전개 계수 계산 방법 제공으로 계산 복잡도 감소
오차 분석: 분수 브라운 운동의 다항식 근사 평균제곱 오차를 정확히 계산할 수 있는 방정식 수립
연속 시간 근사: 이산 시간 점에서만이 아닌 연속 시간에서 분수 브라운 운동의 근사 표현 실현

방법 상세 설명

작업 정의

허스트 지수 $H \in (0,1)$ 를 갖는 분수 브라운 운동 $B_H(t)$ 가 주어졌을 때, 르장드르 다항식 직교 기저에 대한 전개 계수를 찾아 연속 시간 근사를 실현:

$B_H(t) = \sum_{i=0}^{\infty} B_H^i \hat{P}(i,t)$

여기서 $\hat{P}(i,t)$ 는 표준화된 르장드르 다항식이다.

핵심 수학 프레임워크

분수 브라운 운동 정의:

공분산 함수: $R_H(t,\tau) = \frac{t^{2H} + \tau^{2H} - |t-\tau|^{2H}}{2}$
적분 표현: $B_H(t) = \int_0^t k_H(t,\tau) dB(\tau)$

핵 함수: $k_H(t,\tau) = a_H(t-\tau)^{H-1/2} \,_2F_1\left(\frac{1}{2}-H, H-\frac{1}{2}, H+\frac{1}{2}, 1-\frac{t}{\tau}\right)\mathbf{1}(t-\tau)$

여기서 $a_H = \sqrt{\frac{\pi H(1-2H)}{\Gamma(2-2H)\cos\pi H}}$ ( $H \neq 1/2$ 일 때).

주요 이론적 결과

정리 1 (전개 계수의 명시적 표현): 행렬 $K^H$ 의 제 $j$ 열에 대해: $K^H_{*j} = a_H\Gamma(3/2-H)\sqrt{\frac{2j+1}{T}} \sum_{k=0}^j \frac{l_{jk}(3/2-H)_k}{T^k(H+1/2+k)k!} F^{H+1/2+k}$

$H \neq 1/2$ 에 대해, 전개 계수는 다음을 만족: $K^H_{ij} = a_H\Gamma(1/2-H)\sqrt{\frac{2j+1}{T}} F^{H+1/2}_i \sum_{k=0}^j (-1)^{j-k} \frac{1/2-H+k}{H+1/2+k} \Pi^{(ij)}_k$

알고리즘 구현

계수 계산: 재귀 관계식 (12)-(13)을 사용하여 거듭제곱 함수의 전개 계수 계산
행렬 구성: 정리 1을 통해 전개 계수 행렬 $K^H$ 구성
난수 생성: 독립 표준 정규 난수 $V_i$ 생성
근사 구성: $\tilde{B}_H(t) = \sum_{i=0}^{L-1} \tilde{B}^H_i \hat{P}(i,t)$

실험 설정

매개변수 설정

시간 구간: $T = [0,1]$
허스트 지수: $H = 0.1, 0.2, \ldots, 0.9$
절단 길이: $L = 4, 8, 16, 32, 64, 128$

평가 지표

평균제곱 근사 오차: $\varepsilon = E\int_T (B_H(t) - \tilde{B}_H(t))^2 dt$

이론적 공식: $\varepsilon = \|k_H(\cdot)\|^2_{L^2(T^2)} - \|\bar{K}^H\|^2$

여기서 $\|k_H(\cdot)\|^2_{L^2(T^2)} = \frac{T^{2H+1}}{2H+1}$

실험 결과

주요 결과

표 1: 평균제곱 근사 오차 $\varepsilon$

H	L=4	L=8	L=16	L=32	L=64	L=128
0.1	0.3842	0.3229	0.2720	0.2299	0.1950	0.1659
0.3	0.1035	0.0655	0.0423	0.0275	0.0180	0.0118
0.5	0.0357	0.0167	0.0081	0.0040	0.0020	0.0010
0.7	0.0132	0.0049	0.0019	0.0008	0.0003	0.0001
0.9	0.0812	0.0617	0.0469	0.0356	0.0270	0.0205

주요 발견

수렴성 분석:
- $H = 0.5$ 일 때 수렴이 가장 빠름
- $H$ 가 경계값(0 또는 1)에 가까워질수록 수렴 속도가 현저히 감소
- 상대 오차는 $\varepsilon/(2H+1)^{-1}$
방법 비교:
- 새로운 방법과 문헌2의 결과가 질적으로 유사
- 작은 $L$ 값에 대해서도 두 방법의 오차가 매우 가까움
- 행렬 곱 표현이 정확도에 미치는 영향이 미미하지만 계산 안정성 제공
수렴율 분석:
- 거듭제곱 함수 $t^\alpha$ 의 전개 수렴율은 $O(L^{-(2\alpha+1)})$
- $\alpha \to -1/2$ 일 때 수렴율이 0으로 수렴하여 경계 허스트 지수의 낮은 수렴성 설명