3. Lyapunov指数的数学定义

要真正把Lyapunov指数用于混沌判定，必须先厘清它的数学定义。直观上，Lyapunov指数衡量的是相空间中两条无限接近的轨道沿某一方向平均分离或靠拢的速率。它把“初值微小扰动是否被放大”这一混沌的核心特征，转化为一个可计算、可比较的实数。本节从一维映射出发，逐步推广到高维系统的指数谱，并说明指数符号与轨道发散、收敛之间的对应关系。

一、一维映射的Lyapunov指数定义

考虑一维离散映射

x_{n+1} = f(x_n),  x_n ∈ ℝ

设两条轨道的初始点相差一个无穷小量 δx₀。经过一次迭代后，偏差近似为

δx_1 ≈ f'(x_0)·δx_0

连续迭代 n 次，偏差按链式法则累积为

δx_n ≈ |f'(x_0)·f'(x_1)···f'(x_{n-1})|·δx_0

若偏差平均每次迭代按指数规律变化，即 |δx_n| ≈ |δx_0|·e^{λn}，则取对数并令 n→∞，得到一维Lyapunov指数的标准定义：

λ = lim_{n→∞} (1/n) · Σ_{k=0}^{n-1} ln|f'(x_k)|

这个定义有两个要点。第一，它是对轨道上各点导数绝对值的对数取时间平均，而不是对某一点求值，因此 λ 是整条轨道的全局统计量。第二，绝对值保证了我们只关心偏差的放大或缩小幅度，方向翻转不影响分离速率。对一维映射，λ 只有一个值：λ>0 表示相邻轨道指数分离，系统对初值敏感；λ<0 表示轨道收敛到稳定周期或不动点；λ=0 对应临界情形，如倍周期分岔点。

以经典Logistic映射 f(x)=μx(1−x) 为例，其导数 f'(x)=μ(1−2x)。当 μ=4 时，理论Lyapunov指数为 ln2≈0.693，意味着偏差平均每次迭代放大约2倍。实际数值计算时，我们无法取无穷次迭代，通常取 N=10⁴～10⁶ 步，并舍弃前若干步暂态，用样本平均近似：

λ ≈ (1/N) · Σ ln|μ(1 − 2x_k)|

需要注意，若迭代过程中出现 f'(x_k)=0（如 x=0.5），对数为负无穷，实际编程中应跳过该点或加极小量保护，否则会污染平均值。

二、高维系统的Lyapunov指数谱

对 m 维连续系统 ẋ=F(x) 或 m 维映射 x_{n+1}=F(x_n)，扰动不再是一个标量，而是 m 维切空间中的向量。将初始扰动分解到 m 个相互正交的方向上，每个方向对应一个Lyapunov指数，按从大到小排列：

λ₁ ≥ λ₂ ≥ … ≥ λ_m

这组数称为Lyapunov指数谱。其严格定义借助切空间变分方程：设 δx 满足线性化方程 δẋ = J(x(t))·δx，J 为雅可比矩阵。对第 i 个方向，

λ_i = lim_{t→∞} (1/t) · ln( |δx_i(t)| / |δx_i(0)| )

其中 δx_i(t) 是沿第 i 个Lyapunov方向的扰动分量。由于数值上直接迭代会因所有方向迅速对齐到最大扩张方向而丢失信息，实际计算采用Gram-Schmidt正交化：每隔固定时间步对扰动向量组重新正交归一，记录每次的伸缩因子并累加对数，最终除以总时间得到各方向指数。这一流程正是Wolf算法等经典方法的核心。

指数谱提供了比单一指数丰富得多的动力学信息。对三维耗散系统，通常有 λ₁+λ₂+λ₃<0（相空间体积收缩）。若 λ₁>0、λ₂=0、λ₃<0，则对应奇怪吸引子上的混沌运动：一个方向扩张、一个方向沿流向中性、一个方向收缩。若 λ₁=0 且其余为负，则对应极限环。因此判断混沌的常用判据是：至少存在一个正Lyapunov指数，且系统有界（吸引子存在）。

三、指数与轨道发散/收敛的关系

Lyapunov指数的符号直接刻画了轨道沿对应方向的几何行为，这是理解其物理意义的关键。下图以三维系统为例，展示指数谱符号与轨道局部发散、收敛及整体吸引子结构之间的关系。

结合图示可以归纳出清晰的判读规则。当 λ>0 时，沿该方向的微小偏差以 e^{λt} 的速率放大，两条初始邻近的轨道会迅速分离，这正是“对初值敏感依赖”的定量表达；但分离不会无限进行，因为吸引子的有界性和收缩方向（λ₃<0）会把轨道拉回有限区域，于是形成“局部扩张、整体有界”的奇怪吸引子。当所有指数均为负时，任意扰动都被衰减，轨道收敛到稳定不动点或稳定周期。当最大指数为零而其余为负时，沿流向方向既不扩张也不收缩，对应极限环这类中性稳定结构。

在实际计算中，还有几点工程细节值得注意。其一，暂态舍弃：迭代初期轨道尚未进入吸引子，导数统计不代表吸引子上的平均行为，通常舍弃前 10³～10⁴ 步。其二，采样长度：λ 是极限值，有限样本只能给出估计，建议 N≥10⁵ 并观察 λ 随 N 的收敛趋势。其三，连续系统需先离散化或采用变分方程积分，步长过大会引入数值误差，使指数估计偏移。其四，对高维系统务必使用正交化流程，否则所有扰动方向会塌缩到最大扩张方向，导致 λ₂、λ₃ 估计失真。

理解定义之后，Lyapunov指数就不再只是一个抽象公式，而是一把把“轨道分离速率”量化为符号与数值的标尺：一维映射给出单一 λ，高维系统给出按大小排列的谱，而谱中正、零、负的组合，恰好对应发散、中性、收敛三种几何行为，共同构成混沌判定的数学基础。