CN111650552A - 一种基于和差协同阵构建的改进l型阵列设计方法 - Google Patents

一种基于和差协同阵构建的改进l型阵列设计方法 Download PDF

Info

Publication number
CN111650552A
CN111650552A CN202010410839.8A CN202010410839A CN111650552A CN 111650552 A CN111650552 A CN 111650552A CN 202010410839 A CN202010410839 A CN 202010410839A CN 111650552 A CN111650552 A CN 111650552A
Authority
CN
China
Prior art keywords
array
sum
cooperative
shaped
difference
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
CN202010410839.8A
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.)
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 CN202010410839.8A priority Critical patent/CN111650552A/zh
Publication of CN111650552A publication Critical patent/CN111650552A/zh
Pending legal-status Critical Current

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
    • G01S3/00Direction-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/02Direction-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 radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • G01S3/143Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种基于和差协同阵构建的改进L型阵列设计方法,通过对传统L型阵列进行稀疏化排布设计,利用协方差恢复算法能够等效出具有大范围连续段的和差协同阵以用于阵列测向,其实现步骤是:计算阵元间距基本单元,确定L型阵列的阵元数和阵元位置,将选定位置上的阵元移到以原点为中心对称的对应位置上。本发明可以在相同阵元数的条件下,有效增大阵列的自由度,提高波达方向的估计精度,提高空间目标探测能力,减少阵列运行成本。

Description

