一种结合白噪声增益约束的精确阵列响应控制方法及装置制造方法及图纸

本发明专利技术提供了一种结合白噪声增益约束的精确阵列响应控制方法及装置，所述方法适用于时域波束形成器，所述方法包括：初始化滤波器系数，设置初始白噪声增益约束条件及旁瓣级约束条件；选定波束图中需要进行控制的位置，在对应频率与角度添加虚拟干扰源，更新滤波器系数，计算波束形成器的白噪声增益，若不满足白噪声增益约束条件，则修改旁瓣级约束条件重新更新滤波器系数，直至计算得到滤波器系数同时满足白噪声增益约束条件和旁瓣级约束条件或达到设定迭代次数后停止迭代，得到最终的波束图。本发明专利技术的优势在于：可用于时域结构的任意阵型传声器阵列；实现了对波束图中指定位置响应的精确控制，具有灵活的波束设计性能。具有灵活的波束设计性能。具有灵活的波束设计性能。

【技术实现步骤摘要】
一种结合白噪声增益约束的精确阵列响应控制方法及装置


[0001]本专利技术属于声学信号处理领域，具体涉及一种结合白噪声增益约束的精确阵列响应控制方法及装置。

技术介绍

[0002]传声器阵列由按一定规律排列的多个传声器构成，广泛应用于声学信号处理领域，如语音交互、电话会议等。传声器阵列能够提取特定方向的期望信号，同时抑制不同方向的干扰和环境噪声，具有优秀的空间滤波性能。由于语音信号频带较宽，在处理时需要采用宽带波束形成技术。在宽带波束形成的时域与频域实现方法中，时域方法采用滤波后求和的结构，往往应用于要求低处理时延的场景。根据滤波器系数是否依赖接受数据，可将波束形成器分为自适应波束形成器和固定波束形成器。典型的自适应波束形成器包括线性约束最小方差波束形成器(Linearly Constrained Minimum Variance,LCMV)和最小均方无失真响应(Minimum Variance Distortionless Response,MVDR)波束形成器等。固定波束形成器包括延迟求和(Delay
‑
and
‑
Sum,DAS)波束形成器和超指向性(Superdirective,SD)波束形成器等。实际应用中，自适应波束形成器计算量较大且当存在误差时可能损伤目标信号，因此在要求低计算量与高稳健性的应用场景中，固定波束形成器往往得到更多的应用。对于固定波束形成器，如何设计滤波器系数以获得理想的波束图是至关重要的。
[0003]阵列波束设计问题也被称为阵列波束图综合问题。在过去的几十年中，大量的研究致力于阵列波束图综合问题。自1990年，许多学者选择引入自适应波束形成方法中的零陷机制进行窄带波束形成器的设计。Olen和Compton提出了凹槽噪声法，然而该方法在选择虚拟干扰功率时采用迭代的方式，设计效率较低。针对窄带波束设计，张学敬等学者提出了最优精确阵列响应控制算法(Optimal Precise Array Response Control,OPARC)，通过迭代更新窄带虚拟干扰与噪声协方差矩阵(Virtual Interference
‑
plus
‑
Noise Covariance Matrix,VINCM)，能够实现对阵列响应的精确控制。但是上述方法仅适用于窄带信号，无法应用于时域结构的宽带波束形成器。对于时域滤波后求和结构的宽带波束形成器，鄢社锋等提出将阵列波束图综合问题转化为二阶锥规划问题(Second Order Cone Programming,SOCP)问题，并采用基于内点法(Interior
‑
Point Method,IPM)的工具包如CVX求解得到最优滤波器系数。然而，基于优化算法的波束设计方法常使用基于内点法的工具包进行求解，计算时间较长且在不合理的约束条件下可能会出现无解的情况。同时，基于优化算法的波束设计方法预先设定期望波束图，对于其旁瓣级(SideLobe Level,SLL)与白噪声增益(White Noise Gain,WNG)需要通过多次预实验进行确定。因此，如何在无需预实验的条件下进行宽带波束设计，是亟待研究的问题。

技术实现思路

