2025-11-19T00:28:13.831044

Stochastic Stability of Monotone Dynamical Systems. I. The Irreducible Cooperative Systems

Jiang, Sheng, Wang
The current series of papers is concerned with stochastic stability of monotone dynamical systems by identifying the basic dynamical units that can survive in the presence of noise interference. In the first of the series, for the cooperative and irreducible systems, we will establish the stochastic stability of a dynamical order, that is, the zero-noise limit of stochastic perturbations will be concentrated on a simply ordered set consisting of Lyapunov stable equilibria. In particular, we utilize the Freidlin--Wentzell large deviation theory to gauge the rare probability in the vicinity of unordered chain-transitive invariant set on a nonmonotone manifold. We further apply our theoretic results to the stochastic stability of classical positive feedback systems by showing that the zero-noise limit is a convex combination of the Dirac measures on a finite number of asymptotically stable equilibria although such system may possess nontrivial periodic orbits.
academic

Stochastic Stability of Monotone Dynamical Systems. I. The Irreducible Cooperative Systems

基本信息

  • 论文ID: 2412.19977
  • 标题: Stochastic Stability of Monotone Dynamical Systems. I. The Irreducible Cooperative Systems
  • 作者: Jifa Jiang (河南师范大学), Xi Sheng, Yi Wang (中国科学技术大学)
  • 分类: math.DS (动力系统)
  • 发表时间: 2024年12月 (arXiv v2: 2025年11月15日)
  • 论文链接: https://arxiv.org/abs/2412.19977v2

摘要

本系列论文研究单调动力系统的随机稳定性，旨在识别在噪声干扰下能够"存活"的基本动力学单元。在本系列的第一篇中，针对合作且不可约系统，作者建立了动力学序(dynamical order)的随机稳定性，即随机扰动的零噪声极限将集中在由Lyapunov稳定平衡点组成的简单有序集上。作者特别利用Freidlin-Wentzell大偏差理论来估计非单调流形上无序链传递不变集附近的稀有概率。进一步将理论结果应用于经典正反馈系统的随机稳定性，证明尽管系统可能存在非平凡周期轨道，零噪声极限仍是有限个渐近稳定平衡点上Dirac测度的凸组合。

研究背景与动机

研究问题

单调动力系统源于M. W. Hirsch的开创性工作，其特征是存在与状态空间中闭偏序关系(由凸锥诱导)相一致的比较原理。尽管确定性单调系统的前向轨道几乎处处收敛(拓扑意义和测度论意义)，但这些系统在噪声扰动下的行为尚不清楚。

重要性

  1. 实际应用需求：源于微分方程和差分方程的数学模型经常受到噪声扰动，研究噪声如何影响单调系统动力学对理论建模和实际系统行为研究至关重要。
  2. 理论挑战：现有研究主要采用基于轨道的方法(trajectory-based approach)，在随机单调系统框架下研究。然而，这种方法可能导致动力学过度简化(如噪声同步现象)，无法充分揭示噪声扰动的影响。
  3. 方法论创新：基于分布的方法(distribution-based approach)被认为是研究一般噪声扰动的更有效策略，核心问题是识别在噪声干扰下能"存活"的基本统计动力学单元。

现有方法局限性

  1. 轨道方法的简化：随机单调系统的长期行为比确定性对应物简单得多(如噪声同步)，无法提供深入洞察。
  2. 预设目标依赖：先前关于随机稳定性的研究主要集中在预设的动力学目标(如吸引子/排斥子、拟势等价类)，缺乏对系统大尺度结构的分析。
  3. 非典型结构：单调系统既非梯度系统也非公理A系统，长期行为可能无法表示为有限个等价类的并集(由于非单调流形上可能存在无穷多个等价类)。

研究动机

本文采用结构性思想分析单调动力系统的随机稳定性，通过识别零噪声极限的集中位置，建立合作不可约系统的动力学序的随机稳定性。这是首次研究可能具有潜在无限等价类的系统的随机稳定性。

核心贡献

  1. 建立动力学序的随机稳定性：证明合作不可约系统的零噪声极限集中在由Lyapunov稳定平衡点组成的简单有序集(stationary p-arc)上(定理3.1)。
  2. 拟势的关键性质：证明对于非单调流形上的任意无序链传递集K，从K到其上(下)对偶吸引子的拟势恒为零(命题4.2和4.3)，这是技术上的关键突破。
  3. 解析/有限平衡点情形：当系统解析或仅有有限个平衡点时，零噪声极限是有限个渐近稳定平衡点上Dirac测度的凸组合(推论3.2)。
  4. 应用于正反馈系统：证明经典正反馈系统(包括Griffith型生化控制回路)的零噪声极限集中在有限个渐近稳定平衡点上，即使系统可能存在非平凡周期轨道(命题6.1)。
  5. 方法论创新：融合Benäım-Hirsch关于随机逼近算法的思想和Freidlin-Wentzell大偏差理论，直接利用大偏差理论(而非间接的离散时间方法)估计无序不变集附近的稀有概率。

