跳转至

第 9 章 弱非线性振荡器（Weakly-Nonlinear Oscillators）

作者

同 Ch 1-8，Thomas Witelski（Duke University，数学系）和 Mark Bowen（Waseda University，国际理工学中心）合著。本章是 Part II（Solution Techniques）的第 5 章——也是 Ch 6 末段预告的"长时间行为"主题的展开。

本章在全书中的定位：Ch 6 末段"渐近展开在 \(t = O(1/\sqrt{\varepsilon})\) 崩溃"——对振荡系统的"崩溃模式"系统化处理。当 \(\varepsilon\) 小时的长期累积效应（\(\sim \varepsilon t\)、\(\sim \varepsilon^2 t^2\)）会让解发生质变（频率偏移、振幅衰减/增长），简单正则展开无法捕捉。需要频率重正化（Poincaré-Lindstedt）或多时间尺度（MMTS）。

内容概述

本章回答一个核心问题：\(\varepsilon\) 小时的扰动如何在长时间内累积？ —— 通过两种方法： 1. §9.1-9.2 复习与失败：线性振荡器 + 周期驱动；共振（\(\gamma_k = \omega_0\)） 产生久期增长项（secular growth, \(t \cos \omega_0 t\)）。简单正则展开在 \(t = O(1/\varepsilon)\) 崩溃。 2. §9.3 Poincaré-Lindstedt 方法：把频率 \(\Omega(\varepsilon) = \omega_0 + \varepsilon \omega_1 + \ldots\) 也展开为 \(\varepsilon\) 的幂级数，让 \(\omega_1\) 选择消除久期项。有效时间 \(0 \le t \ll 1/\varepsilon^2\)（vs 正则展开的 \(1/\varepsilon\)）。 3. §9.4 多时间尺度法（MMTS, Multiple Scales）：引入慢时间 \(\tau = \varepsilon t\)，让 \(X(t, \tau)\) 是两个独立变量。慢时间 \(\tau\) 上的演化由振幅方程（amplitude equations）给出——耦合 ODE (9.32) 描述 \(A(\tau), B(\tau)\)。

核心思想：振荡系统的"小扰动长期累积"问题，核心是消除久期项。Poincaré-Lindstedt 通过重正化频率消除；MMTS 通过慢时间消除。两者数学等价（在首阶），但 MMTS 更灵活（能处理非周期解如衰减振荡）。

前置知识：Ch 6（扰动方法、主导平衡）、线性 ODE 求解（特征方程 \(r^2 + \omega_0^2 = 0\)）、三角恒等式（Fourier 分解）。

核心方程与概念

1. 线性振荡器复习（§9.1, eq. 9.1-9.4）

模型： $\(\ddot x + \omega_0^2 x = 0 \tag{9.1}\)$ $\(x(t) = A \sin(\omega_0 t) + B \cos(\omega_0 t)\)$

周期驱动（9.4）：\(\ddot x + \omega_0^2 x = C_k \sin(\gamma_k t) + D_k \cos(\gamma_k t)\) - 非共振（\(\gamma_k \neq \omega_0\)）：特解 \(= \frac{C_k}{\omega_0^2 - \gamma_k^2} \sin(\gamma_k t) + \frac{D_k}{\omega_0^2 - \gamma_k^2} \cos(\gamma_k t)\) - 共振（\(\gamma_k = \omega_0\)）：特解 \(= \frac{C_k}{2\omega_0} t \cos(\omega_0 t) + \frac{D_k}{2\omega_0} t \sin(\omega_0 t)\)（久期增长项，9.1 末尾）

久期项（secular terms）\(t \cos \omega_0 t\)、\(t \sin \omega_0 t\)：振幅线性增长，是真正的 ODE 精确解（在 \(\gamma_k = \omega_0\) 时）。但在渐近展开中这种"振幅增长"不是我们想要的——我们想要"\(\varepsilon\) 小时的扰动保持解的有界"。

