CN113406620B - 一种阵列分解的分布式阵列测角方法 - Google Patents

一种阵列分解的分布式阵列测角方法 Download PDF

Info

Publication number
CN113406620B
CN113406620B CN202110467816.5A CN202110467816A CN113406620B CN 113406620 B CN113406620 B CN 113406620B CN 202110467816 A CN202110467816 A CN 202110467816A CN 113406620 B CN113406620 B CN 113406620B
Authority
CN
China
Prior art keywords
array
interference
distributed
weight vector
angle measurement
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
Application number
CN202110467816.5A
Other languages
English (en)
Other versions
CN113406620A (zh
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202110467816.5A priority Critical patent/CN113406620B/zh
Publication of CN113406620A publication Critical patent/CN113406620A/zh
Application granted granted Critical
Publication of CN113406620B publication Critical patent/CN113406620B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/87Combinations of radar systems, e.g. primary radar and secondary radar
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66Radar-tracking systems; Analogous systems
    • G01S13/68Radar-tracking systems; Analogous systems for angle tracking only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/023Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/42Diversity systems specially adapted for radar

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

本发明提供了一种阵列分解的分布式阵列测角方法,基于子阵划分的思想,将阵列分解为NB=K+M个子阵,每个子阵包含NA个阵元,则阵元级子阵列记作阵列A,子阵级子阵列记作阵列B;构成大口径分布式阵列,有效抑制干扰;在抑制干扰的同时,精确测量目标角度;本发明的方法适用于分布式抗干扰阵列,是一种在有效抑制干扰的前提下精确测量目标角度的方法,能够有效解决传统方法在干扰抑制后,阵列和差方向图畸变,无法测角的问题。

Description

一种阵列分解的分布式阵列测角方法
技术领域
本发明涉及角度测量领域,具体涉及一种阵列分解的分布式阵列测角方法。
背景技术
目标角度估计对于目标定位、跟踪有着非常重要的意义,它是制约雷达探测性能的重要因素。单脉冲测角方法利用同时产生的“和”、“差”波束来测量波束宽度范围内的目标角度,原理简单,易于工程实现,已经广泛应用于雷达***中。
随着电子干扰技术的不断发展,雷达目标探测环境变得越来越复杂,自适应单脉冲测角方法有效结合了自适应波束形成技术和单脉冲测角方法,当目标被干扰信号掩盖时,首先采用自适应波束形成技术对干扰信号进行有效抑制,然后使用单脉冲测角方法进行目标角度测量。当干扰信号位于雷达副瓣位置时,该方法具有良好的性能,然而,当干扰信号位于雷达主瓣位置时,自适应波束形成在对消干扰的同时,会造成目标回波能量严重损失,无法完成目标探测。
因此,针对主瓣干扰的场景,可以在主雷达周围增加辅助接收单元形成分布式阵列,通过对多部小孔径单元雷达进行信号级的融合处理,使其性能等效于大孔径雷达。借助大孔径阵列的空间高分辨率优势,在保护目标回波能量的同时对消主瓣干扰。然而,大孔径阵列存在栅瓣,存在测角模糊不利于测角。因此需要对分布式阵列对消主瓣干扰下的目标角度估计算法进行深入研究。
发明内容
有鉴于此,针对上述已有技术的不足,本发明提供了一种阵列分解的分布式阵列测角方法,能够在有效抑制主瓣干扰的同时,精确测量目标角度。
一种阵列分解的分布式阵列测角方法,包括如下步骤:
步骤一、分布式雷达阵列分解:
设分布式雷达包含一个主阵和M个辅助阵,且均为相控阵体制,其中,设主阵包含NM个阵元,阵元间距为λ/2;辅助阵包含NA个阵元,阵元间距为λ/2,且主阵的阵元数与辅助阵的阵元数满足整数倍关系,即NM=K·NA;λ表示信号波长;
则分布式雷达共有(K+M)NA个阵元,将其分解为NB=K+M个子阵,每个子阵包含NA个阵元,将阵元级子阵列记作阵列A,子阵级子阵列记作阵列B;
步骤二、阵列A单脉冲测角:
将阵列A关于天线中心对称放置,天线孔径应该划分为两个象限,每个象限有NA/2个阵元,用该两个象限中阵元完成目标角度测量,由此确定阵列A的权矢量,设为WA(θ);
步骤三、采用阵列B抑制主瓣干扰,由此获得阵列B的权矢量,设为WB(θ);
步骤四、根据如下公式计算分布式雷达权矢量:
Figure BDA0003044796400000021
利用权矢量WC(θ)实现目标角度的测量。
较佳的,阵列B是由NB个阵元构成的线性阵列,该NB个阵元非均匀布阵构成大口径分布式阵列。
较佳的,所述步骤三中,采用自适应波束形成或者主分量求逆方法抑制主瓣干扰。
较佳的,当采用主分量求逆方法抑制主瓣干扰时,具体过程为:
远场处有一个期望信号和一个干扰以平面波入射,到达角度分别为θt和θj,将阵列B接收的快拍数据表示为:
X=a(θt)st+a(θj)sj+n (2)
式中X为NB×1阵列数据向量,n为NB×1维阵列噪声向量;st为目标信源的复包络,sj为干扰信源的复包络,a(θt),a(θj)分别为目标信源和干扰信源的导向矢量;
权矢量为把期望信号导向矢量向与干扰子空间正交的空间投影得到,投影矩阵为:
Figure BDA0003044796400000022
P为一个NB×NB的系数矩阵,其对角线即为抗主瓣干扰权矢量WB(θ),干扰对消后的回波写作:
Figure BDA0003044796400000023
较佳的,所述步骤二中,选定一组加权矢量来构造和、差波束,通过比较完成目标角度测量。
本发明具有如下有益效果:
本发明提供了一种阵列分解的分布式阵列测角方法,基于子阵划分的思想,将阵列分解为NB=K+M个子阵,每个子阵包含NA个阵元,则阵元级子阵列记作阵列A,子阵级子阵列记作阵列B;构成大口径分布式阵列,有效抑制干扰;在抑制干扰的同时,精确测量目标角度;本发明的方法适用于分布式抗干扰阵列,是一种在有效抑制干扰的前提下精确测量目标角度的方法,能够有效解决传统方法在干扰抑制后,阵列和差方向图畸变,无法测角的问题。
附图说明
图1为分布式雷达示意图;
图2为分布式雷达阵列分解示意图;
图3为第一级子阵的和差方向图;
图4为第二级子阵自适应干扰抑制方向图;
图5为分布式阵列雷达全阵和差方向图;
图6为分布式阵列雷达全阵鉴角曲线。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
一种阵列分解的分布式阵列测角方法,包括如下步骤:
步骤一、分布式雷达阵列分解
分布式雷达包含一个主阵和M个辅助阵,均为相控阵体制,其中,主阵包含NM个阵元,阵元间距为λ/2;辅助阵包含NA个阵元,阵元间距同样为λ/2,需要注意的是,主阵阵元数与辅助阵阵元数满足整数倍关系,即NM=K·NA。分布式雷达如图1所示。λ表示信号波长;
因此,分布式雷达共有(K+M)NA个阵元,基于子阵划分的思想,可以将阵列分解为NB=K+M个子阵,每个子阵包含NA个阵元,则阵元级子阵列记作阵列A,子阵级子阵列记作阵列B,阵列的分解如图2所示。
阵列A的相对位置矢量为dA=[dA(1),dA(2),...,dA(NA)]T,阵列B的相对位置矢量为dB=[dB(1),dB(2),...,dB(NB)]T,则合成阵列的相对位置矢量是两个对应子阵列的相对位置矢量的Kronecker和,写作:
Figure BDA0003044796400000041
阵列A的相对导向矢量为aA(θ),阵列B的相对导向矢量为aB(θ),则合成阵列的相对导向矢量是两个对应子阵列的相对导向矢量的Kronecker积,写作:
Figure BDA0003044796400000042
阵列A的权矢量为WA(θ),阵列B的权矢量为WB(θ),则合成阵列的权矢量是两个对应子阵列的权矢量的Kronecker积,写作
Figure BDA0003044796400000043
阵列A的方向图为FA(θ;WA(θ)),阵列B的方向图为FB(θ;WB(θ)),由Kronecker积的混合乘积性质,原分布式雷达的方向图是两个对应子阵列的方向图的乘积,写作
Figure BDA0003044796400000044
因此,经过上述步骤实现了对分布式雷达的阵列分解。
步骤二、阵列A单脉冲测角。
假设阵列A是由NA个阵元构成的线性阵列,为了实现单脉冲角度估计,阵列雷达的天线单元关于天线中心对称放置。NA元线性阵列的天线孔径应该划分为两个象限,每个象限有NA/2个阵元,选定一组加权矢量w,wΔ来构造和、差波束,通过比较完成目标角度测量。
步骤三、阵列B干扰抑制
阵列B是由NB个阵元构成的线性阵列,为了有效抑制阵列干扰,NB个阵元非均匀布阵构成大口径分布式阵列。远场处有一个期望信号和一个干扰以平面波入射(波长为λ),到达角度分别为θt和θj,阵列接收的快拍数据可表示为
X=a(θt)st+a(θj)sj+n (5)
式中X为NB×1阵列数据向量,X=[x1,x2,…,xNB]T。[]T表示矩阵转置,n为NB×1阵列噪声向量,
Figure BDA0003044796400000051
st为目标信源的复包络,sj为干扰信源的复包络,a(θt),a(θj)分别为目标信源和干扰信源的导向矢量。
接下来可用自适应波束形成、主分量求逆等方法抑制主瓣干扰。本文采用主分量求逆法进行分析。主分量求逆法的基本思想是由阵列接收数据估计出强干扰分量,然后从原始数据中减去这些分量(从而抑制干扰),再对剩余数据向量进行空域匹配滤波(以期望信号导向矢量a(θ0)为滤波权矢量)得到阵列输出。换言之,权矢量是把期望信号导向矢量向与干扰子空间正交的空间(噪声子空间)投影得到的。写作
Xc=(I-a(θj)aHj))X (6)
=PX
其中,P=(I-a(θj)aHj))是一个NB×NB的系数矩阵。求得的输出信号Xc已不包含干扰信号。
步骤四、完成干扰抑制下的目标角度测量
阵列A获得单脉冲测角和、差方向图,阵列B获得干扰位置零陷方向图,经过上述处理,分布式雷达能够获得在干扰位置生成零陷的和、差方向图,目标回波能量没有损失,且和差方向图不会畸变。因而,可以在有效抑制干扰的前提下,对目标角度进行测量。
下面给出一个应用本发明的仿真实例,并对角度测量过程进行分析与说明。
分布式阵列雷达由一部主雷达和若干部辅助雷达组成,主雷达和辅助雷达都是相控阵体制,主辅雷达满足两个条件:
主阵和辅阵采用相同的阵元间隔;
主阵尺寸与辅助阵尺寸成整数倍关系;
因此,根据上文分析,该分布式阵列可以分为两级。
仿真参数如下表所示。
表1仿真参数表
Figure BDA0003044796400000061
步骤1:第一级子阵采用相位和差单脉冲测角法,根据和、差波束权矢量,计算和、差天线方向图,如图3所示。
步骤2:第二级子阵采用主分量相消法获得抗干扰权矢量,从而获得在干扰位置产生零陷的自适应干扰抑制方向图,如图4所示。由于阵列稀疏布置,天线方向图有些杂乱,但从图中明显看出,在干扰位置形成了很深的零陷。
步骤3:和、差波束权矢量分别和抗干扰权矢量相乘,获得分布式阵列雷达每个阵元抗干扰下的和、差波束权矢量,从而计算出在干扰位置产生零陷的和、差方向图,如图5所示。
步骤4:构造和差波束,完成分布式阵列雷达存在干扰下的目标角度测量,鉴角曲线如图6所示,从图中可以看出,鉴角曲线没有因为干扰抑制而产生突变,因此,能够精确测量目标角度。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种阵列分解的分布式阵列测角方法,其特征在于,包括如下步骤:
步骤一、分布式雷达阵列分解:
设分布式雷达包含一个主阵和M个辅助阵,且均为相控阵体制,其中,设主阵包含NM个阵元,阵元间距为λ/2;辅助阵包含NA个阵元,阵元间距为λ/2,且主阵的阵元数与辅助阵的阵元数满足整数倍关系,即NM=K·NA;λ表示信号波长;
则分布式雷达共有(K+M)NA个阵元,将其分解为NB=K+M个子阵,每个子阵包含NA个阵元,将阵元级子阵列记作阵列A,子阵级子阵列记作阵列B;
步骤二、阵列A单脉冲测角:
将阵列A关于天线中心对称放置,天线孔径应该划分为两个象限,每个象限有NA/2个阵元,用该两个象限中阵元完成目标角度测量,由此确定阵列A的权矢量,设为WA(θ);
步骤三、采用阵列B抑制主瓣干扰,由此获得阵列B的权矢量,设为WB(θ);
步骤四、根据如下公式计算分布式雷达权矢量:
Figure FDA0003044796390000011
利用权矢量WC(θ)实现目标角度的测量。
2.如权利要求1所述的一种阵列分解的分布式阵列测角方法,其特征在于,阵列B是由NB个阵元构成的线性阵列,该NB个阵元非均匀布阵构成大口径分布式阵列。
3.如权利要求1所述的一种阵列分解的分布式阵列测角方法,其特征在于,所述步骤三中,采用自适应波束形成或者主分量求逆方法抑制主瓣干扰。
4.如权利要求3所述的一种阵列分解的分布式阵列测角方法,其特征在于,当采用主分量求逆方法抑制主瓣干扰时,具体过程为:
远场处有一个期望信号和一个干扰以平面波入射,到达角度分别为θt和θj,将阵列B接收的快拍数据表示为:
X=a(θt)st+a(θj)sj+n (2)
式中X为NB×1阵列数据向量,n为NB×1维阵列噪声向量;st为目标信源的复包络,sj为干扰信源的复包络,a(θt),a(θj)分别为目标信源和干扰信源的导向矢量;
权矢量为把期望信号导向矢量向与干扰子空间正交的空间投影得到,投影矩阵为:
Figure FDA0003044796390000021
P为一个NB×NB的系数矩阵,其对角线即为抗主瓣干扰权矢量WB(θ),干扰对消后的回波写作:
Figure FDA0003044796390000022
5.如权利要求1所述的一种阵列分解的分布式阵列测角方法,其特征在于,所述步骤二中,选定一组加权矢量来构造和、差波束,通过比较完成目标角度测量。
CN202110467816.5A 2021-04-28 2021-04-28 一种阵列分解的分布式阵列测角方法 Active CN113406620B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110467816.5A CN113406620B (zh) 2021-04-28 2021-04-28 一种阵列分解的分布式阵列测角方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110467816.5A CN113406620B (zh) 2021-04-28 2021-04-28 一种阵列分解的分布式阵列测角方法

