Rademacher扰动AR(1)序列的持久性与线性模1变换

《Ergodic Theory and Dynamical Systems》：Persistence of AR($1$ ) sequences with Rademacher innovations and linear mod $1$ transforms

　　

本研究考察一阶自回归（AR(1)）Markov链X_n+1=aX_n+ξ_n+1（其中a∈(0,1)为常数）在长时间内保持非负的概率（持久性概率（persistence probability））的精确渐近行为。假设独立同分布的扰动项ξ_n仅取±1两个值且a≤2/3时，研究人员得到了该概率的精确渐近公式以及X_n在保持非负条件下的弱极限。该极限分布是拟平稳的（quasi-stationary），无原子，且在1/2n=1)=1/2时为[0,3]上的均匀分布。此性质与Bernoulli卷积（Bernoulli convolution）相似。当0n出发，可在反向时间中确定性恢复被杀链的值。基于这一事实，研究人员构造了合适的Banach空间，使得被杀链的转移算子（transition operator）具有紧性（compactness），从而可应用Perron-Frobenius型论证。

当0n=-1}，故P_x(τ>n)=pⁿ，条件分布弱收敛于点1/(1-a)处的δ-测度。

核心结果涵盖a∈(1/2,2/3]情形。定义映射T_a于[0,1/(1-a)]\I_a，其中I_a=((2a-1)/(1-a),1)。该映射具有关键性质：在事件{τ>n}上，X_n-k=T_a^k(X_n)（0≤k≤n），即被杀链在反向时间中完全确定性。令?_a(x)=inf{k≥0: T_a^k(x)∈I_a}，定义示性函数δ_k(x)=1_{{Ta^k(x)≥1}}及计数函数L_k(x)=∑_i=0^k-1δ_i(x)。主要定理表明：存在常数c∈(0,1)，使得一致地有P_x(τ>n) ~ c V(x) λ_aⁿ，其中λ_a>p为方程∑_k=0^?_a δ_k (p/λ)^k+1 (q/p)^L_k =1的唯一正解，V(x)为由轨道跃迁定义的显式函数。条件分布弱收敛于拟平稳概率测度ν_a，满足P_νa(X₁∈A|τ>1)=ν_a(A)。测度ν_a无原子，除a=2/3且p=1/2时为均匀分布外，关于Lebesgue测度奇异。映射a?λ_a连续不减，在[1/2,2/3]\S的每个区间上常值，其中S={a∈[1/2,2/3]: ?_a=∞}为零Lebesgue测度的闭集。

技术方法上，研究人员采用三种途径：一是将状态聚合为有限Markov链（当轨道有限时）；二是基于紧性的算子方法，构造由T_a轨道 jumps 生成的Banach空间U，证明转移算子P在此空间上的拟紧性；三是对a=2/3情形，利用与转移算子（transfer operator）的联系。核心在于建立合适的函数空间使Perron-Frobenius定理适用，并处理p<1/2时需要的一致上界估计。

本研究发表于《Ergodic Theory and Dynamical Systems》，探讨了具有Rademacher扰动（即取±1值的独立同分布随机变量）的一阶自回归（AR(1)）Markov链的持久性问题，即该链在长时间内保持非负的概率渐近行为。

**研究背景与问题**

持久性概率的研究在随机过程理论中具有核心意义。对于一般Markov链，Hinrichs、Kolb和Wachtel此前已在对数尺度下建立了指数衰减率λ_a∈(0,1)的存在性，但其精确尾渐近需要更强的条件：创新项（innovations）的绝对连续性及矩条件。这种对连续分布的依赖源于证明方法所需的转移算子紧性。然而，离散分布创新项的情形同样重要且更具挑战，因为传统紧性论证难以直接套用。本研究针对最简单的离散分布——Rademacher分布（取±1，概率分别为p和q=1-p），且参数a∈(0,2/3]的情形，建立了完整的精确渐近理论。特别地，当a∈(0,1/2]时分析较为直接， persistence probability 为pⁿ；真正的难点在于a∈(1/2,2/3]区间。

