The Heston stochastic-local volatility model, consisting of a asset price process and a Cox--Ingersoll--Ross-type variance process, offers a wide range of applications in the financial industry. The pursuit for efficient model evaluation has been assiduously ongoing and central to which is the numerical simulation of CIR process. Different from the weakly convergent noncentral chi-squared approximation used in 25, this paper considers two strongly convergent and positivity-preserving methods for CIR process under Lamperti transformation, namely, the truncated Euler method and the backward Euler method. It should be noted that these two methods are completely different. The explicit truncated Euler method is computationally effective and remains robust under high volatility, while the implicit backward Euler method provides high computational accuracy and stable performance. Numerical experiments on European call options are presented to show the superiority of different methods.
academic- 論文ID: 2509.24449
- タイトル: Efficient simulation of prices for European call options under Heston stochastic-local volatility model: a comparison of methods
- 著者: Meng Cai, Tianze Li
- 著者所属: 中央財経大学統計数学学院、経済学院
- 分類: q-fin.CP（計算ファイナンス）
- 発表日: 2025年10月16日
- 論文リンク: https://arxiv.org/abs/2509.24449
本論文は、資産価格過程とCox-Ingersoll-Ross型分散過程から構成されるHeston確率-局所ボラティリティモデルを研究している。このモデルは金融業界で広く応用されている。先行研究25で使用された弱収束の非中心カイ二乗近似とは異なり、本論文ではLamperti変換下でCIR過程に対して強収束かつ正性を保つ2つの方法を検討している：切断Euler法と後退Euler法である。明示的切断Euler法は計算効率が高く、高ボラティリティ下でも堅牢性を保つ一方、暗黙的後退Euler法は高い計算精度と安定した性能を提供する。ヨーロッパコールオプションの数値実験を通じて、異なる方法の優位性を示している。
- モデルの必要性：金融工学では完全な校正と正確なダイナミクスのバランスが必要であり、確率-局所ボラティリティ(SLV)モデルがこの解決策を提供する
- 既存モデルの限界：
- 純粋な確率ボラティリティモデルはパラメータが限定的であり、すべての満期と行使価格のインプライドボラティリティを完全に適合させることが困難
- 局所ボラティリティモデルは現在の市場のすべてのオプションを校正できるが、将来のボラティリティ行動の予測に問題があり、フォワードボラティリティスマイルが平坦化する傾向
CIR過程の数値シミュレーションはモデル評価の重要な困難点である。理由は以下の通り：
- 拡散項が無界であり、係数が大域的Lipschitz連続でない
- 従来のEuler-Maruyama法は負の分散を生成する可能性があり、数学的および金融的に無意味
- 既存の非中心カイ二乗近似は弱収束のみを実現
強収束かつ正性を保つ方法を開発し、精密解と近似解の軌跡が互いに接近することを保証し、全体的なオプション価格設定を改善する。
- 2つの新しい数値方法の提案：切断Euler法と後退Euler法、両者とも強収束性と正性を有する
- 理論分析：切断Euler法の強収束性（収束次数1/2）を証明
- 包括的な数値比較：ヨーロッパコールオプション価格設定実験を通じて、異なる方法の性能を体系的に比較
- 実用的ガイダンス：異なる応用シナリオに対する方法選択の推奨を提供
Heston確率-局所ボラティリティモデル下のヨーロッパコールオプション価格設定問題を解く。その中核はCIR分散過程を効率的かつ正確にシミュレートすることである。
dS_t = rS_t dt + \sqrt{V_t} S_t \sigma(t, S_t)\left(\rho dW_t + \sqrt{1-\rho^2} d\tilde{W}_t\right) \\
dV_t = \kappa(\theta - V_t)dt + \gamma\sqrt{V_t} dW_t \\
dW_t d\tilde{W}_t = \rho dt
\end{cases}$$
ここで：
- $S_t$：資産価格過程
- $V_t$：分散過程（CIRダイナミクス）
- $\sigma(t,s)$：局所ボラティリティ関数
- $\kappa, \theta, \gamma$：CIR過程パラメータ
#### Lamperti変換
CIR過程の数値的困難を克服するため、変換$L_t = \sqrt{V_t}$を導入：
$$dL_t = \frac{1}{2}\kappa\left(\frac{\theta}{L_t} - L_t\right)dt + \frac{1}{2}\gamma dW_t$$
### 中核的数値方法
#### 1. 切断Euler法
切断写像を定義：$\pi_\tau(x) = (b\tau^{1/4}) \vee x$
数値スキーム：
$$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{\pi_\tau(L^N_{t_n})} - \pi_\tau(L^N_{t_n})\right) + \frac{1}{2}\gamma\Delta W_n$$
**特徴**：
- 明示的スキーム、計算効率が高い
- 切断を通じて正性を保証
- 高ボラティリティ下で堅牢性を維持
#### 2. 後退Euler法
$$L^N_{t_{n+1}} = L^N_{t_n} + \frac{\kappa}{2}\tau\left(\frac{\theta}{L^N_{t_{n+1}}} - L^N_{t_{n+1}}\right) + \frac{1}{2}\gamma\Delta W_n$$
**特徴**：
- 暗黙的スキーム、反復解法が必要
- 高い計算精度を提供
- 性能は安定だが計算コストが高い
### 技術的革新点
1. **強収束性**：弱収束の非中心カイ二乗法と異なり、軌跡レベルの収束を保証
2. **正性設計**：数学的変換と切断技術を通じて、分散が常に正であることを保証
3. **方法の相補性**：2つの方法は効率と精度において各々の優位性を有し、異なるシナリオに適用可能
## 実験設定
### モデルパラメータ
- **ベンチマークHestonパラメータ**：$\gamma = 0.95$, $\kappa = 1.05$, $\rho = -0.315$, $\theta = 0.0855$, $S_0 = 1$, $r = 0$, $v_0 = 0.0945$, $T = 5$
- **HSLV修正パラメータ**：修正パラメータ$p = 0.25$を導入し、各パラメータを調整
### 評価指標
- 相対誤差百分率
- 計算時間
- 異なる市場条件下での堅牢性
### 比較方法
- 標準Euler法
- AES法（Andersen Exact Simulation）
- 切断Euler法
- 後退Euler法
### 実験設計
1. **条件付き期待値計算**：2D-COS法をベンチマークとして使用、20ビン、10,000パス、ステップサイズ0.001
2. **パラメータ感度分析**：異なる$\theta$および$p$値下での性能をテスト
3. **行使価格比較**：$K = 70\%, 100\%, 150\%$の3つのケース
## 実験結果
### 主要結果
#### 異なるステップ数下での誤差比較
| ステップ数N | K=70% | K=100% | K=150% |
|-------|--------|---------|---------|
| **切断Euler法** |
| 5 | 13.49%(0.24%) | 14.08%(0.24%) | 13.96%(0.85%) |
| 40 | 4.66%(0.53%) | 3.86%(0.74%) | 4.80%(1.35%) |
| **後退Euler法** |
| 5 | 15.49%(0.28%) | 15.97%(0.10%) | 17.32%(1.58%) |
| 40 | 0.06%(0.66%) | 0.46%(0.68%) | 1.44%(1.15%) |
#### 計算時間の比較（40ステップ）
- AES法：0.1114秒
- Euler法：0.0939秒
- 切断Euler法：0.0733秒
- 後退Euler法：137.9476秒
### 主要な発見
1. **ステップ数が少ない場合**：切断Euler法は通常最も低い誤差を示す
2. **ステップ数が多い場合**：後退Euler法は最高の精度を示す
3. **パラメータ堅牢性**：切断Euler法は$\theta$の増加に伴い誤差増加が緩やかであり、優れた堅牢性を示す
4. **計算効率**：切断Euler法が最も高速であり、後退Euler法は反復解法が必要なため計算時間が大幅に増加
### アブレーション実験
#### ボラティリティパラメータ$\theta$の影響
実験は$\theta$の増加（極端な市場条件をシミュレート）に伴い以下を示している：
- AES法と標準Euler法の誤差が著しく増加
- 切断Euler法は堅牢性を保ち、誤差増加率が低い
- これはストレステストと裾部リスク建模において切断法がより適切であることを示唆
#### 修正パラメータ$p$の影響
$p$の増加に伴い、後退Euler法は常に最良の性能を示すが、ステップサイズに対する感度が高い。
## 関連研究
### 主要な研究方向
1. **SLVモデル校正**：Guyon and Henry-Labordère (2011)の先駆的研究
2. **CIR過程の数値方法**：有限差分、Monte Carlo粒子法を含む
3. **深層学習の応用**：Cuchiero等(2020)が生成対抗ネットワークを用いた校正を使用
### 本論文の優位性
既存研究と比較して、本論文は強収束の正性保持法を提供し、弱収束法の不足を補い、実際の応用に対してより信頼性の高い数値ツールを提供する。
## 結論と考察
### 主要な結論
1. **方法の相補性**：切断Euler法は大規模シミュレーションと実時間リスク管理に適し、後退Euler法はストレステストと高精度要件シナリオに適する
2. **堅牢性の優位性**：切断Euler法は極端な市場条件下で優れた性能を示す
3. **精度-効率のトレードオフ**：具体的な応用要件に応じて適切なステップサイズと方法を選択する必要がある
### 限界
1. 後退Euler法の計算コストが高い
2. ステップサイズの選択は精度と効率のバランスが必要
3. 極端なパラメータ組み合わせに対する方法の安定性はさらなる検証が必要
### 今後の方向
1. 適応的ステップサイズ戦略の開発
2. 並列化実装の探索
3. より複雑なSLVモデル変種への拡張
## 深層的評価
### 利点
1. **理論的厳密性**：完全な収束性証明を提供
2. **実用的価値が高い**：2つの方法は各々の優位性を有し、実際の応用に柔軟な選択肢を提供
3. **実験が包括的**：異なる条件下での異なる方法の性能を体系的に比較
4. **記述が明確**：数学的表現が正確で、構造が論理的に明確
### 不足
1. **理論分析の深さ**：後退Euler法の理論分析が相対的に少ない
2. **実際の市場検証**：実際の市場データによる検証が欠如
3. **計算複雑度分析**：詳細な複雑度理論分析が提供されていない
### 影響力
1. **学術的貢献**：CIR過程の数値シミュレーションに新しい強収束法を提供
2. **実用的価値**：金融機関のオプション価格設定とリスク管理に直接適用可能
3. **再現性**：方法の説明が詳細で、実装が容易
### 適用シナリオ
1. **量的取引**：高頻度オプション価格設定とヘッジ戦略
2. **リスク管理**：ストレステストと裾部リスク評価
3. **モデル校正**：SLVモデルのパラメータ推定と検証
## 参考文献
論文は25篇の関連文献を引用しており、SLVモデル理論、CIR過程の数値方法、オプション価格設定などの重要分野の重要な研究を網羅し、研究に堅実な理論的基礎を提供している。