微扰理论:从原理到波形分析

微扰理论的核心思想是将一个无法精确求解的非线性方程,拆解为一系列可以求解的线性方程。 我们将通过一个经典的非线性振荡器模型(Duffing Oscillator 的简化版)来演示这一过程:
ẍ + x + εx³ = 0

第 1 步:问题设定与假设

我们面对的是一个包含微小非线性项 εx³ 的微分方程。我们假设解 x(t) 可以写成关于 ε 的幂级数。

方程: ẍ + x = -ε
假设解: x(t) = x₀ + εx₁ + ε²x₂ + ...

目标:找到 x₀ (零阶解), x₁ (一阶修正), x₂ (二阶修正)。

第 2 步:代入方程

将级数代入原微分方程。这一步看起来很复杂,但我们只需要关注幂次。

(ẍ₀ + εẍ₁...) + (x₀ + εx₁...) = -ε(x₀ + εx₁...)³

利用二项式展开右边的立方项:
(x₀ + εx₁)³ ≈ x₀³ + 3x₀²(εx₁) + ...

第 3 步:分离零阶项 (Order ε⁰)

我们只看那些不包含 ε 的项。这给出了我们熟悉的“无扰动”方程。

ε⁰: ẍ₀ + x₀ = 0 这是标准的简谐振动方程。解为 x₀(t) = A cos(t)。

意义: 这是系统的基准行为(如同弹簧没有非线性时的运动)。

第 4 步:分离一阶项 (Order ε¹)

找出所有系数为 ε¹ 的项。注意右边的源项出现了。

ε¹: ẍ₁ + x₁ = -x₀³ 这变成了一个“受迫振动”方程!驱动力是零阶解的立方。

我们将上一步求出的 x₀ 代入右边,即可解出 x₁。这里会产生 cos(3t) 项。

第 5 步:分离二阶项 (Order ε²)

这一步通常在入门课程中略过,但原理相同:找出系数为 ε² 的项。

ε²: ẍ₂ + x₂ = -3x₀²x₁ 驱动力混合了零阶解 x₀ 和一阶修正 x₁。

总结: 我们不仅仅是在解一个难题,而是在依次解一串简单的题。后一个方程利用前一个方程的解作为输入。

1 / 5

静态波形对比分析

拖动滑块,观察红色微扰波的振幅变化及其对总波形的影响。

0.30
基础解 x₀ (cos t)
微扰项波形 (εx₁)
差异区域
近似总解 x
看图指南:
  • 蓝色虚线:理想的简谐波。
  • 红色实线:这就是“微扰”本身!注意它的频率是基础波的3倍 ($3\omega$)。它展示了我们具体往系统里加了什么。
  • 紫色粗线:蓝线与红线叠加后的结果。
  • 当 ε 为正时,红色波在基波峰值处是凹下去的,导致总解变“平”。