[논문 리뷰] On numerical solutions to stochastic Volterra equations

Anna Karczewska, P. Rozmej|ArXiv.org|2004. 09. 02.

Stochastic processes and financial applications참고 문헌 11인용 수 5

한 줄 요약

이 논문은 잔여 이론과 확률적 컨볼루션 근사법을 사용하여 확률적 및 결정론적 볼테라 적분 방정식을 해결하기 위한 갈레르킨 기반 수치 방법을 제시한다. α=1 및 α=2의 경우 작은 오차로 수렴성을 입증하였으며, 위너 과정 이산화 및 에르미트 다항식 사영을 통한 확률적 적분의 명시적 공식을 제공한다.

ABSTRACT

The aim of the paper is to demonstrate the use of the Galerkin method for some kind of Volterra equations, determininistic and stochastic as well. The paper consists of two parts: the theoretical and numerical one. In the first part we recall some apparently well-known results concerning the Volterra equations under consideration. In the second one we describe a numerical algorithm used and next present some examples of numerical solutions in order to illustrate the pertinent features of the technique used in the paper.

연구 동기 및 목표

• 비마르코프성과 기억 의존성 핵을 갖는 확률적 볼테라 적분 방정식을 해결하기 위한 수치 알고리즘 개발.
• 힐버트 공간 사영을 이용해 잔여 함수군과 함께 갈레르킨 방법을 확률적 볼테라 적분 방정식에 확장.
• 특정 핵 함수(α=1, α=2)에 대해 수치적 방법의 수렴성과 정확도 분석.
• 위너 과정 이산화를 통해 확률적 컨볼루션 적분의 명시적 계산 공식 제공.
• 정규 초기 자료와 가감성 노이즈를 포함한 수치 예제를 통한 방법의 시각화.

제안 방법

• 핵 함수 a(t) 및 연산자 A를 통해 정의된 잔여 함수군 S(t)를 사용하여 확률적 볼테라 적분 방정식의 온전한 해를 수식화.
• L²(ℝ)의 유한차원 부분공간에서 기저 함수 {e_k}를 사용하여 갈레르킨 방법 적용. 방정식을 {e_1, ..., e_n}의 스칼라 곱 공간에 투영.
• 시간 영역을 I개의 간격 τ = T/I로 이산화하고, 독립적인 표준 정규 변수 ζ_i를 사용하여 확률적 적분을 합으로 근사.
• α=1 및 α=2의 경우, 잔여 핵 함수 φ_α(t,x)에 대해 닫힌 형태의 표현식 유도. 이는 확률적 컨볼루션의 해석적 평가 가능하게 한다.
• 잔여 함수의 컨볼루션 성질을 활용하여 확률적 적분을 가우시안의 가중합(α=1) 또는 델타 함수 조합(α=2)으로 계산.
• 갈레르킨 투영으로 유도된 큰 선형 방정식계(1208×1208)를 해결하기 위해 LU 분해를 활용.

실험 결과

연구 질문

• RQ1 기억 핵을 갖는 확률적 볼테라 적분 방정식에 대해 갈레르킨 방법이 효과적으로 적용될 수 있는가?
• RQ2 잔여 함수가 해석적으로 알려진 경우(예: α=1,2) 수치적 해의 정확도는 어떠한가?
• RQ3 분수 또는 다항식 핵을 갖는 경우 확률적 컨볼루션 적분의 거동는 어떠한가?
• RQ4 다양한 핵 유형에 대해 시간 및 공간 이산화에 따른 수치 오차의 척도는 어떻게 되는가?
• RQ5 확률적 적분에 대한 명시적 공식을 도출하고 수치적 방법에 효율적으로 구현할 수 있는가?

주요 결과

• α=1 및 α=2의 경우, 시간 간격 [0,6] 전반에 걸쳐 결정론적 볼테라 적분 방정식의 수치적 해가 비교적 작은 오차를 보인다.
• 확률적 컨볼루션 적분은 표준 정규 변수와 해석적으로 유도된 핵 함수를 포함한 합으로 명시적으로 계산되며, α=1의 경우 가우시안, α=2의 경우 가우시안의 합이다.
• 공간 간격 m=150(=0.2) 및 기저 함수 수 n=8을 사용하여 안정적인 수치적 해를 확보하였으며, 이로 인해 1208×1208 선형 방정식계가 유도된다.
• 확률적 궤적은 결정론적 해와 확률적 적분의 합으로 구성되며, 노이즈 강도 C=0.1이다.
• 잔여 기반 통합과 갈레르킨 투영을 통해 비마르코프 성질을 효과적으로 처리하는 수치적 구현이 성공하였다.
• 확률적 성분을 위한 독립적인 표준 정규 난수 생성을 보장하기 위해 서브루틴 gasdev 및 ran1을 활용.