振动与混沌初步

非线性动力学是经典力学中最活跃的领域之一。从单摆的微小摆动到行星轨道的长期演化，周期运动、拟周期运动和混沌运动构成了力学系统行为的全貌。理解这些行为不仅需要精确计算，更需要对相空间结构有整体性的把握。本篇从周期运动的稳定性出发，逐步引入相图、庞加莱截面和李雅普诺夫指数等工具，最终通过具体模型展示混沌运动的本质特征。

周期运动的稳定性

单摆在小角度摆动时近似做简谐运动，但若推动幅度稍大，摆动周期就开始偏离简谐值；再大一些，摆甚至可能越过顶点做整圈转动。这种对初始条件的敏感性，正是稳定性研究的出发点。

考虑一个一维系统，运动方程为 $\ddot{q} = f(q)$，设系统有一个平衡点 $q_0$ 满足 $f(q_0) = 0$。令 $q = q_0 + \xi$，其中 $\xi$ 是小偏离量，代入运动方程并保留一阶项：

令 $\lambda = f'(q_0)$，则小偏离的演化完全由 $\lambda$ 的符号决定：

单摆有两个平衡点：最低点（$\theta = 0$）和最高点（$\theta = \pi$）。在最低点，$f'(0) = -g/L < 0$，是稳定平衡；在最高点，，是不稳定平衡。

对于周期轨道，稳定性分析需要研究轨道附近的小扰动是否会随时间增长。设系统存在一条周期为 $T$ 的轨道 $q_0(t)$，令 $q(t) = q_0(t) + \xi(t)$，线性化后得到变分方程：

由于 $q_0(t)$ 是周期的，方程中系数 $f'(q_0(t))$ 也以 $T$ 为周期，这类方程称为希尔方程。根据 Floquet 定理，其解具有形式 ，其中 是周期函数， 称为特征指数。若所有特征指数的实部均为零，则周期轨道稳定；若存在正实部，轨道不稳定。

相图的定性分析

将系统状态用广义坐标 $q$ 和广义速度 $\dot{q}$（或广义动量 $p$）在平面上表示，每条运动轨迹就成为相平面中的一条曲线，这种图形称为相图。相图是研究非线性系统最直观的工具，无需求解方程就能判断系统的全局行为。

单摆的哈密顿量为：

其中 $\theta$ 是摆角，$p = mL^2\dot{\theta}$ 是广义动量，$L$ 为摆长，$m$ 为摆球质量，$g\approx 9.8\mathrm{m}\mathrm{/}$。不同能量 对应的等能量曲线就是相图中的轨迹：

• 例一：单摆相图的定性特征

在相平面 $(\theta, p)$ 中，稳定平衡点 $(\theta = 0,\, p = 0)$ 被一族封闭曲线包围，对应有界的振荡运动。不稳定平衡点 $(\theta = \pm\pi,\, p = 0)$ 则是，相轨迹从远处趋近后又沿另一方向散开。将鞍点连接起来的曲线称为，它将相平面分为振荡区域（内部）和转动区域（外部）。分隔线上的运动需要无穷长时间才能到达鞍点，是一种理想极限情况。

保守系统（能量守恒）的相图中，轨迹就是等能量曲线。刘维尔定理保证相体积守恒，不同能量的轨迹不会相交。

对于阻尼系统，相图发生根本变化。以线性阻尼振子 $\ddot{q} + 2\gamma\dot{q} + \omega_0^2 q = 0$ 为例，相轨迹不再是封闭曲线，而是向内螺旋趋向平衡点，因为阻尼不断耗散能量，平衡点从中心型变为（欠阻尼情形）或（过阻尼情形）。

庞加莱截面

对于高维系统，相图难以可视化。庞加莱截面（Poincaré section）通过"采样"的方式将高维问题降维：选取一个低维截面，记录每次轨迹穿越该截面的交点，用这些离散点的分布来反映运动的整体结构。

对于二维哈密顿系统，相空间是四维的 $(q_1, p_1, q_2, p_2)$，能量守恒将运动限制在三维超曲面上。固定 ，只记录 时的穿越点 ，就得到了二维庞加莱截面。

• 例二：庞加莱截面的三种典型模式

庞加莱截面最大的优势在于将连续时间的流"压缩"成离散映射：轨迹在截面上的演化可以用一个庞加莱映射 $\mathbf{x}_{n+1} = P(\mathbf{x}_n)$ 来描述。周期轨道对应映射的不动点，拟周期运动对应映射的不变曲线，混沌运动则对应映射上的随机游走行为。

李雅普诺夫指数

两条初始条件相差 $\delta_0$ 的轨迹，经过时间 $t$ 后间距变为 $\delta(t)$。若轨迹稳定，间距不会无限增长；若运动是混沌的，间距会以指数速率增长：

这里 $\lambda$ 称为李雅普诺夫指数（单位：$\mathrm{s^{-1}}$），是混沌程度的量化判据。严格定义为：

对 $n$ 自由度哈密顿系统，相空间是 $2n$ 维的，共有 $2n$ 个李雅普诺夫指数，从大到小排列为 $\lambda_1 \geq \lambda_2 \geq \cdots \geq \lambda_{2n}$。由于哈密顿系统的时间可逆性，指数成对出现：，故最大指数 与最小指数互为相反数。只要最大李雅普诺夫指数 ，系统就是混沌的。

• 例三：双摆的混沌

双摆由两个串联的单摆组成，是最典型的混沌系统之一。在能量较小时（小摆角），运动接近线性叠加，$\lambda_1$ 接近零，轨迹可预测。当能量超过某阈值，双摆开始做不规则的大幅运动，$\lambda_1$ 变为正值（量级约 $1\,\mathrm{s^{-1}}$），仅几秒后，两条初始位置相差毫米级的轨迹就会完全分离，无法依靠精确计算来做长期预报。

Hénon–Heiles 哈密顿量

Hénon 和 Heiles 于 1964 年提出了一个简单的二维哈密顿模型，用于研究星系势中恒星运动的混沌性质。哈密顿量为：

前两项是各向同性谐振子，后两项是非线性耦合势。势能 $V(x,y) = (x^2 + y^2)/2 + x^2 y - y^3/3$ 具有三重旋转对称性，但不能用变量分离法求解解析解。

以能量 $E$ 为控制参数，系统呈现出从有序到混沌的典型转变：

这个模型的重要意义在于，一个仅含两个自由度的简单确定性系统，随参数变化即可从完全有序过渡到几乎完全混沌。在庞加莱截面上，随 $E$ 增大，原本光滑的闭合曲线逐渐断裂、卷曲，最终分解为随机散点——这一过程直观地展示了混沌的产生机制。

逻辑斯蒂方程与倍周期分岔

逻辑斯蒂方程是理解混沌最简单的离散模型，描述种群在有限资源下的增长：

其中 $x_n \in [0,1]$ 代表第 $n$ 代的种群密度（相对于环境承载量），$r \in (0,4]$ 是增长率参数。

随 $r$ 从 0 增大到 4，系统长期行为发生一系列显著变化：

相邻两次分岔点之间的参数区间之比，趋向一个普适常数，称为费根鲍姆常数：

• 例四：$r = 3.5$ 时的逻辑斯蒂迭代

取 $r = 3.5$，初始值 $x_0 = 0.5$，迭代后序列很快进入 4 周期循环：$\{0.3828,\, 0.8269,\, 0.5009,\, 0.8750,\, 0.3828,\ldots\}$。若改取 （仅差 ），前几十步几乎相同，但约第 50 步后两序列开始出现明显差异，体现了临近分岔时系统对初始值的敏感性已经开始增强。

分形维数

混沌吸引子（耗散系统中混沌运动在相空间留下的轨迹集合）往往具有非整数维数，称为分形。描述分形最直观的工具是盒计数维数（box-counting dimension）：

用边长为 $\varepsilon$ 的小盒子覆盖分形集合，设需要 $N(\varepsilon)$ 个盒子，则盒计数维数定义为：

对光滑曲线（一维对象），$N(\varepsilon) \propto 1/\varepsilon$，故 $d_{\mathrm{f}} = 1$；对平面区域（二维对象），$N(\varepsilon) \propto 1/\varepsilon^2$，故 。分形集合的特征是 的增长介于两者之间， 取非整数值。

• 例五：康托集的分形维数

从区间 $[0,1]$ 出发，每步去掉剩余每段的中间三分之一，无限迭代后的极限集称为康托集。第 $n$ 步剩余 $2^n$ 段，每段长度为 $3^{-n}$，用 $\varepsilon = 3^{-n}$ 大小的盒子覆盖，需要 个：

康托集处处不连续却包含不可数无穷个点，维数在 0 和 1 之间，是最简单的分形之一。

练习题