一种基于和差协同阵构建的改进L型阵列设计方法
技术领域:
本发明属于阵列信号处理技术领域,尤其涉及一种改进L型阵列的构建方法,可用于生成高自由度的 和差协同阵列。
技术背景:
波达方向(Direction of Arrival,DOA)估计是阵列信号处理领域的一个重要研究分支,它在雷达、 通信和射电天文领域有广泛应用。传统DOA估计技术多利用规则排布、阵元间隔小于等于半波长的均匀线 性阵列或均匀矩形阵列。这种均匀排布方式,保证了信号的空间谱形式唯一,空间信息无混叠,并能通过 谱函数唯一辨识空间探测目标,获得有用参数的唯一估计。然而,这些阵列结构中密集排布的阵元会导致 严重的互耦效应,进而严重影响估计精度。此外这些结构也有着极高的计算复杂度和硬件成本。
为改进上述密集天线阵列排布结构,相邻阵元间距应当进一步增大。基于此想法,许多稀疏阵列被相 继提出。其中在二维领域,L型阵列得到了广泛地应用。该阵列在xOy平面的x轴和y轴上分别以原点为 起点布置两条均匀线性阵列。一方面,该种阵列与传统的均匀矩形阵列相比极大降低了密集排布的阵元个 数;另一方面,在相同阵元个数的前提下,L型阵列通过结合差分协同阵列的概念能够获得一个无孔虚拟 阵,进而获得比均匀矩形阵列更高的自由度。但是,由于L型阵列在每条轴线上的子阵列依旧为阵元密集 排布的均匀线性阵列,可以通过设计更为稀疏的阵元排布以获得更小的阵元间互耦效应。同时,目前对于 平面阵列的虚拟阵的研究多是基于差分协同阵概念设计的,而差分协同阵的孔径扩展能力有限。联合利用 差分协同阵与其他类型的协同阵,比如由正和协同阵和负和协同阵并集组成的和协同阵,可以突破差分协 同阵的孔径扩展受限性,进一步增加自由度。
因此,如何设计一种改进的L型阵列,使其基于和差协同阵的概念获得更高自由度(即更长连续段), 并通过阵元稀疏排布降低天线耦合效应,对于实际应用有着重要的意义。
发明内容:
本发明的目的在于针对上述现有技术存在的不足,提出一种L型阵列的优化布阵方案,使其最终产生 的差分协同阵与和协同阵能够组合成具有更大虚拟孔径的协同阵,以有效提高阵列的自由度并降低耦合效 应。
为了解决上述技术问题,本发明是通过以下技术方案实现的,它包括步骤如下:
(1)用d表示阵元间距的基本单元,将其取为半波长,即d=λ/2;λ表示阵列入射信号的波长;
(2)根据阵列给定的总阵元数N,选取两个正整数Nx,Ny作为阵元参数,且N=Nx+Ny+1;
(3)在xOy坐标平面上分别构建L型阵列的两个均匀线性子阵列,每个子阵列中相邻阵元间距均为d; 其中一个子阵列分布在x轴上,由Nx+1个阵元组成,记为X,X=={(nx,0)|nx∈[0,Nx]},其中nx为整数; 另一个子阵列分布在y轴上,由Ny个阵元组成,记为Y,Y=={(0,ny)|ny∈[1,Ny]},其中ny为整数;
(4)分别将L型阵列x轴和y轴方向上两个子阵的偶数位置阵元以原点为中心移动至其中心对称位置 上,即令L性阵列中的x轴上点集{(nx,0)|nx=2,4,6,...}移动到负半轴,成为{(nx,0)|nx=-2,-4,-6,...};类似 地,y轴上点集{(0,ny)|ny=2,4,6,...}移动到负半轴,成为{(0,ny)|ny=-2,-4,-6,...};
此时时对于给定总阵元数N=Nx+Ny+1,且两个正整数参数Nx≥2和Ny≥2,改进L型阵列的阵 元位置可以解析表示为P=X1∪X2∪Y1∪Y2,若设nx和ny为整数,则有
Figure BDA0002493172730000021
此时设移位后x轴上和y轴上的子阵分别表示为X'a和Y'a,则有X'a=X1∪X2,Y'a=Y1∪Y2
(5)由改进L型阵列P的阵元生成和差协同阵列;
首先,由改进L型阵列P的阵元位置坐标相互之间两两求差,得到的一系列差值坐标的集合,构成差 分协同阵列,记为D;
其次,由改进L型阵列P的阵元位置坐标相互之间两两求矢量和,得到一系列和值坐标的集合,同时 将这些和值坐标分别取负值,得到负和值坐标的集合,上述两个集合形成的并集构成和协同阵列,记为S;
最后,由D和S组成的并集构成和差协同阵列,记为DS;
(6)利用和差协同阵列DS所包含的具有最多连续虚拟阵元的均匀矩形面阵,可实施各类波达方向估计 算法,对空间信号的来波方向进行准确估计。
进一步地,生成的差分协同阵列D在xOy第二和第四象限各有一块连续均匀矩形面阵,其虚拟阵元位 置的数学表达式为:{(nx,ny)|nx∈[-Nx,0],ny∈[1,Ny]}∪{(nx,ny)|nx∈[0,Nx],ny∈[-Ny,-1]}。
进一步地,生成的和协同阵列S在xOy第一和第三象限各有一块连续均匀矩形面阵,其虚拟阵元位置 的数学表达式为:{(nx,ny)|nx∈[0,Nx],ny∈[1,Ny]}∪{(nx,ny)|nx∈[-Nx,0],ny∈[-Ny,-1]}。
进一步地,生成的和差协同阵列DS在xOy平面中心有一块连续均匀矩形面阵,其虚拟阵元位置的数 学表达式为:{(nx,ny)|nx∈[-Nx,Nx],ny∈[-Ny,Ny]},该连续范围所提供的阵列自由度为 (2Nx+1)×(2Ny+1)。
进一步地,在上述结构中,间距为d的天线阵列有2对,分别是(0,0)和(0,1),(0,0)和(1,0);间距为
Figure BDA0002493172730000031
的天线阵列有1对,为(1,0)和(0,1)。
本发明的有益效果如下:
根据本发明设计的改进L型阵列所构建的和差协同阵列,包含两个增大了稀疏程度的线性子阵列,同 传统L型阵列或其他二维阵列所构建的差分协同阵列相比,本发明构建的和差协同阵列在相同阵元数目下, 可以提供更高的阵列自由度;反之,在提供相同阵列自由度情况下,本发明的设计方法能够极大地减少对 阵元数目的需求。
附图说明
图1是本发明步骤(4)构建的改进L型阵列的一种示例图,其中Nx=15,Ny=15。
图2是本发明步骤(4)构建的改进L型阵列的和差协同阵列,此时Nx=15,Ny=15。其中蓝色×代 表阵列生成的差分协同阵列位置,红色○代表阵列生成的和协同阵列位置。
具体实施方式
下面结合附图及实例对本发明的内容做进一步的说明。
以给定入射信号频率20GHz,阵元总数N=31为例,对本发明的改进L型阵列(如图1所示)构建进 行详细说明。
步骤(1):计算入射信号波长λ=c/f=3.0×108/20×109=1.5cm;用d表示阵元间距的基本单元, 将其取为半波长,即d=λ/2=7.5mm。
步骤(2):根据阵列给定的总阵元数N=31,选取两个正整数Nx,Ny作为阵元参数,要求 N=Nx+Ny+1;此实例选取Nx=15,Ny=15。
步骤(3):在xOy坐标平面上分别构建L型阵列的两个均匀线性子阵列,每个子阵列中相邻阵元间距 均为d;其中一个子阵列分布在x轴上,由Nx+1=16个阵元组成,记为X,X=={(nx,0)|nx∈[0,15]},其 中nx为整数;另一个子阵列分布在y轴上,由Ny=15个阵元组成,记为Y,Y=={(0,ny)|ny∈[1,15]},其中 ny为整数;
步骤(4):分别将L型阵列x轴和y轴方向上两个子阵的偶数位置阵元以原点为中心移动至其中心对 称位置上。即令L性阵列中的x轴上点集{(nx,0)|nx=2,4,6,...}移动到负半轴,成为{(nx,0)|nx=-2,-4,-6,...}; 类似地,y轴上点集{(0,ny)|ny=2,4,6,...}移动到负半轴,成为{(0,ny)|ny=-2,-4,-6,...}。此时选取 Mx=My=15,子阵X1,X2,Y1,Y2中所选取阵元的归一化位置(省略阵元间距d)可表示为:
Figure BDA0002493172730000041
此时设移位后x轴上和y轴上的子阵分别表示为X'a和Y'a,则有
Figure RE-GDA0002615430410000042
步骤(5):由子阵列X'a和子阵列Y'a的阵元生成和差协同阵列;
首先,由子阵列X'a和子阵列Y'a的阵元位置坐标相互之间两两求差,得到的一系列差值坐标的集合, 构成差分协同阵列,记为D;生成的差分协同阵列D在xOy第二和第四象限各有一块均匀矩形面阵,其虚 拟阵元位置的数学表达式为:{(nx,ny)|nx∈[-15,0],ny∈[1,15]}∪{(nx,ny)|nx∈[0,15],ny∈[-15,-1]}
其次,由子阵列X'a和子阵列Y'a的阵元位置坐标相互之间两两求和,得到的一系列和值坐标集合,同 时将这些和值坐标分别取负值,得到负和值坐标的集合,上述两个集合形成的并集构成和协同阵列,记为 S;生成的和协同阵列S在xOy第一和第三象限各有一块均匀矩形面阵,其虚拟阵元位置的数学表达式为: {(nx,ny)|nx∈[0,15],ny∈[1,15]}∪{(nx,ny)|nx∈[-15,0],ny∈[-15,-1]}。
最后,由D和S组成的并集构成和差协同阵列,记为DS;
如图2实例可见DS中阵元的位置,其中在{(nx,ny)|nx∈[-15,15],ny∈[-15,15]}的范围内形成了一 个无孔的面阵。
进一步地,在此结构中,间距为d的天线阵列有2对,间距为
Figure BDA0002493172730000043
的天线阵列有1对。

