3、阵列天线原理：阵列因子、单元因子、方向图乘积定理、均匀线性阵列（ULA）的波束形成
阵列天线并非简单地将多个辐射单元“捆”在一起，其核心价值在于通过控制各单元信号的幅度与相位，在空间中形成可操控的干涉图样。理解这一过程需要从“单元因子”与“阵列因子”的分离入手，再通过方向图乘积定理将二者统一。本章将系统梳理阵列天线的数学基础，重点剖析均匀线性阵列（ULA）的波束形成机理，并结合工程案例说明如何从方向图指标反推阵列参数。
3.1 核心原理：从单元辐射到阵列合成的数学解耦
阵列天线的分析之所以能够简化，关键在于将总方向图分解为两个独立因子的乘积。这一分解的前提是：阵列中所有单元具有相同的辐射方向图，且单元间的互耦可以忽略或已通过校准补偿。在此条件下，总电场方向图可表示为：
F(θ,φ) = EF(θ,φ) × AF(θ,φ)
其中，单元因子（Element Factor, EF）描述单个辐射单元在自由空间中的固有方向特性，它由单元的类型（如偶极子、贴片、喇叭）、尺寸、工作频率以及所在环境（如无限大地平面、有限接地板）决定。例如，一个沿z轴放置的半波偶极子，其单元因子为 EF(θ) = cos(π/2·cosθ)/sinθ，在H面（θ=90°平面）呈现全向性，在E面呈现“8”字形。
阵列因子（Array Factor, AF）则完全由阵列的几何布局和激励分布决定，与单元类型无关。对于N个各向同性点源组成的阵列，其阵列因子的一般形式为：
AF(θ,φ) = Σ_{n=0}^{N-1} A_n · exp(j·(k·r_n·cosψ_n + φ_n))式中，A_n 为第n个单元的激励幅度，r_n 为位置矢量，ψ_n 为观察方向与位置矢量之间的夹角，φ_n 为激励相位。这一公式揭示了阵列合成的本质：通过调整每个单元的幅度和相位，可以改变空间各方向的干涉叠加结果，从而“塑造”波束。
3.2 均匀线性阵列（ULA）的波束形成：从数学到物理
均匀线性阵列是最基础也最常用的阵列构型。设N个各向同性单元沿x轴等间距d排列，各单元激励幅度相同（A_n = 1），激励相位按等差级数递增（φ_n = n·Δφ）。当观察方向与阵列法线（broadside）夹角为θ时，相邻单元间的空间相位差为 k·d·sinθ，其中 k = 2π/λ 为波数。阵列因子可写成等比级数求和形式：
AF(θ) = Σ_{n=0}^{N-1} exp(j·n·(k·d·sinθ + Δφ))
= sin(N·ψ/2) / (N·sin(ψ/2))其中 ψ = k·d·sinθ + Δφ。当 ψ = 0 时，AF 取得最大值 N，对应主波束方向 θ₀ 满足：
sinθ₀ = -Δφ / (k·d)这一关系式是相控阵波束扫描的核心。通过连续改变相邻单元间的激励相位差 Δφ，即可实现主波束在空间中的连续扫描，而无需任何机械转动。例如，若 d = λ/2，要使波束指向 θ₀ = 30°，则所需相位差为 Δφ = -k·d·sin30° = -π/2，即相邻单元间相位依次滞后90°。
3.3 判断方法与操作流程：从指标到阵列参数设计
在实际工程中，设计一个ULA通常遵循以下流程：
- 确定工作频率与波长：由系统工作频段计算自由空间波长 λ = c/f。
- 选择单元间距 d：为避免栅瓣（grating lobes），需满足 d ≤ λ/(1+|sinθ_max|)，其中 θ_max 为最大扫描角。若要求扫描范围 ±60°，则 d ≤ λ/1.866 ≈ 0.536λ，工程上常取 d = λ/2 作为折中。
- 确定阵元数量 N：根据所需波束宽度（半功率波束宽度 HPBW）估算。对于均匀激励ULA，HPBW ≈ 0.886λ/(N·d·cosθ₀)（弧度）。若要求 θ₀=0° 时 HPBW ≤ 5°，则 N·d ≥ 0.886λ/sin5° ≈ 10.2λ，取 d=λ/2 时 N ≥ 21。
- 计算激励相位分布：按 Δφ = -k·d·sinθ₀ 计算各单元所需相位，通过移相器或数字波束形成（DBF）网络实现。
- 仿真验证与优化：使用全波仿真软件验证方向图、增益、副瓣电平（SLL）等指标，必要时引入幅度锥削（如泰勒分布、切比雪夫分布）以降低副瓣。
| 参数 | 符号 | 典型值/范围 | 对方向图的影响 |
|---|---|---|---|
| 阵元数量 | N | 8~64（常见） | N越大，波束越窄，增益越高 |
| 单元间距 | d | 0.3λ~0.7λ | d过大产生栅瓣，d过小互耦增强 |
| 激励幅度 | A_n | 均匀或锥削 | 锥削降低副瓣，但展宽主瓣 |
| 激励相位 | φ_n | 等差递增 | 控制主波束指向 |
| 工作频率 | f | 由系统决定 | 频率变化导致电尺寸变化 |
3.4 案例拆解：X波段8单元ULA波束扫描设计
背景与约束：某X波段（9.4GHz，λ≈31.9mm）雷达前端需实现方位面±45°电扫，要求波束宽度≤12°，副瓣电平≤-20dB，且阵列尺寸不超过250mm。
初始参数计算：取 d = 0.5λ = 15.95mm，则8单元阵列总长 L = 7d = 111.65mm，满足尺寸约束。均匀激励时，HPBW ≈ 0.886λ/(8×0.5λ) = 0.2215 rad ≈ 12.7°，略超指标。为压缩波束，将单元数增至10，L = 9×15.95 = 143.6mm，HPBW ≈ 10.2°，满足要求。
副瓣抑制设计：均匀激励的ULA第一副瓣电平约为-13.2dB，不满足-20dB要求。采用道尔夫-切比雪夫幅度分布，设定副瓣电平-25dB，计算得到各单元归一化幅度为：
单元序号: 0 1 2 3 4 5 6 7 8 9
幅度: 0.31 0.52 0.78 1.00 1.00 1.00 1.00 0.78 0.52 0.31相位计算：扫描至45°时，相邻单元相位差 Δφ = -k·d·sin45° = -2π×0.5×0.707 ≈ -2.221 rad ≈ -127.3°。各单元相位依次为 0°, -127.3°, -254.6°, ...，实际通过5位移相器量化（步进11.25°）实现。
仿真结果：全波仿真显示，在θ₀=45°时，主波束指向准确，HPBW≈11.8°，副瓣电平-23.5dB，增益较法向下降约3.2dB（由阵列投影面积减小引起）。
风险与代价：切比雪夫分布以展宽主瓣约15%为代价换取副瓣抑制；大扫描角时单元互耦导致阻抗失配，需在功分网络中加入隔离电阻或使用校准算法补偿。
3.5 常见误区与总结复盘
误区一：忽视栅瓣条件。当 d > λ/(1+|sinθ_max|) 时，除主瓣外在可见空间会出现额外的大辐射波瓣，造成能量泄漏和测向模糊。例如 d=0.7λ、扫描至60°时，栅瓣出现在 -90°方向。
误区二：将单元因子与阵列因子混为一谈。单元因子是“个体属性”，阵列因子是“群体行为”。微带贴片单元在H面方向图较宽，即使阵列因子在端射方向（end-fire）有峰值，总方向图也可能因单元因子在端射方向为零而被抑制。
误区三：忽略互耦对相位分布的影响。理论计算假设各单元独立工作，实际中相邻单元的耦合会改变有效激励相位，导致波束指向偏差。工程上需在调试阶段通过近场测量或内置校准网络修正。
误区四：认为增加阵元数总能提高增益。增益提升受限于阵列口径效率，当单元间距过小时，互耦加剧，效率下降，增益反而可能降低。
总结：阵列天线的方向图合成是“单元”与“阵列”两个层面共同作用的结果。方向图乘积定理提供了清晰的分解思路，而ULA的等差相位激励则揭示了波束扫描的本质——通过相位梯度控制等相位面的倾斜角度。掌握这些原理后，工程师便能从系统指标出发，快速完成阵列规模估算、馈电网络设计和性能验证。在实际项目中，务必结合全波仿真与实测数据迭代优化，才能实现理论到工程的可靠落地。