2. 正则展开的失败（§9.2, eq. 9.5-9.12）

例子 1（9.5）：\(\ddot x + x = -\varepsilon x, x(0) = 1, \dot x(0) = 0\)。 - 精确解（9.10）：\(x = \cos(\sqrt{1 + \varepsilon}\, t)\)——频率偏移 - 正则展开（9.8）：\(x \sim \cos t - \frac{1}{2} \varepsilon t \sin t + \frac{1}{8} \varepsilon^2 (t \sin t - t^2 \cos t) + \ldots\)

例子 2（9.6）：\(\ddot x + x = -\varepsilon \dot x, x(0) = 1, \dot x(0) = -\frac{1}{2} \varepsilon\)。 - 精确解（9.11）：\(x = e^{-\varepsilon t/2} \cos(\sqrt{1 - \varepsilon^2/4}\, t)\)——振幅衰减 - 正则展开（9.9）：\(x \sim \cos t - \frac{1}{2} \varepsilon t \cos t + \frac{1}{8} \varepsilon^2 (t \sin t + t^2 \cos t) + \ldots\)

问题（9.12）：渐近阶数 \(\cos t \gg \varepsilon t \sin t\) 需要 \(\varepsilon t \ll 1\)，即 \(t \ll 1/\varepsilon\)。展开的有效时间窗口是 \(0 \le t \ll 1/\varepsilon\)——而精确解的物理意义（如频率偏移、振幅衰减）在 \(t = O(1)\) 已经显现，\(t = O(1/\varepsilon)\) 时已完全失效。

关键观察：虽然两个问题的精确解完全不同（一个频率偏移、一个振幅衰减），但正则展开在 \(t = O(1)\) 上非常相似——这是因为\(O(1)\) 时间内频率偏移和振幅衰减都是 \(O(\varepsilon)\) 量级，无法区分。需要更长时间才能分辨。

3. Poincaré-Lindstedt 方法（§9.3, eq. 9.13-9.20）

核心思想：把频率也展开 \(\Omega(\varepsilon) = \omega_0 + \varepsilon \omega_1 + \ldots\)（9.13-9.14），让 \(\omega_1\) 由消除久期项决定。

变量替换：\(x(t) = \tilde x(\theta), \theta = \Omega t\)，\(\tilde x(\theta) = \tilde x_0(\theta) + \varepsilon \tilde x_1(\theta) + \ldots\)

例子（9.15）：\(\ddot x + 4x = -\varepsilon x, x(0) = 3, \dot x(0) = -8\)。 - \(\omega_0 = 2\)（无扰动频率） - \(O(\varepsilon^0)\)（9.17a）：\(\tilde x_0 = -4 \sin \theta + 3 \cos \theta\) - \(O(\varepsilon)\)（9.17b, 9.19）：\(4 \tilde x_1'' + 4 \tilde x_1 = [16 \omega_1 - 4] \sin \theta + [-12 \omega_1 + 3] \cos \theta\) - 消除久期项（\(\sin \theta, \cos \theta\) 系数为零）：\(\omega_1 = 1/4\) - 最终解（9.20）： $\(x(t) \sim (-2 + \tfrac{1}{2}\varepsilon) \sin[(2 + \tfrac{1}{4}\varepsilon) t] + 3 \cos[(2 + \tfrac{1}{4}\varepsilon) t]\)$ - 有效时间：\(0 \le t \ll 1/\varepsilon^2\)（比正则展开长 \(\varepsilon\) 倍）

关键：Poincaré-Lindstedt 只能找周期解。\(\omega_1\) 唯一由"消除两个久期项"决定——但两个方程 vs 一个未知数，需要系统是简并（强制出现周期解）。对衰减解 (9.11)，没有 \(\omega_1\) 能消除久期项——需要 MMTS。

4. 多时间尺度法（MMTS, §9.4, eq. 9.21-9.32）

