一般波动方程（以经典的机械波/电磁波标量形式为例）的推导，核心思路是：从物理定律出发，得到关于介质中某点位移（或场量）对空间和时间的二阶偏微分方程。下面给出最常见的两种推导路径，并说明其一般形式。

一、一维绳波（机械波）的推导

这是最直观的推导，适合理解“波动方程”的物理来源。

1. 模型假设

  • 一根均匀柔软的细绳，线密度为 \(\mu\)，张力为 \(T\)（假设张力大小恒定，方向沿绳切线）。
  • 绳只在横向做微小振动，位移记为 \(y(x,t)\)。
  • 忽略重力、阻尼等。

2. 取微元受力分析

取绳上 \(x\) 到 \(x+\Delta x\) 的一小段微元。

在两端，张力方向沿切线，与水平方向夹角分别为 \(\theta_1\)、\(\theta_2\)。
微小振动下： \[ \sin\theta \approx \tan\theta = \frac{\partial y}{\partial x} \]

横向合力为： \[ F_y = T\sin\theta_2 - T\sin\theta_1 \approx T\left(\left.\frac{\partial y}{\partial x}\right|_{x+\Delta x} -\left.\frac{\partial y}{\partial x}\right|_x\right) \]

由泰勒展开： \[ \left.\frac{\partial y}{\partial x}\right|_{x+\Delta x} -\left.\frac{\partial y}{\partial x}\right|_x \approx \frac{\partial^2 y}{\partial x^2}\Delta x \]

所以： \[ F_y \approx T\frac{\partial^2 y}{\partial x^2}\Delta x \]

3. 应用牛顿第二定律

微元质量 \(m=\mu\Delta x\)，横向加速度为 \(\frac{\partial^2 y}{\partial t^2}\)，因此： \[ T\frac{\partial^2 y}{\partial x^2}\Delta x = \mu\Delta x \frac{\partial^2 y}{\partial t^2} \]

消去 \(\Delta x\)： \[ \frac{\partial^2 y}{\partial x^2} = \frac{\mu}{T}\frac{\partial^2 y}{\partial t^2} \]

令 \[ v=\sqrt{\frac{T}{\mu}} \] 则得到一维波动方程： \[ \boxed{\frac{\partial^2 y}{\partial x^2} = \frac{1}{v^2}\frac{\partial^2 y}{\partial t^2}} \]

二、一般三维波动方程的形式

把上述思想推广到三维连续介质，设标量场为 \(u(\mathbf r,t)\)，则一般齐次波动方程为： \[ \boxed{\nabla^2 u - \frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}=0} \] 或写成： \[ \boxed{\frac{\partial^2 u}{\partial t^2}=v^2\nabla^2 u} \]

其中：

  • \(\nabla^2\) 是拉普拉斯算子；
  • \(v\) 是波速，由介质性质决定。

在直角坐标系中： \[ \nabla^2 u= \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2}+ \frac{\partial^2 u}{\partial z^2} \]

所以三维波动方程展开为： \[ \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2}+ \frac{\partial^2 u}{\partial z^2}

\frac{1}{v^2}\frac{\partial^2 u}{\partial t^2} \]

三、电磁波波动方程的推导（简述）

对于真空中的电磁波，从麦克斯韦方程组出发。

1. 真空麦克斯韦方程组

\[ \nabla\cdot \mathbf E=0 \] \[ \nabla\cdot \mathbf B=0 \] \[ \nabla\times \mathbf E=-\frac{\partial \mathbf B}{\partial t} \] \[ \nabla\times \mathbf B=\mu_0\varepsilon_0\frac{\partial \mathbf E}{\partial t} \]

2. 对法拉第定律取旋度

\[ \nabla\times(\nabla\times \mathbf E) =-\frac{\partial}{\partial t}(\nabla\times \mathbf B) \]

利用矢量恒等式： \[ \nabla\times(\nabla\times \mathbf E) =\nabla(\nabla\cdot \mathbf E)-\nabla^2\mathbf E \]

在真空中 \(\nabla\cdot \mathbf E=0\)，所以： \[ \nabla\times(\nabla\times \mathbf E)=-\nabla^2\mathbf E \]

又由安培定律： \[ \nabla\times \mathbf B =\mu_0\varepsilon_0\frac{\partial \mathbf E}{\partial t} \]

因此： \[ -\nabla^2\mathbf E =-\mu_0\varepsilon_0\frac{\partial^2 \mathbf E}{\partial t^2} \]

即： \[ \boxed{\nabla^2\mathbf E =\mu_0\varepsilon_0\frac{\partial^2\mathbf E}{\partial t^2}} \]

同理可得： \[ \boxed{\nabla^2\mathbf B =\mu_0\varepsilon_0\frac{\partial^2\mathbf B}{\partial t^2}} \]

令 \[ c=\frac{1}{\sqrt{\mu_0\varepsilon_0}} \] 则： \[ \boxed{\nabla^2\mathbf E =\frac{1}{c^2}\frac{\partial^2\mathbf E}{\partial t^2}} \] \[ \boxed{\nabla^2\mathbf B =\frac{1}{c^2}\frac{\partial^2\mathbf B}{\partial t^2}} \]

这就是电磁波波动方程，\(c\) 为光速。

四、一般波动方程的共性

无论是绳波、声波、电磁波，还是量子力学中的波动方程，推导通常都遵循以下步骤：

  1. 确定研究对象：位移、压强、电场、磁场等。
  2. 取微元或局部区域：分析其受力、通量或场量变化。
  3. 应用基本物理定律：
    • 牛顿第二定律；
    • 胡克定律；
    • 麦克斯韦方程组；
    • 连续性方程等。
  4. 做线性近似：小振动、小扰动、均匀介质等。
  5. 消去中间量，得到关于场量 \(u\) 的二阶偏微分方程。
  6. 整理成标准形式： \[ \frac{\partial^2 u}{\partial t^2}=v^2\nabla^2 u \] 或 \[ \nabla^2 u-\frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}=0 \]

