非线性动力学与混沌入门

许多常见的振动问题，最初往往假定回复力与位移成正比，得到线性谐振子的模型，其解为正弦或余弦函数。实际上，回复力只有在很小偏移时才能近似为线性；一旦角度变大、速度加快，非线性项的影响就变得显著，运动规律也随之改变。

单摆就是一个典型例子：小幅度时周期几乎恒定，振幅增大后周期变长。如果再考虑阻尼、外界驱动或多个系统的耦合，相平面上的轨迹就会从闭合曲线发展为螺旋、极限环，甚至表现出对初始条件极为敏感的长时演化。以下内容将依次介绍单摆的非线性行为、小角近似的误差分析、相图的描述方式、混沌与初值敏感性、参数分岔、以及分形和奇怪吸引子，穿插表格和算例，帮助将抽象概念与可计算公式建立联系。角位移用 $\theta$ 表示，时间导数记为 $\dot{\theta}$，重力加速度为 $g$，摆长为 $l$，均采用国际单位制。

单摆方程：从线性近似到大角度

质量集中于摆球、轻绳长 $l$ 的理想单摆，以竖直向下为参考，角位移 $\theta$ 的动力学方程为

$\ddot{\theta}+\frac{g}{l}\sin\theta=0$

小角度近似 $\sin\theta\approx\theta$（$\theta$ 以弧度计且足够小）把上式化为线性方程 $\ddot{\theta}+(g/l)\theta=0$，角频率 $$，周期 。一旦 不再很小， 与 的差会累积，周期随振幅变化，这就是非线性效应最直接的信号。

例1： 取 $l=1.00\ \mathrm{m}$，$g=9.80\ \mathrm{m}\cdot\mathrm{s}^{-2}$，则 $T_0=2\pi\sqrt{l/g}\approx2.01\ \mathrm{s}$。小角度实验里测得周期接近该值，说明线性近似有效；把最大摆角拉到 ，同样 、 下周期明显大于 ，说明非线性项已起作用。

例2： 无耗散时，系统机械能守恒。以最低点重力势能为零，摆球质量 $m$，则

$E=\frac{1}{2}m l^2\dot{\theta}^2+mgl(1-\cos\theta)$

给定 $E$ 可在相平面 $(\theta,\dot{\theta})$ 上画出等能曲线；能量越大，对应的最大摆角越大，周期越长，与能量积分形式一致。

小角近似：误差从哪一级开始抬头

小角公式 $\sin\theta\approx\theta$ 来自泰勒展开 $\sin\theta=\theta-\theta^3/6+\cdots$，首项误差量级约为 $|\theta|^3/6$（ 以弧度计）。同一摆长下，角振幅越大，非线性对周期的修正越明显；实验上先把振幅压到几度到十几度，测得的 与 接近，再逐步增大振幅，周期偏离即可被秒表分辨。

例9： $\theta=0.50\ \mathrm{rad}$ 时 $\sin\theta\approx0.479$，与 $\theta$ 本身相差约 $4.2\%$；回复力比例偏小，等效“弹簧变软”，周期比 $T_0$ 拉长，与上表趋势一致。课堂演示时把振幅从 增到 ，周期变化往往一眼可见，比死记“非线性”三字更有说服力。

相平面：位置与速度画在同一张图上

把广义坐标与广义速度取为横纵轴（单摆用 $\theta$ 与 $\omega=\dot{\theta}$），每一个瞬时状态对应相平面上的一个点，时间演化描成轨迹。无阻尼单摆能量守恒时，有界运动对应相图上的闭合曲线；小振动近似下，闭合曲线近似椭圆。线性中心平衡点附近轨迹一圈一圈绕开，属于稳定平衡点周围的典型图像。

例3： 仅看线性化 $\ddot{\theta}+\omega_0^2\theta=0$，令 $x=\theta$、$$，则 ，。相轨迹满足 ，积分得 ，即椭圆。把同一初始能量代回大角度方程数值积分，闭合曲线会略胖或略扁，但仍是圈，说明“绕圈”来自周期性，而“是不是正椭圆”才区分线性与否。

例4： 范德波尔振子等模型在适当参数下出现极限环：相平面上存在一条孤立的闭合轨道，附近轨迹随时间卷向该轨道。机械表里擒纵机构与能量补偿配合，宏观上也常抽象成“耗散与补给平衡”的周期吸引，与数学上的极限环概念同类，细节留给专门课程。

初值敏感与确定性混沌

方程完全确定、不含随机项，解仍可能对初值指数分离，长时间预测失效，这称为初值敏感或混沌现象。逻辑斯谛映射

$x_{n+1}=r x_n(1-x_n),\qquad 0\le r\le4,\ x_n\in[0,1]$

是离散时间的一维范例：$r$ 较小时迭代趋于固定点或周期点；$r$ 增大后出现倍周期分岔，进入混沌区后，两条极近的初值轨道会在有限步内明显分开。

