3、波束形成核心思想:空域滤波、导向矢量、波束方向图与主瓣、旁瓣、栅瓣的概念
波束形成的本质,是将一组按特定几何位置排列的传感器(通常是均匀线阵、均匀面阵或圆阵)采集到的信号,通过加权求和的方式合成一个标量输出。这个过程的物理含义,并非简单的“把几路信号加起来”,而是利用各阵元接收信号之间的相位差,在空间中构造出一个具有方向选择性的响应函数。换句话说,波束形成是一种空域滤波器——它允许来自期望方向的信号以低失真通过,同时抑制来自其他方向的干扰和噪声。理解这一思想,需要从三个递进的层次展开:首先明确空域滤波与常规时域滤波的对应关系;其次掌握导向矢量这一描述阵列几何与信号传播关系的数学工具;最后通过波束方向图直观地评估滤波性能,并理解主瓣、旁瓣、栅瓣等关键指标对系统设计的影响。
上图展示了波束形成从物理接收到数学表达再到性能评估的完整链路。阵列接收信号后,通过导向矢量描述各阵元间的相位关系,再经加权求和实现空域滤波,最终以波束方向图作为衡量滤波效果的核心工具。下面逐一深入剖析每个环节。
3.1 空域滤波:从时域滤波到空间维度的扩展
核心原理:传统数字信号处理中的FIR/IIR滤波器,是在时间维度上对信号的采样序列进行加权求和,以保留或抑制特定频率分量。波束形成将这一思想推广到空间维度——阵列中每个阵元相当于空间中的一个“采样点”,各阵元接收信号之间的相位差由信号到达方向(DOA)决定。通过对各阵元输出施加复权值(包含幅度和相位),可以在空间上形成“通带”和“阻带”。例如,当期望信号来自30°方向时,权矢量设计为对该方向信号进行同相叠加,使输出最大;而对来自-60°方向的干扰,则通过相位反相使其相互抵消。
数学表达:设M元均匀线阵,阵元间距为d,窄带信号以角度θ入射。第m个阵元相对于参考阵元的相位延迟为:
τ_m(θ) = (m-1)·d·sinθ / c
其中c为光速。对应的导向矢量(Steering Vector)为:
a(θ) = [1, e^{j·2π·d·sinθ/λ}, e^{j·4π·d·sinθ/λ}, ..., e^{j·(M-1)·2π·d·sinθ/λ}]^T
阵列输出为:
y(t) = w^H · x(t) = Σ_{m=1}^{M} w_m^* · x_m(t)
其中w为权矢量,x(t)为各阵元接收信号矢量。当w = a(θ₀)时,输出在θ₀方向形成最大增益,这就是常规波束形成(CBF)。
3.2 导向矢量:阵列响应的“空间指纹”
核心原理:导向矢量是描述阵列对某一方向入射信号的复数响应向量。它包含了阵列几何结构(阵元位置)、信号频率、传播介质等全部信息。对于均匀线阵,导向矢量具有范德蒙德结构,其相位项线性递增,这为快速DOA估计算法(如ESPRIT)提供了数学基础。导向矢量不仅是波束形成权值设计的依据,也是阵列校准、DOA估计、干扰抑制等算法的核心输入。
关键性质:
- 归一化正交性:不同方向导向矢量的内积|a^H(θ₁)·a(θ₂)|随角度差增大而减小,但并非严格正交,这决定了波束方向图的旁瓣结构。
- 频率依赖性:导向矢量与信号频率相关,宽带信号需采用聚焦变换或子带分解处理。
- 阵列流形:所有可能方向的导向矢量集合构成阵列流形,是阵列信号处理的“字典”。
实际应用中的注意事项:
- 阵元位置误差、幅相不一致会导致导向矢量失配,需定期校准。
- 互耦效应会使实际导向矢量偏离理想模型,高频段尤为明显。
- 对于非均匀阵列(如稀疏阵、圆阵),导向矢量不再具有范德蒙德结构,需数值计算。
设载频f = 10GHz(λ = 3cm),阵元间距d = 1.5cm(λ/2),信号入射角θ = 20°。
相邻阵元相位差:Δφ = 2π·d·sinθ/λ = 2π·0.015·sin20°/0.03 = 2π·0.5·0.342 = 1.074 rad ≈ 61.5°
导向矢量:a(20°) = [1, e^{j61.5°}, e^{j123°}, e^{j184.5°}, e^{j246°}, e^{j307.5°}, e^{j369°}, e^{j430.5°}]
若采用矩形窗权值w = a(20°),则输出在20°方向增益为8倍(18dB),而0°方向的响应为:
|a^H(20°)·a(0°)| = |Σ_{m=0}^{7} e^{-jm·61.5°}| = |sin(8·61.5°/2)/sin(61.5°/2)| ≈ 0.85(约-1.4dB)
这说明矩形窗旁瓣较高,实际系统常采用切比雪夫窗或泰勒窗将旁瓣压低至-30dB以下。
3.3 波束方向图:性能评估的核心工具
核心原理:波束方向图定义为阵列输出幅度(或功率)随入射角变化的函数,数学上为B(θ) = |w^H·a(θ)|。它直观地展示了阵列对不同方向信号的增益响应,是衡量波束形成性能最直接的指标。方向图通常以分贝(dB)表示,归一化后主瓣峰值对应0dB。
关键参数与判断方法:
| 参数 | 定义 | 典型值/设计要求 | 对系统的影响 |
|---|---|---|---|
| 主瓣宽度(3dB波束宽度) | 主瓣峰值两侧功率下降3dB对应的角度范围 | θ_3dB ≈ 0.886λ/(M·d·cosθ₀) | 决定角度分辨率,越窄越好 |
| 第一旁瓣电平(SLL) | 最大旁瓣峰值相对主瓣峰值的比值 | 矩形窗:-13dB;切比雪夫窗:-30~-40dB | 决定抗干扰能力,越低越好 |
| 栅瓣位置 | 与主瓣等高的虚假响应方向 | 出现在sinθ_g = sinθ₀ ± kλ/d | 导致方向模糊,必须避免 |
| 零点位置 | 响应为零的方向 | 由权矢量决定,可自适应控制 | 用于抑制特定方向干扰 |
主瓣、旁瓣、栅瓣的物理意义:
- 主瓣:阵列增益最大的方向,是期望信号的通过区域。主瓣越窄,空间分辨率越高,但代价是旁瓣电平升高(不确定性原理)。
- 旁瓣:主瓣以外的响应峰,是干扰和噪声进入阵列的“后门”。强干扰从旁瓣进入时,即使主瓣增益很高,输出信干噪比(SINR)仍可能严重恶化。
- 栅瓣:当阵元间距d > λ/2时,导向矢量在多个方向出现相同的相位分布,导致方向图出现与主瓣等高的虚假峰。栅瓣会造成DOA估计的模糊,是阵列设计中的“红线”问题。
3.4 案例拆解:从方向图设计到系统性能验证
场景描述:某5G基站采用16元均匀线阵,工作频率3.5GHz(λ≈85.7mm),阵元间距d=42.9mm(λ/2)。系统要求:波束指向30°时,第一旁瓣电平≤-25dB,且不得出现栅瓣。同时,需要在-20°方向形成零点以抑制邻区干扰。
设计过程:
- 初始设计:采用矩形窗权值w = a(30°),仿真得到方向图。主瓣宽度约6.3°,但第一旁瓣电平为-13.2dB,不满足要求。
- 加窗优化:改用64dB切比雪夫窗,旁瓣电平降至-25dB,但主瓣展宽至8.1°,角度分辨率略有下降。权衡后接受该折中。
- 零点约束:在切比雪夫窗基础上,采用线性约束最小方差(LCMV)算法,在-20°方向施加零陷。权矢量更新为:
w_LCMV = R^{-1}·C·(C^H·R^{-1}·C)^{-1}·f
其中C = [a(30°), a(-20°)],f = [1, 0]^T,R为接收信号协方差矩阵(此处用仿真数据估计)。
结果对比:
| 方案 | 主瓣宽度 | 第一旁瓣 | -20°响应 | 栅瓣 |
|---|---|---|---|---|
| 矩形窗 | 6.3° | -13.2dB | -18dB | 无 |
| 切比雪夫窗 | 8.1° | -25dB | -22dB | 无 |
| LCMV零陷 | 8.3° | -25dB | -45dB | 无 |
复盘与代价:LCMV方案在-20°方向实现了45dB的零陷,但主瓣宽度略有增加,且对导向矢量误差敏感——当实际来波方向偏离30°达2°时,主瓣增益下降约1.5dB。工程中需结合DOA估计精度,在零陷深度与稳健性之间取得平衡。
3.5 常见误区与总结复盘
误区一:认为阵元越多越好。阵元数量增加确实能收窄主瓣、提高增益,但代价是硬件通道数、功耗和成本线性增长。且当阵元间距固定时,阵元数增加意味着阵列孔径增大,旁瓣数量增多,若不加窗处理,旁瓣电平反而可能恶化。
误区二:忽视阵元间距与栅瓣的关系。部分设计者为节省阵元数量,采用d > λ/2的稀疏布阵,却未意识到栅瓣风险。正确做法是:若必须稀疏布阵,需通过优化算法(如遗传算法)设计非均匀阵元位置,使栅瓣转化为较低的旁瓣。
误区三:将方向图视为静态指标。实际系统中,阵元幅相误差、互耦、安装平台反射都会使方向图畸变。建议在暗室进行方向图实测,并与仿真对比,必要时引入校准算法(如旋转矢量法)修正权值。
① 空域滤波是波束形成的本质,权矢量设计决定了滤波特性;
② 导向矢量是连接阵列几何与信号处理的桥梁,必须精确建模;
③ 方向图的主瓣宽度、旁瓣电平、栅瓣位置是三大核心指标,设计时需统筹权衡;
④ 工程实现中,务必通过实测验证方向图,并考虑误差鲁棒性。