五、补充：非齐次波动方程

如果存在外力或源项 \(f(\mathbf r,t)\)，则方程变为： \[ \boxed{\nabla^2 u-\frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}=f(\mathbf r,t)} \] 或 \[ \boxed{\frac{\partial^2 u}{\partial t^2}-v^2\nabla^2 u=F(\mathbf r,t)} \]

例如电磁波有电荷、电流源时，就是非齐次波动方程。

总结一句话：
一般波动方程的推导，本质上是把介质或场的局部相互作用，通过牛顿定律、胡克定律或麦克斯韦方程组等写成微分关系，再消去中间变量，最终得到“时间二阶导 = 波速平方 × 空间二阶导”的标准形式。

非微小振动呢？sin不约等于tan试试 如果去掉“微小振动”近似，绳波方程会变成非线性波动方程。核心变化是：不能再把 \(\sin\theta\approx\tan\theta\approx y_x\) 当作线性近似，而必须保留完整的几何非线性。

一、重新建立绳的几何关系

设绳的位移为 \(y(x,t)\)，绳上一点处切线与水平方向夹角为 \(\theta\)，则：

\[ \tan\theta=\frac{\partial y}{\partial x}=y_x \]

而 \[ \sin\theta=\frac{\tan\theta}{\sqrt{1+\tan^2\theta}} =\frac{y_x}{\sqrt{1+y_x^2}} \]

\[ \cos\theta=\frac{1}{\sqrt{1+y_x^2}} \]

这就是不约等于的关键：sin 和 tan 之间差一个 \(\sqrt{1+y_x^2}\) 因子。

二、微元受力分析（不近似）

仍取 \(x\) 到 \(x+\Delta x\) 的一小段绳微元。

1. 张力方向

张力 \(T\) 沿绳的切线方向。设张力大小在绳中可能不是常数，记为 \(T(x,t)\)。

在右端 \(x+\Delta x\)，张力矢量为： \[ \mathbf T_2 = T_2(\cos\theta_2,\sin\theta_2) = T_2\left(\frac{1}{\sqrt{1+y_x^2}},\frac{y_x}{\sqrt{1+y_x^2}}\right)_{x+\Delta x} \]

在左端 \(x\)，张力矢量为： \[ \mathbf T_1 = -T_1(\cos\theta_1,\sin\theta_1) = -T_1\left(\frac{1}{\sqrt{1+y_x^2}},\frac{y_x}{\sqrt{1+y_x^2}}\right)_{x} \]

2. 横向合力

\[ F_y = T_2\frac{y_x}{\sqrt{1+y_x^2}}\bigg|_{x+\Delta x} -T_1\frac{y_x}{\sqrt{1+y_x^2}}\bigg|_{x} \]

如果绳是理想柔软绳且无外力沿绳方向，水平方向合力为零，可以推出张力水平分量沿绳为常数： \[ T\cos\theta = T_0 = \text{常数} \]

于是： \[ T=\frac{T_0}{\cos\theta}=T_0\sqrt{1+y_x^2} \]

代入横向力： \[ T\sin\theta = T_0\sqrt{1+y_x^2}\cdot \frac{y_x}{\sqrt{1+y_x^2}} = T_0 y_x \]

所以横向合力变为： \[ F_y = T_0 y_x|_{x+\Delta x} - T_0 y_x|_x = T_0\left(y_x(x+\Delta x)-y_x(x)\right) \]

泰勒展开： \[ F_y \approx T_0 y_{xx}\Delta x \]

注意：这里横向合力形式上仍然像线性结果，但这只是因为用了“水平张力分量守恒”这一条件。真正的非线性来自另外两个地方：

  1. 弧长与线密度变化：微元实际长度不是 \(\Delta x\)，而是 \(\Delta s=\sqrt{1+y_x^2},\Delta x\)。
  2. 纵向运动：非微小振动时，绳元会沿绳方向运动，不能忽略纵向位移和纵向加速度。

三、只考虑横向运动但保留弧长非线性

如果仍假设绳元只做横向运动，但弧长变化不能忽略，则微元质量为： \[ m=\mu \Delta s=\mu\sqrt{1+y_x^2},\Delta x \]

牛顿第二定律： \[ T_0 y_{xx}\Delta x = \mu\sqrt{1+y_x^2},\Delta x \cdot y_{tt} \]

消去 \(\Delta x\)： \[ \boxed{ T_0 y_{xx} = \mu\sqrt{1+y_x^2}, y_{tt} } \]

整理： \[ \boxed{ y_{tt} = \frac{T_0}{\mu}\frac{y_{xx}}{\sqrt{1+y_x^2}} } \]

这就是一个非线性波动方程。
当 \(|y_x|\ll 1\) 时，\(\sqrt{1+y_x^2}\approx 1\)，退化为经典线性波动方程： \[ y_{tt}=v^2 y_{xx},\qquad v=\sqrt{\frac{T_0}{\mu}} \]

四、更严格：同时考虑纵向位移

非微小振动时，绳上质点不仅有横向位移 \(y\)，还有纵向位移 \(u\)。设质点位置为： \[ \mathbf r(x,t)=\bigl(x+u(x,t),, y(x,t)\bigr) \]

则弧长微元： \[ ds=\sqrt{(1+u_x)^2+y_x^2},dx \]

