CN112630760A - 多目标条件下的强干扰抑制波束形成器设计方法 - Google Patents

多目标条件下的强干扰抑制波束形成器设计方法 Download PDF

Info

Publication number
CN112630760A
CN112630760A CN202011377290.3A CN202011377290A CN112630760A CN 112630760 A CN112630760 A CN 112630760A CN 202011377290 A CN202011377290 A CN 202011377290A CN 112630760 A CN112630760 A CN 112630760A
Authority
CN
China
Prior art keywords
matrix
filter
formula
array
spatial
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011377290.3A
Other languages
English (en)
Inventor
白云
刘海嫚
杨鑫
陆路
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Haiying Enterprise Group Co Ltd
Original Assignee
Haiying Enterprise Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Haiying Enterprise Group Co Ltd filed Critical Haiying Enterprise Group Co Ltd
Priority to CN202011377290.3A priority Critical patent/CN112630760A/zh
Publication of CN112630760A publication Critical patent/CN112630760A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
    • G01S3/802—Systems for determining direction or deviation from predetermined direction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
    • G01S3/86—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves with means for eliminating undesired waves, e.g. disturbing noises

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明涉及水声设备设计与制造领域技术领域，具体涉及一种多目标条件下的强干扰抑制波束形成器设计方法，包括以下步骤：S1、设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，其中矩阵空域滤波原理；S2、采用稀疏超分辨方位估计算法；S3、基于空域矩阵滤波的稀疏超分辨方位估计：利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计；S4、算法的计算机仿真与实测试验数据验证；设计出具有通带和阻带的空间滤波特性的矩阵滤波器，结合高分辨稀疏方位估计算法，从而在信号处理过程中最大程度的实现对干扰的有效抑制，提高声纳对弱目标的检测能力，提高基阵输出信噪比，从而加强声纳的远程探测能力。

Description

