《金利変動モデルと債券価格》

【きんりへんどうもでるとさいけんかかく (interest rate model and bond price)】

　金利変動モデル(interest rate model)を仮定することによって, 債券の理論価格を求めることができる. 以下, $(\Omega ,{\mathcal {F}},P,{\mathcal {F}}_{t})\,$をフィルトレーション付き確率空間, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle T_0 < \infty\, } をtime horizon とし, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(t,T)\, } を時刻構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle T\, } を満期とする割引債の時刻構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\, } での価格とする($0\leq t\leq T\leq T_{0}\,$).

　確率1で, 各t毎に構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle P(t,T)\, } は$T\in (t,T_{0}]\,$について微分可能かつ連続でその微分値は$T\in (t,T_{0}]\,$について有界とする. このとき時刻$t\,$における, 時刻$T\,$からの(瞬間的な)フォワードレート$f(t,T)\,$は, $\displaystyle -{\frac {\partial \log P(t,T)}{\partial T}}\,$で定義される. 逆に, フォワードレート$f(t,T)\,$, $0\leq t\leq T\leq T_{0}\,$が与えられたとき, 債券価格$P(t,T)\,$は, 定義から$\textstyle P(t,T)=\exp(-\int _{t}^{T}f(t,s){\mbox{d}}s)\,$となる.

　Heath, Jarrow, Morton[2]は, フォワードレート$f(t,T)\,$が次の確率微分方程式

${\mbox{d}}f(t,T)=\alpha (t,T){\mbox{d}}t+\sum _{n=1}^{N}\sigma _{n}(t,T){\mbox{d}}B_{t}^{n}\,$

(但し$\alpha (t,T)\,$, $\sigma _{n}(t,T)\,$は${\mathcal {F}}_{t}\,$適合な確率過程, $B_{t}^{n}\,$は標準ブラウン運動($n=1,2,\ldots ,N)\,$) に従うと仮定して金利変動モデルのかなり一般的な枠組みを与えた.

　$\displaystyle r_{t}=\lim _{T\downarrow t}f(t,T)\,$をスポットレートという. スポットレートの変動を確率過程としてモデル化したものはスポットレートモデルと呼ばれ, その下では, 債券価格$P(t,T)\,$は$\textstyle P(t,T)=E^{Q}[\exp(-\int _{t}^{T}r_{s}{\mbox{d}}s)|{\mathcal {F}}_{t}]\,$ となる. ただし, $E[\cdot ]\,$は同値マルチンゲール測度$Q\,$の下での期待値を表す.

　スポットレートモデルの代表例として, Vasicek[4]は, $r_{t}\,$が確率過程

${\mbox{d}}r_{t}=(\varphi -ar_{t}){\mbox{d}}t+\sigma {\mbox{d}}B_{t}\,$

($\varphi \,$, $a\,$, $\sigma \,$は定数)に従うとするモデル(Vasicek モデル)を立てた. またCox, Ingersoll, Ross [1]は, $r_{t}\,$が確率過程

${\mbox{d}}r_{t}=(\varphi -ar_{t}){\mbox{d}}t+\sigma {\sqrt {r_{t}}}{\mbox{d}}B_{t}\,$

に従うとするモデル(C.I.R.モデル)を立てた.

　これらのモデルの下では, 時刻Tを満期とする割引債価格の満たす確率微分方程式は,

${\frac {{\mbox{d}}P(t,T)}{P(t,T)}}=\mu ^{T}(r_{t},t){\mbox{d}}t+\sigma ^{T}(r_{t},t){\mbox{d}}B_{t}\,$

と書けるが, このとき$\displaystyle \lambda _{t}={\frac {\mu ^{T}(r_{t},t)-r_{t}}{\sigma ^{T}(r_{t},t)}}\,$は$T\,$には依らないことが示される. この$\lambda _{t}\,$はリスクの市場価値と呼ばれる. リスクの市場価値$\lambda _{t}\,$を既知とするとき, 同値マルチンゲール測度$Q\,$は

${\frac {{\mbox{d}}Q}{{\mbox{d}}P}}=\exp \left(-\int _{0}^{T_{0}}\lambda _{t}{\mbox{d}}B_{t}-{\frac {1}{2}}\int _{0}^{T_{0}}\lambda _{t}^{2}{\mbox{d}}t\right)\,$

で与えられる. Vasicekモデルの割引債価格は

$P(t,T)\,$	$=\,$	$A(T-t){\mbox{e}}^{-B(T-t)r(t)}\,$
$B(t)\,$	$=\,$	${\frac {1-{\mbox{e}}^{-at}}{a}}\,$
$A(t)\,$	=	$\exp\{-{\frac {(B(t)-t)(a^{2}{\bar {\varphi }}-\sigma ^{2}/2)}{a^{2}}}-{\frac {\sigma ^{2}B^{2}(t)}{4a}}\}\,$
${\bar {\varphi }}\,$	=	${\frac {\varphi -\sigma \lambda }{a}}\,$