单位切向量： \[ \hat{\mathbf t} =\frac{(1+u_x,,y_x)}{\sqrt{(1+u_x)^2+y_x^2}} \]

张力矢量为： \[ \mathbf T=T\hat{\mathbf t} \]

对微元应用牛顿第二定律，分别得到纵向和横向方程：

纵向方程

\[ \mu, u_{tt}

\frac{\partial}{\partial x} \left[ T\frac{1+u_x}{\sqrt{(1+u_x)^2+y_x^2}} \right] \]

横向方程

\[ \mu, y_{tt}

\frac{\partial}{\partial x} \left[ T\frac{y_x}{\sqrt{(1+u_x)^2+y_x^2}} \right] \]

再结合绳的本构关系（张力与应变关系）。若绳近似满足胡克定律： \[ T=T_0+EA\left(\frac{ds}{dx}-1\right) \] 其中 \(E\) 为弹性模量，\(A\) 为截面积，则： \[ T=T_0+EA\left(\sqrt{(1+u_x)^2+y_x^2}-1\right) \]

把 \(T\) 代入上面两个方程，就得到完整的非线性弹性绳方程组。这是一组耦合的非线性偏微分方程，通常没有简单解析解，只能用数值方法或摄动法求解。

五、如果只保留 sin 与 tan 的差异，而不考虑弧长变化

有些人可能会想：只把 \[ \sin\theta\approx y_x \] 换成 \[ \sin\theta=\frac{y_x}{\sqrt{1+y_x^2}} \] 会怎样？

若同时忽略弧长变化、忽略纵向运动，并假设张力大小 \(T\) 为常数，则横向合力为： \[ F_y

T\left( \frac{y_x}{\sqrt{1+y_x^2}}\bigg|_{x+\Delta x}

\frac{y_x}{\sqrt{1+y_x^2}}\bigg|_{x} \right) \]

泰勒展开： \[ F_y \approx T\frac{\partial}{\partial x} \left( \frac{y_x}{\sqrt{1+y_x^2}} \right)\Delta x \]

计算导数： \[ \frac{\partial}{\partial x} \left( \frac{y_x}{\sqrt{1+y_x^2}} \right)

\frac{y_{xx}}{(1+y_x^2)^{3/2}} \]

因此： \[ F_y \approx T\frac{y_{xx}}{(1+y_x^2)^{3/2}}\Delta x \]

质量仍用 \(\mu\Delta x\)，则： \[ \mu\Delta x, y_{tt}

T\frac{y_{xx}}{(1+y_x^2)^{3/2}}\Delta x \]

得到： \[ \boxed{ y_{tt}

\frac{T}{\mu}\frac{y_{xx}}{(1+y_x^2)^{3/2}} } \]

这个方程描述的是几何非线性弦振动，其中非线性完全来自 \(\sin\theta\) 与 \(\tan\theta\) 的差异。

当 \(|y_x|\ll 1\) 时： \[ (1+y_x^2)^{3/2}\approx 1 \] 又回到线性波动方程。

六、两种非线性方程的区别

处理方式 方程 非线性来源
保留弧长变化，假设张力水平分量常数 \(y_{tt}=\frac{T_0}{\mu}\frac{y_{xx}}{\sqrt{1+y_x^2}}\) 质量随弧长变化
保留 sin/tan 差异，假设张力大小常数 \(y_{tt}=\frac{T}{\mu}\frac{y_{xx}}{(1+y_x^2)^{3/2}}\) 张力方向投影非线性
完整纵向+横向耦合 耦合方程组 几何非线性 + 材料非线性

七、物理意义

非微小振动下，会出现：

  1. 波形畸变：不同振幅部分波速不同，大振幅处传播速度改变。
  2. 高次谐波产生：正弦波会逐渐产生二倍频、三倍频等谐波。
  3. 孤子/激波现象：在某些非线性介质中可出现孤立波或陡化。
  4. 不能直接用叠加原理：两个波相遇后不再简单线性叠加。

一句话总结：
去掉 \(\sin\theta\approx\tan\theta\) 后，绳波方程从线性波动方程变成非线性波动方程，典型形式如
\[ y_{tt}=\frac{T_0}{\mu}\frac{y_{xx}}{\sqrt{1+y_x^2}} \quad\text{或}\quad y_{tt}=\frac{T}{\mu}\frac{y_{xx}}{(1+y_x^2)^{3/2}} \] 具体形式取决于你保留哪些物理效应。

把它展开看看前3项 好的，我们来把非线性项做小振幅但非微小的展开，看看前几项是什么。这样能清楚看出：线性波动方程之外，会出现哪些非线性修正项。

一、先看较简单的那个方程

前面得到一个常见形式：

\[ y_{tt}=\frac{T_0}{\mu}\frac{y_{xx}}{\sqrt{1+y_x^2}} \]

令

\[ v^2=\frac{T_0}{\mu} \]

则

\[ y_{tt}=v^2\frac{y_{xx}}{\sqrt{1+y_x^2}} \]

二、展开 \(\frac{1}{\sqrt{1+y_x^2}}\)

利用二项式展开：

\[ (1+z)^{-1/2} =1-\frac12 z+\frac{3}{8}z^2-\frac{5}{16}z^3+\cdots \]

令 \(z=y_x^2\)，得到：

\[ \frac{1}{\sqrt{1+y_x^2}}

1-\frac12 y_x^2+\frac38 y_x^4-\frac{5}{16}y_x^6+\cdots \]

所以方程展开为：

\[ y_{tt}

v^2 y_{xx} \left( 1-\frac12 y_x^2+\frac38 y_x^4-\frac{5}{16}y_x^6+\cdots \right) \]

也就是：

