[논문 리뷰] Efficient Solution of Backward Jump-Diffusion PIDEs with Splitting and Matrix Exponentials

Andrey Itkin|arXiv (Cornell University)|2013. 04. 10.

Differential Equations and Numerical Methods인용 수 4

한 줄 요약

이 논문은 수학적 금융에서 발생하는 후행 점프-확산 PIDE를 해결하기 위한 새로운 조건부 안정성, 2차 정확도를 갖춘 수치적 방법을 제안한다. 이 방법은 확산 및 점프 과정에 대해 연산자 분할을 적용하고, 행렬 지수 또는 패드 근사법을 사용한다. 접근 방식은 비국소적 점프 적분을 국소적 의사미분 연산자로 변환하여, 겹치는 격자 위에서 유한차분 스킴을 통해 효율적이고 양수성 보존 해를 도출한다. 이는 복잡한 모델인 헤스턴+베이글(GVG)과 같은 경우에 FFT 기반 대안보다 뛰어난 수렴성과 성능을 보여준다.

ABSTRACT

We propose a new, unified approach to solving jump-diffusion partial integro-differential equations (PIDEs) that often appear in mathematical finance. Our method consists of the following steps. First, a second-order operator splitting on financial processes (diffusion and jumps) is applied to these PIDEs. To solve the diffusion equation, we use standard finite-difference methods, which for multi-dimensional problems could also include splitting on various dimensions. For the jump part, we transform the jump integral into a pseudo-differential operator. Then for various jump models we show how to construct an appropriate first and second order approximation on a grid which supersets the grid that we used for the diffusion part. These approximations make the scheme to be unconditionally stable in time and preserve positivity of the solution which is computed either via a matrix exponential, or via P{á}de approximation of the matrix exponent. Various numerical experiments are provided to justify these results.

연구 동기 및 목표

• 지속적 확률적 변동성에 점프가 포함된 복잡한 금융 모델에서 발생하는 후행 점프-확산 PIDE를 해결하는 데 있어 계산적 과제를 해결하기 위해.
• 옵션 가격의 양수성을 유지하는 것이 금융 모델링에 매우 중요하므로, 수치적으로 안정적이고 2차 정확도를 갖춘 스킴을 개발하기 위해.
• 메르톤, 쿠, CGMY 등의 다양한 점프 모델을 통합적으로 다룰 수 있도록 비국소적 점프 적분을 국소적 비선형 미분 연산자로 변환하기 위해.
• FFT 기반 방법에서 흔히 발생하는 느린 피카르 반복과 워핑-아웃 오류를 피하기 위해, 행렬 지수 또는 패드 근사법을 활용해 효율적인 계산을 가능하게 하기 위해.
• 특성 함수가 도달 불가능한 모델, 예를 들어 헤스턴에 베이글 점프가 추가된 경우와 같이 변환 방법이 비현실적인 모델로의 적용 가능성을 넓히기 위해.

제안 방법

• PIDE의 확산 및 점프 성분을 분리하기 위해 2차 연산자 분할을 적용하여 시간에 따라 순차적으로 처리한다.
• 확산 부분은 표준 유한차분 스킴을 사용하여 해결하며, 다차원 문제의 경우 차원 분할을 적용한다.
• 점프 적분은 변수 변경을 통해 의사미분 연산자로 변환되어, 확산 격자 외부로 연장된 격자에서 국소적으로 근사할 수 있다.
• 각 점프 모델(Merton, Kou, CGMY)에 대해 1차 및 2차 유한차분 근사를 구성하여 조건부 안정성과 양수성 보존을 보장한다.
• 행렬 지수 또는 행렬 지수의 패드 유리 근사법을 사용하여 시간에 따라 해를 진행시키며, 이는 조건부 안정성과 양수성 보존을 유지한다.
• 핵심 혁신은 겹치는 격자와 보간을 활용하여 비균일 격자에서도 2차 정확도를 유지하면서 FFT 관련 워핑-아웃 오류를 방지하는 것이다.

실험 결과

연구 질문

• RQ1 일반 레비 점프 동역학을 갖는 후행 점프-확산 PIDE에 대해 2차 정확도, 조건부 안정성, 양수성 보존을 동시에 확보한 수치적 스킴을 구성할 수 있는가?
• RQ2 점프 격자와 확산 격자가 일치하지 않을 경우, 특히 FFT에 의존하지 않고 비국소적 점프 적분을 효율적으로 근사할 수 있는가?
• RQ3 행렬 지수 방법이 고정밀 가격 책정에 있어 수렴 속도와 안정성 측면에서 피카르 반복과 같은 반복적 스킴을 능가할 수 있는가?
• RQ4 메르톤, 쿠, CGMY 등의 다양한 점프 모델을 동일한 수치적 프레임워크 내에서 일관된 정확도와 안정성로 통합 다룰 수 있는가?
• RQ5 특성 함수는 알려져 있지만 FFT 기반 방법이 격자 확장과 정확도 간의 상충 관계를 겪는 복잡한 모델, 예를 들어 헤스턴에 베이글 점프가 포함된 경우에서 이 방법은 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 방법은 테스트된 모든 점프 모델, 특히 이전 방법이 어려움을 겪던 α > 1인 CGMY 모델에서도 공간에 대해 2차 수렴성을 보였다.
• M-행렬 및 메츠러 행렬 이론을 활용한 증명을 통해, 이 스킴은 조건부 안정성과 해의 양수성 보존을 확보하였다.
• 헤스턴+VG 모델에서 제안된 방법과 FFT 기반 가격 책정 간의 상대 오차는 0.5% 미만이었으며, 깊은 아웃 오브 더 머니 만기에서는 1.7%로 증가하였다.
• 행렬 지수 계산은 특히 전이 행렬의 스펙트럼 반경이 1에 가까울 경우 피카르 반복보다 수렴 속도가 빠르며, 느린 수렴과 반복 오차 민감도 문제를 피할 수 있다.
• FFT 관련 워핑-아웃 오류를 피하고, 표준 순환 임베딩 기법에서 요구하는 격자 크기 두 배 증가가 필요하지 않다.
• 이 알고리즘은 바리어 옵션, 디지털 옵션과 같은 복잡한 옵션으로 일반화 가능하며, 향후 개선을 통해 미국형/버마디안 옵션 가격 책정에도 적합하다.