第 13 章 · 布朗运动与随机分析基础

13.3 强马尔可夫性与首达问题

Strong Markov Property and Hitting Times
布朗运动在确定时刻之后具有独立增量;若观察结束的时刻由路径决定,还需要强马尔可夫性。本节从通常增广滤过与停时的离散逼近出发,证明停时后的布朗增量独立于过去,再推导反射原理、最大值与单侧首达分布。区间退出概率和平均退出时间则通过有界截停计算,并逐步核对取极限所需的可积性。

13.3.1 通常增广滤过与停时的信息

设 \(B=(B_t)_{t\ge0}\) 是从 \(0\) 出发的一维标准布朗运动,采用上一节的连续路径版本。先把概率空间完备化,令 \(\mathcal N\) 为所有零测集的子集所组成的族,定义自然滤过的通常增广(usual augmentation):

\[\mathcal F_t^0=\sigma(B_s:0\le s\le t),\qquad \mathcal F_t=\bigcap_{u>t}\sigma(\mathcal F_u^0,\mathcal N).\]
(13.3-1)

它满足通常条件:\(\mathcal F_0\) 包含所有零测集,且 \(\mathcal F_t=\bigcap_{u>t}\mathcal F_u\),即滤过右连续。增广后,未来布朗增量仍独立于 \(\mathcal F_t\)。验证这一点时,可取 \(A\in\mathcal F_t\) 和 \(0<u_1<\cdots<u_m\):对 \(0<\varepsilon<u_1\),向量 \((B_{t+u_j}-B_{t+\varepsilon})_j\) 独立于 \(A\)。令 \(\varepsilon\downarrow0\),对有界连续测试函数使用路径连续性及控制收敛,再用单调类论证,便得到所需的独立性。补入零测集不改变这些积分等式。[1,§7.2]

定义 13.3-1 · 停时及停时 σ代数

随机时间 \(\tau:\Omega\to[0,\infty]\) 称为停时(stopping time),若对每个 \(t\ge0\),事件 \(\{\tau\le t\}\in\mathcal F_t\)。停时的信息由 \(\mathcal F_\tau=\{A\in\mathcal F:A\cap\{\tau\le t\}\in\mathcal F_t\text{,对所有 }t\ge0\}\) 表示。约定空集的首达时刻为无穷大。

例如对 \(a>0\),首达时刻 \(\tau_a=\inf\{t\ge0:B_t=a\}\) 是停时,因为 \(\{\tau_a\le t\}=\{\sup_{0\le s\le t}B_s\ge a\}\),而路径连续性允许用可数稠密时刻的上确界表示右侧。它并不要求在时刻 \(t\) 知道未来是否还会再次到达 \(a\)。

若 \(S\le T\) 都是停时,则 \(\mathcal F_S\subseteq\mathcal F_T\);若 \(T_n\downarrow T\),右连续性给出 \(\mathcal F_T=\bigcap_n\mathcal F_{T_n}\)。后一性质可由 \(\{T<t\}=\bigcup_n\{T_n<t\}\) 先处理严格不等式,再让 \(t+1/k\downarrow t\) 得到。配合连续路径的离散逼近,\(\tau\) 与 \(B_\tau\) 在 \(\{\tau<\infty\}\) 上都是 \(\mathcal F_\tau\) 可测的。

13.3.2 强马尔可夫性:从简单停时逼近

定理 13.3-1 · 布朗运动的强马尔可夫性

若 \(\tau\) 是几乎必然有限的停时,则 \(W_u=B_{\tau+u}-B_\tau\),\(u\ge0\),是标准布朗运动,并且整个过程 \(W\) 独立于 \(\mathcal F_\tau\)。特别地,对有界 Borel 函数 \(g:\mathbb R^m\to\mathbb R\) 及确定时刻 \(0\le u_1<\cdots<u_m\),

\[E[g(W_{u_1},\ldots,W_{u_m})\mid\mathcal F_\tau]=E[g(B_{u_1},\ldots,B_{u_m})]\quad\text{几乎必然}.\]
(13.3-2)

这是强马尔可夫性质(strong Markov property)。条件中的“停时”不可删除;有界随机时间也未必符合该条件。

证明

先取有限值简单停时 \(\tau\),其可能值为 \(t_k\)。记 \(V_r=(B_{r+u_j}-B_r)_{j=1}^m\),并令 \(c=E[g(V_0)]\)。对 \(A\in\mathcal F_\tau\),有 \(A\cap\{\tau=t_k\}\in\mathcal F_{t_k}\),所以确定时刻的增量独立性给出