\[ \boxed{ y_{tt}

v^2 y_{xx} -\frac{v^2}{2}y_{xx}y_x^2 +\frac{3v^2}{8}y_{xx}y_x^4 -\cdots } \]

三、前 3 项分别是什么

第 1 项：线性波动项

\[ v^2 y_{xx} \]

对应经典线性波动方程：

\[ y_{tt}=v^2y_{xx} \]

这是微小振动下保留的主项。

第 2 项：三阶非线性项

\[ -\frac{v^2}{2}y_{xx}y_x^2 \]

注意：

  • \(y_{xx}\) 是二阶小量；
  • \(y_x^2\) 是二阶小量；
  • 合起来是四阶小量？这要看振幅量级。

更准确地说，若设振幅为 \(\epsilon\)，则：

\[ y\sim \epsilon,\qquad y_x\sim \epsilon,\qquad y_{xx}\sim \epsilon \]

所以：

\[ y_{xx}y_x^2\sim \epsilon^3 \]

而线性项 \(y_{xx}\sim \epsilon\)。
因此第二项相对于线性项是 \(\epsilon^2\) 阶修正。

所以它是三阶非线性项（在方程中出现的项整体为 \(\epsilon^3\) 量级）。

第 3 项：五阶非线性项

\[ +\frac{3v^2}{8}y_{xx}y_x^4 \]

因为：

\[ y_{xx}y_x^4\sim \epsilon^5 \]

所以它相对于线性项是 \(\epsilon^4\) 阶修正。

四、写成“线性 + 非线性修正”的形式

\[ \boxed{ y_{tt}-v^2y_{xx}

-\frac{v^2}{2}y_{xx}y_x^2 +\frac{3v^2}{8}y_{xx}y_x^4 -\cdots } \]

左边是线性波动算子：

\[ \Box y=y_{tt}-v^2y_{xx} \]

右边是非线性修正项。

五、再看另一个方程

前面还有：

\[ y_{tt}

\frac{T}{\mu}\frac{y_{xx}}{(1+y_x^2)^{3/2}} \]

令

\[ v^2=\frac{T}{\mu} \]

则：

\[ y_{tt}

v^2\frac{y_{xx}}{(1+y_x^2)^{3/2}} \]

利用：

\[ (1+z)^{-3/2}

1-\frac32 z+\frac{15}{8}z^2-\frac{35}{16}z^3+\cdots \]

令 \(z=y_x^2\)，得：

\[ \frac{1}{(1+y_x^2)^{3/2}}

1-\frac32 y_x^2+\frac{15}{8}y_x^4-\frac{35}{16}y_x^6+\cdots \]

所以：

\[ y_{tt}

v^2y_{xx} \left( 1-\frac32 y_x^2+\frac{15}{8}y_x^4-\frac{35}{16}y_x^6+\cdots \right) \]

展开：

\[ \boxed{ y_{tt}

v^2y_{xx} -\frac{3v^2}{2}y_{xx}y_x^2 +\frac{15v^2}{8}y_{xx}y_x^4 -\cdots } \]

即：

\[ \boxed{ y_{tt}-v^2y_{xx}

-\frac{3v^2}{2}y_{xx}y_x^2 +\frac{15v^2}{8}y_{xx}y_x^4 -\cdots } \]

六、两个展开对比

方程形式 线性项 第一个非线性修正 第二个非线性修正
\(y_{tt}=v^2\dfrac{y_{xx}}{\sqrt{1+y_x^2}}\) \(v^2y_{xx}\) \(-\dfrac{v^2}{2}y_{xx}y_x^2\) \(+\dfrac{3v^2}{8}y_{xx}y_x^4\)
\(y_{tt}=v^2\dfrac{y_{xx}}{(1+y_x^2)^{3/2}}\) \(v^2y_{xx}\) \(-\dfrac{3v^2}{2}y_{xx}y_x^2\) \(+\dfrac{15v^2}{8}y_{xx}y_x^4\)

区别只在系数，说明不同物理假设会给出不同强度的非线性修正。

七、这些项为什么叫“非线性项”

因为出现了 \(y_{xx}y_x^2\)、\(y_{xx}y_x^4\) 这种场量及其导数相乘的项。

线性方程中，每一项都只包含 \(y\) 或其导数的一次幂，所以可以叠加：

\[ y=y_1+y_2 \]

仍然满足方程。

但有了 \(y_{xx}y_x^2\) 后，若：

\[ y=y_1+y_2 \]

则：

\[ y_{xx}y_x^2

(y_{1,xx}+y_{2,xx})(y_{1,x}+y_{2,x})^2 \]

会出现交叉项，比如：

\[ y_{1,xx}y_{2,x}^2,\qquad y_{2,xx}y_{1,x}y_{2,x} \]

所以叠加原理失效。

八、如果做多尺度展开，会看到什么

设小参数 \(\epsilon\ll 1\)，令：

\[ y=\epsilon y_1+\epsilon^2 y_2+\epsilon^3 y_3+\cdots \]

代入：

\[ y_{tt}=v^2y_{xx}\left(1-\frac12 y_x^2+\cdots\right) \]

最低阶 \(O(\epsilon)\)：

\[ y_{1,tt}=v^2y_{1,xx} \]

这是线性波。

下一阶 \(O(\epsilon^3)\) 会出现：

\[ y_{3,tt}-v^2y_{3,xx}

-\frac{v^2}{2}y_{1,xx}y_{1,x}^2 \]

右端相当于一个“源项”，会驱动三倍频等谐波。

九、一个具体正弦波例子

设最低阶解为：

\[ y_1=A\cos(kx-\omega t) \]

其中：

\[ \omega^2=v^2k^2 \]

则：

