4、波束形成概述:波束形成的基本概念、波束图、主瓣与旁瓣、阵列增益与分辨率
波束形成(Beamforming)是阵列信号处理的核心技术之一,其本质是通过对阵列各阵元接收信号进行加权、延时或相位补偿,使阵列在特定方向形成指向性响应,从而在空间域实现信号增强与干扰抑制。本章作为波束形成模块的起点,将系统阐述波束形成的基本概念、波束图(Beam Pattern)的物理含义、主瓣与旁瓣的工程权衡,以及阵列增益与分辨率这两个决定系统性能的关键指标。理解这些基础概念,是后续学习自适应波束形成、DOA估计算法(如MUSIC、ESPRIT)的前提。
4.1 波束形成的基本概念与数学表达
波束形成的核心思想源于对空间信号的相干叠加。考虑一个由 M 个全向阵元组成的均匀线阵(ULA),阵元间距为 d。当平面波以角度 θ 入射时,相邻阵元间的传播延时差为 τ = d·sinθ/c,对应相位差为 φ = 2πd·sinθ/λ。若对各阵元输出施加复权重 wm(m = 1, 2, ..., M),则阵列输出为:
y(t) = Σ w_m* · x_m(t) = wᴴ · x(t)
其中 w = [w₁, w₂, ..., w_M]ᵀ 为权向量,x(t) 为阵元接收信号向量,上标 H 表示共轭转置。通过设计权向量 w,可以使来自期望方向的信号同相叠加(增强),而来自其他方向的干扰异相抵消(抑制)。这一过程在频域上等价于空间滤波,在时域上则表现为对阵列孔径的利用。
从工程实现角度看,波束形成可分为三类:固定波束形成(权向量预先计算,不随环境变化)、自适应波束形成(根据接收数据实时调整权向量,如 LCMV、MVDR 算法)和数字波束形成(在基带通过 DSP/FPGA 实现,灵活性最高)。现代相控阵雷达和 5G 大规模 MIMO 系统普遍采用数字波束形成,其优势在于可以同时形成多个波束,且波束指向切换速度仅受限于基带处理时延。
4.2 波束图:主瓣、旁瓣与栅瓣
波束图(Beam Pattern)定义为阵列输出功率对入射角度的响应,数学上表示为 B(θ) = |wᴴ·a(θ)|²,其中 a(θ) 为方向向量(Steering Vector)。对于 ULA,a(θ) = [1, ejφ, ej2φ, ..., ej(M-1)φ]ᵀ。波束图是衡量阵列空间选择性的核心工具,其关键特征包括:
- 主瓣(Main Lobe):波束图的最大响应区域,其中心指向期望信号方向。主瓣宽度(通常用 3dB 波束宽度表示)决定了阵列的角度分辨能力。
- 旁瓣(Side Lobes):主瓣之外的局部极大值区域。旁瓣电平越高,来自非期望方向的干扰越容易泄漏进阵列输出,通常要求旁瓣低于主瓣 20~40dB。
- 栅瓣(Grating Lobes):当阵元间距 d > λ/2 时,在空间上出现的与主瓣等高的虚假响应。栅瓣会导致空间模糊,必须通过约束 d ≤ λ/2 来避免。
- 零陷(Nulls):波束图中的极小值点,可用于抑制特定方向的干扰。
均匀加权(w = [1,1,...,1]ᵀ)时,ULA 的波束图具有解析表达式:
B(θ) = |sin(Mπd·sinθ/λ) / (M·sin(πd·sinθ/λ))|²
该式表明,主瓣宽度与阵列孔径 L = (M-1)·d 成反比,阵元数越多、孔径越大,主瓣越窄,分辨率越高。但均匀加权的第一旁瓣电平约为 -13.3dB,在许多应用中无法满足要求,因此需要引入幅度加权来压低旁瓣,代价是主瓣略微展宽。
| 加权方式 | 第一旁瓣电平 | 3dB 主瓣宽度(相对均匀加权) | 适用场景 |
|---|---|---|---|
| 均匀加权 | -13.3 dB | 1.0× | 最大分辨率,旁瓣要求不高 |
| 汉明窗 | -42.8 dB | 1.47× | 雷达/声呐,需低旁瓣 |
| 泰勒窗(nbar=4) | -30 dB(可设计) | 1.15× | 相控阵雷达,兼顾主瓣与旁瓣 |
| 切比雪夫窗 | 等旁瓣(可设计) | 可调 | 需恒定旁瓣电平的场景 |
4.3 阵列增益与分辨率
阵列增益(Array Gain, AG)定义为阵列输出信噪比(SNR)与单阵元输入 SNR 之比。在各阵元噪声独立同分布(空间白噪声)且信号相干的前提下,均匀加权波束形成的阵列增益为 AG = M(即 10·log₁₀M dB)。这意味着 8 元阵列可获得约 9dB 的 SNR 改善,16 元阵列可获得 12dB。然而,当噪声存在空间相关性(如强干扰)时,阵列增益会下降,此时需要采用自适应波束形成在白噪声与干扰之间取得平衡。
分辨率(Resolution)是指阵列区分两个角度相近信号的能力。瑞利准则给出均匀加权 ULA 的最小可分辨角度约为:
Δθ ≈ 0.886·λ / (M·d·cosθ₀)
其中 θ₀ 为波束指向角。在正侧视(θ₀ = 0°)时,Δθ ≈ 0.886λ/L。例如,一个 32 阵元、半波长间距的 ULA,在 3GHz 频率下(λ = 0.1m),孔径 L = 1.55m,其分辨率约为 3.3°。值得注意的是,瑞利准则仅适用于传统波束形成;超分辨算法(如 MUSIC、ESPRIT)可以突破这一限制,但需要较高的 SNR 和精确的阵列校准。
4.4 案例拆解:8 元均匀线阵的波束设计
为将上述概念串联,我们设计一个典型的 8 元 ULA 波束形成方案。假设工作频率 f = 2.4GHz(λ = 0.125m),阵元间距 d = λ/2 = 0.0625m,期望波束指向 θ₀ = 30°。
初始目标与约束:要求第一旁瓣电平低于 -25dB,同时 3dB 主瓣宽度尽可能窄。系统采用 8 通道接收机,ADC 位宽 12bit,基带处理在 FPGA 中完成。
设计过程:
- 相位加权计算:为使波束指向 30°,各阵元所需相位补偿为 φ_m = -2π·m·d·sin(30°)/λ = -π·m/2(m = 0,1,...,7)。即权向量 w_phase = [1, e-jπ/2, e-jπ, e-j3π/2, ..., e-j7π/2]。
- 幅度加权选择:查表可知,汉明窗可提供 -42.8dB 的旁瓣抑制,但主瓣展宽 1.47 倍;泰勒窗(nbar=4)提供 -30dB 旁瓣,主瓣展宽仅 1.15 倍。鉴于旁瓣要求为 -25dB,选择泰勒窗更为合适。
- 综合权向量:w = w_taylor ⊙ w_phase(逐元素相乘),其中 w_taylor 为泰勒窗系数。
- 仿真验证:在 MATLAB 中计算波束图,确认主瓣指向 30°、第一旁瓣 -28.5dB、3dB 主瓣宽度约 14.2°(均匀加权时为 12.6°)。
- FPGA 实现:将权向量量化至 16bit 定点数,存入 ROM。每个采样周期执行 8 次复数乘加运算,流水线时延 3 个时钟周期,满足 100MSPS 采样率要求。
结果与复盘:实测波束图与仿真吻合良好,旁瓣电平 -27.3dB(略高于仿真值,因量化误差所致)。主瓣宽度增加约 1.5°,在可接受范围内。若需进一步压低旁瓣,可改用 40dB 泰勒窗,但主瓣将展宽至 16.5°,需根据探测场景权衡。
4.5 总结与工程要点
本章从波束形成的基本概念出发,建立了从阵列输入到加权输出、再到波束图与性能指标的完整知识框架。核心要点归纳如下:
- 波束形成的本质是空间匹配滤波,权向量 w 的设计决定了阵列的空间响应。
- 波束图是评估阵列性能的“第一眼”工具,主瓣宽度、旁瓣电平、栅瓣位置必须逐一检查。
- 阵列增益与分辨率是相互制约的指标,需结合任务需求(探测距离、角度精度、干扰环境)进行系统级权衡。
- 工程实现中,通道校准、量化误差、阵元互耦是影响理论性能转化为实际性能的三大关键因素。