CN112630760A - 多目标条件下的强干扰抑制波束形成器设计方法 - Google Patents
多目标条件下的强干扰抑制波束形成器设计方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000013461 design Methods 0.000 title claims abstract description 28
- 230000001629 suppression Effects 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 108
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 46
- 238000001914 filtration Methods 0.000 claims abstract description 45
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000005094 computer simulation Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000013524 data verification Methods 0.000 claims abstract description 4
- 230000014509 gene expression Effects 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011045 prefiltration Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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具有空域滤波特性,即该滤波器能够抑制不感兴趣的方位扇面信号,而让感兴趣的方位扇面信号无失真通过,该矩阵空域滤波器应满足
式(2)中ΘP表示通带扇面,ΘS表示阻带扇面,
由式(1)所示的窄带阵列信号处理模型,矩阵滤波器的滤波输出可以表示为
y(t)=GHx(t)=GHA(θ)s(t)+GHn(t),t=1,…,N (3);
S2、采用稀疏超分辨方位估计算法:
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M,且阵元数M大于信源数N,即Ω>>M>N,
近似最小方差算法就是一种使参数估计协方差矩阵CP最小,即达到其估计协方差矩阵下界的估计方法,其表达式为
令r=vec(R),阵列输出数据协方差矩阵进行向量化运算可以得到
式中Tr(·)为矩阵求迹运算,
式中表示第i次迭代估计出的协方差矩阵,其初始值即为采样协方差矩阵迭代运算的初值可由常规波束形成器输出方位谱给定,迭代截止条件为前后两次迭代估计误差小于某个预设门限,SAMV-0算法和迭代自适应算法的表达式是一致的,如果将代入公式(16)中,可以得到“1”型稀疏近似最小方差算法,简称SAMV-1:
S3、基于空域矩阵滤波的稀疏超分辨方位估计:
利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计;
S4、算法的计算机仿真与实测试验数据验证。
优选的,所述步骤S1中,
还应当满足如下条件,最小均方准则:
在一般的情况下,通带内阵列流形向量经过矩阵变换后不一定要保持完全不变,即可以将原来的阵列流形转化为一个虚拟的阵列流形,此时有
式(4)可以写为
求解矩阵G可以描述为最小均方准则设计问题:
式中||·||F表示Frobenius范数,
优选的,所述步骤S1中,
还应当满足如下条件,阻带约束通带Minimax准则:
分别将通带和阻带方位进行离散化,令θj∈Θp(j=1,…,Np)与θi∈Θs(j=1,…,Ns)分别是通带与阻带扇面离散化的方位点,拟采用阻带约束通带Minimax准则,
阻带约束通带Minimax准则矩阵滤波器优化设计问题可以表述为
式中ξ0为指定的阻带扇面内噪声衰减率,上述问题经过转化可以化为二阶锥规划问题求解。
优选的,所述步骤S1中,
还应当满足如下条件,阻带约束通带最小均方准则:
阻带约束通带最小均方准则矩阵滤波器优化设计问题可以表述为
上述问题可以转化为二阶锥规划的形式,可以采用二阶锥规划方法求解,
优选的,所述步骤S3中,
具体的,还应当满足如下条件,
将式(18)代入(17)式,得到基于空域矩阵滤波的稀疏超分辨方位谱估计最终结果
优选的,所述步骤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具有空域滤波特性,即该滤波器能够抑制不感兴趣的方位扇面信号,而让感兴趣的方位扇面信号无失真通过,该矩阵空域滤波器应满足
式(2)中ΘP表示通带扇面,ΘS表示阻带扇面,
由式(1)所示的窄带阵列信号处理模型,矩阵滤波器的滤波输出可以表示为
y(t)=GHx(t)=GHA(θ)s(t)+GHn(t),t=1,…,N (3);
S2、采用稀疏超分辨方位估计算法:
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M,且阵元数M大于信源数N,即Ω>>M>N,
令r=vec(R),阵列输出数据协方差矩阵进行向量化运算可以得到
式中Tr(·)为矩阵求迹运算,
式中表示第i次迭代估计出的协方差矩阵,其初始值即为采样协方差矩阵迭代运算的初值可由常规波束形成器输出方位谱给定,迭代截止条件为前后两次迭代估计误差小于某个预设门限,参考相关文献可知,SAMV-0算法和迭代自适应算法(Iterative Adaptive Approach,IAA)的表达式是一致的,如果将代入公式(16)中,可以得到“1”型稀疏近似最小方差算法,简称SAMV-1:
S3、基于空域矩阵滤波的稀疏超分辨方位估计:
“1”型稀疏近似最小方差算法可以有效实现超分辨方位估计,当环境中有强干扰源存在时,强干扰会掩蔽目标的存在或者影响目标的定位精度,利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计;
S4、算法的计算机仿真与实测试验数据验证。
进一步的,步骤S1中,
还应当满足如下条件,最小均方准则:
在一般的情况下,通带内阵列流形向量经过矩阵变换后不一定要保持完全不变,即可以将原来的阵列流形转化为一个虚拟的阵列流形,此时有
式(4)可以写为
求解矩阵G可以描述为最小均方准则设计问题:
式中||·||F表示Frobenius范数,
进一步的,步骤S1中,
还应当满足如下条件,阻带约束通带Minimax准则:
分别将通带和阻带方位进行离散化,令θj∈Θp(j=1,…,Np)与θi∈Θs(j=1,…,Ns)分别是通带与阻带扇面离散化的方位点,拟采用阻带约束通带Minimax准则,
阻带约束通带Minimax准则矩阵滤波器优化设计问题可以表述为
式中ξ0为指定的阻带扇面内噪声衰减率,上述问题经过转化可以化为二阶锥规划问题求解。
进一步的,步骤S1中,
还应当满足如下条件,阻带约束通带最小均方准则:
阻带约束通带最小均方准则矩阵滤波器优化设计问题可以表述为
上述问题可以转化为二阶锥规划的形式,可以采用二阶锥规划方法求解,
进一步的,步骤S3中,
具体的,还应当满足如下条件,
将式(18)代入(17)式,得到基于空域矩阵滤波的稀疏超分辨方位谱估计最终结果
进一步的,步骤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具有空域滤波特性,即该滤波器能够抑制不感兴趣的方位扇面信号,而让感兴趣的方位扇面信号无失真通过,该矩阵空域滤波器应满足
式(2)中ΘP表示通带扇面,ΘS表示阻带扇面,
由式(1)所示的窄带阵列信号处理模型,矩阵滤波器的滤波输出可以表示为
y(t)=GHx(t)=GHA(θ)s(t)+GHn(t),t=1,…,N (3);
S2、采用稀疏超分辨方位估计算法:
稀疏处理的必要条件为扫描网格点数Ω远大于阵元数M,且阵元数M大于信源数N,即Ω>>M>N,
令r=vec(R),阵列输出数据协方差矩阵进行向量化运算可以得到
式中Tr(·)为矩阵求迹运算,
式中表示第i次迭代估计出的协方差矩阵,其初始值即为采样协方差矩阵迭代运算的初值可由常规波束形成器输出方位谱给定,迭代截止条件为前后两次迭代估计误差小于某个预设门限,SAMV-0算法和迭代自适应算法的表达式是一致的,如果将代入公式(16)中,可以得到“1”型稀疏近似最小方差算法,简称SAMV-1:
S3、基于空域矩阵滤波的稀疏超分辨方位估计:
利用空域矩阵滤波和“1”型稀疏近似最小方差算法可以实现强干扰环境下的弱目标的超分辨方位估计;
S4、算法的计算机仿真与实测试验数据验证。
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°。
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)
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)
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 | 哈尔滨工程大学 | 基于空域矩阵滤波和干扰对消的强干扰抑制方法 |
-
2020
- 2020-11-30 CN CN202011377290.3A patent/CN112630760A/zh active Pending
Patent Citations (5)
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)
Title |
---|
杨龙: "水下目标方位估计与稳健自适应波束形成方法研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》 * |
Cited By (4)
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 |
---|---|---|
CN104181552B (zh) | 一种动态gnss接收机的抗干扰正态零陷加宽的方法 | |
CN109116337B (zh) | 一种基于矩阵滤波的稀疏近似最小方差方位估计方法 | |
CN109143190B (zh) | 一种零陷展宽的宽带稳健自适应波束形成方法 | |
CN106569181A (zh) | 基于协方差矩阵重构稳健Capon波束形成的算法 | |
CN107153178B (zh) | 外辐射源雷达参考信号含有多径干扰时的目标检测方法 | |
CN111273237B (zh) | 基于空域矩阵滤波和干扰对消的强干扰抑制方法 | |
CN109212500A (zh) | 一种基于稀疏重构的ka-stap杂噪协方差矩阵高精度估计方法 | |
CN103984676A (zh) | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 | |
CN109459744B (zh) | 一种实现多干扰抑制的稳健自适应波束形成方法 | |
Hossain et al. | Efficient robust broadband antenna array processor in the presence of look direction errors | |
CN112612006B (zh) | 基于深度学习的机载雷达非均匀杂波抑制方法 | |
CN112630760A (zh) | 多目标条件下的强干扰抑制波束形成器设计方法 | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
CN110261826A (zh) | 一种零陷展宽的相干干扰抑制方法 | |
CN115372925A (zh) | 一种基于深度学习的阵列稳健自适应波束形成方法 | |
CN112051552A (zh) | 一种基于多站的主瓣抗干扰方法及装置 | |
CN108152796B (zh) | 一种基于灰色卡尔曼滤波的主瓣移动干扰消除方法 | |
CN107167809B (zh) | 一种基于信号子空间聚焦的宽带阻塞阵波束形成方法 | |
CN113238194B (zh) | 基于分数域-频域处理的宽带相控阵雷达抗诱骗干扰方法 | |
CN107728132B (zh) | 一种改善干扰阻塞算法输出信噪比的方法 | |
Van Lancker | Acoustic goniometry: a spatio-temporal approach | |
CN112649787B (zh) | 一种基于低频圆环阵的目标方位估计方法 | |
CN114578311B (zh) | 天波超视距雷达特征域抗杂波和干扰方法与装置 | |
CN113094627B (zh) | 加权空域处理拖线阵声纳平台辐射噪声零响应抑制方法 | |
CN113094634B (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 |