[논문 리뷰] A Milstein Scheme for SPDEs

Arnulf Jentzen, Michael Roeckner|arXiv (Cornell University)|2010. 01. 15.

Stochastic processes and financial applications인용 수 6

한 줄 요약

이 논문은 흔적-클래스 노이즈를 가진 반선형 스토크라스틱 편미분방정식(SPDEs)에 대해 미르스타인의 수치적 방법의 무한차원 버전을 제안한다. 이는 특정 클래스의 반선형 SPDEs에서 성립하는 교환 조건을 활용한 것으로, 이로 인해 기존 알고리즘 대비 더 낮은 계산 비용과 더 적은 난수 변수를 사용하여 더 효율적인 시뮬레이션을 가능하게 한다. 수치 결과는 스토크라스틱 열 방정식과 반응-확산 방정식에 대해 이와 같은 성능 향상을 입증한다.

ABSTRACT

This article studies an infinite dimensional analog of Milstein's scheme for finite dimensional stochastic ordinary differential equations (SODEs). The Milstein scheme is known to be impressively efficient for SODEs which fulfill a certain commutativity type condition. This article introduces the infinite dimensional analog of this commutativity type condition and observes that a certain class of semilinear stochastic partial differential equation (SPDEs) with multiplicative trace class noise naturally fulfills the resulting infinite dimensional commutativity condition. In particular, a suitable infinite dimensional analog of Milstein's algorithm can be simulated efficiently for such SPDEs and requires less computational operations and random variables than previously considered algorithms for simulating such SPDEs. The analysis is supported by numerical results for a stochastic heat equation and stochastic reaction diffusion equations showing signifficant computational savings.

연구 동기 및 목표

• 유한차원에서의 밀스타인 방법을 무한차원 SPDEs로 확장한다.
• 고차수 약한 수렴을 보장하는 무한차원 교환 조건의 무한차원 해석을 규명한다.
• 다중성 흔적-클래스 노이즈를 가진 SPDEs를 시뮬레이션하기 위한 수치적으로 효율적인 알고리즘을 개발한다.
• 실제 SPDE 시뮬레이션에서 기존 방법 대비 뚜렷한 계산 비용 절감 효과를 입증한다.

제안 방법

• 유한차원에서의 교환 조건을 힐버트 공간에 값이 있는 SDEs로 확장하여 무한차원 밀스타인 방법의 버전을 유도한다.
• 무한차원 설정에서, 반복 스토크라스틱 적분이 필요 없이 강수렴을 달성할 수 있도록 보장하는 교환 조건을 정의한다.
• 자연스럽게 제안된 교환 조건을 만족하는 흔적-클래스 노이즈에 의해 구동되는 반선형 SPDEs에 이 방법을 적용한다.
• 무한차원 상태 공간을 다루기 위해 스펙트럼 콜로케이션 또는 갈레르킨 근사법을 사용하여 알고리즘을 구현한다.
• 스토크라스틱 열 방정식과 스토크라스틱 반응-확산 방정식에 대해 수치적 검증을 수행한다.
• 표준 옆러-마루야마 방법 및 기타 기존 방법과의 비교를 통해 연산 횟수와 난수 변수 수 측면에서의 계산 효율성을 평가한다.

실험 결과

연구 질문

• RQ1 밀스타인의 고차수 수치 방법은 다중성 노이즈를 가진 무한차원 SPDEs로 일반화될 수 있는가?
• RQ2 유한차원 SODEs에서 관찰된 효율성 향상과 동일한 성과를 달성하기 위해 필요한 무한차원 조건은 무엇인가?
• RQ3 흔적-클래스 노이즈를 가진 SPDEs는 이 일반화된 교환 조건을 만족하는가?
• RQ4 제안된 방법은 SPDEs에 대한 기존 수치적 방법보다 계산 비용이 얼마나 줄어들 수 있는가?
• RQ5 대표적인 SPDE 모델에 대해 실질적으로 얼마 정도의 계산 비용 절감 효과를 기대할 수 있는가?

주요 결과

• 제안된 무한차원 밀스타인 방법은 흔적-클래스 노이즈를 가진 SPDEs에 대해 표준 옆러 방법보다 더 높은 약한 수렴 차수를 달성한다.
• 교환 조건은 스토크라스틱 열 방정식과 반응-확산 방정식을 포함한 광범위한 반선형 SPDEs 클래스에서 자연스럽게 만족된다.
• 이전 알고리즘 대비 더 적은 난수 변수와 계산 연산을 요구하여 뚜렷한 효율성 향상을 이룬다.
• 수치 실험 결과는 이론적 예측을 확인하며, 스토크라스틱 열 방정식과 반응-확산 시스템 모두에서 계산 비용의 측정 가능한 감소를 보여준다.
• 반복 스토크라스틱 적분이 필요 없이도 안정적이고 정확한 SPDE 시뮬레이션을 가능하게 한다.