\[ y_{1,x}=-Ak\sin(kx-\omega t) \]

\[ y_{1,xx}=-Ak^2\cos(kx-\omega t) \]

于是非线性源项：

\[ -\frac{v^2}{2}y_{1,xx}y_{1,x}^2

-\frac{v^2}{2} \left(-Ak^2\cos\phi\right) \left(A^2k^2\sin^2\phi\right) \]

\[

\frac{v^2A^3k^4}{2}\cos\phi\sin^2\phi \]

利用：

\[ \sin^2\phi=\frac{1-\cos2\phi}{2} \]

\[ \cos\phi\sin^2\phi

\frac12\cos\phi-\frac12\cos\phi\cos2\phi \]

再用：

\[ \cos\phi\cos2\phi=\frac12(\cos3\phi+\cos\phi) \]

得到：

\[ \cos\phi\sin^2\phi

\frac14\cos\phi-\frac14\cos3\phi \]

所以源项包含：

\[ \cos\phi \quad\text{和}\quad \cos3\phi \]

也就是说：
一个纯正弦波，由于非线性，会产生三次谐波。

十、总结成一句话

非线性绳波方程展开后，前几项是：

\[ \boxed{ y_{tt}

v^2y_{xx} -\frac{v^2}{2}y_{xx}y_x^2 +\frac{3v^2}{8}y_{xx}y_x^4 -\cdots } \]

或者另一种假设下：

\[ \boxed{ y_{tt}

v^2y_{xx} -\frac{3v^2}{2}y_{xx}y_x^2 +\frac{15v^2}{8}y_{xx}y_x^4 -\cdots } \]

它们的含义是：

  1. 第一项：线性波动；
  2. 第二项：最低阶非线性修正，导致波形畸变、谐波产生；
  3. 第三项：更高阶修正，通常更弱，但在大振幅下不可忽略。

这是否和微扰类似 是的，非常类似，但严格说它更接近非线性摄动 / 多尺度展开，而不是量子力学里常见的“微扰论”。两者在数学结构上同源，但物理含义和展开对象有区别。

一、相似之处：都是“小参数展开”

微扰论的核心思想是：

\[ \text{精确解}=\text{已知简单解}+\text{小修正} \]

比如量子力学中：

\[ H=H_0+\lambda H’ \]

波函数和能量按 \(\lambda\) 展开：

\[ E=E^{(0)}+\lambda E^{(1)}+\lambda^2 E^{(2)}+\cdots \]

\[ |\psi\rangle=|\psi^{(0)}\rangle+\lambda|\psi^{(1)}\rangle+\lambda^2|\psi^{(2)}\rangle+\cdots \]

我们刚才做的非线性绳波展开：

\[ y_{tt}=v^2y_{xx}\left(1-\frac12 y_x^2+\cdots\right) \]

设振幅小参数 \(\epsilon\)，令：

\[ y=\epsilon y_1+\epsilon^2 y_2+\epsilon^3 y_3+\cdots \]

最低阶：

\[ y_{1,tt}=v^2y_{1,xx} \]

下一阶出现非线性源项：

\[ y_{3,tt}-v^2y_{3,xx}

-\frac{v^2}{2}y_{1,xx}y_{1,x}^2 \]

这和微扰论的结构完全一样：

  • 零阶：已知线性解；
  • 高阶：由低阶解产生的“源项”驱动修正。

所以从数学框架看，是的，这就是一种摄动展开。

二、关键区别：微扰 vs 非线性摄动

方面 量子力学微扰论 非线性波摄动
展开对象 哈密顿量 \(H=H_0+\lambda H’\) 解 \(y=\epsilon y_1+\epsilon^2 y_2+\cdots\)
小参数来源 相互作用强度 \(\lambda\) 振幅 \(\epsilon\)
方程类型 线性方程 非线性方程
叠加原理 零阶仍可叠加 高阶出现交叉项，叠加失效
修正项 能量、波函数修正 谐波、波形畸变、波速修正
主要困难 久期项、简并 久期项、谐波共振、激波形成

最重要的区别是：

微扰论通常处理线性方程；非线性摄动处理的是非线性方程，非线性本身来自解的高阶项。

三、为什么说它“类似微扰”

因为操作步骤几乎一样：

1. 找零阶解

\[ y_1=A\cos(kx-\omega t) \]

2. 把零阶解代入非线性项

得到：

\[ -\frac{v^2}{2}y_{1,xx}y_{1,x}^2

\frac{v^2A^3k^4}{2}\cos\phi\sin^2\phi \]

3. 展开成谐波

\[ \cos\phi\sin^2\phi

\frac14\cos\phi-\frac14\cos3\phi \]

于是高阶方程变成：

\[ y_{3,tt}-v^2y_{3,xx}

\frac{v^2A^3k^4}{8}\cos\phi -\frac{v^2A^3k^4}{8}\cos3\phi \]

4. 解高阶方程

  • \(\cos3\phi\) 项产生三次谐波；
  • \(\cos\phi\) 项与齐次解同频，会产生久期项，需要重新定义频率或波速。

这和量子微扰论中处理共振、简并的方式非常像。

四、但有一个本质不同：久期项与重normalization

在量子微扰论中，如果微扰导致能级移动，通常要把修正吸收到能量里。

在非线性波中，同样的事情发生：

\[ y_{3,tt}-v^2y_{3,xx}

\frac{v^2A^3k^4}{8}\cos\phi \]

右端 \(\cos\phi\) 与齐次解同频，直接求解会出现：

\[ t\cos\phi \]

这种随时间线性增长的久期项，说明零阶频率 \(\omega=vk\) 需要修正。

正确做法是引入非线性色散关系：

