Heston Model 总结

Heston Model 总结

---

Heston 模型基本方程为，
$$\begin{array}{rl}d{S}_t& =r{S}_{t}dt+\sqrt{V_t}{S}_t[\rho d{W}_t^{(1)}+\sqrt{1-{\rho }^2}d{W}_t^{(2)}],\\ d{V}_t& =\kappa (\theta -V_t)dt+{\sigma }_v\sqrt{V_t}d{W}_t^{(1)}.\end{array}$$
Broadie 和 Kaya (2006) 给出了 exact simulation，但由于其计算量过大，很难实际应用，后续有众多基于精确模拟算法的近似方法，在尽可能保留精确度的基础上（超过 Euler 格式和 Milstein 格式），尽量降低计算量。

---

精确模拟算法

Heston 模型的精确模拟分为三个主要步骤：（1）方差过程的模拟；（2）方差积分过程的模拟；（3）价格过程的模拟。其中方差过程的模拟涉及非中心卡方分布，一般的模拟方法是拒绝接受采样；方差积分过程涉及特征函数的计算和傅里叶变换；价格过程的模拟较为简单。

方差过程的模拟

$V_t|V_u$ 服从非中心卡方分布，
$$V_t=\frac{{\sigma }_v^2(1-e^{-\kappa (t-u)})}{4\kappa }{\chi }_d^{\prime 2}\left(\frac{4\kappa e^{-\kappa (t-u)}}{{\sigma }_v^2(1-e^{-\kappa (t-u)})}{V}_u\right),d=\frac{4\theta \kappa }{{\sigma }_v^2}$$

方差积分过程的模拟

计算方差积分过程的特征函数，
$$\begin{array}{rl}\Psi (a,v(s),v(t))=& E\left[\exp \left(ia{\int }_s^{t}v(u)du\right)|v(s),v(t)\right]\\ =& \frac{\gamma (a)e^{\frac{1}{2}(\gamma (a)-\kappa )(t-s)}(1-\exp (-\kappa (t-s)))}{\kappa (1-e^{-\gamma (a)(t-s)})}\\ & \times \exp \left[\frac{v(s)+v(t)}{{\xi }^2}\left(\frac{\kappa (1+e^{-\kappa (t-s)})}{1-e^{-\kappa (t-s)}}-\frac{\gamma (a)(1+e^{-\gamma (a)(t-s)})}{1-e^{-\gamma (a)(t-s)}}\right)\right]\\ & \times I_{\frac{1}{2}d-1}[\sqrt{v(s)v(t)}4\gamma (a)e^{-\frac{\kappa }{2}(t-s)}/{\xi }^2(1\end{array}$$