9. Multi-Dimensional Brownian Motion
9. Multi-Dimensional Brownian Motion
本章导读 本章把布朗运动推广到 \(d\) 维。§9.1 定义(Def 9.1 由 \(d\) 个独立一维 BM 拼成 \(\mathbf B_t\);Def 9.2 等价四条件定义;Fact 9.1 密度;Def 9.3 带漂移多维 BM)。§9.2 Dirichlet 问题(Def 9.4/9.5 调和函数的两个定义——球面均值性质 = Laplace 方程 \(\Delta f=0\);Thm 9.1 最大值原理;Dirichlet 问题及 Prop 9.1 解的唯一性)。§9.3 常返 vs 暂留:9.3.1 退出概率(球壳上首达外球的概率 \(f\) 是 Dirichlet 问题的唯一解,球对称假设 \(\varphi(|\mathbf x|^2)\),解 ODE 得三维情形公式 (9.4));9.3.2 定义(Def 9.6 点常返、Def 9.7 邻域常返、Def 9.8 暂留);9.3.3 一维=点常返(Thm 9.2);9.3.4 二维=邻域常返但非点常返(Thm 9.3);9.3.5 三维及以上=暂留(Thm 9.4)。§9.4 Green 函数:9.4.1 转移密度 \(P_t(\mathbf x,\mathbf y)\) (9.8);9.4.2 热方程(Def 9.9、Prop 9.2);9.4.3 Green 函数 \(G(\mathbf x,\mathbf y)=\int_0^\infty P_t\,dt=C_d|\mathbf x-\mathbf y|^{2-d}\) (9.10);9.4.4 条件 Green 函数(嵌入停时,(9.11))。无图。
9. Multi-Dimensional Brownian Motion
Overview This chapter extends Brownian motion to \(d\) dimensions. §9.1 definitions (Def 9.1 \(\mathbf B_t\) assembled from \(d\) independent one-dimensional BMs; Def 9.2 the equivalent four-condition definition; Fact 9.1 density; Def 9.3 multi-dimensional BM with drift). §9.2 Dirichlet problem (Def 9.4/9.5 two definitions of a harmonic function — the spherical mean-value property $=$ Laplace's equation \(\Delta f=0\); Thm 9.1 maximum principle; the Dirichlet problem and Prop 9.1 uniqueness of its solution). §9.3 recurrence vs transience: 9.3.1 exit probability (the probability \(f\) of hitting the outer sphere first on a spherical shell is the unique Dirichlet solution; under the spherical-symmetry ansatz \(\varphi(|\mathbf x|^2)\), solving the ODE gives the dimension-dependent formula (9.4)); 9.3.2 definitions (Def 9.6 point recurrent, Def 9.7 neighborhood recurrent, Def 9.8 transient); 9.3.3 1-D is point recurrent (Thm 9.2); 9.3.4 2-D is neighborhood- but not point-recurrent (Thm 9.3); 9.3.5 3-D and higher are transient (Thm 9.4). §9.4 Green's function: 9.4.1 transition density \(P_t(\mathbf x,\mathbf y)\) (9.8); 9.4.2 heat equation (Def 9.9, Prop 9.2); 9.4.3 Green's function \(G(\mathbf x,\mathbf y)=\int_0^\infty P_t\,dt=C_d|\mathbf x-\mathbf y|^{2-d}\) (9.10); 9.4.4 conditional Green's function (embedding a stopping time, (9.11)). No figures.
9.1 定义 / Definitions
9.1 Definitions
定义 9.1、9.2(多维布朗运动)/ Definitions 9.1, 9.2
定义 9.1(多维布朗运动):设 \(\{B_{1,t}\},\{B_{2,t}\},\dots,\{B_{d,t}\}\) 是 \(d\) 个独立的标准一维布朗运动,则定义 \(\{\mathbf B_t\}\)(\(\forall t\),\(\mathbf B_t=(B_{1,t},B_{2,t},\dots,B_{d,t})\))为 \(d\) 维标准布朗运动。定义 9.2(等价定义):设 \(\{\mathbf B_t\}\),\(\mathbf B_t\) 是 \(\mathbb R^d\) 中 \(d\) 维向量,满足:(1) 从原点出发 \(\mathbf B_0=\mathbf 0\);(2) 增量独立:\(rDefinition 9.1 (Multi-dimensional Brownian motion): let \(\{B_{1,t}\},\{B_{2,t}\},\dots,\{B_{d,t}\}\) be \(d\) independent standard one-dimensional Brownian motions; then define \(\{\mathbf B_t\}\) (for all \(t\), \(\mathbf B_t=(B_{1,t},B_{2,t},\dots,B_{d,t})\)) to be a \(d\)-dimensional standard Brownian motion. Definition 9.2 (Alternative definition): suppose \(\{\mathbf B_t\}\) with \(\mathbf B_t\) a \(d\)-dimensional vector in \(\mathbb R^d\) satisfies: (1) starts from the origin \(\mathbf B_0=\mathbf 0\); (2) independent increments: for \(r
事实 9.1(密度)与定义 9.3(带漂移)/ Fact 9.1 (density) and Definition 9.3 (with drift) 事实 9.1:标准 \(d\) 维布朗运动的增量 \(\mathbf B_t-\mathbf B_s\) 的密度为(\(\mathbf x=(x_1,\dots,x_d)\in\mathbb R^d\))Fact 9.1: the density of the increment \(\mathbf B_t-\mathbf B_s\) of a standard \(d\)-dimensional Brownian motion is (for \(\mathbf x=(x_1,\dots,x_d)\in\mathbb R^d\))
$$f(\mathbf x)=\prod_{j=1}^d\frac1{\sqrt{2\pi(t-s)}}e^{-\frac{x_j^2}{2(t-s)}}=\frac1{[2\pi(t-s)]^{\frac d2}}e^{-\frac{|\mathbf x|^2}{2(t-s)}}.$$
定义 9.3(带漂移的多维布朗运动):设 \(\{\mathbf B_t\}\) 是标准 \(d\) 维布朗运动,若 \(\mathbf Y_t=\mathbf A\mathbf B_t+t\boldsymbol\mu\)(\(\mathbf A\) 是 \(k\times d\) 矩阵,\(\boldsymbol\mu\) 是 \(k\) 维向量),则 \(\{\mathbf Y_t\}\) 是 \(k\) 维带漂移布朗运动,\(\forall t\),\(\mathbf Y_t\sim\mathcal N(t\boldsymbol\mu,\mathbf A\mathbf A')\)。Definition 9.3 (Multi-dimensional Brownian motion with drift): let \(\{\mathbf B_t\}\) be a standard \(d\)-dimensional Brownian motion; if \(\mathbf Y_t=\mathbf A\mathbf B_t+t\boldsymbol\mu\) (\(\mathbf A\) a \(k\times d\) matrix, \(\boldsymbol\mu\) a \(k\)-dimensional vector), then \(\{\mathbf Y_t\}\) is a \(k\)-dimensional Brownian motion with drift, and for all \(t\), \(\mathbf Y_t\sim\mathcal N(t\boldsymbol\mu,\mathbf A\mathbf A')\).
9.2 Dirichlet 问题 / Dirichlet Problem
定义 9.4、9.5(调和函数)/ Definitions 9.4, 9.5 (Harmonic function) 定义 9.4(调和函数·球面均值性质):函数 \(f:D\to\mathbb R\) 是调和的,若它局部可积(脚注 9.1:\(f\) 局部可积指 \(\forall\mathbf x\in D\),\(\exists\varepsilon>0\) 使 \(\int_{|\mathbf y-\mathbf x|<\varepsilon}f(\mathbf y)\,d\mathbf y<\infty\))且满足球面均值性质(脚注 9.2:任意极小球面上所有点的均值等于其球心处的值):当 \(\mathbf x\in D\) 且 \(\varepsilon<\mathrm{distance}(\mathbf x,\partial D)\) 时,\(f(\mathbf x)=MV(f;\mathbf x,\varepsilon)=\int_{|\mathbf y-\mathbf x|<\varepsilon}f(\mathbf y)\,ds_\varepsilon(\mathbf y)\),其中 \(s_\varepsilon(\cdot)\) 是归一化的表面测度(脚注 9.3),\(\int_{|\mathbf y-\mathbf x|<\varepsilon}ds_\varepsilon(\mathbf y)=1\)。定义 9.5(调和函数·Laplace 方程):\(f:D\to\mathbb R\)(\(D\) 是 \(\mathbb R^d\) 的开子集)是调和的,若它满足 Laplace 方程,即 \(\forall(x_1,\dots,x_d)\in D\),\(\Delta f\equiv\dfrac{\partial^2 f}{(\partial x_1)^2}+\dfrac{\partial^2 f}{(\partial x_2)^2}+\dots+\dfrac{\partial^2 f}{(\partial x_d)^2}=0\)。Definition 9.4 (Harmonic function, mean-value form): a function \(f:D\to\mathbb R\) is harmonic if it is locally integrable (footnote 9.1: \(f\) is locally integrable if for all \(\mathbf x\in D\) there is \(\varepsilon>0\) with \(\int_{|\mathbf y-\mathbf x|<\varepsilon}f(\mathbf y)\,d\mathbf y<\infty\)) and satisfies the spherical mean-value property (footnote 9.2: the mean of all points on the surface of any extremely small ball equals the value at its center): when \(\mathbf x\in D\) and \(\varepsilon<\mathrm{distance}(\mathbf x,\partial D)\), \(f(\mathbf x)=MV(f;\mathbf x,\varepsilon)=\int_{|\mathbf y-\mathbf x|<\varepsilon}f(\mathbf y)\,ds_\varepsilon(\mathbf y)\), where \(s_\varepsilon(\cdot)\) is the normalized surface measure (footnote 9.3), \(\int_{|\mathbf y-\mathbf x|<\varepsilon}ds_\varepsilon(\mathbf y)=1\). Definition 9.5 (Harmonic function, Laplace form): \(f:D\to\mathbb R\) (\(D\) an open subset of \(\mathbb R^d\)) is harmonic if it satisfies Laplace's equation, i.e. for all \((x_1,\dots,x_d)\in D\), \(\Delta f\equiv\dfrac{\partial^2 f}{(\partial x_1)^2}+\dfrac{\partial^2 f}{(\partial x_2)^2}+\dots+\dfrac{\partial^2 f}{(\partial x_d)^2}=0\).
定理 9.1(最大值原理)/ Theorem 9.1 (Maximum principle on a bounded domain \(D\)) 若 \(f\) 是调和函数,则 \(f\) 没有局部极大值,除非 \(f\) 是常数。If \(f\) is a harmonic function, then \(f\) has no local maximum unless \(f\) is constant.
定理 9.1 证明 / Proof of Theorem 9.1 反证:若不然,存在 \(\mathbf x_0\in D\) 使得对某 \(\varepsilon>0\),\(f(\mathbf x_0)\ge f(\mathbf x)\) 对 \(\forall\mathbf x:|\mathbf x-\mathbf x_0|<\varepsilon\) 成立且至少在某处严格大于。则 \(f(\mathbf x_0)>\int_{|\mathbf x-\mathbf x_0|<\varepsilon}f(\mathbf x)\,ds_\varepsilon(\mathbf x)=MV(f;\mathbf x_0,\varepsilon)\),与球面均值性质矛盾。\(\blacksquare\)By contradiction: if not, there is \(\mathbf x_0\in D\) such that for some \(\varepsilon>0\), \(f(\mathbf x_0)\ge f(\mathbf x)\) for all \(\mathbf x:|\mathbf x-\mathbf x_0|<\varepsilon\) with strict inequality at least somewhere. Then \(f(\mathbf x_0)>\int_{|\mathbf x-\mathbf x_0|<\varepsilon}f(\mathbf x)\,ds_\varepsilon(\mathbf x)=MV(f;\mathbf x_0,\varepsilon)\), contradicting the mean-value property. \(\blacksquare\)
Dirichlet 问题与命题 9.1(唯一性)/ The Dirichlet problem and Proposition 9.1 (uniqueness) Dirichlet 问题:设 \(D\subset\mathbb R^d\) 是定义域(开且连通),\(F:\partial D\to\mathbb R\) 连续且有界。求一个有界函数 \(f:\bar D\to\mathbb R\)(脚注 9.4:\(\bar D=D\cup\partial D\) 是 \(D\) 的闭包)使得:\(f\) 在 \(\bar D\) 上连续;\(f(\mathbf x)=F(\mathbf x)\) 对 \(\forall\mathbf x\in\partial D\);\(f\) 在 \(D\) 上调和。命题 9.1:Dirichlet 问题的解 \(f\) 是唯一的。The Dirichlet problem: let \(D\subset\mathbb R^d\) be the domain (open and connected) and \(F:\partial D\to\mathbb R\) continuous and bounded. Find a bounded function \(f:\bar D\to\mathbb R\) (footnote 9.4: \(\bar D=D\cup\partial D\) is the closure of \(D\)) such that: \(f\) is continuous on \(\bar D\); \(f(\mathbf x)=F(\mathbf x)\) for all \(\mathbf x\in\partial D\); \(f\) is harmonic on \(D\). Proposition 9.1: the solution \(f\) to the Dirichlet problem is unique.
命题 9.1 证明 / Proof of Proposition 9.1 先证:若 \(F(\mathbf x)=0\) 对 \(\forall\mathbf x\in\partial D\),则 \(f(\mathbf x)=0\) 对 \(\forall\mathbf x\in\bar D\)。这是因为:由调和性、连续性与零边界条件,\(f\) 在 \(\bar D\) 上不能取严格单调(沿任意路径),故 \(f\) 在 \(\bar D\) 内部有局部极大值;由最大值原理(Thm 9.1)\(f\) 是常数,即 \(f(\mathbf x)=0\)。再设有两个解 \(f,g\) 满足 \(f(\mathbf x)=F(\mathbf x)=g(\mathbf x)\) 于 \(\partial D\),则 \(h(\mathbf x)\equiv f(\mathbf x)-g(\mathbf x)=0\) 于 \(\partial D\);\(h\) 是调和函数的线性组合,故也调和;由上半部分,\(h(\mathbf x)=0\) 对 \(\forall\mathbf x\in\bar D\),即 \(f,g\) 处处相同。\(\blacksquare\)First: if \(F(\mathbf x)=0\) for all \(\mathbf x\in\partial D\), then \(f(\mathbf x)=0\) for all \(\mathbf x\in\bar D\). This is because, by harmonicity, continuity and the zero boundary condition, \(f\) cannot be strictly monotonic along any path in \(\bar D\), so \(f\) has a local maximum in the interior of \(\bar D\); by the maximum principle (Thm 9.1), \(f\) is constant, i.e. \(f(\mathbf x)=0\). Now suppose there exist two solutions \(f,g\) satisfying \(f(\mathbf x)=F(\mathbf x)=g(\mathbf x)\) on \(\partial D\); then \(h(\mathbf x)\equiv f(\mathbf x)-g(\mathbf x)=0\) on \(\partial D\); \(h\) is a linear combination of harmonic functions, hence also harmonic; by the first half, \(h(\mathbf x)=0\) for all \(\mathbf x\in\bar D\), i.e. \(f\) and \(g\) are the same everywhere. \(\blacksquare\)
9.3 性质:常返 vs 暂留 / Properties: Recurrence vs Transience
9.3.1 退出概率 / Exit Probability
球壳上的退出概率 / Exit probability on a spherical shell
\(\{\mathbf B_t\}\) 是标准 \(d\) 维布朗运动。设 \(D=\{\mathbf x\in\mathbb R^d:r<|\mathbf x|
\(f\) 是 Dirichlet 问题的唯一解 / \(f\) is the unique Dirichlet solution
\(f(\mathbf x_0)\) 自然是连续的调和函数(最终可验证)。它是如下 Dirichlet 问题的解:\(D=\{r<|\mathbf x|
求解 \(\varphi\):化为 ODE 并代入边界条件 / Solving \(\varphi\): reducing to an ODE and applying boundary conditions 由 \(f\) 调和,\(\Delta f=\sum_{i}\partial^2 f/(\partial x_i)^2=0\),即 \(\Delta\varphi(|\mathbf x|^2)=0\) (9.1)。对 \(j=1,\dots,d\):\(\dfrac{\partial\varphi(|\mathbf x|^2)}{\partial x_j}=2x_j\varphi'(|\mathbf x|^2)\),\(\dfrac{\partial^2\varphi(|\mathbf x|^2)}{(\partial x_j)^2}=4x_j^2\varphi''(|\mathbf x|^2)+2\varphi'(|\mathbf x|^2)\) (9.2)。代入 (9.1):Since \(f\) is harmonic, \(\Delta f=\sum_i\partial^2 f/(\partial x_i)^2=0\), i.e. \(\Delta\varphi(|\mathbf x|^2)=0\) (9.1). For \(j=1,\dots,d\): \(\dfrac{\partial\varphi(|\mathbf x|^2)}{\partial x_j}=2x_j\varphi'(|\mathbf x|^2)\), \(\dfrac{\partial^2\varphi(|\mathbf x|^2)}{(\partial x_j)^2}=4x_j^2\varphi''(|\mathbf x|^2)+2\varphi'(|\mathbf x|^2)\) (9.2). Substituting into (9.1):
$$\sum_{j=1}^d\!\left(4x_j^2\varphi''+2\varphi'\right)=4|\mathbf x|^2\varphi''+2d\varphi'=0\;\Rightarrow\;4y\varphi''(y)+2d\varphi'(y)=0\;\Rightarrow\;\varphi''(y)=-\frac d{2y}\varphi'(y)\quad(y=|\mathbf x|^2).$$
该 ODE 关于 \(\varphi'\) 的解为 \(\varphi'(y)=Cy^{-d/2}\)。分维数积分并代入边界 \(\varphi(R^2)=1\)、\(\varphi(r^2)=0\):This ODE in \(\varphi'\) has solution \(\varphi'(y)=Cy^{-d/2}\). Integrating per dimension and applying the boundaries \(\varphi(R^2)=1\), \(\varphi(r^2)=0\):
- \(d=1\):\(\varphi(y)=C_1|\mathbf x|+D_1\),\(C_1=\dfrac1{R-r}\),\(D_1=\dfrac{-r}{R-r}\),故 \(f(\mathbf x)=\dfrac{|\mathbf x|-r}{R-r}\)。\(d=1\): \(\varphi(y)=C_1|\mathbf x|+D_1\), \(C_1=\dfrac1{R-r}\), \(D_1=\dfrac{-r}{R-r}\), so \(f(\mathbf x)=\dfrac{|\mathbf x|-r}{R-r}\).
- \(d=2\):\(\varphi(y)=2C_2\ln|\mathbf x|+D_2\),\(2C_2=\dfrac1{\ln R-\ln r}\),\(D_2=\dfrac{-\ln r}{\ln R-\ln r}\),故 \(f(\mathbf x)=\dfrac{\ln|\mathbf x|-\ln r}{\ln R-\ln r}\)。\(d=2\): \(\varphi(y)=2C_2\ln|\mathbf x|+D_2\), \(2C_2=\dfrac1{\ln R-\ln r}\), \(D_2=\dfrac{-\ln r}{\ln R-\ln r}\), so \(f(\mathbf x)=\dfrac{\ln|\mathbf x|-\ln r}{\ln R-\ln r}\).
- \(d\ge3\):\(\varphi(y)=C_3|\mathbf x|^{2-d}+D_3\),\(C_3=\dfrac1{R^{2-d}-r^{2-d}}\),\(D_3=\dfrac{-r^{2-d}}{R^{2-d}-r^{2-d}}\),故 \(f(\mathbf x)=\dfrac{|\mathbf x|^{2-d}-r^{2-d}}{R^{2-d}-r^{2-d}}\)。\(\blacksquare\)\(d\ge3\): \(\varphi(y)=C_3|\mathbf x|^{2-d}+D_3\), \(C_3=\dfrac1{R^{2-d}-r^{2-d}}\), \(D_3=\dfrac{-r^{2-d}}{R^{2-d}-r^{2-d}}\), so \(f(\mathbf x)=\dfrac{|\mathbf x|^{2-d}-r^{2-d}}{R^{2-d}-r^{2-d}}\). \(\blacksquare\)
退出概率公式 (9.4) / Exit probability formula (9.4)
综上,对 \(\mathbf x_0\in D=\{r<|\mathbf x|
$$f(\mathbf x_0)=\mathbb P\{|\mathbf B_T|=R\mid\mathbf B_0=\mathbf x_0\}=\begin{cases}\dfrac{|\mathbf x_0|-r}{R-r}&\text{if }d=1\\[2mm]\dfrac{\ln|\mathbf x_0|-\ln r}{\ln R-\ln r}&\text{if }d=2\\[2mm]\dfrac{|\mathbf x_0|^{2-d}-r^{2-d}}{R^{2-d}-r^{2-d}}&\text{if }d\ge3\end{cases}\tag{9.4}$$
它确实球对称、连续、调和(脚注 9.5:调和性来自我们在求解 \(f\) 时施加了 \(f\) 调和的限制)。which is indeed spherically symmetric, continuous, and harmonic (footnote 9.5: harmonicity follows as we imposed the restriction of \(f\) being harmonic when solving for \(f\)).
9.3.2 常返与暂留的定义 / Definitions of Recurrent and Transient
定义 9.6–9.8 / Definitions 9.6–9.8 定义 9.6(点常返):\(\{\mathbf B_t\}\) 点常返,若 \(\forall\mathbf x\in\mathbb R^d\),存在随机序列 \(t_n\uparrow\infty\) 使得以概率 1,\(\mathbf B_{t_n}=\mathbf x\) 对 \(\forall n\in\mathbb N_+\)。定义 9.7(邻域常返):\(\{\mathbf B_t\}\) 邻域常返,若 \(\forall\mathbf x\in\mathbb R^d\)、\(\forall\varepsilon>0\),存在随机序列 \(t_n\uparrow\infty\) 使得以概率 1,\(|\mathbf B_{t_n}-\mathbf x|<\varepsilon\) 对 \(\forall n\in\mathbb N_+\)。定义 9.8(暂留):\(\{\mathbf B_t\}\) 暂留,若以概率 1,\(\lim_{t\to\infty}|\mathbf B_t|=\infty\)。Definition 9.6 (Point recurrent): \(\{\mathbf B_t\}\) is point recurrent if for all \(\mathbf x\in\mathbb R^d\) there is a random sequence \(t_n\uparrow\infty\) such that with probability 1, \(\mathbf B_{t_n}=\mathbf x\) for all \(n\in\mathbb N_+\). Definition 9.7 (Neighborhood recurrent): \(\{\mathbf B_t\}\) is neighborhood recurrent if for all \(\mathbf x\in\mathbb R^d\) and all \(\varepsilon>0\) there is a random sequence \(t_n\uparrow\infty\) such that with probability 1, \(|\mathbf B_{t_n}-\mathbf x|<\varepsilon\) for all \(n\in\mathbb N_+\). Definition 9.8 (Transient): \(\{\mathbf B_t\}\) is transient if with probability 1, \(\lim_{t\to\infty}|\mathbf B_t|=\infty\).
9.3.3 一维布朗运动:点常返 / One-Dimensional BM: Point Recurrent
定理 9.2 / Theorem 9.2 \(\{B_t\}\) 是标准一维布朗运动,则它是点常返的。\(\{B_t\}\) is a standard one-dimensional Brownian motion; then it is point recurrent.
定理 9.2 证明 / Proof of Theorem 9.2
对 \(0
9.3.4 二维布朗运动:邻域常返但非点常返 / Two-Dimensional BM: Neighborhood Recurrent But Not Point Recurrent
定理 9.3 / Theorem 9.3 \(\{\mathbf B_t\}\) 是标准二维布朗运动,则它是邻域常返但非点常返的。\(\{\mathbf B_t\}\) is a standard two-dimensional Brownian motion; then it is neighborhood recurrent but not point recurrent.
定理 9.3 证明 / Proof of Theorem 9.3
邻域常返:对 \(0
9.3.5 三维及以上:暂留 / Three or Higher Dimensional BM: Transient
定理 9.4 / Theorem 9.4 \(\{\mathbf B_t\}\) 是标准三维或更高维布朗运动,则它是暂留的。\(\{\mathbf B_t\}\) is a standard three- or higher-dimensional Brownian motion; then it is transient.
定理 9.4 证明 / Proof of Theorem 9.4
对 \(0
$$\lim_{R\to\infty}\mathbb P\{|\mathbf B_T|=r\mid\mathbf B_0=\mathbf x_0\}=\left(\frac r{|\mathbf x_0|}\right)^{d-2}.\tag{9.7}$$
故 \(\{\mathbf B_t\}\) 以概率 \(\left(r/|\mathbf x_0|\right)^{d-2}\) 在跑向无穷前到达半径 \(r\) 的小球。一旦到达,由定理 8.1,以概率 1 \(\{\mathbf B_t\}\) 会在足够大的 \(t\) 回到某 \(\mathbf x_1\)(\(|\mathbf x_1|=|\mathbf x_0|\));再由增量独立性,第二次到达小球的概率仍为 \(\left(r/|\mathbf x_0|\right)^{d-2}\)。记在跑向无穷前到达小球 \(N\) 次的事件为 \(A_N\),概率 \(p_N=\left(r/|\mathbf x_0|\right)^{N(d-2)}\),则 \(\lim_{N\to\infty}p_N=0\)。故其补事件 \(A_N^c\)(在到达 \(N\) 次小球前就偏离到无穷)概率为 1(\(\lim_{N\to\infty}1-p_N=1\))。这对任意 \(r\) 成立。故对足够大的 \(t\),\(\{\mathbf B_t\}\) 以概率 1 偏离到无穷,即 \(\{\mathbf B_t\}\) 暂留。\(\blacksquare\)So \(\{\mathbf B_t\}\) has probability \(\left(r/|\mathbf x_0|\right)^{d-2}\) of reaching the small ball of radius \(r\) before going to infinity. Once it reaches that ball, by Theorem 8.1, with probability 1 \(\{\mathbf B_t\}\) returns to some \(\mathbf x_1\) with \(|\mathbf x_1|=|\mathbf x_0|\) at sufficiently large \(t\); then by increment independence, the second time it has probability \(\left(r/|\mathbf x_0|\right)^{d-2}\) of reaching the ball. Denoting by \(A_N\) the event of reaching the small ball \(N\) times before going to infinity, with probability \(p_N=\left(r/|\mathbf x_0|\right)^{N(d-2)}\), we have \(\lim_{N\to\infty}p_N=0\). So its complement \(A_N^c\) (deviating to infinity before reaching the ball \(N\) times) has probability 1 (\(\lim_{N\to\infty}1-p_N=1\)). This is true for any \(r\). So for large enough \(t\), \(\{\mathbf B_t\}\) deviates to infinity with probability 1, i.e. \(\{\mathbf B_t\}\) is transient. \(\blacksquare\)
9.4 Green 函数 / Green's Function
9.4.1 转移密度 / Transition Density
转移密度 \(P_t(\mathbf x,\mathbf y)\) (9.8) 与注 9.1 / Transition density and Remark 9.1 由事实 9.1,标准 \(d\) 维 BM 增量 \(\mathbf B_t-\mathbf B_s\) 的密度为 \(f(\mathbf x)=\dfrac1{[2\pi(t-s)]^{d/2}}e^{-\frac{|\mathbf x|^2}{2(t-s)}}\)。可等价地定义转移密度 \(P_t(\mathbf x,\mathbf y)\):\(\mathbf x\) 是起点(\(\mathbf B_0=\mathbf x\)),\(P_t(\mathbf x,\mathbf y)\) 是时刻 \(t\) 处 \(\mathbf y\) 的密度,即By Fact 9.1, the density of the increment \(\mathbf B_t-\mathbf B_s\) of a standard \(d\)-dimensional BM is \(f(\mathbf x)=\dfrac1{[2\pi(t-s)]^{d/2}}e^{-\frac{|\mathbf x|^2}{2(t-s)}}\). We can equivalently define the transition density \(P_t(\mathbf x,\mathbf y)\): \(\mathbf x\) is the start (\(\mathbf B_0=\mathbf x\)) and \(P_t(\mathbf x,\mathbf y)\) is the density of \(\mathbf y\) at time \(t\), i.e.
$$P_t(\mathbf x,\mathbf y)=\frac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}.\tag{9.8}$$
注 9.1:\(P_t(\mathbf x,\mathbf y)\) 可直觉理解为 \(\{\mathbf B_t\}\) 从 \(\mathbf x\) 出发、于时刻 \(t\) 到达 \(\mathbf y\) 的概率。Remark 9.1: \(P_t(\mathbf x,\mathbf y)\) can be intuitively interpreted as the probability of \(\{\mathbf B_t\}\) starting at \(\mathbf x\) and reaching \(\mathbf y\) at time \(t\).
9.4.2 热方程 / Heat Equation
定义 9.9(热方程)与命题 9.2 / Definition 9.9 (Heat equation) and Proposition 9.2 定义 9.9(热方程):函数 \(f(\mathbf x,\mathbf y,t)\) 满足热方程,若 \(\dfrac\partial{\partial t}f(\mathbf x,\mathbf y,t)=\alpha\big(\Delta_{\mathbf x,\mathbf y}f(\mathbf x,\mathbf y,t)\big)\),其中 \(\Delta_{\mathbf x,\mathbf y}\) 是关于 \(\mathbf x,\mathbf y\) 的 Laplace 算子。命题 9.2:(9.8) 中的转移密度满足热方程,即 \(\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf x}P_t(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf y}P_t(\mathbf x,\mathbf y)\)。Definition 9.9 (Heat equation): a function \(f(\mathbf x,\mathbf y,t)\) satisfies the heat equation if \(\dfrac\partial{\partial t}f(\mathbf x,\mathbf y,t)=\alpha\big(\Delta_{\mathbf x,\mathbf y}f(\mathbf x,\mathbf y,t)\big)\), where \(\Delta_{\mathbf x,\mathbf y}\) is the Laplace operator w.r.t. \(\mathbf x,\mathbf y\). Proposition 9.2: the transition density in (9.8) satisfies the heat equation, i.e. \(\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf x}P_t(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf y}P_t(\mathbf x,\mathbf y)\).
命题 9.2 证明 / Proof of Proposition 9.2 对 (9.8) 关于 \(t\) 求导:\(\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)=\left(\dfrac{|\mathbf x-\mathbf y|^2}{2t^2}-\dfrac d{2t}\right)\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}\)。对 \(\mathbf x\) 中 \(x_i\) 求二阶导:\(\dfrac{\partial^2}{\partial x_i^2}P_t(\mathbf x,\mathbf y)=\left(\dfrac{(x_i-y_i)^2}{t^2}-\dfrac1t\right)\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}\)。对 \(i\) 求和:Differentiate (9.8) w.r.t. \(t\): \(\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)=\left(\dfrac{|\mathbf x-\mathbf y|^2}{2t^2}-\dfrac d{2t}\right)\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}\). The second derivative w.r.t. \(x_i\) in \(\mathbf x\): \(\dfrac{\partial^2}{\partial x_i^2}P_t(\mathbf x,\mathbf y)=\left(\dfrac{(x_i-y_i)^2}{t^2}-\dfrac1t\right)\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}\). Summing over \(i\):
$$\Delta_{\mathbf x}P_t(\mathbf x,\mathbf y)=\sum_{i=1}^d\frac{\partial^2}{\partial x_i^2}P_t=\frac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x-\mathbf y|^2}{2t}}\left(\frac{|\mathbf x-\mathbf y|^2}{t^2}-\frac dt\right)\;\Rightarrow\;\frac12\Delta_{\mathbf x}P_t=\frac\partial{\partial t}P_t.$$
由 \(\mathbf x,\mathbf y\) 的对称性,同样有 \(\dfrac12\Delta_{\mathbf y}P_t(\mathbf x,\mathbf y)=\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)\)。\(\blacksquare\)By the symmetry of \(\mathbf x,\mathbf y\), likewise \(\dfrac12\Delta_{\mathbf y}P_t(\mathbf x,\mathbf y)=\dfrac\partial{\partial t}P_t(\mathbf x,\mathbf y)\). \(\blacksquare\)
9.4.3 Green 函数 / Green's Function
定义 9.10(Green 函数)与事实 9.2 / Definition 9.10 (Green's function) and Fact 9.2 定义 9.10(Green 函数):Green 函数 \(G(\mathbf x,\mathbf y)\) 定义为从 \(\mathbf x\) 出发、在整个未来访问 \(\mathbf y\) 的期望次数,即 \(G(\mathbf x,\mathbf y)=\int_0^\infty P_t(\mathbf x,\mathbf y)\,dt\)。事实 9.2:由 \(P_t(\mathbf x,\mathbf y)\) 中 \(\mathbf x,\mathbf y\) 的对称性,\(G(\mathbf x,\mathbf y)=G(\mathbf y,\mathbf x)\);由 \(P_t\) 的平移性质,\(G(\mathbf x,\mathbf y)=G(\mathbf y-\mathbf x)=G(\mathbf x-\mathbf y)\);特别地 \(G(\mathbf x)=G(\mathbf 0,\mathbf x)\)。Definition 9.10 (Green's function): the Green's function \(G(\mathbf x,\mathbf y)\) is defined by the expected number of visits to \(\mathbf y\) in the entire future starting from \(\mathbf x\), i.e. \(G(\mathbf x,\mathbf y)=\int_0^\infty P_t(\mathbf x,\mathbf y)\,dt\). Fact 9.2: by the symmetry of \(\mathbf x,\mathbf y\) in \(P_t(\mathbf x,\mathbf y)\), \(G(\mathbf x,\mathbf y)=G(\mathbf y,\mathbf x)\); by the translation property of \(P_t\), \(G(\mathbf x,\mathbf y)=G(\mathbf y-\mathbf x)=G(\mathbf x-\mathbf y)\); in particular \(G(\mathbf x)=G(\mathbf 0,\mathbf x)\).
计算 \(G(\mathbf x)\):化为 Gamma 积分 / Computing \(G(\mathbf x)\): reduction to a Gamma integral \(G(\mathbf x)=G(\mathbf 0,\mathbf x)=\displaystyle\int_0^\infty P_t(\mathbf 0,\mathbf x)\,dt=\int_0^\infty\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x|^2}{2t}}\,dt\) (9.9)。令 \(u=\dfrac{|\mathbf x|^2}{2t}\),则 \(t=\dfrac{|\mathbf x|^2}{2u}\)、\(dt=-\dfrac{|\mathbf x|^2}{2u^2}\,du\),(9.9) 化为\(G(\mathbf x)=G(\mathbf 0,\mathbf x)=\displaystyle\int_0^\infty P_t(\mathbf 0,\mathbf x)\,dt=\int_0^\infty\dfrac1{(2\pi t)^{d/2}}e^{-\frac{|\mathbf x|^2}{2t}}\,dt\) (9.9). Let \(u=\dfrac{|\mathbf x|^2}{2t}\), then \(t=\dfrac{|\mathbf x|^2}{2u}\), \(dt=-\dfrac{|\mathbf x|^2}{2u^2}\,du\), and (9.9) becomes
$$G(\mathbf x)=\int_0^\infty\frac1{2\pi^{d/2}}|\mathbf x|^{-d}\,u^{\frac d2-1}e^{-u}\,du\,|\mathbf x|^{2}=|\mathbf x|^{2-d}\frac1{2\pi^{d/2}}\int_0^\infty u^{\frac d2-1-1}e^{-u}\,du=|\mathbf x|^{2-d}\frac{\Gamma\!\left(\frac d2-1\right)}{2\pi^{d/2}},\tag{9.10}$$
末式由 Gamma 函数定义 \(\Gamma(z)\equiv\int_0^\infty u^{z-1}e^{-u}\,du\)。故一般地 \(G(\mathbf x,\mathbf y)=C_d|\mathbf x-\mathbf y|^{2-d}\),常数 \(C_d=\dfrac{\Gamma\!\left(\frac d2-1\right)}{2\pi^{d/2}}\) 与维数 \(d\) 有关。\(\blacksquare\)where the last line uses the Gamma function definition \(\Gamma(z)\equiv\int_0^\infty u^{z-1}e^{-u}\,du\). So in general \(G(\mathbf x,\mathbf y)=C_d|\mathbf x-\mathbf y|^{2-d}\) with the dimension-dependent constant \(C_d=\dfrac{\Gamma\!\left(\frac d2-1\right)}{2\pi^{d/2}}\). \(\blacksquare\)
9.4.4 条件 Green 函数 / Conditional Green's Function
条件转移概率 \(P_t^D\) 及其性质 / Conditional transition probability \(P_t^D\) and its properties 把停时嵌入转移概率与 Green 函数。设 \(D\) 是定义域,标准 \(d\) 维 BM \(\{\mathbf B_t\}\) 从 \(\mathbf B_0=\mathbf x\in D\) 出发,停时 \(\tau=\inf\{t:\mathbf B_t\notin D\}\)。定义条件转移概率 \(P_t^D(\mathbf x,\mathbf y)\) 为 \(\{\mathbf B_t\}\) 从 \(\mathbf x\) 出发、在离开 \(D\) 之前于时刻 \(t\) 到达 \(\mathbf y\) 的概率。由定义 \(\mathbb P\{t<\tau\mid\mathbf B_0=\mathbf x\in D\}=\int_D P_t^D(\mathbf x,\mathbf y)\,d\mathbf y\)。其性质:(i) 与无条件的关系 \(P_t^D(\mathbf x,\mathbf y)=P_t(\mathbf x,\mathbf y)-\mathbb E[P_{t-\tau}(\mathbf B_\tau,\mathbf y)\mid t>\tau\text{ and }\mathbf B_0=\mathbf x\in D]\)(直觉:从总路径中减去"先离开再返回"的不合格路径);(ii) 对称性 \(P_t^D(\mathbf x,\mathbf y)=P_t^D(\mathbf y,\mathbf x)\);(iii) 热方程 \(\dfrac\partial{\partial t}P_t^D(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf x}P_t^D(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf y}P_t^D(\mathbf x,\mathbf y)\)。Embed a stopping time into the transition probability and Green's function. Let \(D\) be the domain, a standard \(d\)-dimensional BM \(\{\mathbf B_t\}\) start at \(\mathbf B_0=\mathbf x\in D\), and \(\tau=\inf\{t:\mathbf B_t\notin D\}\). Define the conditional transition probability \(P_t^D(\mathbf x,\mathbf y)\) as the probability of \(\{\mathbf B_t\}\) starting at \(\mathbf x\) and reaching \(\mathbf y\) at time \(t\) before ever leaving \(D\). By definition \(\mathbb P\{t<\tau\mid\mathbf B_0=\mathbf x\in D\}=\int_D P_t^D(\mathbf x,\mathbf y)\,d\mathbf y\). Properties: (i) relation to the unconditional \(P_t^D(\mathbf x,\mathbf y)=P_t(\mathbf x,\mathbf y)-\mathbb E[P_{t-\tau}(\mathbf B_\tau,\mathbf y)\mid t>\tau\text{ and }\mathbf B_0=\mathbf x\in D]\) (intuitively: subtract the "leave then return" unqualified paths from the total); (ii) symmetry \(P_t^D(\mathbf x,\mathbf y)=P_t^D(\mathbf y,\mathbf x)\); (iii) heat equation \(\dfrac\partial{\partial t}P_t^D(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf x}P_t^D(\mathbf x,\mathbf y)=\dfrac12\Delta_{\mathbf y}P_t^D(\mathbf x,\mathbf y)\).
条件 Green 函数 (9.11) / Conditional Green's function (9.11) 条件 Green 函数 \(G^D(\mathbf x,\mathbf y)\) 是从 \(\mathbf x\) 出发、在离开 \(D\) 之前于整个未来访问 \(\mathbf y\) 的期望次数:The conditional Green's function \(G^D(\mathbf x,\mathbf y)\) is the expected number of visits to \(\mathbf y\) in the entire future starting from \(\mathbf x\) before ever leaving \(D\):
$$G^D(\mathbf x,\mathbf y)=\int_0^\infty P_t^D(\mathbf x,\mathbf y)\,dt=G(\mathbf x,\mathbf y)-\mathbb E[G(\mathbf B_\tau,\mathbf y)\mid\mathbf B_0=\mathbf x\in D].\tag{9.11}$$
直觉:条件 Green 函数 = 期望总访问次数 \(G(\mathbf x,\mathbf y)\) 减去不合格(先离开再返回)的期望访问次数 \(\mathbb E[G(\mathbf B_\tau,\mathbf y)\mid\mathbf B_0=\mathbf x\in D]\)。Intuitively: the conditional Green's function $=$ the expected total visits \(G(\mathbf x,\mathbf y)\) minus the expected unqualified (leave-and-return) visits \(\mathbb E[G(\mathbf B_\tau,\mathbf y)\mid\mathbf B_0=\mathbf x\in D]\).