3. 周期轨道的数学刻画

在上一节中，我们从几何直观上认识了“轨道会绕回来”这一现象。本节要把这种直观转化为可计算、可判定的数学语言。核心问题是：给定一个动力系统，如何用方程和条件精确地说明“某条解是周期的”？这个刻画一旦建立，后续的数值延拓、稳定性分析和分岔追踪才有统一的出发点。

周期解的定义与基本记号

考虑一般的一阶常微分方程组

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

若存在常数 T > 0，使得某个解 x(t) 满足对一切 t 都有

x(t + T) = x(t)

则称 x(t) 为周期解，T 为一个周期。使上式成立的最小正数 T 称为基本周期（最小周期）。注意定义中“对一切 t 成立”是关键：只在一个时刻相等不构成周期性，必须是整条解曲线在时间平移下不变。

从相空间看，周期解对应一条闭合曲线。若 x(t) 是周期解，则它自身在相空间中画出一条闭轨；反过来，落在闭轨上的解必然是周期解。这一点在自治系统中尤其干净，后文会专门讨论。

周期与频率：一对互为倒数的参数

周期 T 描述“绕一圈要多久”，频率则描述“单位时间绕几圈”。二者关系为

f = 1 / T,    ω = 2π / T

其中 f 是普通频率（赫兹），ω 是角频率（弧度/秒）。在数值延拓中，我们通常把 T 作为未知量直接求解，因为周期本身出现在边界条件 x(0) = x(T) 中；而 ω 更多用于频域分析和与实验谱线对照。举例来说，若某振荡器的周期测得 T = 0.628 s，则 ω ≈ 10 rad/s，f ≈ 1.59 Hz。这三个量只是同一事实的不同表述，选择哪个取决于后续分析的需要。

需要提醒的是，周期解的基本周期未必等于系统中某个参数的“自然周期”。非线性系统里，周期往往依赖于振幅和参数，这正是延拓要追踪的对象。

自治系统与非自治系统

若 f 不显含 t，即 dx/dt = f(x)，系统称为自治系统；若 f 显含 t，则称为非自治系统。这一区分对周期轨道的刻画影响很大。

在自治系统中，解具有时间平移不变性：若 x(t) 是解，则 x(t + t₀) 也是解。因此周期解对应的闭轨是一条“轨道”，其上任意一点都可以作为起点，周期相同。自治系统的闭轨不会自相交（否则违反解的唯一性），这给数值计算带来便利：我们只需在一条横截面上寻找回归点。

非自治系统则不同。若 f(x, t) 关于 t 以 T₀ 为周期，即 f(x, t + T₀) = f(x, t)，系统可能产生与激励周期相关的周期响应，也可能产生周期为 T₀ 整数倍甚至有理数倍的次谐波响应。此时“周期解”的周期不一定是 T₀，需要单独判定。非自治系统的周期解在相空间中不一定闭合，因为相空间没有把时间作为坐标；若把时间也纳入，扩展相空间中的轨迹才是闭合的。这一点在数值实现中必须明确，否则容易把“看起来回来了”误判为周期。

周期轨道与闭轨的关系

二者的关系可以概括为：周期解是时间域的概念，闭轨是相空间的概念，在自治系统中二者一一对应。具体来说：

• 自治系统中，每条闭轨都对应一簇周期解，它们仅相差时间平移；基本周期相同。
• 闭轨上不能有平衡点，否则解会在有限时间停住，无法继续绕行。
• 闭轨是孤立曲线还是成族出现，取决于系统结构：保守系统（如无阻尼谐振子）的闭轨连续成族，周期随能量变化；耗散系统（如范德波尔振子）通常只有孤立闭轨，即极限环。
• 非自治系统中，相空间轨迹可以不闭合，但扩展相空间（x, t）中的轨迹闭合，周期判定应在此空间进行。

这一区分直接指导数值策略：对自治系统，可在相空间中取庞加莱截面，把周期轨道转化为截面上不动点问题；对非自治系统，则常以激励周期为采样间隔构造频闪映射。

两个标准例子：谐振子与范德波尔振子

无阻尼谐振子是最简单的周期系统：

dx/dt = v
dv/dt = -ω₀² x

其解为 x(t) = A cos(ω₀t + φ)，周期 T = 2π/ω₀，与振幅 A 无关。相空间中是同心椭圆族，每条都是闭轨，但都不是孤立的——任意小的扰动仍落在另一条闭轨上。这类系统没有极限环，周期由参数 ω₀ 唯一决定，数值延拓时若以 ω₀ 为参数，周期随之线性变化，可作为验证代码正确性的基准算例。

范德波尔振子则展示了非线性与耗散的共同作用：

dx/dt = v
dv/dt = μ(1 - x²)v - x

当 μ > 0 时，原点是不稳定焦点，而系统存在一条孤立的极限环。小振幅扰动被推向外侧，大振幅扰动被拉向内侧，最终都收敛到同一条闭轨。其周期不再由某个线性参数直接给出，而随 μ 变化：μ 较小时 T ≈ 2π（接近谐振子），μ 增大后周期单调增大，μ 很大时进入张弛振荡区，T 约按 μ 增长。这条极限环正是后续延拓章节的主要研究对象——我们需要追踪它在参数 μ 变化下如何移动、何时失稳。