\[ \omega

vk\left(1+\alpha A^2+\cdots\right) \]

其中 \(\alpha\) 由非线性项系数决定。

这相当于量子微扰论中：

\[ E_n=E_n^{(0)}+\lambda E_n^{(1)}+\lambda^2 E_n^{(2)}+\cdots \]

只是这里修正的是频率/波速，而不是能级。

五、更准确的名字：非线性摄动 / 多尺度方法

这类展开通常叫：

  1. 非线性摄动法
  2. 多尺度展开法
  3. Stokes 展开
  4. 弱非线性波理论

标准做法是引入多个时间尺度：

\[ T_0=t,\qquad T_1=\epsilon t,\qquad T_2=\epsilon^2 t,\cdots \]

设：

\[ y=\epsilon y_1+\epsilon^2 y_2+\epsilon^3 y_3+\cdots \]

对时间的导数变为：

\[ \frac{\partial}{\partial t}

\frac{\partial}{\partial T_0} +\epsilon\frac{\partial}{\partial T_1} +\epsilon^2\frac{\partial}{\partial T_2} +\cdots \]

这样可以把久期项消掉，得到缓慢变化的振幅和频率修正。

六、和量子微扰的对照表

步骤 量子微扰论 非线性波摄动
零阶 \(H_0\psi^{(0)}=E^{(0)}\psi^{(0)}\) \(y_{1,tt}=v^2y_{1,xx}\)
小参数 \(\lambda\) \(\epsilon=A\)
一阶方程 \((H_0-E^{(0)})\psi^{(1)}=(E^{(1)}-H’)\psi^{(0)}\) \(y_{2,tt}-v^2y_{2,xx}=0\) 或源项
高阶方程 由低阶解产生源项 由 \(y_{1,xx}y_{1,x}^2\) 产生源项
共振/久期 简并微扰、能级移动 同频源项、频率修正
最终修正 能级、波函数 波速、谐波、波形

七、一句话总结

是的，这和微扰非常类似，本质上是“非线性摄动展开”。

区别在于：

  • 微扰论通常展开的是线性方程中的小相互作用；
  • 非线性波摄动展开的是非线性方程中的小振幅解；
  • 两者都会遇到久期项，都需要把低阶效应吸收到零阶参数中；
  • 非线性波摄动还会产生高次谐波和波形畸变，这是线性微扰论中通常没有的。

所以你可以把它理解成：

\[ \boxed{ \text{非线性波摄动}

\text{微扰思想} + \text{非线性效应} + \text{多尺度处理} } \]

非线性波动方程时，y_x«1是什么意思 好问题。这里容易混淆的是：\(y_x\ll 1\) 到底是在说“振幅小”，还是“斜率小”？ 答案更准确地说，是斜率小，但它通常也隐含了振幅小和波长足够长的条件。

一、\(y_x\) 的几何意义

\[ y_x=\frac{\partial y}{\partial x} \]

它表示绳的局部斜率，也就是绳在该点处切线与水平方向夹角的正切：

\[ y_x=\tan\theta \]

所以：

\[ y_x\ll 1 \quad\Longleftrightarrow\quad |\tan\theta|\ll 1 \quad\Longleftrightarrow\quad |\theta|\ll 1 \]

也就是说，绳上每一点的切线都几乎水平，绳的倾斜非常小。

二、它不等于“振幅小”

这是最关键的一点。

振幅小，指的是：

\[ |y|\ll L \]

其中 \(L\) 是某个特征长度，比如绳长或波长。

但 \(y_x\ll 1\) 说的是：

\[ \left|\frac{\partial y}{\partial x}\right|\ll 1 \]

这两者不是一回事。

反例：振幅很小但斜率很大

假设：

\[ y(x)=A\sin(kx) \]

则：

\[ y_x=Ak\cos(kx) \]

所以：

\[ |y_x|_{\max}=Ak \]

即使 \(A\) 很小，如果 \(k\) 很大，也就是波长很短，那么：

\[ Ak\gg 1 \]

仍然会有 \(y_x\gg 1\)。

这时振幅虽然小，但绳在很短距离内上下剧烈起伏，斜率很大，微小振动近似失效。

正例：振幅大但斜率小

反过来，如果波长很长：

\[ k\ll 1 \]

即使 \(A\) 不算特别小，只要：

\[ Ak\ll 1 \]

那么：

\[ y_x\ll 1 \]

仍然成立。

此时绳整体起伏可能很大，但局部倾斜很小，线性波动方程仍然近似成立。

三、对于正弦波，条件具体是什么

设：

\[ y=A\sin(kx-\omega t) \]

则：

\[ y_x=Ak\cos(kx-\omega t) \]

所以：

\[ |y_x|_{\max}=Ak \]

因此：

\[ y_x\ll 1 \quad\Longleftrightarrow\quad Ak\ll 1 \]

即：

\[ \boxed{ A\ll \frac{1}{k}=\frac{\lambda}{2\pi} } \]

其中 \(\lambda=2\pi/k\) 是波长。

也就是说：

振幅远小于波长除以 \(2\pi\)。

更粗略地说：

\[ \boxed{ A\ll \lambda } \]

即振幅远小于波长。

四、为什么这个条件能让非线性项变小

回到展开：

\[ y_{tt}

v^2y_{xx} \left( 1-\frac12 y_x^2+\frac38 y_x^4-\cdots \right) \]

非线性修正项相对于线性项的比值为：

\[ \frac{\text{非线性项}}{\text{线性项}} \sim y_x^2 \]

所以：

\[ y_x\ll 1 \quad\Longrightarrow\quad y_x^2\ll 1 \]

非线性修正就很小，可以忽略，方程退化为：

\[ y_{tt}\approx v^2y_{xx} \]