\[E[\mathbf1_Ag(V_\tau)]=\sum_k E[\mathbf1_{A\cap\{\tau=t_k\}}g(V_{t_k})]=\sum_k P(A\cap\{\tau=t_k\})c=P(A)c.\]
(13.3-3)

当取值可数时同样成立,因为有界 \(g\) 使上面的级数绝对可求和。对一般的几乎必然有限停时,取严格从上逼近的二进制网格:

\[\tau_n=2^{-n}\bigl(\lfloor2^n\tau\rfloor+1\bigr),\qquad \tau<\tau_n\le\tau+2^{-n},\qquad \tau_n\downarrow\tau.\]
(13.3-4)

每个 \(\tau_n\) 都是停时:在网格点 \(k2^{-n}\) 上,\(\{\tau_n\le k2^{-n}\}=\{\tau<k2^{-n}\}\in\mathcal F_{k2^{-n}}\)。若先把 \(\tau\) 截在确定上界内,这些逼近都是有限值简单停时;一般情形可直接用已证明的可数取值情形。

由于 \(\mathcal F_\tau\subseteq\mathcal F_{\tau_n}\),式 (13.3-3) 对 \(\tau_n\) 及任意 \(A\in\mathcal F_\tau\) 成立。先令 \(g\) 有界连续,路径连续性给出 \(V_{\tau_n}\to V_\tau\),控制收敛使该积分等式传到极限。再用单调类定理扩展到有界 Borel 函数,得到式 (13.3-2)。过程 \(W\) 的路径连续,有限维分布与 \(B\) 相同;连续路径空间的 Borel σ代数由有理时刻坐标生成,因此有限维独立性也给出整个过程与 \(\mathcal F_\tau\) 的独立性。[2]

13.3.3 反射原理的完整推导

固定 \(a>0\)、\(t>0\),令 \(M_t=\max_{0\le s\le t}B_s\)。我们尚未证明 \(\tau_a\) 有限,因此先在有界停时 \(\sigma=\tau_a\wedge t\) 处反射:

\[\widetilde B_u=\begin{cases}B_u,&0\le u\le\sigma,\\2B_\sigma-B_u,&u>\sigma.\end{cases}\]
(13.3-5)
反射保持路径分布

由强马尔可夫性,\((B_{\sigma+v}-B_\sigma)_{v\ge0}\) 是独立于 \(\mathcal F_\sigma\) 的标准布朗运动。将这一整段增量变号,其分布不变,且仍独立于此前的信息。停时、停止前的路径与停止位置都包含在 \(\mathcal F_\sigma\) 中,故将变号后的增量接回原路径,所得 \(\widetilde B\) 与 \(B\) 有相同的路径分布。

取 \(b<a\)。对事件 \(\{M_t\ge a,\ B_t\le b\}\) 中的路径,反射前后首次触及 \(a\) 的时刻相同,而终点变为 \(2a-B_t\ge2a-b\)。反过来,任何终点至少为 \(2a-b>a\) 的连续路径必先触及 \(a\),再作同一反射就回到原事件。反射是自身的逆映射,且保持概率,因此

\[P(M_t\ge a,\ B_t\le b)=P(B_t\ge2a-b),\qquad a>0,\quad b<a.\]
(13.3-6)

令 \(b\uparrow a\),利用 \(P(B_t=a)=0\),左侧趋于 \(P(M_t\ge a,B_t<a)\),右侧趋于 \(P(B_t\ge a)\)。另一方面,\(\{B_t\ge a\}\subseteq\{M_t\ge a\}\)。把终点低于 \(a\) 与不低于 \(a\) 的两类路径相加,得到反射原理(reflection principle):

\[P(M_t\ge a)=P(\tau_a\le t)=2P(B_t\ge a)=2\left[1-\Phi\!\left(\frac a{\sqrt t}\right)\right].\]
(13.3-7)

其中 \(\Phi\) 是标准正态分布函数。整个推导只使用了有界停时 \(\sigma\) 的强马尔可夫性,没有预先使用首达时间的有限性。[1,例 7.4.2]

在首次触及水平a之后反射,低于a的终点b变为2a减b 反射示意:首次触及 a 之前不变,之后的增量变号 a00τat时间 b2a−b 原路径反射路径
图 13.3-1 折线仅作路径反射示意。对已经在 \(t\) 前触及 \(a\) 的路径,终点 \(b<a\) 与终点 \(2a-b>a\) 成对出现;路径连续性保证上方终点对应的路径必已触及该水平。

