15 单因子利率模型蒙卡模拟
15.1 简介
15.1.1 单因子模型蒙卡抽样
在单因子模型下,瞬时无风险利率的变化过程为
d r ( t ) = ( θ ( t ) − a r ( t ) ) d t + σ d z . dr(t) = (\theta(t)-ar(t))dt+\sigma dz \; . dr(t)=(θ(t)−ar(t))dt+σdz.
其中 a a a和 σ \sigma σ为常数, θ ( t ) \theta(t) θ(t)由初始零息利率期限结构确定。
如果我们考虑在 0 0 0至 T T T时间段内,利率 r ( t ) r(t) r(t)的变化过程被离散化为 s t e p s steps steps次变化,在不同时间区间 r ( t ) r(t) r(t)的取值为 r ( t i ) , t i = i × Δ t , i = 0 , 1 , . . . , s t e p s , Δ t = T / ( s t e p s + 1 ) r(t_i),\,t_i=i\times\Delta t,\, i=0,1,...,steps,\, \Delta t=T/(steps+1) r(ti),ti=i×Δt,i=0,1,...,steps,Δt=T/(steps+1)。 r ( t i ) r(t_i) r(ti)由瞬时无风险利率变为 t i t_i ti至 t i = t i + Δ t t_i=t_i+\Delta t ti=ti+Δt之间的短期无风险利率。离散化后的 r ( t i ) r(t_i) r(ti)的变化过程可以表示为
r ( t i + Δ t ) = r ( t i ) + ( θ i − a r ( t i ) ) Δ t + σ Δ t ε . r(t_i+\Delta t) = r(t_i)+(\theta_i-ar(t_i))\Delta t+\sigma\sqrt{\Delta t}\varepsilon\; . r(ti+Δt)=r(ti)+(θi−ar(ti))Δt+σΔtε.
ε \varepsilon ε为一标准正态分布随机数, θ i \theta_i θi为由初始零息利率表示的变量,下面说明如何确定 θ i \theta_i θi。
由《期权、期货及其他衍生产品》书中第31章内容,我们知道在利率 r ( t ) r(t) r(t)连续变化的情况下,
θ ( t ) = F t ( 0 , t ) + a F ( 0 , t ) + σ 2 2 a ( 1 − e − 2 a t ) . \theta(t) = F_t(0, t)+aF(0,t)+\frac{\sigma^2}{2a}\left(1-e^{-2at}\;\;\right)\;. θ(t)=Ft(0,t)+aF(0,t)+2aσ2(1−e−2at).
其中 F ( 0 , t ) F(0,t) F(0,t)和 F t ( 0 , t ) F_t(0,t) Ft(0,t)分别为 t t t时刻的瞬时远期利率和 t t t时刻瞬时远期利率的导数。
然后为了清晰方便起见,我们先统一将各种利率都用 R ( t 0 , t 1 , t 2 ) R(t_0, t_1, t_2) R(t0,t1,t2)表示,其中 t 0 t_0 t0为观察时刻, t 1 t_1 t1为利率作用开始时刻, t 2 t_2 t2为结束时刻。此时初始零息利率 R ( t ) = R ( 0 , 0 , t ) R(t)=R(0,0,t) R(t)=R(0,0,t),离散化后的短期无风险利率 r ( t i ) = R ( t i , t i , t i + Δ t ) r(t_i)=R(t_i,t_i,t_i+\Delta t) r(ti)=R(ti,ti,ti+Δt),短期远期利率 F ( 0 , t i ) = R ( 0 , t i , t i + Δ t ) F(0,t_i)=R(0,t_i,t_i+\Delta t) F(0,ti)=R(0,ti,ti+Δt)。考虑在连续复利时,
exp [ − R ( 0 , 0 , t i ) × t i ] ⋅ exp [ − R ( 0 , t i , t i + 1 ) × Δ t ] = exp [ − R ( 0 , 0 , t i + 1 ) × t i + 1 ] , \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}\,]}\;, exp[−R(0,0,ti)×ti]⋅exp[−R(0,ti,ti+1)×Δt]=exp[−R(0,0,ti+1)×ti+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 ) ×