となり, C.I.R.モデルの割引債価格は

$P(t,T)\,$	$=\,$	$A(T-t){\mbox{e}}^{-B(T-t)r(t)}\,$
$B(t)\,$	$=\,$	${\frac {2({\mbox{e}}^{-\gamma t}-1)}{(a+\gamma )({\mbox{e}}^{-\gamma t}-1)+2\gamma }}\,$
$A(t)\,$	$=\,$	$\left[{\frac {2\gamma {\mbox{e}}^{(a+\gamma )t/2}}{(a+\gamma )({\mbox{e}}^{-\gamma t}-1)+2\gamma }}\right]^{2a{\bar {\varphi }}/\sigma ^{2}}\,$
$\gamma \,$	$=\,$	${\sqrt {a^{2}+2\sigma ^{2}}}\,$
${\bar {\varphi }}\,$	$=\,$	構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\frac {\varphi -\sigma {\sqrt {r_{t}}}\lambda }{a}}\,}

となる.

　HullとWhite[3]はVasicek モデルにおけるパラメータ$\varphi \,$を時刻$t\,$の確定的な関数として拡張することにより, 市場で観測される [3] イールドカーブ(yield curve)に整合するモデルを提唱した. 具体的には次の通り. $y(t),0\leq t\leq T_{0}\,$を割引債利回りのカーブとし, $t\,$について2階微分可能と仮定する. このモデルでは, $r_{t}\,$は同値マルチンゲール測度$Q\,$の下で

${\mbox{d}}r_{t}\,$	$=\,$	構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle (\varphi (t)-a(t)r_{t}){\mbox{d}}t+\sigma (t){\mbox{d}}B_{t}^{Q}\,}
		$\;\;B^{Q}\,$は$Q\,$ブラウン運動, 構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \varphi (t),a(t),\sigma (t)\,} は確定的な関数
$\varphi (t)\,$	$=\,$	構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_0^t \sigma^2(u) \mbox{e}^{-2\int_u^t a(v) \mbox{d}v}\mbox{d}u +a(t)(ty(t))' +(ty(t))'' \, }

を満たし, 割引債価格は

$P(t,T)\,$	$=\,$	$H_{1}(t,T){\mbox{e}}^{-H_{2}(t,T)r(t)}\,$
構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle H_{2}(t,T)\,}	$=\,$	$\int _{t}^{T}{\mbox{e}}^{-\int _{t}^{s}a(v){\mbox{d}}v}{\mbox{d}}s\,$
構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle H_{1}(t,T)\,}	$=\,$	$\exp\{{\frac {1}{2}}\int _{t}^{T}\sigma ^{2}(u)H_{2}^{2}(u,T){\mbox{d}}u-\int _{t}^{T}\varphi (u)H_{2}(u,T){\mbox{d}}u\}\,$

で得られる. これを変形すると,

$P(t,T)\,$	$=\,$	構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\frac {P(0,T)}{P(0,t)}}\exp \left\{(\theta (t)-r(t))H_{2}(t,T){\frac {}{}}\right.\,}
		$\left.+{\frac {1}{2}}\int _{t}^{T}\sigma ^{2}(u)H_{2}^{2}(u,T){\mbox{d}}u-{\frac {1}{2}}\int _{0}^{T}\sigma ^{2}(u)H_{2}^{2}(u,T){\mbox{d}}u+{\frac {1}{2}}\int _{0}^{t}\sigma ^{2}(u)H_{2}^{2}(u,t){\mbox{d}}u\right\}\,$
構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \theta (t)\, }	$=\,$	構文解析に失敗 (Conversion error. Server ("https://en.wikipedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \int _{0}^{t}\sigma ^{2}(u)H_{2}(u,t){\mbox{e}}^{-\int _{u}^{t}a(v){\mbox{d}}v}{\mbox{d}}u+(tY(t))\,}

となる.

---

参考文献

[1] J. C. Cox, J. E. Ingersoll and S. A. Ross, "A Theory of the Term Structure of Interest Rates," Econometrica, 53 (1985), 385-407.

[2] D. Heath, R. Jarrow and A. Morton, "Bond Pricing and the Term Structure of Interest Rates: A New Methodology for Contingent Claims Valuation," Econometrica, 60 (1992), 77-105.

[3] J. Hull and A. White, "Pricing Interest-rate-derivative Securities," The Review of Financial Studies, 3 (1990), 573-592.

[4] Vasicek, O.A., "An equilibrium characterization of the term structure," Journal of Financial Economics, 5 (1977), 177-188.