平面电磁波在介质中的传播与波动方程推导
当人们谈论无线信号穿透墙壁、光在光纤中传输，或雷达探测远距离目标时，其背后统一的物理图景正是——电磁波在介质中的传播。这一现象的数学根基，并非来自某种经验公式，而是深植于一百多年前麦克斯韦建立的那组简洁而深刻的方程之中。
从这些方程出发，我们不仅能理解电磁波为何存在，还能精确描述它如何在不同材料中行进、衰减甚至发生折射与反射。本文将沿着经典电磁理论的路径，逐步揭示：平面电磁波是如何从麦克斯韦方程组中自然涌现的？损耗性介质又如何影响其传播特性？
无源介质中的基本方程体系
研究电磁波传播的第一步，是明确其所处环境的基本假设。考虑一个线性、均匀且各向同性的电介质，其宏观响应由三个参数刻画：介电常数 $\varepsilon$、磁导率 $\mu$ 和电导率 $\sigma$。若该区域不含自由电荷（$\rho = 0$）和外加电流源，则麦克斯韦方程组在时域下简化为：
$$
\nabla \cdot \mathbf{E} = 0 \tag{1}
$$
$$
\nabla \cdot \mathbf{B} = 0 \tag{2}
$$
$$
\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} \tag{3}
$$
$$
\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t} \tag{4}
$$
结合本构关系 $\mathbf{D} = \varepsilon \mathbf{E}$、$\mathbf{B} = \mu \mathbf{H}$ 以及欧姆定律形式的传导电流 $\mathbf{J} = \sigma \mathbf{E}$，可将 (4) 式改写为：
$$
\nabla \times \mathbf{H} = \sigma \mathbf{E} + \varepsilon \frac{\partial \mathbf{E}}{\partial t} \tag{5}
$$
这组方程构成了分析电磁波行为的基础框架。特别地，当 $\sigma = 0$ 时，介质为理想绝缘体，能量无耗散；而一旦 $\sigma > 0$，则意味着介质具有导电性，电磁波在其中传播将伴随焦耳热损耗。
电场波动方程的系统推导
为了获得电场 $\mathbf{E}$ 的传播动力学，需对方程进行进一步操作。对法拉第定律 (3) 两边取旋度：
$$
\nabla \times (\nabla \times \mathbf{E}) = -\nabla \times \left( \frac{\partial \mathbf{B}}{\partial t} \right)
$$
利用矢量恒等式：
$$
\nabla \times (\nabla \times \mathbf{A}) = \nabla(\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A}
$$
并代入高斯定律 (1) 中 $\nabla \cdot \mathbf{E} = 0$，左边化为：
$$
\nabla \times (\nabla \times \mathbf{E}) = -\nabla^2 \mathbf{E}
$$
右边交换空间与时间导数顺序（假设场足够光滑），有：
$$
-\nabla \times \left( \frac{\partial \mathbf{B}}{\partial t} \right) = -\frac{\partial}{\partial t} (\nabla \times \mathbf{B}) = -\mu \frac{\partial}{\partial t} (\nabla \times \mathbf{H})
$$
再将安培-麦克斯韦方程 (5) 代入：
$$
-\mu \frac{\partial}{\partial t} \left( \sigma \mathbf{E} + \varepsilon \frac{\partial \mathbf{E}}{\partial t} \right) = -\mu \sigma \frac{\partial \mathbf{E}}{\partial t} - \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2}
$$
联立左右两侧结果，最终得到电场满足的偏微分方程：
$$
\nabla^2 \mathbf{E} = \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2} + \mu \sigma \frac{\partial \mathbf{E}}{\partial t} \tag{6}
$$
这个方程被称为广义电磁波方程，适用于包含能量损耗的一般介质。它清楚地表明：电场的时空演化不仅受二阶时间导数驱动（对应波动行为），还受到一阶时间导数项的影响（对应阻尼效应）。当 $\sigma = 0$ 时，该项消失，退化为标准波动方程：
$$
\nabla^2 \mathbf{E} = \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2} \tag{7}
$$
此时解代表以速度 $v = 1/\sqrt{\mu\varepsilon}$ 匀速传播而不衰减的波。
磁场的对称性波动方程
类似过程可用于推导磁场 $\mathbf{H}$ 的传播规律。从修正后的安培定律 (5) 出发，取旋度：
$$
\nabla \times (\nabla \times \mathbf{H}) = \nabla \times \left( \sigma \mathbf{E} + \varepsilon \frac{\partial \mathbf{E}}{\partial t} \right)
$$
左边仍用矢量恒等式，并利用 $\nabla \cdot \mathbf{H} = 0$（因 $\mathbf{B} = \mu \mathbf{H}, \nabla \cdot \mathbf{B} = 0$），得：
$$
\nabla \times (\nabla \times \mathbf{H}) = -\nabla^2 \mathbf{H}
$$
右边展开为：
$$
\nabla \times (\sigma \mathbf{E}) + \varepsilon \frac{\partial}{\partial t} (\nabla \times \mathbf{E}) = \sigma (\nabla \times \mathbf{E}) + \varepsilon \frac{\partial}{\partial t} \left( -\frac{\partial \mathbf{B}}{\partial t} \right)
= \sigma (-\mu \frac{\partial \mathbf{H}}{\partial t}) - \varepsilon \mu \frac{\partial^2 \mathbf{H}}{\partial t^2}
$$
因此：
$$
-\nabla^2 \mathbf{H} = -\mu \sigma \frac{\partial \mathbf{H}}{\partial t} - \mu \varepsilon \frac{\partial^2 \mathbf{H}}{\partial t^2}
$$
整理后即得：
$$
\nabla^2 \mathbf{H} = \mu \varepsilon \frac{\partial^2 \mathbf{H}}{\partial t^2} + \mu \sigma \frac{\partial \mathbf{H}}{\partial t} \tag{8}
$$
可见，磁场也遵循与电场完全相同的波动方程结构。两者互为激发源，在空间中共振前行，形成自洽的横电磁波（TEM）模式。
自由空间中的极限情形
在真空中，所有材料参数取其基准值：$\varepsilon = \varepsilon_0$, $\mu = \mu_0$, $\sigma = 0$。此时波动方程进一步简化为：
$$
\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}, \quad
\nabla^2 \mathbf{H} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{H}}{\partial t^2}
$$
定义真空光速：
$$
c = \frac{1}{\sqrt{\mu_0 \varepsilon_0}} \approx 3 \times 10^8~\mathrm{m/s}
$$
于是方程可写作更紧凑的形式：
$$
\nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = 0
$$
这不仅是电磁波存在的直接证据，也是爱因斯坦狭义相对论中“光速不变原理”的理论起点。现代通信技术，无论是卫星链路还是移动网络，本质上都是对此类方程解的应用与操控。
单频激励下的复数表示：亥姆霍兹方程
工程实践中，多数系统工作于稳态正弦激励下。设电场随时间作简谐变化：
$$
\mathbf{E}(\mathbf{r}, t) = \mathrm{Re} \left{ \tilde{\mathbf{E}}(\mathbf{r}) e^{j\omega t} \right}
$$
其中 $\tilde{\mathbf{E}}(\mathbf{r})$ 是空间相关的复相量，$\omega$ 为角频率。
将其代入含损耗的波动方程 (6)，注意：
$$
\frac{\partial \mathbf{E}}{\partial t} \to j\omega \tilde{\mathbf{E}} e^{j\omega t}, \quad
\frac{\partial^2 \mathbf{E}}{\partial t^2} \to -\omega^2 \tilde{\mathbf{E}} e^{j\omega t}
$$
消去公共因子 $e^{j\omega t}$ 后，空间部分满足：
$$
\nabla^2 \tilde{\mathbf{E}} + \omega^2 \mu \varepsilon \left(1 + \frac{j\sigma}{\omega \varepsilon}\right) \tilde{\mathbf{E}} = 0
$$
引入复介电常数：
$$
\tilde{\varepsilon} = \varepsilon \left(1 + \frac{j\sigma}{\omega \varepsilon}\right) = \varepsilon - j\frac{\sigma}{\omega}
$$
则上式变为：
$$
\nabla^2 \tilde{\mathbf{E}} + \omega^2 \mu \tilde{\varepsilon} \tilde{\mathbf{E}} = 0 \tag{9}
$$
这就是著名的亥姆霍兹方程，广泛应用于天线设计、波导分析、光学成像等领域。它的解直接给出电磁场的空间分布特征，尤其适合处理单频辐射与散射问题。
平面波解及其横向传播特性
考虑在无限大均匀介质中沿 $+z$ 方向传播的平面波，其电场可设为：
$$
\mathbf{E}(z,t) = \mathbf{E}_0 \cos(\omega t - k z + \phi)
$$
其中波数 $k = \omega \sqrt{\mu \varepsilon}$。在复域中写作：
$$
\tilde{\mathbf{E}}(z) = \mathbf{E}_0 e^{-j k z}, \quad k = \omega \sqrt{\mu \tilde{\varepsilon}}
$$
代入亥姆霍兹方程 (9)，验证：
$$
\frac{d^2}{dz^2} \tilde{\mathbf{E}} + k^2 \tilde{\mathbf{E}} = (-k^2 + k^2)\tilde{\mathbf{E}} = 0
$$
成立。这说明该函数确实是方程的一个特解。
更重要的是，电场方向垂直于传播方向（$z$ 轴），表现出典型的横波性质。由法拉第定律可推出对应的磁场：
$$
\mathbf{H} = \frac{1}{\eta} \hat{k} \times \mathbf{E}
$$
其中 $\eta = \sqrt{\mu / \tilde{\varepsilon}}$ 称为介质的本征阻抗，决定了电场与磁场之间的幅值比和相位关系。
对于非磁性材料（$\mu \approx \mu_0$），$\eta$ 主要由介电响应决定。例如在空气中约为 $377~\Omega$，而在水中由于高介电常数和一定电导率，阻抗显著降低，导致电磁波难以深入穿透。
复折射率与振幅衰减机制
在光学和高频电磁学中，常用折射率来表征介质对波的响应：
$$
n = \frac{c}{v_p} = \sqrt{\frac{\varepsilon_r \mu_r}{1}} \approx \sqrt{\varepsilon_r} \quad (\text{当 } \mu_r \approx 1)
$$
但在有损耗的情况下，必须推广至复折射率：
$$
\tilde{n} = n - j\kappa
$$
其中实部 $n$ 控制相位传播速度，虚部 $\kappa$（消光系数）描述振幅衰减。
其与复介电常数的关系为：
$$
\tilde{n}^2 = \varepsilon_r \left(1 + j \frac{\sigma}{\omega \varepsilon_0 \varepsilon_r}\right)
$$
令 $\varepsilon’ = \varepsilon$, $\varepsilon’’ = \sigma / \omega$，则可通过分离实虚部求得：
$$
n = \sqrt{ \frac{|\tilde{\varepsilon}| + \varepsilon’}{2} }, \quad
\kappa = \sqrt{ \frac{|\tilde{\varepsilon}| - \varepsilon’}{2} }
$$
电场在传播过程中呈指数衰减：
$$
|\mathbf{E}(z)| \propto e^{-\alpha z}, \quad \alpha = \omega \sqrt{\frac{\mu \varepsilon}{2}} \left[ \sqrt{1 + \left(\frac{\sigma}{\omega \varepsilon}\right)^2} - 1 \right]^{1/2}
$$
当满足低损耗条件 $\sigma \ll \omega \varepsilon$ 时，近似有：
$$
\alpha \approx \frac{\sigma}{2} \sqrt{\frac{\mu}{\varepsilon}}
$$
这意味着，即使微弱的电导率也会导致信号随距离迅速减弱，这也是为什么海水对无线电波几乎是“不透明”的原因。
边界行为与横波结构的保持条件
在两种介质交界面，电磁场必须满足连续性边界条件：
- 电场切向分量连续：$E_{t1} = E_{t2}$
- 磁场切向分量连续（无表面电流）：$H_{t1} = H_{t2}$
- 电位移法向分量连续（无自由面电荷）：$D_{n1} = D_{n2}$
- 磁感应强度法向分量连续：$B_{n1} = B_{n2}$
这些条件共同决定了电磁波在界面处的反射与透射行为，是菲涅耳公式的推导基础。
在均匀各向同性介质中，平面波的电场 $\mathbf{E}$、磁场 $\mathbf{H}$ 与波矢 $\mathbf{k}$ 相互垂直，构成右手正交系：
$$
\mathbf{E} \perp \mathbf{H}, \quad \mathbf{E} \perp \mathbf{k}, \quad \mathbf{H} \perp \mathbf{k}
$$
坡印廷矢量 $\mathbf{S} = \mathbf{E} \times \mathbf{H}$ 指示能量流动方向，通常与 $\mathbf{k}$ 一致。
然而，若介质为各向异性（如石英晶体或铁氧体），$\mathbf{D}$ 与 $\mathbf{E}$ 不再平行，导致 $\mathbf{E}$ 不一定垂直于 $\mathbf{k}$，出现所谓的“非常规波”（extraordinary wave），此时波矢方向与能流方向发生分离，带来复杂的双折射效应。
| 物理情形 | 对应波动方程 |
|---|---|
| 一般耗损介质 | $\nabla^2 \mathbf{E} = \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2} + \mu \sigma \frac{\partial \mathbf{E}}{\partial t}$ |
| 无损介质 | $\nabla^2 \mathbf{E} = \mu \varepsilon \frac{\partial^2 \mathbf{E}}{\partial t^2}$ |
| 自由空间 | $\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}$ |
| 时谐场（复数域） | $\nabla^2 \tilde{\mathbf{E}} + \omega^2 \mu \tilde{\varepsilon} \tilde{\mathbf{E}} = 0$ |
| 复折射率定义 | $\tilde{n}^2 = \varepsilon_r (1 + j \frac{\sigma}{\omega \varepsilon})$ |
这套从麦克斯韦方程组出发的完整推导链条，展示了电磁波作为场的动力学实体，如何在数学上被严格构造出来。无论是在光纤中穿行的光脉冲，还是穿过人体组织的毫米波信号，其本质都可以归结为此处所呈现的波动模型。
这种基于第一性原理的分析方法，不仅提供了预测能力，更为新材料设计、高性能天线开发以及复杂电磁环境建模奠定了坚实的理论基础。