13.3.4 最大值、首达密度与无穷均值

反射原理给出 \(P(M_t\le a)=2\Phi(a/\sqrt t)-1\),\(a\ge0\),因而 \(M_t\overset d=|B_t|\),且 \(EM_t=\sqrt{2t/\pi}\)。这是一维分布相同,不能据此把两个随机变量逐路径视为相等。

把式 (13.3-7) 看成关于时间 \(t\) 的分布函数并求导,利用标准正态密度 \(\varphi\),得到首达时间密度(first-passage density):

\[\begin{aligned}f_{\tau_a}(t)&=\frac a{t^{3/2}}\varphi\!\left(\frac a{\sqrt t}\right)=\frac a{\sqrt{2\pi}\,t^{3/2}}\exp\!\left(-\frac{a^2}{2t}\right),\qquad t>0,\\P(\tau_a>t)&=2\Phi\!\left(\frac a{\sqrt t}\right)-1\sim a\sqrt{\frac2\pi}\,t^{-1/2}\quad(t\to\infty).\end{aligned}\]
(13.3-8)

尾概率趋于零,所以 \(P(\tau_a<\infty)=1\);但由 Tonelli 定理,\(E\tau_a=\int_0^\infty P(\tau_a>t)\,dt=\infty\)。最后的渐近式来自 \(\Phi(z)-1/2\sim z/\sqrt{2\pi}\),而 \(t^{-1/2}\) 在无穷远处不可积。首达几乎必然发生与平均等待无穷并不矛盾。缩放还给出 \(\tau_a\overset d=a^2\tau_1\);对负水平,利用 \(B\) 与 \(-B\) 同分布即可。[3]

例 13.3-1 · 曾越过阈值,终点却回到下方

标准布朗运动在时间 \(4\) 前触及水平 \(2\) 的概率为 \(P(\tau_2\le4)=2[1-\Phi(2/\sqrt4)]=2[1-\Phi(1)]\approx0.3173\)。若再要求终点不大于 \(0\),式 (13.3-6) 给出 \(P(M_4\ge2,B_4\le0)=P(B_4\ge4)=1-\Phi(2)\approx0.0228\)。后一个事件包含“先到达、再回落”的路径,单看终点低于阈值无法排除它。

因此时间 \(4\) 前的最大值均值为 \(EM_4=\sqrt{8/\pi}\),而第一次触及 \(2\) 的等待时间均值仍为无穷大。前者是固定时间窗内的空间最大值,后者是可延续任意久的时间变量。

13.3.5 区间退出:先截停,再证明可积

用布朗独立增量展开条件期望,可直接验证 \(B_t\) 与 \(B_t^2-t\) 都是鞅。连续时间可选停止不能只引用离散结论而略过极限;先建立本页所需的有界停时版本。

引理 13.3-1 · 有界截停的两个等式

若 \(\theta\le K<\infty\) 是停时,则

\[EB_\theta=0,\qquad EB_\theta^2=E\theta.\]
(13.3-9)
证明及极限的依据

用式 (13.3-4) 的网格从上逼近 \(\theta\),所得 \(\theta_m\) 受同一个整数 \(L>K+1\) 控制。在网格上对 \(B\) 与 \(B^2-t\) 使用离散有界停时定理,得到 \(EB_{\theta_m}=0\) 和 \(EB_{\theta_m}^2=E\theta_m\)。

Doob 的 \(L^2\) 最大不等式先用于 \([0,L]\) 的有限二进制网格,再由路径连续性与单调收敛传到稠密网格,给出 \(E[\sup_{0\le u\le L}|B_u|^2]\le4EB_L^2=4L\)。因此 \(B_{\theta_m}\) 和 \(B_{\theta_m}^2\) 分别被可积的上确界及其平方控制。路径连续性给出几乎必然收敛,控制收敛遂证明式 (13.3-9)。

定理 13.3-2 · 有界区间的退出分布与期望

令 \(X_t=x+B_t\),其中 \(\ell<x<r\) 且两个端点有限;取 \(T=\inf\{t\ge0:X_t\notin(\ell,r)\}\)。则 \(T\) 是可积停时,退出位置只可能是端点,并且

\[P(X_T=r)=\frac{x-\ell}{r-\ell},\qquad P(X_T=\ell)=\frac{r-x}{r-\ell},\qquad ET=(x-\ell)(r-x).\]
(13.3-10)
完整证明

连续路径使首次退出等于首次触及两端点之一,故 \(T\) 为停时。对每个整数 \(n\),只在有界停时 \(T\wedge n\) 上使用引理。由 \(EX_{T\wedge n}=x\),有

