日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

从布朗运动到Black–Scholes

發布時間:2023/12/16 编程问答 42 豆豆
生活随笔 收集整理的這篇文章主要介紹了 从布朗运动到Black–Scholes 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

我們先看Shreve給出的布朗運動的定義

簡單來說布朗運動就是一個“無記憶”的隨機過程。每次的增量服從正態分布,這個增量與之前的路徑無關,增量的variance等于這個增量對應時間差(Var[W(t+Δt)?W(t)]=ΔtVar[W(t+\Delta t)-W(t)]=\Delta tVar[W(t+Δt)?W(t)]=Δt).

現在我們考慮更為復雜的隨機過程。在微積分中,我們會有微分的概念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). 如果我們簡單的有dx=μdt,x(0)=0dx=\mu dt, x(0)=0dx=μdt,x(0)=0, 我們就可以推出x(t)=μtx(t) = \mu tx(t)=μt. 我們現在運用剛學的布朗運動,在原有的微分的基礎上加入隨機擾動 dW=W(t+dt)?W(t)dW = W(t+dt)-W(t)dW=W(t+dt)?W(t). 用我們剛學得關于布朗運動的性質,這個隨機繞動 dWdWdW 就是一個均值為0,variance為 dtdtdt 的正態隨機量。在微分中加入隨機擾動 σdW\sigma dWσdW, 這樣我們得到 dx=μdt+σdWdx = \mu dt + \sigma dWdx=μdt+σdW. 此時 x(t)x(t)x(t) 就不在是一個確定的函數而變成了一個隨機變量。我們可以類似的求得x(t)=μt+σW(t)x(t) = \mu t + \sigma W(t)x(t)=μt+σW(t).

現在我們來看一個略微復雜的隨機過程,dS=μS dt+σS dWdS = \mu S \,dt + \sigma S \,dWdS=μSdt+σSdW. 這是Black–Scholes模型中假設股票SSS服從的隨機過程。為了理解這個過程,我們來看一下Ito公式。假設函數C(t,x)C(t,x)C(t,x)有連續的二階導數,那么我們有
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公式的一個簡單方法就是用泰勒展開。但我們必須記住保留二階微分項(dS)2(dS)^2(dS)2因為(dW)2=dt(dW)^2=dt(dW)2=dt.通過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. 我們考慮一個歐式看漲期權,期權的價格應該與時間和股票價格有關。我們假設這個期權在 ttt 時刻的價格是 C(t,S)C(t, S)C(t,S). 我們現在來看一個這個期權加上 ?CS-C_S?CS? share 股票的一個投資組合。當時間變化了 dtdtdt 這個投資組合的變化量應該為 dC(t,S)?CSdSdC(t, S) - C_S dSdC(t,S)?CS?dS. 應用 Ito 公式,我們可以算出這個變化量為 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. 可以看到這個投資組合并無任何隨機項 (dWdWdW項)。對于這樣一個無風險的投資組合,它的回報率只能為無風險利率 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. 結合上面的結果,我們推導出 (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. 我們可以做變量代換把這個PDE轉化成熱方程從而得出期權的公式。這里我們不繼續深入討論這個做法。下面我們要介紹用中性(risk neutral)概率的方法來得到期權公式。概率的做法的好處是我們極其容易將它轉換成一套可以用蒙特卡洛方法的到期權價格的算法。

我們先來引入一個概率的概念鞅(martingale): 簡而言之,鞅XXX就是滿足條件期望不增不減的隨機過程 (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&lt;ts &lt; ts<t). 如果期權價格C(t,S)C(t,S)C(t,S)是鞅的話,期權定價將變得極其容易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)], 而期權在最終時刻的價格就是C(T,S(T))=(S(T)?K)+C(T,S(T))=(S(T)-K)^+C(T,S(T))=(S(T)?K)+. 不幸的是期權價格C(t,S)C(t,S)C(t,S)并不是鞅。幸運的是,它幾乎等于鞅。我們從上面的推導看到期權加上 ?CS-C_S?CS? share 股票的投資組合的回報率等于rrr. 這樣我們很容易可以看出期權等同于CSC_SCS? share 股票加上把剩余的前存進無風險銀行(或者從銀行借錢)的組合(CSC_SCS? share of S and the rest of money in money account M where dM=rMdtdM=rMdtdM=rMdt). 從股票的公式dS=μS&ThinSpace;dt+σS&ThinSpace;dWdS = \mu S \,dt + \sigma S \,dWdS=μSdt+σSdW我們看到它的平均回報率是 μ\muμ,這個回報率一般高于rrr來reward投資者接受的風險σ\sigmaσ (μ?rσ\frac{\mu-r}{\sigma}σμ?r?就是夏普比Sharpe ratio)。Girsanov定理告訴我們,我們可以通過改變概率分布來改變隨機過程的drift(平均回報率項)。我們現在取一個中性概率空間,使得股票在此概率空間中的平均回報率剛好等于rrr,即dS=rS&ThinSpace;dt+σS&ThinSpace;dW~dS = r S \,dt + \sigma S \,d\tilde{W}dS=rSdt+σSdW~. 在此空間中,我們剛才的股票加銀行賬戶的平均回報率就是rrr, 因而期權的平均回報率也就等于rrr. 在此概率空間,我們簡單的有dC=rCdt+...dW~dC = rC dt + ... d\tilde{W}dC=rCdt+...dW~. 于是de?rtC=...dW~de^{-rt}C=...d\tilde{W}de?rtC=...dW~. 既然我們只有布朗運動項,e?rtCe^{-rt}Ce?rtC這個過程就必定是鞅。所以我們就有期權在中性概率下的定價公式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)].

未完待續。接下來我們將講如何用上面的定價公式和蒙特卡洛方法來算歐式期權價格。

總結

以上是生活随笔為你收集整理的从布朗运动到Black–Scholes的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。