例5： 取 $r=3.9$，初值 $x_0=0.300000$ 与 $x_0=0.300001$，用同一递推式各算 步，会发现后期数值差异远大于 。步数再多，两条序列几乎毫无关联，而方程本身没有写入任何随机数，敏感完全来自非线性迭代结构。

例6： 洛伦兹气象模型是三维连续动力系统的名例，参数合适时出现著名“蝴蝶形”吸引子。其意义在于：同一组方程，初值差一丁点，两条轨迹不久分道扬镳，长期天气式预测难度与初值精度绑在一起，“蝴蝶效应”是这种敏感的形象说法，不是文学夸张。

参数扫过临界点：从不动点到倍周期

映射 $x_{n+1}=f(x_n)$ 的不动点满足 $x^\ast=f(x^\ast)$。在逻辑斯谛映射里 ，平凡解 总存在；另一支在 时出现 。线性稳定性看导数 ：小于 时邻近点迭代被吸向该点，大于 则失稳。 从小变大，不动点失稳后常先出现周期 轨道（两点轮流跳），再出现周期 、、，称为倍周期序列；参数继续增大，系统进入混沌窗与周期窗交替的区域。全程不需要手算高次方程，数值实验画 横轴、长时间 纵轴的分岔图，就能把“参数一挪，吸引子改头换面”看在眼里。

a

例8： 取 $r=2.8$，从 $x_0=0.2$ 迭代足够多步，数值上会向 $x^\ast\approx0.643$ 靠拢，验证 $$ 的稳定吸引。改 ，同一初值经暂态后会看到两点交替，周期 占主导。两次实验同一递推式、只改 ，相空间形态完全不同，这就是分岔在实验桌面的缩影。

分形维数与奇怪吸引子：尺度里的自相似

通常曲线维数看作 $1$，曲面为 $2$。奇怪吸引子上点集的“有效维数”往往不是整数：放大后仍出现复杂层次，称为自相似或统计自相似。入门阶段不必死记多种维数定义，只需把握：用边长为 $\varepsilon$ 的小盒子去覆盖吸引子，所需盒子数 $N(\varepsilon)$ 随 $\varepsilon$ 变小的幂律若写 $N(\varepsilon)\propto\varepsilon^{-D}$，指数 可取非整数值，作为分形维的一种度量思路。

例7： 海岸线的长度随测量标尺变短而变长，标尺越细，弯折越多，总长大趋势上升，与光滑曲线不同；这类几何在科普与教材里常作为分形直觉入口，与物理里吸引子维数共享同一套“尺度律”语言。

练习题

选择题

1. 理想单摆大角度方程为 $\ddot{\theta}+(g/l)\sin\theta=0$，与小角度线性方程 $\ddot{\theta}+(g/l)\theta=0$ 相比，关于周期 与振幅关系的说法，正确的是（　　）

A. 大角度下 $T$ 恒等于 $2\pi\sqrt{l/g}$，与振幅无关

B. 大角度下 $T$ 随振幅增大而变长（有阻尼时另论）

C. 大角度下 $T$ 一定随振幅增大而变短

D. 非线性方程没有周期运动

2. 用相平面 $(\theta,\omega)$ 描述无阻尼单摆的周期运动，下列定性描述最合适的是（　　）

A. 轨迹沿直线远离原点，永不回头

B. 轨迹为伸向无穷的抛物线

C. 有界运动对应闭合相轨迹

D. 相轨迹必为严格正圆

3. 关于混沌的初值敏感性，下列说法正确的是（　　）

A. 初值敏感意味着方程里必须显含随机噪声项

B. 初值敏感意味着系统没有确定性方程

C. 确定性非线性系统也可能对初值指数分离

D. 初值敏感只出现在三维以上连续系统，一维映射不可能

4. 关于奇怪吸引子与分形维，下列说法最恰当的是（　　）

A. 奇怪吸引子一定是空间中的光滑曲面

B. 分形维数只能是正整数

C. 奇怪吸引子往往具有复杂几何，分形维可非整数

D. 有吸引子就一定混沌

计算题

5. 单摆摆长 $l=0.80\ \mathrm{m}$，取 $g=9.80\ \mathrm{m}\cdot\mathrm{s}^{-2}$。小角度近似下求角频率 $\omega_0$ 与周期 （单位分别用 与 ）。

6. 逻辑斯谛映射 $x_{n+1}=r x_n(1-x_n)$，取 $$，。求 、（保留三位有效数字）。

$T_0=\dfrac{2\pi}{\omega_0}=\dfrac{2\pi}{3.50}\approx1.80\ \mathrm{s}$。

（亦可用 $T_0=2\pi\sqrt{l/g}$ 直接得到同一结果。）

$x_2=r x_1(1-x_1)=2.5\times0.400\times(1-0.400)=2.5\times0.400\times0.600=0.600$。

故 $x_1=0.400$，$x_2=0.600$。