从布朗运动到Black–Scholes
我們先看Shreve給出的布朗運(yùn)動(dòng)的定義
簡(jiǎn)單來(lái)說(shuō)布朗運(yùn)動(dòng)就是一個(gè)“無(wú)記憶”的隨機(jī)過(guò)程。每次的增量服從正態(tài)分布,這個(gè)增量與之前的路徑無(wú)關(guān),增量的variance等于這個(gè)增量對(duì)應(yīng)時(shí)間差(Var[W(t+Δt)?W(t)]=ΔtVar[W(t+\Delta t)-W(t)]=\Delta tVar[W(t+Δt)?W(t)]=Δt).
現(xiàn)在我們考慮更為復(fù)雜的隨機(jī)過(guò)程。在微積分中,我們會(huì)有微分的概念dx=f′(t)dt(x=f(t),dt=t+dt?t)dx = f'(t)dt (x = f(t), dt = t+dt - t)dx=f′(t)dt(x=f(t),dt=t+dt?t). 如果我們簡(jiǎn)單的有dx=μdt,x(0)=0dx=\mu dt, x(0)=0dx=μdt,x(0)=0, 我們就可以推出x(t)=μtx(t) = \mu tx(t)=μt. 我們現(xiàn)在運(yùn)用剛學(xué)的布朗運(yùn)動(dòng),在原有的微分的基礎(chǔ)上加入隨機(jī)擾動(dòng) dW=W(t+dt)?W(t)dW = W(t+dt)-W(t)dW=W(t+dt)?W(t). 用我們剛學(xué)得關(guān)于布朗運(yùn)動(dòng)的性質(zhì),這個(gè)隨機(jī)繞動(dòng) dWdWdW 就是一個(gè)均值為0,variance為 dtdtdt 的正態(tài)隨機(jī)量。在微分中加入隨機(jī)擾動(dòng) σdW\sigma dWσdW, 這樣我們得到 dx=μdt+σdWdx = \mu dt + \sigma dWdx=μdt+σdW. 此時(shí) x(t)x(t)x(t) 就不在是一個(gè)確定的函數(shù)而變成了一個(gè)隨機(jī)變量。我們可以類(lèi)似的求得x(t)=μt+σW(t)x(t) = \mu t + \sigma W(t)x(t)=μt+σW(t).
現(xiàn)在我們來(lái)看一個(gè)略微復(fù)雜的隨機(jī)過(guò)程,dS=μS dt+σS dWdS = \mu S \,dt + \sigma S \,dWdS=μSdt+σSdW. 這是Black–Scholes模型中假設(shè)股票SSS服從的隨機(jī)過(guò)程。為了理解這個(gè)過(guò)程,我們來(lái)看一下Ito公式。假設(shè)函數(shù)C(t,x)C(t,x)C(t,x)有連續(xù)的二階導(dǎo)數(shù),那么我們有
d C(t,S)=Ct dt+CS dS+12CSS (dS)2=Ct dt+CS(μS dt+σS dW)+12CSS (μS dt+σS dW)2=Ct dt+CS(μS dt+σS dW)+12CSSσ2S2 dt=(Ct+μSCS+σ22S2CSS) dt+σSCS dW.d\,C(t,S) = C_t\,dt + C_S\,dS + \frac{1}{2} C_{SS}\,(dS)^2 = C_t\,dt + C_S (\mu S\,dt + \sigma S\,dW) + \frac{1}{2} C_{SS}\,(\mu S\,dt + \sigma S\,dW)^2 = C_t\,dt + C_S(\mu S\,dt + \sigma S\,dW) + \frac{1}{2} C_{SS}\sigma^2 S^2\,dt=(C_t + \mu S C_S + \frac{\sigma^2}{2} S^2 C_{SS})\,dt + \sigma S C_S \, dW.dC(t,S)=Ct?dt+CS?dS+21?CSS?(dS)2=Ct?dt+CS?(μSdt+σSdW)+21?CSS?(μSdt+σSdW)2=Ct?dt+CS?(μSdt+σSdW)+21?CSS?σ2S2dt=(Ct?+μSCS?+2σ2?S2CSS?)dt+σSCS?dW.
理解Ito公式的一個(gè)簡(jiǎn)單方法就是用泰勒展開(kāi)。但我們必須記住保留二階微分項(xiàng)(dS)2(dS)^2(dS)2因?yàn)?span id="vt6mr5x" class="katex--inline">(dW)2=dt(dW)^2=dt(dW)2=dt.通過(guò)Ito公式,我們可以得到
d(lnS)=(μ?σ22)dt+σdWd(ln S) = (\mu-\frac{\sigma^2}{2})dt+\sigma dWd(lnS)=(μ?2σ2?)dt+σdW
因此
S=e(μ?σ22)t+σW(t).S = e^{(\mu-\frac{\sigma^2}{2})t+\sigma W(t)}.S=e(μ?2σ2?)t+σW(t).
Excel for the simulation of the above: https://github.com/mathengineer/BrownianMotionToAmericanOption
從Ito公式,我們很容易能得到Black-Scholes PDE. 我們考慮一個(gè)歐式看漲期權(quán),期權(quán)的價(jià)格應(yīng)該與時(shí)間和股票價(jià)格有關(guān)。我們假設(shè)這個(gè)期權(quán)在 ttt 時(shí)刻的價(jià)格是 C(t,S)C(t, S)C(t,S). 我們現(xiàn)在來(lái)看一個(gè)這個(gè)期權(quán)加上 ?CS-C_S?CS? share 股票的一個(gè)投資組合。當(dāng)時(shí)間變化了 dtdtdt 這個(gè)投資組合的變化量應(yīng)該為 dC(t,S)?CSdSdC(t, S) - C_S dSdC(t,S)?CS?dS. 應(yīng)用 Ito 公式,我們可以算出這個(gè)變化量為 dC(t,S)?CSdS=(Ct+σ22S2CSS) dtdC(t, S) - C_S dS = (C_t + \frac{\sigma^2}{2} S^2 C_{SS})\,dtdC(t,S)?CS?dS=(Ct?+2σ2?S2CSS?)dt. 可以看到這個(gè)投資組合并無(wú)任何隨機(jī)項(xiàng) (dWdWdW項(xiàng))。對(duì)于這樣一個(gè)無(wú)風(fēng)險(xiǎn)的投資組合,它的回報(bào)率只能為無(wú)風(fēng)險(xiǎn)利率 rrr, 所以 dC(t,S)?CSdS=r(C?CSS)dtdC(t, S) - C_S dS = r (C-C_S S)dtdC(t,S)?CS?dS=r(C?CS?S)dt. 結(jié)合上面的結(jié)果,我們推導(dǎo)出 (Ct+σ22S2CSS) dt=r(C?CSS) dt(C_t + \frac{\sigma^2}{2} S^2 C_{SS})\,dt = r (C - C_S S)\,dt(Ct?+2σ2?S2CSS?)dt=r(C?CS?S)dt. 這就是Black-Scholes PDE. 我們可以做變量代換把這個(gè)PDE轉(zhuǎn)化成熱方程從而得出期權(quán)的公式。這里我們不繼續(xù)深入討論這個(gè)做法。下面我們要介紹用中性(risk neutral)概率的方法來(lái)得到期權(quán)公式。概率的做法的好處是我們極其容易將它轉(zhuǎn)換成一套可以用蒙特卡洛方法的到期權(quán)價(jià)格的算法。
我們先來(lái)引入一個(gè)概率的概念鞅(martingale): 簡(jiǎn)而言之,鞅XXX就是滿(mǎn)足條件期望不增不減的隨機(jī)過(guò)程 (E[X(t)∣F(s)]=X(s)E[X(t)|\mathcal{F}(s)] = X(s)E[X(t)∣F(s)]=X(s) for all s<ts < ts<t). 如果期權(quán)價(jià)格C(t,S)C(t,S)C(t,S)是鞅的話(huà),期權(quán)定價(jià)將變得極其容易C(t,S(t))=E[C(T,S(T))∣F(t)]C(t,S(t))=E[C(T,S(T))|\mathcal{F}(t)]C(t,S(t))=E[C(T,S(T))∣F(t)], 而期權(quán)在最終時(shí)刻的價(jià)格就是C(T,S(T))=(S(T)?K)+C(T,S(T))=(S(T)-K)^+C(T,S(T))=(S(T)?K)+. 不幸的是期權(quán)價(jià)格C(t,S)C(t,S)C(t,S)并不是鞅。幸運(yùn)的是,它幾乎等于鞅。我們從上面的推導(dǎo)看到期權(quán)加上 ?CS-C_S?CS? share 股票的投資組合的回報(bào)率等于rrr. 這樣我們很容易可以看出期權(quán)等同于CSC_SCS? share 股票加上把剩余的前存進(jìn)無(wú)風(fēng)險(xiǎn)銀行(或者從銀行借錢(qián))的組合(CSC_SCS? share of S and the rest of money in money account M where dM=rMdtdM=rMdtdM=rMdt). 從股票的公式dS=μS dt+σS dWdS = \mu S \,dt + \sigma S \,dWdS=μSdt+σSdW我們看到它的平均回報(bào)率是 μ\muμ,這個(gè)回報(bào)率一般高于rrr來(lái)reward投資者接受的風(fēng)險(xiǎn)σ\sigmaσ (μ?rσ\frac{\mu-r}{\sigma}σμ?r?就是夏普比Sharpe ratio)。Girsanov定理告訴我們,我們可以通過(guò)改變概率分布來(lái)改變隨機(jī)過(guò)程的drift(平均回報(bào)率項(xiàng))。我們現(xiàn)在取一個(gè)中性概率空間,使得股票在此概率空間中的平均回報(bào)率剛好等于rrr,即dS=rS dt+σS dW~dS = r S \,dt + \sigma S \,d\tilde{W}dS=rSdt+σSdW~. 在此空間中,我們剛才的股票加銀行賬戶(hù)的平均回報(bào)率就是rrr, 因而期權(quán)的平均回報(bào)率也就等于rrr. 在此概率空間,我們簡(jiǎn)單的有dC=rCdt+...dW~dC = rC dt + ... d\tilde{W}dC=rCdt+...dW~. 于是de?rtC=...dW~de^{-rt}C=...d\tilde{W}de?rtC=...dW~. 既然我們只有布朗運(yùn)動(dòng)項(xiàng),e?rtCe^{-rt}Ce?rtC這個(gè)過(guò)程就必定是鞅。所以我們就有期權(quán)在中性概率下的定價(jià)公式e?rtC(t,S(t))=E[e?rTC(T,S(T))∣S(t)]e^{-rt}C(t, S(t)) = E[e^{-rT}C(T, S(T))|S(t)]e?rtC(t,S(t))=E[e?rTC(T,S(T))∣S(t)].
未完待續(xù)。接下來(lái)我們將講如何用上面的定價(jià)公式和蒙特卡洛方法來(lái)算歐式期權(quán)價(jià)格。
總結(jié)
以上是生活随笔為你收集整理的从布朗运动到Black–Scholes的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 11行代码实现滚动公告栏
- 下一篇: php公告滚动源码,如何利用JavaSc