由于金融市场潜在的随机本质(假设),随机分析在金融过程的数据建模中非常重要。
规则:抛出正面赢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,因此条件期望等于第五次的现金数量
随机变量 S i S_i Si的所有过去时间的条件分布,只依赖于上一时刻的值 S i − 1 S_{i-1} Si−1
鞅属性的定义 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时刻的条件期望
二次变分的定义 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
修改抛硬币实验的规则:
二次变分 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 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
带随机漂移的简单布朗运动
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,都指的是均方极限语境