摘 要 假設股票價格遵循指數OU過程,利用隨機分析中的鞅方法,得到了具有隨機波動率的歐式期權的定價公式,推廣了BS模型.
關鍵詞 期權定價;隨機波動率;OU過程
中圖分類號 F830.9 文獻標識碼:A
1 引 言
期權定價理論是現在金融理論的核心內容,主要的革命性成果是Black和Scholes(1973)的BlackScholes期權定價公式[1].標的資產的價格的波動率是用來度量一定時間內,標的資產價格變動的不確定性,它是影響期權價格的重要因素,是期權定價模型中重要的變量.BS期權定價公式假定波動率是常數,然而大量實證研究表明波動率是一個隨機變量.Scott[2]、Wiggins[3]、Renault和Touzi[4]等人已經把BS模型推廣到隨機波動率的情形,但這些情形大多數沒有給出解析解,并且需要使用數值分析方法得到期權價格.1987年,Hull和White[5]提出著名的波動率平方服從幾何布朗運動,利用Taylor展式,近似地給出了期權定價公式,在理論上產生了較大的影響,但還是沒有給出精確的解析解.陳俊霞、蹇明[6]在他們二人模型的假設下, 利用鞅方法, 推導出歐式期權定價的解析解.本文作者討論了股票價格遵循指數OU過程,具有隨機波動率的歐式期權定價的問題.
2 模型及若干引理
給定具有濾子的概率空間{Ω,F,Ftt≥0,P},其中濾子流Ftt≥0滿足通常條件(完備、單調遞增、右連續),F0是平凡的且FT=F.設市場上有兩種資產,一種為無風險資產,稱為債券,其價格過程p(t)滿足微分方程:
dp(t)=p(t)r(t)dt,p(0)=1,
其中r(t)是瞬時利率(是無風險利率),則p(t)=e-∫t0r(s)ds.
另一種為風險資產(如股票),其價格過程S=(St)0≤t≤T滿足隨機微分方程
dSt=(μ(t)-αln St)Stdt+σ(t)StdBt,(1)
dσ2(t)=uσ2(t)dt+ξ σ2(t)dWt,(2)
其中S0=S>0,u,ξ為常數,μ(t)為時間t的確定性函數,σ(t)為股票的波動率,α為非負常數, B=(Bt)0≤t≤T和W=(Wt)0≤t≤T是Ft適應的、相互獨立的布朗運動.
注1 當α=0時式(1)即為BlackScholes模型中的風險資產價格方程,當α>0時意味著當股票價格St上升到一定高度后,它使St有下降的趨勢.此時模型稱為指數OU過程[7].
注2 由于波動率是不可交易資產,毫無疑問這里的市場模型是不完備的,因而等價鞅測度以及相應的價格過程是不唯一的,即不可能得到與投資者風險態度無關的唯一定價.關于如何選取一個恰當的鞅測度,可參閱文獻[8].本文考慮當市場中鞅測度選定以后,如何進行期權定價.
令Xt=Ste-∫t0r(s)ds, 則X=(Xt)0≤t≤T為風險資產貼現價格過程,過程X具有微分
dXt=(μ(t)-r(t)-αln St)Xtdt+σ(t)XtdBt.
下面旨在給出等價鞅測度,定義過程Z=(Zt)0≤t≤T,這里
Zt=exp (-∫t0μ(s)-r(s)-αln Ssσ(s)dBs-12∫t0(μ(s)-r(s)-αln Ssσ(s))2ds).
引理1 E(ZT)=1,因而過程Z是P鞅.
證明 注意到過程B與W的獨立性,由Liptser Shiryayev[9]的第6章第2節例題4可證E(ZT)=1,引理的第二個結論是顯然的.
定義測度Q,使得dQ/dP=ZT,于是由引理1得出,Q是P的等價鞅測度.記BQt=Bt+∫t0μ(s)-r(s)-αln Ssσ(s)ds,由Girsanov定理可知:BQt是關于Q的布朗運動.若當前時刻為0時刻,由It公式知:St=Sexp {∫t0(r(s)-12σ2(s))ds+∫t0σ(t)dBQs}.(3)
引理2[6] ∫t0σ(s)dBQs服從N(0,∫t0σ2(s)ds).
引理3[10]記p(t,a,b;x)為隨機變量∫t0exp (au+bBu)du的密度函數,對于x>0,有
p(t,a,b;x)=M(x)∫
SymboleB@ 0N(v)[∫
SymboleB@ 0y2ab2exp (-2b2x(y2+2ycosh(v)+1))dy]dv,
經 濟 數 學第 27 卷
第2期朱利芝等:隨機波動率下股價服從OU過程的期權定價
其中
M(x)=8(πb3x22πt)-1exp (4π2-a2t22b2t),N(v)=sin (4πvb2t)sinh(v)exp (-2v2b2t).
3 期權定價公式
定理1 股票價格服從式(1)和式(2),執行價格為K,到期日為T的歐式看漲期權在0時刻的價格V0為:V0=C(0,S)=∫
SymboleB@ 0SΦ(12σ2(0)x-ln(KS)+∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx
-Ke-∫T0r(s)ds∫
SymboleB@ 0Φ(-12σ2(0)x-ln (KS)+∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx.
證明由期權定價的鞅方法,執行價格為K,到期日為T的歐式看漲期權在0時刻的價格為
V0=C(0,S)=EQ[e-∫T0r(s)ds(ST-K)+]
=EQ[e-∫T0r(s)dsSTI{ST>K}]-Ke-∫T0r(s)dsEQ[I{ST>K}].(4)
而 ST>K∫T0σ(s)dBQs>12∫T0σ2(s)ds+ln (KS)-∫T0r(s)ds.
記Δ=∫T0σ2(s)ds,δ=12Δ+ln (KS)-∫T0r(s)ds,則
EQ[e-∫T0r(s)dsSTI{ST>K}σ(WS)]=∫+
SymboleB@ 0Sexp (-12Δ+x)12πΔexp (-x22Δ)dx=SΦ(Δ-δΔ)
由式(2)可得,σ2(t)=σ2(0)exp {(u-12ξ)t+ξdWt},從而由引理3,Δσ2(0)的密度函數為p(T,u-12ξ,ξ;x),因此,由全期望公式得
EQ[e-∫T0r(s)dsSTI{ST>K}]=EQ[EQ[e-∫T0r(s)dsSTI{ST>K}σ(WS)]]
=∫+
SymboleB@ 0SΦ(12σ2(0)x-ln (KS)+∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx.(5)
同理可得
EQ[I{ST>K}]=∫+
SymboleB@ 0Φ(-12σ2(0)x-ln (KS)+∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx. (6)
將式(5)和式(6)代入式(4)即得定理.
推論1 當ξ=0,u=0時,σ(t)=σ為常數,再假定r(t)=r為常數,可得到與定理相同條件下歐式看漲期權在0時刻的價格V0為:
V0=SΦ(12σ2T-ln (KS)+rTσT)-Ke-rTΦ(-12σ2T-ln (KS)+rTσT).
此公式即為BS公式.
推論2 股票價格服從式(1)和式(2),執行價格為K,到期日為T的歐式看跌期權在0時刻的價格P(0,S)為:
P(0,S)=EQ[e-∫T0r(s)ds(K-ST)+]
=Ke-∫T0r(s)ds∫
SymboleB@ 0Φ(12σ2(0)x+ln (KS)-∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx
-∫
SymboleB@ 0SΦ(-12σ2(0)x+ln (KS)-∫T0r(s)dsσ(0)x)p(T,u-12ξ,ξ;x)dx.
推論3 (歐式看漲期權與看跌期權的平價關系)執行價格為K,到期日為T的歐式看漲期權與看跌期權在0時刻的價格之間有關系式:
C(0,S)-P(0,S)=S-Ke-∫T0r(s)ds.
參考文獻
[1] BIACK F,SCHOLES M. The pricing of options and corporate liabilities [J]. Political Economy,1973,81(3):637-659.
[2] SCOTTL O. Option pricing when the variance changes randomly: theory, estimation, and application [J]. Journal of Financial and Quantitative Analysis,1987,22:419-438.
[3] WIGGINS J.B. Option values under stochastic volatility [J]. Journal of Financial Economics,1987,15:351-372.
[4] RENAULT E, TOUZI N. Optional hedging and implied volatilities in a stochastic volatility model [J]. Mathematical Finance,1996,6:279-302.
[5] HULL J A,WHITE A. The pricing of options on assets with stochastic volatilities[J]. The Journal of Finance,1987,42:281-300.
[6] 陳俊霞, 蹇明.隨機波動率情形下期權定價的解析解[J].統計與決策,2007,(8):21-22.
[7] 林建華, 王福昌,馮敬海.股價波動的指數OU過程模型[J].經濟數學,2000,17(4):29-32.
[8] FOLLMER H, SCHWEIZERr M. Hedging of contingent claims under incomplete information [M].New York: Gordon and Breach, 1991.
[9] LIPTSER R S,A N SHIRYAYEU.Statistics of random process[M].New York :Springverlag ,1977.
[10]YANG Z J, HUANG L H.MA C Q. Explicit expression for the valuation and hedging of the arithmetic asian option [J]. Journal of Systems and Complexity,2003,16:557-561.
Option Pricing under Stochastic Volatility and Stock PriceDirven by Ornsteinhlenback Process
ZHU lizhi,MA Peng,YU Junwu
(College of Mathematics and Computational Science,Hunan university of Science and Technology,Xiangtan,Hunan 411201,China)
Abstract Undertheassumption that stock price process is driven by OrnsteinUhlenback process,the pricing formula of European option with Stochastic Volatility was derived by using martingale method,which generalizes theBS model.
Keywordsoptionpricing ;stochastic volatility;OrnsteinUhlenback process