📚 Wave Equation: Basic Forms and Solution Methods | 波动方程的基本形式与解法
The wave equation is one of the most fundamental partial differential equations in mathematics and physics. It describes how waves propagate through a medium, from sound waves in air to light waves in electromagnetic fields. Understanding its basic forms and solution methods is essential for students preparing for advanced mathematics examinations.
波动方程是数学和物理学中最基本的偏微分方程之一。它描述了波如何在介质中传播，从空气中的声波到电磁场中的光波。理解其基本形式和求解方法对于准备高级数学考试的学生至关重要。
1. What is the Wave Equation? | 什么是波动方程？
The one-dimensional wave equation is typically written as:
一维波动方程通常写成以下形式：
∂²u/∂t² = c² ∂²u/∂x²
Here, u(x,t) represents the displacement of the wave at position x and time t. The constant c represents the wave speed. This equation states that the acceleration of the wave at any point is proportional to the curvature of the wave at that point.
其中u(x,t)表示波在位置x和时间t处的位移，常数c表示波速。该方程表明波在任何一点的加速度与该点的曲率成正比。
For higher dimensions, the wave equation generalizes to:
对于更高维度，波动方程推广为：
∂²u/∂t² = c² ∇²u
where ∇² is the Laplacian operator. In three dimensions, this becomes ∂²u/∂t² = c²(∂²u/∂x² + ∂²u/∂y² + ∂²u/∂z²).
其中∇²是拉普拉斯算子。在三维情况下，它变成∂²u/∂t² = c²(∂²u/∂x² + ∂²u/∂y² + ∂²u/∂z²)。
2. Physical Derivation | 物理推导
The wave equation can be derived from physical principles. Consider a vibrating string with linear density ρ and tension T. Applying Newton’s second law to a small segment of the string yields:
波动方程可以从物理原理推导出来。考虑一根线密度为ρ、张力为T的振动弦。对弦的一小段应用牛顿第二定律可得：
∂²u/∂t² = (T/ρ) ∂²u/∂x²
Thus, the wave speed for a string is c = √(T/ρ). This shows that the wave speed increases with tension and decreases with density.
因此，弦上的波速为c = √(T/ρ)。这表明波速随张力增大而增大，随密度增大而减小。
Similarly, for sound waves in a gas, the wave speed is c = √(γP/ρ), where γ is the adiabatic index and P is the pressure. For light waves in a vacuum, c is the speed of light, approximately 3 × 10⁸ m/s.
类似地，对于气体中的声波，波速为c = √(γP/ρ)，其中γ是绝热指数，P是压强。对于真空中的光波，c是光速，约为3 × 10⁸ m/s。
3. D’Alembert’s Solution | 达朗贝尔解
One of the most elegant solution methods for the one-dimensional wave equation is D’Alembert’s formula. It uses a change of variables to transform the wave equation into a simpler form.
一维波动方程最优雅的求解方法之一是达朗贝尔公式。它通过变量变换将波动方程转化为更简单的形式。
Define new variables ξ = x – ct and η = x + ct. Using the chain rule, the wave equation becomes:
定义新变量ξ = x – ct和η = x + ct。利用链式法则，波动方程变为：
∂²u/∂ξ∂η = 0
This equation can be integrated directly. The general solution is:
这个方程可以直接积分。其通解为：
u(x,t) = f(x – ct) + g(x + ct)
where f and g are twice differentiable arbitrary functions. The function f(x – ct) represents a wave traveling to the right at speed c, while g(x + ct) represents a wave traveling to the left at speed c.
其中f和g是任意二阶可微函数。函数f(x – ct)表示以速度c向右传播的波，而g(x + ct)表示以速度c向左传播的波。
4. Initial Value Problem | 初值问题
For the wave equation on an infinite string, we impose initial conditions:
对于无限长弦上的波动方程，我们施加初始条件：
u(x,0) = φ(x), ∂u/∂t(x,0) = ψ(x)
Here, φ(x) is the initial displacement and ψ(x) is the initial velocity. D’Alembert’s formula for this initial value problem is:
其中φ(x)是初始位移，ψ(x)是初始速度。该初值问题的达朗贝尔公式为：
u(x,t) = ½[φ(x – ct) + φ(x + ct)] + (1/2c)∫₍ₓ₋꜀ₜ₎₍ₓ₊꜀ₜ₎ ψ(s) ds
This formula shows that the solution consists of two parts: the average of the initial displacement at x – ct and x + ct, plus an integral of the initial velocity over the interval [x – ct, x + ct]. This interval is called the domain of dependence.
该公式表明解由两部分组成：初始位移在x – ct和x + ct处的平均值，加上初始速度在区间[x – ct, x + ct]上的积分。这个区间称为依赖区间。
Let us verify this satisfies the wave equation. Computing the second derivatives with respect to t and x and substituting, we find that both sides are equal to c² times the second partial derivative with respect to x. This verification is a common exam question.
我们验证这个解满足波动方程。计算对t和对x的二阶偏导数并代入，我们发现两边都等于c²乘以对x的二阶偏导数。这种验证是常见的考试题目。
5. Separation of Variables | 分离变量法
For boundary value problems on finite intervals, the method of separation of variables is the standard approach. Consider a string of length L fixed at both ends:
对于有限区间上的边值问题，分离变量法是标准方法。考虑一根长度为L、两端固定的弦：
u(0,t) = u(L,t) = 0
We assume a solution of the form u(x,t) = X(x)T(t). Substituting into the wave equation gives:
我们假设解具有形式u(x,t) = X(x)T(t)。代入波动方程得到：
X”/X = (1/c²)(T”/T) = -λ
where λ is a separation constant. The boundary conditions imply X(0) = X(L) = 0. Solving the spatial equation X” + λX = 0 gives eigenvalues:
其中λ是分离常数。边界条件意味着X(0) = X(L) = 0。求解空间方程X” + λX = 0得到特征值：
λₙ = (nπ/L)², Xₙ(x) = sin(nπx/L), n = 1, 2, 3, …
The temporal equation T” + c²λT = 0 has solutions:
时间方程T” + c²λT = 0的解为：
Tₙ(t) = Aₙ cos(cnπt/L) + Bₙ sin(cnπt/L)
The general solution is a superposition of normal modes:
通解是简正模式的叠加：
u(x,t) = Σₙ₌₁^∞ [Aₙ cos(cnπt/L) + Bₙ sin(cnπt/L)] sin(nπx/L)
6. Fourier Coefficients | 傅里叶系数
The coefficients Aₙ and Bₙ are determined by the initial conditions. For initial displacement u(x,0) = φ(x):
系数Aₙ和Bₙ由初始条件确定。对于初始位移u(x,0) = φ(x)：
Aₙ = (2/L) ∫₀ᴸ φ(x) sin(nπx/L) dx
For initial velocity ∂u/∂t(x,0) = ψ(x):
对于初始速度∂u/∂t(x,0) = ψ(x)：
Bₙ = (2/(cnπ)) ∫₀ᴸ ψ(x) sin(nπx/L) dx
These are exactly the Fourier sine series coefficients. In practice, one computes these integrals for the given initial functions.
这些正是傅里叶正弦级数的系数。实际中，我们针对给定的初始函数计算这些积分。
Example: Suppose L = 1, φ(x) = x(1 – x), and ψ(x) = 0. Then Aₙ = (2/L)∫₀¹ x(1 – x) sin(nπx) dx = 4(1 – cos(nπ))/(n³π³) = 8/(n³π³) for n odd, 0 for n even. The solution is u(x,t) = Σₙ₌₁,₃,₅… [8/(n³π³)] cos(nπt) sin(nπx).
例如：设L = 1，φ(x) = x(1 – x)，ψ(x) = 0。那么Aₙ = (2/L)∫₀¹ x(1 – x) sin(nπx) dx = 4(1 – cos(nπ))/(n³π³)，当n为奇数时为8/(n³π³)，n为偶数时为0。解为u(x,t) = Σₙ₌₁,₃,₅… [8/(n³π³)] cos(nπt) sin(nπx)。
7. The Wave Equation in Two and Three Dimensions | 二维和三维波动方程
In higher dimensions, solutions become more complex. For a circular membrane (drum), we use polar coordinates. The Laplacian in polar coordinates is:
在更高维度中，解变得更加复杂。对于圆形膜（鼓），我们使用极坐标。极坐标中的拉普拉斯算子为：
∇²u = (1/r) ∂/∂r (r ∂u/∂r) + (1/r²) ∂²u/∂θ²
Applying separation of variables in the form u(r,θ,t) = R(r)Θ(θ)T(t) leads to Bessel’s equation for the radial part:
采用分离变量法u(r,θ,t) = R(r)Θ(θ)T(t)后，径向部分满足贝塞尔方程：
r² R” + r R’ + (k²r² – m²)R = 0
The solutions are Bessel functions Jₘ(kr). The allowed values of k are determined by the boundary condition R(a) = 0 for a drum of radius a.
其解为贝塞尔函数Jₘ(kr)。k的允许值由边界条件R(a) = 0（半径为a的鼓）确定。
In three dimensions, spherical coordinates lead to spherical harmonics and radial equations involving Bessel functions. These topics are typically covered in advanced courses and are essential for understanding electromagnetic waves and quantum mechanics.
在三维空间中，球坐标引入球谐函数和涉及贝塞尔函数的径向方程。这些主题通常在高级课程中涉及，对于理解电磁波和量子力学至关重要。
8. Energy Conservation | 能量守恒
An important property of the wave equation is energy conservation. For the one-dimensional wave equation with appropriate boundary conditions, define the total energy:
波动方程的一个重要性质是能量守恒。对于具有适当边界条件的一维波动方程，定义总能量：
E(t) = ½ ∫₀ᴸ [(∂u/∂t)² + c² (∂u/∂x)²] dx
Taking the derivative with respect to t and using integration by parts, we find dE/dt = 0, provided the boundary terms vanish. Thus, energy is conserved.
对t求导并使用分部积分，我们发现dE/dt = 0，前提是边界项消失。因此，能量守恒。
This energy conservation principle is used in many exams to check the validity of solutions or to derive qualitative behavior without solving the equation explicitly.
这个能量守恒原理在许多考试中用于检查解的合理性，或在不解方程的情况下推导定性行为。
9. Boundary Conditions and their Physical Meaning | 边界条件及其物理意义
There are several types of boundary conditions commonly encountered:
有几种常见的边界条件类型：
- Dirichlet condition: u(0,t) = 0. This fixes the displacement at the boundary, e.g., a fixed end of a string.
- 诺伊曼条件： ∂u/∂x(0,t) =
Published by TutorHao | Mathematics Revision Series | aleveler.com
更多咨询请联系16621398022（同微信）
屏轩国际教育cambridge primary/secondary checkpoint, cat4, ukiset,ukcat,igcse,alevel,PAT,STEP,MAT, ibdp,ap,ssat,sat,sat2课程辅导，国外大学本科硕士研究生博士课程论文辅导