[衍生商品] 淺談 Black-Scholes Model 的性質 (0)

這次要跟大家介紹衍生商品市場的 Black-Scholes Model (B-S model)，此 Formula 是由 Professor Fisher Black, Myron Scholes 與 Robert Merton 在選擇權定價領域的重大突破，此模型亦指引了衍生商品該如何透過無套利機會來獲得合理價格的研究大門。

事實上， B-S model 本質上是用來作為進行 歐式選擇權 (European Option) 定價公式。

在介紹之前，我們需先知道使用 B-S model 的一些假設：

對於股價分布的假設

1. 股價服從連續複利 Log-normal 分布
2. 波動度(Volatility)為已知常數
3. 未來的股息已知

對於市場的假設

1. 無交易手續費用、無稅收
2. 證卷交易為連續進行
3. 短期無風險利率 $r$ 為已知常數
4. 可以基於無風險利率執行 short sell or borrow
5. 不存在無風險套利機會 

Comment:
1. 股價的 Log-normal 分布假設： (此 Comment 需要隨機過程與隨機分析的背景知識，有興趣的讀者請參考BLOG中相關隨機分析的文章)

考慮股價模型為 Geometric Brownian Motion，故可寫做下列隨機微分方程 SDE

$$dS_t := \mu S_t dt + \sigma S_t dB_t, \ 0 \leq t \leq T$$ 其中 $S_t$ 為時刻 $t$ 的股價， $\mu :=$ 股價每年的收益率期望值 (或稱 drift rate)，$\sigma:=$ 股價每年的波動度 (volatility)，$B_t$ 為標準布朗運動。注意到這邊我們假設 $\mu, \sigma$ 為固定常數。

