CN109507643B - 旁瓣相消中展宽零陷的方法 - Google Patents
旁瓣相消中展宽零陷的方法 Download PDFInfo
- Publication number
- CN109507643B CN109507643B CN201811293682.4A CN201811293682A CN109507643B CN 109507643 B CN109507643 B CN 109507643B CN 201811293682 A CN201811293682 A CN 201811293682A CN 109507643 B CN109507643 B CN 109507643B
- Authority
- CN
- China
- Prior art keywords
- main
- array
- aux
- auxiliary
- training sample
- 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.)
- Active
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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
-
- 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2925—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by using shape of radiation pattern
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明揭示了一种旁瓣相消中展宽零陷的方法,包括如下步骤:S1、分别获得雷达接收的主阵列训练样本和辅助阵列训练样本;S2、确定零陷需要展宽的宽度;S3、计算主阵列的扰动向量和辅助阵列的扰动向量:S4、计算获得扰动后的主通道训练样本和辅助通道训练样本;S5、对扰动后的主通道训练样本进行波束形成,得到主通道期望输出;S6、计算协方差矩阵和互相关向量;S7、计算展宽后的旁瓣相消权值,并进行旁瓣相消。本发明的方法在旁瓣相消自适应权值训练阶段对训练快拍进行扰动,使得旁瓣相消的合成方向图在干扰方向形成宽零陷,从而提高了雷达在非平稳干扰场景下的抗干扰能力。
Description
技术领域
本发明涉及一种旁瓣相消中展宽零陷的方法,属于雷达技术领域
背景技术
旁瓣相消技术是目前雷达系统中所采用的主要抗干扰手段之一,其需要在主阵列的周围增加若干辅助天线,组成辅助阵列,主阵列做普通波束形成、辅助阵列以接收数据为驱动,通过自适应算法动态调节各辅助天线的复权值。主辅阵列的输出相减,最终合成的方向图既能够保证期望方向的来波信号无失真通过,又能自动在干扰方向形成零陷,抑制旁瓣方向的干扰。
在常规的旁瓣相消过程中,自适应权值是通过从回波中选择一批训练数据、再使用自适应算法学习而得到的。随后将自适应权值用于辅助阵列中接收数据的滤波,完成相消操作,则能够最终实现干扰抑制。如果训练数据和应用数据中干扰的角度一致,则学习得到的自适应权值就能有效地实现应用数据中的干扰抑制。但是在实际的应用过程中,由于天线平台的转动或抖动、干扰源与雷达平台之间的相对运动或干扰传播路径的非平稳性等多方面的原因,干扰的角度在接收数据中常常会呈现出非平稳性。此时,训练数据和应用权值的数据会出现失配,应用数据中的干扰将从旁瓣相消合成方向图中的尖锐零陷中移出,进而导致干扰信号不能得到有效抑制。
针对上述问题,目前业界的主要解决方法是进行零陷展宽,通过对训练样本或协方差矩阵进行处理,可以增加合成方向图中零陷的宽度,从而抑制非平稳干扰。但现有技术中的零陷展宽方法大多不能用于旁瓣相消技术。为此,也有部分研究者提出将基于加窗的零陷展宽方法应用于旁瓣相消技术中。尽管这样的想法能够达到一定的效果,但这样一来,每一个自由度均需要进行一次加窗,无疑大大增加了系统整体的复杂度。
综上所述,如何在现有技术的基础上提出一种新的旁瓣相消中展宽零陷的方法,既能够保证干扰信号得到有效抑制,又不至于导致系统复杂度的过度增加,也就成为了业内技术人员新的研究方向。
发明内容
鉴于现有技术存在上述缺陷,本发明提出了一种旁瓣相消中展宽零陷的方法,包括如下步骤:
S1、在旁瓣相消阶段中分别获得雷达接收的N个主阵列训练样本xmain(n)和辅助阵列训练样本xaux(n);
S2、根据雷达工作参数及实际探测情况,确定零陷需要展宽的宽度Δ;
S3、根据零陷需要展宽的宽度Δ,确定主通道训练样本xmain(n)需要的N个L×1维扰动向量Emain(n)和辅助通道训练样本xaux(n)需要的N个M×1维扰动向量Eaux(n);
S7、利用协方差矩阵R和互相关向量P计算展宽后的旁瓣相消权值w,并进行旁瓣相消。
优选地,S1中所述主阵列训练样本xmain(n)的表达式为
xmain(n)=[xmain,1(n),xmain,2(n),…,xmain,L(n)]T,
所述辅助阵列训练样本xaux(n)的表达式为
xaux(n)=[xaux,1(n),xaux,2(n),…,xaux,M(n)]T,
其中,n表示训练样本编号,n=1,2,…,N,L表示主阵列阵元数,l=1,2,…,L,M表示辅助阵列阵元数,m=1,2,…,M,[·]T表示转置。
优选地,S2中所述零陷需要展宽的宽度Δ的表达式为
其中,d为主阵列阵元间距,λ为雷达工作波长,θ1与θ2分别为所需宽零陷的起始角度与终止角度。
优选地,S3具体包括如下步骤:
S31、根据零陷需要展宽的宽度Δ,按照均匀分布产生主阵列随机变量φmain(n)和辅助阵列随机变量φaux(n),
S32、利用主阵列随机变量φmain(n)和辅助阵列随机变量φaux(n),产生主阵列的扰动向量Emain(n)和辅助阵列的扰动向量Eaux(n),
其中,d1,d2,…,dM为辅助阵列的阵元相对于主阵列相位中心的距离。
优选地,S5具体包括如下步骤:
S51、根据系统预设的主波束指向,得到主阵列的波束形成权值WCBF,
其中,θ0表示主波束指向。
其中,[·]H表示共轭转置。
优选地,S6具体包括如下步骤:
其中,[·]*表示共轭。
优选地,S7中所述展宽后的旁瓣相消权值w的表达式为
w=R-1P,
其中,[·]-1表示矩阵求逆。
本发明的优点主要体现在以下几个方面:
本发明的方法在旁瓣相消自适应权值训练阶段对训练快拍进行扰动,使得旁瓣相消的合成方向图在干扰方向形成宽零陷,从而提高了雷达在非平稳干扰场景下的抗干扰能力。
同时,由于本发明通过对训练样本扰动的形式展宽零陷,与现有技术相比,只需在扰动后进行一次主阵列合成,方法的实现方式更为简单、所需要的运算资源更少,更适合工程实践使用。
此外,本发明也为同领域内的其他相关问题提供了参考,可以以此为依据进行拓展延伸,运用于同领域内其他旁瓣相消技术的设计方案中,具有十分广阔的应用前景。
以下便结合实施例附图,对本发明的具体实施方式作进一步的详述,以使本发明技术方案更易于理解、掌握。
附图说明
图1是本发明的实现流程图;
图2是使用本发明的方法与常规旁瓣相消方法获得的阵列方向图的比较图;
图3是使用本发明的方法与现有方法获得的阵列方向图的比较图。
具体实施方式
如图1所示,本发明揭示了一种旁瓣相消中展宽零陷的方法,包括如下步骤:
S1、在旁瓣相消阶段中分别获得雷达接收的N个主阵列训练样本xmain(n)和辅助阵列训练样本xaux(n)。
所述主阵列训练样本xmain(n)的表达式为
xmain(n)=[xmain,1(n),xmain,2(n),…,xmain,L(n)]T,
所述辅助阵列训练样本xaux(n)的表达式为
xaux(n)=[xaux,1(n),xaux,2(n),…,xaux,M(n)]T,
其中,n表示训练样本编号,n=1,2,…,N,L表示主阵列阵元数,l=1,2,…,L,M表示辅助阵列阵元数,m=1,2,…,M,[·]T表示转置。
S2、根据雷达工作参数及实际探测情况,确定零陷需要展宽的宽度Δ。
所述零陷需要展宽的宽度Δ的表达式为
其中,d为主阵列阵元间距,λ为雷达工作波长,θ1与θ2分别为所需宽零陷的起始角度与终止角度。
S3、根据零陷需要展宽的宽度Δ,确定主通道训练样本xmain(n)需要的N个L×1维扰动向量Emain(n)和辅助通道训练样本xaux(n)需要的N个M×1维扰动向量Eaux(n)。
S3具体包括如下步骤:
S31、根据零陷需要展宽的宽度Δ,按照均匀分布产生主阵列随机变量φmain(n)和辅助阵列随机变量φaux(n),
S32、利用主阵列随机变量φmain(n)和辅助阵列随机变量φaux(n),产生主阵列的扰动向量Emain(n)和辅助阵列的扰动向量Eaux(n),
其中,d1,d2,…,dM为辅助阵列的阵元相对于主阵列相位中心的距离。
S5具体包括如下步骤:
S51、根据系统预设的主波束指向,得到主阵列的波束形成权值WCBF,
其中,θ0表示主波束指向。
其中,[·]H表示共轭转置。
S6具体包括如下步骤:
其中,[·]*表示共轭。
S7、利用协方差矩阵R和互相关向量P计算展宽后的旁瓣相消权值w,并进行旁瓣相消。
所述展宽后的旁瓣相消权值w的表达式为
w=R-1P,
其中,[·]-1表示矩阵求逆。
以下结合仿真实验,详细介绍本发明的使用效果。
实验场景如下:
以一维均匀线阵为例,主阵列阵元数L=20,阵元间距为半波长,辅助阵列阵元数M=8,以主阵列相位中心为参考,主阵列阵元间距为单位,辅助阵列的位置为[-11,-12,-14,-15,11,13,14,17],干扰入射方向为-29°和37°,训练样本数N=300,展宽宽度Δ=0.1。
实验内容与结果如下:
实验1、为证明本发明方法的有效性,利用本发明方法对仿真数据计算合成方向图,并与常规旁瓣相消技术的合成方向图进行比较,比较结果如如图2所示,其中实线表示本发明方法得到的合成方向图,虚线表示常规旁瓣相消方法得到的合成方向图。
从图2中可以看出利用本发明方法可以得到设计的宽零陷,能够有效对抗非平稳干扰。
实验2、为比较本发明方法与现有旁瓣相消下零陷展宽方法,分别用本发明方法和现有方法处理仿真数据,画出合成方向图进行比较,如图3所示,其中,实线表示本发明方法得到的合成方向图,虚线表示现有方法得到的合成方向图。
从图3中可以看出利用本发明方法得到的合成方向图和现有方法得到的合成方向图在各方面性能相同,但本发明方法需要的运算资源更少,更适合工程实践使用。
本发明的方法在旁瓣相消自适应权值训练阶段对训练快拍进行扰动,使得旁瓣相消的合成方向图在干扰方向形成宽零陷,从而提高了雷达在非平稳干扰场景下的抗干扰能力。
同时,由于本发明通过对训练样本扰动的形式展宽零陷,与现有技术相比,只需在扰动后进行一次主阵列合成,方法的实现方式更为简单、所需要的运算资源更少,更适合工程实践使用。
此外,本发明也为同领域内的其他相关问题提供了参考,可以以此为依据进行拓展延伸,运用于同领域内其他旁瓣相消技术的设计方案中,具有十分广阔的应用前景。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神和基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内,不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (7)
1.一种旁瓣相消中展宽零陷的方法,其特征在于,包括如下步骤:
S1、在旁瓣相消阶段中分别获得雷达接收的N个主阵列训练样本xmain(n)和辅助阵列训练样本xaux(n);
S2、根据雷达工作参数及实际探测情况,确定零陷需要展宽的宽度Δ;
S3、根据零陷需要展宽的宽度Δ,确定主通道训练样本xmain(n)需要的N个L×1维扰动向量Emain(n)和辅助通道训练样本xaux(n)需要的N个M×1维扰动向量Eaux(n);
S7、利用协方差矩阵R和互相关向量P计算展宽后的旁瓣相消权值w,并进行旁瓣相消。
2.根据权利要求1所述的旁瓣相消中展宽零陷的方法,其特征在于:S1中所述主阵列训练样本xmain(n)的表达式为
xmain(n)=[xmain,1(n),xmain,2(n),…,xmain,L(n)]T,
所述辅助阵列训练样本xaux(n)的表达式为
xaux(n)=[xaux,1(n),xaux,2(n),…,xaux,M(n)]T,
其中,n表示训练样本编号,n=1,2,…,N,L表示主阵列阵元数,l=1,2,…,L,M表示辅助阵列阵元数,m=1,2,…,M,[·]T表示转置。
7.根据权利要求1所述的旁瓣相消中展宽零陷的方法,其特征在于:S7中所述展宽后的旁瓣相消权值w的表达式为
w=R-1P,
其中,[·]-1表示矩阵求逆。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811293682.4A CN109507643B (zh) | 2018-11-01 | 2018-11-01 | 旁瓣相消中展宽零陷的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811293682.4A CN109507643B (zh) | 2018-11-01 | 2018-11-01 | 旁瓣相消中展宽零陷的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109507643A CN109507643A (zh) | 2019-03-22 |
CN109507643B true CN109507643B (zh) | 2022-10-14 |
Family
ID=65747407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811293682.4A Active CN109507643B (zh) | 2018-11-01 | 2018-11-01 | 旁瓣相消中展宽零陷的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109507643B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110188406B (zh) * | 2019-05-09 | 2022-12-23 | 西安电子科技大学 | 基于旁瓣相消器的自适应零陷展宽方法 |
CN110412552B (zh) * | 2019-07-20 | 2023-02-17 | 中国船舶重工集团公司第七二四研究所 | 一种干扰样本自适应选择的副瓣对消方法 |
CN110868248A (zh) * | 2019-11-21 | 2020-03-06 | 南京邮电大学 | 一种短突发碰撞信号分离方法 |
CN111313948B (zh) * | 2020-02-14 | 2021-05-28 | 北京邮电大学 | 一种信号传输方法、装置及电子设备 |
CN111736119B (zh) * | 2020-06-05 | 2024-01-30 | 西安电子科技大学 | 一种相控阵雷达抗干扰处理的设计方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5134416A (en) * | 1990-01-08 | 1992-07-28 | Cafarelli Nicholas J | Scanning antenna having multipath resistance |
CN101425832B (zh) * | 2008-11-10 | 2012-09-26 | 中兴通讯股份有限公司 | 一种基于零陷展宽的自适应多用户波束成形方法 |
CN101482610B (zh) * | 2008-12-30 | 2012-11-21 | 中国人民解放军空军雷达学院 | 相控阵雷达自适应抑制特殊干扰方法 |
CN102508213B (zh) * | 2011-11-04 | 2013-09-25 | 西安电子科技大学 | 基于零陷加宽的波束域自适应干扰抑制方法 |
CN104345299A (zh) * | 2014-11-03 | 2015-02-11 | 西安电子科技大学 | 基于简化ec的机载mimo雷达空时自适应处理方法 |
US10764102B2 (en) * | 2016-04-25 | 2020-09-01 | Idac Holdings, Inc. | Apparatus and methods for non-systematic complex coded discrete fourier transform spread orthogonal frequency division multiplexing |
CN106842237B (zh) * | 2017-01-18 | 2019-05-07 | 南京理工大学 | 快速任意形状方向图主瓣保形自适应波束形成方法 |
-
2018
- 2018-11-01 CN CN201811293682.4A patent/CN109507643B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109507643A (zh) | 2019-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109507643B (zh) | 旁瓣相消中展宽零陷的方法 | |
Fan et al. | MIMO radar waveform design for quasi-equiripple transmit beampattern synthesis via weighted $ l_p $-minimization | |
CN109143190B (zh) | 一种零陷展宽的宽带稳健自适应波束形成方法 | |
Yang et al. | Robust wideband adaptive beamforming with null broadening and constant beamwidth | |
CN106301498A (zh) | 子带处理方法及频空级联的宽带自适应波束获取方法 | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
KR20160012284A (ko) | 송신신호 파형의 직교성을 이용한 재머 신호 억제 및 원신호 도래각 추정 방법 및 장치 | |
CN109116334A (zh) | 基于超波束加权的声纳波束形成方法及系统 | |
Yang et al. | Robust adaptive beamforming using an iterative FFT algorithm | |
Battista et al. | IRLS based inverse methods tailored to volumetric acoustic source mapping | |
CN111162878B (zh) | 一种基于子阵降维带约束的多域联合抗干扰方法 | |
Godara et al. | Convolution constraints for broadband antenna arrays | |
CN113593596A (zh) | 一种基于子阵划分的鲁棒自适应波束形成定向拾音方法 | |
Al-Sadoon et al. | A more efficient AOA method for 2D and 3D direction estimation with arbitrary antenna array geometry | |
CN113608179B (zh) | 一种雷达信号干扰抑制方法 | |
Pannert | Spatial smoothing for localized correlated sources–Its effect on different localization methods in the nearfield | |
CN113805139A (zh) | 基于聚焦变换的宽带信号稀疏表示波达方向估计方法 | |
Li et al. | A null broadening beamforming method of virtual antenna array | |
CN113391301A (zh) | 基于时间调制线性阵列的超波束形成方法及系统 | |
Wang et al. | Optimal array phase center study for frequency-domain constrained space-time broadband beamforming | |
Lu et al. | Virtual reconstruction-based robust adaptive beamforming for distributed digital subarray antennas | |
Kong et al. | A new reduced-dimension GSC for target tracking and interference suppression | |
Jiang et al. | An efficient ADBF algorithm based on Keystone transform for wideband array system | |
Mengzhong et al. | A method for suppressing grating lobes of wideband reconnaissance DBF | |
CN115453503B (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 | ||
CB02 | Change of applicant information |
Address after: Room 201, building 2, phase II, No.1 Kechuang Road, Yaohua street, Qixia District, Nanjing City, Jiangsu Province Applicant after: NANJING University OF POSTS AND TELECOMMUNICATIONS Address before: 210003 Gulou District, Jiangsu, Nanjing new model road, No. 66 Applicant before: NANJING University OF POSTS AND TELECOMMUNICATIONS |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |