2. 离散动力系统与连续动力系统

在计算Lyapunov指数之前，必须先明确所研究的系统属于哪一类动力系统。混沌现象既可以出现在离散迭代中，也可以出现在连续时间演化中，但两者的数学描述、数值实现和指数计算方式存在本质差异。本节围绕“映射”与“流”这两个基本概念展开，并介绍三个最常用的混沌系统示例，最后说明状态空间与轨迹的几何含义。

一、离散系统：映射与迭代

离散动力系统用差分方程描述，其一般形式为

x_{n+1} = F(x_n),  x_n ∈ R^m

其中 F 称为映射（map），n 为离散时间步。系统的演化不是连续进行的，而是通过一次次迭代“跳跃”产生。给定初值 x_0，反复作用 F 就得到轨道 {x_0, x_1, x_2, …}。离散系统的核心特征是：时间只取整数，状态更新完全由映射决定，没有微分方程中的时间导数。

最经典的例子是 Logistic 映射：

x_{n+1} = r · x_n · (1 - x_n),  x_n ∈ [0, 1]

参数 r 控制系统行为。当 r 从 0 增大时，系统依次经历不动点、倍周期分岔，最终在 r ≈ 3.57 之后进入混沌区。以 r = 4 为例，取初值 x_0 = 0.2，迭代可得 x_1 = 0.64、x_2 = 0.9216、x_3 ≈ 0.2890……序列看似随机，却完全由确定性规则生成。这种“确定性中的不可预测性”正是混沌的直观体现。Logistic 映射是一维映射，计算其 Lyapunov 指数只需沿轨道累加 ln|F'(x_n)|，实现非常直接。

二、连续系统：流与微分方程

连续动力系统用常微分方程描述：

dx/dt = f(x),  x ∈ R^m

其解 x(t) = φ_t(x_0) 称为流（flow），φ_t 表示从初值 x_0 出发、经过时间 t 后到达的状态。与映射的离散跳跃不同，流是随时间连续演化的，时间 t 可取任意实数。数值上通常用四阶 Runge–Kutta 等方法以固定步长 h 积分，把连续流离散化为近似迭代，但要注意步长会直接影响轨迹精度和指数估计的可靠性。

连续混沌系统的代表是 Lorenz 系统：

dx/dt = σ(y - x)
dy/dt = x(ρ - z) - y
dz/dt = xy - βz

取经典参数 σ = 10、ρ = 28、β = 8/3，系统呈现著名的“蝴蝶”双翼吸引子：轨迹在两个不稳定平衡点附近来回切换，永不重复也不发散。另一个常用例子是 Rössler 系统：

dx/dt = -y - z
dy/dt = x + a·y
dz/dt = b + z(x - c)

取 a = 0.2、b = 0.2、c = 5.7 时，轨迹呈现单卷的螺旋带状结构，其拓扑比 Lorenz 更简单，常被用作教学与算法验证对象。连续系统的 Lyapunov 指数需要同时演化切向量（变分方程）或采用相邻轨道重正化方法，比一维映射复杂，但思想仍是衡量相邻轨道的平均分离速率。

三、状态空间与轨迹

无论离散还是连续，系统的状态都可以看作状态空间（相空间）中的一个点。状态空间的维数等于独立状态变量的个数：Logistic 映射是一维的，状态空间就是区间 [0,1]；Lorenz 和 Rössler 系统是三维的，状态空间为 R³。系统随时间演化时，状态点在状态空间中划出的曲线称为轨迹（轨道）。

轨迹的几何形态直接反映系统性质。稳定不动点对应轨迹收敛到一点；周期运动对应闭合曲线；而混沌系统的轨迹被限制在一个有界区域内，却永不闭合、永不重复，这个区域称为奇异吸引子。Lorenz 吸引子的分形维数约为 2.06，介于二维曲面与三维体积之间，说明它既不是简单曲线也不是实体，而是具有自相似结构的复杂集合。

理解状态空间对计算 Lyapunov 指数至关重要：指数衡量的是状态空间中相邻轨迹沿不同方向的平均分离或收缩速率。三维连续系统有三个 Lyapunov 指数，按大小排列为 λ₁ ≥ λ₂ ≥ λ₃。混沌的判据是最大指数 λ₁ > 0；同时，耗散系统满足指数之和为负（对应相体积收缩）。例如 Lorenz 系统在经典参数下，典型结果约为 λ₁ ≈ 0.9、λ₂ ≈ 0、λ₃ ≈ -14.6，三者之和约等于 -(σ+1+β) ≈ -13.67，与散度一致。这一关系可作为数值计算的校验依据。

四、两类系统的实现差异

在代码层面，离散系统只需一个循环反复调用映射函数；连续系统则需要选择积分器和步长。以 Lorenz 为例，用步长 h = 0.01 的 RK4 积分，每步同时更新状态向量和切向量矩阵，每若干步做一次 Gram–Schmidt 重正化，避免切向量因指数增长而溢出。步长过大将引入积分误差，使指数估计偏离；步长过小则计算量上升。实践中可先用 h = 0.01 与 h = 0.005 各算一次，若结果差异在可接受范围内，说明步长已足够。

综上，映射与流是描述混沌的两套语言：前者以迭代跳跃刻画离散演化，后者以连续轨迹刻画时间流。Logistic、Lorenz、Rössler 分别代表一维离散、三维连续双翼与三维连续单卷三类典型对象。掌握状态空间与轨迹的几何图像，是后续正确计算并解读 Lyapunov 指数的前提。