15 单因子利率模型蒙卡模拟

15 单因子利率模型蒙卡模拟

• 15.1 简介
  • 15.1.1 单因子模型蒙卡抽样
  • 15.1.2 由抽样结果计算债券期权价格
• 15.2 蒙卡模拟步骤
  • 15.2.1 蒙卡抽样利率变化路径
  • 15.2.2 计算零息债券期权价格
• 15.3 步骤Python代码实现
• 15.4 计算示例
• 15.5 参考资料

15.1 简介

15.1.1 单因子模型蒙卡抽样

        在单因子模型下，瞬时无风险利率的变化过程为
$$dr(t)=(\theta (t)-ar(t))dt+\sigma dz.$$
其中$a$和$\sigma$为常数，$\theta (t)$由初始零息利率期限结构确定。

        如果我们考虑在$0$至$T$时间段内，利率$r(t)$的变化过程被离散化为$steps$次变化，在不同时间区间$r(t)$的取值为$r(t_i),t_i=i\times \Delta t,i=0,1,...,steps,\Delta t=T/(steps+1)$。$r(t_i)$由瞬时无风险利率变为$t_i$至$t_i=t_i+\Delta t$之间的短期无风险利率。离散化后的$r(t_i)$的变化过程可以表示为
$$r(t_i+\Delta t)=r(t_i)+({\theta }_i-ar(t_i))\Delta t+\sigma \sqrt{\Delta t}\epsilon .$$
$\epsilon$为一标准正态分布随机数，${\theta }_i$为由初始零息利率表示的变量，下面说明如何确定${\theta }_i$。

        由《期权、期货及其他衍生产品》书中第31章内容，我们知道在利率$r(t)$连续变化的情况下，
$$\theta (t)=F_t(0,t)+aF(0,t)+\frac{{\sigma }^2}{2a}\left(1-e^{-2at}\right).$$
其中$F(0,t)$和$F_t(0,t)$分别为$t$时刻的瞬时远期利率和$t$时刻瞬时远期利率的导数。

        然后为了清晰方便起见，我们先统一将各种利率都用$R(t_0,t_1,t_2)$表示，其中$t_0$为观察时刻，$t_1$为利率作用开始时刻，$t_2$为结束时刻。此时初始零息利率$R(t)=R(0,0,t)$，离散化后的短期无风险利率$r(t_i)=R(t_i,t_i,t_i+\Delta t)$，短期远期利率$F(0,t_i)=R(0,t_i,t_i+\Delta t)$。考虑在连续复利时，
$$\exp [-R(0,0,t_i)\times t_i]\cdot \exp [-R(0,t_i,t_{i+1})\times \Delta t]=\exp [-R(0,0,t_{i+1})\times t_{i+1}],$$
得
F ( 0 , t i ) = R ( 0 , t i , t i + Δ t ) = 1 Δ t ( R ( 0 , 0 , t i + Δ t ) × ( t i + Δ t ) − R ( 0 , 0 , t i ) ×