Claims (5)

1.一种基于和差协同阵构建的改进L型阵列设计方法,其特征在于,包括步骤如下:
(1)用d表示阵元间距的基本单元,将其取为半波长,即d=λ/2;λ表示阵列入射信号的波长;
(2)根据阵列给定的总阵元数N,选取两个正整数Nx,Ny作为阵元参数,则有N=Nx+Ny+1;
(3)在xOy坐标平面上分别构建L型阵列的两个均匀线性子阵列,每个子阵列中相邻阵元间距均为d;其中一个子阵列分布在x轴上,由Nx+1个阵元组成,记为X,X=={(nx,0)|nx∈[0,Nx]},其中nx为整数;另一个子阵列分布在y轴上,由Ny个阵元组成,记为Y,Y=={(0,ny)|ny∈[1,Ny]},其中ny为整数;
(4)分别将L型阵列x轴和y轴方向上两个子阵的偶数位置阵元以原点为中心移动至其中心对称位置上,即令L性阵列中的x轴上点集{(nx,0)|nx=2,4,6,...}移动到负半轴,成为{(nx,0)|nx=-2,-4,-6,...};类似地,y轴上点集{(0,ny)|ny=2,4,6,...}移动到负半轴,成为{(0,ny)|ny=-2,-4,-6,...};
此时对于给定总阵元数N=Nx+Ny+1,且两个正整数参数为Nx≥2和Ny≥2,改进L型阵列的阵元位置可以解析表示为P=X1∪X2∪Y1∪Y2,若设nx和ny为整数,则有
Figure FDA0002493172720000011
此时设移位后x轴上和y轴上的子阵分别表示为X'a和Y'a,则有X'a=X1∪X2,Y'a=Y1∪Y2
(5)由改进L型阵列P的阵元生成和差协同阵列;
首先,由改进L型阵列P的阵元位置坐标相互之间两两求矢量差,得到的一系列差值坐标的集合构成差分协同阵列,记为D;
其次,由改进L型阵列P的阵元位置坐标相互之间两两求矢量和,得到一系列和值坐标的集合,同时将这些和值坐标分别取负值,得到负和值坐标的集合,上述两个集合形成的并集构成和协同阵列,记为S;
最后,由D和S组成的并集构成和差协同阵列,记为DS;
(6)利用和差协同阵列DS所包含的具有最多连续虚拟阵元的均匀矩形面阵,可实施各类波达方向估计算法,对空间信号的来波方向进行准确估计。
2.根据权利要求1所述的一种基于和差协同阵构建的改进L型阵列设计方法,其特征在于,生成的差分协同阵列D在xOy第二和第四象限各有一块连续均匀矩形面阵,其虚拟阵元位置的数学表达式为:{(nx,ny)|nx∈[-Nx,0],ny∈[1,Ny]}∪{(nx,ny)|nx∈[0,Nx],ny∈[-Ny,-1]}。
3.根据权利要求1所述的一种基于和差协同阵构建的改进L型阵列设计方法,其特征在于,生成的和协同阵列S在xOy第一和第三象限各有一块连续均匀矩形面阵,其虚拟阵元位置的数学表达式为:{(nx,ny)|nx∈[0,Nx],ny∈[1,Ny]}∪{(nx,ny)|nx∈[-Nx,0],ny∈[-Ny,-1]}。
4.根据权利要求1所述的一种基于和差协同阵构建的改进L型阵列设计方法,其特征在于,生成的和差协同阵列DS在xOy平面中心有一块连续均匀矩形面阵,其虚拟阵元位置的数学表达式为:{(nx,ny)|nx∈[-Nx,Nx],ny∈[-Ny,Ny]},该连续范围所提供的阵列自由度为(2Nx+1)×(2Ny+1)。
5.根据权利要求1所述的一种基于和差协同阵构建的改进L型阵列设计方法,其特征在于,在该结构中,间距为d的天线阵列有2对,分别是(0,0)和(0,1),(0,0)和(1,0);间距为
Figure FDA0002493172720000021
的天线阵列有1对,为(1,0)和(0,1)。
CN202010410839.8A 2020-05-15 2020-05-15 一种基于和差协同阵构建的改进l型阵列设计方法 Pending CN111650552A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010410839.8A CN111650552A (zh) 2020-05-15 2020-05-15 一种基于和差协同阵构建的改进l型阵列设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010410839.8A CN111650552A (zh) 2020-05-15 2020-05-15 一种基于和差协同阵构建的改进l型阵列设计方法

