[논문 리뷰] Cubature methods to solve BSDEs: Error expansion and complexity control

Jean-François Chassagneux, Camilo Andrés García Trillos|arXiv (Cornell University)|2017. 02. 03.

Simulation Techniques and Applications참고 문헌 30인용 수 3

한 줄 요약

이 논문은 위너 공간 위의 큐버처를 사용하여 마코프성 백분산 확률미분방정식(BSDEs)를 풀기 위한 고차수 수치적 방법을 제안한다. 이 방법은 Richardson-Romberg 보간을 통한 수렴 가속화를 가능하게 하는 명시적 오차 전개를 포함한다. 희소 격자 보간 투영을 도입함으로써 계산 복잡도를 크게 감소시키면서도 정확도를 유지하며, 특히 유한 격자 위에서의 선형 보간을 통해 뛰어난 성능을 발휘한다.

ABSTRACT

We obtain an explicit error expansion for the solution of Backward Stochastic Differential Equations (BSDEs) using the cubature on Wiener spaces method. The result is proved under a mild strengthening of the assumptions needed for the application of the cubature method. The explicit expansion can then be used to construct implementable higher order approximations via Richardson-Romberg extrapolation. To allow for an effective efficiency improvement of the interpolated algorithm, we introduce an additional projection on finite grids through interpolation operators. We study the resulting complexity reduction in the case of the linear interpolation.

연구 동기 및 목표

• 마코프성 BSDEs에 큐버처 방법을 적용할 때, 약간의 정규성 조건을 강화한 가정 하에서 명시적 오차 전개를 유도하는 것.
• 유도된 오차 전개를 활용해 Richardson-Romberg 보간을 통해 BSDE 근사의 고차수 수렴을 달성하는 것.
• 보간 연산자를 통한 유한 격자 투영을 도입하여 큐버처 기반 BSDE 해법기의 계산 복잡도를 감소시키는 것.
• 특히 선형 및 희소 격자 보간을 사용할 때, 해법의 복잡도-성능 트레이드오���을 분석하고 최적화하는 것.
• BSDE 및 관련 비선형 편미분방정식의 수치적 해법에서 이론적 및 실용적 효율성 향상을 확립하는 것.

제안 방법

• 가정 1.2에 따라 정규성 조건을 강화한 상황에서, 위너 공간 위의 큐버처 방법을 사용한 BSDE의 약한 근사에 대한 명시적 오차 전개를 유도한다.
• 오차 전개에 Richardson-Romberg 보간을 적용하여, 기초 이산화 단계를 늘리지 않고도 고차수 수렴 스킴을 구성한다.
• 해결책을 유한 격자에 투영하기 위해 보간 연산자를 도입하여 계산 비용을 감소시키면서도 정확도를 유지한다.
• 유한 격자 위에서 선형 및 희소 격자 보간(순서 $p$의 희소 격자 노드 공간을 통해)을 사용하여 근사화에 필요한 격점 수를 최소화한다.
• 큐버처 측도 $\hat{\mathbb{Q}}$ 하에서 조건부 기대값을 사용하여 격점 $D_i$에서 근사해 $\check{u}_i(x)$ 및 $\check{v}_i(x)$를 계산하기 위해 후행 스킴을 수정한다.
• 다음과 같은 재귀적 갱신을 통해 희소 큐버처 후행 스킴을 정의한다: $\check{u}_i(x) = \pi_{\mathcal{V}_{p_i - d + 1}}^{D_i}[ \mathbb{E}^{\hat{\mathbb{Q}}}[\check{u}_{i+1}(\hat{X}^{t_i,x}_{t_{i+1}})] + h_i f(\cdot, \check{u}_i, \check{v}_i) ](x)$, $\check{v}_i(x)$ 는 마팅게일 차분 항을 통해 계산된다.

실험 결과

연구 질문

• RQ1 마코프성 BSDEs에 큐버처 방법을 적용할 때, 정규성 조건을 약간 강화한 상황에서 명시적 오차 전개를 도출할 수 있는가?
• RQ2 유도된 오차 전개를 어떻게 활용하여 Richardson-Romberg 보간을 통해 고차수 수렴 스킴을 구성할 수 있는가?
• RQ3 유한 격자에 해를 보간 투영하는 것이 큐버처 기반 BSDE 해법기의 전체 복잡도와 정확도에 어떤 영향을 미치는가?
• RQ4 전체 격자 보간 대비 희소 격자 보간을 사용할 때 계산 복잡도는 어떻게 감소하는가? 동시에 수렴 속도는 유지되는가?
• RQ5 큐버처 방법과 보간 기반 격자 투영을 결합함으로써 BSDE 근사에서 이론적 및 실용적 복잡도 감소는 어느 정도 이루어지는가?

주요 결과

• 가정 1.2 하에서 위너 공간 위의 큐버처 방법에 대한 명시적 오차 전개를 도출하였으며, 이는 Richardson-Romberg 보간을 통한 고차수 수렴 가능성을 보장한다.
• 오차 전개는 계수, 시험 함수, 생성자 구조에 대한 의존성을 드러내며, 체계적인 수렴 가속화의 기초를 제공한다.
• 특히 선형 및 희소 격자 보간을 사용한 유한 격자 위의 보간은 정확도를 떨어뜨리지 않으면서도 필요한 평가 점 수를 크게 감소시킨다.
• 희소 격자 구현은 $d$차원 문제에서 전체 격자 복잡도 $\mathcal{O}(n^d)$ 대비 복잡도 감소 비율 $\mathcal{O}(n^{d-1})$을 달성하며, 오차 저하도 최소화된다.
• 이론적 분석에 따르면, 희소 격자 보간 오차는 $\mathcal{O}(n^{-\gamma(\beta^* - \beta + 1) + 1})$ 으로 제어되며, 여기서 $\gamma$ 및 $\beta$ 는 해와 시험 함수의 정규성과 관련된다.
• 최종 알고리즘은 $\cup_{i=0}^{n-1}\{t_i\} \times D_i^{\blacksquare}$ 와 같은 더 큰 영역에서 정의되며, 재계산 없이도 중간 점에서의 효율적 평가를 가능하게 한다.