CN106785486A - 一种广义互质面阵天线结构及角度估计方法 - Google Patents
一种广义互质面阵天线结构及角度估计方法 Download PDFInfo
- Publication number
- CN106785486A CN106785486A CN201710014939.7A CN201710014939A CN106785486A CN 106785486 A CN106785486 A CN 106785486A CN 201710014939 A CN201710014939 A CN 201710014939A CN 106785486 A CN106785486 A CN 106785486A
- Authority
- CN
- China
- Prior art keywords
- array
- relatively prime
- antenna
- broad sense
- axis direction
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q21/00—Antenna arrays or systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种广义互质面阵(Generalized coprime planar array,GCPA)天线结构,该结构由两个矩形均匀面阵组成,天线数分别为N1×M1和N2×M2,其中N1,N2为x轴方向上的天线数,M1、M2为y轴方向上的天线数,N1、N2和M1、M2为两组互质数;两个子阵只在原点处有一个阵元重合。本发明天线结构进行角度估计的方法为:首先,根据MUSIC或ESPRIT算法得到存在模糊的角度估计值;然后,利用互质特性寻找模糊值重合的角度估计值,即为对应理论值的真实值。本发明的天线结构优势及优点在于拓展了天线孔径,天线阵列布局更加灵活,且能获得相对较高的自由度。两个子阵保持均匀的特性,经典的信号到达角估计算法能够直接应用,利用互质特性,角度模糊问题可以被消除。
Description
技术领域
本发明涉及一种广义互质面阵天线结构及角度估计方法,属于阵列天线布局的领域。
背景技术
传统天线布局受到半波长的限制,即阵元间距要小于半波长,因此,在天线数一定的情况下,天线的孔径受到了限制。
当前,一种被称为互质阵的天线布局受到了广泛得关注。其阵元间距突破了半波长得限制,使得天线孔径得到极大得拓展,能获得很好的角度估计性能的提升。但是多数研究还停留在一维线阵的研究,不能满足现实需求。由两个方形面阵构成的互质面阵的可行性已经得到证明,但是其自由度受到较少阵元数的子阵的限制。
发明内容
本发明所要解决的技术问题是提供一种广义互质面阵(Generalized coprimeplanar array,GCPA)天线结构,该GCPA结构,由两个矩形子面阵构成,且子阵列的阵元数为总阵元数的一半,从而使得GCPA能够获得更好的自由度;并且在总的天线数相同的情况,我们证明GCPA能够获得比方形互质面阵更好的角度估计性能。本发明的广义互质面阵具有更灵活的天线布局,更高的自由度,且更好的角度估计性能的优势。可以应用与无线通信、声呐、定位等领域中。
本发明为解决上述技术问题采用以下技术方案:
一方面,本发明提供一种广义互质面阵天线结构,该结构由两个矩形均匀面阵组成,两个矩形均匀面阵的天线数分别为N1×M1和N2×M2,其中,N1、N2为x轴方向上的天线数,M1、M2为y轴方向上的天线数,且N1、N2和M1、M2为两组互质数;两个矩形均匀面阵仅在原点处有一个阵元重合。
作为本发明的进一步优化方案,该结构的总阵元数为Tgcpa=N1M1+N2M2-1。
作为本发明的进一步优化方案,阵元数为N1×M1的矩形均匀面阵:x轴方向上的阵元间距为dx1=N2λ/2,y轴方向上的阵元间距为dy1=M2λ/2,其中,λ为载波波长。
作为本发明的进一步优化方案,阵元数为N2×M2的矩形均匀面阵:x轴方向上的阵元间距为dx2=N1λ/2,y轴方向上的阵元间距为dy2=M1λ/2,其中,λ为载波波长。
另一方面,本发明还提供一种上述广义互质面阵天线结构进行角度估计的方法,该方法具体为:首先,根据MUSIC或ESPRIT算法得到存在模糊的角度估计值;然后,利用互质特性寻找模糊值重合的角度估计值,即为对应理论值的真实值。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
①充分利用了阵元,从而获得更高的自由度;
②具有天线布局灵活的特点;
③该天线布局能获得更好的角度估计性能。
附图说明
图1是广义互质面阵GCPA天线结构示意图;
图2是当9个信号入射到GCPA结构时,利用MUSIC算法得到的角度估计散点图;
图3是ESPRIT和MUSIC算法在GCPA和方形互质面阵下的角度估计性能对比示意图,其中,(a)为仰角的性能对比图,(b)为方位角的性能对比图;
图4是ESPRIT和MUSIC算法在GCPA和方形互质面阵下的角度估计性能随着快拍数变化的对比示意图,其中,(a)为仰角的性能对比图,(b)为方位角的性能对比图。
具体实施方式
下面结合附图以及具体实施例对本发明的技术方案做进一步的详细说明:
本发明设计一种广义互质面阵(Generalized coprime planar array,GCPA)天线结构,该结构由两个矩形均匀面阵组成,天线数分别为N1×M1和N2×M2,其中N1、N2为x轴方向上的天线数,M1、M2为y轴方向上的天线数,两个子阵只在原点处有一个阵元重合,所以总的阵元数为Tgcpa=N1M1+N2M2-1。N1,N2和M1,M2为两组互质数。阵元数为N1×M1的子阵,其x轴方向上的阵元间距为dx1=N2λ/2,y轴方向上的阵元间距为dy1=M2λ/2,其中λ为载波波长;阵元数为N2×M2的子阵,其x轴方向上的阵元间距为dx2=N1λ/2,y轴方向上的阵元间距为dy2=M1λ/2,其中λ为载波波长。
一、数据模型
如图1所示的是一个GCPA天线结构的例子,其中,N1=5,M1=3,N2=6,M2=4。
假设K个来自{(θk,φk)|k=1,2,…,K}的非相干窄带信号入射到GCPA上,其中,θk和φk分别为第k个信号的仰角与方位角(K<min{N1M1,N2M2},θk∈(0,90°),φk∈(0,180°)),定义uk=sinθk cosφk∈(-1,1)和vk=sinθk sinφk∈(0,1)。为了简化,在接下来的部分用Ni×Mi(i=1,2)来表示GCPA两个子阵的阵元数。
接收信号可以表示为
Xi=AiS+Ni (i=1,2)
其中,S=[s1,s2,…,sK]T为信号矩阵,sk=[sk(1),sk(2),…,sk(l),…,sk(L)],L为快拍数,sk(l)为对第k个信号的l次采样,l=1,2,…,L;为加性高斯白噪声,均值为零、方差为 为阵元数为Ni×Mi(i=1,2)的第i个子阵的方向矩阵,并且axi(uk)和ayi(vk)为第i个子阵在x方向和y轴方向的方向矢量,可以表示为
axi(uk)=[1,exp(-j2πdxiuk/λ),…,exp(-j2π(Ni-1)dxiuk/λ)]T
ayi(vk)=[1,exp(-j2πdyivk/λ),…,exp(-j2π(Mi-1)dyivk/λ)]T。
实际工程中,由于采样都是在有限快拍数下进行,信号矩阵的协方差矩阵由计算得到。特征分解后可以得到其中,Esi为由对应K个最大特征值的特征矢量构成的信号子空间,Dsi是由K个最大特征值构成的对角矩阵,Eni为其余特征矢量构成的噪声子空间,Dni是由其它所有特征值构成的对角矩阵。
二、自由度(Degree Of Freedom,DOF)
在相同天线数的情况下,GCPA结构可以获得比方形互质面阵更高的自由度。
假设方形互质面阵的两个子阵天线数分别为P1×P1和P2×P2,其中P1,P2为互质数,天线总数为T=P1×P1+P2×P2-1。假设GCPA结构的两个子阵天线数分别为N1×M1和N2×M2,其中N1,N2和M1,M2为两组互质数,其天线总数也为T=P1×P1+P2×P2-1=N1×M1+N2×M2-1。不失一般性,假设P1<P2,N1×M1≤N2×M2,N2×M2=N1×M1+δ,其中δ≥0是一个整数。那么,方形互质面阵的自由度为DOFcpa=P1×P1-1,GCPA的自由度为DOFgcpa=N1×M1-1。
T=N1×M1+N2×M2-1=2N1×M1+δ-1
所以,GCPA结构的最大自由度为当两个子阵的天线数分别为和
其中因为P1<P2,所以GCPA结构最多可以比方形互质面阵多的自由度为
三、克拉美罗界
定义GCPA的方向矩阵为:
其中,Axi=[axi(u1),axi(u2),…,axi(uK)],axi(uk)=[1,exp(-j2πdxiuk/λ),…,exp(-j2π(Ni-1)dxiuk/λ)]T,Ayi=[ayi(v1),ayi(v2),…,ayi(vK)],ayi(vk)=[1,exp(-j2πdyivk/λ),…,exp(-j2π(Mi-1)dyivk/λ)]T,i=1,2;At由Ax2的第二行到最后一行构成;Dm(·)为取矩阵的第m行构成一个对角矩阵;信号协方差矩阵可以由Ps=SSH/L得到。
角度估计的克拉美罗界可由下式计算得到:
其中,ak为A的第k列,为N1M1+N2M2-1维的单位矩阵,k=1,2,…,K,
四、角度估计方法
1、经典二维角度估计算法
根据噪声子空间与信号子空间正交的特性,我们可以构造二维多重信号分类(Multiple Signal Classification,MUSIC)算法的空间谱函数为
其中方向矢量Eni为噪声子空间;
axi(u)=[1,exp(-j2πdxiu/λ),…,exp(-j2π(Ni-1)dxiu/λ)]T,
ayi(v)=[1,exp(-j2πdyiv/λ),…,exp(-j2π(Mi-1)dyiv/λ)]T。
对谱函数进行二维谱峰搜索,我们可以得到存在角度模糊的波达角参数。另外,由于GCPA的两个子阵的方向矩阵依然具有范德蒙矩阵特性,所以经典的ESPRIT算法也适用于GCPA阵列。
2、模糊消除
假设只有一个来自(θt,φt)方向的信号入射到GCPA阵列上。我们假设(θa,φa)为(θt,φt)的一个模糊角度。由于以自然常数为基的指数函数具有2π周期性,我们可以得到真实值与模糊值之间的关系:
2πdxi(ut-ua)/λ=2kuiπ
2πdyi(vt-va)/λ=2kviπ
其中,ut=sinθt cosφt,vt=sinθt sinφt,ua=sinθa cosφa,va=sinθa sinφa,为整数,dxi=Njλ/2,dyi=Mjλ/2(i,j=1,2;i≠j)。
根据假设,我们可以得到ua∈(-1,1),va∈(0,1),这两个参数不仅需要分别满足要求,还需要满足
根据互质数的特性,我们可以消除角度模糊问题,得到角度估计的真实值。我们假设(θp,φp)为理论角度(θt,φt)在两个子阵中出现的模糊角度值。
ut-up=2kuiπ/Nj,vt-vp=2kviπ/Mj
其中,up=sinθp cosφp,vp=sinθp sinφp,kui∈(-Nj,Nj),kvi∈(-Mj,Mj)(i,j=1,2;i≠j)。
ku1/N2=ku2/N1,kv1/M2=kv2/M1
因为N1,N2和M1,M2为两组互质数,所以只有ku1=ku2=0和kv1=kv2=0可以使得上式成立,这意味着(θp,φp)即为对应(θt,φt)的真实值。
图2当9个信号入射到GCPA结构时,利用MUSIC算法得到的角度估计散点图,其中方形互质面阵的两个子阵的阵元数分别为3×3和5×5,此时基于该阵列的子空间算法已经失效;GCPA的两个子阵的阵元数分别为4×4和3×5,虽然GCPA比方形互质面阵少3根天线,但是并不影响结果,由图可以看到GCPA此时依然能够很好的工作(L=100,SNR=8dB);
图3是ESPRIT和MUSIC算法在GCPA和方形互质面阵下的角度估计性能对比,其中方形互质面阵的两个子阵阵元数分别为5×5和8×8,GCPA结构的两个子阵的阵元数为5×9和4×11,两个阵列的天线数相同,都为88根,(θ1,φ1)=(20°,20°),(θ2,φ2)=(40°,40°),L=100,其中CPA代表方形互质面阵。
图4是ESPRIT和MUSIC算法在GCPA和方形互质面阵下的角度估计性能随着快拍数变化的对比,其中方形互质面阵的两个子阵阵元数分别为5×5和8×8,GCPA结构的两个子阵的阵元数为5×9和4×11,两个阵列的天线数相同,都为88根,SNR=0dB,(θ1,φ1)=(20°,20°),(θ2,φ2)=(40°,40°),其中CPA代表方形互质面阵。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种广义互质面阵天线结构,其特征在于,该结构由两个矩形均匀面阵组成,两个矩形均匀面阵的天线数分别为N1×M1和N2×M2,其中,N1、N2为x轴方向上的天线数,M1、M2为y轴方向上的天线数,且N1、N2和M1、M2为两组互质数;两个矩形均匀面阵仅在原点处有一个阵元重合。
2.根据权利要求1所述的一种广义互质面阵天线结构,其特征在于,该结构的总阵元数为Tgcpa=N1M1+N2M2-1。
3.根据权利要求1所述的一种广义互质面阵天线结构,其特征在于,阵元数为N1×M1的矩形均匀面阵:x轴方向上的阵元间距为dx1=N2λ/2,y轴方向上的阵元间距为dy1=M2λ/2,其中,λ为载波波长。
4.根据权利要求1所述的一种广义互质面阵天线结构,其特征在于,阵元数为N2×M2的矩形均匀面阵:x轴方向上的阵元间距为dx2=N1λ/2,y轴方向上的阵元间距为dy2=M1λ/2,其中,λ为载波波长。
5.一种如权利要求1所述的广义互质面阵天线结构进行角度估计的方法,其特征在于,该方法具体为:首先,根据MUSIC或ESPRIT算法得到存在模糊的角度估计值;然后,利用互质特性寻找模糊值重合的角度估计值,即为对应理论值的真实值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710014939.7A CN106785486B (zh) | 2017-01-09 | 2017-01-09 | 一种广义互质面阵天线结构 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710014939.7A CN106785486B (zh) | 2017-01-09 | 2017-01-09 | 一种广义互质面阵天线结构 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106785486A true CN106785486A (zh) | 2017-05-31 |
CN106785486B CN106785486B (zh) | 2020-02-18 |
Family
ID=58948619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710014939.7A Active CN106785486B (zh) | 2017-01-09 | 2017-01-09 | 一种广义互质面阵天线结构 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106785486B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107611624A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 低旁瓣的基于互质思想的分子阵布阵方法 |
CN109239651A (zh) * | 2018-07-24 | 2019-01-18 | 南京航空航天大学 | 互质面阵下的二维doa跟踪方法 |
CN109327249A (zh) * | 2018-11-13 | 2019-02-12 | 钟祥博谦信息科技有限公司 | Mimo系统中的到达角估计方法 |
CN109738853A (zh) * | 2018-11-21 | 2019-05-10 | 南京航空航天大学 | 一种电磁矢量互质阵基于旋转不变性的角度与极化估计方法 |
CN110736959A (zh) * | 2019-10-25 | 2020-01-31 | 北京理工大学 | 一种基于和差协同阵构建的平面互质阵列设计方法 |
JP2020085461A (ja) * | 2018-11-15 | 2020-06-04 | 株式会社東芝 | 電子機器及びイメージング方法 |
CN111239679A (zh) * | 2020-02-12 | 2020-06-05 | 南京航空航天大学 | 一种用于互质面阵下相干信源doa估计的方法 |
CN113156362A (zh) * | 2021-03-15 | 2021-07-23 | 北京邮电大学 | 波达方向的确定方法、装置及信号获取方法、装置 |
CN113484821A (zh) * | 2021-07-06 | 2021-10-08 | 北京邮电大学 | 一种新型虚拟阵列结构及其doa估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912791A (zh) * | 2016-04-15 | 2016-08-31 | 山东农业大学 | 虚拟互质阵列中基于局部搜索的doa估计方法 |
-
2017
- 2017-01-09 CN CN201710014939.7A patent/CN106785486B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105912791A (zh) * | 2016-04-15 | 2016-08-31 | 山东农业大学 | 虚拟互质阵列中基于局部搜索的doa估计方法 |
Non-Patent Citations (2)
Title |
---|
QIHUI WU等: "Two-Dimensional Direction-of-Arrival Estimation for Co-Prime Planar Arrays: A Partial Spectral Search Approach", 《IEEE SENSORS JOURNAL》 * |
WANG ZHENG等: "Generalized Coprime Planar Array Geometry for 2-D DOA Estimation", 《IEEE COMMUNICATIONS LETTERS》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107611624A (zh) * | 2017-08-24 | 2018-01-19 | 电子科技大学 | 低旁瓣的基于互质思想的分子阵布阵方法 |
CN109239651A (zh) * | 2018-07-24 | 2019-01-18 | 南京航空航天大学 | 互质面阵下的二维doa跟踪方法 |
CN109239651B (zh) * | 2018-07-24 | 2023-06-20 | 南京航空航天大学 | 互质面阵下的二维doa跟踪方法 |
CN109327249A (zh) * | 2018-11-13 | 2019-02-12 | 钟祥博谦信息科技有限公司 | Mimo系统中的到达角估计方法 |
JP7086824B2 (ja) | 2018-11-15 | 2022-06-20 | 株式会社東芝 | 電子機器及びイメージング方法 |
JP2020085461A (ja) * | 2018-11-15 | 2020-06-04 | 株式会社東芝 | 電子機器及びイメージング方法 |
US11677159B2 (en) | 2018-11-15 | 2023-06-13 | Kabushiki Kaisha Toshiba | Electronic apparatus and imaging method |
CN109738853A (zh) * | 2018-11-21 | 2019-05-10 | 南京航空航天大学 | 一种电磁矢量互质阵基于旋转不变性的角度与极化估计方法 |
CN110736959A (zh) * | 2019-10-25 | 2020-01-31 | 北京理工大学 | 一种基于和差协同阵构建的平面互质阵列设计方法 |
CN111239679A (zh) * | 2020-02-12 | 2020-06-05 | 南京航空航天大学 | 一种用于互质面阵下相干信源doa估计的方法 |
CN113156362A (zh) * | 2021-03-15 | 2021-07-23 | 北京邮电大学 | 波达方向的确定方法、装置及信号获取方法、装置 |
CN113156362B (zh) * | 2021-03-15 | 2024-01-09 | 北京邮电大学 | 波达方向的确定方法、装置及信号获取方法、装置 |
CN113484821A (zh) * | 2021-07-06 | 2021-10-08 | 北京邮电大学 | 一种新型虚拟阵列结构及其doa估计方法 |
CN113484821B (zh) * | 2021-07-06 | 2024-04-12 | 北京邮电大学 | 一种新型虚拟阵列结构及其doa估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106785486B (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106785486A (zh) | 一种广义互质面阵天线结构及角度估计方法 | |
CN108120967B (zh) | 一种平面阵列doa估计方法及设备 | |
Yang et al. | DOA estimation with sub-array divided technique and interporlated esprit algorithm on a cylindrical conformal array antenna | |
CN108344967A (zh) | 基于互质面阵的二维波达方向快速估计方法 | |
Yang et al. | A unified array geometry composed of multiple identical subarrays with hole-free difference coarrays for underdetermined DOA estimation | |
Wei et al. | Pair-matching method by signal covariance matrices for 2D-DOA estimation | |
Li et al. | Two-dimensional direction finding with parallel nested arrays using DOA matrix method | |
CN110244273A (zh) | 一种基于均匀分布式阵列的目标角度估计方法 | |
IHEDRANE et al. | Direction of arrival estimation using MUSIC, ESPRIT and maximum-likelihood algorithms for antenna arrays | |
CN106443571A (zh) | 一种保角变换旋转不变子空间波达方向估计方法 | |
CN102142609A (zh) | 具有低旁瓣特性的子阵级自适应数字波束形成器 | |
CN110531312A (zh) | 一种基于稀疏对称阵列的doa估计方法和系统 | |
Yue et al. | Closed-form two-dimensional DOA and polarization joint estimation using parallel non-collocated sparse COLD arrays | |
Nechaev et al. | Evaluation of the influence of directivity factor of directive elements of conformal antenna arrays on the performances of azimuth-elevation DOA estimation | |
Zhang et al. | A novel high degree of freedom sparse array with displaced multistage cascade subarrays | |
CN104656055A (zh) | 一种基于大规模多天线系统的单一信号到达角估计方法 | |
Peng et al. | DOA estimation using MUSIC algorithm on a cylinderical conformal array | |
Peshkov et al. | The estimation of radio direction-finding performance in volume antenna arrays with directive radiators by music method | |
Nechaev et al. | Accuracy researching of direction-of-arrival estimation via music for circular, octagonal, hexagonal and rectangular antenna arrays | |
Liu et al. | Grade nested array with increased degrees of freedom for quasi-stationary signals | |
Liu et al. | Triple two-level nested array with improved degrees of freedom | |
Bermudez et al. | Simulation study on DOA estimation using ESPRIT algorithm | |
Liu et al. | DOA estimation with compressing sensing based array interpolation technique in multi-carrier frequency MIMO HFSWR | |
Yang et al. | Fast 2-d DOA and polarization estimation using arbitrary conformal antenna array | |
Taillefer et al. | ESPRIT algorithm for a seven-element regular-hexagonal shaped ESPAR antenna |
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 |