CQF笔记M1L4随机分析和伊藤引理

CQF笔记M1L4随机分析和伊藤引理

      • Module 1 Building Blocks of Quant Finance
        • Lecture 4 Stochastic Calculus and Itˆo’s Lemma
          • 抛硬币实验
          • 马可夫属性 Markov property
          • 鞅属性 martingale property
          • 二次变分 Quadratic variation
          • 布朗运动 Brownian motion
          • 随机变量的函数
          • 泰勒级数和伊藤引理
          • 随机微分方程

Module 1 Building Blocks of Quant Finance

Lecture 4 Stochastic Calculus and Itˆo’s Lemma

由于金融市场潜在的随机本质(假设),随机分析在金融过程的数据建模中非常重要。

抛硬币实验

规则:抛出正面赢1元,抛出反面输1元,概率都是 1 2 \frac{1}{2} 21,任意抛硬币的行为相互独立

R i = { 1 , i f   H − 1 , i f   T \begin{aligned} R_i= \begin{cases} 1, if \ H\\ -1, if \ T\\ \end{cases} \end{aligned} Ri={1,if H1,if T

E ( R i ) = 1 × 1 2 + ( − 1 ) × 1 2 = 0 V ( R i ) = 1 2 × 1 2 + ( − 1 ) 2 × 1 2 = 1 E ( R i R j ∣ i ≠ j ) = 0 \begin{aligned} & E(R_i) = 1 \times \frac{1}{2}+(-1) \times \frac{1}{2} = 0 \\ & V(R_i) = 1^2 \times \frac{1}{2}+(-1)^2 \times \frac{1}{2} = 1 \\ & E(R_iR_j|i \neq j) = 0 \\ \end{aligned} E(Ri)=1×21+(1)×21=0V(Ri)=12×21+(1)2×21=1E(RiRji=j)=0

考察抛 i i i次硬币之后的现金总和
S i = ∑ j = 1 i R j , S 0 = 0 \begin{aligned} S_i = \sum_{j=1}^{i} R_j, S_0 = 0 \end{aligned} Si=j=1iRj,S0=0

E ( S i ) = 0 V ( S i ) = E ( S i 2 ) = ∑ m = 1 i ∑ n = 1 i E ( R m R n ) = ∑ n = 1 i E ( R n 2 ) = i \begin{aligned} & E(S_i) = 0 \\ & V(S_i) = E(S_i^2) = \sum_{m=1}^{i}\sum_{n=1}^{i}E(R_mR_n) = \sum_{n=1}^{i}E(R_n^2) = i \\ \end{aligned} E(Si)=0V(Si)=E(Si2)=m=1in=1iE(RmRn)=n=1iE(Rn2)=i

条件期望 E ( S 6 ∣ R 1 , ⋯   , R 5 ) = S 5 \begin{aligned} & E(S_6 | R_1, \cdots, R_5) = S_5 \\ \end{aligned} E(S6R1,,R5)=S5
备注:第五次抛出硬币后,第六次自身的期望为0,因此条件期望等于第五次的现金数量

马可夫属性 Markov property

随机变量 S i S_i Si的所有过去时间的条件分布,只依赖于上一时刻的值 S i − 1 S_{i-1} Si1

  • Markov property表示随机游走没有记忆
  • 不要求随机变量 S i S_i Si的期望值与上一时刻 S i − 1 S_{i-1} Si1的期望相同
  • 更一般化一些,对于 1 ≤ j < j 1 \le j \lt j 1j<j, 只有含有信息的最大的 j j j对应的 S j S_j Sj,对估计 S i S_i Si是有用的(有信息的)
鞅属性 martingale property

鞅属性的定义 E ( S i ∣ S j , j < i ) = S j \begin{aligned} & E(S_i | S_j, j < i) = S_j \\ \end{aligned} E(SiSj,j<i)=Sj

站在 j j j时刻看未来的 i i i时刻, i i i时刻的条件期望

二次变分 Quadratic variation

二次变分的定义 Q = ∑ j = 1 i ( S j − S j − 1 ) 2 \begin{aligned} & Q = \sum_{j=1}^{i} (S_j - S_{j-1})^2 \\ \end{aligned} Q=j=1i(SjSj1)2

由于抛硬币的 ( S j − S j − 1 ) 2 = 1 (S_j - S_{j-1})^2 = 1 (SjSj1)2=1, 所以 Q = i Q=i Q=i

布朗运动 Brownian motion

修改抛硬币实验的规则:

  • 抛6次,总时长为 t t t, 则单次实验的时长为 t / 6 t/6 t/6
  • 赌注从 1 1 1改为 t / 6 \sqrt{t/6} t/6

二次变分 Q = ∑ j = 1 6 ( S j − S j − 1 ) 2 = 6 ( t 6 ) 2 = t \begin{aligned} & Q = \sum_{j=1}^{6} (S_j - S_{j-1})^2 = 6 (\sqrt{\frac{t}{6}})^2 = t\\ \end{aligned} Q=j=16(SjSj1)2=6(6t )2=t

假设抛 n n n次,总时长仍然是 t t t, 赌注改为 t / n \sqrt{t/n} t/n , 二次变分 Q Q Q仍然等于 t t t

当逐渐增加 n n n, 步长以 n − 1 n^{-1} n1的速度下降, 而赌注以 n − 1 / 2 n^{-1/2} n1/2的速度下降

CQF笔记M1L4随机分析和伊藤引理_第1张图片

n = ∞ n=\infin n=
E [ S ( t ) ] = E [ lim ⁡ n → ∞ ∑ i = 1 n R i ] = 0 V [ S ( t ) ] = E [ S ( t ) 2 ] = t \begin{aligned} & E[S(t)] = E[\lim_{n \to \infty} \sum_{i=1}^{n}R_i] = 0 \\ & V[S(t)] = E[S(t)^2] = t \end{aligned} E[S(t)]=E[nlimi=1nRi]=0V[S(t)]=E[S(t)2]=t

E [ S ( t ) 2 ] = lim ⁡ n → ∞ ∑ i = 1 n ∑ j = 1 n E ( R i R j ) = lim ⁡ n → ∞ ∑ i = 1 n E ( R i 2 ) = lim ⁡ n → ∞ n × t n = t \begin{aligned} E[S(t)^2] = \lim_{n \to \infty} \sum_{i=1}^{n}\sum_{j=1}^{n}E(R_iR_j) = \lim_{n \to \infty} \sum_{i=1}^{n}E(R_i^2) = \lim_{n \to \infty} n \times \frac{t}{n} = t \end{aligned} E[S(t)2]=nlimi=1nj=1nE(RiRj)=nlimi=1nE(Ri2)=nlimn×nt=t

这种随机游走称为布朗运动,记为 X ( t ) X(t) X(t), 布朗运动有如下性质:

  • X ( 0 ) = x 0 = 0 X(0) = x_0 = 0 X(0)=x0=0
  • Finiteness: 赌注随时间的平方根增长
  • Continuity: 路径连续
  • Markov: 在给定的 τ \tau τ时刻, τ < t \tau \lt t τ<t, X ( t ) X(t) X(t)的条件分布仅依赖 X ( τ ) X(\tau) X(τ)
  • Martingale: 在给定的 τ \tau τ时刻, τ < t \tau \lt t τ<t, X ( t ) X(t) X(t)的条件期望是 X ( τ ) X(\tau) X(τ)
  • Normality: 在有限时间 t i − 1 t_{i-1} ti1 t i t_i ti, X ( t i ) − X ( t i − 1 ) ∼ N ( 0 , t i − t i − 1 ) X(t_i) - X(t_{i-1}) \sim N(0, t_i - t_{i-1}) X(ti)X(ti1)N(0,titi1), d X ∼ N ( 0 , d t ) dX \sim N(0, dt) dXN(0,dt)

X t X_t Xt看做随机游走的最终结果
X t ∼ N ( 0 , t ) X_t \sim N(0, t) XtN(0,t), 即 X t = ϕ t , ϕ   N ( 0 , 1 ) X_t = \phi \sqrt{t}, \phi ~ N(0, 1) Xt=ϕt ,ϕ N(0,1)

X X X的增量
d X ∼ N ( 0 , d t ) dX \sim N(0, dt) dXN(0,dt), 即 d X = ϕ d t dX = \phi \sqrt{dt} dX=ϕdt

随机变量的函数

F ( X ) = X 2 F(X)=X^2 F(X)=X2
CQF笔记M1L4随机分析和伊藤引理_第2张图片

随机变量的微分

随机世界的两个变量:时间 t t t和布朗运动 X X X

d X dX dX d t \sqrt{dt} dt 的同阶无穷小,因此远大于 d t dt dt

考虑梯度、斜率、微分、敏感性时,要非常小心,因为这些会使 d t dt dt趋向于0

因此在随机世界中,使用随机微分方程: d F = ⋯ d t + ⋯ d X dF=\cdots dt + \cdots dX dF=dt+dX

普通的微积分法则不适用于随机世界
假设 F ( X ) = X 2 F(X)=X^2 F(X)=X2, 则 d F ≠ 2 X d X dF \neq 2XdX dF=2XdX

泰勒级数和伊藤引理

简单泰勒展开,忽略高阶项
F ( X + d X ) = F ( X ) + d F d X d X + 1 2 d 2 F d X 2 d X 2 \begin{aligned} F(X+dX)=F(X)+\frac{dF}{dX}dX + \frac{1}{2}\frac{d^2F}{dX^2}dX^2 \end{aligned} F(X+dX)=F(X)+dXdFdX+21dX2d2FdX2

移项,得到
d F = d F d X d X + 1 2 d 2 F d X 2 d X 2 \begin{aligned} dF=\frac{dF}{dX}dX + \frac{1}{2}\frac{d^2F}{dX^2}dX^2 \end{aligned} dF=dXdFdX+21dX2d2FdX2

d X 2 dX^2 dX2在时间步长逐渐变小为 d t dt dt时,等于其均值 d t dt dt, 也就是说没有随机性。
这里似乎能说得通,但是没有严格证明,按照同样的逻辑 d X dX dX应该等于 d t \sqrt{dt} dt

RULE OF THUMB一般指拇指规则。拇指规则(RULE OF THUMB) ,中文又译为“大拇指规则 ”,又叫”经验法则“,是一种可用于许多情况的简单的,经验性的,探索性的但不是很准确的原则。

这里得到 d X 2 = d t dX^2=dt dX2=dt是Rule of thumb
似乎这就是伊藤引理的结论

对于 F = X 2 F=X^2 F=X2, 带入泰勒级数
d F d X = 2 X d 2 F d X 2 = 2 d F = d t + 2 X d X \begin{aligned} & \frac{dF}{dX} = 2X \\ & \frac{d^2F}{dX^2} = 2 \\ & dF = dt + 2XdX \\ \end{aligned} dXdF=2XdX2d2F=2dF=dt+2XdX
这是随机微分方程的一个例子

随机微分方程

随机微分方程用于建模随机量,例如股价。
这种随机量,例如股价,可以分解为确定性(可预测)和随机性两部分。

确定性部分用 d t dt dt建模,随机性部分用 d X dX dX建模: d S = f ( S , t )   d t + g ( S , t )   d X dS = f(S,t) \ dt + g(S,t) \ dX dS=f(S,t) dt+g(S,t) dX
其中 f ( S , t ) f(S,t) f(S,t)称为growth rate或者drift, g ( S , t ) g(S,t) g(S,t)与S的波动率有关。

带漂移的简单布朗运动

d S = μ d t + σ d X dS = \mu dt + \sigma dX dS=μdt+σdX

CQF笔记M1L4随机分析和伊藤引理_第3张图片
这种形式的S可能为负值

带随机漂移的简单布朗运动

d S = μ S d t + σ S d X dS = \mu S dt + \sigma S dX dS=μSdt+σSdX

CQF笔记M1L4随机分析和伊藤引理_第4张图片
如果 S S S的初值为正数,则 S S S永远不可能变为负数,因为 S S S越接近 0 0 0时, d S dS dS的漂移项越接近0

假设 σ = 0 \sigma = 0 σ=0, 变成普通的微分方程,求解得到 S t = S 0 e μ t S_t = S_0 e^{\mu t} St=S0eμt

均值复归的随机游走

d S = ( ν − μ S ) d t + σ d X dS = (\nu - \mu S)dt + \sigma dX dS=(νμS)dt+σdX

CQF笔记M1L4随机分析和伊藤引理_第5张图片

S S S变大到 d t dt dt的系数为负时, S S S的均值会变小

r r r替换 S S S, 就是短期利率的 Vasicek model

短期利率的Cox, Ingersoll & Ross model

d S = ( ν − μ S ) d t + σ S 1 / 2 d X dS = (\nu - \mu S)dt + \sigma S^{1/2} dX dS=(νμS)dt+σS1/2dX

当S接近0时,随机性会变小
这个微分方程在后面很重要

系数都是 t t t的函数

d W = g ( t ) d t + f ( t ) d W dW = g(t)dt + f(t)dW dW=g(t)dt+f(t)dW

这个SDE是下面随机积分的shorthand

W ( t ) = ∫ 0 t g ( τ ) d τ + ∫ 0 t f ( τ ) d X ( τ ) W(t) = \int_0^t g(\tau)d\tau + \int_0^t f(\tau)dX(\tau) W(t)=0tg(τ)dτ+0tf(τ)dX(τ)

mean square limit

用在随机积分中

考察随机量
E [ ( ∑ j = 1 n ( X ( t j ) − X ( t j − 1 ) ) 2 − t ) 2 ]   w h e r e   t j = j t / n = E [ ∑ j = 1 n ( X ( t j ) − X ( t j − 1 ) ) 4 + 2 ∑ i = 1 n ∑ j < i n ( X ( t i ) − X ( t i − 1 ) ) 2 ( X ( t j ) − X ( t j − 1 ) ) 2 − 2 t ∑ j = 1 n ( X ( t j ) − X ( t j − 1 ) ) 2 + t 2 ] = n 3 t 2 n 2 + n ( n − 1 ) t 2 n 2 − 2 t n t n + t 2 ∼ O ( 1 n ) \begin{aligned} &E[(\sum_{j=1}^{n}(X(t_j) - X(t_{j-1}))^2 - t)^2] \ where \ t_j = jt/n \\ = &E[\sum_{j=1}^{n}(X(t_j) - X(t_{j-1}))^4 + 2 \sum_{i=1}^{n} \sum_{j \lt i}^{n} (X(t_i) - X(t_{i-1}))^2 (X(t_j) - X(t_{j-1}))^2 - 2t\sum_{j=1}^{n}(X(t_j) - X(t_{j-1}))^2 + t^2] \\ = & n \frac{3t^2}{n^2} + n(n-1)\frac{t^2} {n^2}-2tn\frac{t}{n}+t^2 \\ \sim &O(\frac{1}{n}) \end{aligned} ==E[(j=1n(X(tj)X(tj1))2t)2] where tj=jt/nE[j=1n(X(tj)X(tj1))4+2i=1nj<in(X(ti)X(ti1))2(X(tj)X(tj1))22tj=1n(X(tj)X(tj1))2+t2]nn23t2+n(n1)n2t22tnnt+t2O(n1)

推导过程中有个隐含前提 X ( t j ) − X ( t j − 1 ) ∼ N ( 0 , t / n ) X(t_j) - X(t_{j-1}) \sim N(0, t/n) X(tj)X(tj1)N(0,t/n), 其二阶矩(方差)为 t / n t/n t/n, 四阶矩为 3 t 2 / n 2 3t^2/n^2 3t2/n2

n → ∞ n \to \infin n时前面计算的随机量等于0

在mean square limit 均方极限语境下,认为 ∑ j = 1 n ( X ( t j ) − X ( t j − 1 ) 2 = t \sum_{j=1}^{n}(X(t_j) - X(t_{j-1})^2 = t j=1n(X(tj)X(tj1)2=t, 记做: ∫ 0 t ( d X ) 2 = t \int_0^t (dX)^2 = t 0t(dX)2=t

后续的证明中,提到相等equality,都指的是均方极限语境

你可能感兴趣的:(CQF)