[0004]本专利技术的目的在于克服现有技术进行宽带阵列波束设计时需要进行多次预实验的缺陷。
[0005]为了实现上述目的，本专利技术提出了一种结合白噪声增益约束的精确阵列响应控制方法，所述方法适用于时域波束形成器，所述时域波束形成器在每个传声器之后接有一组有限脉冲响应滤波器；
[0006]所述方法包括：
[0007]步骤S1：初始化滤波器系数，设置初始白噪声增益约束条件及旁瓣级约束条件；
[0008]步骤S2：选定波束图中需要进行控制的位置，在对应频率与角度添加虚拟干扰源，更新滤波器系数，计算波束形成器的白噪声增益，若不满足白噪声增益约束条件，则修改旁瓣级约束条件重新更新滤波器系数，直至计算得到滤波器系数同时满足白噪声增益约束条件和旁瓣级约束条件或达到设定迭代次数后停止迭代，得到最终的滤波器系数。
[0009]作为上述方法的一种改进，所述初始化滤波器系数，具体包括：
[0010]构建初始宽带虚拟干扰与噪声协方差矩阵为单位阵；计算初始滤波器系数h0：
[0011][0012]其中，表示目标方向约束矩阵，θ0表示声源方向，表示取实部操作，表示取虚部操作，f
k
，k＝1,
…
,K表示受到控制的第k个频点对应的频率，K表示所控制的频率数；T表示矩阵转置；u(f
K
,θ0)表示对于位于方向θ0、频率为f
K
的远场信号，时域波束形成器的导向向量：
[0013][0014]其中，表示克罗内克积；e(f
k
)表示L维傅里叶变换因子，L表示有限脉冲响应滤波器的阶数；a(f
k
,θ0)表示M维频域导向向量，M表示传声器个数；
[0015]e(f
k
)＝[1,exp(
‑
j2πf
k
/f
s
),...,exp(
‑
j2π(L
‑
1)f
k
/f
s
)]T
[0016]a(f
k
,θ0)＝exp[
‑
j2πf
k
sin(θ0)d/c][0017]其中，M维向量d包含每个传声器与参考位置间距离，c表示声速，f
s
表示采样频率；
[0018]g＝[1,...,1,0,...,0]由K个1和K个0组成；表示ML
×
ML单位矩阵。
[0019]作为上述方法的一种改进，所述选定波束图中需要进行控制的位置，具体为：
[0020][0021]其中，f
j
，θ
j
分别为第j次迭代控制点的频率与角度；B为控制频带；Θ为控制角度范围；P
d
(f,θ)为期望阵列响应；
[0022]归一化后的阵列功率响应P
j
(f,θ)为：
[0023][0024]其中，h
j
表示第j次迭代时的滤波器系数；θ0表示声源方向；f0表示归一化参考频率。
[0025]作为上述方法的一种改进，所述更新滤波器系数，具体为：
[0026]h
j
＝H
j
[(S
j
)4,(S
j
)3,(S
j
)2,S
j
,1]T
[0027]其中，S
j
表示虚拟干扰的功率，由下式计算得到：
[0028][(S
j
)4,(S
j
)3,(S
j
)2,S
j
,1]G
j
[(S
j
...

【技术保护点】

【技术特征摘要】
1.一种结合白噪声增益约束的精确阵列响应控制方法，所述方法适用于时域波束形成器，所述时域波束形成器在每个传声器之后接有一组有限脉冲响应滤波器；所述方法包括：步骤S1：初始化滤波器系数，设置初始白噪声增益约束条件及旁瓣级约束条件；步骤S2：选定波束图中需要进行控制的位置，在对应频率与角度添加虚拟干扰源，更新滤波器系数，计算波束形成器的白噪声增益，若不满足白噪声增益约束条件，则修改旁瓣级约束条件重新更新滤波器系数，直至计算得到滤波器系数同时满足白噪声增益约束条件和旁瓣级约束条件或达到设定迭代次数后停止迭代，得到最终的滤波器系数。2.根据权利要求1所述的结合白噪声增益约束的精确阵列响应控制方法，其特征在于，所述初始化滤波器系数，具体包括：构建初始宽带虚拟干扰与噪声协方差矩阵为单位阵；计算初始滤波器系数h0：其中，表示目标方向约束矩阵，θ0表示声源方向，表示取实部操作，表示取虚部操作，f
k
，k＝1,
…
,K表示受到控制的第k个频点对应的频率，K表示所控制的频率数；T表示矩阵转置；u(f
K
,θ0)表示对于位于方向θ0、频率为f
K
的远场信号，时域波束形成器的导向向量：其中，表示克罗内克积；e(f
k
)表示L维傅里叶变换因子，L表示有限脉冲响应滤波器的阶数；a(f
k
,θ0)表示M维频域导向向量，M表示传声器个数；e(f
k
)＝[1,exp(
‑
j2πf
k
/f
s
),
…
,exp(
‑
j2π(L
‑
1)f
k
/f
s
)]
T
a(f
k
,θ0)＝exp[
‑
j2πf
k
sin(θ0)d/c]其中，M维向量d包含每个传声器与参考位置间距离，c表示声速，f
s
表示采样频率；g＝[1,
…
,1,0,
…
,0]由K个1和K个0组成；表示ML
×
ML单位矩阵。3.根据权利要求2所述的结合白噪声增益约束的精确阵列响应控制方法，其特征在于，所述选定波束图中需要进行控制的位置，具体为：其中，f
j
，θ
j
分别为第j次迭代控制点的频率与角度；B为控制频带；Θ为控制角度范围；P
d
(f,θ)为期望阵列响应；归一化后的阵列功率响应P
j
(f,θ)为：其中，h
j
表示第j次迭代时的滤波器系数；θ0表示声源方向；f0表示归一化参考频率。4.根据权利要求3所述的结合白噪声增益约束的精确阵列响应控制方法，其特征在于，所述更新滤波器系数，具体为：h
j
＝H
j
[(S
j
)4,(S
j
)3,(S
j
)2,S
j
,1]
T
其中，S
j
表示...

【专利技术属性】
技术研发人员：厉剑，徐嘉懿，孟维鑫，郑成诗，李晓东，
申请(专利权)人：中国科学院声学研究所，
类型：发明
国别省市：