**核心研究发现**

研究人员的首要突破是揭示了被杀AR(1)链（在首次离开[0,∞)时被终止）与经典动力系统之间的深刻联系。定义分段线性映射T_a，其关键性质在于：在反向时间中，被杀链的值可由该动力系统的轨道确定性恢复，即X_n-k=T_a^k(X_n)。这一发现将随机问题转化为确定性动力系统分析。

基于上述联系，研究人员构造了适当的Banach空间，证明被杀链的转移算子P具有拟紧性（quasi-compactness），从而应用Perron-Frobenius型论证。主要定理表明：存在常数c∈(0,1)和函数V(x)，使得对所有初始点x∈[0,1/(1-a)]，一致地有P_x(τ>n) ~ cV(x)λ_aⁿ。衰减率λ_a>p由涉及动力系统轨道的超越方程唯一确定，且条件分布弱收敛于拟平稳分布（quasi-stationary distribution）ν_a。

该拟平稳分布展现出与Bernoulli卷积相似的奇异性质：当1/2a无原子且关于Lebesgue测度奇异；仅当a=2/3且p=1/2的特殊情形下，ν_a为[0,3]上的均匀分布。映射a?λ_a连续不减，在零测闭集S外为常值，其Lebesgue-Stieltjes测度无原子且奇异。

**关键技术方法**

研究人员采用了三种互补的技术途径。第一种是状态聚合（lumping）方法：当零点的T_a-轨道有限时，可将原链聚合为有限状态Markov链，直接应用经典Perron-Frobenius定理。第二种是核心的紧性方法：通过在由T_a轨道jumps生成的函数空间U上建立算子P的拟紧性，适用于所有a∈(1/2,2/3]。该空间中的函数由其在轨道点上的跃度完全刻画，P的作用可通过显式的线性算子A在?₁空间上表示。第三种针对临界情形a=2/3，利用与加权转移算子（transfer operator）的直接联系，此时P=(3/4)P_T，其中P_T为T_2/3的Perron-Frobenius算子。

对于p<1/2的困难情形，研究人员发展了对(1/a,1/2)-展开中零频率的一致上界估计，这是保证无穷级数收敛性的关键。具体而言，通过比较零点与任意点的轨道返回时间，建立了L_n(x)/n的有效控制。

**研究结果详解**

在确定性反向时间动力学部分，研究人员严格证明了：当a∈[1/2,2/3)时，T_a将X_n的轨道唯一确定地反向恢复，且创新项可由(1/a,1/2)-展开的数字编码。这一性质直接导出了持久性概率的奇异连续特征。

在算子谱分析部分，通过将P限制于U空间，研究人员证明其等价于?₁上的线性算子A。A的领先特征值即为λ_a，对应的左右特征向量分别给出函数V(x)和归一化常数c。特别地，当轨道最终周期时，A退化为有限维矩阵，但分析方法仍然统一。

对于λ_a的性质，研究人员证明了其作为a的函数的连续性与单调性，并刻画了其支撑集的拓扑结构。利用动力系统轨道的遍历性论证，表明S集合（λ_a非常值点集）是零Lebesgue测度的闭集，这与β-变换的正规数理论密切相关。

**结论部分总结**

本研究完整地解决了Rademacher扰动AR(1)序列的持久性渐近问题，将概率方法与现代动力系统理论有机结合。核心结论包括：精确持久性渐近P_x(τ>n) ~ cV(x)λ_aⁿ；拟平稳分布ν_a的存在性、唯一性及奇异特征；以及参数λ_a的精细结构。当a=2/3且p=1/2时出现的均匀分布例外情形，反映了动力系统临界性态与概率对称性之间的深刻关联。研究人员还指出，方法可推广至更一般的双值扰动情形，核心在于保持反向时间的确定性结构。