上述隨機微分方程可解得 (Proof: omitted，有興趣讀者請參閱
[隨機分析] How to solve SDE practically (4) - Geometric Brownian Motion
$$S_t = S_0 \exp \left\{ {\left( {\mu  - \frac{{{\sigma ^2}}}{2}} \right)t + \sigma {B_t}} \right\}$$ 改寫上式
$$\Rightarrow \ln \frac{{{S_t}}}{{{S_0}}} = \left( {\mu  - \frac{{{\sigma ^2}}}{2}} \right)t + \sigma {B_t}$$ 由於 $B_t$ 為標準布朗運動，由定義可知標準布朗運動為 Gaussian process with mean 0, variance $t$ 與 covariance $\min{(s,t)}$，故我們可推論

$$\begin{array}{l} \ln \frac{{{S_T}}}{{{S_0}}}\sim {{\cal N}}((\mu  - \frac{{{\sigma ^2}}}{2})T,{\sigma ^2}T)\\  \Rightarrow \ln {S_T}\sim {{\cal N}}(\ln{S_0} + (\mu  - \frac{{{\sigma ^2}}}{2})T,{\sigma ^2}T) \end{array}$$ 其中 $S_T$ 是未來時間 $T$ 時的股價，$S_0$ 是時間 $0$ 時的股價。上式表明 $\ln {S_T}$ 服從 normal distribution，故 $S_T$ 為Log-normal (亦即 取Log之後為normal)，

2. 關於 波動度 (Volatility) $\sigma$
股價的波動度 $\sigma$ 用於測量股價收益的不確定性。一般而言介於 $15 \% \sim 50 \%$。
B-S formula 假設：給定任意 $K,S,T,r..$， $\sigma$ 均為常數不變。

另外值得一提的是波動度有兩種，一種是 歷史資料波動度 (Historical Volatility) 與 隱含波動度 (Implied Volatility)，其中歷史波動度是由歷史資料股價計算收益再由此歷史收益計算標準差將其定為波動度。但是隱含波動度則是透過 Black-Scholes Formula 反推而得。此波動度會在之後再作介紹。

3. 關於 Black-Scholes Formula 本質：
透過 購買/賣出 股票與債卷  來複製 選擇權的收益。透過此法可建構一組無風險投資組合，且此組合收益僅為無風險利率 (亦及無套利機會)，相關推導請參閱此文： [隨機分析] Black-Scholes PDE for European Call option (0)

Black-Scholes-Metron Formula: 

B-S Model For Stock Option with Constant Dividend: 
以下為  對 European Call option 與 European put option 的 B-S formula

$$\left\{ \begin{array}{l} c = {S_0}{e^{ - qT}}N\left( {{d_1}} \right) - K{e^{ - rT}}N\left( {{d_2}} \right)\\ p = K{e^{ - rT}}N\left( {{-d_2}} \right) - {S_0}{e^{ - qT}}N\left( {{-d_1}} \right) \end{array} \right.$$ 其中 $c$ 為 European Call option 的價格、$p$ 為 European Put option 價格、$q$ 為連續複利的固定股息、$r$ 為連續複利的無風險利率、 $S_0$ 為現時股價、$T$ 到期時間(以年為單位)、$K$ 為執行價格、$\sigma$ 為股價波動度；且 $N(\cdot)$ 為 Standard Normal cumulative distirbuiton function，其中

$$\left\{ \begin{array}{l} {d_1} = \frac{{\ln \left( {{S_0}/K} \right) + \left( {r - q + {\sigma ^2}/2} \right)T}}{{\sigma \sqrt T }}\\ {d_2} = {d_1} - \sigma \sqrt T \end{array} \right.$$ 關於上式的證明有興趣的讀者請參考
 [隨機分析] Black-Scholes PDE for European Call option (1)

如果是要計算的話，讀者可以使用 MATLAB 指令 blsprice 協助計算，在此不贅述

Comment:
1. 仔細觀察 B-S Formula，我們可發現需要的參數有
$S_0, K, T, r, q, \sigma$ ，亦即選擇權 $f$ 為
$$f(S,K,T,r,q,\sigma)$$ 注意到其中只有 $\sigma$ 無法被直接觀測獲得，其餘參數都可直接獲得。但是事實上選擇權的報價是知道的 $c,p$已知，故我們可以透過代入選擇權價格以及其餘已知資訊來反推 $\sigma$，此法求得的 $\sigma$ 即為之前所提及的隱含波動度( Implied Volatility)

2. 關於上述 B-S formula的性質：
首先考慮 $S_0$ 很大 (deep in the money for call option, and out of money for put option) 的時候，
則上述 $d_1, d_2 \rightarrow \infty$ 故 $N(d_1), N(d_2) \rightarrow 1$，$N(-d_1), N(-d_2) \rightarrow 0$亦即
$$\left\{ \begin{array}{l} c = {S_0}{e^{ - qT}}N\left( {{d_1}} \right) - K{e^{ - rT}}N\left( {{d_2}} \right) \to {S_0}{e^{ - qT}} - K{e^{ - rT}}\\ p = K{e^{ - rT}}N\left( { - {d_2}} \right) - {S_0}{e^{ - qT}}N\left( { - {d_1}} \right) \to 0 \end{array} \right.$$ 可發現此情況確實為 $c, p$ 的收益價格。

當 $S_0$ 非常小的時候 (deep in the money for put option, and out of money for call option)
則上述 $d_1, d_2 \rightarrow -\infty$ 故 $N(d_1), N(d_2) \rightarrow 0$，$N(-d_1), N(-d_2) \rightarrow 1$亦即
$$\left\{ \begin{array}{l} c = {S_0}{e^{ - qT}}N\left( {{d_1}} \right) - K{e^{ - rT}}N\left( {{d_2}} \right) \to 0\\ p = K{e^{ - rT}}N\left( { - {d_2}} \right) - {S_0}{e^{ - qT}}N\left( { - {d_1}} \right) \to K{e^{ - rT}} - {S_0}{e^{ - qT}} \end{array} \right.$$ 可發現此情況確實為 $c, p$ 的收益價格。

當 $T \rightarrow 0$ 的時候 ：
如果 $S_0 \geq K$ 則 ${{d_1}} \rightarrow \infty \Rightarrow d_2 \rightarrow \infty$，故
$$N(d_1) \rightarrow 1, N(d_2) \rightarrow 1$$ 且
$$\begin{array}{l} \left\{ {\begin{array}{*{20}{l}} {c = {S_0}{e^{ - qT}}N\left( {{d_1}} \right) - K{e^{ - rT}}N\left( {{d_2}} \right)}\\ {p = K{e^{ - rT}}N\left( { - {d_2}} \right) - {S_0}{e^{ - qT}}N\left( { - {d_1}} \right)} \end{array}} \right.  \Rightarrow \left\{ {\begin{array}{*{20}{l}} {c = {S_0} - K}\\ {p = 0} \end{array}} \right. \end{array}$$
如果 $S_0 < K$ 則 ${{d_1}} \rightarrow -\infty \Rightarrow d_2 \rightarrow -\infty$，故
$$N(d_1) \rightarrow 0, N(d_2) \rightarrow 0$$ 且
$$\begin{array}{l} \left\{ {\begin{array}{*{20}{l}} {c = {S_0}{e^{ - qT}}N\left( {{d_1}} \right) - K{e^{ - rT}}N\left( {{d_2}} \right)}\\ {p = K{e^{ - rT}}N\left( { - {d_2}} \right) - {S_0}{e^{ - qT}}N\left( { - {d_1}} \right)} \end{array}} \right.  \Rightarrow \left\{ {\begin{array}{*{20}{l}} {c = 0}\\ {p = K - {S_0}} \end{array}} \right. \end{array}$$
故總結上述兩種情況，我們得到 當 $T \rightarrow 0$ ，選擇權價格為
$$\left\{ {\begin{array}{*{20}{l}} {c = \max \left\{ {{S_0} - K,0} \right\}}\\ {p = \max \left\{ {K - {S_0},0} \right\}} \end{array}} \right.$$


3. 現在如果考慮B-S model 用到其他商品上的情況

B-S Model For Currency Option
此時 $S_0$ 為現時匯率；$q$ 為外幣利率(foreign currency interest rate)。

B-S Model For Future Option

此時 $S_0 = F_0$ 為現時期貨價格、$q = r$ 

以下我們看兩個例子：

Example 1:

使用 B-S model 來計算 European Put Option on the euro currencies 

考慮 現時匯率 $1.05/euro$、執行價格 $1.1/ euro$、歐元利率 $3.1 \%$、美金利率 $5.5 \%$、波動度 $10 \%$、到期時間 $4$個月

Solution

首先改寫已知資訊 $S=1.05, K=1.1, q=3.1\%, r=5.5\%, \sigma=10\%, T=4/12$，

接著計算 $d_1, d_2$：
$$\left\{ \begin{array}{l} {d_1} = \frac{{\ln \left( {1.05/1.1} \right) + \left( {5.5\%  - 3.1\%  + {{\left( {0.1} \right)}^2}/2} \right)\left( {4/12} \right)}}{{\left( {0.1} \right)\sqrt {(4/12)} }} =  - 0.6383\\ {d_2} =  - 0.6961 \end{array} \right.$$ 帶入 B-S formula 可得 $p=0.05/euro$. $\square$

Example 2:
使用 B-S formula 計算 Put Options on Futures
考慮 期貨現時價格 $\$ 31$、執行價格 $31$、利率 $2 \%$、波動度 $30 \%$、到期時間 $6$個月

Solution

首先改寫已知資訊 $S=31, K=31, q=r=2\%, \sigma=30\%, T=6/12$，(注意到在期貨中 $q=r$ !!)

接著計算 $d_1, d_2$：

$$\left\{ \begin{array}{l} {d_1} = \frac{{\ln \left( {31/31} \right) + \left( {2\%  - 2\%  + {{\left( {0.3} \right)}^2}/2} \right)\left( {6/12} \right)}}{{\left( {0.3} \right)\sqrt {(6/12)} }} =  0.1061\\ {d_2} =  - 0.1060 \end{array} \right.$$ 帶入 B-S formula 可得 Put Future Option 價格 $p= \$ 2.59$. $\square$

ref: John C. Hull, Options, Futures and Other Derivatives 7th.