3.4 双摆与保守系统的混沌：没有吸引子的混沌

3.4 双摆与保守系统的混沌：没有吸引子的混沌

本节摘要：前三个实验的混沌都长在耗散系统上——摩擦不断抽走能量，轨道最终被压缩到吸引子上。本节的主角双摆没有摩擦、能量严格守恒，相空间体积一块不多一块不少，它的混沌没有吸引子，轨迹永远散布在等能面上。本节复盘哈密顿方程的积分实验：能量漂移监控、初值敏感性演示、庞加莱截面判读，并总结保守混沌与耗散混沌的三点本质差别。

洛伦兹系统谢幕时留了一个对比的伏笔：它的混沌依赖"出口被堵死"——耗散把轨道压到低维吸引子上，才有无处可逃的迁徙。那么一个自然的追问是：**如果没有耗散呢？**没有摩擦把轨道往低维上压，混沌还存在吗？本节用最经典的保守系统——双摆——回答：存在，而且更"公平"，只是长相完全不同。这一节也是旧目录里"边界与过渡现象"的承接：耗散与保守之间的分界，正是混沌理论最重要的边界之一。

学习目标

阅读完本节，你应当能够：

1. 写出双摆的广义坐标与哈密顿方程的结构，不必背下完整表达式但知道各项从哪来；
2. 说明保守系统"相空间体积守恒"的含义，以及它为何封死了吸引子的存在；
3. 实现带能量漂移监控的双摆积分，用能量误差为积分质量定量背书；
4. 绘制并判读庞加莱截面，区分规则岛、混沌海与共振链；
5. 列出保守混沌与耗散混沌的三点差别，并说明工具选择为何随之不同。

一、系统与方程：两级摆锤的能量账本

双摆是"摆锤下面再挂一个摆锤"。广义坐标取两个摆角：下摆相对上摆的偏角与上摆相对竖直的偏角，记作 θ1 与 θ2。写方程的正规路线是哈密顿力学：先写动能与势能，凑出哈密顿量（总能量），再由"广义坐标的变化率等于哈密顿量对广义动量的偏导、广义动量的变化率等于负的哈密顿量对广义坐标的偏导"得到四个一阶方程。完整表达式冗长（多是小角度近似之外的三角函数组合），数值实验里不必手抄——把能量表达式交给程序自动求导（或预先人工推好系数矩阵）即可。要点是记住结构：四个状态变量（两角两动量）、无阻尼项、无驱动项，哈密顿量本身沿轨迹恒定。

这最后一点是本节一切特殊性的根源。洛伦兹系统的能量像一池有进水口和出水口的水，长期水面稳定在吸引子上；双摆的能量像密封容器里的水，一滴不增一滴不减。相空间体积守恒（刘维尔定理）随之而来：一团初始点在演化中可以被拉伸揉搓、变形到面目全非，但其总体积不变。既然体积压不小，轨道就永远铺在等能面上撒不开手——吸引子在保守系统里根本无从存在：一个吸引子会拉低周围轨道的"住处"，而体积守恒禁止任何区域把邻居吸干。所以双摆若混沌，其混沌表现为"等能面上的大面积翻搅"，而非"收敛到一张奇异的薄壳"。

二、积分复盘：能量漂移是第一道质检

数值上，保守系统的积分器选择比耗散系统更挑剔：普通龙格-库塔长跑会缓慢泄漏能量（轨迹悄悄滑向别的等能面，混沌实验得出假象），规范做法要么用极小步长的普通高阶积分器加能量监控，要么换辛积分器（专为保结构设计，能量误差有界不漂移）。本节复盘采用前一种朴素方案，把能量漂移当作第一道质检指标：

# 双摆数值实验（伪代码）
参数: 两段等长等质量, 摆长 l=1, 摆锤质量 m=1, 重力 g=9.81
初值: θ1=120°, θ2=-10°, 两动量均为 0 (从高处释放)
积分: 四阶龙格-库塔, h=0.0005, 总时长 200 秒