Publications (1)

Publication Number Publication Date
CN111650552A true CN111650552A (zh) 2020-09-11

Family

ID=72349235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010410839.8A Pending CN111650552A (zh) 2020-05-15 2020-05-15 一种基于和差协同阵构建的改进l型阵列设计方法

Country Status (1)

Country Link
CN (1) CN111650552A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113435027A (zh) * 2021-06-23 2021-09-24 哈尔滨工程大学 一种高自由度低耦合稀疏线性阵列排布方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140152504A1 (en) * 2012-12-02 2014-06-05 Khalifa University of Science, Technology & Research (KUSTAR) Method and system for measuring direction of arrival of wireless signal using circular array displacement
CN104280732A (zh) * 2014-09-16 2015-01-14 电子科技大学 一种基于等效协同阵的穿墙雷达建筑布局成像方法
CN109581276A (zh) * 2018-11-26 2019-04-05 电子科技大学 一种基于求和求差嵌套阵的doa估计方法
CN109917329A (zh) * 2019-04-15 2019-06-21 南京邮电大学 一种基于协方差匹配准则的l型阵列波达方向估计方法
CN109917328A (zh) * 2019-04-15 2019-06-21 南京邮电大学 一种基于原子范数最小化的l型阵列波达方向估计方法
CN109932680A (zh) * 2019-04-04 2019-06-25 哈尔滨工程大学 一种基于平移互质阵列的非圆信号波达方向估计方法
CN110736959A (zh) * 2019-10-25 2020-01-31 北京理工大学 一种基于和差协同阵构建的平面互质阵列设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140152504A1 (en) * 2012-12-02 2014-06-05 Khalifa University of Science, Technology & Research (KUSTAR) Method and system for measuring direction of arrival of wireless signal using circular array displacement
CN104280732A (zh) * 2014-09-16 2015-01-14 电子科技大学 一种基于等效协同阵的穿墙雷达建筑布局成像方法
CN109581276A (zh) * 2018-11-26 2019-04-05 电子科技大学 一种基于求和求差嵌套阵的doa估计方法
CN109932680A (zh) * 2019-04-04 2019-06-25 哈尔滨工程大学 一种基于平移互质阵列的非圆信号波达方向估计方法
CN109917329A (zh) * 2019-04-15 2019-06-21 南京邮电大学 一种基于协方差匹配准则的l型阵列波达方向估计方法
CN109917328A (zh) * 2019-04-15 2019-06-21 南京邮电大学 一种基于原子范数最小化的l型阵列波达方向估计方法
CN110736959A (zh) * 2019-10-25 2020-01-31 北京理工大学 一种基于和差协同阵构建的平面互质阵列设计方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113435027A (zh) * 2021-06-23 2021-09-24 哈尔滨工程大学 一种高自由度低耦合稀疏线性阵列排布方法