Publications (2)

Publication Number Publication Date
CN113406620A CN113406620A (zh) 2021-09-17
CN113406620B true CN113406620B (zh) 2022-08-19

Family

ID=77678060

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110467816.5A Active CN113406620B (zh) 2021-04-28 2021-04-28 一种阵列分解的分布式阵列测角方法

Country Status (1)

Country Link
CN (1) CN113406620B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102064892A (zh) * 2010-12-28 2011-05-18 哈尔滨工业大学 子阵级自适应单脉冲的两级干扰抑制方法
CN103018730B (zh) * 2012-11-27 2015-04-08 西安电子科技大学 分布式子阵波达方向估计方法
JP6584714B1 (ja) * 2018-08-28 2019-10-02 三菱電機株式会社 レーダ装置及び目標測角方法
CN109212526B (zh) * 2018-10-17 2022-07-01 哈尔滨工业大学 用于高频地波雷达的分布式阵列目标角度测量方法
CN110244273B (zh) * 2019-06-18 2023-06-30 西安电子科技大学 一种基于均匀分布式阵列的目标角度估计方法
CN110988835B (zh) * 2019-11-27 2023-03-31 中国船舶重工集团公司第七二四研究所 一种分布式相参雷达测角方法
CN111427022A (zh) * 2020-05-08 2020-07-17 北京理工大学重庆创新中心 一种基于极大似然估计的阵列雷达测角方法

