四、阵列天线理论：均匀线阵、平面阵、共形阵的数学模型与波束形成

阵列天线通过多个辐射单元的相干叠加，在空间特定方向形成高增益波束，同时抑制干扰和杂波。本章从均匀线阵（ULA）出发，建立阵列因子与方向图乘积定理的数学框架，进而推广至平面阵与共形阵，最后给出波束形成的基本公式与工程实现路径。理解这些模型是掌握相控阵波束扫描、低副瓣设计、自适应置零等高级技术的基础。

4.1 均匀线阵的数学模型与阵列因子

4.1.1 核心原理：从阵元到阵列的叠加

均匀线阵是最基础的阵列构型，N 个各向同性阵元沿直线等间距排列，间距为 d。设第 n 个阵元的复加权为 wₙ = Aₙ·e^{j·φₙ}，入射平面波方向与阵列法线夹角为 θ，则相邻阵元间的空间相位差为 Δψ = k·d·sinθ（以法线为参考）或 k·d·cosθ（以轴线为参考）。阵列远区辐射场为各阵元贡献的相干叠加：

E_total(θ) = Σ_{n=0}^{N-1} wₙ · e^{j·n·k·d·sinθ} · F_element(θ)

其中 F_element(θ) 为单元方向图。当各阵元方向图一致时，总方向图可分解为单元因子与阵列因子的乘积，即方向图乘积定理。阵列因子仅取决于阵元位置、加权系数和观察方向，与单元形式无关，这为独立设计阵列波束提供了便利。

4.1.2 阵列因子计算与波束指向

对于等幅激励（wₙ = 1），均匀线阵的阵列因子可写成等比级数求和形式：

AF(θ) = Σ_{n=0}^{N-1} e^{j·n·k·d·sinθ} = sin(N·ψ/2) / sin(ψ/2) · e^{j·(N-1)·ψ/2}

其中 ψ = k·d·sinθ。当 ψ = 2πm（m 为整数）时，AF 取得最大值 N。若引入线性相位递变 φₙ = -n·k·d·sinθ₀，则波束最大值指向 θ₀，实现电扫描。工程中常用相移量 Δφ = k·d·sinθ₀ 来设定波束指向，当 θ₀ 从 -60° 扫至 +60° 时，所需相移范围随 d/λ 增大而增大。

4.1.3 案例拆解：8 元线阵的波束扫描设计

某 X 波段（10 GHz，λ=30 mm）雷达采用 8 元均匀线阵，d=15 mm（0.5λ）。要求波束在 ±45° 范围内扫描，副瓣电平低于 -20 dB。设计步骤如下：

1. 确定相移量：对 θ₀=30°，Δφ = k·d·sin30° = (2π/30)·15·0.5 = π/2，即相邻阵元相位差 90°。
2. 选择幅度加权：采用 -35 dB 泰勒加权，查表得幅度分布为 [0.28, 0.62, 0.92, 1.0, 1.0, 0.92, 0.62, 0.28]。
3. 计算阵列因子：代入公式得主瓣宽度约 12.8°，第一副瓣 -35 dB。
4. 验证栅瓣条件：d=0.5λ 时，栅瓣出现在 |sinθ|≥2 处，扫描范围内无栅瓣。

实测方向图与仿真吻合良好，但在大扫描角（45°）时主瓣增益下降约 1.8 dB，这是由阵元投影面积减小所致，可通过单元方向图补偿。

4.2 平面阵的二维波束形成

4.2.1 核心原理：矩形栅格平面阵

平面阵将阵元分布在二维平面上，可同时控制方位角和俯仰角。以 M×N 矩形栅格为例，阵元位于 (m·dₓ, n·d_y) 处，m=0~M-1，n=0~N-1。观察方向用方向余弦 (u, v) 表示：u = sinθ·cosφ，v = sinθ·sinφ。阵列因子为：

AF(u,v) = Σ_{m=0}^{M-1} Σ_{n=0}^{N-1} wₘₙ · e^{j·k·(m·dₓ·u + n·d_y·v)}

当加权可分离时，即 wₘₙ = wₘ·wₙ，阵列因子可分解为两个一维阵列因子的乘积：AF(u,v) = AFₓ(u)·AF_y(v)，这大大简化了设计。但实际中为获得圆形等高性能波束，常采用不可分离的二维加权。

4.2.2 波束指向与栅瓣抑制

设期望波束指向 (θ₀, φ₀)，则所需相移为：

φₘₙ = -k·(m·dₓ·u₀ + n·d_y·v₀)

其中 u₀ = sinθ₀·cosφ₀，v₀ = sinθ₀·sinφ₀。为避免栅瓣，需满足 dₓ/λ ≤ 1/(1+|u₀|) 且 d_y/λ ≤ 1/(1+|v₀|)。对于全空域扫描（|u₀|≤1），要求 dₓ, d_y ≤ 0.5λ。实际系统常取 0.5λ~0.55λ 以平衡成本与性能。

4.2.3 案例拆解：8×8 平面阵的波束赋形

某相控阵通信终端采用 8×8 平面阵，dₓ=d_y=0.5λ，工作频率 28 GHz。要求波束指向 (θ₀=30°, φ₀=45°)，副瓣低于 -25 dB。计算过程：

• 计算方向余弦：u₀ = sin30°·cos45° = 0.3536，v₀ = sin30°·sin45° = 0.3536
• 计算相移：φₘₙ = -2π·(m·0.5·0.3536 + n·0.5·0.3536) = -π·0.3536·(m+n)
• 采用 -30 dB 二维泰勒加权，生成 8×8 幅度矩阵
• 仿真验证：主瓣宽度约 12.5°，第一副瓣 -27.3 dB，指向误差 <0.1°