多目标条件下的强干扰抑制波束形成器设计方法
技术领域
本发明涉及水声设备设计与制造领域技术领域，具体涉及一种多目标条件下的强干扰抑制波束形成器设计方法。
背景技术
强干扰环境下的超分辨方位估计是声纳等阵列信号处理领域的一个重要研究课题。例如拖线阵声纳中拖船噪声就是一个强干扰源，它的存在会掩蔽水下目标或者影响水下目标的定位精度。提高强干扰环境条件下的目标方位分辨能力是目前的研究重点。矩阵滤波技术是实现已知方位强干扰抑制的有效手段。矩阵滤波器首先在频率滤波领域被提出和使用，Vaccaro给出了频域矩阵滤波器的设计准则，并指出其设计方法是一种凸优化问题。Macinnes利用伪逆求解矩阵滤波设计问题，得到矩阵滤波设计的最小二乘解，该方法无法给出准确的旁瓣约束和过渡带范围。鄢社峰将矩阵滤波设计的凸优化问题转换为二阶锥优化问题，并利用sedumi软件进行有效求解。
然而在多目标环境下进行的方位估计时，常规的定位方法如常规波束形成方法(CBF)、最小方差无失真波束形成方法(MVDR)虽具有一定的空间分辨能力，但对于空间邻近分布和相关声源的定位问题则束手无策。近些年，稀疏方位估计算法已被广泛使用，如M-FOCUSS、-SVD和SpSF算法等。当方位扫描网格数远大于基阵接收的信号数时，该信号模型具备稀疏性。利用信号模型的稀疏性的先验信息可以提高空间目标的方位估计能力。
重要航道(港口)，复杂作战区海面目标多、干扰强度大、覆盖范围宽等因素严重影响和威胁对目标的远程预警探测能力。
基于此种背景，研究如何在多目标、强干扰的背景下通过高增益干扰抑制波束形成器的设计。
发明内容
针对现有技术的不足，本发明提供了一种多目标条件下的强干扰抑制波束形成器设计方法，本发明利用矩阵滤波技术，设计出具有通带和阻带的空间滤波特性的矩阵滤波器，结合高分辨稀疏方位估计算法，从而在信号处理过程中最大程度的实现对干扰的有效抑制，提高声纳对弱目标的检测能力，提高基阵输出信噪比，从而加强声纳的远程探测能力。对现役岸基声呐改进升级或新型声呐型号研制具有十分重要的现实意义。
本发明通过以下技术方案予以实现：
一种多目标条件下的强干扰抑制波束形成器设计方法，包括以下步骤：
S1、设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，其中矩阵空域滤波原理，
应当满足如下条件：
对于一个由M个阵元组成的已知任意几何形状基阵，假设有D(D<M)个远场窄带平面波信号源从D个方向ΘD＝[θ1,…θd…θD]入射到该基阵，这M个阵元接收到的窄带阵列信号模型可表示为
x(t)＝A(ΘD)s(t)+n(t) (1)；
式(1)中A(ΘD)＝[a(θ1),a(θ2),…,a(θD)]是M×D维阵列流形矩阵，s(t)＝[s1(t),s2(t),…,sD(t)]T是D×1维信号源向量，n(t)是M×1维噪声向量，
对于式(1)设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，即该滤波器能够抑制不感兴趣的方位扇面信号，而让感兴趣的方位扇面信号无失真通过，该矩阵空域滤波器应满足
Figure BDA0002808494540000021
式(2)中ΘP表示通带扇面，ΘS表示阻带扇面，
由式(1)所示的窄带阵列信号处理模型，矩阵滤波器的滤波输出可以表示为
y(t)＝GHx(t)＝GHA(θ)s(t)+GHn(t),t＝1,…,N (3)；
S2、采用稀疏超分辨方位估计算法：
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M，且阵元数M大于信源数N，即Ω＞＞M＞N，
根据MVDR功率谱估计算法可知
Figure BDA0002808494540000022
Figure BDA0002808494540000023
的计算，可以采用迭代方式实现参数估计，在此应用近似最小方差算法完成迭代计算，
近似最小方差算法就是一种使参数估计协方差矩阵CP最小，即达到其估计协方差矩阵下界的估计方法，其表达式为
Figure BDA0002808494540000024
代价函数
Figure BDA0002808494540000025
和扩展不变准则(EXIP)的表达式一致，而EXIP可以等效为大样本条件下的近似最大似然估计，
令r＝vec(R)，阵列输出数据协方差矩阵进行向量化运算可以得到
Figure BDA0002808494540000026
式中vec(·)表示矩阵向量化算子，
Figure BDA0002808494540000027
表示Kronecker乘积，上标“*”为共轭运算符，矩阵
Figure BDA0002808494540000028
和向量
Figure BDA0002808494540000029
分别为
Figure BDA0002808494540000031
根据公式(11)，对于第k个参数
Figure BDA0002808494540000032
设
Figure BDA0002808494540000033
则
Figure BDA0002808494540000034
令
Figure BDA0002808494540000035
将
Figure BDA0002808494540000036
代入上式，可以得到参数
Figure BDA0002808494540000037
的估计值为
Figure BDA0002808494540000038
利用矩阵向量化运算
Figure BDA0002808494540000039
和Kronecker乘积
Figure BDA00028084945400000310
的性质，可以化简为
Figure BDA00028084945400000311
式中Tr(·)为矩阵求迹运算，
将其代入公式中，实现
Figure BDA00028084945400000312
的参数估计，即得到“0”型稀疏近似最小方差算法，简称SAMV-0：
Figure BDA00028084945400000313
式中
Figure BDA00028084945400000314
表示第i次迭代估计出的协方差矩阵，其初始值
Figure BDA00028084945400000315
即为采样协方差矩阵
Figure BDA00028084945400000316
迭代运算的初值可由常规波束形成器输出方位谱给定，迭代截止条件为前后两次迭代估计误差小于某个预设门限，SAMV-0算法和迭代自适应算法的表达式是一致的，如果将
Figure BDA00028084945400000317
代入公式(16)中，可以得到“1”型稀疏近似最小方差算法，简称SAMV-1：
Figure BDA0002808494540000041
S3、基于空域矩阵滤波的稀疏超分辨方位估计：
利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计；
S4、算法的计算机仿真与实测试验数据验证。
优选的，所述步骤S1中，
还应当满足如下条件，最小均方准则：
在一般的情况下，通带内阵列流形向量经过矩阵变换后不一定要保持完全不变，即可以将原来的阵列流形转化为一个虚拟的阵列流形，此时有
Figure BDA0002808494540000042
式中
Figure BDA0002808494540000043
是
Figure BDA0002808494540000044
维虚拟的阵列流形向量，或称期望响应向量，G为
Figure BDA0002808494540000045
维矩阵，这里
Figure BDA0002808494540000046
不一定等于M。
式(4)可以写为
Figure BDA0002808494540000047
令
Figure BDA0002808494540000048
求解矩阵G可以描述为最小均方准则设计问题：
Figure BDA0002808494540000049
式中||·||F表示Frobenius范数，
该最小均方问题的解为
Figure BDA00028084945400000410
式中(·)+表示伪逆。
优选的，所述步骤S1中，
还应当满足如下条件，阻带约束通带Minimax准则：
分别将通带和阻带方位进行离散化，令θj∈Θp(j＝1,…,Np)与θi∈Θs(j＝1,…,Ns)分别是通带与阻带扇面离散化的方位点，拟采用阻带约束通带Minimax准则，
阻带约束通带Minimax准则矩阵滤波器优化设计问题可以表述为
Figure BDA00028084945400000411
式中ξ0为指定的阻带扇面内噪声衰减率，上述问题经过转化可以化为二阶锥规划问题求解。
优选的，所述步骤S1中，
还应当满足如下条件，阻带约束通带最小均方准则：
阻带约束通带最小均方准则矩阵滤波器优化设计问题可以表述为
Figure BDA0002808494540000051
上述问题可以转化为二阶锥规划的形式，可以采用二阶锥规划方法求解，
Figure BDA0002808494540000052
优选的，所述步骤S3中，
具体的，还应当满足如下条件，
设在给定频率f0和空间滤波特性的条件下，使用公式(7)所示的空域矩阵滤波器
Figure BDA0002808494540000053
离散化方位角区域形成Ω个方位扫描网格，形成扫描方位集合ΘΩ，并由公式(16)可知，频域滤波输出后的数据协方差矩阵可以表示为
Figure BDA0002808494540000054
将式(18)代入(17)式，得到基于空域矩阵滤波的稀疏超分辨方位谱估计最终结果
Figure BDA0002808494540000055
优选的，所述步骤S4中，
具体的，计算机仿真还应当满足如下条件：
考虑一个36元均匀直线阵，使用公式(7)所示准则设计空域矩阵滤波器，其空域滤波特性为：通带角度扇区为[-15°,15°]，阻带角度扇区为[-90°,-20°]∪[20°,90°]，直线阵阵元间距为2.5m，在半波长间距条件下对应的处理频率f0为300Hz。阻带衰减量为-25dB，空间方位-40°方向存在一个强干扰，其信号功率为15dB，同时位于-10°和10°方向各存在一个目标信号，其中前者为运动声源，在分析时段内从-10°方向运动到20°方向，后者为固定声源，信号的功率为0dB，噪声功率为6dB,系统采样率为2000Hz，总处理时长为5min。
优选的，所述步骤S4中，
具体的，实测试验数据验证采用湖上试验对大坝发电站强干扰的抑制效果，湖试中采用7个阵元的等间隔线列阵，目标方位125°。
本发明的有益效果为：
本发明设计的多目标条件下的干扰抑制波束形成算法，利用矩阵空域滤波技术对阵列数据进行空域滤波，滤波输出仍保留阵元域数据特性，作为空域预滤波器，能起到净化数据，提高信噪比，改善后续处理性能的目的。矩阵空域滤波技术与高分辨稀疏方位估计技术相结合，具有以下优点：
1)强干扰环境下，有效抑制已知方位的强干扰信号，无失真保留感兴趣的目标信号能量；
2)强干扰环境下，提高了空间目标的方位分辨能力；
3)实现了强干扰环境下的弱目标的超分辨方位估计，提高弱目标检测性能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明中矩阵滤波器空域响应幅度曲线图；
图2为本发明中干扰抑制前方位历程图；
图3为本发明中干扰抑制后方位历程图；
图4为本发明中干扰抑制前某时刻处理结果；
图5为本发明中干扰抑制后某时刻处理结果。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一：
本发明具体公开了一种多目标条件下的强干扰抑制波束形成器设计方法的技术方案，包括以下步骤：
S1、设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，其中矩阵空域滤波原理，
应当满足如下条件：
对于一个由M个阵元组成的已知任意几何形状基阵，假设有D(D<M)个远场窄带平面波信号源从D个方向ΘD＝[θ1,…θd…θD]入射到该基阵，这M个阵元接收到的窄带阵列信号模型可表示为
x(t)＝A(ΘD)s(t)+n(t) (1)；
式(1)中A(ΘD)＝[a(θ1),a(θ2),…,a(θD)]是M×D维阵列流形矩阵，s(t)＝[s1(t),s2(t),…,sD(t)]T是D×1维信号源向量，n(t)是M×1维噪声向量，
对于式(1)设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，即该滤波器能够抑制不感兴趣的方位扇面信号，而让感兴趣的方位扇面信号无失真通过，该矩阵空域滤波器应满足
Figure BDA0002808494540000071
式(2)中ΘP表示通带扇面，ΘS表示阻带扇面，
由式(1)所示的窄带阵列信号处理模型，矩阵滤波器的滤波输出可以表示为
y(t)＝GHx(t)＝GHA(θ)s(t)+GHn(t),t＝1,…,N (3)；
S2、采用稀疏超分辨方位估计算法：
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M，且阵元数M大于信源数N，即Ω＞＞M＞N，
根据MVDR功率谱估计算法可知
Figure BDA0002808494540000072
Figure BDA0002808494540000073
的计算，可以采用迭代方式实现参数估计，在此应用近似最小方差算法完成迭代计算，
近似最小方差算法就是一种使参数估计协方差矩阵
Figure BDA0002808494540000074
最小，即达到其估计协方差矩阵下界的估计方法，其表达式为
Figure BDA0002808494540000075
代价函数
Figure BDA0002808494540000076
和扩展不变准则(EXIP)的表达式一致，而EXIP可以等效为大样本条件下的近似最大似然估计，
令r＝vec(R)，阵列输出数据协方差矩阵进行向量化运算可以得到
Figure BDA0002808494540000077
式中vec(·)表示矩阵向量化算子，
Figure BDA0002808494540000078
表示Kronecker乘积，上标“*”为共轭运算符，矩阵
Figure BDA0002808494540000079
和向量
Figure BDA00028084945400000710
分别为
Figure BDA00028084945400000711
根据公式(11)，对于第k个参数
Figure BDA00028084945400000712
设
Figure BDA00028084945400000713
则
Figure BDA0002808494540000081
令
Figure BDA0002808494540000082
将
Figure BDA0002808494540000083
代入上式，可以得到参数
Figure BDA0002808494540000084
的估计值为
Figure BDA0002808494540000085
利用矩阵向量化运算
Figure BDA0002808494540000086
和Kronecker乘积
Figure BDA0002808494540000087
的性质，可以化简为
Figure BDA0002808494540000088
式中Tr(·)为矩阵求迹运算，
将其代入公式中，实现
Figure BDA0002808494540000089
的参数估计，即得到“0”型稀疏近似最小方差算法，简称SAMV-0：
Figure BDA00028084945400000810
式中
Figure BDA00028084945400000811
表示第i次迭代估计出的协方差矩阵，其初始值
Figure BDA00028084945400000812
即为采样协方差矩阵
Figure BDA00028084945400000813
迭代运算的初值可由常规波束形成器输出方位谱给定，迭代截止条件为前后两次迭代估计误差小于某个预设门限，参考相关文献可知，SAMV-0算法和迭代自适应算法(Iterative Adaptive Approach,IAA)的表达式是一致的，如果将
Figure BDA00028084945400000814
代入公式(16)中，可以得到“1”型稀疏近似最小方差算法，简称SAMV-1：
Figure BDA00028084945400000815
S3、基于空域矩阵滤波的稀疏超分辨方位估计：
“1”型稀疏近似最小方差算法可以有效实现超分辨方位估计，当环境中有强干扰源存在时，强干扰会掩蔽目标的存在或者影响目标的定位精度，利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计；
S4、算法的计算机仿真与实测试验数据验证。
进一步的，步骤S1中，
还应当满足如下条件，最小均方准则：
在一般的情况下，通带内阵列流形向量经过矩阵变换后不一定要保持完全不变，即可以将原来的阵列流形转化为一个虚拟的阵列流形，此时有
Figure BDA0002808494540000091
式中
Figure BDA0002808494540000092
是
Figure BDA0002808494540000093
维虚拟的阵列流形向量，或称期望响应向量，G为
Figure BDA0002808494540000094
维矩阵，这里
Figure BDA0002808494540000095
不一定等于M。
式(4)可以写为
Figure BDA0002808494540000096
令A＝[A(Θp),A(Θs)]，
Figure BDA0002808494540000097
求解矩阵G可以描述为最小均方准则设计问题：
Figure BDA0002808494540000098
式中||·||F表示Frobenius范数，
该最小均方问题的解为
Figure BDA0002808494540000099
式中(·)+表示伪逆。
进一步的，步骤S1中，
还应当满足如下条件，阻带约束通带Minimax准则：
分别将通带和阻带方位进行离散化，令θj∈Θp(j＝1,…,Np)与θi∈Θs(j＝1,…,Ns)分别是通带与阻带扇面离散化的方位点，拟采用阻带约束通带Minimax准则，
阻带约束通带Minimax准则矩阵滤波器优化设计问题可以表述为
Figure BDA00028084945400000910
式中ξ0为指定的阻带扇面内噪声衰减率，上述问题经过转化可以化为二阶锥规划问题求解。
进一步的，步骤S1中，
还应当满足如下条件，阻带约束通带最小均方准则：
阻带约束通带最小均方准则矩阵滤波器优化设计问题可以表述为
Figure BDA0002808494540000101
上述问题可以转化为二阶锥规划的形式，可以采用二阶锥规划方法求解，
Figure BDA0002808494540000102
进一步的，步骤S3中，
具体的，还应当满足如下条件，
设在给定频率f0和空间滤波特性的条件下，使用公式(7)所示的空域矩阵滤波器
Figure BDA0002808494540000103
离散化方位角区域形成Ω个方位扫描网格，形成扫描方位集合ΘΩ，并由公式(16)可知，频域滤波输出后的数据协方差矩阵可以表示为
Figure BDA0002808494540000104
将式(18)代入(17)式，得到基于空域矩阵滤波的稀疏超分辨方位谱估计最终结果
Figure BDA0002808494540000105
进一步的，步骤S4中，
具体的，计算机仿真还应当满足如下条件：
考虑一个36元均匀直线阵，使用公式(7)所示准则设计空域矩阵滤波器，其空域滤波特性为：通带角度扇区为[-15°,15°]，阻带角度扇区为[-90°,-20°]∪[20°,90°]，直线阵阵元间距为2.5m，在半波长间距条件下对应的处理频率f0为300Hz。阻带衰减量为-25dB，空间方位-40°方向存在一个强干扰，其信号功率为15dB，同时位于-10°和10°方向各存在一个目标信号，其中前者为运动声源，在分析时段内从-10°方向运动到20°方向，后者为固定声源，信号的功率为0dB，噪声功率为6dB,系统采样率为2000Hz，总处理时长为5min。
将设计得到的矩阵空域滤波器的空域幅度响应绘制于图1中，从图中可以看出，设计得到的矩阵空域滤波器的空域响应满足设计要求，在通带内空域滤波响应为0dB，阻带部分则严格控制在-25dB以下，且各分析频点的空域响应保持一致。图2是CBF算法下干扰抑制前的结果，图中存在三个目标方位，由于在-40°方向存在强干扰，导致位于正横方向附近的两个能量相对较弱的目标信号容易被背景噪声淹没。图3是CBF算法下干扰抑制后的方位历程图，可见位于-40°方向的强干扰被完全抑制。由于矩阵滤波器滤除了能量较强的干扰，因此弱目标信号的相对能量值变高，在BTR图上的运动角度轨迹更加凸显。
进一步的，步骤S4中，
具体的，实测试验数据验证采用湖上试验对大坝发电站强干扰的抑制效果，湖试中采用7个阵元的等间隔线列阵，目标方位125°。
图4为干扰抑制前的某时刻的波束形成结果，图中可以看到位于60°方位存在强干扰目标为大坝发电站，其能量高于目标能量。图5为应用干扰抑制后的某时刻的波束形成结果，60°方向强干扰被完全抑制，弱目标信号相应增强。同时矩阵空域滤波结合超分辨估计算法，与CBF算法比较，提高了信号处理增益和方位分辨力。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，包括以下步骤：
S1、设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，其中矩阵空域滤波原理，
应当满足如下条件：
对于一个由M个阵元组成的已知任意几何形状基阵，假设有D(D<M)个远场窄带平面波信号源从D个方向ΘD＝[θ1,…θd…θD]入射到该基阵，这M个阵元接收到的窄带阵列信号模型可表示为
x(t)＝A(ΘD)s(t)+n(t) (1)；
式(1)中A(ΘD)＝[a(θ1),a(θ2),…,a(θD)]是M×D维阵列流形矩阵，s(t)＝[s1(t),s2(t),…,sD(t)]T是D×1维信号源向量，n(t)是M×1维噪声向量，
对于式(1)设计一个滤波矩阵G∈CM×M，使得滤波器输出y(t)＝Gx(t),t＝1,…,N具有空域滤波特性，即该滤波器能够抑制不感兴趣的方位扇面信号，而让感兴趣的方位扇面信号无失真通过，该矩阵空域滤波器应满足
Figure FDA0002808494530000011
式(2)中ΘP表示通带扇面，ΘS表示阻带扇面，
由式(1)所示的窄带阵列信号处理模型，矩阵滤波器的滤波输出可以表示为
y(t)＝GHx(t)＝GHA(θ)s(t)+GHn(t),t＝1,…,N (3)；
S2、采用稀疏超分辨方位估计算法：
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M，且阵元数M大于信源数N，即Ω＞＞M＞N，
根据MVDR功率谱估计算法可知
Figure FDA0002808494530000012
Figure FDA0002808494530000013
的计算，可以采用迭代方式实现参数估计，在此应用近似最小方差算法完成迭代计算，
近似最小方差算法就是一种使参数估计协方差矩阵
Figure FDA0002808494530000014
最小，即达到其估计协方差矩阵下界的估计方法，其表达式为
Figure FDA0002808494530000015
代价函数
Figure FDA0002808494530000016
和扩展不变准则(EXIP)的表达式一致，而EXIP可以等效为大样本条件下的近似最大似然估计，
令r＝vec(R)，阵列输出数据协方差矩阵进行向量化运算可以得到
Figure FDA0002808494530000017
Figure FDA0002808494530000021
式中vec(·)表示矩阵向量化算子，
Figure FDA0002808494530000022
表示Kronecker乘积，上标“*”为共轭运算符，矩阵
Figure FDA0002808494530000023
和向量
Figure FDA0002808494530000024
分别为
Figure FDA0002808494530000025
根据公式(11)，对于第k个参数
Figure FDA0002808494530000026
设
Figure FDA0002808494530000027
则
Figure FDA0002808494530000028
令
Figure FDA0002808494530000029
将
Figure FDA00028084945300000210
代入上式，可以得到参数
Figure FDA00028084945300000211
的估计值为
Figure FDA00028084945300000212
利用矩阵向量化运算
Figure FDA00028084945300000213
和Kronecker乘积
Figure FDA00028084945300000214
的性质，可以化简为
Figure FDA00028084945300000215
式中Tr(·)为矩阵求迹运算，
将其代入公式中，实现
Figure FDA00028084945300000216
的参数估计，即得到“0”型稀疏近似最小方差算法，简称SAMV-0：
Figure FDA00028084945300000217
式中
Figure FDA00028084945300000218
表示第i次迭代估计出的协方差矩阵，其初始值
Figure FDA00028084945300000219
即为采样协方差矩阵
Figure FDA0002808494530000031
迭代运算的初值可由常规波束形成器输出方位谱给定，迭代截止条件为前后两次迭代估计误差小于某个预设门限，SAMV-0算法和迭代自适应算法的表达式是一致的，如果将
Figure FDA0002808494530000032
代入公式(16)中，可以得到“1”型稀疏近似最小方差算法，简称SAMV-1：
Figure FDA0002808494530000033
S3、基于空域矩阵滤波的稀疏超分辨方位估计：
利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计；
S4、算法的计算机仿真与实测试验数据验证。
2.根据权利要求1所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S1中，
还应当满足如下条件，最小均方准则：
在一般的情况下，通带内阵列流形向量经过矩阵变换后不一定要保持完全不变，即可以将原来的阵列流形转化为一个虚拟的阵列流形，此时有
Figure FDA0002808494530000034
式中
Figure FDA0002808494530000035
是
Figure FDA0002808494530000036
维虚拟的阵列流形向量，或称期望响应向量，G为
Figure FDA0002808494530000037
维矩阵，这里
Figure FDA0002808494530000038
不一定等于M，式(4)可以写为
Figure FDA0002808494530000039
令A＝[A(Θp),A(Θs)]，
Figure FDA00028084945300000310
求解矩阵G可以描述为最小均方准则设计问题：
Figure FDA00028084945300000311
式中||·||F表示Frobenius范数，
该最小均方问题的解为
Figure FDA00028084945300000312
式中(·)+表示伪逆。
3.根据权利要求1所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S1中，
还应当满足如下条件，阻带约束通带Minimax准则：
分别将通带和阻带方位进行离散化，令θj∈Θp(j＝1,…,Np)与θi∈Θs(j＝1,…,Ns)分别是通带与阻带扇面离散化的方位点，拟采用阻带约束通带Minimax准则，
阻带约束通带Minimax准则矩阵滤波器优化设计问题可以表述为
Figure FDA0002808494530000041
s.t.||GHa(θi)||≤ξ0,θi∈Θs,i＝1,…,Ns (7)；
式中ξ0为指定的阻带扇面内噪声衰减率，上述问题经过转化可以化为二阶锥规划问题求解。
4.根据权利要求1所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S1中，
还应当满足如下条件，阻带约束通带最小均方准则：
阻带约束通带最小均方准则矩阵滤波器优化设计问题可以表述为
Figure FDA0002808494530000042
上述问题可以转化为二阶锥规划的形式，可以采用二阶锥规划方法求解，
Figure FDA0002808494530000043
5.根据权利要求2、3或4任一项所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S3中，
具体的，还应当满足如下条件，
设在给定频率f0和空间滤波特性的条件下，使用公式(7)所示的空域矩阵滤波器
Figure FDA0002808494530000044
离散化方位角区域形成Ω个方位扫描网格，形成扫描方位集合ΘΩ，并由公式(16)可知，频域滤波输出后的数据协方差矩阵可以表示为
Figure FDA0002808494530000045
将式(18)代入(17)式，得到基于空域矩阵滤波的稀疏超分辨方位谱估计最终结果
Figure FDA0002808494530000046
6.根据权利要求5所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S4中，
具体的，计算机仿真还应当满足如下条件：
考虑一个36元均匀直线阵，使用公式(7)所示准则设计空域矩阵滤波器，其空域滤波特性为：通带角度扇区为[-15°,15°]，阻带角度扇区为[-90°,-20°]∪[20°,90°]，直线阵阵元间距为2.5m，在半波长间距条件下对应的处理频率f0为300Hz，阻带衰减量为-25dB，空间方位-40°方向存在一个强干扰，其信号功率为15dB，同时位于-10°和10°方向各存在一个目标信号，其中前者为运动声源，在分析时段内从-10°方向运动到20°方向，后者为固定声源，信号的功率为0dB，噪声功率为6dB,系统采样率为2000Hz，总处理时长为5min。
7.根据权利要求6所述的一种多目标条件下的强干扰抑制波束形成器设计方法，其特征在于，所述步骤S4中，
具体的，实测试验数据验证采用湖上试验对大坝发电站强干扰的抑制效果，湖试中采用7个阵元的等间隔线列阵，目标方位125°。
CN202011377290.3A 2020-11-30 2020-11-30 多目标条件下的强干扰抑制波束形成器设计方法 Pending CN112630760A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011377290.3A CN112630760A (zh) 2020-11-30 2020-11-30 多目标条件下的强干扰抑制波束形成器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011377290.3A CN112630760A (zh) 2020-11-30 2020-11-30 多目标条件下的强干扰抑制波束形成器设计方法