方法详解

任务定义

研究合作常微分方程系统的随机扰动：

  • 确定性系统：dxdt=b(x)\frac{dx}{dt} = b(x)，其中b:Rr→Rrb: \mathbb{R}^r \to \mathbb{R}^r连续可微，满足合作性(∂bi∂xj≥0,i≠j\frac{\partial b_i}{\partial x_j} \geq 0, i \neq j)和不可约性(Jacobi矩阵不可约)。
  • 随机扰动：dXtε=b(Xtε)dt+εσ(Xtε)dWtdX^\varepsilon_t = b(X^\varepsilon_t)dt + \varepsilon\sigma(X^\varepsilon_t)dW_t，其中ε\varepsilon是小扰动参数，WtW_t是标准rr维Wiener过程。
  • 目标：确定不变测度με\mu_\varepsilon的零噪声极限μ\mu(当ε→0\varepsilon \to 0时με\mu_\varepsilon的弱极限)的支撑集。

理论框架

1. 单调动力系统基础

  • 偏序关系：x≤y⇔xi≤yi,∀ix \leq y \Leftrightarrow x_i \leq y_i, \forall i；x≪y⇔xi<yi,∀ix \ll y \Leftrightarrow x_i < y_i, \forall i
  • 强单调性：合作不可约系统的解流Φt\Phi_t满足：x<y⇒Φt(x)≪Φt(y),∀t>0x < y \Rightarrow \Phi_t(x) \ll \Phi_t(y), \forall t > 0
  • 链传递集的二分性(引理4.1)：链传递集KK要么无序，要么是静态p-arc(平衡点弧)

2. Freidlin-Wentzell大偏差理论

速率能量函数：

