摘 要 運用混合分數布朗運動的It公式,將幾何平均亞式期權定價化成一個偏微分方程求解問題,通過偏微分方程求解獲得了幾何平均型亞式看漲期權的定價公式.
關鍵詞 混合分數布朗運動; 幾何平均型亞式期權; Black-Scholes偏微分方程
中圖分類號 O211.6;F830.9 文獻標識碼 A
Asian Option Pricing Model in Mixed Fractional Brownian Motion Environment
SUN Yu-dong, SHI Yi-min
(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an,Shannxi 710072,China)
AbstractIn order to overcome the shortcomings, the geometric average Asian option was changed into the question of solving partial differential equation by mixed It formula. The pricing formula of the geometric average Asian call and put option were obtained by partial differential equation theory.
Key words mixed fractional Brownian motion; geometric average Asian option; Black-Scholes partial differential equation
1 引 言
亞式期權作為一種路徑依賴期權,是一種新型期權,其執行與否取決于合同期內股票平均價格的高低.由于采用平均值可以減少價格波動的所帶來的影響,從而亞式期權比類似的常規期權更便宜,因此在貨幣和商品市場中比較流行.亞式期權是一種依賴標的資產價格路徑的期權, 它在到期日的收益依賴于期權整個有效期內標的資產的平均價格,根據取平均的方法亞式期權可以分為算術平均型和幾何平均型兩種類型.筆者主要考慮幾何平均亞式期權.文獻[1-2] 用偏微分方程的方法,研究了幾何平均型亞式期權定價模型,在求解方程時,為了使得方程能夠求解,假定所涉及的參數均為常數;文獻[3-4] 在風險中性測度下采用鞅方法,討論了幾何平均亞式期權定價問題.但是上訴文獻都是布朗運動下的結論,近些年來,在實際金融市場中由于股票價格對過去價格具有依賴性;一些學者考慮用分數布朗運動刻畫股票價格的變化[5],甚至有部分學者開始用混合布朗運動研究歐式期權[6] .本文假定股票價格遵循分數驅動的隨機微分方程,利用分數布朗運動隨機分析理論,得到了混合分數布朗運動環境下幾何平均型亞式期權價格所滿足的Black-Scholes偏微分方程, 并通過求解該偏微分方程給出了亞式看漲期權的定價公式.
2 金融市場數學模型
考慮存在兩種可自由、連續進行交易的資產的無套利和無摩擦市場, Bt表示t時刻無風險資產的價格,St表示風險資產的價格,如股票.本文假設它們在物理測度下仍滿足隨機微分方程
dBt=rtBtdt,(1)
dSt=St((μt-qt)dt+σ1dWHt+σ2dWt),(2)
其中參數rt,μt都為時間的函數,他們分別表示銀行短期利率、股票期望回報率,波動項系數σ1和σ2為非負常數,Wt是布朗運動, WHt為參數為H的分數布朗運動,Wt是與WHt相互獨立.這里假定作為原生資產的股票支付紅利,紅利率為qt. 設在t時刻形成投資組合θ=(θ0t,θ1t), 則投資組合的回報為Vt=θ0tBt+θ1tSt,并且由自融資策略(即短時間不改變市場份額)可知V滿足
dVt=θ1tdSt+θ0tdBt+θ1tqtStdt.
令Jt=exp {1t∫t0ln Sτdτ}表示股票價格的路徑變量, 則路徑變量Jt對t求導可得
dJtdt=Jt(ln St-ln Jtt).(3)
定理1 股票價格所滿足隨機微分方程(2)的解為
St=S0exp {∫t0(μτ-qτ)dτ-σ22t-σ22t2H+σ1dWHt+σ2dWt}. (4)
證明 記f(t,x,y)=S0exp {∫t0(μτ-qτ)dτ-λσ22t-σ22t2H+σ1x+σ2y}, 由混合分數布朗運動It公式[6]可得
dSt=df(t,WHt,Wt)=ft(t,WHt,Wt)dt+
fx(t,WHt,Wt)dWHt+fy(t,WHt,Wt)dWt+
Ht2H-1fxx(t,WHt,Wt)dt+12fyy(t,WHt,Wt)dt
=(μt-qt-Hσ21t2H-1-0.5σ22)f(t,WHt,Wt)dt+
σ1f(t,WHt,Wt)dWHt+σ2f(t,WHt,Wt)dWt+
Hσ21t2H-1f(t,WHt,Wt)dt+0.5σ22f(t,WHt,
Wt)dt=f(t,WHt,Wt)((μt-qt)dt+σ1dWHt+
σ2dWt)=St((μt-qt)dt+σ1dWHt+σ2dWt).
證畢.
定理2 混合分數布朗環境下幾何平均亞式期權定價可以化為偏微分方程定解問題:
Ft+(rt-qt)xFx+y(ln x-ln yt)Fy+
(Hσ21t2H-1+σ222)x22Fx2-rtF=0,F(T,x,y)=f(x,y),(t,x,y)∈[0,T]×
R+×R+.
證明 像其他的普通的歐式未定權益一樣,幾何平均亞式看漲期權都時間和股票價格St,但是它作為一個強路徑依賴期權,與普通的歐式期權有著很大的不同,幾何平均亞式看漲期權在當前時刻的價值還有它過去的價格平均值有關,即它還依賴路徑變量Jt.因此,未定權益的價格過程可表示為F=F(t,St,Jt). 由式(4)和混合分數布朗運動It公式有
dF(t,St,Jt)=(Ft(t,St,Jt)+(μt-qt-λσ22-
Hσ21t2H)StFx(t,St,Jt))dt+lnSt-lnJttJt#8226;
Fy(t,St,Jt)dt+σ1StFx(t,St,Jt)dWHt+
σ2StFx(t,St,Jt)dWt+
(Hσ21t2H-1+
σ222)(S2t2Fx2(t,St,Jt)+StFx(t,St,Jt))dt
=(Ft(t,St,Jt)+(μt-
qt)StFx(t,St,Jt)+
lnSt-lnJttJtFy(t,St,Jt))dt+
(Hσ21t2H-1+σ222)S2t2Fx2(t,St,Jt)dt+
σ1StFx(t,St,Jt)dWHt+σ2StFx(t,St,Jt)dWt.
(5)
因為
dVt=θ0tdBt+θ1tdSt+qtθ1tStdt
=θ0trtBtdt+θ1t(μtStdt+σ1StdWHt+σ2StdWt)
=(θ0trtBt+θ1tμtSt)dt+θ1t(σ1StdWHt+σ2StdWt). (6)
那么根據自融資交易策略Vt=F(t,St,Jt)以及式(5)、式(6)可得
θ1t=Fx(t,St,Jt),
θ0t=(rtBt)-1(Ft(t,St,Jt)-qtStFx(t,St,Jt)+
ln x-ln JttJtFx(t,St,Jt)+(Hσ21t2H-1+
σ222)S2t2Fx2(t,St,Jt)).
再由Vt=F(t,St,Jt)=θ0tBt+θ1tSt知θ0t=B-1t(F(t,St,Jt)-StFx(t,St,Jt)),因此得二階偏微分方程
Ft+(rt-qt)xFx+ln x-ln ytyFy+
(Hσ21t2H-1+σ222)x22Fx2-rtF=0,
由F(T,ST,JT)=f(ST,JT)知,F(t,x,y)=f(x,y),(t,x,y)∈[0,T]×R+×R+, 所以有結論成立.
3 模型求解
定理3 幾何平均亞式看漲期權在t時刻價格可以表示為
F(t,St,Jt)=(JttST-tt)1/Texp{r*(T-t)-
∫Ttr(τ)dτ+(σ*H)22(T2H-t2H)+(σ*)22(T-
t)}Φ(d1)-Kexp{-∫Ttr(τ)dτ}Φ(d2),
其中
d1=
1TlnJttST-ttKT+r*(T-t)+(σ*)2(T-t)+(σ*H)2(T2H-t2H)(σ*)2(T-t)+(σ*H)2(T2H-t2H),
d2=d1-(σ*)2(T-t)+(σ*H)2(T2H-t2H),
Jt=exp{1t∫t0lnSτdτ}r*=∫Tt(rτ-qτ)T-τTdτT-t-
σ22(T-t)4T-σ21(T2H-t2H)2(T-t)+Hσ21(T2H+1-t2H+1)(2H+1)T(T-t),
σ*H=[1-4H(T2H+1-t2H+1)(2H+1)T(T2H-t2H)+
H(T2H+2-t2H+2)(H+1)T2(T2H-t2H)]1/2σ1,σ*=T-t3Tσ2.
證明 因為此時考慮具有固定執行價格的幾何平均亞式期權,所以取f(x,y)=(y-K)+,令
ξ=tln y+(T-t)ln xT,F=U(t,ξ),
則
Ft=Ut+ln y-ln xTUξ,
Fx=T-tTxUξ,
2Fx2=(T-tTx)22Uξ2-T-tTx2Uξ,
Fy=tTyUξ.
從而混合分數布朗運動環境下Black-Scholes偏微分方程化為
Ut+(Hσ21t2H-1+λσ222)(FtT)22Uξ2+(rt-
qt-Hσ21t2H-1-σ222)(T-tT)Uξ-rtU=0U(T,ξ)=(eξ-K)+.
再令
C(s,z)=U(t,ξ)exp {∫Ttrτdτ},z=ξ+∫Tt(rτ-
qτ)T-τTdτ-σ224T(T-t)2-σ212(T2H-
t2H)+Hσ21T(2H+1)(T2H+1-t2H+1),
s=σ226T2(T-t)3-σ212(T2H-t2H)+
2Hσ21T(2H+1)(T2H+1-t2H+1)-
Hσ212(H+1)T2(T2H+2-t2H+2),
則分數跳擴散環境下Black-Scholes偏微分方程可化為
Cs=2Cz2,C0,z=(ez-K).
根據熱傳導方程經典解理論,其唯一強解C(s,z)=12πs∫
-
u0(τ)e(τ-z)24sdτ, 將邊界條件帶入,可得
C(s,z)=ez+sΦ(z+2s-ln K2s)-
KΦ(z-ln K2s),
對上述變換進行逆變換可得定理證明.
定理4 當σ2=0時可得分數布朗運動環境下幾何平均亞式看漲期權在t時刻價格
C(t,St,Jt)=(JttST-tt)1/Texp{r*(T-t)-
∫Ttr(τ)dτ+(σ*H}22(T2H-t2H)}Φ(d1)-
Kexp{-∫Ttr(τ)dτ}Φ(d2),
其中
d1=1TlnJttST-ttKT+r*(T-t)+(σ*H)2(T2H-t2H)(σ*)2(T-t)+(σ*H)2(T2H-t2H),
d2=d1-(σ*H)2(T2H-t2H),
Jt=exp{1t∫t0ln Sτdτ},r*=∫Tt(rτ-qτ)T-τTdτT-t-
σ21(T2H-t2H)2(T-t)+Hσ21(T2H+1-t2H+1)(2H+1)T(T-t),
σ*H=[1-4H(T2H+1-t2H+1)(2H+1)T(T2H-t2H)+
H(T2H+2-t2H+2)(H+1)T2(T2H-t2H)]1/2σ1.
特別地,當H=0.5可得文獻[1-4]中的結果.
參考文獻
[1] 王志明,朱芳芳.幾何平均型亞式期權定價模型的新解法[J].武漢科技大學學報:自然科學版,2008, 32(2): 321-322.
[2] 姜禮尚. 期權定價的數學模型和方法[M].2版. 北京:高等教育出版社,2008.
[3] 羅慶紅, 楊向群. 幾何型亞式期權定價研究[J].湖南文理學院學報,2007, 19(1): 5-7.
[4] 王國強, 馬德全, 宋華. 亞式期權的一種定價方法[J].數學理論與應用, 2007, 27(3): 410-411.
[5] 薛紅, 王拉省. 分數布朗運動環境中最值期權定價[J].工程數學學報, 2008, 25(5): 843-850.
[6] R J Elliott,JHoek.A general fractional white noise theory and applications to finance[J]. Mathematical Finance, 2003, 13 (2): 301-330.
注:本文中所涉及到的圖表、注解、公式等內容請以PDF格式閱讀原文