Publications (1)

Publication Number Publication Date
CN112630760A true CN112630760A (zh) 2021-04-09

Family

ID=75306991

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011377290.3A Pending CN112630760A (zh) 2020-11-30 2020-11-30 多目标条件下的强干扰抑制波束形成器设计方法

Country Status (1)

Country Link
CN (1) CN112630760A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113176558A (zh) * 2021-04-20 2021-07-27 哈尔滨工程大学 一种矢量舷侧阵稳健波束形成方法
CN113673419A (zh) * 2021-08-19 2021-11-19 西北工业大学 适用于强干扰环境的波束域快速稀疏贝叶斯方位估计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109116337A (zh) * 2018-07-30 2019-01-01 西北工业大学 一种基于矩阵滤波的稀疏近似最小方差方位估计方法
CN110673119A (zh) * 2019-10-30 2020-01-10 海鹰企业集团有限责任公司 基于压缩感知的非正则化方位估计方法及系统
CN110687528A (zh) * 2019-10-30 2020-01-14 海鹰企业集团有限责任公司 自适应波束形成器生成方法及系统
CN110824428A (zh) * 2019-11-06 2020-02-21 哈尔滨工程大学 一种垂直矢量阵水下声线匹配被动定位方法
CN111273237A (zh) * 2019-05-29 2020-06-12 哈尔滨工程大学 基于空域矩阵滤波和干扰对消的强干扰抑制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109116337A (zh) * 2018-07-30 2019-01-01 西北工业大学 一种基于矩阵滤波的稀疏近似最小方差方位估计方法
CN111273237A (zh) * 2019-05-29 2020-06-12 哈尔滨工程大学 基于空域矩阵滤波和干扰对消的强干扰抑制方法
CN110673119A (zh) * 2019-10-30 2020-01-10 海鹰企业集团有限责任公司 基于压缩感知的非正则化方位估计方法及系统
CN110687528A (zh) * 2019-10-30 2020-01-14 海鹰企业集团有限责任公司 自适应波束形成器生成方法及系统
CN110824428A (zh) * 2019-11-06 2020-02-21 哈尔滨工程大学 一种垂直矢量阵水下声线匹配被动定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨龙: "水下目标方位估计与稳健自适应波束形成方法研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113176558A (zh) * 2021-04-20 2021-07-27 哈尔滨工程大学 一种矢量舷侧阵稳健波束形成方法
CN113176558B (zh) * 2021-04-20 2023-09-29 哈尔滨工程大学 一种矢量舷侧阵稳健波束形成方法
CN113673419A (zh) * 2021-08-19 2021-11-19 西北工业大学 适用于强干扰环境的波束域快速稀疏贝叶斯方位估计方法
CN113673419B (zh) * 2021-08-19 2024-05-28 西北工业大学 适用于强干扰环境的波束域快速稀疏贝叶斯方位估计方法

