电磁波基础:麦克斯韦方程简介、波动方程、平面波与球面波
阵列天线的工作本质是控制电磁波在空间中的叠加与干涉。要理解阵列如何形成波束、如何扫描、如何出现栅瓣,首先必须清楚单个辐射单元向外辐射的到底是什么——是电磁波,以及电磁波在自由空间中如何传播。本节从麦克斯韦方程组出发,推导出波动方程,再引入两种最基础的波前形态:平面波与球面波。它们是后续分析阵列方向图、相位差与口径分布的数学基石。
核心原理:麦克斯韦方程组与波动方程的由来
电磁场的时空演化由四个方程统一描述。在无源、线性、各向同性的均匀介质中,麦克斯韦方程组可写为:
∇ × E = -μ ∂H/∂t (法拉第定律)
∇ × H = ε ∂E/∂t + J (安培-麦克斯韦定律,无源时 J=0)
∇ · E = ρ/ε (高斯定律,无源时 ρ=0)
∇ · H = 0 (磁通连续性)
对第一个方程两边取旋度,并利用矢量恒等式 ∇×(∇×E) = ∇(∇·E) - ∇²E,在无源区域(ρ=0, J=0)中消去 H,即可得到关于电场 E 的齐次波动方程:
∇²E - με ∂²E/∂t² = 0
同理,磁场 H 也满足相同形式的方程。其中 μ 为磁导率,ε 为介电常数。在自由空间中,μ₀ = 4π×10⁻⁷ H/m,ε₀ ≈ 8.854×10⁻¹² F/m,由此得到电磁波在真空中的传播速度:
c = 1/√(μ₀ε₀) ≈ 2.998×10⁸ m/s
波动方程的意义在于:它揭示了电磁场可以脱离源而以波的形式独立传播,且传播速度仅由介质参数决定。阵列天线设计中最关心的频率 f 与波长 λ 的关系,正是由波动方程的解所决定:λ = c/f。在介质中,波长会缩短为 λ = c/(f√εᵣ),其中 εᵣ 为相对介电常数——这一点在微带贴片天线设计中尤为关键。
关键推论:波动方程是线性的,因此多个源产生的电磁场可以线性叠加。这正是阵列天线“多个单元合成波束”的数学基础。若方程是非线性的,阵列方向图合成将不复存在。
平面波:远场分析的理想化工具
波动方程最简单的解是平面波。设波沿 +z 方向传播,电场仅沿 x 方向极化,则解可写为:
E_x(z,t) = E₀ cos(ωt - kz + φ₀)
其中 ω = 2πf 为角频率,k = 2π/λ 为波数,φ₀ 为初始相位。平面波的核心特征是:在任意垂直于传播方向的平面上,场强幅度和相位处处相同。也就是说,等相位面是无限大的平面,且等相位面垂直于传播方向。
平面波是理想化的模型——现实中不存在无限大的等相位面。但在阵列天线的远场区域(通常取 r > 2D²/λ,D 为天线最大尺寸),观察点距离足够远时,球面波的局部波前可以近似为平面波。这一近似极大简化了阵列方向图的计算:各单元到远场观察点的路径差可以用简单的几何投影关系表达,从而直接推导出阵列因子。
工程直觉:平面波没有距离衰减(幅度恒定),因此它描述的是“已经传播到无穷远”的波。在阵列分析中,我们关心的是远场方向图,所以将观察点视为处于平面波区域是合理的。
球面波:单点源的真实辐射形态
一个理想的点源(无穷小偶极子或各向同性辐射元)向空间均匀辐射时,其等相位面是以源点为球心的同心球面。在球坐标系中,波动方程的解为:
E(r,t) = (E₀/r) cos(ωt - kr + φ₀)
球面波有两个显著区别于平面波的特征:
- 幅度随距离衰减:场强与 r 成反比,功率密度与 r² 成反比(满足能量守恒)。
- 等相位面为球面:在距离源 r 的球面上,各点相位相同。
阵列天线的每个单元在近场区域辐射的都是球面波。当多个单元的球面波在空间叠加时,由于各单元到观察点的距离不同,会产生路径差 Δr,进而引起相位差 Δφ = k·Δr = 2πΔr/λ。阵列波束形成的本质,就是通过控制各单元的激励相位来补偿或利用这些路径差,使目标方向上的场同相叠加,其他方向上异相抵消。
常见误区:初学者常将“球面波”与“各向同性辐射”混淆。球面波描述的是波前形状(等相位面为球面),而各向同性辐射描述的是幅度方向图(所有方向辐射强度相同)。一个实际天线的球面波在近场成立,但其幅度方向图并非均匀——例如偶极子在端射方向辐射为零。
平面波与球面波的对比与工程意义
| 特征 | 平面波 | 球面波 |
|---|---|---|
| 等相位面 | 无限大平面 | 同心球面 |
| 幅度变化 | 恒定,无衰减 | 与 r 成反比衰减 |
| 源模型 | 无穷远点源 | 有限距离点源 |
| 适用场景 | 远场方向图分析、阵列因子计算 | 近场辐射、单元互耦分析 |
| 相位差来源 | 仅由传播方向投影决定 | 由球面几何路径差决定 |
| 阵列分析中的角色 | 观察点处的入射波近似 | 各单元实际辐射形态 |
在阵列天线设计中,两种波的角色分工明确:单元设计阶段关注球面波(近场互耦、阻抗匹配),阵列综合阶段采用平面波近似(方向图乘积定理)。理解从球面波到平面波的过渡条件(即远场判据 r > 2D²/λ),是正确使用方向图乘积定理的前提。
案例:从球面波到平面波的过渡距离计算
设某阵列天线最大尺寸 D = 0.5 m,工作频率 f = 10 GHz(λ = 0.03 m)。则远场起始距离为:
r_far = 2D²/λ = 2 × (0.5)² / 0.03 ≈ 16.7 m
这意味着在 16.7 m 以内,观察点处的波前仍具有明显的球面弯曲,各单元到观察点的路径差不能简单用平面波投影计算。若在此距离内测量方向图,测得的波束宽度会比理论值偏宽,副瓣电平也会偏离设计值。实际测试中,若暗室尺寸无法满足远场条件,则需要使用紧缩场(反射面将球面波校准为平面波)或近场扫描测量。
设计启示:对于 77 GHz 车载雷达(λ ≈ 3.9 mm),若阵列尺寸 D = 0.1 m,远场距离仅约 5.1 m,很容易满足。但对于低频段(如 300 MHz,λ = 1 m)的大型阵列(D = 10 m),远场距离高达 200 m,几乎无法在室内直接测量——此时必须依赖近场测量或外场测试。
误区与复盘:平面波近似失效的场景
平面波近似虽然强大,但在以下场景中会失效,需要特别警惕:
- 近场测量:探头距阵列太近时,波前弯曲明显,直接套用平面波假设会导致方向图畸变。
- 大型阵列的近距离观察:即使满足 r > 2D²/λ,若阵列本身存在幅度锥削,近场区的幅度分布与远场方向图之间并非简单傅里叶变换关系。
- 相控阵波束扫描时的近场效应:大角度扫描时,阵列口径投影缩小,等效 D 减小,远场距离变短,但若仍按法向远场距离计算,可能低估近场效应。
- 互耦分析:单元之间的互耦发生在极近的距离(亚波长量级),此时必须使用完整的球面波展开或全波数值方法,平面波假设完全不适用。
复盘上述内容,核心脉络可归纳为:麦克斯韦方程组 → 波动方程 → 平面波与球面波解 → 远场近似条件 → 阵列分析中的角色分工。下一节将基于平面波模型,推导二元阵的干涉公式,正式进入阵列方向图合成的数学表达。