[논문 리뷰] Linear multi-step schemes for BSDEs

Jean-François Chassagneux|arXiv (Cornell University)|2013. 06. 24.

Numerical methods for differential equations참고 문헌 2인용 수 6

한 줄 요약

이 논문은 안정성 분석과 절삭 오차 분석을 활용하여 임의의 수렴 차수를 가진 후행 확률 미분 방정식(BSDEs)을 해결하기 위한 선형 다단계 스킴을 제안한다. 오일러 형 방법을 초월하는 고차수 수치적 방법을 개발하고, 고전적 선형 다단계 방법(예: 애덤스, 니스트롬, 밀린)을 일반적인 확산 계수와 Z에 의존하는 드라이버를 갖는 BSDEs로 확장하며, 시간 이산화에서 모멘트 매칭을 통해 수렴 속도를 m차수까지 증명한다.

ABSTRACT

We study the convergence rate of a class of linear multi-step methods for BSDEs. We show that, under a sufficient condition on the coefficients, the schemes enjoy a fundamental stability property. Coupling this result to an analysis of the truncation error allows us to design approximation with arbitrary order of convergence. Contrary to the analysis performed in \cite{zhazha10}, we consider general diffusion model and BSDEs with driver depending on $z$. The class of methods we consider contains well known methods from the ODE framework as Nystrom, Milne or Adams methods. We also study a class of Predictor-Correctot methods based on Adams methods. Finally, we provide a numerical illustration of the convergence of some methods.

연구 동기 및 목표

• 오일러 형 방법을 초월하는 고차수 수치적 스킴을 BSDEs에 개발하여 수렴 속도의 한계를 해결한다.
• 고전적 선형 다단계 방법(예: 애덤스, 밀린, 니스트롬)을 일반적인 확산 계수와 Z에 의존하는 드라이버를 갖는 BSDEs의 확률적 설정으로 확장한다.
• 후행 SDEs의 맥락에서 안정성과 절삭 오차 분석을 통해 수렴의 이론적 기반을 확립한다.
• 정확도와 안정성을 향상시키기 위해 애덤스 스킴을 기반으로 한 예측-수정 방법을 설계한다.
• 부드러움과 리프시츠 조건 하에서 엄밀한 수렴 분석을 제공하며, 시간 간격과 모멘트 매칭 차수에 따라 명시적인 오차 한계를 제시한다.

제안 방법

• r단계 초기화를 갖는 후행 재귀를 사용하여 BSDEs에 대한 선형 다단계 스킴의 클래스를 수립하며, 계수 $ a_j, b_{i,j}, \alpha_j, \beta_{i,j} $ 가 사전 일致성 및 유계 조건을 만족하도록 정의한다.
• 필터레이션 $ \mathcal{F}_{t_{i+j}} $ 에 적응하는 랜덤 가중치 $ H^{Y}_{i,j}, H^{f}_{i,j} $ 를 도입하여 조건부 기대값이 0이고 2차 모멘트가 유계임을 보장한다.
• 시간 격자 $ \pi $ 에서 $ X_t $ 의 전방 과정을 근사하기 위해 모멘트를 순서 $ K $ 까지 매칭하는 이산 마코프 과정 $ \widehat{W}^{x,0}_{t} $ 을 구성한다.
• 조건부 기대값을 사용하여 스킴의 예측 및 수정 단계를 정의한다: $ Y_i = \mathbb{E}_{t_i}[\sum a_j Y_{i+j} + h \sum b_{i,j} f(Y_{i+j}, Z_{i+j})] $, 마찬가지로 $ Z_i $ 에 대해서도 정의한다.
• 정확한 해 $ (Y_t, Z_t) $ 와 이산 스킴 $ (\widehat{Y}_i, \widehat{Z}_i) $ 를 비교하기 위해 펌터베이션 분석을 적용하며, 오차를 시간 이산화 오차와 공간 이산화 오차 성분으로 분해한다.
• 와이너 증분 $ W^{0,t_i}_{t_{i+j}} $ 와 그 이산 대체물 $ \widehat{W}^{0,t_i}_{t_{i+j}} $ 간의 모멘트 매칭을 통해 공간 이산화 오차를 순서 $ K $ 까지 제어한다.

실험 결과

연구 질문

• RQ1 일반적인 확산 계수와 Z에 의존하는 드라이버를 갖는 BSDEs에 대해, ODEs의 선형 다단계 방법을 고차수 수렴을 유지하면서 적응시킬 수 있는가?
• RQ2 리프시츠 조건 하에서 선형 다단계 스킴의 수렴을 보장하는 안정성 조건은 무엇인가?
• RQ3 시간 이산화에서의 모멘트 매칭 차수는 스킴의 전반적 수렴 속도에 어떤 영향을 미치는가?
• RQ4 제안된 스킴에서 시간 이산화 오차와 공간 이산화 오차 간의 상호작용은 어떠한가?
• RQ5 애덤스 방법을 기반으로 한 예측-수정 변형 스킴은 BSDE 근사에서 더 높은 정확도와 안정성을 달성할 수 있는가?

주요 결과

• 제안된 선형 다단계 스킴은 모멘트 매칭 차수 $ K $ 가 $ K = 2m + 1 $ 를 만족할 경우, 임의의 수렴 차수 $ m $ 를 달성한다.
• 계수에 대한 충분한 조건이 만족되면, 스킴는 수렴 분석을 가능하게 하는 기본적인 안정성 성질을 갖는다.
• 시간 이산화 오차는 차수 $ m $ 이며, 공간 이산화 오차는 $ C|\pi|^{(K-1)/2} $ 로 유계이므로, 전반적 오차 한계는 $ |\widetilde{Y}_0 - \widehat{Y}_0| \leq C|\pi|^m + C|\pi|^{(K-1)/2} $ 로 나타난다.
• 만약 $ K = 2m + 1 $ 이면, 전반적 오차는 차수 $ m $ 으로 수렴하며, 이는 스킴가 원하는 임의의 수렴 차수를 달성할 수 있음을 확인한다.
• 이 방법은 애덤스, 니스트롬, 밀린와 같은 잘 알려진 ODE 스킴을 특수 케이스로 포함하며, 이들의 적용 범위를 BSDEs로 확장한다.
• 수치적 예시는 이론적 수렴 속도를 확인하며, 제안된 스킴의 고차수 성능을 검증한다.