Similar Documents

Publication Publication Date Title
CN109116337B (zh) 一种基于矩阵滤波的稀疏近似最小方差方位估计方法
CN109143190B (zh) 一种零陷展宽的宽带稳健自适应波束形成方法
CN107153178B (zh) 外辐射源雷达参考信号含有多径干扰时的目标检测方法
CN111273237B (zh) 基于空域矩阵滤波和干扰对消的强干扰抑制方法
CN109212500A (zh) 一种基于稀疏重构的ka-stap杂噪协方差矩阵高精度估计方法
CN103984676A (zh) 一种基于协方差矩阵重构的正交投影自适应波束形成方法
CN109459744B (zh) 一种实现多干扰抑制的稳健自适应波束形成方法
CN112612006B (zh) 基于深度学习的机载雷达非均匀杂波抑制方法
CN109669182B (zh) 无源双基地sar动/静目标联合稀疏成像方法
Hossain et al. Efficient robust broadband antenna array processor in the presence of look direction errors
CN112630760A (zh) 多目标条件下的强干扰抑制波束形成器设计方法
CN112612005B (zh) 基于深度学习的雷达抗主瓣干扰方法
CN113835068B (zh) 一种基于独立成分分析的盲源分离实时抗主瓣干扰方法
CN111693971B (zh) 一种用于弱目标检测的宽波束干扰抑制方法
CN110531311A (zh) 一种基于矩阵重组的lte外辐射源雷达doa估计方法
CN115372925A (zh) 一种基于深度学习的阵列稳健自适应波束形成方法
CN112051552A (zh) 一种基于多站的主瓣抗干扰方法及装置
CN108152796B (zh) 一种基于灰色卡尔曼滤波的主瓣移动干扰消除方法
CN113238194B (zh) 基于分数域-频域处理的宽带相控阵雷达抗诱骗干扰方法
CN105334502A (zh) 基于fir滤波器的雷达空时二维自适应处理方法
Van Lancker Acoustic goniometry: a spatio-temporal approach
CN112649787B (zh) 一种基于低频圆环阵的目标方位估计方法
CN114578311B (zh) 天波超视距雷达特征域抗杂波和干扰方法与装置
CN113094634B (zh) 加权空域处理拖线阵声纳平台辐射噪声限定响应抑制方法
CN116520249A (zh) 基于空域矩阵预滤波的反卷积聚焦定位方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210409