《金利変動モデルと債券価格》のソースを表示

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

　[[金利変動モデル]](interest rate model)を仮定することによって, 債券の理論価格を求めることができる. 以下, <math>(\Omega, {\mathcal F}, P, {\mathcal F}_t )\, </math>をフィルトレーション付き確率空間, <math>T_0 < \infty\, </math>をtime horizon とし, <math>P(t,T)\, </math>を時刻<math>T\, </math>を満期とする割引債の時刻<math>t\, </math>での価格とする(<math>0 \leq t \leq T \leq T_0\, </math>). 

　確率1で, 各<math>t \, </math>毎に<math>P(t,T)\, </math>は<math>T \in (t, T_0]\, </math>について微分可能かつ連続でその微分値は<math>T\in (t, T_0]\, </math>について有界とする. このとき時刻<math>t\, </math>における, 時刻<math>T\, </math>からの(瞬間的な)フォワードレート<math>f(t,T)\, </math>は, <math>\displaystyle - \frac{\partial \log P(t,T)}{\partial T}\, </math>で定義される. 逆に, フォワードレート<math>f(t,T)\, </math>, <math>0 \leq t \leq T \leq T_0\, </math>が与えられたとき, 債券価格<math>P(t,T)\, </math>は, 定義から<math>\textstyle P(t,T)=\exp (-\int_t^T f(t,s)\mbox{d}s)\, </math>となる.

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


<center><math>\mbox{d}f(t,T)=\alpha(t,T)\mbox{d}t + \sum_{n=1}^{N} \sigma_n(t,T)\mbox{d}B^n_t\, </math></center>


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

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

　スポットレートモデルの代表例として, Vasicek[4]は, <math>r_t\, </math>が確率過程


<center><math>\mbox{d}r_t = (\varphi - ar_t)\mbox{d}t + \sigma \mbox{d}B_t\, </math></center>


(<math>\varphi\, </math>, <math>a\, </math>, <math>\sigma\, </math>は定数)に従うとするモデル(Vasicek モデル)を立てた. またCox, Ingersoll, Ross [1]は, <math>r_t\, </math>が確率過程


<center><math>\mbox{d}r_t = (\varphi - ar_t)\mbox{d}t + \sigma \sqrt{r_t} \mbox{d}B_t\, </math></center> 


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

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


<center><math>\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\, </math></center> 


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


<center>
<math>\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)\, </math>
</center> 


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


<table align="center">
<tr>
<td align="right"><math>P(t,T)\, </math></td>
	<td><math>=\, </math></td>
		<td><math>A(T-t) \mbox{e}^{-B(T-t)r(t)}\, </math> </td>
</tr>
<tr>
<td align="right"><math>B(t)\, </math></td>
	<td><math>=\, </math></td>
		<td><math>\frac{1-\mbox{e}^{-at}}{a}\, </math> </td>
</tr>
<tr>
<td align="right"><math>A(t)\, </math> </td>
	<td>=</td>
		<td><math>\exp \{-\frac{(B(t)-t)(a^2\bar{\varphi}-\sigma^2/2)}{a^2}
                                 -\frac{\sigma^2B^2(t)}{4a} \}\, </math></td>
</tr>
<tr>
<td align="right"><math>\bar{\varphi}\, </math></td>
	<td>=</td>
		<td><math>\frac{\varphi -\sigma\lambda}{a}\, </math> </td>
</tr>
</table>


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


<table align="center">
<tr>
<td align="right"><math>P(t,T) \, </math></td>
	<td><math> = \, </math></td>
			<td><math> A(T-t) \mbox{e}^{-B(T-t)r(t)} \, </math></td>
</tr>
<tr>
<td align="right"><math> B(t) \, </math></td>
	<td><math> = \, </math></td>
			<td><math>\frac{ 2(\mbox{e}^{-\gamma t}-1)}{(a+\gamma)(\mbox{e}^{-\gamma t}-1)+2\gamma} \, </math></td>
</tr>
<tr>
<td align="right"><math>A(t) \, </math></td>
	<td><math> = \, </math></td>
			<td><math> \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}\, </math></td>
</tr>
<tr>
<td align="right"><math> \gamma \, </math></td>
	<td><math> = \, </math></td>
			<td><math>\sqrt{a^2 + 2\sigma^2} \, </math></td>
</tr>
<tr>
<td align="right"><math>\bar{\varphi} \, </math></td>
	<td><math> = \, </math></td>
			<td><math>\frac{\varphi -\sigma \sqrt{r_t}\lambda}{a} \, </math></td>
</tr>
</table>


となる. 

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


<table align="center">
<tr>
<td align="right"><math>\mbox{d}r_t \, </math></td>
	<td><math> =\, </math></td>
			<td><math> (\varphi (t)-a(t)r_t)\mbox{d}t+\sigma(t) \mbox{d}B^Q_t\, </math></td>
</tr>
<tr>
<td align="right"></td>
	<td></td>
			<td><math>\; \; B^Q \, </math>は<math>Q\, </math>ブラウン運動, 
<math>\varphi (t),a(t),\sigma(t)\, </math>は確定的な関数</td></tr>
<tr>
<td align="right"><math>\varphi(t)\, </math></td>
	<td><math>= \, </math></td>
			<td><math>\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))'' \, </math></td>
</tr>
</table>


を満たし, 割引債価格は


<table align="center">
<tr>
<td align="right"><math> P(t,T)\, </math></td>
	<td><math>= \, </math></td>
			<td><math>H_1(t,T) \mbox{e}^{-H_2(t,T)r(t)} \, </math></td>
</tr>
<tr>
<td align="right"><math>H_2(t,T) \, </math></td>
	<td><math> =\, </math></td>
			<td><math> \int _t^T \mbox{e}^{-\int _t^s a(v)\mbox{d}v}\mbox{d}s \, </math></td>
</tr>
<tr>
<td align="right"><math> H_1(t,T)\, </math></td>
	<td><math>= \, </math></td>
			<td><math> \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\}\, </math></td>
</tr>
</table>


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


<table align="center">
<tr>
<td align="right"><math> P(t,T)\, </math></td>
	<td><math> =\, </math></td>
			<td><math> \frac{P(0,T)}{P(0,t)}\exp \left\{(\theta(t)-r(t))H_2(t,T)\frac{}{}\right.\, </math></td>
</tr>
<tr>
<td align="right"></td>
	<td></td>
			<td><math>\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\} \, </math></td>
</tr>
<tr>
<td align="right"><math> \theta (t)\, </math></td>
	<td><math> =\, </math></td>
			<td><math>\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)) \, </math></td>
</tr>
</table>


となる. 



----
'''参考文献'''

[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.

[[category:ファイナンス|きんりへんどうもでるとさいけんかかく]]