中国科学院大学数学科学学院, 北京 100049
2015年12月30日 收稿; 2016年05月20日 收修改稿
基金项目: 国家自然科学基金(11171342)资助
通信作者: 胡晓予, E-mail:xyhu@ucas.ac.cn
摘要: 考虑变化环境中分枝树上的广义有偏随机游动,分别得到随机游动为暂留的、正常返的和零常返的充分条件,为进一步研究这类随机游动的中心极限定理等性质做了铺垫。
关键词: 树上随机游动变化环境中分枝过程常返性暂留性
Classification of states for a biased random walk on a branching tree in the varying environment
ZHANG Zhiyang, HU Xiaoyu
School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: In this study we investigate biased random walk on branching trees in the varying environment. We have obtained the sufficient conditions for transience, positive recurrence, and null recurrence. Our main result enables the further study on the central limit theory for this kind of random walks.
Key words: random walk on a treebranching process in the varying environmentrecurrencetransience
随机游动是随机过程中的一个经典模型。早在概率论发展初期,就有人曾经考虑过直线上的简单随机游动和有偏随机游动等问题[1]。如今这方面的研究对象已经从直线上的随机游动扩展到随机图和网络上的随机游动, 以及随机环境中的随机游动和各种代数结构上的随机游动[2-4]。
随机游动研究中的一个基本问题是状态分类问题。历史上Kolmogorov, Polya, Kakutani等都曾经讨论过这一问题, Polya还给出了
2006年Collevecchio[9]利用另一种方法, 在Lyons工作的基础上得到经典分枝树上的λ-有偏随机游动的状态分类结果。我们则考虑变化环境中的分枝树上的广义有偏随机游动, 得到其状态分类的结果,文献[9]中的结果恰为我们的结论的特例。
1 背景知识设(Ω,
${c_1}g\left( x \right) \le f\left( x \right) \le {c_2}g\left( x \right), $ |
定义1.1?? (树)设U为由如下形式的有限正整数序列所构成的集合: < i1, i2, …, in > , 其中ik∈
(a)空序列?∈T;
(b)若对某个j∈N*, < i1, i2, …, ik, j > ∈T,则必有 < i1, …, ik > ∈T;
(c)若u= < i1, i2, …, ik > ∈T, νu(T)是依赖于u和T的某正整数(也称νu(T)为T上节点u的分枝数),则对于所有的1≤j≤νu(T),总有 < i1, i2, …, ik,j > ∈T;
则称T为一棵树,空序列称为树的根, 以o记之。称T中全体序列长度的上确界为树T的高度。若v= < i > ∈T,则称
${T^v}:{\rm{ = }}\{ < i, {i_1}, {i_2}, \cdots, {i_n} > < i, {i_1}, {i_2}, \cdots, {i_n} > \in T{\rm{\} }}$ |
$T{|_n}{\rm{: = \{ }} < {i_1}, {i_2}, \cdots, {i_k} > \in T;k \le n\} .$ |
记
定义1.2(G-W过程和分枝树)令Z0=1且对于任意的n≥0,
${z_{n + 1}} = \left\{ {\begin{array}{*{20}{c}}{\sum\limits_{j = 1}^{{Z_n}} {{X_{n, j}}, } }&{{Z_n} \ne 0, }\\{0, \;\;\;\;\;\;\;\;}&{{Z_n} = 0, }\end{array}} \right.$ |
定义1.3 ??(G-W测度)设{p(k), k≥0}为某Galton-Watson过程的分布律,
$\begin{array}{c}GW\left\{ {T \in F:{v_j} \in T, 1 \le j \le n, \mathop \cap \limits_{j = 1}^n \left\{ {{d_{{v_j}}} = {k_j}} \right\}} \right\}\\ = \prod\limits_{j = 1}^n p \left( {{k_j}} \right).\end{array}$ |
类似地,可以定义变化环境中的分枝过程,以及变化环境中的分枝树。
定义1.4 ??令Y0=1且对于任意的n≥0,
${Y_{n + 1}} = \left\{ {\begin{array}{*{20}{c}}{\sum\limits_{j = 1}^{{Y_n}} {{U_{n + 1, j}}, \left( x \right)} }&{{y_n} \ne 0;}\\{0, \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;}&{{Y_n} = 0, }\end{array}} \right.$ |
${p^{(n)}}(0) = 0, 0 < {p^{(n)}}(1) < 1, n \ge 1.$ |
$\begin{array}{c}VGW\left\{ {T \in F:{v_j} \in T, 1 \le j \le n, \;\mathop \cap \limits_{j = 1}^n \left\{ {{d_{{v_j}}} = {k_j}} \right\}} \right\}\\ = \prod\limits_{j = 1}^n {{p^{\left( {\left| {{v_j}} \right|} \right)}}\left( {{k_j}} \right).} \end{array}$ |
定义1.6 ??给定一棵由上临界的G-W树T,其分布律为{p(k), k≥0}, 设{Xn:n≥0}为定义在
${P_T}\left( {{X_{n + 1}} = w|{X_n} = v} \right) = \left\{ \begin{array}{l}\frac{\lambda }{{\lambda + {d_v}}}, v = \overrightarrow w, \\\frac{1}{{\lambda + {d_v}}}, w = \overrightarrow v, v \ne o, \\\frac{1}{d}, v = o, w = \overrightarrow o .\end{array} \right.$ |
类似地,可以定义变化环境中分枝树上的{λn, n≥1}-有偏随机游动如下。
定义1.7 ??给定一棵上临界的变化环境中的分枝树T,其后代分布律为{p(n)(j), n≥1, j≥0},又给定一列正实数{λl, l≥1},若{Vn:n≥0}为定义在
${P_T}\left( {{V_{n + 1}} = w|{V_n} = v} \right) = \left\{ \begin{array}{l}\frac{{{\lambda _{\left| v \right|}}}}{{{\lambda _{\left| v \right|}} + {d_v}}}, v = \overrightarrow w, \\\frac{1}{{{\lambda _{\left| v \right|}} + {d_v}}}, w = \overrightarrow v, v \ne o, \\\frac{1}{{{d_o}}}, v = o, w = \overrightarrow {o.} \end{array} \right.$ |
根据Lyons[7]关于树上有偏随机游动状态分类的结论可得
对于一个个体平均生殖个数为m>1的G-W树上的λ-有偏随机游动{Xn:n≥0}而言:当λ>m时,对几乎所有的GW-树T,λ-有偏随机游动{Xn:n≥0}是正常返的;当λ=m时,对几乎所有的GW-树T, λ-有偏随机游动{Xn:n≥0}是零常返的;当λ < m时,对几乎所有的GW-树T, λ-有偏随机游动{Xn:n≥0}是暂留的。
关于变化环境中分枝树上的{λn, n≥1}-有偏随机游动,我们的主要结果如下:
定理1.1 ??设{Yn, n≥0}是一个上临界的变化环境中的分枝过程,设{mn, n≥1}如前所定义。给定一列实数λn>0(n≥1)及一棵由{Yn, n≥0}决定的分枝树T,{Xn, n≥0}是取值于T的{λn, n≥1}-有偏随机游动,则有如下结论:
1)若{λn, n≥1}与{mn, n≥1}满足:存在正整数n>1,使得
$\frac{{{m_2} \cdots {m_n}}}{{1 + {\lambda _1} + {\lambda _1}{\lambda _2} + \cdots + {\lambda _1} \cdots {\lambda _{n-1}}}} > 1$ |
$\frac{{{m}_{(k-1)n+1}}-{{m}_{kn}}}{1+\sum\limits_{l=1}^{n}{\prod\limits_{j=1}^{l}{{{\lambda }_{\left( k-1 \right)n+j-1}}}}}>1,k\ge 2,$ | (1) |
2)若{λn, n≥1}与{mn, n≥1}满足
$\sum\limits_{n=2}^{\infty }{\frac{{{m}_{1}}\cdots {{m}_{n}}}{1\text{+}\sum\limits_{l=1}^{n-1}{\prod\limits_{j=1}^{l}{{{\lambda }_{j}}}}} < \infty ,}$ | (2) |
3)若对VGW-a.s.T,随机游动{Xn, n≥0}关于PT是常返的,而且
$\sum\limits_{n=2}^{\infty }{\frac{\left| {{D}_{n}}(T) \right|}{1+\sum\nolimits_{l=1}^{n-1}{(\prod\nolimits_{j=1}^{l}{{{\lambda }_{j}}})}}=\infty ,}$ | (3) |
2 主要定理的证明命题2.1 ??VGW可唯一扩张为
证明 ??考虑乘积空间
$\times _{i = 1}^n{q^{\left( {{u_i}} \right)}}(\{ < {i_1}, {i_2}, \cdots, {i_n} > \} ) = \prod\limits_{j = 1}^n {{p^{\left( {{k_i} + 1} \right)}}\left( j \right)} .$ |
给定u∈U,设
${j_{k + 1}} \le v_{ < {j_1}, \cdots, {j_k} > }^*({w^*}), 0 \le k \le p.$ |
通篇假定对于任意的n≥1,
${m_n} = {\sum\limits_{k \ge 1} {kp} ^{(n)}}\left( k \right) > 1, $ |
受文献[9]中方法的启发,我们用变化环境中的分枝树的生殖规律与有偏随机游动的有偏参数的性质讨论状态分类。沿用文献[11]中使用的测度定义体系, 严格给出了状态分类定理的证明。
为证明定理1.1,首先需要下面一个引理。因为我们没有在文献中找到这一结果的证明,所以将详细的证明陈述如下。
引理2.1 ??设{St:t≥0}是取值于
$P\left( {{S_{t + 1}} = 1|{S_t} = 0} \right) = 1 = {p_0}, $ |
$P\left( {{S_{t + 1}} = n + 1|{S_t} = n} \right) = {p_n}, $ |
$P\left( {{S_{t + 1}} = n + 1|{S_t} = n} \right) = 1-{p_n}: = {q_n}.$ |
${\tau _B} = {\rm{inf}}\left\{ {k \ge 0:{S_k} \in B} \right\}, $ |
${\tau _j} = {\tau _{\left\{ j \right\}}}, j \in \mathbb{N}, $ |
${\delta _j} = \frac{{{q_j}}}{{{p_j}}}, j \ge 1, $ |
$P\left( {{\tau _0}.{\tau _{n + 1}}|{S_0} = 1} \right) = 1/[\sum\limits_{k = 1}^n {(\prod\limits_{j = 1}^k {{\delta _j}} ) + 1}], n \ge 1, $ |
证明
记
记ak=P(τ0>τn+1|S0=k), k≥0。
首先有
${a_0} = P({\tau _0} > {\tau _{n + 1}}|{S_0} = 0) = 0, $ |
${a_0} = P({\tau _0} > {\tau _{n + 1}}|{S_0} = 0) = 0, \;\;{a_{n + 1}} = P\left( {{\tau _0} > {\tau _{n + 1}}|{S_0} = n + 1} \right) = 1, n \ge 1.$ |
$P({\tau _0} > {\tau _{n + 1}}|{S_0} = k) = P\left( {\tau _0^{(1)} > \tau _{n + 1}^{(1)}|{S_1} = k} \right), $ |
$\begin{array}{l}{a_k} = P({\tau _0} > {\tau _{n + 1}}|{S_1} = k + 1)P\left( {{S_1} = k + 1|{S_0}{\rm{ = }}k} \right) + \\P({\tau _0} > {\tau _{n + 1}}|{S_1} = k-1)P({S_1} = k-1|{S_0} = k) = \\{P_k}{a_{k + 1}} + {q_k}{a_{k-1}}, \end{array}$ |
$\begin{array}{c}{a_{k + 1}}-{a_k} = \frac{{{q_k}}}{{{p_k}}}\left( {{a_k}-{a_{k-1}}} \right) = {\delta _k}({a_k} - {a_{k - 1}}), \\0 < k < n + 1.\end{array}$ |
${a_{k + 1}}-{a_k} = (\prod\limits_{j = 1}^k {{\delta _j}} ){a_1}, 0 < k < n + 1, $ |
$1 = {a_{n + 1}}-{a_0} = \sum\limits_{k = 0}^n {({a_{k + 1}}-{a_k}) = \sum\limits_{k = 1}^n {(\prod\limits_{j = 1}^k {{\delta _j}} ){a_1} + {a_1}, } } $ |
$[\sum\limits_{k = 1}^n {(\prod\limits_{j = 1}^k {{\delta _j}} ) + 1}]{a_1} = 1$ |
定理1.1的证明:
先证明1)。设:
${\overline X _n}(T, \omega ):(\mathscr{T} \times \Omega, \mathscr{B}(\mathscr{T}) \times \mathscr{T}) \to (U, \mathscr{G}), n \ge 1, $ |
$\tau (T, \omega ) = inf\{ j \ge 1:{\overline X _j}(T, \omega ) = o\} ;$ |
${{T}_{v}}(T,\omega )=\left\{ \begin{matrix} \inf \{k\ge 0:{{\overline{X}}_{k}}=v\}, & v\in T, \\ \infty ,\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ & v\notin T; \\\end{matrix} \right.$ |
${{H}_{v}}(T,\omega )=\left\{ \begin{matrix} \inf \{k\ge {{T}_{v}}:{{\overline{X}}_{k}}=\overleftarrow{v}\}, & v\in T, \\ \infty ,\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ & v\notin T. \\\end{matrix} \right.$ |
$\bar P(A) = \int\limits_\mathscr{F} {{P_T}(A(T))dVGW(T)}, $ |
为证结论1), 只需证明:对于VGW-a.s.T, {Xn, n≥0}是暂留的(关于PT),亦即,对VGW-a.s.T, {Xn, n≥0}在时间τ(T)之前以正概率访问了无穷多个T上的节点。
谬设P(τ < ∞)=1,此即对VGW-a.s.T, PT(τ(T) < ∞)=1。
设n是满足(1)的正整数,n>1。给定树T,满足PT(τ(T) < ∞)=1。对于树T上第n层的任一节点v,若过程{ Xn(T), n≥0}于时刻τ(T)之前访问它, 则将之涂成白色,以Γ1记第n层全体白点的个数。对于k≥2,若T上(k-1)n层的白点依次标记为v1, v2, …, vj,以Vk, i记vi,位于第kn层的白色后代的个数(vi位于第kn层的某后代被涂为白色当且仅当{ Xn(T), n≥0}于时刻Hvi之前访问它)。我们还假定T上其他层的节点绝不涂成白色。
由上述规则知道, kn层的白色节点的位于(k-1)n层的祖先一定是白色的(k≥2)。以Γk记第kn层全体白点的个数。于是
${\mathit{\Gamma }_{k + 1}} = \sum\limits_{j-1}^{_{{\mathit{\Gamma }_k}}} {{V_{k + 1, j}}, k \ge 0, } $ |
上式定义的{Γk, k≥0}是
${P_T}({H_v}(T) < \infty |{T_v}(T) < \infty ) = 1.$ |
${p_i} + {q_i} = 1, $ |
${p_i} = \frac{1}{{\lambda + {d_{{u_i}}}}}(1 + ({d_{{u_i}}}-1)p({u_i}, {{\bar u}_j})), $ |
${q_i} = \frac{\lambda }{{\lambda + {d_{{u_i}}}}}(1 + ({d_{{u_i}}}-1)p({u_i}, {{\bar u}_j})), $ |
以u0记根节点o, {Xn(T), n≥0}沿u0u1…un移动可以看成一个广义的赌徒输光问题, 此时引理2.1所定义的{Sn, n≥0}从i移至i+1可看成{Xn(T), n≥0}从ui有效下行至ui+1, {Sn, n≥0}从i移至i-1可看成{Xn(T), n≥0}从ui有效上行至ui-1,而输光概率P(τ0>τn|S0=1)即是{Xn(T), n≥0}于Hv之前击中u的概率, 故由引理2.1知
PT[{Xn(T)}于Hv前击中u|Tv < ∞]=
以|Dn(T, v)|记v的位于T上第n层的全体后代个数, 则
EP(v在第n层的后代个数|Tv < ∞)=
$\begin{array}{l}{E_{VGW}}(\left| {{D_n}(T, v) \cdot {\eta _n}} \right|) = \\\frac{{{m_2} \cdots {m_n}}}{{1 + {\lambda _1} + {\lambda _1}{\lambda _2} + \cdots + {\lambda _1} \cdots {\lambda _{n-1}}}}.\end{array}$ |
${{E}_{{\bar{P}}}}\left( {{\mathit{\Gamma }}_{1}} \right)\ge \frac{{{m}_{2}}\cdots {{m}_{n}}}{1+{{\lambda }_{1}}+{{\lambda }_{1}}{{\lambda }_{2}}+\cdots +{{\lambda }_{1}}\cdots {{\lambda }_{n-1}}}.$ | (4) |
$\frac{{{m}_{(k-1)n+1}}\cdots {{m}_{kn}}}{1+\sum\limits_{l=1}^{n}{\prod\limits_{j=1}^{l}{{{\lambda }_{\left( k-1 \right)n+j-1}}}}},$ | (5) |
现在根据1)可知,EP(Γ1)>1, EP(Vk, j)>1, ?k, ?j.这说明过程{Γk, k≥0}是一个上临界的依赖于代的分枝过程, 因此
$\bar P\left( {\mathop {\lim }\limits_{k \to \infty } {\mathit{\Gamma }_k} = \infty } \right) > 0,$ | (6) |
${{P}_{T}}(B)>0,$ |
其次证明2)。给定k≥2,对于第k层的任一节点v,事件{Tv < τ}发生,必然使得X1为v之祖先且{Xn(T), n≥0}于τ之前击中v,根据前面的计算有
${P_T}({T_v} < \tau ) \le \frac{1}{{1 + {\lambda _1} + {\lambda _1}{\lambda _2} + \cdots + {\lambda _1} \cdots {\lambda _{k-1}}}}, $ |
$\frac{{{m_1} \cdots {m_2}}}{{1 + {\lambda _1} + {\lambda _1}{\lambda _2} + \cdots + {\lambda _1} \cdots {\lambda _{k - 1}}}}.$ |
${E_T}(L|\tau (T)) = \frac{1}{2}\tau (T), $ |
${E_T}(L|\tau (T)) \le {E_T}(\{ {X_n}, n \ge 0\} 在\tau (T)前击中的节点数|\tau (T)), $ |
$\begin{array}{l}{E_{\bar p}}\left( \tau \right) \le 2{E_{\bar p}}(\{ {{\bar X}_n}, n \ge 0\} 在\tau 前击中的节点数)\\ \le 2{m_1} + 2\sum\limits_{k = 2}^\infty {\frac{1}{{1 + {\lambda _1} + {\lambda _1}{\lambda _2} + \cdots + {\lambda _1} \cdots {\lambda _{k-1}}}}, } \end{array}$ |
最后证明3)。事实上要证3), 只需要证明对VGW-a.s.T, ET[τ(T)]=∞。根据第一部分的证明知:对于VGW-a.s.T,在τ之前第k层(k≥2)被{Xn, n≥0}击中的节点个数的期望(关于PT)是
$\frac{{\left| {{D_k}(T)} \right|}}{{1 + \sum\nolimits_{l = 1}^{k-1} {(\prod\nolimits_{j = 1}^l {{\lambda _j}} )} }}, $ |
${E_T}(\tau (T)) \ge \sum\limits_{k = 2}^\infty {\frac{{\left| {{D_k}(T)} \right|}}{{1 + \sum\nolimits_{l = 1}^{k-1} {(\prod\nolimits_{j = 1}^l {{\lambda _j}} )} }}, } $ |
注2.1 ??当
$\begin{align} & VGW(\left| {{D}_{k}}\left( T \right) \right|)\le \frac{{{m}_{1}}\cdots {{m}_{k}}}{{{k}^{\frac{1}{2}\left( 1+\epsilon \right)}}}\le \\ & \frac{1}{{{k}^{1+\epsilon}}}{{E}_{VGW}}{{\left[\frac{{{m}_{1}}\cdots {{m}_{k}}}{\left| {{D}_{k}}(T) \right|} \right]}^{2}}, \le \frac{C}{{{k}^{1+\epsilon}}}, \\ \end{align}$ |
$\begin{align} & \sum\limits_{n=2}^{\infty }{\frac{\left| {{D}_{n}}(T) \right|}{1+\sum\nolimits_{l=1}^{n-1}{(\prod\nolimits_{j=1}^{l}{{{\lambda }_{j}}})}}}\ge \\ & K\sum\limits_{n=2}^{\infty }{\left( \frac{{{m}_{1}}\cdots {{m}_{n}}}{1+{{m}_{1}}+\cdots +{{m}_{1}}\cdots {{m}_{n-1}}}\cdot \frac{1}{{{n}^{\frac{1}{2}\left( 1+\epsilon \right)}}} \right),} \\ \end{align}$ | (7) |
$\begin{align} & \sum\limits_{n=2}^{\infty }{\frac{\left| {{D}_{n}}\left( T \right) \right|}{1+\sum\nolimits_{l=1}^{n-1}{\left( \prod\nolimits_{j=1}^{l}{{{\lambda }_{j}}} \right)}}}\ge \\ & \tilde{K}\sum\limits_{n=2}^{\infty }{{{n}^{\frac{1}{2}\left( 1+\epsilon \right)}}},=\infty , \\ \end{align}$ | (8) |
注2.2 ??考虑经典分枝树上的λn-有偏随机游动{Xn, n≥0},假设分枝树的平均生殖个数为m∈(1,
注2.3 ??根据注2.1和注2.2, 当λn=mn+1(n≥1)时或λn↓m时,可以考虑变化环境中分枝树上的中心极限定理。我们猜测有类似于文献[11]中的结论。
参考文献
[1] | 威廉·费勒.概率论及其应用[M]. 3版.胡迪鹤译.北京:人民邮电出版社, 2006. |
[2] | Zeitouni O. Random walks in random environment, XXXI Summer school in probability, St Flour[J].Lecture notes in Math, 2004, 1837:193–312. |
[3] | Révész P. Random Walk in Random and Non-Random Environments[M]. 2nd ed. Hackensack and New York:World Scientific Publishing Co Pte Ltd, 2005. |
[4] | Wolfgang W. Random walk on infinite graphs and groups[M].Cambridge: Cambridge University Press, 2000. |
[5] | Polya G. Vber eine Anfgabe der Wahrscheinlichkeitsrechnung betref fend die Irrfahrt im Strassennetz[J].Mathematische Annalen, 1921, 84:149–160.DOI:10.1007/BF01458701 |
[6] | Doyle P, Snell J L. Random walks and electric networks[J].American Mathematical Monthly, 2000, 26(4):595–599. |
[7] | Lyons R. Random walks and percolation on trees[J].Ann Probab, 1990, 18:931–958.DOI:10.1214/aop/1176990730 |
[8] | Peres Y, Zeitouni O. A central limit theorem for biased random walks on Galton-Watson trees[J].Probab Theory Realated Fields, 2008, 140(3/4):595–629. |
[9] | Collevecchio A. On the Transience of processes defined on Galton-Watson trees[J].Ann Probab, 2006, 34:870–878.DOI:10.1214/009117905000000837 |
[10] | Neveu J. Arbres et processus de Galton-Watson[J].Ann Inst H Poincaré Probab Stastist, 1968, 22:199–207. |
[11] | Lyons R, Peres Y. Probability on trees and networks[M].Cambridge: Cambridge University Press, 2016. |