这学期会时不时更新一下伊曼纽尔·德曼(Emanuel Derman) 教授与迈克尔B.米勒(Michael B. Miller)的《The Volatility Smile》这本书,本意是协助导师课程需要,发在这里有意的朋友们可以学习一下,思路不一定够清晰且由于分工原因我是从书本第13章写起,还请大家见谅。

第15章 局部波动率的影响

局部波动率的DUPIRE公式

在第11章中曾推导过Breeden-Lizenberger公式:
p ( S , t , K , T ) = e r ( T − t ) ∂ 2 C ( S , t , K , T ) ∂ K 2 p(S,t,K,T)=e^{r(T-t)}\frac{\partial^2C(S,t,K,T)}{\partial K^2} p(S,t,K,T)=er(T−t)∂K2∂2C(S,t,K,T)​
其中, C ( S , t , K , T ) C(S,t,K,T) C(S,t,K,T) 表示行权价为 K K K、到期日为 T T T 的标准看涨期权在 t t t 时刻的市场价格, ∂ 2 C / ∂ K 2 \partial^2C/\partial K^2 ∂2C/∂K2 表示价差无限小的蝶式期权组合, p ( S , t , K , T ) p(S,t,K,T) p(S,t,K,T) 表示在 t t t 时刻、标的股票价格为 S S S 时,对应的在到期日 T T T 股票价格变为 K K K 的风险中性概率密度函数

下面给出Breeden-Lizenberger公式的推导:

假设有一个欧式衍生品,标的资产是 S,到期日为 T T T,只要我们知道在 T T T 时刻所有 S T S_T ST​ 价格下伪概率 p ( S , t , S T , T ) p(S,t,S_T,T) p(S,t,ST​,T) 的值,那么我们就可以计算在到期日前的任意 t t t 时点该衍生品的价值,下面计算 p ( S , t , S T , T ) p(S,t,S_T,T) p(S,t,ST​,T) :

假设有一个标准欧式看涨期权,在到期日 T T T 的损益为 max ⁡ ( S T − K , 0 ) \max(S_T-K,0) max(ST​−K,0),其在 t t t 时点的价格表示为 C ( S , t , K , T ) C(S,t,K,T) C(S,t,K,T),于是有:
C ( S t , t , K , T ) = e − r ( T − t ) ∫ K ∞ p ( S , t , S T , T ) ( S T − K ) d S T C(S_t,t,K,T)=e^{-r(T-t)}\int_K^\infin p(S,t,S_T,T)(S_T-K)dS_T C(St​,t,K,T)=e−r(T−t)∫K∞​p(S,t,ST​,T)(ST​−K)dST​
两边对 K K K 求导可得:
∂ C ( S , t , K , T ) ∂ K = − e − r ( T − t ) ∫ K ∞ p ( S , t , S T , T ) d S T \frac{\partial C(S,t,K,T)}{\partial K}=-e^{-r(T-t)}\int_K^\infin p(S,t,S_T,T)dS_T ∂K∂C(S,t,K,T)​=−e−r(T−t)∫K∞​p(S,t,ST​,T)dST​
再次对 K K K 进行二次微分求导可得:
∂ 2 C ( S , t , K , T ) ∂ K 2 = − e − r ( T − t ) p ( S , t , K , T ) p ( S , t , K , T ) = e r ( T − t ) ∂ 2 C ( S , t , K , T ) ∂ K 2 \frac{\partial^2C(S,t,K,T)}{\partial K^2}=-e^{-r(T-t)}p(S,t,K,T)\\ p(S,t,K,T)=e^{r(T-t)}\frac{\partial^2C(S,t,K,T)}{\partial K^2} ∂K2∂2C(S,t,K,T)​=−e−r(T−t)p(S,t,K,T)p(S,t,K,T)=er(T−t)∂K2∂2C(S,t,K,T)​

Breeden-Lizenberger公式不需要任何模型假设条件

相类似的,将要推导的Dupire公式描述的就是上一章提到的局部波动率与标准期权价格对行权价 K K K 和到期日 T T T 的偏微分之间的关系

假设无风险利率和股息率等于0,Dupire公式可以表示为:
σ 2 ( K , T ) 2 = ∂ C ( S , t , K , T ) ∂ T K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 \frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C(S,t,K,T)}{\partial T}}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}} 2σ2(K,T)​=K2∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​​
这里的 σ ( K , T ) \sigma(K,T) σ(K,T) 是指在未来某时刻 T T T,股票价格为 K K K 时的局部波动率,该局部波动率对应的是更早一些的时刻 t t t,股票价格为 S S S 时的期权价格

如果无风险利率 r ≠ 0 r\neq0 r=0,那么Dupire公式可以表示为:
σ 2 ( K , T ) 2 = ∂ C ( S , t , K , T ) ∂ T + r K ∂ C ( S , t , K , T ) ∂ K K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 \frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C(S,t,K,T)}{\partial T}+rK\frac{\partial C(S,t,K,T)}{\partial K}}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}} 2σ2(K,T)​=K2∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​+rK∂K∂C(S,t,K,T)​​

理解公式

考虑无风险利率等于0的情况,此时
p ( S , t , K , T ) = ∂ 2 C ( S , t , K , T ) ∂ K 2 p(S,t,K,T)=\frac{\partial^2C(S,t,K,T)}{\partial K^2} p(S,t,K,T)=∂K2∂2C(S,t,K,T)​
∂ 2 C ( S , t , K , T ) ∂ K 2 \dfrac{\partial^2C(S,t,K,T)}{\partial K^2} ∂K2∂2C(S,t,K,T)​ 就是Breeden-Lizenberger公式中的风险中性概率密度函数 p ( S , t , K , T ) p(S,t,K,T) p(S,t,K,T),这个密度函数可以看作一个价差无限小的蝶式期权组合价格极限值

同时,由于:
∂ C ( S , t , K , T ) ∂ T = lim ⁡ d T → 0 C ( S , t , K , T + d T ) − C ( S , t , K , T ) d T \frac{\partial C(S,t,K,T)}{\partial T}=\lim_{dT\to0}\frac{C(S,t,K,T+dT)-C(S,t,K,T)}{dT} ∂T∂C(S,t,K,T)​=dT→0lim​dTC(S,t,K,T+dT)−C(S,t,K,T)​
所以 ∂ C ( S , t , K , T ) ∂ T \dfrac{\partial C(S,t,K,T)}{\partial T} ∂T∂C(S,t,K,T)​ 可以看作 1 / d T 1/dT 1/dT 份时间价格无限小的跨期期权组合价格极限值

下面分析跨期期权组合的变动情况, C ( S , t , K , T + d T ) − C ( S , t , K , T ) C(S,t,K,T+dT)-C(S,t,K,T) C(S,t,K,T+dT)−C(S,t,K,T),其中在时刻 T T T,近期期权到期,此时股票价格 S = S T S=S_T S=ST​。有三种可能性: S T ≤ K , S T ≈ K , S T ≥ K S_T\leq K,S_T\approx K,S_T\geq K ST​≤K,ST​≈K,ST​≥K

  1. S T ≤ K S_T\leq K ST​≤K,在到期日近期期权的价值等于0,且对于一个无限小的时间间隔 d T dT dT,远期期权的价值也接近于0
  2. S T ≥ K S_T\geq K ST​≥K,相似的,两个期权都处于实值状态,且价值相等

