CQF笔记M1L4随机分析和伊藤引理
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 H−1,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(RiRj∣i=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=1∑iRj,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=1∑in=1∑iE(RmRn)=n=1∑iE(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(S6∣R1,⋯,R5)=S5
备注:第五次抛出硬币后,第六次自身的期望为0,因此条件期望等于第五次的现金数量
马可夫属性 Markov property
随机变量 S i S_i Si的所有过去时间的条件分布,只依赖于上一时刻的值 S i − 1 S_{i-1} Si−1
- Markov property表示随机游走没有记忆
- 不要求随机变量 S i S_i Si的期望值与上一时刻 S i − 1 S_{i-1} Si−1的期望相同
- 更一般化一些,对于 1 ≤ j < j 1 \le j \lt j 1≤j<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(Si∣Sj,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=1∑i(Sj−Sj−1)2
由于抛硬币的 ( S j − S j − 1 ) 2 = 1 (S_j - S_{j-1})^2 = 1 (Sj−Sj−1)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=1∑6(Sj−Sj−1)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} n−1的速度下降, 而赌注以 n − 1 / 2 n^{-1/2} n−1/2的速度下降
当
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[n→∞limi=1∑nRi]=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]=n→∞limi=1∑nj=1∑nE(RiRj)=n→∞limi=1∑nE(Ri2)=n→∞limn×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} ti−1到 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(ti−1)∼N(0,ti−ti−1), d X ∼ N ( 0 , d t ) dX \sim N(0, dt) dX∼N(0,dt)
将
X
t
X_t
Xt看做随机游走的最终结果
X
t
∼
N
(
0
,
t
)
X_t \sim N(0, t)
Xt∼N(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)
dX∼N(0,dt), 即
d
X
=
ϕ
d
t
dX = \phi \sqrt{dt}
dX=ϕdt
随机变量的函数
F
(
X
)
=
X
2
F(X)=X^2
F(X)=X2
随机变量的微分
随机世界的两个变量:时间 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
这种形式的S可能为负值
带随机漂移的简单布朗运动
d S = μ S d t + σ S d X dS = \mu S dt + \sigma S dX dS=μSdt+σSdX
如果
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
当 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=1∑n(X(tj)−X(tj−1))2−t)2] where tj=jt/nE[j=1∑n(X(tj)−X(tj−1))4+2i=1∑nj<i∑n(X(ti)−X(ti−1))2(X(tj)−X(tj−1))2−2tj=1∑n(X(tj)−X(tj−1))2+t2]nn23t2+n(n−1)n2t2−2tnnt+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(tj−1)∼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(tj−1)2=t, 记做: ∫ 0 t ( d X ) 2 = t \int_0^t (dX)^2 = t ∫0t(dX)2=t
后续的证明中,提到相等equality,都指的是均方极限语境