实测中由于阵元互耦，副瓣抬升约 2 dB，通过加装金属隔离墙和优化馈电网络后达标。

4.3 共形阵的建模与计算

4.3.1 核心原理：曲面阵列的通用模型

共形阵将阵元共形于载体表面（如导弹弹头、机翼前缘、卫星蒙皮），具有低雷达散射截面（RCS）和宽角覆盖优势。其数学模型需考虑每个阵元的空间位置矢量 rₙ 和单元方向图 Fₙ(θ,φ)：

AF(θ,φ) = Σ_{n=1}^{N} wₙ · Fₙ(θ,φ) · e^{j·k·(rₙ·â)}

其中 â 为观察方向的单位矢量，rₙ·â 为阵元到观察点的空间相位差。由于阵元朝向不同，Fₙ 各不相同，方向图乘积定理不再适用，必须逐元计算。

4.3.2 圆柱共形阵的阵列因子

以圆柱阵为例，阵元分布在半径为 R 的圆柱面上，第 n 个阵元位置为 (R·cosφₙ, R·sinφₙ, zₙ)。观察方向 (θ, φ) 的单位矢量为 â = (sinθ·cosφ, sinθ·sinφ, cosθ)。则空间相位为：

rₙ·â = R·sinθ·cos(φ-φₙ) + zₙ·cosθ

阵列因子为：

AF(θ,φ) = Σ wₙ · Fₙ(θ,φ) · e^{j·k·[R·sinθ·cos(φ-φₙ) + zₙ·cosθ]}

由于阵元方向图 Fₙ 随 φₙ 变化，即使等幅激励，阵列因子也无法写成闭式解析式，需数值计算。

4.3.3 案例拆解：圆柱阵的宽角扫描

某导引头采用 32 元圆柱共形阵，R=0.8λ，阵元沿圆周均匀分布（每 11.25° 一个），分 4 层，层间距 0.5λ。要求波束在方位 360°、俯仰 ±60° 内扫描。设计要点：

1. 单元选择：采用微带贴片，单元波束宽度约 120°，覆盖相邻阵元间隔。
2. 加权计算：对每个扫描角，计算各阵元的空间相位并取共轭，得到最大比合成加权。
3. 方向图补偿：由于阵元朝向差异，需对每个阵元方向图进行归一化补偿，否则扫描增益波动达 ±3 dB。
4. 互耦修正：实测互耦矩阵 S，用有源方向图替代孤立单元方向图，修正后副瓣改善 4~6 dB。

实测结果：方位全向扫描增益波动 ±1.2 dB，俯仰 60° 处增益下降 2.5 dB，满足导引头指标。

4.4 波束形成的基本公式与工程实现

4.4.1 核心原理：从阵列因子到波束形成器

波束形成的本质是确定加权矢量 w，使阵列输出在期望方向增益最大、干扰方向增益最小。基本公式为：

y(t) = wᴴ · x(t)

其中 x(t) 为阵列接收信号矢量，w 为复加权矢量。常规波束形成（CBF）取 w = a(θ₀)，a(θ₀) 为期望方向的导向矢量。最小方差无失真响应（MVDR）波束形成器为：

w_MVDR = R⁻¹·a(θ₀) / (aᴴ(θ₀)·R⁻¹·a(θ₀))

其中 R 为接收数据的协方差矩阵。MVDR 可在干扰方向自动形成零陷，但需要估计 R 且对阵列误差敏感。

4.4.2 数字波束形成（DBF）架构

现代相控阵普遍采用 DBF，每个阵元通道包含 LNA、下变频、ADC，数字域完成加权求和。典型流程：

1. 阵元接收 → 模拟下变频 → ADC 采样
2. 数字下变频 (DDC) → I/Q 数据
3. 计算导向矢量 a(θ₀,φ₀)
4. 生成加权 w = a ⊙ w_win（或自适应算法）
5. 加权求和 → 波束输出

DBF 的优势在于可同时形成多个波束、自适应置零和校准补偿。但 ADC 的采样率、位数和通道一致性直接影响系统性能。

4.4.3 案例拆解：自适应波束形成抑制干扰

某 16 元线阵接收信号，期望信号来自 10°，两个干扰分别来自 -30° 和 40°，干噪比 40 dB。采用 MVDR 算法：

• 采集 1000 个快拍，估计协方差矩阵 R̂ = (1/K)·Σ x·xᴴ
• 计算加权 w_MVDR，在 -30° 和 40° 形成 -45 dB 零陷
• 输出 SINR 从 CBF 的 -5 dB 提升至 22 dB
• 对角加载（加载量 10⁻³·trace(R)）提高稳健性，零陷深度变为 -38 dB，但输出 SINR 仅下降 0.8 dB

实际系统中，通道幅相误差会导致零陷抬升，需在算法前进行通道校准，校准精度要求幅度 <0.1 dB、相位 <1°。

4.4.4 总结复盘

本章从均匀线阵的解析模型出发，建立了阵列因子与方向图乘积定理的完整框架，并推广至平面阵和共形阵。核心要点归纳如下：

波束形成的基本公式 y = wᴴ·x 统一了常规和自适应方法，工程实现需综合考虑阵元互耦、通道误差、量化效应和计算复杂度。建议在系统设计初期就建立包含单元方向图和互耦的完整阵列模型，避免后期返工。