阻带总体响应约束空域矩阵滤波器设计方法.pdfVIP

水声学 阻带总体响应约束空域矩阵滤波器设计方法 韩东，李启虎 中国科学院声学研究所综合声纳实验室 摘要：针对空域矩阵滤波技术用于实时阵列信号处理必须解决的多阵元、高效率、宽带 等问题，本文提出了双边阻带响应误差整体约束空域矩阵滤波器设计方法。通过建立两 个最优化问题设计滤波器，实现滤波器对左右阻带整体响应的独立约束。并直接给出两 个最优化问题的最优解，本文所提出的滤波器设计方法设计效率高，适用于多阵元情况。 关键词：空域矩阵滤波；阻带响应约束；目标方位估计：匹配场处理； 1引言 空域矩阵滤波技术通过对接收传感器阵列方向向量设置不同的响应，实现对空 间区域入射信号的抑制或无失真通过。由于传感器阵列的方向向量与时域信号傅 立叶变换的展开向量具有类似的形式，所以矩阵滤波器也可被用于做信号滤波， 与传统的FIR滤波器相比，不需要经过一定长度的训练序列就可以直接实现滤波 的收敛效果，尤其在短数据滤波中更具优势【1,2】。现有的矩阵滤波器设计技术主要 包含两大类，一类是与接收阵列数据无关的常规设计方法，另一类是与之相关的 自适应设计方法。 在常规的设计方法中，Maclnnes等提出了通带和阻带响应误差之和最小的最 小二乘矩阵滤波器设计方法，并使用常规方法实现了超出阵元数的目标检测能力 【3】。Vaccaro将滤波器对通带响应误差、阻带响应和过渡带振幅作为目标函数或约 束条件，构造了一系列的最优化问题设计矩阵滤波器，并应用于水声信号匹配场 处理，在有效抑制水面强干扰的情况下，实现了水下弱目标运动轨迹的估计[41。鄢 社锋等人构造了凸规划问题，并将凸规划问题转化为二阶锥规划问题，使用最优 化软件Sedumi实现了平面波目标方位估计及匹配场处理的空域矩阵滤波器设计 【5，6J。Zhu等设计了半无限最优化问题，并将滤波器用于MUSIC波束形成算法，验 证了空域矩阵滤波技术的应用使MUSIC算法的均方根误差接近于克拉美罗界【．71。 韩东等提出了零点约束、通带响应无失真、阻带整体响应约束、通带响应误差约 束、本舰辐射噪声抑制和宽带空域矩阵滤波器等一系列设计方法【89】，将滤波器设 计问题转化为矩阵运算或非线性方程求解问题，并给出了滤波器的闭式最优解， 具有较高的设计效率。自适应空域矩阵滤波方法由Hassanien提出，所构造最优化 问题的目标函数中包含阵列数据协方差矩阵，也即包含输入信号在各频率的能量 信息，利用此信息实现自适应滤波【l叭。 本文提出了双边阻带响应误差整体约束空域矩阵滤波器设计方法，实现了滤波 ．．79．． 器对左右阻带整体响应的独立约束，通过构造两个最优化问题实现。利用Lagrange 乘子理论，推导得出两个最优化问题的最优解，并利用两个最优解之间的关联， 给出最优化问题的最优值。 2双边阻带响应误差整体约束空域矩阵滤波器 考虑由无指向性水听器构成的均匀线列阵，假设阵元数目为Ⅳ，具有相同角频 COo的D个窄带平面波信号源，从D个方向0=b，岛，…，铭]入射到基阵，则接收 阵采样信号可表示为： x(t，COo)=A(e，‰必O，‰)+刀O，04)o) (1) 式中x(t，‰)=[五(f，‰)，…，h 个具有相同角频率‰的信号源，拧(f，‰)=[啊(f，‰)，…，‰t，‰)]1是阵列接收数据 的背景噪声，A(O，‰)=[口(B，‰)，…，口(够，‰)，…，口(铭，‰)]∈口M。是阵列流形矩 是信号源到接收阵的方向向量，这里(．)1表示转置，△代表阵元间隔。 假设通带和左右阻带方向向量构成的阵列流形分别为％∈口M尸、rs。∈口肌墨和 K2∈口M＆，且 咋2 l，pl，．一，，矿…，％j，1≤i≤P 噱1=h，…，％，…，，埚I，1≤，≤S 略2=l叱I，．一，‰，…，kI，1≤k≤是 其中l，，，、％和，砖分别是通带及阻带离散化后的第i、第／和第七个方向向量，P、 S和是分别为相应的离散化方向向量数目。 设