电磁波基础:麦克斯韦方程简介、波动方程、平面波与球面波

阵列天线的工作本质是控制电磁波在空间中的叠加与干涉。要理解阵列如何形成波束、如何扫描、如何出现栅瓣,首先必须清楚单个辐射单元向外辐射的到底是什么——是电磁波,以及电磁波在自由空间中如何传播。本节从麦克斯韦方程组出发,推导出波动方程,再引入两种最基础的波前形态:平面波与球面波。它们是后续分析阵列方向图、相位差与口径分布的数学基石。

核心原理:麦克斯韦方程组与波动方程的由来

电磁场的时空演化由四个方程统一描述。在无源、线性、各向同性的均匀介质中,麦克斯韦方程组可写为:

∇ × 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 减小,远场距离变短,但若仍按法向远场距离计算,可能低估近场效应。
  • 互耦分析:单元之间的互耦发生在极近的距离(亚波长量级),此时必须使用完整的球面波展开或全波数值方法,平面波假设完全不适用。

复盘上述内容,核心脉络可归纳为:麦克斯韦方程组 → 波动方程 → 平面波与球面波解 → 远场近似条件 → 阵列分析中的角色分工。下一节将基于平面波模型,推导二元阵的干涉公式,正式进入阵列方向图合成的数学表达。

电磁波传播形态与阵列分析的关系 麦克斯韦方程组 ∇×E = -μ∂H/∂t ∇×H = ε∂E/∂t 波动方程 ∇²E - με ∂²E/∂t² = 0 线性 → 叠加原理成立 平面波 E = E₀ cos(ωt - kz + φ₀) 等相位面:无限大平面 幅度恒定,无距离衰减 球面波 E = (E₀/r) cos(ωt - kr + φ₀) 等相位面:同心球面 幅度 ∝ 1/r,功率 ∝ 1/r² 源 等相位面 (球面) r 传播方向 z 等相位面 (平面) |E| 恒定 远场条件:r > 2D²/λ 球面波 → 局部近似为平面波 阵列方向图分析 路径差 → 相位差 → 波束形成 单元设计与互耦 近场球面波展开、全波仿真 阵列天线远场分析采用平面波近似,近场互耦分析必须使用球面波模型