核心思想：引入慢时间 \(\tau = \varepsilon t\)，让 \(X(t, \tau)\) 是两个独立变量。快时间 \(t\) 描述"周期振荡"；慢时间 \(\tau\) 描述"振幅 / 相位的缓慢演化"。

链式法则（9.22-9.23）： $\(\frac{dx}{dt} = X_t + \varepsilon X_\tau, \quad \frac{d^2 x}{dt^2} = X_{tt} + 2\varepsilon X_{t\tau} + \varepsilon^2 X_{\tau\tau}\)$

ODE → PDE（9.24-9.25）：弱非线性振荡器 $\(X_{tt} + X = \varepsilon [f(X, X_t + \varepsilon X_\tau) - 2 X_{t\tau}] - \varepsilon^2 X_{\tau\tau}\)$

为什么把 ODE 变 PDE 反而更好？——因为 PDE 允许 \(X\) 随 \(\tau\) 变化，自由度多。ODE 只有 \(t\) 一个变量，解的形式被ODE 的局部结构严格限制；PDE 多了 \(\tau\)，让"振幅 \(A(\tau)\)"成为待定函数。

例子（9.27-9.32）：Duffing 振荡器 \(\ddot x + x = -\varepsilon \kappa \dot x + \varepsilon x^3\)。 - \(O(\varepsilon^0)\)：\(X_0(t, \tau) = A(\tau) \sin t + B(\tau) \cos t\)（系数是 \(\tau\) 的函数——这是 MMTS 与 Poincaré-Lindstedt 的关键区别） - \(O(\varepsilon)\)（9.31）：\(X_{1tt} + X_1 = \text{resonant terms (}\cos t, \sin t\text{)} + \text{non-resonant terms (}\cos 3t, \sin 3t\text{)}\) - 消除共振项给出振幅方程（9.32）： $\(\begin{aligned} -\kappa A + \frac{3}{4} A^2 B + \frac{3}{4} B^3 - 2 \frac{dA}{d\tau} &= 0 \\ \kappa B + \frac{3}{4} A^3 + \frac{3}{4} A B^2 + 2 \frac{dB}{d\tau} &= 0 \end{aligned}\)$ 这是关于 \(A(\tau), B(\tau)\) 的耦合 ODE——它们描述振幅和相位的慢演化。

最终解（9.34 极坐标形式）：\(X_0(t, \tau) = R(\tau) \sin(\omega_0 t + \Phi(\tau))\)，其中 \(R, \Phi\) 由极坐标振幅方程（9.32 转极坐标）描述。

5. 振幅方程的可解性条件（§9.5 末段 + §9.4 末尾）

作者明确指出（§9.5）：消除共振项 = Fredholm alternative 定理的应用。

Fredholm alternative（1900s）：线性方程 \(A x = f\) 有解 \(\iff\) \(f\) 与 \(A^*\) 的零空间正交（即 \(\langle f, v \rangle = 0\) 对所有 \(A^* v = 0\)）。对振荡方程，\(A = d^2/dt^2 + \omega_0^2\)（自伴），零空间是 \(\sin \omega_0 t, \cos \omega_0 t\)。所以 \(f\) 与这两个函数正交就是"消除共振项"。

这是一个一般原理：扰动展开中每阶的可解性条件 = Fredholm alternative。

6. 与 Ch 10 的衔接

Ch 9 是"多时间尺度"的弱非线性版本——振幅 \(A(\tau), B(\tau)\) 的演化是"慢"过程。 Ch 10 是"多时间尺度"的强非线性版本——van der Pol 振荡器、Michaelis-Menten 反应（\(\varepsilon\) 不在时间导数，而在非线性项系数）。 Ch 8 是"多空间尺度"——长波近似 + 边界层。

三章构成"多尺度方法"的完整图景： - Ch 8：空间多尺度（边界层） - Ch 9：时间多尺度（慢时间） - Ch 10：函数依赖性的多尺度（van der Pol、酶动力学）

