基于空域矩阵滤波和干扰对消的强干扰抑制方法技术

本发明专利技术是基于空域矩阵滤波和干扰对消的强干扰抑制方法。本发明专利技术主要用于抑制观测扇面外和观测扇面内的强干扰。通过子空间矩阵滤波，减小观测扇面外干扰影响，进而获得观测扇面内强干扰方位，并将其应用于阻塞矩阵的设计，构建新的不降低数据维数的阻塞阵。通过阻塞矩阵和空域矩阵滤波器对阵列接收数据进行处理，最后采用MUSIC谱进行方位估计。本发明专利技术在抑制观测扇面内强干扰的同时，保留了邻近方位的弱目标信息，实现强干扰条件下的弱目标方位估计。本发明专利技术属于一种水声阵列信号处理方法，可应用于阵列信号处理、弱目标方位探测等领域。

A strong interference suppression method based on spatial matrix filtering and interference cancellation

【技术实现步骤摘要】
基于空域矩阵滤波和干扰对消的强干扰抑制方法
本专利技术涉及水声工程
，是一种基于空域矩阵滤波和干扰对消的强干扰抑制方法。
技术介绍
基于阵列信号处理的DOA估计方法广泛应用于声纳、雷达等领域，其中，强干扰背景下的目标方位估计首先要对干扰进行抑制。空域矩阵滤波器应用于DOA估计，可以抑制阻带区域的干扰和噪声，提高了低信噪比下的方位估计精度。但是，当有一个强干扰位于观测扇面内距离目标方位比较近时，矩阵空域滤波器失效。传统的阻塞矩阵，也可以用于去除强干扰影响，但是此技术需要预知强干扰的方位，并且会降低阻塞矩阵维数。
技术实现思路
本专利技术为解决上述现有技术存在的问题，提供了一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，本专利技术提供了以下技术方案：一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，所述方法包括如下步骤：步骤一：建立阵元信号接收模型，设定矩阵滤波器通带，确定通带内强干扰方位；步骤二：根据估计的强干扰方位，构造阻塞矩阵，对阻塞矩阵接收数据进行处理，去除通带内的强干扰；步骤三：去除阻塞矩阵对协方差矩阵的影响；步骤四：对协方差矩阵进行空域滤波处理，获得MUSIC空间谱。优选地，所述步骤一具体为：第一步：建立阵元信号接收模型，所述模型由N个阵元组成的均匀线列阵，有P个远场窄带平面波信号从不同方向入射到所述模型，所述阵元间距为信号频率的半波长，计算得到阵元接收信号，所述阵元接收信号通过下式表示：n(t)＝[n1(t),n2(t),...,nN(t)]T(2)其中，x(t)为阵元接收信号，si(t)为源信号，P为远场窄带平面波信号数量，a(θi)为第i个源信号的阵列流形向量，n(t)为噪声向量，[·]T为转置，N为阵元数；第二步：设定各源信号之间以及源信号与噪声之间互不相干，各阵元的背景噪声为相互独立的零均值高斯白噪声，建立阵列接收数据的协方差矩阵，所述协方差矩阵通过下式表示：其中，Rx为阵列接收数据的协方差矩阵，K表示快拍总数，k表示快拍数，pi表示第i个源信号的能量，I为单位矩阵，[·]H表示共轭转置,表示噪声能量,λi为特征值，ηi表示第i个特征值并且满足η1＞η2＞···＞ηN，ui是ηi对应的特征向量；第三步：设定矩阵滤波器的通带为Θ，阻带为对通带Θ外的强干扰采用空域矩阵滤波器对接收的阵列信号数据进行滤波，使得通带Θ内的信号无失真通过，所述空域矩阵滤波器预期响应通过下式表示：其中，T为空域矩阵滤波器，a(θ)为指定角度域的导向矢量，为预期响应；第四步：设定空域矩阵滤波器，对通带响应进行约束优化，使得阻带响应最大误差最小化，根据所述约束优化得到下式：minr其中，||·||2表示向量的2-范数。为通带响应误差约束，r为阻带约束；第五步：估计通道Θ内的强干扰方位，所述通带Θ内的强干扰方位估计通过下式表示：UI＝[u1,u2,...,uP-1](7)其中，为通带Θ内的子空间方位谱，a(θm)为通带内方位的导向矢量，θm∈Θ，为通带内强干扰方位的导向矢量的估计值，为通带内强干扰的估计方位，UI为Rx中的强干扰子空间。优选地，所述步骤二具体为：第一步：根据估计出的构建阻塞矩阵，通过下式表示阻塞矩阵：其中，B为N×N的阻塞矩阵，d为阵元间距，λ为源信号的波长；第二步：源信号经过阻塞矩阵处理，处理后的源信号通过下式表示：y(t)＝Bx(t)(9)其中，y(t)为经过阻塞矩阵处理后的阵列数据。优选地，所述步骤三具体为：第一步：采用阻塞矩阵处理协方差矩阵，处理后的协方差矩阵通过下式表示：其中，为采用阻塞矩阵处理的协方差矩阵；第二步：减小阻塞矩阵对协方差矩阵中高斯白噪声的影响，通过下式表示减小阻塞矩阵影响的协方差矩阵：其中，为减小阻塞矩阵影响的协方差矩阵，为的估计值。优选地，所述通过下式进行估计：优选地，所述步骤四具体为：第一步：采用矩阵滤波器处理协方差矩阵，通过下式表示处理后的协方差矩阵：其中，为矩阵滤波器处理后的协方差矩阵，γi表示特征值，ei是与γi对应的特征矢量(i＝1,...,N)；第二步：根据下式求解MUSIC空间谱：其中，a(θk)为搜索导向矢量，Pmusic为MUSIC空间谱，En＝[eP,eP+1,...,eN]为噪声空间。优选地，特征值γi满足γ1≥γ2≥···≥γP≥···γN，且Es＝[e1,e2,...,eP-1]，Es为协方差矩阵中去除通带内强干扰后，包含目标方位的信号空间。本专利技术具有以下有益效果：(1)在多个强干扰下，本专利技术方法不仅可以抑制观测扇面外的强干扰，还可以抑制观测扇面内目标方位较近的强干扰。(2)与传统的阻塞矩阵相比，本专利技术构建了不降低维数的新阻塞矩阵，更完整地保留了目标信息。(3)与传统DOA估计方法相比，本专利技术在抑制了强干扰之后，提高了弱目标DOA估计能力。附图说明图1是基于空域矩阵滤波和阻塞矩阵处理的强干扰抑制算法流程图。图2是强干扰下MUSIC谱估计结果。图3是本专利技术在强干扰下目标方位估计结果。具体实施方式以下结合具体实施例，对本专利技术进行了详细说明。具体实施例一：根据图1所示，本专利技术提供一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，所述方法包括如下步骤：步骤一：建立阵元信号接收模型，设定矩阵滤波器通带，确定通带内强干扰方位；步骤二：根据估计的强干扰方位，构造阻塞矩阵，对阻塞矩阵接收数据进行处理，去除通带内的强干扰；步骤三：去除阻塞矩阵对协方差矩阵的影响；步骤四：对协方差矩阵进行空域滤波处理，获得MUSIC空间谱。步骤一：通带内强干扰方位的确定。首先建立阵列信号接收模型，对于一个由N阵元组成的均匀线列阵，假设有P个远场窄带平面波信号从不同方向入射到该基阵，阵元间距为信号频率的半波长，阵列接收信号可表示为其中，si(t)为源信号，包括P-1个强干扰信号和1个目标信号。本专利技术假设：一个干扰位于目标观测扇面内，其他干扰位于观测扇面外，干扰强度高于目标。第i个源信号的阵列流形向量为a(θi)，噪声向量为n(t)＝[n1(t),n2(t),...,nN(t)]T，表示各阵元接收的背景复噪声，[·]T表示转置。这里假设各源信号之间以及源信号与噪声之间互不相关，且各阵元的背景噪声为相互独立的零均值高斯白噪声，方差为阵列接收数据的协方差矩阵表示为其中pi表示第i个源信号的能量，K表示快拍总数，k表示快拍数，表示噪声能量，I为单位矩阵。ηi为特征值，ui是ηi对应的特征向量，特征值满足η1≥η2≥···≥ηP≥···ηN，此处假设信源数本文档来自技高网...

