导读：Duffing 方程 $\ddot x + 2\zeta\dot x + x + \alpha x^3 = F\cos\omega t$ 是非线性振动的果蝇：软/硬弹簧弯曲、跳跃现象、多解并存一应俱全。
渐近法与幅频曲线
谐波平衡（一阶谐波近似）：响应幅值 $A$ 满足 $\left[(1 + \frac{3\alpha A^2}{4} - \omega^2)^2 + (2\zeta\omega)^2\right]A^2 = F^2$。$ ilde\alpha>0$ 硬弹簧曲线右弯，$ ilde\alpha<0$ 软弹簧左弯——固有频率随幅值漂移是非线性的指纹。
右弯曲线在扫频中出现跳跃：升频与降频路径不同，振幅在拐点突跳——多解并存的直接观测证据。中间那支（负斜率）不稳定，实验不可达。
从幅频到混沌
激励加大后三次谐波、次谐波响应出现；再进入混沌吸引子（Poincaré 截面 + Lyapunov 指数判别，呼应本站混沌文章）。Duffing 的价值：一个方程串起非线性振动的全部现象谱系。
工程对应：隔振器橡胶大变形（软弹簧）、磁悬浮与结构大变形（硬弹簧）、船舶横摇（软硬切换）——幅频弯曲与跳跃在这些系统中有直接安全含义。
数值实践
打靶法/延续法追踪幅频曲线的全部解支（含不稳定支）；直接时域积分验证。数值线性代数中的延续法与矩阵计算（本站矩阵计算文章）在此会师。
常见误区
- 线性化后用固有频率外推大幅响应（幅值越大误差越大）。
- 只做升频扫频，漏掉跳跃与滞后路径。
吃透自检
- 能用谐波平衡推导幅频方程并画出三支曲线
- 能解释跳跃现象与不稳定支
- 能数值复现 Duffing 的周期与混沌响应