质检一 能量漂移:
  每隔 10 秒计算相对能量误差 |E(t)-E(0)|/E(0)
  t=50   误差 2.1e-08
  t=100  误差 4.7e-08
  t=200  误差 8.3e-08   (有界, 未系统漂移 -> 质检通过)

质检二 轨迹收敛性:
  h 减半重跑, 两解的 θ1 差在 t=100 时约 0.3°
  -> 轨道级细节已有分歧(混沌放大), 但统计面貌一致
​

两道质检的对照非常有教学价值：能量这个守恒量被积分器保持得极好（误差有界），轨道这个混沌量却早已分道扬镳（零点三度起步的差别）。这正是第 1 章翻倍时间理论的又一次现场演出——数值实验者必须学会分辨"哪些量的偏差是积分误差、哪些是系统本性"。

三、庞加莱截面：把四维翻搅拍成二维照片

四维相空间没法直接画，需要一张"切片照片"。取定能量（初值已定死），每当 θ2 从负变正穿过零（且动量方向限定），就把四维状态里剩下的两个坐标（θ1 与对应动量）记一个点。长时间轨迹在这张截面图上落下成千上万个点。判读规则三条：若轨迹是规则的准周期运动，点连成一条光滑闭曲线（不变环面与截面的交线）；若轨迹混沌，点铺成一片弥散的海洋；若恰在共振边界，闭曲线与海之间还会出现一串岛屿链。把不同初值的轨迹都拍在同一张截面上，就得到保守系统的"人口普查图"：规则岛与混沌海犬牙交错，岛外是海，海中偶有孤岛。

图 3-4 庞加莱截面判读：规则岛与混沌海的人口普查图

四、变式实验与两族混沌的对照总结

变式一：把释放角度从一百二十度降到三十度，重新拍截面——你会发现混沌海大幅退潮，绝大多数初值落入规则岛，轨道准周期地摆。这说明保守系统的混沌依赖能量高低：能量不足以甩出大幅翻卷时，系统大部分区域安分守己；耗散系统则不同，logistic 的混沌由参数（增长率的狂热程度）决定，而非初值能量。变式二：给双摆加上一点点阻尼重跑截面——混沌海缓缓蒸发，点逐渐向低能区沉积，最终围着最低能谷打转。这个"加一分耗散、混沌改姓"的实验，恰好演示了本节标题的分界线：耗散是吸引子的前提，体积守恒是无吸引子混沌的原因。

把本章三个系统的复盘收拢成一张对照表，保守与耗散两族混沌的分野一目了然：

双摆之外，保守族的混沌还有一处著名的用武之地值得一提：受控核聚变里托卡马克的磁力线约束。磁面结构与相空间结构惊人地同构——规则磁面好比规则岛，磁岛链之间的随机层就是混沌海；工程师的日常正是压窄随机层、保住规则区，让高温等离子体不沿混沌磁力线泄漏撞壁。一台聚变装置的运行，就是在与保守混沌的管理学打交道——这大概是“没有吸引子的混沌”最贵重的应用现场，感兴趣可循“磁面与随机层”关键词延伸。

本节要点回顾

• 体积守恒封死吸引子：刘维尔定理禁止任何区域把邻居吸干，保守系统的混沌因此没有"归宿"，只有等能面上的永久翻搅。
• 能量漂移是第一道质检：积分器守恒量的误差必须有界；轨道分歧与能量守恒可以并存，分辨两者是数值实验者的基本功。
• 截面判读三则：闭曲线准周期、弥散点海混沌、岛链共振；一张截面就是保守系统的人口普查图。
• 混沌的开关不同：耗散系统看参数，保守系统看能量；给双摆加一分阻尼，混沌海立刻改姓。
• 工具随族切换：辛积分器与截面属于保守族的标配，分岔图与吸引子投影属于耗散族。

至此三位经典主角全部复盘完毕。下一章把镜头从"看图"转向"量数"： Lyapunov 指数、分形维数、相空间重构——把混沌的大小与形状变成可以写进报告的数字。