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过程强收敛且保正性的方法：截断Euler方法和后向Euler方法。显式截断Euler方法计算有效且在高波动率下保持稳健，而隐式后向Euler方法提供高计算精度和稳定性能。通过欧式看涨期权的数值实验展示了不同方法的优越性。
- 模型需求：金融工程中需要在完美校准和正确动态之间实现平衡，随机-局部波动率(SLV)模型为此提供了解决方案
- 现有模型局限性：
- 纯随机波动率模型参数有限，难以完美拟合所有期限和执行价格的隐含波动率
- 局部波动率模型虽能校准当前市场所有期权，但预测未来波动率行为存在问题，导致前向波动率微笑趋于平坦
CIR过程的数值模拟是模型评估的关键难点，因为：
- 扩散项无界且系数非全局Lipschitz连续
- 传统Euler-Maruyama方法可能产生负方差，在数学和金融上都无意义
- 现有非中心卡方近似只能实现弱收敛
开发强收敛且保正性的方法，确保精确解和近似解的轨迹彼此接近，从而改进整体期权定价。
- 提出两种新的数值方法：截断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. **方法互补性**：两种方法在效率和精度上各有优势，适用于不同场景
## 实验设置
### 模型参数
- **基准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个bins，10,000条路径，步长0.001
2. **参数敏感性分析**：测试不同$\theta$和$p$值下的表现
3. **执行价格比较**：$K = 70\%, 100\%, 150\%$三种情况
## 实验结果
### 主要结果
#### 不同步数下的误差比较
| 步数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和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. **实用价值高**：两种方法各有优势，为实际应用提供灵活选择
3. **实验全面**：系统比较了不同方法在各种条件下的表现
4. **写作清晰**：数学表述准确，结构逻辑清晰
### 不足
1. **理论分析深度**：对后向Euler方法的理论分析相对较少
2. **实际市场验证**：缺乏真实市场数据的验证
3. **计算复杂度分析**：未提供详细的复杂度理论分析
### 影响力
1. **学术贡献**：为CIR过程数值模拟提供新的强收敛方法
2. **实用价值**：直接适用于金融机构的期权定价和风险管理
3. **可复现性**：方法描述详细，易于实现
### 适用场景
1. **量化交易**：高频期权定价和对冲策略
2. **风险管理**：压力测试和尾部风险评估
3. **模型校准**：SLV模型的参数估计和验证
## 参考文献
论文引用了25篇相关文献，涵盖了SLV模型理论、CIR过程数值方法、期权定价等关键领域的重要工作，为研究提供了坚实的理论基础。