Phase 6 — Brownian Motion 基础

1. 从 SRW 到 Brownian Motion

Brownian Motion 不是凭空定义的。它是 Symmetric Random Walk 在"时间和空间同时缩放"下的极限。

设 $X_n$ 是 SRW:$X_n = \sum_{i=1}^n \varepsilon_i$,$\varepsilon_i = \pm 1$ w.p. $1/2$。

固定一个时间 $t > 0$,考虑 $\frac{1}{\sqrt{n}} X_{\lfloor nt \rfloor}$。这是"走了 $\lfloor nt \rfloor$ 步之后的位置,除以 $\sqrt{n}$"。

由 Central Limit Theorem:

$$\frac{1}{\sqrt{n}} X_{\lfloor nt \rfloor} = \frac{1}{\sqrt{n}} \sum_{i=1}^{\lfloor nt \rfloor} \varepsilon_i \xrightarrow{d} N(0, t)$$

($\lfloor nt \rfloor$ 个均值 0、方差 1 的 i.i.d. 求和,除以 $\sqrt{n}$,方差是 $\lfloor nt \rfloor / n \to t$)

更一般地,对多个时间点 $0 < t_1 < t_2 < \ldots < t_k$,增量

$$\frac{1}{\sqrt{n}}(X_{\lfloor nt_1 \rfloor}),\quad \frac{1}{\sqrt{n}}(X_{\lfloor nt_2 \rfloor} - X_{\lfloor nt_1 \rfloor}),\quad \ldots$$

是独立的(因为它们是不重叠区间上的 $\varepsilon_i$ 的和),每个都趋向于正态分布。

Brownian Motion 就是这个极限对象——“连续时间、连续空间的 SRW”。

2. Brownian Motion 的定义

$(B_t)_{t \geq 0}$ 是 Brownian Motion(也叫 Wiener Process),如果满足以下四个条件:

条件 1: $B_0 = 0$

条件 2: 对任意 $0 < t_1 < t_2 < \ldots < t_k$,增量

$$B_{t_1},\quad B_{t_2} - B_{t_1},\quad B_{t_3} - B_{t_2},\quad \ldots,\quad B_{t_k} - B_{t_{k-1}}$$

互相独立的。

条件 3: 对任意 $s < t$,$B_t - B_s \sim N(0, t - s)$。

增量的分布只取决于时间差 $t - s$,不取决于起点 $s$。方差等于时间长度。

条件 4: $t \mapsto B_t$ 是连续函数(a.s.)。

路径是连续的——没有跳跃。

从条件 2 和 3 可以推出: $(B_{t_1}, B_{t_2}, \ldots, B_{t_k})$ 是 multivariate normal(多元正态)。满足这个性质的过程叫做 Gaussian Process

3. BM 的基本性质

性质 A:BM 是 Markov Process

$$\text{给定 } B_t \text{ 的值,}(B_r)_{r > t} \text{ 的条件分布不依赖于 } (B_r)_{r \leq t} \text{ 中 } B_t \text{ 以外的信息}$$

为什么?因为 $B_r - B_t$($r > t$)独立于 $(B_s)_{s \leq t}$(由条件 2)。所以未来增量完全不受过去影响,只需要知道当前位置 $B_t$。

这和离散时间的 Markov property 完全类比。

性质 B:BM 有 Strong Markov Property

如果 $T$ 是 stopping time,那么给定 $T < \infty$ 和 $B_T$ 的值,过程 $(B_{T+s} - B_T)_{s \geq 0}$ 是一个全新的 BM,独立于 $T$ 时刻之前的信息。

这和离散 MC 的 Strong Markov Property 完全类比。后面 Reflection Principle 的证明会用到这个。

性质 C:BM 处处连续但处处不可微

路径是连续的(条件 4),但:

$$\frac{B_{t+\Delta t} - B_t}{\Delta t} \sim N\left(0, \frac{1}{\Delta t}\right)$$

当 $\Delta t \to 0$ 时,这个"导数"的方差趋向无穷。所以 $B_t$ 在任何点都没有导数(a.s.)。

这就是为什么后面需要 Itô calculus——普通微积分依赖可微性,BM 没有这个性质。

4. BM 作为 Martingale

$(B_t)_{t \geq 0}$ 是 martingale。

验证:对 $s < t$,

$$E[B_t | \mathcal{F}_s] = E[B_t - B_s + B_s | \mathcal{F}_s]$$

$B_s$ 在 $\mathcal{F}_s$ 中已知,提出来:

$$= B_s + E[B_t - B_s | \mathcal{F}_s]$$

$B_t - B_s$ 独立于 $\mathcal{F}_s$(独立增量性质),所以条件期望等于无条件期望:

