[논문 리뷰] Multi-level Picard approximations of high-dimensional semilinear parabolic differential equations with gradient-dependent nonlinearities

Martin Hutzenthaler, Thomas Kruse|arXiv (Cornell University)|2017. 11. 03.

Stochastic processes and financial applications참고 문헌 10인용 수 41

한 줄 요약

이 논문은 고차원 비선형 포아송-파라볼릭 PDE를 그라디언트에 의존하는 비선형성을 가진 경우에 대해 다중 수준 피카르 반복법과 가우스-레지앙드르 정적분을 조합하여 근사하는 방법을 제안한다. 적절한 정규성 조건 하에서 계산 복잡도가 차원 d와 정밀도 역수 ε에 대해 다항식적으로 증가함을 증명하였으며, 특히 임의의 δ > 0에 대해 O(dε^{-(4+δ)})로 유계임을 보였다. 이는 고차원 PDE 근사 분야에서 오랫동안 남아 있던 열린 문제를 해결한 것이다.

ABSTRACT

Parabolic partial differential equations (PDEs) and backward stochastic differential equations (BSDEs) have a wide range of applications. In particular, high-dimensional PDEs with gradient-dependent nonlinearities appear often in the state-of-the-art pricing and hedging of financial derivatives. In this article we prove that semilinear heat equations with gradient-dependent nonlinearities can be approximated under suitable assumptions with computational complexity that grows polynomially both in the dimension and the reciprocal of the accuracy.

연구 동기 및 목표

• 고차원 비선형 포아송-파라볼릭 PDE의 근사에 있어 계산적 과제를 해결하기 위해.
• 이러한 PDE가 차원 d와 정밀도 역수 ε에 대해 계산 복잡도가 다항식적으로 증가하는지 여부라는 열린 문제를 해결하기 위해.
• 기존의 피카르 기반 방법을 확장하여 이전 방법이 일반적인 종단 조건 하에서 다룰 수 없었던 그라디언트에 의존하는 비선형성을 처리할 수 있도록 하기 위해.
• 제안된 방법에 대해 엄밀한 복잡도 상한을 확립하여 고차원에서의 확장 가능성을 보장하기 위해.

제안 방법

• 비선형 포아송-파라볼릭 PDE의 해를 반복적으로 근사하기 위해 다중 수준 피카르 반복을 적용한다.
• 피카르 반복 과정에서 시간 적분을 효율적으로 근사하기 위해 가우스-레지앙드르 정적분 규칙을 통합하며, 부드러운 피적분 함수에 대해 스펙트럼 수렴 특성을 활용한다.
• 다중 수준 구조의 각 수준에서의 전반적 오차를 제어하기 위해 재귀적 오차 상한(부등식 54)을 사용한다.
• 그라디언트에 의존하는 비선형성 존재 하에서 오차 전파를 관리하기 위해 이산 그론발라 유형 부등식과 준노름 추정을 적용한다.
• 재귀적 오차 추정을 반복하여 비재귀적 전반적 오차 상한(부등식 65)을 유도한다.
• 후행 확률 미분 방정식(스모드-디피에이션 방정식, BSDE)을 통한 PDE 해의 확률적 표현과 몬테카를로 샘플링을 활용하여 수치적 평가를 수행한다.

실험 결과

연구 질문

• RQ1 고차원 비선형 포아송-파라볼릭 PDE가 그라디언트에 의존하는 비선형성을 가질 경우, 계산 복잡도가 차원 d와 정밀도 역수 ε에 대해 다항식적으로 증가하는 방식으로 근사 가능한가?
• RQ2 가우스-레지앙드르 정적분을 사용하는 다중 수준 피카르 방법이 비선형성이 해의 그라디언트에 의존할 경우에도 다항식 복잡도를 유지하는가?
• RQ3 일반적인 해의 정규성 가정 하에서 제안된 방법의 계산 비용이 임의의 δ > 0에 대해 O(dε^{-(4+δ)})로 유계가 되는가?
• RQ4 지난 분기 확산 방법과 달리, 종단 조건이 작을 필요 없이 일반적인 종단 조건을 다룰 수 있는가?
• RQ5 그라디언트에 의존하는 경우에 있어 계산 복잡도가 차원 d와 정밀도 ε에 대해 정확히 어떻게 의존하는가?

주요 결과

• 다중 수준 피카르 근사의 계산 복잡도는 임의의 δ > 0에 대해 O(dε^{-(4+δ)})로 유계이며, 여기서 d는 차원이고 ε는 지정된 정밀도이다.
• 정확한 해에 대해 적절한 정규성 조건을 만족할 경우, 이 복잡도 상한을 달성한다. 이는 해 u∞의 C∞ 연속성 포함한다.
• 가우스-레지앙드르 정적분의 사용으로 시간 적분의 수렴 속도가 빨라지며, 이는 낮은 복잡도를 유지하는 데 핵심적이다.
• 오차 분석을 통해 비재귀적 전반적 오차 상한(부등식 65)을 확립하였으며, 이는 복잡도 추정의 기초가 된다.
• 이 방법은 일반적인 종단 조건에 적용 가능하며, 비선형성이 작을 필요가 없어 이전의 분기 확산 방법과는 다르다.
• 결과적으로, 제안된 다중 수준 피카르 접근법을 통해 이 클래스의 PDE에 대해 차원의 저주가 극복됨을 확인한다.