\[E(T\wedge n)=EX_{T\wedge n}^{\,2}-x^2\le C^2-x^2,\qquad C=\max(|\ell|,|r|).\]
(13.3-11)

不论此时是否已退出,路径连续性均保证 \(X_{T\wedge n}\in[\ell,r]\)。单调收敛于是先给出 \(ET\le C^2-x^2<\infty\),同时证明 \(T<\infty\) 几乎必然。这一步没有使用待求的期望公式。

现在可令 \(n\to\infty\):\(X_{T\wedge n}\to X_T\in\{\ell,r\}\),且位置和位置平方都有确定的可积上界。控制收敛给出 \(EX_T=x\) 与 \(ET=EX_T^2-x^2\)。记 \(p=P(X_T=r)\),则 \(x=\ell(1-p)+rp\),解得 \(p=(x-\ell)/(r-\ell)\)。最后代入 \(ET=\ell^2(1-p)+r^2p-x^2=(x-\ell)(r-x)\)。[1,定理 7.5.1–7.5.5]

标准布朗运动从区间零到二内部出发的上端退出概率和平均退出时间 区间 (0,2):上端退出概率与平均退出时间 012012 1010xx p(x)=x/2m(x)=x(2−x) (1,1/2)(1,1) 越靠近上端,从上端退出的概率越大 平均等待在中点达到最大值
图 13.3-2 标准布朗运动从 \(x\in(0,2)\) 出发时的精确函数图。端点处按 \(T=0\) 延拓,平均退出时间为零;上端退出概率随起点线性变化,平均退出时间为二次函数。
例 13.3-2 · 改变扩散尺度

过程 \(X_t=1+2B_t\) 首次离开 \((-1,4)\) 的时刻记为 \(T\)。按空间尺度转换,等价于 \(B\) 首次离开 \((-1,3/2)\)。所以从上端退出的概率为 \(1/(1+3/2)=2/5\),平均退出时间为 \(1\cdot(3/2)=3/2\)。

一般地,\(X_t=x+\sigma B_t\)、\(\sigma>0\) 时,上端退出概率仍是 \((x-\ell)/(r-\ell)\),而 \(ET=(x-\ell)(r-x)/\sigma^2\)。空间比例不受 \(\sigma\) 影响,时间尺度则除以方差增长率 \(\sigma^2\)。这里没有漂移项,加入漂移后不能继续套用这两个公式。

单侧首达不能直接代入无界可选停止

对 \(a>0\),虽然 \(B_{\tau_a}=a\) 几乎必然,但每个 \(n\) 都有 \(EB_{\tau_a\wedge n}=0\)。这些变量几乎必然收敛到 \(a\),却满足 \(E|B_{\tau_a\wedge n}-a|=a\),因为停止前位置始终不大于 \(a\)。它们甚至有 \(E|B_{\tau_a\wedge n}|\le2a\),仍然不一致可积;若一致可积,Vitali 定理会推出 \(L^1\) 收敛,与上述等式矛盾。这正是直接写 \(EB_{\tau_a}=EB_0\) 时遗漏的极限条件。

13.3.6 习题与详解

习题 13.3-1 · 相邻水平之间的首达时间

设 \(0<a<b\)。证明 \(\tau_b-\tau_a\) 独立于 \(\mathcal F_{\tau_a}\),且与 \(\tau_{b-a}\) 同分布。再求 \(P(\tau_1\le1,\ \tau_2-\tau_1\le1)\)。

解

由反射原理已知这些首达时刻均几乎必然有限。路径连续性保证 \(\tau_a<\tau_b\);令 \(W_u=B_{\tau_a+u}-a\),则 \(\tau_b-\tau_a=\inf\{u\ge0:W_u=b-a\}\)。强马尔可夫性使整个 \(W\) 独立于 \(\mathcal F_{\tau_a}\),并具有标准布朗运动的分布,对它取首达时间这一可测函数便得到结论。取 \(a=1,b=2\),两段等待独立且都与 \(\tau_1\) 同分布,因此所求概率为 \(\{2[1-\Phi(1)]\}^2\)。这并不表示两段等待的期望有限。

习题 13.3-2 · 给定终点后,路径是否始终低于阈值?

固定 \(a>0,t>0\)。利用反射原理求 \(P(M_t<a\mid B_t=y)\) 的一个条件概率版本,并说明对概率为零的事件 \(B_t=y\) 怎样理解该表达式。

解