Similar Documents

Publication Publication Date Title
CN107102291B (zh) 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法
CN107329108B (zh) 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法
CN107315160B (zh) 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法
CN108344967B (zh) 基于互质面阵的二维波达方向快速估计方法
CN106407723B (zh) 面向低副瓣的稀疏排布阵列天线激励电流幅度的确定方法
CN107015190A (zh) 基于虚拟阵列协方差矩阵稀疏重建的互质阵列波达方向估计方法
CN112115595A (zh) 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法
WO2021068496A1 (zh) 基于结构化虚拟域张量信号处理的互质面阵二维波达方向估计方法
CN107703478B (zh) 基于互相关矩阵的扩展孔径二维doa估计方法
CN110954861A (zh) 一种基于增强型嵌套阵列的doa估计方法
CN108872930B (zh) 扩展孔径二维联合对角化doa估计方法
CN111983554A (zh) 非均匀l阵下的高精度二维doa估计
Si et al. Two-dimensional DOA estimation via a novel sparse array consisting of coprime and nested subarrays
CN113189538A (zh) 一种基于互质稀疏排列的三元阵列及其空间谱估计方法
CN111650552A (zh) 一种基于和差协同阵构建的改进l型阵列设计方法
CN112731277A (zh) 一种基于改进混合基线的任意阵列干涉仪测向方法
CN110232228B (zh) 一种多联骨牌非规则子阵阵面优选设计方法
CN112711000B (zh) 基于最小化准则的电磁矢量互质面阵张量功率谱估计方法
CN113075649B (zh) 一种适用于分布式网络化雷达的信号级直接定位方法
CN110726967B (zh) 用于一维阵列测向的单边稀疏嵌套阵设计方法
CN111323742A (zh) 一种基于曲线阵列的相位干涉仪及其测向方法
CN110736959B (zh) 一种基于和差协同阵构建的平面互质阵列设计方法
CN115616478A (zh) 一种基于混合远近场的新型对称稀疏阵列的设计方法
CN112710983B (zh) 基于乘性张量波束扫描的电磁矢量互质面阵多维参数估计方法
CN112100837A (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200911

WD01 Invention patent application deemed withdrawn after publication