在上述两种情况下,当 d T dT dT 趋近于0时,跨期期权组合的价值等于0

而当 S T = K S_T=K ST​=K 时,近期期权在到期日处于平值状态,价值等于0,而远期期权价值为正,并且局部波动率 σ ( K , T ) \sigma(K,T) σ(K,T) 会引导股票价格向实值区间变动,于是跨期期权组合价值为正。因此,在 d T → 0 dT\to0 dT→0 时,只有当 S T ≈ K S_T\approx K ST​≈K 时, 1 / d T 1/dT 1/dT 份跨期期权组合的价值为正,决定该组合价值的波动率就是局部波动率 σ ( K , T ) \sigma(K,T) σ(K,T)

我们将简短证明,间隔期间无限小的跨期期权组合的价值实际上是与局部波动率的平方 σ 2 ( K , T ) \sigma^2(K,T) σ2(K,T)​ 成比例关系。先假设这一结论成立,我们就有
C ( S , t , K , T + d T ) − C ( S , t , K , T ) ∝ p ( S , t , K , T ) σ 2 ( K , T ) d T C ( S , t , K , T + d T ) − C ( S , t , K , T ) ∝ ∂ 2 C ( S , t , K , T ) ∂ K 2 σ 2 ( K , T ) d T C ( S , t , K , T + d T ) − C ( S , t , K , T ) d T ∝ ∂ 2 C ( S , t , K , T ) ∂ K 2 σ 2 ( K , T ) d T → 0 : ∂ C ( S , t , K , T ) ∂ T ∝ ∂ 2 C ( S , t , K , T ) ∂ K 2 σ 2 ( K , T ) σ 2 ( K , T ) ∝ ∂ C ( S , t , K , T ) ∂ T ∂ 2 C ( S , t , K , T ) ∂ K 2 C(S,t,K,T+dT)-C(S,t,K,T)\propto p(S,t,K,T)\sigma^2(K,T)dT\\ C(S,t,K,T+dT)-C(S,t,K,T)\propto \frac{\partial^2C(S,t,K,T)}{\partial K^2}\sigma^2(K,T)dT\\ \frac{C(S,t,K,T+dT)-C(S,t,K,T)}{dT}\propto \frac{\partial^2C(S,t,K,T)}{\partial K^2}\sigma^2(K,T)\\dT\to0:\frac{\partial C(S,t,K,T)}{\partial T}\propto \frac{\partial^2C(S,t,K,T)}{\partial K^2}\sigma^2(K,T)\\ \sigma^2(K,T)\propto \frac{\frac{\partial C(S,t,K,T)}{\partial T}}{\frac{\partial^2C(S,t,K,T)}{\partial K^2}} C(S,t,K,T+dT)−C(S,t,K,T)∝p(S,t,K,T)σ2(K,T)dTC(S,t,K,T+dT)−C(S,t,K,T)∝∂K2∂2C(S,t,K,T)​σ2(K,T)dTdTC(S,t,K,T+dT)−C(S,t,K,T)​∝∂K2∂2C(S,t,K,T)​σ2(K,T)dT→0:∂T∂C(S,t,K,T)​∝∂K2∂2C(S,t,K,T)​σ2(K,T)σ2(K,T)∝∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​​
这是局部波动率模型一个非常有吸引力的特点,那就是局部波动率与两个最常用的期权交易策略的比值密切相关。如果市场上有足够的跨期期权组合和蝶式期权组合的信息,我们就可以确定局部波动率曲面

若无风险利率不等于0:
σ 2 ( K , T ) 2 = ∂ C ( S , t , K , T ) ∂ T + r K ∂ C ( S , t , K , T ) ∂ K K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 ∂ C ( S , t , K , T ) ∂ T + r K ∂ C ( S , t , K , T ) ∂ K − σ 2 ( K , T ) 2 K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 = 0 \frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C(S,t,K,T)}{\partial T}+rK\frac{\partial C(S,t,K,T)}{\partial K}}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}}\\ \frac{\partial C(S,t,K,T)}{\partial T}+rK\frac{\partial C(S,t,K,T)}{\partial K}-\frac{\sigma^2(K,T)}{2}K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}=0 2σ2(K,T)​=K2∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​+rK∂K∂C(S,t,K,T)​​∂T∂C(S,t,K,T)​+rK∂K∂C(S,t,K,T)​−2σ2(K,T)​K2∂K2∂2C(S,t,K,T)​=0
上式类似于BSM等式,但是更加严格,其局部波动率模型只对普通欧式看涨期权或者看跌期权成立,将在下节做出证明

Dupire公式的价值在于,它可以告诉你如何通过一个标准期权的市场价格得到某特定的局部波动率函数 σ ( K , T ) \sigma(K,T) σ(K,T)。知道了所有 K K K 和 T T T 的局部波动率函数 σ ( K , T ) \sigma(K,T) σ(K,T),就可以建立一个包含这些局部波动率函数的隐含(波动率)树形图,并以此来计算奇异期权的价值,并且对标准期权进行对冲

DUPIRE公式的二叉树推导式

σ 2 ( K , T ) 2 = ∂ C ( S , t , K , T ) ∂ T K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 \frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C(S,t,K,T)}{\partial T}}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}} 2σ2(K,T)​=K2∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​​

在推导过程中用到了二叉树模型的相关理论框架,并且假设利率和股息率等于0

树形图

首先建立一个Jarrow-Rudd树形图,每次上行下行的概率都等于 1 / 2 1/2 1/2。在这个树形图中,令每个阶段的时间间隔都等于 d T / 2 dT/2 dT/2,这会使得跨期期权的估值更加简便,假设有个节点满足对应条件 ( K , T ) (K,T) (K,T),如下图所示:

Dupire公式的分子项与一个跨期期权组合价值成比例关系,该跨期期权包括买入一个行权价和期限分别为 ( K , T + d T ) (K,T+dT) (K,T+dT) 的看涨期权以及卖出一个行权价和期限分别为 ( K , T ) (K,T) (K,T) 的看涨期权。用 p K , T p_{K,T} pK,T​ 表示从条件 ( S , t ) (S,t) (S,t) 移动到 ( K , T ) (K,T) (K,T) 的风险中性概率

