가중치 비앵커드 Sobolev 공간에서 무한 영역 적분을 위한 중위값 QMC 방법

초록

본 논문은 다변량 적분 I₍ϕ₎(f)=∫_{ℝ^s}f(y)ϕ(y)dy 를 위해, 가중치 비앵커드 Sobolev 공간 F 에 정의된 함수에 대해, 무작위 이동 격자 규칙을 k=O(log N)번 독립적으로 수행한 뒤 그 중위값을 취하는 새로운 QMC 추정기를 제안한다. 평균 절대 오차는 O(N^{-r+ε}) (ε∈(0,r−½])를 만족하며, 이는 기존 CBC 기반 격자 규칙과 동일한 수렴률이면서 가중치 구조에 대한 사전 지식이 필요하지 않다. 실험 결과는 이 방법이 CBC와 비슷한 정확도를 보이며 단순한 Monte Carlo보다 우수함을 확인한다.

상세 분석

본 연구는 무한 영역 ℝ^s 위의 적분 문제를 QMC 방법으로 해결하고자 하는데, 핵심 아이디어는 두 단계의 무작위화를 결합한 ‘중위값 QMC’이다. 첫 번째 무작위화는 전통적인 랜덤 시프트(lattice shift)이며, 두 번째는 생성 벡터 z∈G_{s,N} 을 균등하게 무작위 선택하는 것이다. 각각의 (z,Δ)쌍에 대해 평균값 Q_{P_{z,Δ}}(f)=\frac1N∑_{n=0}^{N-1}f∘Φ^{-1}({nz/N+Δ}) 을 계산하고, k개의 독립 추정값 중 중위값 M_k(f) 을 최종 추정치로 채택한다.

이때 함수 공간 F는 Nichols와 Kuo가 제안한 가중치 비앵커드 Sobolev 공간으로, 노름
‖f‖F^2 = ∑{u⊆{1,…,s}} γ_u^{-1} ∫{ℝ^{|u|}} |∂^{|u|}f/∂y_u|^2 ∏{j∈u} ψ_j^2(y_j) dy_u
에 의해 정의된다. 여기서 γ_u는 가중치, ψ_j는 밀도 ϕ 에 의존하는 보조 가중치 함수이며, 강한 조건(2.1)을 만족하면 F 는 RKHS이면서 L^2_ϕ 에 포함된다.

논문은 먼저 랜덤 시프트 격자 규칙의 평균 최악오차 e_{sh,s,N}(z) 를 Fourier 계수 b_{θ_j}(h) 와 가중치 γ_u 의 조합으로 표현한다(식 3.2). 핵심 정리 3.2는 각 b_{θ_j}(h) 가 |h|^{-2r_j} 꼴의 상한을 갖는 경우, 평균값의 λ‑모멘트(λ∈(1/(2r),1])가 γ_u와 ζ(2r_jλ) 의 곱으로 제한됨을 보인다. 이를 통해 임의의 z 에 대한 평균 최악오차가 O(N^{-r}) 수렴을 보이며, r>½는 θ_j 의 스무스니스 파라미터에 의해 결정된다.

중위값 추정기의 확률적 경계는 정리 3.4와 그 뒤의 코롤라리 3.3을 통해 도출된다. 구체적으로, ε>0와 ρ∈(0,1) 에 대해
P{ |M_k(f)-I_ϕ(f)| ≥ c ρ^{-½−r+ε} N^{-r+ε} } ≤ ρ^{k+1/2}/4
가 성립한다. 여기서 k=O(log N)이면 위 확률이 N에 대해 초지수적으로 감소한다. 평균 절대 오차(Mean Absolute Error, MAE)는 이 확률 경계와 L^2_ϕ 에 대한 포함 관계를 이용해
E