|
双曲局部波动率Hull-White模型在我们的第二个示例中,我们考虑一个由dS(t)St=r(t)dt+σH(St)dW(t),S(0)=S,dr(t)=a(θ(t)- r(t))dt+σ(ρdW(t)+p1- ρdW(t)),r(0)=r,(73),其中σH(St)=νn(1- β + β)β+(β - 1) βStqSt+β(1- St)- βo、 (74)ν>0表示波动性水平和β∈ [0,1]显示倾斜参数。文献[39]中介绍的该模型与CEV模型的行为非常相似,并已用于文献[40,41]中的数值试验。它的优点是避免零成为可达到的边界,然后允许避免一些数值不稳定性,如基础S为0时的CEV模型中所示(参见例[42])。它对应于β=1的Black-Scholes模型,并且当β6=1时,波动率表面呈现出倾斜。图15说明了参数β对波动率曲面倾斜的影响。我们观察到,随着β值的降低,偏斜显著增加。例如,当ν=0.2,β=0.2时,50%和100%冲击之间的波动性差异约为20%。我们通过分别考虑负相关ρ=-30%,正相关ρ=30%。表2中选择了其他模型参数。对于这两个测试,我们在ds=0.012、dr=0.002、dt=0.0099的统一网格中求解P Z至成熟度T的PDEin方程(31)。然后,我们利用PDE网格点和解,通过数值积分,对不同的敲打进行成熟度为T的欧式看涨期权定价。我们还通过对具有dt=和100万条路径的SDE(73)使用ingueler离散化进行蒙特卡罗模拟定价(参见例如[37]或[43])。对于集合1(分别针对集合2),定价结果如图16(分别在图18中)所示,其差异如图17(分别在图19中)所示。
|