下面考虑,在时点 T T T,树形图中的各个节点对 ( K , T + d T ) (K,T+dT) (K,T+dT) 看涨期权以及 ( K , T ) (K,T) (K,T) 看涨期权的价值分别贡献了多少

  1. 在 T T T 时刻,只要节点股票价格 S T ′ S_T' ST′​ 低于 K K K 水平,那么 ( K , T ) (K,T) (K,T) 看涨期权就是都处于虚值状态。在此基础上,到了 T + d T T+dT T+dT 时刻,节点股票价格变为 S T + d T ′ S'_{T+dT} ST+dT′​,于是 ( K , T + d T ) (K,T+dT) (K,T+dT) 看涨期权的价值等于0。因此,在 T T T 时刻所有股票价格低于 K K K 水平的节点,对跨期期权组合的价值都没有影响
  2. 在到期日 T T T 时刻,股票价格 S T ′ S_T' ST′​ 高于 K K K 水平的任意一个节点,对 ( K , T ) (K,T) (K,T) 看涨期权的价值贡献都是 ( S T ′ − K ) (S_T'-K) (ST′​−K)。在 T + d T T+dT T+dT 时刻,该节点会衍生出3个平值或者实值状态的节点,对应的股票价格为 S T + d T ′ S'_{T+dT} ST+dT′​。因为所有的这些节点都是平值或者实值状态,于是我们可以证明,在 T + d T T+dT T+dT 时, ( K , T + d T ) (K,T+dT) (K,T+dT) 看涨期权在这3个节点对应的价值 ( S T + d T ′ − K ) (S'_{T+dT}-K) (ST+dT′​−K) 的预期折现值就等于在 T T T 时, ( K , T ) (K,T) (K,T) 看涨期权在节点 S T ′ S_T' ST′​ 的价值 ( S T ′ − K ) (S_T'-K) (ST′​−K),节点 S T ′ S_T' ST′​ 也是后3个节点的起点。因此,在 T T T 时刻,所有股票价格低于 K K K 水平的节点,对于跨期期权组合的价值都没有影响
  3. 考虑节点 ( K , T ) (K,T) (K,T) 本身。该节点本身就代表一个行权价为 K K K,期限为 T T T 的平值看涨期权,因此对于 ( K , T ) (K,T) (K,T) 看涨期权的价值没有影响。两个时间间隔阶段之后,在 T + d T T+dT T+dT 时点,该节点衍生出3个节点,如上图所示。对于行权价为 K K K,期限为 ( T + d T ) (T+dT) (T+dT) 的看涨期权来说,较低的两个节点分别处于虚值状态和平值状态,价值等于0。在3个节点中最高的节点,也就是股价等于 K e σ 2 Δ T Ke^{\sigma\sqrt{2\Delta T}} Keσ2ΔT ​ 时,行权价为 K K K 的看涨期权的价值就为正数。到达该节点的概率等于 1 / 4 × p K , T 1/4\times p_{K,T} 1/4×pK,T​

从 ( K , T ) (K,T) (K,T) 节点出发,沿着树形图,后续节点中有的会对 ( K , T + d T ) (K,T+dT) (K,T+dT) 看涨期权的价值有影响,而这些节点恰恰对 ( K , T ) (K,T) (K,T) 看涨期权的价值没有影响。除此之外,其他所有节点对两个看涨期权价值的影响都是相等的。因此跨期期权的价值就完全取决于上图中的粗线段所描述的转换路径。我们可以计算这个值是多少。每个时间间隔阶段等于 d T / 2 dT/2 dT/2,于是上行参数就等于:
u = σ d T 2 u=\sigma\sqrt{\frac{dT}{2}} u=σ2dT​ ​
其中 σ \sigma σ 是局部波动率 σ ( K , T ) \sigma(K,T) σ(K,T) 的简写,从节点 ( K , T ) (K,T) (K,T) 出发,连续向上运动两次最高节点对应的股票价格等于:
S = K e 2 u = K e σ 2 d T S=Ke^{2u}=Ke^{\sigma\sqrt{2dT}} S=Ke2u=Keσ2dT ​
跨期期权组合在 ( T + d T ) (T+dT) (T+dT) 时刻的价值就等于:
d K ≡ K e σ 2 d T − K ≈ K σ 2 d T dK\equiv Ke^{\sigma\sqrt{2dT}}-K\approx K\sigma\sqrt{2dT} dK≡Keσ2dT ​−K≈Kσ2dT ​
在时点 t t t,假设无风险利率等于0,那么跨期期权组合的风险中性价值就等于该损益乘以实现该损益的风险中性概率,于是:
C ( S , t , K , T + d T ) − C ( S , t , K , T ) = 1 4 p K , T d K ∂ C ( S , t , K , T ) ∂ T ≈ C ( S , t , K , T + d T ) − C ( S , t , K , T ) d T = 1 4 p K , T d K d T C(S,t,K,T+dT)-C(S,t,K,T)=\frac{1}{4}p_{K,T}dK\\ \frac{\partial C(S,t,K,T)}{\partial T}\approx\frac{C(S,t,K,T+dT)-C(S,t,K,T)}{dT}=\frac{1}{4}p_{K,T}\frac{dK}{dT} C(S,t,K,T+dT)−C(S,t,K,T)=41​pK,T​dK∂T∂C(S,t,K,T)​≈dTC(S,t,K,T+dT)−C(S,t,K,T)​=41​pK,T​dTdK​

蝶式价差期权组合

Dupire公式的分母与价差无限小的蝶式期权组合价值成比例关系:
C ( S , t , K − d K , T ) − 2 C ( S , t , K , T ) + C ( S , t , K + d K , T ) C(S,t,K-dK,T)-2C(S,t,K,T)+C(S,t,K+dK,T) C(S,t,K−dK,T)−2C(S,t,K,T)+C(S,t,K+dK,T)

上图展示了这个蝶式价差期权组合的损益情况,在 T T T 时刻,如果股票价格恰好等于节点 K K K,那么组合的收益就等于 d K dK dK,在相邻节点,组合的损益等于0

之前曾经定义过离散概率 p K , T p_{K,T} pK,T​,该概率约定,在节点 ( K , T ) (K,T) (K,T),组合将支付1美元,因此该概率就等于 1 / d K 1/dK 1/dK 份价差无限小的蝶式期权组合的价值,于是,在极限条件 d K → 0 dK\to0 dK→0 时:
p K , T = C ( S , t , K − d K , T ) − 2 C ( S , t , K , T ) + C ( S , t , K + d K , T ) d K = C ( S , t , K + d K , T ) − C ( S , t , K , T ) d K − C ( S , t , K , T ) − C ( S , t , K − d K , T ) d K ≈ ∂ C ( S , t , K , T ) ∂ K − ∂ C ( S , t , K − d K , T ) ∂ K ≈ ∂ 2 C ( S , t , K , T ) ∂ K 2 d K p_{K,T}=\frac{C(S,t,K-dK,T)-2C(S,t,K,T)+C(S,t,K+dK,T)}{dK}\\=\frac{C(S,t,K+dK,T)-C(S,t,K,T)}{dK}-\frac{C(S,t,K,T)-C(S,t,K-dK,T)}{dK}\\\approx\frac{\partial C(S,t,K,T)}{\partial K}-\frac{\partial C(S,t,K-dK,T)}{\partial K}\\\approx\frac{\partial^2C(S,t,K,T)}{\partial K^2}dK pK,T​=dKC(S,t,K−dK,T)−2C(S,t,K,T)+C(S,t,K+dK,T)​=dKC(S,t,K+dK,T)−C(S,t,K,T)​−dKC(S,t,K,T)−C(S,t,K−dK,T)​≈∂K∂C(S,t,K,T)​−∂K∂C(S,t,K−dK,T)​≈∂K2∂2C(S,t,K,T)​dK
于是上一节最后的等式可以写成:
∂ C ( S , t , K , T ) ∂ T = 1 4 p K , T d K d T = 1 4 ∂ 2 C ( S , t , K , T ) ∂ K 2 d K 2 d T \frac{\partial C(S,t,K,T)}{\partial T}=\frac{1}{4}p_{K,T}\frac{dK}{dT}=\frac{1}{4}\frac{\partial^2C(S,t,K,T)}{\partial K^2}\frac{dK^2}{dT} ∂T∂C(S,t,K,T)​=41​pK,T​dTdK​=41​∂K2∂2C(S,t,K,T)​dTdK2​
又因为 d K = K σ 2 d T dK=K\sigma\sqrt{2dT} dK=Kσ2dT ​,于是有:
∂ C ( S , t , K , T ) ∂ T = 1 2 σ 2 K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 → σ 2 ( K , T ) 2 = ∂ C ( S , t , K , T ) ∂ T K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 \frac{\partial C(S,t,K,T)}{\partial T}=\frac{1}{2}\sigma^2K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}\\\to\frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C(S,t,K,T)}{\partial T}}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}} ∂T∂C(S,t,K,T)​=21​σ2K2∂K2∂2C(S,t,K,T)​→2σ2(K,T)​=K2∂K2∂2C(S,t,K,T)​∂T∂C(S,t,K,T)​​
这就是当无风险利率等于0时,局部波动率的Dupire公式

