2025-11-10T02:50:53.145942

Two approaches to average stochastic perturbations of integrable systems

Kuksin

We discuss two approaches to study the long-time behaviour and infinite-time behaviour of solutions for integrable hamiltonian systems under small stochastic perturbations. Then we compare these results with those for deterministic perturbations of integrable systems.

academic

Two approaches to average stochastic perturbations of integrable systems

基本信息

论文ID: 2509.17159
标题: Two approaches to average stochastic perturbations of integrable systems
作者: Sergei Kuksin
分类: math.DS (Dynamical Systems)
发表时间: 2025年10月11日 (arXiv v2)
论文链接: https://arxiv.org/abs/2509.17159

摘要

本文讨论了研究可积哈密顿系统在小随机扰动下解的长时间行为和无穷时间行为的两种方法，并将这些结果与可积系统确定性扰动的结果进行了比较。

研究背景与动机

问题背景

经典问题: 研究可积哈密顿系统的小扰动在大时间尺度上的行为是动力系统的经典问题，可追溯到18世纪末
传统方法: 现有研究主要关注两种情况：
- A) 扰动是小的光滑向量场
- B) 扰动是小的哈密顿场
新的挑战: 本文考虑第三种情况——随机扰动，这在理论和应用上都具有重要意义

研究动机

随机扰动比确定性扰动更能反映实际物理系统中的噪声影响
现有的平均化理论主要针对确定性系统，随机情况下的理论相对不完善
需要发展新的数学工具来处理随机扰动下的长时间动力学行为

核心贡献

提出了两种不同的方法来研究可积系统随机扰动的平均化问题
建立了快慢系统方法，通过作用-角坐标将问题转化为慢变量的研究
发展了有效方程方法，直接在原坐标系中构造平均化方程
证明了收敛性定理，给出了随机扰动下作用变量的极限行为
提供了与确定性扰动结果的详细比较，揭示了随机性的稳定化作用

方法详解

基本设置

考虑辛空间 $(R^{2n}_{x,y}, dx \wedge dy)$ 中的Birkhoff可积哈密顿系统，其哈密顿函数只依赖于作用变量 $I_j = (x_j^2 + y_j^2)/2$ 。引入复坐标 $v_j = x_j + iy_j$ ，系统可写为：

$\frac{\partial}{\partial t}v_k(t) = i\nabla_k H(I)v_k, \quad k = 1, \ldots, n$

ε-小随机扰动

考虑如下随机扰动系统（慢时间 $\tau = εt$ ）：

$\dot{v}_k(\tau) = iε^{-1}\nabla_k H(I)v_k + P_k(v) + \sum_{j=1}^{n_1} B_{kj}(v)\dot{\beta}^c_j(\tau)$

其中 $\{\beta^c_k(t)\}$ 是标准独立复Wiener过程。

基本假设

(A1) Anosov非退化条件: 对几乎所有 $I \in R^n$ ，频率向量 $\nabla H(I)$ 的分量在整数上线性无关 (A2) 局部Lipschitz条件和多项式增长条件 (A3) 方程适定性和矩的一致有界性 (A4) 噪声非退化：矩阵 $B(v)$ 的秩等于 $n$

第一种方法：快慢系统

作用-角坐标变换

引入作用-角坐标 $(I, \phi)$ ：

$I_j = \frac{1}{2}|v_j|^2$
$\phi_j = \arg v_j$

慢变量方程

作用变量满足： $\dot{I}^ε(\tau) = F(v^ε) + G(v^ε)\dot{\beta}(\tau)$

这是一个快慢系统，其中作用 $I^ε_k$ 是慢变量，角度 $\phi^ε_k$ 是快变量。

平均化方程

通过在快变量上平均得到： $\dot{I}(\tau) = \langle\langle F \rangle\rangle(I) + \langle\langle G \rangle\rangle(I)\dot{\beta}(\tau)$

其中平均化定义为： $\langle\langle F \rangle\rangle(I) = \frac{1}{(2\pi)^n} \int_{T^n} F(I, \phi) d\phi$

主要结果

定理 2.1: 解 $v^ε(\tau)$ 的作用曲线 $I^ε(\tau)$ 在分布意义下预紧，每个极限过程都是平均化方程的弱解。

局限性: 当矩阵 $B(v)$ 不是常数时，平均化系数可能只有 Hölder-1/2 光滑性，导致解的唯一性问题。

第二种方法：有效方程

直接平均化

受Krylov-Bogolyubov平均化启发，直接在复坐标中构造有效方程：

$\dot{a}(\tau) = \langle\langle P \rangle\rangle(a) + \langle\langle B \rangle\rangle(a)\dot{\beta}^c(\tau)$

其中平均化通过旋转算子定义： $\langle\langle P \rangle\rangle(a) = \frac{1}{(2\pi)^n} \int_{T^n} \Phi_\omega \circ P(\Phi_{-\omega}a) d\omega$

主要优势

避免了作用-角坐标的奇异性
保持原始方程的光滑性
总是有唯一解

核心定理

定理 3.1: 在假设A1)-A4)下，有效方程有唯一解，当 $ε \to 0$ 时， $I^ε(\tau)$ 在分布意义下收敛到 $a(\tau)$ 的作用。

定理 3.3: 如果另外假设A5)（长时间估计），则收敛性在时间上一致成立。

定理 3.4: 如果 $P(v) = P_1(v) + P_2(v)$ ，其中 $P_2$ 是哈密顿场，则有效方程可修正为只平均化非哈密顿部分。

实验设置与应用

应用场景

5.1 振子链: 考虑非线性振子链： $\ddot{q}_k = -Q(q_k), \quad k = 1, \ldots, n$ 其中势函数 $Q$ 光滑、凸且奇，在零点附近 $Q(q) = αq + βq^3 + o(|q|^4)$ 。

5.2 阻尼/驱动哈密顿系统: 考虑方程： $\dot{v}_k(\tau) = iε^{-1}\nabla_k H(I(v))v_k + 2i\frac{\partial}{\partial \bar{v}_k}h(v) - ν_k v_k + b_k \dot{\beta}^c_k(\tau)$