• 1. 关于相图中稳定平衡点和不稳定平衡点的特征，下列说法正确的是：

A. 稳定平衡点在相图中表现为鞍点，不稳定平衡点表现为中心

B. 稳定平衡点在相图中被封闭轨迹包围（中心型），不稳定平衡点处有两条轨迹趋近、两条轨迹离去（鞍点型）

C. 两类平衡点在相图中的表现完全相同，只能通过求解运动方程来区分

D. 保守系统中不存在稳定平衡点，所有平衡点均为鞍点

• 2. 逻辑斯蒂方程 $x_{n+1} = r\,x_n(1-x_n)$ 中，当参数 从 1 增大到 4 时，系统经历倍周期分岔进入混沌。关于这一过程，下列说法正确的是：

A. 第一次倍周期分岔发生在 $r = 1$，此时系统从灭绝转为稳定平衡

B. 倍周期分岔是指周期从 $T$ 变为 $T/2$，即频率加倍

C. 相邻分岔点之间的间距之比趋于费根鲍姆常数 $\delta \approx 4.669$，这一常数对一大类映射具有普适性

D. 进入混沌区域后，系统不再存在任何周期轨道

• 3. 李雅普诺夫指数 $\lambda$ 描述相邻轨迹的指数分离速率。关于李雅普诺夫指数，下列说法正确的是：

A. 混沌运动要求系统所有的李雅普诺夫指数均大于零

B. 哈密顿系统中，只要最大李雅普诺夫指数 $\lambda_1 > 0$，系统就是混沌的

C. 李雅普诺夫指数越大，表示系统越稳定

D. 周期运动的最大李雅普诺夫指数一定等于 $-1$

• 4. 关于分形维数，下列说法正确的是：

A. 分形维数必须是整数，不能取非整数值

B. 用边长 $\varepsilon$ 的盒子覆盖分形集合，所需盒子数 $N(\varepsilon)$ 满足 $N(\varepsilon) \propto \varepsilon^{d_{\mathrm{f}}}$

C. 康托集的分形维数为 $\ln 2/\ln 3 \approx 0.631$，介于 0 维点和 1 维曲线之间

D. 混沌吸引子的分形维数总等于其所在相空间的整数维数

• 5. 逻辑斯蒂方程 $x_{n+1} = r\,x_n(1-x_n)$ 在 时存在稳定不动点 。（一）由 求出非零不动点的表达式；（二）令 ，对不动点处线性化，推导 满足的递推关系，并给出稳定性条件；（三）解释为何在 处不动点失稳并发生倍周期分岔。

• 6. 用盒计数法计算科赫曲线的分形维数。科赫曲线的构造规则是：将线段三等分，用两段等边三角形斜边替换中间三分之一（每段因此变为 4 段，长度缩为原来的 $1/3$）。（一）经过 $n$ 步后，共有多少段？每段长度是多少？（二）用边长 $\varepsilon = (1/3)^n$ 的盒子覆盖，需要多少个？（三）由此求出科赫曲线的分形维数。

$x^* = 0$ 是灭绝状态（$r > 1$ 时不稳定）；$x^* = (r-1)/r$ 是非零不动点，在 $1 < r < 3$ 时稳定。

（二）令 $x_n = x^* + \xi_n$，代入 $x_{n+1} = r\,x_n(1-x_n)$，保留一阶项：

利用 $x^* = r\,x^*(1-x^*)$ 消去常数项，得：

将 $x^* = (r-1)/r$ 代入：$r(1 - 2x^*) = r\bigl(1 - 2(r-1)/r\bigr) = 2 - r$，故：

稳定性条件为 $|2 - r| < 1$，即 $1 < r < 3$，与预期一致。

（三）在 $r = 3$ 时，线性化乘子恰为 $|2 - 3| = 1$，不动点处于稳定与不稳定的临界点。当 $r$ 稍超过 3，乘子 $|2-r| > 1$，小偏离不再收缩，不动点变为不稳定。取而代之的是 2 周期轨道：$x^{(1)} \to x^{(2)} \to x^{(1)} \to \cdots$，这正是倍周期分岔的发生机制——一个稳定不动点分裂为一对稳定两点循环。

（三）由盒计数维数公式：

科赫曲线的分形维数约为 1.262，介于 1 维曲线和 2 维面之间，反映了其"比曲线更弯折、但尚未填满平面"的几何特征。科赫雪花由三条科赫曲线围成，其周长无穷大，但面积有限，是分形几何中最经典的例子之一。