记 \(p_t(y)=(2\pi t)^{-1/2}e^{-y^2/(2t)}\)。对 \(y<a\),式 (13.3-6) 给出事件 \(\{M_t\ge a\}\) 上终点的次概率密度 \(p_t(2a-y)\)。因此始终低于阈值时的终点密度是 \(p_t(y)-p_t(2a-y)\)。除以严格为正的终点密度,得到 \(1-p_t(2a-y)/p_t(y)=1-\exp\{-2a(a-y)/t\}\)。当 \(y\ge a\) 时可取零。这个函数是由联合密度与边缘密度之比确定的正则条件概率版本;条件期望只要求它对 \(B_t\) 的分布几乎处处成立,不是对零概率事件作普通概率除法。

习题 13.3-3 · 有界但使用未来的随机时间

令 \(S=1\) 当 \(B_2\ge0\),否则 \(S=2\)。证明 \(S\) 不是停时,并计算 \(EB_S\),说明“时间有界”本身不足以保证可选停止。

解

\(\{S\le1\}=\{B_2\ge0\}\),而 \(P(B_2\ge0\mid\mathcal F_1)=\Phi(B_1)\) 严格介于 \(0\) 与 \(1\) 之间,不能等于该事件的指示量,故事件不属于 \(\mathcal F_1\)。由联合高斯分布,\(E[B_1\mid B_2]=B_2/2\)。又 \(B_2\sim N(0,2)\),所以 \(E[B_2\mathbf1_{\{B_2\ge0\}}]=1/\sqrt\pi\)。于是 \(EB_S=\frac12E[B_2\mathbf1_{\{B_2\ge0\}}]+E[B_2\mathbf1_{\{B_2<0\}}]=-1/(2\sqrt\pi)\)。虽然 \(S\le2\),它不满足可选停止定理的停时前提。

习题 13.3-4 · 首达时间的拉普拉斯变换

对 \(a>0,q>0\),用指数鞅求 \(E[e^{-q\tau_a}]\),写明截停和取极限的依据。提示:取 \(\theta=\sqrt{2q}\)。

解

独立高斯增量给出 \(Z_t=\exp(\theta B_t-\theta^2t/2)\) 是鞅,且对每个有限 \(L\),有 \(EZ_L^2=e^{\theta^2L}<\infty\)。因此同引理 13.3-1 的网格论证,Doob 的 \(L^2\) 界保证在有界停时 \(\tau_a\wedge n\) 上 \(EZ_{\tau_a\wedge n}=1\)。由于 \(\theta>0\) 且 \(B_{\tau_a\wedge n}\le a\),有 \(0\le Z_{\tau_a\wedge n}\le e^{\theta a}\)。首达时间已经证明几乎必然有限,控制收敛于是给出 \(1=E[e^{\theta a-q\tau_a}]\),即 \(E[e^{-q\tau_a}]=e^{-a\sqrt{2q}}\)。该推导不要求 \(E\tau_a<\infty\),使用的是停止后指数变量的确定上界。

关键术语

  • 通常增广(usual augmentation):补入零测集并使滤过右连续。
  • 停时 σ代数(stopping-time σ-algebra):描述在随机停止时刻可确定的事件。
  • 强马尔可夫性(strong Markov property):在停时之后重新开始的布朗增量独立于该停时的信息。
  • 反射原理(reflection principle):在首达水平处反射后续路径以建立事件之间的等概率对应。
  • 首达时间(first-passage time):过程首次到达指定水平的时间。
  • 退出时间(exit time):过程首次离开指定区域的时间。
  • 有界截停(bounded stopping):用 \(T\wedge n\) 先限制观察时长,再检查极限。
  • 一致可积(uniform integrability):可把依概率收敛提升为 \(L^1\) 收敛的尾部控制条件。

参考文献

  1. Durrett, R. Probability: Theory and Examples, 5th ed., 2019。§7.2 讨论增广信息下的马尔可夫性质;定理 7.3.9 为强马尔可夫定理,例 7.4.2 为反射原理,§7.5 讨论布朗鞅、可选停止与退出时间。作者公开版本。
  2. Pitman, J.,记录者 Yan, D. The Strong Markov Property,UC Berkeley Stat 205B,2003,Lecture 16。停时 σ代数、从上离散逼近及强马尔可夫定理。大学官方讲义。
  3. Lalley, S. P. Brownian Motion,University of Chicago 课程讲义。反射原理、首达密度与 §3 “Wald Identities for Brownian Motion” 的有界截停计算。作者课程讲义。