$$= B_s + E[B_t - B_s] = B_s + 0 = B_s \quad \checkmark$$

$(B_t^2 - t)_{t \geq 0}$ 也是 martingale。

这和离散情况下 $X_n^2 - n$ 是 martingale 完全对应。验证:

$$E[B_t^2 - t | \mathcal{F}_s] = E[(B_s + (B_t - B_s))^2 - t | \mathcal{F}_s]$$

展开平方:$(B_s + (B_t - B_s))^2 = B_s^2 + 2B_s(B_t - B_s) + (B_t - B_s)^2$

逐项取条件期望:

  • $E[B_s^2 | \mathcal{F}_s] = B_s^2$(已知)
  • $E[2B_s(B_t - B_s) | \mathcal{F}_s] = 2B_s \cdot E[B_t - B_s] = 2B_s \cdot 0 = 0$($B_s$ 提出来,增量独立且均值零)
  • $E[(B_t - B_s)^2 | \mathcal{F}_s] = E[(B_t - B_s)^2] = t - s$(增量独立,方差 $= t - s$)

合起来:

$$= B_s^2 + 0 + (t - s) - t = B_s^2 - s \quad \checkmark$$

5. Continuous-Time Stopping Time 和 OST

Stopping Time 的连续时间定义

离散时间:${T = n} \in \mathcal{F}_n$。

连续时间的类比:随机时间 $T$ 是 stopping time,如果对所有 $t \geq 0$:

$${T \leq t} \in \mathcal{F}_t$$

注意这里用的是 ${T \leq t}$ 而不是 ${T = t}$。在连续时间中,$P(T = t) = 0$ 对大多数 $t$ 成立(因为 $T$ 通常是连续分布的),所以用 ${T = t}$ 不方便。${T \leq t}$ 问的是"到时间 $t$ 为止,$T$ 是否已经发生了"——这和离散情况的精神是一样的。

离散 stopping time 的所有性质都自然延伸:min、max 保持 stopping time,hitting time 是 stopping time,等等。

Continuous-Time OST

定理: 设 $(X_t)_{t \geq 0}$ 是连续时间 martingale(连续路径),$T$ 是 stopping time,$P(T < \infty) = 1$。如果:

(i) $E[|X_T|] < \infty$

(ii) $\lim_{t \to \infty} E[|X_t| \cdot \mathbf{1}(T > t)] = 0$

那么 $E[X_T] = E[X_0]$。

形式和离散版本完全一样,只是把 $n$ 换成 $t$。同样地,如果过程在 $T$ 之前有界,或者过程是 UI,OST 也成立。

6. 应用:Continuous-Time Gambler’s Ruin

问题: $(B_t)_{t \geq 0}$ 是 BM。$a, b > 0$。$T = \inf{t \geq 0 : B_t = -a \text{ or } B_t = b}$。求 $P(B_T = b)$ 和 $E[T]$。

Part 1:求 $P(B_T = b)$

用 $B_t$ 这个 martingale。

在 $T$ 之前,$-a \leq B_t \leq b$(还没 exit),所以 $|B_t| \leq \max(a, b)$,有界。OST 适用。

$$E[B_T] = E[B_0] = 0$$

$B_T$ 只能取 $b$ 或 $-a$:

$$b \cdot P(B_T = b) + (-a) \cdot P(B_T = -a) = 0$$

又 $P(B_T = b) + P(B_T = -a) = 1$,设 $p = P(B_T = b)$:

$$bp - a(1 - p) = 0$$

$$bp - a + ap = 0$$

$$p(a + b) = a$$

$$p = \frac{a}{a + b}$$

Part 2:求 $E[T]$

用 $M_t = B_t^2 - t$ 这个 martingale。

需要验证 OST 条件。这比 Part 1 更微妙,因为 $M_t$ 在 $T$ 之前不是有界的($B_t^2$ 有界,但 $-t$ 可以任意负)。

课上的处理方式:先证 $E[T] < \infty$。

怎么证?考虑一个固定的时间 $T_0$。从 $(-a, b)$ 区间的任意起点出发,BM 在时间 $T_0$ 内到达 $-a$ 或 $b$ 的概率有一个正的下界 $\eta > 0$(这是因为 $B_{T_0} \sim N(B_0, T_0)$,正态分布的尾概率始终为正)。

所以每过 $T_0$ 时间,有至少 $\eta$ 的概率停下来。这意味着:

$$P(T > kT_0) \leq (1 - \eta)^k$$

所以 $T$ 有指数衰减的尾概率,$E[T] < \infty$。

