[논문 리뷰] Error estimates on ergodic properties of discretized Feynman-Kac semigroups

Grégoire Ferré, Gabriel Stoltz|arXiv (Cornell University)|2017. 12. 11.

Stochastic processes and financial applications참고 문헌 71인용 수 13

한 줄 요약

이 논문은 확산 과정과 관련된 파인만–카프 반군의 시간 이산화에 대한 엄밀한 오차 추정을 제공하며, 생성자에 대한 불변 측도와 주 eigen값의 편향을 분석한다. 탈라예–투바로 유형의 분석을 통해 컴팩트 다양체 위의 이산화 동역학에 대한 수렴 속도를 확립하며, 확산 몽테카를로 및 관련 분야의 수치적 방법에 대한 이론적 근거를 제시한다.

ABSTRACT

We consider the numerical analysis of the time discretization of Feynman-Kac semigroups associated with diffusion processes. These semigroups naturally appear in several fields, such as large deviation theory, Diffusion Monte Carlo or non-linear filtering. We present errors estimates a la Talay-Tubaro on their invariant measures when the underlying continuous stochastic differential equation is discretized; as well as on the leading eigenvalue of the generator of the dynamics, which corresponds to the rate of creation of probability. This provides criteria to construct efficient integration schemes of Feynman-Kac dynamics, as well as a mathematical justification of numerical results already observed in the Diffusion Monte Carlo community. Our analysis is illustrated by numerical simulations.

연구 동기 및 목표

• 확산 과정과 관련된 파인만–카프 반군의 시간 이산화가 유도하는 편향을 분석하기 위해.
• 시간 이산화 하에서 생성자의 불변 측도와 주 고유값에 대한 오차를 정량화하기 위해.
• 확산 몽테카를로 및 비선형 필터링에 사용되는 수치적 방법에 대한 이론적 근거를 제공하기 위해.
• 재가중치를 포함한 비선형, 측도 기반 파인만–카프 동역학에 탈라예–투바로 오차 추정을 확장하기 위해.
• 콤���트 영역에서 파인만–카프 동역학의 효율적 수치적 통합을 위한 수렴 기준을 설정하기 위해.

제안 방법

• 가우시안 노이즈가 첨가된 확률적 미분 방정식과 가중치 함수 $W$를 통해 정의된 연속 시간 파인만–카프 반군을 분석한다.
• 유한한 시간 간격 $\Delta t$를 사용하여 SDE를 이산화하고, 시간 적분을 구적법으로 근사한다.
• 결과적으로 얻어진 이산 시간 마르코프 체인 $q^n$과 그에 대응하는 반군 $\Phi_{\Delta t,n}^W$를 연구한다.
• 탈라예–투바로 유형의 오차 분석을 적용하여 불변 측도 $\nu_{W,\Delta t}$가 연속 극한 $\nu_W$로 수렴하는 수렴 속도를 유도한다.
• 도브루시닌의 에르고딕 계수와 수축 노름을 사용하여 이산 반군의 혼합 성질과 수렴 속도를 근사한다.
• 기대값의 비율에 대한 균일한 상한을 통해 주 고유값 $\lambda_{\Delta t} = \tilde{\lambda}_{\Delta t}$의 등식을 증명하여 고유값 근사의 일致성을 입증한다.

실험 결과

연구 질문

• RQ1 시간 이산화가 파인만–카프 반군의 불변 측도에 어떤 영향을 미치는가?
• RQ2 이산화된 불변 측도가 진짜 불변 측도로 수렴하는 속도는 얼마인가?
• RQ3 이산화 오차가 확률 생성 속도를 결정하는 생성자의 주 고유값에 어떤 영향을 미치는가?
• RQ4 탈라예–투바로 프레임워크는 파인만–카프 과정과 같은 비선형, 측도 기반 반군에 어떻게 적용될 수 있는가?
• RQ5 이산화된 반군이 연속 동역학의 에르고딕 성질을 유지하기 위한 조건은 무엇인가?

주요 결과

• 논문은 이산화된 파인만–카프 반군의 불변 측도 $\nu_{W,\Delta t}$가 $\Delta t \to 0$일 때 진짜 불변 측도 $\nu_W$로 수렴하며, 탈라예–투바로 분석을 통해 이에 대한 오차 한계를 도출함을 입증한다.
• 이산화된 생성자의 주 고유값 $\lambda_{\Delta t}$는 진짜 고유값 $\lambda$로 수렴하며, 기대값의 비율에 대한 균일한 상한을 통해 $\lambda_{\Delta t} = \tilde{\lambda}_{\Delta t}$가 증명되어 고유값 근사의 일致성을 확보한다.
• 이산 반군의 수축 노름은 $1 - \alpha(Q)$로 유계이며, 여기서 $\alpha(Q)$는 도브루시닌의 에르고딕 계수이므로, 불변 측도로의 지수 수렴이 보장된다.
• 불변 측도의 오차는 벡터장 $b$, 가중치 함수 $W$, 그리고 이산화 방법의 약한 순서에 의해 제어된다.
• 수치 시뮬레이션은 이론적 오차 추정을 확인하고 실용적 환경에서의 수렴 속도를 검증한다.
• 분석을 통해 확산 몽테카를로에서 지상 상태 에너지 추정의 편향을 줄이기 위해 고차수 이산화 방법을 사용하는 것이 타당함을 입증한다.