关键结论

  1. 久期项 = 长期累积的"指示器"（§9.1-9.2）：当 \(t \cos \omega_0 t\) 这种项出现在渐近展开中，它告诉我们："\(\varepsilon\) 小时的扰动在 \(t = O(1/\varepsilon)\) 时间产生 \(O(1)\) 效应"。简单展开无法处理——需要重新参数化。
  2. Poincaré-Lindstedt = 频率重正化（§9.3）：把 \(\Omega = \omega_0 + \varepsilon \omega_1\) 看作"扰动修正的频率"，\(\omega_1\) 由"消除久期项"决定。有效时间延长到 \(1/\varepsilon^2\)——因为频率偏移是 \(O(\varepsilon)\)，振幅累积是 \(O(\varepsilon t)\)，但累积的累积是 \(O(\varepsilon^2 t^2)\)，需要 \(t \ll 1/\varepsilon^2\) 才能让二阶项可忽略。
  3. MMTS = 慢时间上的振幅演化（§9.4）：振幅 \(A(\tau), B(\tau)\) 是慢变量，由耦合 ODE (9.32) 描述。这等价于 Poincaré-Lindstedt 对周期解，但对衰减/增长解（如 (9.11)）也适用——振幅方程自动给出 \(A(\tau) \sim e^{-\kappa \tau/2}\) 形式的衰减。
  4. Poincaré-Lindstedt vs. MMTS 的对比：
  5. Poincaré-Lindstedt：只处理周期解；一个自由度 \(\omega_1\) 消除两个共振项需要简并系统
  6. MMTS：处理所有解（周期、衰减、增长）；两个自由度 \(A(\tau), B(\tau)\) 独立演化
  7. 振幅方程的物理意义（9.32）：
  8. \(-\kappa A - 2 dA/d\tau\)：线性阻尼 + 振幅演化
  9. \(+\frac{3}{4} A^2 B + \frac{3}{4} B^3\)：非线性立方刚度（Duffing 振子）
  10. 耦合 ODE 描述"非线性振子从初始振幅演化到稳态"
  11. Fredholm alternative 的一般性（§9.5 末段）：作者明说"the mathematical theory underpinning the solvability conditions is the Fredholm alternative theorem"。这是数学上严格的理由——为什么"消除久期项"是必须的：否则 ODE 无解。这种"必须满足的可解性条件"在 Ch 3 变分法、Ch 8 绝缘导线（8.25 末段）、Ch 10 快/慢系统中反复出现。
  12. 非线性效应 vs. 阻尼效应的耦合（9.32）：Duffing 振子 \(\ddot x + x = -\varepsilon \kappa \dot x + \varepsilon x^3\) 同时含阻尼（\(-\varepsilon \kappa \dot x\)）和非线性刚度（\(\varepsilon x^3\)）。振幅方程 (9.32) 给出两者相互作用：阻尼让 \(R\) 衰减，立方刚度让 \(R\) 增长——稳态 \(dR/d\tau = 0\) 给出极限环振幅 \(R_* = \sqrt{4\kappa/3}\)（van der Pol 极限环）。这是非线性动力学的"自维持振荡"。
  13. 与 Ch 10 van der Pol 的衔接（§9.3 末段 + §9.4 末段）：Ch 9 的振幅方程自动导出 van der Pol 极限环（当阻尼是非线性的 \(f(x) = (x^2 - 1) \dot x\)）。Ch 10 详细处理 van der Pol 振荡器。

