RESEARCH / PRACTICE / NOTES

常规波束形成:从时延补偿到方位谱

从均匀线阵的传播时延和导向向量推导延迟求和波束形成,并用双目标仿真说明孔径、角度间隔和噪声对方位分辨的影响。

#波束形成#阵列信号处理
关键词DOA估计
连续狭缝的影像在一点汇聚
本文目录 · 3 个章节
  1. 从传播时延到加权向量
  2. 从输出功率到方位谱
  3. 双目标仿真

波束形成是阵列信号处理中最基础的一类方法,直到现在仍是水声与语音阵列系统的组成部分。它一方面把观察方向上的信号相干叠加,以改善输出信噪比;另一方面,通过扫描不同方向并比较输出功率,得到方位谱,用于估计信号的到达方向。两种用途使用的是同一套空间加权思想,但评价标准并不完全相同。

从传播时延到加权向量

对远场窄带平面波,阵列上各阵元接收的是同一源信号的不同延迟版本。设 NN 个特性相同的阵元组成均匀线阵,相邻间距为 dd,入射角 θ\theta 从阵列法线(正横方向)起算,波速为 cc,则相邻阵元的几何时延幅度为

τ(θ)=dsin⁡θc.\tau(\theta)=\frac{d\sin\theta}{c}.

若频率为 ff、角频率为 ω=2πf\omega=2\pi f、波长为 λ=c/f\lambda=c/f,相邻阵元的相位差可写成 ωτ=2π(d/λ)sin⁡θ\omega\tau=2\pi(d/\lambda)\sin\theta。τ\tau 是时间,不能把波长和时间直接相乘。时延的正负号取决于阵元编号与来波方向约定;下文固定一种约定,把导向向量写成

a(θ)=[1e−jωτ(θ)⋯e−jω(N−1)τ(θ)] ⁣T.\mathbf a(\theta)= \begin{bmatrix} 1 & e^{-j\omega\tau(\theta)} & \cdots & e^{-j\omega(N-1)\tau(\theta)} \end{bmatrix}^{\!\mathrm T}.

以 θ0\theta_0 为观察方向,均匀延迟求和权向量取 w(θ0)=a(θ0)/N\mathbf w(\theta_0)=\mathbf a(\theta_0)/N。若接收向量为 x\mathbf x,输出是

y=wH(θ0)x=∑m=1Nwm∗(θ0)xm.y=\mathbf w^{\mathrm H}(\theta_0)\mathbf x =\sum_{m=1}^{N}w_m^*(\theta_0)x_m.

此时 wH(θ0)a(θ0)=1\mathbf w^{\mathrm H}(\theta_0)\mathbf a(\theta_0)=1:观察方向上的信号得到相位补偿后相干相加。求和式里的权值应取共轭,因为前面采用的是 wHx\mathbf w^{\mathrm H}\mathbf x 记法。若改用 ∑mwmxm\sum_m w_mx_m,就必须同时修改权向量定义,不能把两种约定混在一起。

从宽带信号看,“补偿时延”可以是时域分数延时,也可以在每个频点施加相位权值。固定频率的导向向量只描述窄带情形;把它直接用于整个语音带宽,会使不同频率上的指向性不同。

从输出功率到方位谱

设阵列协方差矩阵为 R=E[xxH]\mathbf R=\mathbb E[\mathbf x\mathbf x^{\mathrm H}],则扫描方向 θ\theta 的常规波束形成功率为

PCBF(θ)=wH(θ)Rw(θ).P_{\mathrm{CBF}}(\theta)= \mathbf w^{\mathrm H}(\theta)\mathbf R\mathbf w(\theta).

实际只有有限个快拍 x1,…,xL\mathbf x_1,\ldots,\mathbf x_L,通常用样本协方差代替理论期望:

R^=1L∑ℓ=1LxℓxℓH=1LXXH,X=[x1,…,xL].\widehat{\mathbf R}= \frac{1}{L}\sum_{\ell=1}^{L}\mathbf x_\ell\mathbf x_\ell^{\mathrm H} =\frac{1}{L}\mathbf X\mathbf X^{\mathrm H}, \qquad \mathbf X=[\mathbf x_1,\ldots,\mathbf x_L].

在独立同分布的零均值复高斯样本等特定模型下,它也与协方差的极大似然估计一致;一般情况下称为样本协方差更准确。用它代入 PCBFP_{\mathrm{CBF}},便得到可以在离散角度网格上计算的方位谱。谱峰说明该扫描方向具有较高输出功率,但不保证每个真实声源都对应一个可分开的峰。

常规波束形成的分辨力主要受阵列孔径、信号频率、角度位置和噪声条件限制。阵元数不变时,增大孔径通常能缩窄主瓣,但如果相邻间距超过半波长,视场内可能出现空间混叠。快拍数增加可使协方差估计更稳定,却不能把由孔径造成的主瓣宽度无限缩小。

双目标仿真

考虑 16 个各向同性阵元组成的均匀线阵,阵元间距为 λ/2\lambda/2,两个等强度窄带信号频率为 2 kHz,背景为高斯白噪声,信号与噪声互不相关,快拍数为 1024。以阵列法线为 0∘0^\circ,改变两个目标的角度间隔和输入信噪比,可得到下列方位谱。

两个目标位于负 5 度与正 5 度,15 dB 信噪比时的方位谱

两目标分别在 −5∘-5^\circ 与 5∘5^\circ,间隔 10∘10^\circ,谱中可以看到明显分开的两个峰。

两个目标位于负 5 度与正 2 度,15 dB 信噪比时的方位谱

把间隔收窄到 7∘7^\circ,两个峰开始靠近,仍能在这组条件下区分。

两个目标位于负 5 度与正 1 度,15 dB 信噪比时的方位谱

间隔进一步收窄至 6∘6^\circ,两个目标落进一片合并的主瓣里,常规波束形成无法从这张谱上稳定读出两个独立方位。

两个目标位于负 5 度与正 2 度,10 dB 信噪比时的方位谱

保持 7∘7^\circ 间隔而降低信噪比,谱峰相对噪声和旁瓣的区分度下降。这里的 6∘6^\circ、7∘7^\circ 只是在上述阵列和仿真条件下的观察,不是常规波束形成的一条固定分辨率界限;改动孔径、频率、目标方位、信号相关性或快拍数,结论都会变化。

从这些图也能看出,扫描方位谱与提高输出信噪比虽然共用同一个波束形成器,却要关注不同问题。用于拾音时,要看目标方向的增益、通带宽度和干扰抑制;用于 DOA 时,要看相邻目标能否形成可区分的谱峰。只追求主瓣变窄而忽略旁瓣、扫描范围与麦克风失配,往往不能得到稳定的方位估计。

交流与讨论

评论加载中…

随便读一篇 →

站内搜索

输入关键词开始搜索。

搜索内容不会上传 Esc 关闭

订阅明月入怀

将下面的网址添加到你的 RSS 阅读器。