Also Published As

Publication number Publication date
CN113406620A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
US6661366B2 (en) Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation
CN108508423B (zh) 基于异型阵的子阵数字和差单脉冲测角方法
CN109669178B (zh) 一种星载三阵元单脉冲两维测向方法
EP1348978B1 (en) Radar processing system and method for detecting and maintaining targets
Su et al. Adaptive beamforming for nonstationary HF interference cancellation in skywave over-the-horizon radar
CN110673086A (zh) 一种基于数字阵列雷达的二维角度超分辨方法
US6720910B2 (en) Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter
EP1167995B1 (en) Matrix monopulse ratio radar processor for two target azimuth and elevation angle determination
Zhu et al. Combined constrained adaptive sum and difference beamforming in monopulse angle estimation
CN110736976B (zh) 一种任意阵形的声纳波束形成器性能估计方法
CN110208757B (zh) 一种抑制主瓣干扰的稳健自适应波束形成方法及装置
US20030184473A1 (en) Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation
CN108490428B (zh) 一种用于抗主瓣干扰的降维子阵比相跟踪测角方法
CN112230215B (zh) 未知极化参数的自适应单脉冲测向方法
CN105974377A (zh) 一种对数字阵列雷达自适应调零技术的干扰方法
CN113406620B (zh) 一种阵列分解的分布式阵列测角方法
CN109669172B (zh) 基于主瓣内强干扰抑制的弱目标方位估计方法
CN109061594B (zh) 一种强杂波环境下基于子阵空时自适应处理的测角方法
CN113820654B (zh) 基于波束域降维的s波段雷达目标低仰角doa估计方法
CN111830458B (zh) 一种平行线阵单快拍二维测向方法
Zhang et al. DOA estimation method for underwater acoustic signals based on two-dimensional power distribution (TPD) for few element array
Yu et al. Methods to combine deterministic nulling and adaptive nulling
Yeom et al. Performance analysis of beamspace MUSIC with beamforming angle
Yu Mainlobe cancellation, orthogonal nulling and product patterns
CN111175694A (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
GR01 Patent grant
GR01 Patent grant