【技术保护点】
1.一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，其特征是：所述方法包括如下步骤：/n步骤一：建立阵元信号接收模型，设定矩阵滤波器通带，确定通带内强干扰方位；/n步骤二：根据估计的强干扰方位，构造阻塞矩阵，对阻塞矩阵接收数据进行处理，去除通带内的强干扰；/n步骤三：去除阻塞矩阵对协方差矩阵的影响；/n步骤四：对协方差矩阵进行空域滤波处理，获得MUSIC空间谱。/n

【技术特征摘要】
1.一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，其特征是：所述方法包括如下步骤：
步骤一：建立阵元信号接收模型，设定矩阵滤波器通带，确定通带内强干扰方位；
步骤二：根据估计的强干扰方位，构造阻塞矩阵，对阻塞矩阵接收数据进行处理，去除通带内的强干扰；
步骤三：去除阻塞矩阵对协方差矩阵的影响；
步骤四：对协方差矩阵进行空域滤波处理，获得MUSIC空间谱。


2.根据权利要求1所述的一种基于空域矩阵滤波和干扰对消的强干扰抑制方法，其特征是：所述步骤一具体为：
第一步：建立阵元信号接收模型，所述模型由N个阵元组成的均匀线列阵，有P个远场窄带平面波信号从不同方向入射到所述模型，所述阵元间距为信号频率的半波长，计算得到阵元接收信号，所述阵元接收信号通过下式表示：