\int_0^T L(\phi(t), \dot{\phi}(t))dt, & \phi \in AC_T \\ \infty, & \text{otherwise} \end{cases}$$ 其中$L(u,\beta) = \frac{1}{2}(\beta - b(u))^T a^{-1}(u)(\beta - b(u))$，$a(x) = \sigma(x)\sigma^T(x)$。 **拟势**：从$x$到$y$的拟势定义为 $$V(x,y) = \inf_{T>0} \inf_{\phi \in AC^x_T} \{S^x_T(\phi): \phi(0)=x, \phi(T)=y\}$$ **大偏差原理**：在假设(H2)-(H3)下，系统满足Freidlin-Wentzell一致大偏差原理，包括： - (C) 紧性条件 - (L) 下界估计：$P\{\rho_T(X^{\varepsilon,x}_\cdot, \phi) < \delta\} \geq \exp\{-\frac{S^x_T(\phi)+\gamma}{\varepsilon^2}\}$ - (U) 上界估计：$P\{\rho_T(X^{\varepsilon,x}_\cdot, F^x_T(s)) \geq \delta\} \leq \exp\{-\frac{s-\gamma}{\varepsilon^2}\}$ ### 技术创新点 #### 1. 拟势的零值性质(核心突破) **命题4.2**：设$K$为无序链传递集，则存在吸引子$A = A(K)$满足$K \cap A = \emptyset$且 $$V(y,A) = 0, \quad \forall y \in K$$ **证明思路**： 1. 构造上对偶吸引子：$A = \Lambda \cap \bigcap_{t>0} \Phi_t[p,+\infty]]$，其中$p = \sup K$ 2. 定义$H^+(K) = \{x: y \ll \Phi_s(x) \text{ for some } y \in K, s \geq 0\}$ 3. 证明$H^+(K) \subset \text{Basin}(A)$ 4. 利用线性插值函数(LIF)构造能量接近零的路径 **关键洞察**：无序链传递集位于非单调流形上，从该集合到其对偶吸引子的转移几乎不需要能量消耗。 #### 2. 稀有事件概率估计(引理5.1) 对于无序链传递集或具有不稳定端点的静态p-arc $K$，存在吸引子$A$和常数$\kappa_2 > \kappa_1 > 0$使得： - (i) $\liminf_{\varepsilon \to 0} \varepsilon^2 \log P\{X^{\varepsilon,z}_T \in A_\eta\} \geq -\kappa_1$ - (ii) $\limsup_{\varepsilon \to 0} \varepsilon^2 \log P\{X^{\varepsilon,z}_T \notin A_\eta\} \leq -\kappa_2$ 这表明从不稳定集逃逸到吸引子的概率远大于从吸引子逃逸的概率。 #### 3. 主定理证明策略 通过不等式链： $$\mu_\varepsilon(B_\delta(y)) \leq I_1 \exp\{\frac{\kappa_1+\gamma}{\varepsilon^2}\}, \quad \mu_\varepsilon(A_\eta) \geq I_2 \exp\{\frac{\kappa_2-\gamma}{\varepsilon^2}\}$$ 其中$I_1 \leq I_2$(由不变测度性质)，得到 $$\mu_\varepsilon(B_\delta(y)) \leq \exp\{-\frac{(\kappa_2-\kappa_1)-2\gamma}{\varepsilon^2}\} \to 0$$ 从而$\mu(K) = 0$，排除无序集和不稳定p-arc。 ## 实验设置 ### 应用模型：生化控制回路 研究Griffith型正反馈系统： $$\begin{cases} \dot{x}_1 = f(x_r) - \alpha_1 x_1 \\ \dot{x}_j = x_{j-1} - \alpha_j x_j, \quad 2 \leq j \leq r \end{cases}$$ 其中$f(z) = \text{sgn}(z)\frac{|z|^m}{1+|z|^m}$（参数$m \geq 1$）。 ### 随机扰动模型 $$\begin{cases} dx_1 = (f(x_r) - \alpha_1 x_1)dt + \varepsilon\sigma_1(x_1)dW^1_t \\ dx_j = (x_{j-1} - \alpha_j x_j)dt + \varepsilon\sigma_j(x_j)dW^j_t \end{cases}$$ ### 验证条件 构造Lyapunov函数$V(x) = x^T Bx$（$B$为正定矩阵），验证： 1. 系统耗散性(H1) 2. 合作性和不可约性(或强单调性) 3. 条件(H2)和(H3) ## 实验结果 ### 主要结果：零噪声极限分类 **命题6.1**：在Griffith型假设下，零噪声极限$\mu$是有限个渐近稳定平衡点上Dirac测度的凸组合。 #### 情形1：$m = 1$（线性反馈） 定义$\varphi = \prod_{i=1}^r \alpha_i$，$h(z) = \frac{1}{1+z}$： - 若$\varphi \geq 1$：$\mu = \delta_O(\cdot)$（原点唯一稳定） - 若$0 < \varphi < 1$：$\mu = \lambda_1 \delta_{h^{-1}(\varphi)V_0}(\cdot) + \lambda_2 \delta_{-h^{-1}(\varphi)V_0}(\cdot)$ #### 情形2：$m > 1$（非线性反馈） 定义$\varphi_m = \max_{z>0} h(z) = \frac{m-1}{m} \cdot m^{-\frac{1}{m}}$： - 若$\varphi \geq \varphi_m$：$\mu = \delta_O(\cdot)$（原点全局吸引） - 若$0 < \varphi < \varphi_m$：存在$z_2 > 0$使得 $$\mu = \lambda_1 \delta_O(\cdot) + \lambda_2 \delta_{z_2 V_0}(\cdot) + \lambda_3 \delta_{-z_2 V_0}(\cdot)$$ ### 关键发现 1. **周期轨道的排除**：对于$m > 1$和$0 < \varphi < \varphi_m$，系统可能存在Hopf分岔产生的周期轨道(Selgrade, 1982)，但零噪声极限仍集中在平衡点上，不会集中在周期轨道上。 2. **参数依赖性**：零噪声极限的结构完全由参数$\varphi$和$m$决定，呈现明确的分岔结构。 3. **对称性**：由于$f$的奇对称性，非零平衡点总是成对出现($\pm z V_0$)。 ## 相关工作 ### 随机单调系统的轨道方法 - **Arnold & Chueshov (1998)**：随机单调系统的长期动力学简化 - **Flandoli et al. (2017)**：噪声同步现象 - **Chueshov & Scheutzow (2004)**：吸引子和不变测度结构 ### 基于分布的方法 - **Freidlin & Wentzell (1970, 2012)**：大偏差理论和拟势，等价类的随机稳定性 - **Huang et al. (2015-2018)**：Fokker-Planck方程的稳态，梯度系统的随机稳定性 - **Young (1986, 2002)**：双曲吸引子的随机稳定性，SRB测度作为零噪声极限 - **Kifer (1974, 1988)**：公理A系统的随机扰动 ### 单调动力系统理论 - **Hirsch (1984-1988)**：单调系统的收敛性理论 - **Smith (1995, 2017)**：合作系统专著 - **Mierczyński (1994)**：p-arcs理论 - **Poláčik (1989-1992)**：抛物方程的单调流 ### 本文优势 1. 首次处理可能具有无限等价类的系统 2. 不依赖预设动力学目标，基于系统大尺度结构 3. 直接利用连续时间大偏差理论(非离散化方法) 4. 提供正反馈系统的完整分类 ## 结论与讨论 ### 主要结论 1. **理论贡献**：建立了合作不可约系统的零噪声极限集中在Lyapunov稳定平衡点组成的简单有序集上，证明了动力学序的随机稳定性。 2. **技术突破**：证明了从非单调流形上无序链传递集到其对偶吸引子的拟势为零，这是克服计算高维拟势困难的关键。 3. **应用价值**：完整分类了Griffith型正反馈系统的零噪声极限，证明即使存在周期轨道，零噪声极限仍集中在有限个稳定平衡点上。 ### 局限性 1. **假设条件**： - 需要系统耗散性(H1) - 需要Lyapunov函数满足技术条件(H2)-(H3) - 不可约性(或强单调性)是必需的 2. **维数限制**：虽然理论适用于任意有限维，但Lyapunov函数的构造在高维可能困难。 3. **非合作系统**：方法不直接适用于非单调或竞争系统。 4. **权重系数未知**：虽然确定了零噪声极限的支撑集，但凸组合的系数$\lambda_i$未给出显式公式。 ### 未来方向 1. **系列后续工作**：本文是系列的第一篇，后续将研究： - 更一般的单调系统(非不可约情形) - 竞争系统的随机稳定性 - 时间依赖(非自治)单调系统 2. **定量分析**： - 计算凸组合系数$\lambda_i$ - 估计收敛速率 - 研究大偏差速率函数的精细结构 3. **应用拓展**： - 更复杂的生物化学网络 - 神经网络动力学 - 生态系统模型 4. **计算方法**：开发数值算法计算拟势和零噪声极限。 ## 深度评价 ### 优点 1. **理论深度**： - 首次系统研究可能具有无限等价类系统的随机稳定性 - 巧妙融合单调动力系统理论和大偏差理论 - 拟势零值性质的证明具有重要理论意义 2. **方法创新**： - 直接利用连续时间Freidlin-Wentzell理论，避免离散化复杂性 - 不依赖预设动力学目标，基于系统结构性质 - 构造性证明提供了清晰的几何图像 3. **数学严谨性**： - 假设条件明确且可验证 - 证明完整且逻辑清晰 - 引理和命题层次分明 4. **应用价值**： - 正反馈系统的完整分类具有实际意义 - 排除周期轨道的结果具有重要生物学含义 - 为生化网络的随机建模提供理论基础 ### 不足 1. **技术复杂性**： - 需要多个技术假设(H1)-(H3)的验证 - Lyapunov函数的构造在一般情况下可能困难 - 对非专家读者门槛较高 2. **定量信息有限**： - 未给出凸组合系数的显式表达 - 缺少收敛速率的估计 - 数值实验验证缺失 3. **适用范围**： - 限于合作不可约系统 - 耗散性假设排除了某些重要模型 - 非自治情形未涉及 4. **可计算性**： - 拟势的实际计算仍然困难 - 缺少数值算法和实现 ### 影响力 1. **学术贡献**： - 开创单调系统随机稳定性的系统研究 - 为后续研究提供框架和工具 - 连接动力系统和随机分析两大领域 2. **实用价值**： - 为生物化学系统、神经网络等提供理论指导 - 预测噪声环境下系统的长期行为 - 指导鲁棒性设计 3. **可复现性**： - 理论结果可验证 - 正反馈系统的例子可重现 - 需要补充数值实现代码 ### 适用场景 1. **生物系统**： - 基因调控网络(正/负反馈回路) - 细胞信号传导 - 种群动力学 2. **化学系统**： - 催化反应网络 - 振荡反应(如Belousov-Zhabotinsky反应) 3. **工程系统**： - 控制系统设计 - 神经网络动力学 - 经济系统建模 4. **理论研究**： - 单调系统的随机扰动 - 大偏差理论的应用 - 非平衡统计物理 ## 参考文献（关键文献） 1. **Freidlin & Wentzell (2012)**: Random Perturbations of Dynamical Systems - 大偏差理论基础 2. **Hirsch (1985, 1988)**: 单调动力系统的开创性工作 3. **Smith (1995)**: Monotone Dynamical Systems专著 4. **Benäım & Hirsch (1999)**: 随机逼近算法的合作系统 5. **Huang et al. (2015-2018)**: Fokker-Planck方程和随机稳定性系列工作 6. **Mierczyński (1994)**: p-arcs理论 7. **Selgrade (1979-1982)**: 正反馈系统的渐近行为和Hopf分岔 --- **总体评价**：这是一篇高质量的理论数学论文，在单调动力系统的随机稳定性研究中取得重要突破。通过巧妙结合动力系统理论和大偏差理论，解决了长期存在的技术难题。虽然技术性较强，但为理解噪声环境下复杂系统的长期行为提供了深刻洞察，具有重要的理论和应用价值。建议后续工作补充数值实验和更多实际应用案例。