挑战和开放性问题

  1. Poincaré-Lindstedt 的"二阶"展开（§9.3）：作者只做到 \(O(\varepsilon)\)（确定 \(\omega_1\)）。\(O(\varepsilon^2)\) 需要确定 \(\omega_2\)——但两个共振项只能用一个自由度 \(\omega_2\) 消除？ 这是 Poincaré-Lindstedt 的根本限制：简并性要求系统具有某种对称性。对非对称系统（"周期解不一定存在"），需要Lindstedt-Poincaré 的修正版或Newton's method**（迭代找周期解）。
  2. MMTS 的多尺度选择（9.21 末段）：作者固定 \(\tau = \varepsilon t\)。对更复杂情形，需要 \(\tau_2 = \varepsilon^2 t, \tau_3 = \varepsilon^3 t, \ldots\) 多套时间尺度——对 \(n\) 阶精度需要 \(n\) 套时间尺度。这与重正化群（renormalization group）的"对数修正"有相似性。
  3. 振幅方程的稳定性（9.32）：Duffing 振子 (9.27) 的振幅方程 (9.32) 是非线性耦合 ODE——它的固定点稳定性决定解的长期行为。\(R_* = \sqrt{4\kappa/3}\) 是稳定极限环（吸引子），但反相位 \(\theta_* = 0, \pi\) 等是鞍点。这种稳定性分析是非线性动力学的核心，作者没展开。
  4. 缺失：混沌与非线性振荡。当非线性是强（\(\varepsilon = O(1)\)）时，MMTS 失效——需要整个系统的几何分析（KAM 理论、Melnikov 方法）。本书不涉及混沌——这是 Part II 的一个重要限制。
  5. §9.5 末段提到的"averaging method"：作者说"the method of averaging [77, 92, 102] and near-identity transformations [58, 73] are related"。Averaging method（Krylov-Bogoliubov 1947）把 \(x(t)\) 在一个周期上平均得到"慢变量"——与 MMTS 数学等价但几何直观不同。近恒等变换（near-identity transformation）通过重新参数化消除久期项。
  6. §9.5 提到的"harmonic balance"（工程方法）：假设 \(x(t) = A \cos(\omega t) + B \sin(\omega t)\)，代入 ODE 平衡一次谐波系数。这是工程中常用的"近似方法"——但只有对弱非线性才精确。当非线性含 \(\sin(2t)\) 等高次谐波时，harmonic balance 失效。
  7. Poincaré-Lindstedt 与 Lindstedt-Poincaré：作者在 §9.3 用了"Poincaré-Lindstedt"，但实际工程文献常用"Lindstedt-Poincaré"。两种名称反映历史优先权争议（Poincaré 1892, Lindstedt 1882）。本书沿用"Poincaré-Lindstedt"——但读者应知道两种叫法都对。
  8. §9.5 末段的"Fredholm alternative"作者只在末段提到名字——严格定义在习题 9.13。这对 ODE 理论是核心——Riesz 投影、自伴算子谱理论、Sturm-Liouville 理论都是它的应用。

个人反思与批判性分析