有了 $E[T] < \infty$,可以证明 $E[|M_t| \cdot \mathbf{1}(T > t)] \to 0$(因为 $|M_t| \leq \max(a,b)^2 + t$,而 $P(T > t)$ 指数衰减),OST 条件满足。

应用 OST:

$$E[M_T] = E[M_0] = B_0^2 - 0 = 0$$$$E[B_T^2 - T] = 0$$$$E[B_T^2] = E[T]$$

计算 $E[B_T^2]$:

$$E[B_T^2] = b^2 \cdot P(B_T = b) + (-a)^2 \cdot P(B_T = -a)$$$$= b^2 \cdot \frac{a}{a+b} + a^2 \cdot \frac{b}{a+b}$$$$= \frac{ab^2 + a^2 b}{a+b} = \frac{ab(a+b)}{a+b} = ab$$

所以 $E[T] = ab$。

注意这里用了两个不同的 martingale 解同一个问题的两个不同方面: $B_t$ 求 hitting probability,$B_t^2 - t$ 求 expected hitting time。这和离散 Gambler’s Ruin 中用 $X_n$ 和 $X_n^2 - n$ 的思路完全一样。

7. Reflection Principle

这是 Lecture 8 中一个重要的新技巧,离散时间没有直接对应物。

问题: $T_a = \inf{t > 0 : B_t \geq a}$(首次到达 $a$ 的时间),$a > 0$。求 $P(T_a \leq t)$。

关键观察:

$$P(T_a \leq t) = P\left(\max_{0 \leq s \leq t} B_s \geq a\right)$$

这是因为"在时间 $t$ 之前到达过 $a$“等价于"路径的最大值 $\geq a$"。

推导: 把事件 ${B_t \geq a}$ 按 $T_a$ 是否 $\leq t$ 来分:

$$P(B_t \geq a) = P(T_a \leq t) \cdot P(B_t \geq a ,|, T_a \leq t)$$

(如果 $T_a > t$,即在时间 $t$ 之前从未到过 $a$,那 $B_t \geq a$ 不可能发生。因为 BM 路径连续,不到 $a$ 就不可能超过 $a$。所以 $P(B_t \geq a, T_a > t) = 0$。)

现在算 $P(B_t \geq a ,|, T_a \leq t)$。

给定 $T_a \leq t$,BM 在时间 $T_a$ 时到达了 $a$。由 Strong Markov Property,$(B_{T_a + s} - B_{T_a})_{s \geq 0}$ 是一个全新的 BM,独立于 $T_a$ 之前的信息。

所以:

$$P(B_t \geq a ,|, T_a \leq t) = P(B_t - B_{T_a} \geq 0 ,|, T_a \leq t)$$

$B_t - B_{T_a}$ 是新 BM 在时间 $t - T_a$ 的值,分布是 $N(0, t - T_a)$。正态分布关于 0 对称,所以:

$$P(B_t - B_{T_a} \geq 0 ,|, T_a \leq t) = \frac{1}{2}$$

代回去:

$$P(B_t \geq a) = P(T_a \leq t) \cdot \frac{1}{2}$$

所以:

$$P(T_a \leq t) = 2 \cdot P(B_t \geq a) = 2\Phi\left(-\frac{a}{\sqrt{t}}\right)$$

(其中 $\Phi$ 是标准正态 CDF,$P(B_t \geq a) = P\left(\frac{B_t}{\sqrt{t}} \geq \frac{a}{\sqrt{t}}\right) = 1 - \Phi\left(\frac{a}{\sqrt{t}}\right) = \Phi\left(-\frac{a}{\sqrt{t}}\right)$)

Reflection Principle 的一个推论

$$P\left(\max_{0 \leq s \leq t} B_s \geq a\right) = P(T_a \leq t) = 2P(B_t \geq a) = P(|B_t| \geq a)$$

所以 $\max_{0 \leq s \leq t} B_s$ 和 $|B_t|$ 有相同的分布。

这个结论很漂亮:BM 路径在 $[0, t]$ 上的最大值的分布,完全由终点 $B_t$ 的绝对值决定。

8. Phase 6 小结

概念要点
BM 定义$B_0=0$、独立增量、$B_t - B_s \sim N(0, t-s)$、连续路径
BM 作为 MG$B_t$ 是 MG,$B_t^2 - t$ 也是 MG
Continuous-time OST形式和离散版本相同,条件相同
Gambler’s Ruin$B_t$ → hitting prob $= a/(a+b)$;$B_t^2 - t$ → $E[T] = ab$
Reflection Principle$P(T_a \leq t) = 2P(B_t \geq a)$,$\max B_s \sim$
处处不可微导数的方差 $\to \infty$ → 需要新的 calculus