DUPIRE公式的严格证明

下面给出更为严格的证明

首先给出风险中性条件下股票价格的随机偏微分方程如下:
d S S = r d t + σ ( S , t , . . . ) d Z \frac{dS}{S}=rdt+\sigma(S,t,...)dZ SdS​=rdt+σ(S,t,...)dZ
其中 σ ( S , t , . . . ) \sigma(S,t,...) σ(S,t,...) 表示股票价格的瞬时波动率,除了 S , t S,t S,t 之外可能还受其他变量的影响

在 t t t 时刻,一个行权价为 K K K,期限为 T T T​ 的普通欧式看涨期权的价值为
C ( S , t , K , T ) = e − r ( T − t ) E [ ( S T − K ) + ] C(S,t,K,T)=e^{-r(T-t)}E[(S_T-K)_+] C(S,t,K,T)=e−r(T−t)E[(ST​−K)+​]

引入 Heaviside 方程 H ( x ) H(x) H(x) 以及 Dirac delta 方程 δ ( x ) \delta(x) δ(x):
H ( x ) = { 1 , x ≥ 0 0 , x < 0 H(x)=\begin{cases}1,x\geq0\\0,x<0 \end{cases} H(x)={1,x≥00,x<0​
Dirac delta 方程 δ ( x ) \delta(x) δ(x) 可看作是 Heaviside 函数的广义导数

于是看涨期权的价值可以表示为:
C ( S , t , K , T ) = e − r ( T − t ) E [ ( S T − K ) H ( S T − K ) ] C(S,t,K,T)=e^{-r(T-t)}E[(S_T-K)H(S_T-K)] C(S,t,K,T)=e−r(T−t)E[(ST​−K)H(ST​−K)]
对上式分别对 K K K 求一阶导和二阶导:
∂ C ( S , t , K , T ) ∂ K = − e − r ( T − t ) E [ H ( S T − K ) ] ∂ 2 C ( S , t , K , T ) ∂ K 2 = e − r ( T − t ) E [ δ ( S T − K ) ] \frac{\partial C(S,t,K,T)}{\partial K}=-e^{-r(T-t)}E[H(S_T-K)]\\ \frac{\partial^2C(S,t,K,T)}{\partial K^2}=e^{-r(T-t)}E[\delta(S_T-K)] ∂K∂C(S,t,K,T)​=−e−r(T−t)E[H(ST​−K)]∂K2∂2C(S,t,K,T)​=e−r(T−t)E[δ(ST​−K)]
在 K K K 不变而 T T T 增加的情况下,要找到 C ( S , t , K , T ) C(S,t,K,T) C(S,t,K,T) 对 T T T 的总导数,我们需要考虑 d T dT dT 变动对 C C C 的影响,以及随机变量 S T S_T ST​ 变动对 C C C 的影响。因此:
d T C ∣ K = E [ ∂ C ∂ T d T + ∂ C ∂ S T d S T + 1 2 ∂ 2 C ∂ S T 2 ( d S T ) 2 ] = E [ − r C d T + e − r ( T − t ) H ( S T − K ) d S T + 1 2 e − r ( T − t ) δ ( S T − K ) ( d S T ) 2 ] = e − r ( T − t ) E [ − r e r ( T − t ) C d T + H ( S T − K ) d S T + 1 2 δ ( S T − K ) ( d S T ) 2 ] = e − r ( T − t ) E [ − r e r ( T − t ) C d T + H ( S T − K ) d S T + 1 2 δ ( S T − K ) σ 2 ( S T , T , . . . ) S T 2 d T ] d_TC|_K=E[\frac{\partial C}{\partial T}dT+\frac{\partial C}{\partial S_T}dS_T+\frac{1}{2}\frac{\partial^2C}{\partial S_T^2}(dS_T)^2]\\=E[-rCdT+e^{-r(T-t)}H(S_T-K)dS_T+\frac{1}{2}e^{-r(T-t)}\delta(S_T-K)(dS_T)^2]\\=e^{-r(T-t)}E[-re^{r(T-t)}CdT+H(S_T-K)dS_T+\frac{1}{2}\delta(S_T-K)(dS_T)^2]\\=e^{-r(T-t)}E[-re^{r(T-t)}CdT+H(S_T-K)dS_T+\frac{1}{2}\delta(S_T-K)\sigma^2(S_T,T,...)S_T^2dT] dT​C∣K​=E[∂T∂C​dT+∂ST​∂C​dST​+21​∂ST2​∂2C​(dST​)2]=E[−rCdT+e−r(T−t)H(ST​−K)dST​+21​e−r(T−t)δ(ST​−K)(dST​)2]=e−r(T−t)E[−rer(T−t)CdT+H(ST​−K)dST​+21​δ(ST​−K)(dST​)2]=e−r(T−t)E[−rer(T−t)CdT+H(ST​−K)dST​+21​δ(ST​−K)σ2(ST​,T,...)ST2​dT]
在最后一项中用 S T S_T ST​ 代替 K K K,可以得到:
d T C ∣ K = e − r ( T − t ) E [ − r e r ( T − t ) C d T + H ( S T − K ) d S T + 1 2 δ ( S T − K ) σ 2 ( S T , T , . . . ) K 2 d T ] d_TC|_K=e^{-r(T-t)}E[-re^{r(T-t)}CdT+H(S_T-K)dS_T+\frac{1}{2}\delta(S_T-K)\sigma^2(S_T,T,...)K^2dT] dT​C∣K​=e−r(T−t)E[−rer(T−t)CdT+H(ST​−K)dST​+21​δ(ST​−K)σ2(ST​,T,...)K2dT]
将 C ( S , t , K , T ) = e − r ( T − t ) E [ ( S T − K ) H ( S T − K ) ] C(S,t,K,T)=e^{-r(T-t)}E[(S_T-K)H(S_T-K)] C(S,t,K,T)=e−r(T−t)E[(ST​−K)H(ST​−K)] 代入得:
d T C ∣ K = e − r ( T − t ) E [ − r ( S T − K ) H ( S T − K ) d T + H ( S T − K ) d S T + 1 2 δ ( S T − K ) σ 2 ( S T , T , . . . ) K 2 d T ] d_TC|_K=e^{-r(T-t)}E[-r(S_T-K)H(S_T-K)dT+H(S_T-K)dS_T\\+\frac{1}{2}\delta(S_T-K)\sigma^2(S_T,T,...)K^2dT] dT​C∣K​=e−r(T−t)E[−r(ST​−K)H(ST​−K)dT+H(ST​−K)dST​+21​δ(ST​−K)σ2(ST​,T,...)K2dT]
根据股票价格的随机偏微分方程,可以将 d S T dS_T dST​ 进行如下代替:
d T C ∣ K = e − r ( T − t ) E [ − r ( S T − K ) H ( S T − K ) d T + H ( S T − K ) ( S T r d T + σ ( S T , T , . . . ) d Z ) + 1 2 δ ( S T − K ) σ 2 ( S T , T , . . . ) K 2 d T ] d T C ∣ K = e − r ( T − t ) E [ r K H ( S T − K ) d T + 1 2 δ ( S T − K ) σ 2 ( S T , T , . . . ) K 2 d T ] d_TC|_K=e^{-r(T-t)}E[-r(S_T-K)H(S_T-K)dT+H(S_T-K)(S_TrdT+\sigma(S_T,T,...)dZ)\\+\frac{1}{2}\delta(S_T-K)\sigma^2(S_T,T,...)K^2dT]\\ d_TC|_K=e^{-r(T-t)}E[rKH(S_T-K)dT+\frac{1}{2}\delta(S_T-K)\sigma^2(S_T,T,...)K^2dT] dT​C∣K​=e−r(T−t)E[−r(ST​−K)H(ST​−K)dT+H(ST​−K)(ST​rdT+σ(ST​,T,...)dZ)+21​δ(ST​−K)σ2(ST​,T,...)K2dT]dT​C∣K​=e−r(T−t)E[rKH(ST​−K)dT+21​δ(ST​−K)σ2(ST​,T,...)K2dT]
我们之前有:
∂ C ( S , t , K , T ) ∂ K = − e − r ( T − t ) E [ H ( S T − K ) ] ∂ 2 C ( S , t , K , T ) ∂ K 2 = e − r ( T − t ) E [ δ ( S T − K ) ] \frac{\partial C(S,t,K,T)}{\partial K}=-e^{-r(T-t)}E[H(S_T-K)]\\ \frac{\partial^2C(S,t,K,T)}{\partial K^2}=e^{-r(T-t)}E[\delta(S_T-K)] ∂K∂C(S,t,K,T)​=−e−r(T−t)E[H(ST​−K)]∂K2∂2C(S,t,K,T)​=e−r(T−t)E[δ(ST​−K)]
将上述左式作替换代入可得:
d T C ∣ K = − r K ∂ C ( S , t , K , T ) ∂ K d T + 1 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 E [ σ 2 ( K , T , . . . ) ] K 2 d T ∂ C ∂ T ∣ K = − r K ∂ C ( S , t , K , T ) ∂ K + 1 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 E [ σ 2 ( K , T , . . . ) ] K 2 E [ σ 2 ( K , T , . . . ) ] 2 = ∂ C ∂ T ∣ K + r K ∂ C ( S , t , K , T ) ∂ K ∣ T K 2 ∂ 2 C ( S , t , K , T ) ∂ K 2 ∣ T d_TC|_K=-rK\frac{\partial C(S,t,K,T)}{\partial K}dT+\frac{1}{2}\frac{\partial^2C(S,t,K,T)}{\partial K^2}E[\sigma^2(K,T,...)]K^2dT\\ \frac{\partial C}{\partial T}|_K=-rK\frac{\partial C(S,t,K,T)}{\partial K}+\frac{1}{2}\frac{\partial^2C(S,t,K,T)}{\partial K^2}E[\sigma^2(K,T,...)]K^2\\ \frac{E[\sigma^2(K,T,...)]}{2}=\frac{\frac{\partial C}{\partial T}|_K+rK\frac{\partial C(S,t,K,T)}{\partial K}|_T}{K^2\frac{\partial^2C(S,t,K,T)}{\partial K^2}|_T} dT​C∣K​=−rK∂K∂C(S,t,K,T)​dT+21​∂K2∂2C(S,t,K,T)​E[σ2(K,T,...)]K2dT∂T∂C​∣K​=−rK∂K∂C(S,t,K,T)​+21​∂K2∂2C(S,t,K,T)​E[σ2(K,T,...)]K22E[σ2(K,T,...)]​=K2∂K2∂2C(S,t,K,T)​∣T​∂T∂C​∣K​+rK∂K∂C(S,t,K,T)​∣T​​
上式也是Dupire公式。现在可以用未来在 K K K 价格水平和 T T T 时刻瞬时方差的平均数来表示广义局部方差 σ 2 ( K , T ) \sigma^2(K,T) σ2(K,T),该表达式还可以包括影响方差值的所有其他变量,即:
σ 2 ( K , T ) = E [ σ 2 ( K , T , . . . ) ] \sigma^2(K,T)=E[\sigma^2(K,T,...)] σ2(K,T)=E[σ2(K,T,...)]

局部波动率和隐含波动率之间的严格关系式及相关应用

当无风险利率和股息率都等于0的时候有如下关系:
σ 2 ( K , T ) 2 = ∂ C ∂ T ∣ K K 2 ∂ 2 C ∂ K 2 ∣ T \frac{\sigma^2(K,T)}{2}=\frac{\frac{\partial C}{\partial T}|_K}{K^2\frac{\partial^2C}{\partial K^2}|_T} 2σ2(K,T)​=K2∂K2∂2C​∣T​∂T∂C​∣K​​
若用 B S M BSM BSM 隐含波动率 Σ \Sigma Σ 来对期权进行估值,则有:
C ( S , t , K , T ) = C B S M ( S , t , K , T , Σ ( S , t , K , T ) ) C(S,t,K,T)=C_{BSM}(S,t,K,T,\Sigma(S,t,K,T)) C(S,t,K,T)=CBSM​(S,t,K,T,Σ(S,t,K,T))
在微分过程中运用链式法则,并且代人 BSM 模型中的各希腊字母参数,我们可以得到:
σ 2 ( K , τ ) = 2 ∂ Σ ∂ τ + Σ τ K 2 [ ∂ 2 Σ ∂ K 2 − d 1 τ ( ∂ Σ ∂ K ) 2 + 1 Σ ( 1 K τ + d 1 ∂ Σ ∂ K ) 2 ] τ = T − t , d 1 = 1 Σ τ ln ⁡ ( S K ) + 1 2 Σ τ \sigma^2(K,\tau)=\frac{2\frac{\partial\Sigma}{\partial\tau}+\frac{\Sigma}{\tau}}{K^2[\frac{\partial^2\Sigma}{\partial K^2}-d_1\sqrt{\tau}(\frac{\partial\Sigma}{\partial K})^2+\frac{1}{\Sigma}(\frac{1}{K\sqrt{\tau}}+d_1\frac{\partial\Sigma}{\partial K})^2]}\\ \tau=T-t,d_1=\frac{1}{\Sigma\sqrt{\tau}}\ln(\frac{S}{K})+\frac{1}{2}\Sigma\sqrt{\tau} σ2(K,τ)=K2[∂K2∂2Σ​−d1​τ ​(∂K∂Σ​)2+Σ1​(Kτ ​1​+d1​∂K∂Σ​)2]2∂τ∂Σ​+τΣ​​τ=T−t,d1​=Στ ​1​ln(KS​)+21​Στ ​
上式的局部波动率曲面公式 σ ( K , T ) \sigma(K,T) σ(K,T) 是用BSM隐含波动率来表示的,而不是用期权价格

在不存在斜度的条件下,隐含方差等于局部方差在期权有效期内的平均值

如果 Σ \Sigma Σ 跟行权价 K K K 无关,那么 σ ( K , τ ) \sigma(K,\tau) σ(K,τ) 也跟行权价无关,上式可以简化为:
σ 2 ( τ ) = 2 ∂ Σ ∂ τ + Σ τ K 2 [ 1 Σ ( 1 K τ ) 2 ] = 2 τ Σ ∂ Σ ∂ τ + Σ 2 = ∂ ∂ τ ( τ Σ 2 ) \sigma^2(\tau)=\frac{2\frac{\partial\Sigma}{\partial\tau}+\frac{\Sigma}{\tau}}{K^2[\frac{1}{\Sigma}(\frac{1}{K\sqrt{\tau}})^2]}\\=2\tau\Sigma\frac{\partial\Sigma}{\partial\tau}+\Sigma^2\\=\frac{\partial}{\partial\tau}(\tau\Sigma^2) σ2(τ)=K2[Σ1​(Kτ ​1​)2]2∂τ∂Σ​+τΣ​​=2τΣ∂τ∂Σ​+Σ2=∂τ∂​(τΣ2)
对等式两侧进行积分可以得到:
∫ 0 τ σ 2 ( u ) d u = τ Σ 2 ( τ ) \int_0^\tau\sigma^2(u)du=\tau\Sigma^2(\tau) ∫0τ​σ2(u)du=τΣ2(τ)
从一个更广义的角度来看,当 Σ \Sigma Σ 跟行权价格无关的时候,远期方差的合计数就等于总的隐含方差

两倍定律回顾

假设局部波动率只受行权价格的影响,而不受到期日的影响:
Σ = Σ ( K ) ∂ Σ ∂ τ = 0 \Sigma=\Sigma(K)\\ \frac{\partial\Sigma}{\partial\tau}=0 Σ=Σ(K)∂τ∂Σ​=0
同时,假设斜度跟 K K K 之间存在弱线性相关关系,这样我们就只需要保留与 ∂ C / ∂ K \partial C/\partial K ∂C/∂K 成比例关系的变量,而更高阶的 ( ∂ C / ∂ K ) 2 , ∂ 2 C / ∂ K 2 (\partial C/\partial K)^2,\partial^2C/\partial K^2 (∂C/∂K)2,∂2C/∂K2 就可以忽略不计。于是:
σ 2 ( K , T ) = Σ τ K 2 [ 1 Σ ( 1 K τ + d 1 ∂ Σ ∂ K ) 2 ] = Σ 2 ( 1 + d 1 K τ ∂ Σ ∂ K ) 2 σ ( K , τ ) = Σ ( K ) 1 + d 1 K τ ∂ Σ ∂ K \sigma^2(K,T)=\frac{\frac{\Sigma}{\tau}}{K^2[\frac{1}{\Sigma}(\frac{1}{K\sqrt{\tau}}+d_1\frac{\partial\Sigma}{\partial K})^2]}=\frac{\Sigma^2}{(1+d_1K\sqrt{\tau}\frac{\partial\Sigma}{\partial K})^2}\\ \sigma(K,\tau)=\frac{\Sigma(K)}{1+d_1K\sqrt{\tau}\frac{\partial\Sigma}{\partial K}} σ2(K,T)=K2[Σ1​(Kτ ​1​+d1​∂K∂Σ​)2]τΣ​​=(1+d1​Kτ ​∂K∂Σ​)2Σ2​σ(K,τ)=1+d1​Kτ ​∂K∂Σ​Σ(K)​
当期权接近于平值状态,也就是当 K = S + Δ K K=S+\Delta K K=S+ΔK 时,我们就有:
d 1 ≈ ln ⁡ ( S K ) Σ τ ≈ − Δ K S ( Σ τ ) ≈ − Δ K K ( Σ τ ) d_1\approx\frac{\ln(\frac{S}{K})}{\Sigma\sqrt{\tau}}\approx-\frac{\Delta K}{S(\Sigma\sqrt{\tau})}\approx-\frac{\Delta K}{K(\Sigma\sqrt{\tau})} d1​≈Στ ​ln(KS​)​≈−S(Στ ​)ΔK​≈−K(Στ ​)ΔK​
于是对于领头阶就有:
σ ( K ) ≈ Σ ( K ) 1 − Δ K Σ ∂ Σ ∂ K ≈ Σ ( K ) ( 1 + Δ K Σ ∂ Σ ∂ K ) ≈ Σ ( K ) + Δ K ∂ Σ ∂ K \sigma(K)\approx\frac{\Sigma(K)}{1-\frac{\Delta K}{\Sigma}\frac{\partial\Sigma}{\partial K}}\approx\Sigma(K)(1+\frac{\Delta K}{\Sigma}\frac{\partial\Sigma}{\partial K})\approx\Sigma(K)+\Delta K\frac{\partial\Sigma}{\partial K} σ(K)≈1−ΣΔK​∂K∂Σ​Σ(K)​≈Σ(K)(1+ΣΔK​∂K∂Σ​)≈Σ(K)+ΔK∂K∂Σ​
其中 K = S + Δ K K=S+\Delta K K=S+ΔK, Δ K \Delta K ΔK 表示股价与平值期权行权价的差异。因此:
σ ( S + Δ K ) ≈ Σ ( S + Δ K ) + Δ K ∂ Σ ( S + Δ K ) ∂ K \sigma(S+\Delta K)\approx\Sigma(S+\Delta K)+\Delta K\frac{\partial\Sigma(S+\Delta K)}{\partial K} σ(S+ΔK)≈Σ(S+ΔK)+ΔK∂K∂Σ(S+ΔK)​
在 S S S 进行泰勒展开可得:
σ ( S ) + ∂ σ ( S ) ∂ S Δ K ≈ Σ ( S ) + 2 ∂ Σ ( S ) ∂ S Δ K \sigma(S)+\frac{\partial\sigma(S)}{\partial S}\Delta K\approx\Sigma(S)+2\frac{\partial\Sigma(S)}{\partial S}\Delta K σ(S)+∂S∂σ(S)​ΔK≈Σ(S)+2∂S∂Σ(S)​ΔK
于是有 σ ( S ) ≈ Σ ( S ) \sigma(S)\approx\Sigma(S) σ(S)≈Σ(S),且
∂ ∂ S σ ( S ) ≈ 2 ( ∂ Σ ∂ K ) ∣ K = S \frac{\partial}{\partial S}\sigma(S)\approx2(\frac{\partial\Sigma}{\partial K})|_{K=S} ∂S∂​σ(S)≈2(∂K∂Σ​)∣K=S​
即:局部波动率随股票价格的变动速度是隐含波动率随行权价格的变动速度的两倍

当期限较短时,隐含波动率等于当前股价和行权价之间的局部波动率的调和平均数

对于一个正数序列, x 1 , x 2 , . . . , x n x_1,x_2,...,x_n x1​,x2​,...,xn​,其调和平均数 μ H \mu_H μH​ 的定义是:
μ H = n 1 x 1 + 1 x 2 + . . . + 1 x n \mu_H=\frac{n}{\frac{1}{x_1}+\frac{1}{x_2}+...+\frac{1}{x_n}} μH​=x1​1​+x2​1​+...+xn​1​n​
假设有一个连续随机变量,其正态密度函数为 f ( x ) f(x) f(x),那么其调和平均数就等于:
μ H = 1 ∫ 0 ∞ 1 x f ( x ) d x \mu_H=\frac{1}{\int_0^\infin\frac{1}{x}f(x)dx} μH​=∫0∞​x1​f(x)dx1​
考虑之前推导过的下式:
σ 2 ( K , τ ) = 2 ∂ Σ ∂ τ + Σ τ K 2 [ ∂ 2 Σ ∂ K 2 − d 1 τ ( ∂ Σ ∂ K ) 2 + 1 Σ ( 1 K τ + d 1 ∂ Σ ∂ K ) 2 ] \sigma^2(K,\tau)=\frac{2\frac{\partial\Sigma}{\partial\tau}+\frac{\Sigma}{\tau}}{K^2[\frac{\partial^2\Sigma}{\partial K^2}-d_1\sqrt{\tau}(\frac{\partial\Sigma}{\partial K})^2+\frac{1}{\Sigma}(\frac{1}{K\sqrt{\tau}}+d_1\frac{\partial\Sigma}{\partial K})^2]} σ2(K,τ)=K2[∂K2∂2Σ​−d1​τ ​(∂K∂Σ​)2+Σ1​(Kτ ​1​+d1​∂K∂Σ​)2]2∂τ∂Σ​+τΣ​​
假设将等式右侧分子分母同时乘以 τ \tau τ,可得:
σ 2 ( K , τ ) = 2 τ ∂ Σ ∂ τ + Σ K 2 [ τ ∂ 2 Σ ∂ K 2 − d 1 τ τ ( ∂ Σ ∂ K ) 2 + 1 Σ ( 1 K + d 1 τ ∂ Σ ∂ K ) 2 ] \sigma^2(K,\tau)=\frac{2\tau\frac{\partial\Sigma}{\partial\tau}+\Sigma}{K^2[\tau\frac{\partial^2\Sigma}{\partial K^2}-d_1\tau\sqrt{\tau}(\frac{\partial\Sigma}{\partial K})^2+\frac{1}{\Sigma}(\frac{1}{K}+d_1\sqrt{\tau}\frac{\partial\Sigma}{\partial K})^2]} σ2(K,τ)=K2[τ∂K2∂2Σ​−d1​ττ ​(∂K∂Σ​)2+Σ1​(K1​+d1​τ ​∂K∂Σ​)2]2τ∂τ∂Σ​+Σ​
当 τ → 0 \tau\to0 τ→0 时:
lim ⁡ τ → 0 σ 2 ( K , τ ) = Σ K 2 [ 1 Σ ( 1 K + d 1 τ ∂ Σ ∂ K ) 2 ] = Σ 2 ( 1 + d 1 τ K ∂ Σ ∂ K ) 2 \lim_{\tau\to0}\sigma^2(K,\tau)=\frac{\Sigma}{K^2[\frac{1}{\Sigma}(\frac{1}{K}+d_1\sqrt{\tau}\frac{\partial\Sigma}{\partial K})^2]}=\frac{\Sigma^2}{(1+d_1\sqrt{\tau}K\frac{\partial\Sigma}{\partial K})^2} τ→0lim​σ2(K,τ)=K2[Σ1​(K1​+d1​τ ​∂K∂Σ​)2]Σ​=(1+d1​τ ​K∂K∂Σ​)2Σ2​
其中:
d 1 τ = 1 Σ ln ⁡ ( S K ) + 1 2 Σ τ lim ⁡ τ → 0 d 1 τ = 1 Σ ln ⁡ ( S K ) d_1\sqrt{\tau}=\frac{1}{\Sigma}\ln(\frac{S}{K})+\frac{1}{2}\Sigma\tau\\ \lim_{\tau\to0}d_1\sqrt{\tau}=\frac{1}{\Sigma}\ln(\frac{S}{K}) d1​τ ​=Σ1​ln(KS​)+21​Σττ→0lim​d1​τ ​=Σ1​ln(KS​)
代入得:
lim ⁡ τ → 0 σ 2 ( K , τ ) = Σ 2 ( 1 + 1 Σ ln ⁡ ( S K ) K ∂ Σ ∂ K ) 2 σ ( K ) = lim ⁡ τ → 0 σ ( K , τ ) = Σ 1 + K Σ ln ⁡ ( S K ) d Σ d K \lim_{\tau\to0}\sigma^2(K,\tau)=\frac{\Sigma^2}{(1+\frac{1}{\Sigma}\ln(\frac{S}{K})K\frac{\partial\Sigma}{\partial K})^2}\\ \sigma(K)=\lim_{\tau\to0}\sigma(K,\tau)=\frac{\Sigma}{1+\frac{K}{\Sigma}\ln(\frac{S}{K})\frac{d\Sigma}{dK}} τ→0lim​σ2(K,τ)=(1+Σ1​ln(KS​)K∂K∂Σ​)2Σ2​σ(K)=τ→0lim​σ(K,τ)=1+ΣK​ln(KS​)dKdΣ​Σ​
令 x = ln ⁡ ( K / S ) x=\ln(K/S) x=ln(K/S),于是有:
d Σ d x = d Σ / ( S K × 1 S d K ) = K d Σ d K \frac{d\Sigma}{dx}=d\Sigma/(\frac{S}{K}\times\frac{1}{S}dK)=K\frac{d\Sigma}{dK} dxdΣ​=dΣ/(KS​×S1​dK)=KdKdΣ​
于是, σ ( K ) \sigma(K) σ(K) 可以被改写为:
σ ( K ) = Σ 1 + K Σ ln ⁡ ( S K ) d Σ d K = Σ 1 − x Σ d Σ d x \sigma(K)=\frac{\Sigma}{1+\frac{K}{\Sigma}\ln(\frac{S}{K})\frac{d\Sigma}{dK}}=\frac{\Sigma}{1-\frac{x}{\Sigma}\frac{d\Sigma}{dx}} σ(K)=1+ΣK​ln(KS​)dKdΣ​Σ​=1−Σx​dxdΣ​Σ​
令 V = 1 / Σ V=1/\Sigma V=1/Σ,于是:
d Σ d x = − 1 V 2 d V d x σ ( K ) = Σ 1 − x Σ d Σ d x = 1 / V 1 + x V 1 V 2 d V d x = 1 V [ 1 + x V d V d x ] V + x d V d x = 1 σ ( K ) d d x ( x V ) = 1 σ ( K ) \frac{d\Sigma}{dx}=-\frac{1}{V^2}\frac{dV}{dx}\\ \sigma(K)=\frac{\Sigma}{1-\frac{x}{\Sigma}\frac{d\Sigma}{dx}}=\frac{1/V}{1+xV\frac{1}{V^2}\frac{dV}{dx}}=\frac{1}{V[1+\frac{x}{V}\frac{dV}{dx}]}\\ V+x\frac{dV}{dx}=\frac{1}{\sigma(K)}\\ \frac{d}{dx}(xV)=\frac{1}{\sigma(K)} dxdΣ​=−V21​dxdV​σ(K)=1−Σx​dxdΣ​Σ​=1+xVV21​dxdV​1/V​=V[1+Vx​dxdV​]1​V+xdxdV​=σ(K)1​dxd​(xV)=σ(K)1​
此前,局部波动率 σ ( K ) \sigma(K) σ(K) 的表达式中,用到的参数是股票价格 K K K,现在我们将该参数转化为函数 ln ⁡ ( K / S ) \ln(K/S) ln(K/S)。因此,当我们提到 σ ( x ) \sigma(x) σ(x) 的时候,我们指的仍然是跟之前一样的波动率数值,只是我们将参数换成了函数 x = ln ⁡ ( K / S ) x=\ln(K/S) x=ln(K/S)。注意,在隐含波动率函数中,变量 K K K 表示行权价,而在局部波动率函数中, K K K 表示的是股票价格,这点非常重要。于是上式变成:
d d x ( x V ) = 1 σ ( x ) \frac{d}{dx}(xV)=\frac{1}{\sigma(x)} dxd​(xV)=σ(x)1​
对等式两侧同时积分,范围是从 x = 0 x=0 x=0 到 x = ln ⁡ ( K / S ) x=\ln(K/S) x=ln(K/S)​:
x V ∣ 0 ln ⁡ ( K / S ) = ∫ 0 ln ⁡ ( K / S ) 1 σ ( x ) d x ln ⁡ ( K S ) V ( S , K ) = ∫ 0 ln ⁡ ( K S ) 1 σ ( x ) d x ln ⁡ ( K S ) Σ ( K S ) = ∫ 0 ln ⁡ ( K S ) 1 σ ( x ) d x xV|_0^{\ln(K/S)}=\int_0^{\ln(K/S)}\frac{1}{\sigma(x)}dx\\ \ln(\frac{K}{S})V(S,K)=\int_0^{\ln(\frac{K}{S})}\frac{1}{\sigma(x)}dx\\ \frac{\ln(\frac{K}{S})}{\Sigma(\frac{K}{S})}=\int_0^{\ln(\frac{K}{S})}\frac{1}{\sigma(x)}dx xV∣0ln(K/S)​=∫0ln(K/S)​σ(x)1​dxln(SK​)V(S,K)=∫0ln(SK​)​σ(x)1​dxΣ(SK​)ln(SK​)​=∫0ln(SK​)​σ(x)1​dx
假设当前股票价格为 S S S,某标准期权的行权价为 K K K,随着到期日临近,期权的隐含波动率函数 ln ⁡ ( K / S ) \ln(K/S) ln(K/S) 就等于 S , K S,K S,K 之间所有股票价格 S ′ S' S′ 对应的局部波动率的调和平均数,其中局部波动率函数可以表示为 ln ⁡ ( S ′ / S ) \ln(S'/S) ln(S′/S),范围是0到 ln ⁡ ( K / S ) \ln(K/S) ln(K/S)

第15章 局部波动率的影响相关推荐

  1. 第16章 局部波动率模型——对冲比率及奇异期权估值

    这学期会时不时更新一下伊曼纽尔·德曼(Emanuel Derman) 教授与迈克尔B.米勒(Michael B. Miller)的<The Volatility Smile>这本书,本意是 ...

  2. 第17章 关于局部波动率的一些总结

    这学期会时不时更新一下伊曼纽尔·德曼(Emanuel Derman) 教授与迈克尔B.米勒(Michael B. Miller)的<The Volatility Smile>这本书,本意是 ...

  3. 一文弄懂元学习 (Meta Learing)(附代码实战)《繁凡的深度学习笔记》第 15 章 元学习详解 (上)万字中文综述

    <繁凡的深度学习笔记>第 15 章 元学习详解 (上)万字中文综述(DL笔记整理系列) 3043331995@qq.com https://fanfansann.blog.csdn.net ...

  4. 第 15 章 垃圾回收相关算法

    第 15 章 垃圾回收相关算法 1.标记阶段:引用计数器 1.1.标记阶段的目的 垃圾标记阶段:判断对象是否存活 在堆里存放着几乎所有的Java对象实例,在GC执行垃圾回收之前,首先需要区分出内存中哪 ...

  5. 第19章 随机波动率模型入门

    这学期会时不时更新一下伊曼纽尔·德曼(Emanuel Derman) 教授与迈克尔B.米勒(Michael B. Miller)的<The Volatility Smile>这本书,本意是 ...

  6. 第21章 随机波动率模型:无相关性时的微笑曲线

    这学期会时不时更新一下伊曼纽尔·德曼(Emanuel Derman) 教授与迈克尔B.米勒(Michael B. Miller)的<The Volatility Smile>这本书,本意是 ...

  7. 第15章 MiniVGGNet:更深的CNNs

    第15章 MiniVGGNet:更深的CNNs VGGNet(或简称为VGG),第一次在文献<Very Deep Learning Convolutional Neural Networks f ...

  8. 第15章习题解答(一)——《x86汇编语言:从实模式到保护模式》读书笔记40

    1. 第15章代码修改 先不说习题,说一说我对源码的修改.从运行结果来看,主要是增加了颜色支持.不过把我的代码与配书代码相比较的话,还是有很多不同的.这些修改是怎么来的,可以参考我之前的博文. 运行效 ...

  9. 《Real-Time Rendering 4th Edition》全文翻译 - 第15章 非真实感渲染(下)15.3 ~ 15.5

    连更两篇,冲鸭! 业余翻译,若有不周到之处,还请多多指教! 实时渲染(第四版)Real-Time Rendering (Fourth Edition) 第15章 非真实感渲染  Chapter 15  ...

最新文章

  1. Java项目:药品管理系统(java+swing+Gui+mysql)
  2. 根据经纬度批量计算多个点到多个点之间的距离
  3. 微型计算机的典型应用场景,单片机有哪些类型和应用场景?-MCU解决方案
  4. mysql 复制表结构到一个新表_MySQL复制表数据或表结构到新表中
  5. Linux是不是共享软件,linux – 是否有可能在应用程序之间共享Cuda上下文?
  6. LeetCode # Array # Easy # 217. Contains Duplicate
  7. 2021牛客暑期多校训练营4 E-Tree Xor(异或+思维+区间交 or Trie树)
  8. win2008 php mysql zend phpmyadmin_Windows2008 最新版Apache2.PHP5.MySQL6.PHPMyadmin.ZendOptimizer安装图解...
  9. mysql [ERR] 1273 - Unknown collation: ‘utf8mb4_0900_ai_ci‘
  10. leetcode390(2022 1.2)
  11. ubuntu系统安装python hello_ubuntu 下python安装及hello world
  12. 常用简体汉字unicode编码(3500个)
  13. 联众打码写滑动_如何能够对接联众打码平台
  14. 163邮箱邮件被退回,什么企业邮箱可以接收国外邮件?
  15. 「云渲染」渲染农场的架构原理、特点
  16. 用计算机唱歌 丑八怪乐谱,《成都》计算器谱子 抖音计算器按出的音乐乐谱大全...
  17. 王者荣耀会成为下一个微信吗?
  18. 树莓派4接7寸液晶屏竖屏设置
  19. 单片机并行口控制蜂鸣器播放音乐
  20. spire.office,Spire.Office 5.3.7组件的最新版本

热门文章

  1. LC-799.香槟塔
  2. Hadoop空闲时无法提交任务
  3. 服务器不装系统可以看配置么,怎么查看电脑服务器配置
  4. 选择示波器要考虑的十五大因素
  5. 【win10】磁盘保护清除
  6. 表达心情的的C语言程序,今天心情好来分享一段代码,嘎嘎嘎
  7. vue教程入门视频,vue入门视频教程
  8. 新零售风口之下,大数据及人工智能有哪些应用?
  9. A5s 150M便携式智能无线路由器
  10. vm眼影真的好用吗_便宜的眼影盘有什么坏处?