引言：从一根振动的吉他弦开始
想象一下，你拨动吉他的一根弦。弦开始振动，发出优美的声音。如果你用高速摄像机拍摄这个过程，会看到弦的形状随时间不断变化：向上弯曲，向下弯曲，再向上弯曲……这种运动有什么规律？
更具体地说，如果已知某个时刻弦的形状，你能预测下一时刻它的形状吗？这个问题看似简单，但它引领我们走向数学物理中最重要的方程之一——波动方程。
在 18 世纪，几位伟大的数学家——达朗贝尔（d’Alembert）、欧拉（Euler）和伯努利（Bernoulli）——都在思考这个问题。他们的答案不仅解释了弦振动，还为声学、光学、地震学甚至量子力学奠定了基础。
让我们从这根弦开始，一步步揭开波动方程的面纱。
第一章：波动的物理本质
什么是波？
在开始推导方程之前，我们需要明确：什么是波？
波是振动在空间中的传播。当某个点的物理量（如位移、压力、电场等）随时间振动时，这种振动会影响周围的点，并传播出去。波不需要物质的长距离移动，它传播的是能量和信息。
想象一下水面上的波纹。当你往平静的水面投一块石子，水并没有整体移动，但波纹会一圈圈扩散开来——这就是波的传播。
波的分类
波可以分为两大类：
横波（Transverse Wave）：振动方向与传播方向垂直
- 例子：吉他弦振动、光波
- 特点：弦上下的振动，波沿弦的方向传播
纵波（Longitudinal Wave）：振动方向与传播方向平行
- 例子：声波（空气分子的振动）
- 特点：空气分子沿声音传播方向前后振动
波的基本性质
描述波的几个关键参数：
- 频率 $f$：单位时间内振动的次数（单位：赫兹 Hz）
- 周期 $T = \frac{1}{f}$：完成一次振动所需的时间
- 波长 $\lambda$：波完成一个周期在空间中传播的距离
- 波速 $c$：波传播的速度，满足 $c = f\lambda$
- 振幅 $A$：波偏离平衡位置的最大值
这些参数不是孤立的，它们通过波动方程联系在一起。
第二章：一维波动方程的诞生
牛顿第二定律与弦的振动
考虑一根均匀的弦，两端固定（比如吉他弦）。设弦的线密度（单位长度的质量）为 $\rho$，张力为 $T_0$。弦在平衡时是一条直线。
当弦发生微小振动时，设弦上位置 $x$、时间 $t$ 的横向位移为 $u(x, t)$。我们的目标是推导 $u(x, t)$ 满足的方程。
取弦上从 $x$ 到 $x + \Delta x$ 的一小段。这一段的长度约为 $\Delta x$，质量为 $\rho \Delta x$。
根据牛顿第二定律（$F = ma$），这一小段的运动方程为：
$$ \rho \Delta x \frac{\partial^2 u}{\partial t^2} = F_{\text{net}} $$
其中 $F_{\text{net}}$ 是作用在这段弦上的净力。
张力的作用
弦上每一点都受到张力。张力沿着弦的切线方向。考虑这一小段两端：
- 在 $x$ 处，张力为 $T_0$，与水平方向的夹角为 $\theta_1$
- 在 $x + \Delta x$ 处，张力为 $T_0$，与水平方向的夹角为 $\theta_2$
假设振动很小，角度 $\theta$ 也很小。此时：
- 水平方向的张力分量：$T_0 \cos\theta \approx T_0$（因为 $\cos\theta \approx 1$）
- 垂直方向的张力分量：$T_0 \sin\theta \approx T_0 \tan\theta = T_0 \frac{\partial u}{\partial x}$（因为 $\sin\theta \approx \tan\theta$）
垂直方向的净力为：
$$ F_{\text{net}} = T_0 \sin\theta_2 - T_0 \sin\theta_1 \approx T_0 \left( \frac{\partial u}{\partial x}(x + \Delta x, t) - \frac{\partial u}{\partial x}(x, t) \right) $$
注意到右边括号内就是 $\frac{\partial u}{\partial x}$ 在 $[x, x + \Delta x]$ 上的变化量，可以写成：
$$ \frac{\partial u}{\partial x}(x + \Delta x, t) - \frac{\partial u}{\partial x}(x, t) = \frac{\partial^2 u}{\partial x^2} \Delta x $$
因此：
$$ F_{\text{net}} = T_0 \frac{\partial^2 u}{\partial x^2} \Delta x $$
波动方程
将净力代入牛顿第二定律：
$$ \rho \Delta x \frac{\partial^2 u}{\partial t^2} = T_0 \frac{\partial^2 u}{\partial x^2} \Delta x $$
两边除以 $\rho \Delta x$，令 $\Delta x \to 0$：
$$ \frac{\partial^2 u}{\partial t^2} = \frac{T_0}{\rho} \frac{\partial^2 u}{\partial x^2} $$
定义波速 $c = \sqrt{\frac{T_0}{\rho}}$，就得到了著名的一维波动方程：
$$ \boxed{\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}} $$
这个方程告诉我们：位移对时间的二阶导数（加速度）与位移对空间的二阶导数（曲率）成正比。弦越弯曲的地方，加速度越大。
第三章：达朗贝尔公式
特征线法
1746 年，法国数学家达朗贝尔（Jean le Rond d’Alembert）发现了一个优雅的方法来解这个方程。他的思路是：找到一组新的变量，使得方程变得更容易处理。
做变量替换：
$$ \xi = x - ct, \quad \eta = x + ct $$
其中 $c$ 是波速。这两个新的变量沿着"特征线"变化：
- $\xi = \text{常数}$：右行波的特征线（向右传播的波）
- $\eta = \text{常数}$：左行波的特征线（向左传播的波）
坐标变换
计算偏导数：
$$ \frac{\partial}{\partial x} = \frac{\partial \xi}{\partial x} \frac{\partial}{\partial \xi} + \frac{\partial \eta}{\partial x} \frac{\partial}{\partial \eta} = \frac{\partial}{\partial \xi} + \frac{\partial}{\partial \eta} $$
$$ \frac{\partial}{\partial t} = \frac{\partial \xi}{\partial t} \frac{\partial}{\partial \xi} + \frac{\partial \eta}{\partial t} \frac{\partial}{\partial \eta} = -c \frac{\partial}{\partial \xi} + c \frac{\partial}{\partial \eta} $$
二阶导数：
$$ \frac{\partial^2 u}{\partial x^2} = \frac{\partial^2 u}{\partial \xi^2} + 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} $$
$$ \frac{\partial^2 u}{\partial t^2} = c^2 \left( \frac{\partial^2 u}{\partial \xi^2} - 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} \right) $$
简化的方程
将二阶导数代入波动方程：
$$ c^2 \left( \frac{\partial^2 u}{\partial \xi^2} - 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} \right) = c^2 \left( \frac{\partial^2 u}{\partial \xi^2} + 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} \right) $$
化简：
$$ \frac{\partial^2 u}{\partial \xi^2} - 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} = \frac{\partial^2 u}{\partial \xi^2} + 2 \frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} $$
$$ -4 \frac{\partial^2 u}{\partial \xi \partial \eta} = 0 \quad \Rightarrow \quad \frac{\partial^2 u}{\partial \xi \partial \eta} = 0 $$
达朗贝尔公式
积分这个方程：
$$ \frac{\partial u}{\partial \eta} = f(\eta) $$
再对 $\eta$ 积分：
$$ u(\xi, \eta) = f(\eta) + g(\xi) $$
其中 $f$ 和 $g$ 是任意函数。换回原变量：
$$ \boxed{u(x, t) = f(x + ct) + g(x - ct)} $$
这就是达朗贝尔公式！
物理意义
这个公式的物理意义非常深刻：
$g(x - ct)$：右行波，以速度 $c$ 向右传播
- 在 $t=0$ 时，形状为 $g(x)$
- 在 $t$ 时刻，形状相同，但向右平移了 $ct$
$f(x + ct)$：左行波，以速度 $c$ 向左传播
总波是两个行波的叠加。
例如，如果初始时刻弦被拨动成某个形状 $u(x, 0) = \phi(x)$，并且初始速度为零 $u_t(x, 0) = 0$，那么解为：
$$ u(x, t) = \frac{1}{2} \phi(x + ct) + \frac{1}{2} \phi(x - ct) $$
初始形状分裂成两个波，一个向左传播，一个向右传播，振幅各减半。
第四章：分离变量法与傅里叶级数
达朗贝尔公式适用于无限长的弦。对于两端固定的弦（如吉他弦），我们需要考虑边界条件。
边界条件
设弦的两端分别固定在 $x=0$ 和 $x=L$：
$$ u(0, t) = 0, \quad u(L, t) = 0, \quad \forall t > 0 $$
分离变量法
假设解可以写成空间部分和时间部分的乘积：
$$ u(x, t) = X(x) \cdot T(t) $$
代入波动方程 $\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}$：
$$ X(x) T’’(t) = c^2 X’’(x) T(t) $$
两边除以 $c^2 X(x) T(t)$：
$$ \frac{T’’(t)}{c^2 T(t)} = \frac{X’’(x)}{X(x)} $$
左边只依赖 $t$，右边只依赖 $x$，必须都等于同一个常数。设这个常数为 $-\lambda$：
$$ \frac{T’’(t)}{c^2 T(t)} = \frac{X’’(x)}{X(x)} = -\lambda $$
这给了我们两个常微分方程：
空间方程： $$ X’’(x) + \lambda X(x) = 0 $$
时间方程： $$ T’’(t) + c^2 \lambda T(t) = 0 $$
空间方程的解
空间方程的解取决于 $\lambda$ 的符号。为了得到有物理意义的解，我们取 $\lambda > 0$。令 $\lambda = k^2$，则：
$$ X(x) = A \cos(kx) + B \sin(kx) $$
应用边界条件 $X(0) = 0$：
$$ A \cos(0) + B \sin(0) = A = 0 \quad \Rightarrow \quad A = 0 $$
因此 $X(x) = B \sin(kx)$。
应用边界条件 $X(L) = 0$：
$$ B \sin(kL) = 0 $$
要得到非零解，必须有 $\sin(kL) = 0$，即：
$$ kL = n\pi \quad \Rightarrow \quad k_n = \frac{n\pi}{L}, \quad n = 1, 2, 3, \ldots $$
因此特征值和特征函数为：
$$ \lambda_n = \left(\frac{n\pi}{L}\right)^2, \quad X_n(x) = \sin\left(\frac{n\pi x}{L}\right) $$
时间方程的解
时间方程为：
$$ T’’(t) + c^2 \lambda_n T(t) = 0 \quad \Rightarrow \quad T’’(t) + c^2 \left(\frac{n\pi}{L}\right)^2 T(t) = 0 $$
这是简谐振动的方程，解为：
$$ T_n(t) = C_n \cos\left(\frac{n\pi c t}{L}\right) + D_n \sin\left(\frac{n\pi c t}{L}\right) $$
叠加原理
由于方程是线性的，一般解是这些特解的叠加：
$$ u(x, t) = \sum_{n=1}^{\infty} \sin\left(\frac{n\pi x}{L}\right) \left[ A_n \cos\left(\frac{n\pi c t}{L}\right) + B_n \sin\left(\frac{n\pi c t}{L}\right) \right] $$
驻波与固有频率
每一项都代表一个驻波（standing wave）：
$$ u_n(x, t) = \sin\left(\frac{n\pi x}{L}\right) \cos\left(\omega_n t\right) $$
其中 $\omega_n = \frac{n\pi c}{L}$ 是第 $n$ 个固有频率。
驻波的特点是：波不传播，而是原地振动。弦上有一些点始终静止，称为节点（nodes）；有一些点振动幅度最大，称为腹点（antinodes）。
对于吉他弦：
- $n=1$：基频，声音最低
- $n=2$：第一泛音，频率是基频的 2 倍
- $n=3$：第二泛音，频率是基频的 3 倍
音乐中的泛音就是这些驻波！
第五章：扩展到多维空间
二维波动方程
对于薄膜（如鼓皮），波动方程扩展到二维：
$$ \frac{\partial^2 u}{\partial t^2} = c^2 \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} \right) = c^2 \nabla^2 u $$
其中 $\nabla^2 = \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}$ 是二维拉普拉斯算子。
圆形鼓膜的固有频率不再是简单的整数倍，而是与贝塞尔函数（Bessel functions）有关。这就是为什么鼓的声音不如弦乐器"纯粹"。
三维波动方程
在三维空间中，波动方程为：
$$ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u = c^2 \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right) $$
这是声波方程（acoustic wave equation），描述了声音在空气、水等介质中的传播。
球对称情况：球面波
对于球对称的情况（如点声源），使用球坐标系，波动方程简化为：
$$ \frac{\partial^2 u}{\partial t^2} = c^2 \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial u}{\partial r} \right) $$
做变量替换 $v(r, t) = r u(r, t)$，得到：
$$ \frac{\partial^2 v}{\partial t^2} = c^2 \frac{\partial^2 v}{\partial r^2} $$
这与一维波动方程形式相同！解为：
$$ v(r, t) = f(r - ct) + g(r + ct) $$
因此：
$$ u(r, t) = \frac{f(r - ct)}{r} + \frac{g(r + ct)}{r} $$
球面波的振幅随距离衰减，与 $1/r$ 成正比。这就是为什么远处传来的声音会越来越小。
第六章：应用与推广
1. 声学
声波是三维波动方程的经典应用。从音乐厅的声学设计到降噪技术，波动方程无处不在。
- 声速：在空气中约为 340 m/s，取决于温度和气压
- 多普勒效应：当声源和观察者相对运动时，频率发生变化
- 声学共振：建筑物、乐器中的共振现象
2. 光学与电磁波
麦克斯韦方程组（Maxwell’s equations）推导出的电磁波方程也是波动方程：
$$ \nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0 $$
$$ \nabla^2 \mathbf{B} - \frac{1}{c^2} \frac{\partial^2 \mathbf{B}}{\partial t^2} = 0 $$
其中 $\mathbf{E}$ 是电场，$\mathbf{B}$ 是磁场，$c$ 是光速。
光波本质上是电磁波！
3. 地震学
地震波传播用波动方程描述。主要有两种类型：
- P 波（纵波）：速度快，先到达
- S 波（横波）：速度慢，后到达
通过分析地震波的传播路径，可以探测地球内部结构。
4. 量子力学
薛定谔方程（Schrödinger equation）虽然不是波动方程，但形式类似：
$$ i\hbar \frac{\partial \psi}{\partial t} = -\frac{\hbar^2}{2m} \nabla^2 \psi + V \psi $$
其中 $\psi$ 是波函数，描述粒子的量子态。在自由空间（$V=0$）中，薛定谔方程描述了物质波的传播。
5. 波的反射与折射
当波遇到介质界面时，会发生反射和折射。通过波动方程和边界条件，可以推导出：
- 反射定律：入射角等于反射角
- 折射定律（斯涅尔定律）：$\frac{\sin\theta_1}{v_1} = \frac{\sin\theta_2}{v_2}$
6. 干涉与衍射
波的叠加原理导致了许多有趣的现象：
- 干涉：两个波相遇时，振幅叠加。相长干涉（波峰对波峰）和相消干涉（波峰对波谷）
- 衍射：波遇到障碍物时会"绕过"障碍物
这些现象可以用波动方程精确描述。
第七章：数值方法简介
对于复杂的几何形状或非均匀介质，解析解很难找到，这时需要数值方法。
有限差分法
将时间和空间离散化，用差分近似导数：
$$ \frac{\partial^2 u}{\partial t^2} \approx \frac{u_i^{n+1} - 2u_i^n + u_i^{n-1}}{\Delta t^2} $$
$$ \frac{\partial^2 u}{\partial x^2} \approx \frac{u_{i+1}^n - 2u_i^n + u_{i-1}^n}{\Delta x^2} $$
代入波动方程，得到显式格式：
$$ u_i^{n+1} = 2u_i^n - u_i^{n-1} + r^2 (u_{i+1}^n - 2u_i^n + u_{i-1}^n) $$
其中 $r = \frac{c \Delta t}{\Delta x}$ 是 Courant 数。
稳定性条件
为了保证数值稳定，必须满足 CFL 条件（Courant-Friedrichs-Lewy condition）：
$$ r = \frac{c \Delta t}{\Delta x} \leq 1 $$
时间步长不能太大，否则数值解会发散。
结语：从弦振动到宇宙波动
波动方程的伟大之处在于，它用简洁的数学语言统一描述了形形色色的波动现象。从吉他弦的振动，到声波的传播；从光波的干涉，到地震的探测；从量子世界的粒子波，到宇宙早期的引力波——都遵循着相同的数学规律。
从达朗贝尔在 18 世纪的发现，到今天在 5G 通信、医学成像、地震预警中的应用，波动方程已经走过了近三百年的历史。它告诉我们：自然界的波动现象虽然看起来千差万别，但背后有着深刻的统一性。
下次当你听到音乐、看到光波、感受到声波时，你可以自豪地说：“我知道这背后的方程——它描述的不仅仅是波，而是宇宙最基本的规律之一。”
参考资料
- d’Alembert, J. (1747). Recherches sur la courbe que forme une corde tendue mise en vibration. Histoire de l’Académie Royale des Sciences et Belles-Lettres de Berlin, 3, 214-219.
- Strauss, W. A. (2007). Partial Differential Equations: An Introduction (2nd ed.). Hoboken, NJ: Wiley.
- Evans, L. C. (2010). Partial Differential Equations (2nd ed.). Providence, RI: American Mathematical Society.
- Courant, R., & Hilbert, D. (1962). Methods of Mathematical Physics, Volume II: Partial Differential Equations. New York: Interscience Publishers.