3.3 周期轨道与极限环：自持振荡从哪里来

3.3 周期轨道与极限环：自持振荡从哪里来

本节摘要：定义先行——极限环是相空间中"孤立"的闭合轨道：附近的轨道从里从外都向它收拢（或从它逃开），振荡的振幅由方程本身决定。它与线性中心的根本区别在于有无"能量账本"：中心的振荡不结账、一碰就散，极限环每个周期自动打平收支（小振幅注入能量、大振幅消耗能量），因此是结构稳定的自持振荡。本节以范德波尔方程为主角，把账本、判法定理与数值实验一次配齐。

本节在图鉴中位置关键：它是不动点与混沌之间的中间层——分岔理论里霍普夫分岔的产物，倍周期瀑布的起点。也是"周期轨道与极限环"这组对照的主场：同样是来回振荡，线性与非线性的身份天差地别。

谁在给钟表上发条

线性无阻尼振子（2.1 节的"中心"）也振荡，但它像一部不用发条的行星钟：振幅完全由你拨多大决定，拨多少转多少圈，永远转下去。这份"永动"是幻觉——刚度或阻尼的一丝偏差就让它衰减或发散，且任何真实器件都有损耗，没人能给行星钟上发条。工程与自然里的振荡器——钟摆钟、心脏、电子哨、神经元——全都不是中心：它们在摩擦与损耗的世界里维持固定振幅，秘密是每个周期都有一本打平的能量账。

范德波尔方程把这本账写成了三项式：

• x'' − μ(1 − x²)·x' + x = 0，μ>0

读法：等效阻尼是 −μ(1−x²)，振幅小于一为负（注入能量，振荡壮大），大于一为正（消耗能量，振荡收敛）。于是振幅被夹逼到一个中间值附近：从任意初值出发，最终都落在同一条闭合轨道上——极限环。1920 年代范德波尔在真空管电路里发现它，1927 年就用它模拟过心脏；今天神经元模型、激光速率方程、气动颤振分析里都是同一个骨架。

线性中心对非线性极限环

图：线性中心与非线性的稳定极限环

判定定理：什么时候该找极限环

二维连续系统有一条著名的"存在担保"——庞加莱-本迪克松定理：若一条轨道被困在一个不发散的紧区域内，且区域内平衡点只有有限个（比如一个不稳定焦点），则轨道的极限集必含极限环。直觉版本：出不去、又无处可停，只好转圈。判"没有"也有工具（杜拉克判据、梯度系统无环），但工程上更常用的正是定理的直觉版：证一个平衡点不稳定、再画一个囚笼（常由零倾线围出），环的存在即板上钉钉。

展开案例：范德波尔的数值体检

背景：验证"初值不同、终态同圈"，并量一量极限环的振幅与周期。

操作：取 μ=1，从四个天差地别的初值出发积分，末段测量振幅与周期：

import numpy as np
from scipy.integrate import solve_ivp

mu = 1.0
def vdp(_t, s):
    x, v = s
    return [v, mu * (1 - x * x) * v - x]

for x0, v0 in [(0.1, 0.0), (3.0, 0.0), (0.0, 4.0), (-2.5, 1.0)]:
    sol = solve_ivp(vdp, [0, 80], [x0, v0], rtol=1e-9, max_step=0.05,
                    dense_output=True)
    ts = np.linspace(60, 80, 4000)             # 丢弃暂态，只看末段
    x = sol.sol(ts)[0]
    amp = 0.5 * (x.max() - x.min())
    peaks = ts[1:][np.diff((x[1:] > 1.8).astype(int)) == 1]
    T = np.mean(np.diff(peaks)) if len(peaks) > 1 else float("nan")
    print(f"初值({x0},{v0}) -> 振幅 {amp:.3f} 周期 {T:.3f}")
​

结果（典型输出）：四个初值的末段振幅都收敛到约 2.01，周期约 6.6——初值只改变了到达路径，没有改变终点。

解读：振幅 2.01 恰在方程"账本"预测的位置：x>1 后阻尼转正，动能开始归还。变式一：把 μ 调到 5 重跑，波形从近正弦变成"缓慢充电、瞬间放电"的弛豫振荡，周期拉长到约 8（大 μ 时周期≈1.61μ≈8.1）——神经元发放与滴水的龙头都属于这个形态。变式二：加强迫项 A·sin(ωt)，扫 A 会遇到锁频：振荡频率被拖到与强迫成有理数比——这就是第六章心律失常案例里"起搏器拖带"的数学内核。

💡 关键直觉：判断真实系统里的振荡是"共振"还是"自持"，做一个减法实验即可——撤掉外部激励。共振的振荡随激励消失（线性响应），自持的振荡毫发无损（极限环是系统自己的属性）。桥的颤振属于后者，所以危险得多。

⚠️ 常见坑：极限环只在二维及以上的连续流中存在（一维流单调，2.2 节的老规矩）；而且"有振荡波形"不等于"有极限环"——强迫系统里可能只是被拖着走的受迫响应，要靠撤激励实验或庞加莱截面（5.3 节）甄别。

问题：实验里怎么确认"这是极限环"而不是别的振荡？

按证据强度递进做四件事。其一，测振幅稳定性：撤掉扰动后振幅回到同一水平（初值无关）是极限环的第一签名，对照线性共振"激励多大振幅多大"。其二，做撤激励实验：外部输入归零后振荡依旧，排除受迫响应。其三，画相轨迹或李萨如图：稳定极限环在相平面上收拢为固定的闭合曲线，线宽由噪声决定。其四，扫参数看振幅如何出生：若振幅从零按平方根曲线连续长出，即霍普夫分岔的标准姿态（3.2 节）。四条证据里前两条在现场十分钟可做，足以把"自持"与"受迫"分开——这一步分清，后续控制策略（断供能还是改频率）才不会开错药方。

本节要点回顾

• 极限环 = 孤立闭合轨道，振幅由方程决定；线性中心闭合轨道成族、振幅由初值决定。
• 能量账本是自持的全部秘密：范德波尔的负阻尼注入与正阻尼消耗在每个周期打平，账平则环稳。
• 庞加莱-本迪克松定理：囚笼加有限平衡点即可担保极限环存在，工程排查按这个顺序走。
• 弛豫振荡与锁频：大 μ 的充放电形态与外加激励的频率拖带，分别在神经元与心律场景反复出现。

振荡有了，下一步更狂野：轨道永不闭合、又永不发散的混沌，见下节。