这就是线性波动方程。

如果：

\[ y_x\sim 1 \]

那么：

\[ y_x^2\sim 1 \]

非线性项和线性项同量级，绝对不能忽略，必须用非线性波动方程。

五、物理图像

\(y_x\ll 1\) 时

  • 绳几乎水平；
  • 张力方向几乎沿水平；
  • \(\sin\theta\approx\tan\theta\approx\theta\)；
  • 弧长近似等于水平投影： \[ ds=\sqrt{1+y_x^2},dx\approx dx \]
  • 微元质量近似为 \(\mu dx\)；
  • 线性波动方程成立。

\(y_x\) 不再远小于 1 时

  • 绳局部倾斜明显；
  • \(\sin\theta\) 和 \(\tan\theta\) 差别不可忽略；
  • 弧长明显大于水平投影： \[ ds=\sqrt{1+y_x^2},dx>dx \]
  • 微元质量随位置变化；
  • 可能出现纵向运动；
  • 非线性效应显现。

六、和“小振幅”的关系

通常教材说“微小振动”时，实际上包含两个条件：

  1. 振幅小： \[ |y|\ll L \]

  2. 斜率小： \[ |y_x|\ll 1 \]

其中第二个条件才是线性化能否成立的关键。

因为波动方程里的非线性项来自：

\[ \sqrt{1+y_x^2} \]

而不是直接来自 \(y\) 本身。

所以更准确的说法是：

线性波动方程要求的是斜率小，而不只是振幅小。

七、总结

\[ \boxed{ y_x\ll 1 } \]

意思是：

  1. 绳的局部斜率很小，切线几乎水平；
  2. 等价于 \(\theta\ll 1\)，即 \(\tan\theta\approx\sin\theta\approx\theta\)；
  3. 对于正弦波 \(y=A\sin(kx-\omega t)\)，等价于： \[ Ak\ll 1 \] 即： \[ A\ll \frac{\lambda}{2\pi} \] 也就是振幅远小于波长；
  4. 它保证了： \[ \sqrt{1+y_x^2}\approx 1 \] 从而非线性项 \(y_x^2\)、\(y_x^4\) 可以忽略，方程退化为线性波动方程。

一句话：

\(y_x\ll 1\) 是“斜率小”的条件，物理上要求振幅相对于波长足够小，而不是单纯振幅小。

|y’|«1的含义，是否太苛刻，也就是几乎不存在垂直振动 \(\lvert y’\rvert\ll 1\) 指的是斜率小，不是振幅 \(y\) 本身接近 0。垂直振动完全可以很大，只要波长更长。

\(y’\) 是什么

\[ y’=\frac{\partial y}{\partial x}=\tan\theta\approx\theta \]

（弧度）。它是相邻两点的高度差除以水平距离，也就是波形有多陡。

单色波 \(y=A\cos(kx-\omega t)\) 时

\[ \lvert y’\rvert_{\max}=kA=\frac{2\pi A}{\lambda}. \]

所以条件是

\[ \frac{A}{\lambda}\ll\frac{1}{2\pi}\approx0.16, \]

即振幅比波长小很多，不是振幅比“几乎为零”还小。

并不苛刻的例子

\(A\) \(\lambda\) \(kA\) 小斜率？
1 mm 1 m \(0.006\) 是
1 cm 1 m \(0.06\) 是
10 cm 2 m \(0.31\) 勉强
10 cm 20 cm \(3.1\) 否，必须用 \(\sin\theta=y_x/\sqrt{1+y_x^2}\)

吉他弦、绳子上的普通波包、水槽里不太陡的涟漪，\(A/\lambda\) 通常是百分之几，线性横波够用。
只有陡前沿、大甩动、快要打圈的弦，才谈得上“几乎不存在可线性化的竖直振动”。

和“几乎没有竖直运动”的差别

  • \(\lvert y\rvert\ll a\) 或 \(\lvert y\rvert\ll L\)：振幅相对格距或全长很小——这才是“振动很弱”。
  • \(\lvert y’\rvert\ll 1\)：局部像一条缓坡。同一条缓坡可以是“整根弦抬高 10 cm、但 2 m 才完成一个起伏”。

几何缩短 \(u_x\approx-\frac12 y_x^2\) 是 \(O((A/\lambda)^2)\)。\(kA=0.1\) 时这项约 \(0.005\)，线性理论丢掉它很合理；\(kA\sim 1\) 时不能丢。

一句话

\(\lvert y’\rvert\ll 1\) 是“波形平缓”，不是“不要竖直振动”。竖直位移可以明显看见，只要起伏铺在足够长的一段弦上。真正苛刻的是短波大振幅（尖峰、甩鞭），那时必须用上一问的非线性方程。

那上面的非线性方程是否符合洛伦兹变换 好问题。答案是：

\[ \boxed{ \text{上面那个非线性绳波方程通常不满足洛伦兹协变性。} } \]

也就是说，它不是在洛伦兹变换下保持形式不变的。原因在于它的非线性项含有 \(y_x\)，而 \(y_x\) 在洛伦兹变换下会混入时间导数，导致方程形式改变。

下面具体看。

一、回顾非线性方程

我们之前得到的一个典型形式是：

\[ y_{tt}=c^2\frac{y_{xx}}{\sqrt{1+y_x^2}} \]

其中：

\[ y_x=\frac{\partial y}{\partial x} \]

这里 \(c^2=T_0/\mu\) 是线性波速。

这个方程来自绳的几何非线性，是一种非相对论性模型。

二、洛伦兹变换下的关键问题

洛伦兹变换：

\[ x’=\gamma(x-vt) \]

\[ t’=\gamma\left(t-\frac{v}{c^2}x\right) \]