本章是 Ch 6 扰动方法在长时间行为上的精致推广——把"\(\varepsilon\) 小时的扰动"和"\(\varepsilon\) 时间的累积"结合。几个值得反思的点：

  1. "久期项"是"渐近展开语言中的崩溃信号"。当 \(t \cos \omega_0 t\) 出现在展开中，它不是"误差"——它是物理：当 \(\varepsilon\) 小时的扰动以 \(\omega_0\) 频率共振累积时，解的振幅必须线性增长。这种"线性增长"在长时间是不真实的（物理振荡器的振幅不会无限增长），所以应该被新参数（频率偏移或阻尼）吸收。久期项 = 物理遗漏的信号。
  2. Poincaré-Lindstedt 与 MMTS 的"等价但不同角度"。Poincaré-Lindstedt 把 \(\omega_1\) 看作"频率的扰动"——\(\omega(\varepsilon) = \omega_0 + \varepsilon \omega_1 + \ldots\)，只有 \(\omega_0 + \varepsilon \omega_1\) 的精确度。MMTS 把 \(A(\tau), B(\tau)\) 看作"振幅和相位的演化"——任意长时间。对周期解两者等价（MMTS 的 \(A(\tau), B(\tau)\) 在稳态时为常数，等价于 Poincaré-Lindstedt 的固定频率）；对非周期解只有 MMTS 适用。
  3. "慢时间"\(\tau = \varepsilon t\) 是物理洞察。作者在 §9.4 开头用一个比喻：\(t\) 是"昼夜节律"（24 小时），\(\tau\) 是"年份"（累积变化）。这种"宏观-微观"分离是科学中普遍的模式——流体力学中的"湍流平均场 + 涨落"、量子力学中的"绝热演化 + Berry 相"、金融学中的"日内波动 + 长期趋势"——都是"快-慢"分离。
  4. "可解性条件"是 Ch 8-9-10 的共同主线。作者在 §8.3 末段（绝缘导线）、§9.4 末段（MMTS）、Ch 10（van der Pol）反复使用"消除久期项 = Fredholm alternative"思想。这是线性算子谱理论的应用——对自伴算子 \(A = d^2/dt^2 + \omega_0^2\)，零空间是 \(\sin, \cos\)，正交性就是"消除"。
  5. 作者没用群论语言。MMTS 实质上是"标度对称性"：\(\tau = \varepsilon t\) 是"小 \(\varepsilon\) 下的标度变换"——\(\tau\) 是 \(t\) 在 \(\varepsilon\) 标度下的"重正化时间"。这是 Ch 4 + Ch 5 + Ch 9 的统一——重正化群方法（Kadanoff, Wilson）把"短波长"和"长波长"用标度变换连接，MMTS 把"快时间"和"慢时间"用标度变换连接。本质都是"找一个标度群不变性，把多尺度问题约化为单尺度问题"。作者没明说——但 Ch 9 是群论视角在工程中的应用。
  6. Duffing 振荡器 (9.27) 的振幅方程 (9.32) 是"非线性科学"的入口。当 \(R\) 演化到固定点 \(R_*\)时，系统进入稳态振荡——这是Hopf 分岔的典型行为。当 \(R\) 继续演化（如变 \(R\) 的初值），系统可能"锁定"在某个相位——这是锁频（mode locking）。作者没展开——但Ch 1 末段"分岔"和 Ch 9 振幅方程是连贯的。
  7. 缺失主题：KAM 理论。当 \(\varepsilon\) 不是小量时，振幅方程 (9.32) 失效——系统进入准周期或混沌。KAM 理论（Kolmogorov 1954, Arnold 1963, Moser 1962）描述"小扰动下可积系统的准周期轨道大部分仍然存在"。本书不涉及——这是 Ch 9 的"高级"内容。

重要参考文献

  • [X1] Ali Hasan Nayfeh, Introduction to Perturbation Techniques, Wiley, 1981. (扰动技术的经典)
  • [X2] Ali Hasan Nayfeh and Dean T. Mook, Nonlinear Oscillations, Wiley, 1979. (非线性振荡)
  • [X3] Richard H. Rand, Lecture Notes on Nonlinear Vibrations, Cornell University, 2012. (非线性振动的免费讲义)
  • [X4] Jan A. Sanders and Ferdinand Verhulst, Averaging Methods in Nonlinear Dynamical Systems, Springer, 1985. (平均法)
  • [X5] Earl A. Coddington and Norman Levinson, Theory of Ordinary Differential Equations, McGraw-Hill, 1955. (ODE 理论经典)
  • [X6] Michael Tabor, Chaos and Integrability in Nonlinear Dynamics, Wiley, 1989. (混沌与可积性)
  • [X7] Steven H. Strogatz, Nonlinear Dynamics and Chaos, Westview, 1994. (非线性动力学入门)
  • [X8] Lawrence Perko, Differential Equations and Dynamical Systems, Springer, 2001. (ODE 与动力系统)
  • [X9] Carmen Chicone, Ordinary Differential Equations with Applications, Springer, 2006. (ODE 应用)
  • [X10] Jack K. Hale, Ordinary Differential Equations, Wiley, 1969. (ODE 的标准教材)