设$ T $是一个局部有限的无限树, $ \sigma $, $ t\; (\sigma\neq t) $是$ T $中任两个顶点, 则存在唯一的从$ \sigma $到$ t $的路径, $ \sigma = z_1, z_2, \cdots, z_m = t $, 其中$ z_1, z_2, \cdots, z_m $互不相同且$ z_i, z_{i+1} $为相邻两顶点, $ m-1 $称为$ \sigma $到$ t $的距离.为给$ T $中的顶点编号, 选定一个顶点作为根顶点(简称根), 记为$ o $.如果一个顶点$ \sigma $位于根$ o $到顶点$ t $的唯一路径上, 则记$ \sigma\leq t $.若$ \sigma, t $为$ T $上不同的两个顶点, 记$ \sigma\wedge t $为同时满足$ \sigma\wedge t\leq t $和$ \sigma\wedge t\leq \sigma $且离根$ o $最远的顶点.
若$ t $为$ T $中的任一顶点, 记$ |t| $为顶点$ t $到根$ o $的距离.若$ |t| = n $, 称$ t $位于树的第$ n $层.记$ T^{(n)} $表示从根$ o $到第$ n $层所有顶点的子图, $ L_n $表示第$ n $层所有顶点的集合, $ L_m^n $表示含有$ T $的从$ m $层到$ n $层所有顶点的集合.对于任一个顶点$ t $, 从根$ o $到顶点$ t $的路径上存在唯一一个离顶点$ t $最近的顶点称为$ t $的父代, 记为$ 1_t $, 且称$ t $为$ 1_t $的子代.令$ X^S = \{X_t, t\in S\} $, $ S\subset T, $ $ x^S $为$ X^S $的实现, 且记$ |S| $为$ S $中顶点的个数.如果树图$ T $的根顶点有$ N $个相邻顶点, 而其它顶点有$ N +1 $个相邻顶点, 即$ T $的每个顶点都有$ N $个子代, 则称此树为Cayley树, 记为$ T_{C, N} $ (见图 1).
树指标随机过程是新兴的概率论研究方向. Benjamini和Peres [1]给出了树指标马氏链的定义并研究了其常返性和射线常返性; 陈晓雪和杨卫国[2]等研究了树指标马氏链的等价定义; Berger和叶中行[3]研究了齐次树图上平稳随机场熵率的存在性; 叶中行和Berger[4, 5]利用Pemantle在文献[6]中的结果及组合方法, 在依概率收敛意义下研究了齐次树图上PPG不变和遍历随机场的Shannon -McMillan定理; 杨卫国和刘文[7]研究了齐次树图上马氏链场(这实际上是树指标马氏链和PPG不变随机场的特殊情形)状态发生频率的强大数定律.近年来, 杨卫国[8, 9]研究了树上马氏链的强大数定律和渐近均分割性; 黄辉林和杨卫国[10]研究了一致有界树上马氏链的强大数定律和Shannon-McMillian定理; 石志岩和杨卫国[11]研究了树上非齐次马氏链随机转移概率调和平均的强极限性质; 王豹等[12]研究了Cayley树指标可列马氏链的强大数定律; 党慧等[13]给出了离散状态下二叉树上非齐次分支马氏链的定义, 并研究其等价性及存在性, 同时研究了有限状态下二叉树上非齐次分支马氏链的强大数定理和熵遍历定理.
随机环境中马氏链的研究已有相当长的历史. Nawrotzki[14, 15]建立了该主题的一般理论; Cogburn[16-18]构造了Hopf -链, 利用Hopf -链理论深入研究了平稳环境中马氏链的遍历理论、中心极限定理、直接收敛和转移函数的周期性关系以及不变概率测度的存在性; 胡迪鹤[19, 20]对连续时间参数的随机环境中的马氏过程的存在性、等价性、$ q $ -过程的存在唯一性进行了研究; 李应求[21, 22]利用完善的鞅差理论来研究随机环境中的马氏链, 在假设马氏双链遍历的条件下, 得到了马氏环境中马氏链的强大数定律成立的充分条件以及马氏环境中若干强极限定理; 石志岩等[23]给出了离散状态下随机环境中树指标马氏链的定义, 并证明了该定义在概率空间中可以实现; 黄辉林[24]研究了有限i.i.d.随机环境下齐次树指标马氏链的强大数定律和Shannon-McMillian定理.
本文首先给出离散状态下随机环境中树指标马氏链的定义, 并且证明了该定义在概率空间中可以实现, 给出了马氏环境下树指标马氏链与树指标马氏双链的等价定义, 该部分重述了文献[23]中的一些结果.最后研究了有限状态下马氏环境中树指标马氏链的随机转移概率调和平均的强极限性质.
设$ \Theta = \{0, 1, 2, \cdots\}, \chi = \{0, 1, 2, \cdots \} $为可列状态空间, $ \xi^{T} = \{\xi_{t}, t\in T\} $和$ X^{T} = \{X_{t}, t\in T\} $分别是$ (\Omega, {\cal F}, P) $上取值于$ \Theta $和$ \chi $的随机变量族.假定$ p_\theta = \{p(\theta;x), x\in \chi\}, \theta \in \Theta $是关于参数$ \theta $的一个分布, 且$ P_\theta = \{ p(\theta;x, y), x, y\in \chi \}, \theta \in \Theta $是定义在$ \chi^2 $上的关于参数$ \theta $的一个转移矩阵.
定义 1 [2] 设$ T $为一树图, $ \chi = \{0, 1, 2, \cdots\} $是可列状态空间, $ \{X_{t}, t\in T\} $是定义在概率空间$ (\Omega, {\cal F}, P ) $上在$ \chi $中取值的树指标变量族.设$ p = \{p(x), x\in \chi\} $是$ \chi $上一概率分布, $ P = \{p(x, y), x, y\in \chi\} $是定义在$ \chi^{2} $上的转移矩阵.如果$ \forall t\in T, \forall n\geq 1, $有
且
则称$ \{X_{t}, t\in T\} $为具有初始分布$ p $与转移矩阵$ P $在$ \chi $中取值的树指标马氏链.
注 1 文献[2]介绍了树指标马氏链的各种等价定义, 定义1仅仅是其中的一种形式, 对于其他的定义形式, 读者可参阅文献[2].
生物学家研究杆状菌的分裂时, 总结出杆状菌分裂的规律, 即一个杆状菌在分裂时, 从中间断开, 这样就分裂成两个新杆状菌, 这两个新的杆状菌为原来杆状菌的后代.如果我们把每一次分裂中的杆状菌看成一个顶点, 那么杆状菌的分裂过程就可以抽象为一个二叉树结构(分支马氏链) (见文献[25]).如果杆状菌分裂时受到周边环境的影响, 这样杆状菌的分裂过程就可以抽象为一个随机环境中二叉树模型, 因此研究随机环境中树指标过程不仅具有理论意义, 更具有应用价值.类似于定义$ 1 $中的树指标马氏链的定义, 结合随机环境中马氏链的定义, 给出随机环境中树指标马氏链的定义.
定义 2 [23] 设$ T $为树, $ X^{T} = \{X_{t}, t\in T\}, \xi^{T} = \{\xi_{t}, t\in T\} $分别在$ \chi, \Theta $中取值的随机变量族.设$ p_{\theta} = \{p(\theta;x), x\in \chi\}, \theta\in \Theta $是$ \chi $上一含参数的分布, $ P_{\theta} = \{p(\theta;x, y), x, y\in \chi\}, \theta\in \Theta $是定义在$ \chi^{2} $上的含参数的转移矩阵.若
则称$ (X^{T}, \xi^{T}) $为由含参数的分布$ p_{\theta} $与含参数的转移矩阵$ P_{\theta} $确定的随机环境中树指标马氏链, 其中$ \xi^{T} $为随机环境.若$ \xi^{T} $为树指标马氏链, 称$ (X^{T}, \xi^{T}) $为马氏环境中树指标马氏链.
注 2 当$ \xi^{T} $取常数时, 随机环境中树指标马氏链就是一般树指标马氏链.如果树的每一顶点只有一个子代, 则随机环境中树指标马氏链即为随机环境中马氏链.因此随机环境中树指标马氏链是树指标马氏链和随机环境中马氏链的推广.
引理 1 [23] $ (X^{T}, \xi^{T}) $为定义2定义的由含参数的分布$ p_{\theta} $与含参数的转移矩阵$ P_{\theta} $确定的随机环境中的树指标马氏链的充要条件是
(i)
(ii) 当$ k\geq n-1 $时,
引理 2 [23]引理$ 1 $中$ (5) $和$ (6) $式成立的充要条件是
注 3 设$ \overline{\Omega} = \chi^{T}\times \Theta^{T}, {\cal F}' $为$ \overline{\Omega} $中所有有限维柱集生成的$ \sigma $ -代数.定义$ (\overline{\Omega}, {\cal F}') $上随机环境中树指标随机过程如下$ \omega = (x^{T}, \theta^{T})\in \overline{\Omega}. $定义$ X^{T}(x^{T}, \theta^{T}) = x^{T}, \xi^{T}(x^{T}, \theta^{T}) = \theta^{T}. $定义$ \mu_{P} $在柱集$ (X^{T^{(n)}} = x^{T^{(n)}}, \xi^{T^{(m)}} = \theta^{T^{(m)}}) $上的值如$ (7) $式的右端, 由Kolmogorov存在定理(见文献[26]), $ \mu_{P} $可以扩张为整个空间$ (\overline{\Omega}, {\cal F}') $上的概率测度, 于是$ (X^{T}, \xi^{T}) $在$ \mu_{P} $下为随机环境中的树指标马氏链, 这样就得到了随机环境中树指标马氏链的实现.
注 4 若$ \xi^{T} $为初始分布为$ p^{'}(\theta), $转移矩阵族为$ K_t = (K_t(\theta, a)) $的非齐次树指标马氏链.由$ (7) $式和定义1知
若令$ Q_t(x, \theta;y, a) = p(\theta;x, y)K_t(\theta, a) $, $ q(\theta_{o}, x_{o}) = p(\theta_{o};x_{o})p^{'}(\theta_{o}) $, 则
由此可知马氏环境下树指标马氏链与树指标非齐次马氏双链是等价的.此时, $ (X^{T}, \xi^{T}) $为初始分布为$ q(\theta_{o}, x_{o}) $, 转移矩阵族为$ \{Q_t(x, \theta;y, a), t\in T\} $的树指标非齐次马氏双链.
以后总假定$ \xi^T $是马氏环境, 则$ (X^T, \xi^T) $是具有初始分布为$ q(\theta_{o}, x_{o}) $, 转移矩阵族为$ Q_t = \{Q_t(x, \theta;y, a), t\in T\} $的树指标非齐次马氏双链.在下节中, 将给出有限状态下马氏环境中树指标马氏链的随机转移概率的调和平均的极限性质.
引理 3 [24] 设$ (X^{T}, \xi^{T}) $为马氏环境中的树指标马氏链, $ \{g_{t}(x, \theta;y, \alpha), t\in T\} $是定义在$ (\chi \times\Theta)^{2} $上的函数族.令$ L_{o} = \{o\} $, $ {\cal F}_{n} = \sigma(X^{T^{(n)}}, \xi^{T^{(n)}}) $,
其中$ \lambda $为实数, 则$ \{t_{n}(\lambda, \omega), {\cal F}_{n} , n \geq1 \} $是非负鞅.
定理 1 设$ \chi = \{1, 2, \cdots, M\}, \Theta = \{1, 2, \cdots, N\} $, 且$ {\xi}^T $的初始分布为$ p^{'}(\theta) $, 转移概率族为$ K_t = \{K_t(\theta, \alpha), t\in T\} $.设$ (X^{T}, \xi^{T}) $为马氏环境中的树指标马氏链, 其初始分布和转移概率族满足
令
若存在常数$ h>0 $, 使得
则随机转移概率$ \{Q_{t}(X_{1_{t}}, \xi_{1_{t}};X_{t}, \xi_{t}), t\in T^{(n)}\backslash\{o\}\} $的调和平均$ {\rm a.e.} $收敛于$ \frac{1}{MN}, $即
注 5 假设$ \forall t\in T, b_t\equiv\frac{1}{2} $, 此时
因此等式(12)是可以实现的.
证 在引理$ 3 $中令$ g_{t}(x, \theta;y, \alpha) = Q _{t}(x, \theta;y, \alpha)^{-1}, $则由引理$ 3 $知
为非负鞅, 由Doob鞅收敛定理得
故
由式$ (14) $和$ (16) $得
由式$ (17) $和不等式$ \ln x\leq x-1 (x>0), 0\leq e^{x}-1-x \leq \frac{x^{2}}{2}e^{|x|} $得
易知当$ 0<\lambda<1 $时有
当$ 0<\lambda<h $时, 由式$ (11), (12), (18) $与$ (19) $得
令$ \lambda \to 0^{+} $, 由式$ (20) $有
当$ -h<\lambda<0 $时, 将式$ (18) $两边同除以$ \lambda $并利用式$ (11), (12) $与$ (19) $得
令$ \lambda \to 0^{-}, $由式$ (22) $有
由式$ (21) $与$ (23) $即可以知式$ (13) $成立.