在这个变换下：

\[ \partial_x

\gamma\left(\partial_{x’}-\frac{v}{c^2}\partial_{t’}\right) \]

\[ \partial_t

\gamma\left(\partial_{t’}-v\partial_{x’}\right) \]

因此：

\[ y_x

\gamma\left(y_{x’}-\frac{v}{c^2}y_{t’}\right) \]

也就是说：

在洛伦兹变换下，\(y_x\) 不再只是空间导数，它会混入时间导数 \(y_{t’}\)。

这是关键。

三、非线性项会变成什么样

原方程中的非线性因子：

\[ \frac{1}{\sqrt{1+y_x^2}} \]

在 \(S’\) 系中变成：

\[ \frac{1}{\sqrt{1+\gamma^2\left(y_{x’}-\frac{v}{c^2}y_{t’}\right)^2}} \]

它不再只是 \(y_{x’}\) 的函数，而是同时依赖 \(y_{x’}\) 和 \(y_{t’}\)。

所以方程在 \(S’\) 系中会变成：

\[ y_{t’t’}=c^2 \frac{y_{x’x’}}{\sqrt{1+\gamma^2\left(y_{x’}-\frac{v}{c^2}y_{t’}\right)^2}} +\text{其他项} \]

这显然和原来的形式：

\[ y_{t’t’}=c^2\frac{y_{x’x’}}{\sqrt{1+y_{x’}^2}} \]

不一样。

四、为什么线性波动方程可以，非线性不行

线性波动方程：

\[ y_{tt}=c^2y_{xx} \]

在洛伦兹变换下，二阶导数组合：

\[ y_{tt}-c^2y_{xx} \]

是一个洛伦兹标量（更准确地说，是达朗贝尔算符作用在标量场上）：

\[ \Box y

\frac{\partial^2 y}{\partial t^2} -c^2\frac{\partial^2 y}{\partial x^2} \]

它在洛伦兹变换下保持不变：

\[ \Box y=\Box’ y \]

所以线性波动方程是洛伦兹协变的。

但非线性方程中出现了：

\[ y_x \]

单独一个一阶空间导数，不是洛伦兹协变的。

在洛伦兹变换下，\(y_x\) 会和时间导数混合，因此任何只依赖 \(y_x\) 的函数都不可能保持形式不变。

五、什么样的非线性方程才可能洛伦兹协变

要让非线性方程洛伦兹协变，非线性项必须由洛伦兹标量构成。

例如：

  1. 相对论性标量场方程：

\[ \Box\phi+\frac{dV}{d\phi}=0 \]

其中：

\[ \Box=\partial_t^2-c^2\nabla^2 \]

是非线性项来自势能 \(V(\phi)\)，而 \(\phi\) 本身是洛伦兹标量。

  1. sine-Gordon 方程：

\[ \phi_{tt}-c^2\phi_{xx}+\sin\phi=0 \]

这里 \(\sin\phi\) 是 \(\phi\) 的函数，而 \(\phi\) 是标量，所以洛伦兹协变。

  1. \(\phi^4\) 模型：

\[ \Box\phi+\lambda\phi^3=0 \]

也是洛伦兹协变的。

关键区别：

非线性项必须是洛伦兹标量的函数，而不是某个非协变导数的函数。

六、绳波方程为什么不是洛伦兹协变

绳波方程来自非相对论力学：

  • 牛顿第二定律；
  • 微小振动近似；
  • 弦的张力；
  • 几何非线性。

这些都不是相对论性的。

具体地：

  1. 绳的微元质量： \[ dm=\mu,ds \] 这里的 \(\mu\) 是静止线密度，但在运动参考系中会有长度收缩，而绳波方程没有考虑这一点。

  2. 张力方向： \[ \sin\theta=\frac{y_x}{\sqrt{1+y_x^2}} \] 这个几何关系是在经典欧几里得空间中写的，不是洛伦兹协变的。

  3. 时间： \[ t \] 是牛顿绝对时间，不是洛伦兹变换下的坐标时。

所以：

\[ \boxed{ \text{绳波非线性方程是伽利略协变或非协变的，不是洛伦兹协变的。} } \]

七、如果强行要求洛伦兹协变，方程会变成什么

如果把绳看作相对论性弦，需要从相对论性作用量出发：

\[ S=-\int T\sqrt{1-\frac{1}{c^2}\left(\frac{\partial y}{\partial t}\right)^2} ,dx,dt \]

得到的方程会包含：

\[ \frac{\partial}{\partial t} \left( \frac{y_t}{\sqrt{1-y_t^2/c^2}} \right)

c^2\frac{\partial}{\partial x} \left( \frac{y_x}{\sqrt{1+y_x^2}} \right) \]

这也不是我们之前那个绳波方程。

所以：

经典非线性绳波方程不是相对论性方程，因此不满足洛伦兹协变。

八、总结

方程 洛伦兹协变
\(y_{tt}=c^2y_{xx}\) 是
\(y_{tt}=c^2\dfrac{y_{xx}}{\sqrt{1+y_x^2}}\) 否
\(y_{tt}=c^2\dfrac{y_{xx}}{(1+y_x^2)^{3/2}}\) 否
\(\Box\phi+\sin\phi=0\) 是
\(\Box\phi+\lambda\phi^3=0\) 是

关键原因：

\[ \boxed{ \text{非线性绳波方程中的 }y_x\text{ 不是洛伦兹标量，} } \]

在洛伦兹变换下它会和时间导数混合，因此方程形式改变。

所以：

\[ \boxed{ \text{经典非线性绳波方程不满足洛伦兹协变性。} } \]

它属于非相对论性非线性波动方程。