Skip to main content
QUICK REVIEW

[논문 리뷰] Generalised multilevel Picard approximations

Michael B. Giles, Arnulf Jentzen|arXiv (Cornell University)|2019. 11. 08.
Advanced Numerical Methods in Computational Mathematics참고 문헌 95인용 수 13
한 줄 요약

이 논문은 고차원 비선형 포아송-파라볼릭 PDE를 해결하기 위한 일반화된 전기역사 재귀 다수준 피카르(MLP) 근사 방법의 클래스를 제안한다. 이 방법들이 차원의 고통을 극복함을 증명하며, 계산 비용이 차원 $d$와 $ε^{-1}$에 대해 다항식적으로 증가함을 보이며, 일반적인 시간 구간과 초기 조건에서 MLP 방법에 대해 처음으로 체계적인 복잡도 결과를 확립한다.

ABSTRACT

It is one of the most challenging problems in applied mathematics to approximatively solve high-dimensional partial differential equations (PDEs). In particular, most of the numerical approximation schemes studied in the scientific literature suffer under the curse of dimensionality in the sense that the number of computational operations needed to compute an approximation with an error of size at most $ \varepsilon > 0 $ grows at least exponentially in the PDE dimension $ d \in \mathbb{N} $ or in the reciprocal of $ \varepsilon $. Recently, so-called full-history recursive multilevel Picard (MLP) approximation methods have been introduced to tackle the problem of approximately solving high-dimensional PDEs. MLP approximation methods currently are, to the best of our knowledge, the only methods for parabolic semi-linear PDEs with general time horizons and general initial conditions for which there is a rigorous proof that they are indeed able to beat the curse of dimensionality. The main purpose of this work is to investigate MLP approximation methods in more depth, to reveal more clearly how these methods can overcome the curse of dimensionality, and to propose a generalised class of MLP approximation schemes, which covers previously analysed MLP approximation schemes as special cases. In particular, we develop an abstract framework in which this class of generalised MLP approximations can be formulated and analysed and, thereafter, apply this abstract framework to derive a computational complexity result for suitable MLP approximations for semi-linear heat equations. These resulting MLP approximations for semi-linear heat equations essentially are generalisations of previously introduced MLP approximations for semi-linear heat equations.

연구 동기 및 목표

  • 고차원 비선형 포아송-파라볼릭 PDE를 해결하는 데 있어 고전적 방법이 차원의 고통으로 인해 계산적으로 비가능한 문제에 대응한다.
  • 유한 차분법이나 딥러닝 방법과 같은 기존 수치적 방법의 한계를 극복한다. 이들 방법은 기하급수적인 비용 증가를 겪거나 엄밀한 오차 분석이 부족하다.
  • 기존의 MLP 방법을 통합하고 확장할 수 있는 전기역사 재귀 다수준 피카르(MLP) 근사의 일반화된 프레임워크를 개발한다.
  • 일반적인 초기 조건과 시간 구간에서 제안된 일반화된 MLP 근사 방법이 차원의 고통을 극복할 수 있음을 엄밀히 증명하는 복잡도 분석을 제공한다.
  • 차원 $d$와 역오차 허용 오차 $ε^{-1}$에 대해 다항식적 의존성을 보여주는 계산 복잡도 상한을 수립한다.

제안 방법

  • 힐베르트 공간에서의 편향-분산 분해와 일반화된 분산 항등식을 사용하여 일반화된 MLP 근사의 추상적 프레임워크를 수립한다.
  • 모든 수준에서 브라운 운동 경로의 몬테카를로 샘플링을 통해 해를 정밀화하는 전기역사 재귀 기반의 반복적 근사 방법을 정의한다.
  • 역방향 확률미분방정식(BSDE)을 통한 PDE 해의 확률적 표현을 사용하고, 브라운 운동의 함수형을 반복적으로 근사하여 해를 도출한다.
  • 경로 시뮬레이션과 기대값 계산을 위해 독립적인 랜덤 변수 $U^\theta$와 브라운 운동 $W^{d,\theta}$를 포함하는 재귀적 공식을 도입한다.
  • 각 수준당 몬테카를로 샘플 수와 계산 연산 수를 분석하여 비용 상한을 유도하며, 비용이 $ (M_j)^n $의 형태로 증가함을 보인다. 여기서 $M_j$는 샘플 수를 제어한다.
  • 재귀적 방법이 잘 정의되고 수렴하기 위해 필요한 적분 가능성과 가측성 조건을 확립한다.

실험 결과

연구 질문

  • RQ1고차원 PDE에 대한 기존의 MLP 방법을 통합하고 확장할 수 있는 일반화된 다수준 피카르 근사 방법의 클래스를 구성할 수 있는가?
  • RQ2일반적인 초기 조건과 시간 구간에서 비선형 열 방정식의 맥락에서 제안된 일반화된 MLP 방법이 차원의 고통을 어떻게 극복하는가?
  • RQ3차원 $d$와 오차 허용 오차 $ε$에 대해 일반화된 MLP 근사 방법의 계산 복잡도는 어떻게 되는가?
  • RQ4초기 자료와 비선형성에 대해 어떤 조건에서 일반화된 MLP 방법이 다항식 비용 증가를 유지하는가?
  • RQ5추상적 프레임워크를 구체화하여 비선형 열 방정식에 대해 명시적인 복잡도 상한을 도출할 수 있는가?

주요 결과

  • 제안된 일반화된 MLP 근사 방법은 일반적인 시간 구간과 초기 조건에서 포아송-비선형 PDE에 대해 차원의 고통을 극복함을 엄밀히 증명한 최초의 방법이다.
  • 비선형 열 방정식의 경우, MLP 근사의 계산 비용은 $ \mathrm{Cost}_{d,N_{d,\varepsilon},N_{d,\varepsilon}} \leq C_\delta \, d^{1+p(1+\delta)} \varepsilon^{-2(1+\delta)} $ 를 만족하며, 이는 $d$와 $ε^{-1}$에 대해 다항식적 증가를 증명한다.
  • 근사의 오차는 $ \mathbb{E}[|Y_{n,n}^{d,0}(0,\xi_d) - y_d(0,\xi_d)|^2]^{1/2} \leq \varepsilon $ 를 만족하며, $n \geq N_{d,\varepsilon}$ 일 때 진짜 해로의 수렴과 제어 가능한 오차를 보장한다.
  • 이 프레임워크는 일반적인 초기 조건 $\xi_d$와 리프시츠 비선형성 $f$를 허용하며, $g_d$의 성장 조건이 해와 그 모멘트의 유계성을 보장한다.
  • 분석을 통해 비용이 $d$와 $ε^{-1}$에 대해 최대 다항식적으로 증가함을 입증하였으며, 지수는 종단 조건과 비선형성의 다항식 성장률 $p$에 따라 달라진다.
  • 초기 자료와 계수에 대해 온건한 적분 가능성과 가측성 가정 하에 결과가 성립하며, 상수는 $L$, $T$, $K$, $p$에 따라 명시적으로 결정된다.

더 나은 연구,지금 바로 시작하세요

논문 읽기부터 검토까지, 연구 시간을 획기적으로 줄여보세요.

카드 등록 없음 · 무료 플랜 제공

이 리뷰는 AI가 만들고, 인간 에디터가 검토했습니다.