n(t)＝[n1(t),n2(t),...,nN(t)]T(2)
其中，x(t)为阵元接收信号，si(t)为源信号，P为远场窄带平面波信号数量，a(θi)为第i个源信号的阵列流形向量，n(t)为噪声向量，[·]T为转置，N为阵元数；
第二步：设定各源信号之间以及源信号与噪声之间互不相干，各阵元的背景噪声为相互独立的零均值高斯白噪声，建立阵列接收数据的协方差矩阵，所述协方差矩阵通过下式表示：



其中，Rx为阵列接收数据的协方差矩阵，K表示快拍总数，k表示快拍数，pi表示第i个源信号的能量，I为单位矩阵，[·]H表示共轭转置,表示噪声能量,λi为特征值，ηi表示第i个特征值并且满足η1＞η2＞···＞ηN，ui是ηi对应的特征向量；
第三步：设定矩阵滤波器的通带为Θ，阻带为对通带Θ外的强干扰采用空域矩阵滤波器对接收的阵列信号数据进行滤波，使得通带Θ内的信号无失真通过，所述空域矩阵滤波器预期响应通过下式表示：



其中，T为空域矩阵滤波器，a(θ)为指定角度域的导向矢量，为预期响应；
第四步：设定空域矩阵滤波器，对通带响应进行约束优化，使得阻带响应最大误差最小化，根据所述约束优化得到下式：



其中，||·||2表示向量的2-范数。为通带响应误差约束，r为阻带约束；
第五步：估计通道Θ内的强干扰方位，所述通带Θ内的强干扰方位估计通过下式表示：



UI＝...

【专利技术属性】
技术研发人员：邹男，柳国龙，付进，梁国龙，齐滨，邱龙皓，王逸林，孙思博，李晨牧，李雄辉，张文琪，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙;23