CN109597020A - 一种使用互质线阵进行低复杂度角度估计的方法 - Google Patents
一种使用互质线阵进行低复杂度角度估计的方法 Download PDFInfo
- Publication number
- CN109597020A CN109597020A CN201910083523.XA CN201910083523A CN109597020A CN 109597020 A CN109597020 A CN 109597020A CN 201910083523 A CN201910083523 A CN 201910083523A CN 109597020 A CN109597020 A CN 109597020A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- relatively prime
- linear array
- virtual
- 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
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
- 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/02—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 radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
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
本发明公开了一种使用互质线阵进行低复杂度角度估计的方法,首先,对阵列接收信号的协方差矩阵进行向量化处理,删除冗余行并截取连续阵元部分,得到阵元间距为λ/2的虚拟均匀线阵的接收信号;然后,对虚拟接收信号进行离散傅里叶变换(Discrete Fourier Transform,DFT)并找出K个最大峰值(K为信源数)的位置,由此可得波达方向的初始估计;最后,对经过相位旋转的虚拟接收信号进行DFT变换,在小范围内对旋转相位η进行搜索,分别找出使得K个峰值收缩至最大的旋转相位,并由此得到波达方向的精确估计。本发明将DFT变换与互质阵列的DOA估计问题相结合,不需要对虚拟接收信号进行解相关处理,能够获得相对较高的空间自由度,同时具有较低的计算复杂度和较高的估计精度。
Description
技术领域
本发明涉及空间谱估计领域,特别是涉及一种使用互质线阵进行低复杂度角度估计的方法。
背景技术
当前,一种被称为互质阵的阵列天线布局得到了关注,该阵列结构突破了传统天线阵阵元间距半波长的限制,使得天线孔径得到极大的扩展,能够获得角度估计性能的提升。
使用互质线阵进行角度估计的方法有解模糊和虚拟化两类。其中,解模糊方法能够获得比均匀线阵更好的估计性能,但使用两个子阵分别进行角度估计的方法使得空间自由度大大降低。而虚拟化方法则能够获得比实际阵元数更高的空间自由度,极大地增加了可探测的信源数。由于虚拟化后等效的虚拟信号是单快拍的相干信号,基于虚拟化的角度估计算法通常需要使用空间平滑技术对接收信号进行解相关处理,这一过程降低了阵列的空间自由度。
发明内容
本发明所要解决的技术问题是提供一种使用互质线阵进行角度估计的方法,该方法能够获得较大的空间自由度,同时具有较低的计算复杂度和较高的估计精度。可以应用于无线通信、声呐、定位等领域中。
本发明为解决上述技术问题采用以下技术方案:提出一种结合虚拟化和离散傅里叶变换(Discrete Fourier Transform,DFT)的互质线阵角度估计方法,称为DFT算法,该方法具体为:
首先,对互质线阵接收信号的协方差矩阵进行向量化处理,然后将获得的向量按相位进行排序并删除冗余行,截取此时向量的中间2MN+2M-1行元素作为虚拟阵元接收信号。
然后,对虚拟接收信号进行DFT变换并找出K个最大峰值(K为信源数)的位置,获得波达方向初始估计。
最后,构造相位旋转矩阵,对经过相位旋转的虚拟接收信号进行DFT变换,在小范围内对旋转相位η进行搜索,分别找出使得K个峰值收缩至最大的旋转相位,得到波达方向的精确估计。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
①充分利用了阵列结构的互质特性,且避免了空间平滑过程,能够获得更高的自由度;
②具有计算复杂度较低的特点;
③该方法能够获得很好的角度估计性能,例如互质阵相较于嵌套阵,互质阵的阵元分布更为稀疏,阵元间的互耦效应较弱,因此在强互耦环境下具有更好的估计性能。
附图说明
图1是本发明所使用的互质线阵结构示意图;
图2是互质线阵虚拟阵列结构示意图;
图3是当10个信号入射到互质线阵时,利用DFT算法得到的角度估计散点图;
图4是DFT算法初始估计和精确估计的性能对比;
图5是不同快拍数下DFT算法性能随信噪比变化趋势的比较;
图6是不同信源数下DFT算法性能比较。
具体实施方式
下面结合附图以及具体实施例对本发明的技术方案做进一步的详细说明:
本发明提出一种使用互质线阵进行角度估计的方法,该方法使用的阵列天线结构由两个阵元数分别为2M和N的均匀线阵组成,阵元间距分别为Nλ/2和Mλ/2,其中M和N为互质数且M<N,λ为载波波长,两个子阵只在原点处有一个阵元重合。要求阵元数是2M和N是因为这样产生的连续虚拟阵元的个数和范围有闭式解,同时可以获取更多的连续阵元。
一、数据模型
如图1所示的是一个可使用本发明的互质线阵例子,其中,M=3,N=4。
假设K个来自θk(k=1,2,…,K)的窄带信号入射到如图1所示的互质线阵上,则阵列接收信号可表示为
X=AS+N
其中,S=[s1,s2,…,sK]T为信号矩阵,sk=[sk(1),sk(2),…,sk(L)],L为快拍数,sk(l)为对第k个信号的第l次采样结果,l=1,2,…,L;N为阵列的加性高斯白噪声,均值为零、方差为A=[a(θ1),a(θ2),…,a(θK)]为阵列的方向矩阵,a(θk)为θk方向上的方向向量。两个子阵的方向向量可分别表示为
其中,d=λ/2。故使用整个阵列进行DOA估计时,方向向量可表示为
其中,a21(θk)是a2(θk)去除第一行后的向量。
二、角度估计方法
1、虚拟化
首先,求互质线阵接收信号X的协方差矩阵。实际工程中,由于信号采样都是在有限快拍数下进行,接收信号的协方差矩阵由计算得到。将进行向量化处理得到
其中,可看作一个长虚拟阵列的方向矩阵,为单快拍信号向量,和分别为噪声功率和第k个信号的信号功率,I为单位矩阵,vec(·)表示向量化操作,[·]H表示矩阵或向量的共轭转置,[·]T表示矩阵或者向量转置,[·]*代表矩阵或向量的共轭操作。
由于互质阵列的虚拟阵列由一段连续的均匀线阵和一些不连续的阵元组成,可以证明,均匀线阵的范围为[-(MN+M-1)d,(MN+M-1)d],即虚拟阵列的中间2M(N+1)-1个阵元是连续分布的。如图2所示是M=3,N=4时的虚拟阵列。由于DFT算法通常要求阵列为均匀线阵,所以我们删除矩阵中重复的行,并截取连续阵元部分,得到由连续虚拟阵列接收到的信号
其中,为连续虚拟阵列的方向矩阵,T=2M(N+1)-1为连续虚拟阵元总数,矩阵B的第(l,k)个元素为向量仅有第(T+1)/2元素为1,其余元素均为零。接下来我们对连续虚拟阵列的接收信号z使用DFT算法。
2、初始估计
首先构造归一化的DFT矩阵即
其中,矩阵F的第(p,q)个元素为对b(θk)进行DFT变换,得到 的第q个元素为
从该式可以看出,若qk=Tsinθk/2为整数,则有且仅有第qk个元素不为零。此时,θk的值可由非零元素的位置估计得出。当qk=Tsinθk/2不为整数时,的第round{Tsinθk/2}个元素及其邻近的元素均为非零值(round{·}为四舍五入操作),但由于仍然是稀疏的向量(非零元素集中在round{T sinθk/2}附近),所以峰值的位置依旧可用于对θk进行初始估计。
在实际应用中,虽然方位角未知无法计算方向向量,但我们可对虚拟化后的虚拟接收信号z进行DFT变换。令变换后的信号为yini=Fz,yini的第q个元素为
记它K个最大峰值的位置为由此可得信号波达方向的初始估计
3、精确估计
由上述分析可知,当Tsinθk/2不为整数时,角度估计的精确性无法提升。为了进一步提升该算法的估计精度,我们引入相位旋转来对误差进行补偿。
定义相位旋转矩阵Φ(η)为
其中,偏移相位η∈(-π/T,π/T)。对经过相位旋转的方向向量进行DFT变换,得到它的第q个元素为
显然,一定存在一个偏移相位ηk∈(-π/T,π/T),使得等式成立,此时有且只有一个非零元素,解决了初始估计中非零元素不唯一导致估计精度降低的问题。ηk的值可通过在(-π/T,π/T)内小范围搜索得出,即
其中,为矩阵F的第行。由此我们可以获得信号波达方向的精确估计
三、性能分析
1、空间自由度(Degree of freedom,DOF)
由上述分析可知,DFT算法可利用互质阵列的全部连续虚拟阵元,得到的空间自由度为DOFdft=2MN+2M-1。相比之下,常用的空间平滑类算法需要将连续虚拟阵元划分为若干子阵,因此空间自由度仅为DOFss=MN+M,几乎是DFT算法的一半。由此可见,在相同天线数情况下DFT算法能够获得更高的空间自由度。
2、计算复杂度
以复乘次数作为复杂度评判标准,则DFT算法的计算复杂度主要包括:计算接收信号协方差矩阵需要O{(2M+N-1)2L},对虚拟接收信号进行DFT变换需要O{T2},精确估计过程需要O{GTK},其中G表示精确估计过程中的相位搜索次数。DFT算法的总复杂度为O{(2M+N-1)2L+T2+GTK}。相比之下,DFT算法是一种低复杂度的算法。
图3当10个信号入射到互质阵列时,利用DFT算法得到的角度估计散点图。此时互质线阵的阵元数为M=7,N=9,L=500,SNR=0dB。从图中可以看出,该算法能够有效估计出波达方向。
图4是DFT算法初始估计和精确估计的性能比较,此时互质线阵的阵元数为M=7,N=9,信号的方位角为(10°,20°),L=500。其中case1是G=10的情况,case2是G=100的情况。
图5是不同快拍数下算法性能比较,此时互质线阵的阵元数为M=7,N=9,信号的方位角为(10°,20°)。
图6是不同信源数下DFT算法性能比较,此时互质线阵的阵元数为M=7,N=9,L=500。当信源数K=2时,方位角为(10°,20°);K=3时,方位角为(10°,20°,30°);K=4时,(10°,20°,30°,40°)。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (3)
1.一种使用互质线阵进行低复杂度角度估计的方法,其特征在于,包括如下步骤:
步骤1:使用互质线阵结构的阵列天线接收信号,得到阵列接收信号矩阵;
步骤2:计算所述接收信号的协方差矩阵;
步骤3:将所述协方差矩阵向量化处理,将获得的向量按相位进行排序并删除冗余行,得到虚拟阵列接收信号;
步骤4:截取所述虚拟阵列接收信号的连续阵元部分,得到连续虚拟阵列接收信号z;
步骤5:对所述连续虚拟阵列接收信号z进行DFT运算,得到接收信号谱yini,接收信号谱yini的K个最大峰值的位置为K为信号源数目;
步骤6:对所述连续虚拟阵列接收信号z进行相位旋转,之后进行DFT变换,得到旋转接收信号谱zro,在区域(-π/T,π/T)内搜索旋转接收信号谱zro的K个峰值以获得最佳相移ηk;T为连续虚拟阵元总数;
步骤7:利用最佳相移ηk和计算精确的DOA估计。
2.根据权利要求1所述的一种使用互质线阵进行低复杂度角度估计的方法,其特征在于,步骤1中,所述互质线阵由两个阵元数分别为2M和N的均匀线阵组成,阵元间距分别为Nλ/2和Mλ/2,其中,M和N为互质数且M<N,λ为载波波长;两个均匀线阵仅在原点处有一个阵元重合,阵元总数为2M+N-1。
3.根据权利要求2所述的一种使用互质线阵进行低复杂度角度估计的方法,其特征在于,步骤4中,将所述虚拟阵列接收信号按相位进行排序并删除冗余行之后截取此时向量的中间2MN+2M-1行元素作为连续虚拟阵元接收信号z。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810840884.XA CN108931758A (zh) | 2018-07-27 | 2018-07-27 | 一种使用互质线阵进行低复杂度角度估计的方法 |
CN201810840884X | 2018-07-27 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109597020A true CN109597020A (zh) | 2019-04-09 |
Family
ID=64445101
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810840884.XA Withdrawn CN108931758A (zh) | 2018-07-27 | 2018-07-27 | 一种使用互质线阵进行低复杂度角度估计的方法 |
CN201910083523.XA Pending CN109597020A (zh) | 2018-07-27 | 2019-01-29 | 一种使用互质线阵进行低复杂度角度估计的方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810840884.XA Withdrawn CN108931758A (zh) | 2018-07-27 | 2018-07-27 | 一种使用互质线阵进行低复杂度角度估计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN108931758A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927658A (zh) * | 2019-12-04 | 2020-03-27 | 南京理工大学实验小学 | 一种互质线阵中互质数的优选方法 |
CN111965591A (zh) * | 2020-07-21 | 2020-11-20 | 南京航空航天大学 | 一种基于四阶累积量矢量化dft的测向估计方法 |
CN112164898A (zh) * | 2020-09-09 | 2021-01-01 | 南京航空航天大学 | 一种双洞互质阵列天线结构 |
CN112698264A (zh) * | 2020-12-10 | 2021-04-23 | 南京航空航天大学 | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 |
CN112731280A (zh) * | 2020-12-24 | 2021-04-30 | 南京航空航天大学 | 互质阵列混合噪声环境下的esprit-doa估计方法 |
CN113296049A (zh) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN114167354A (zh) * | 2021-11-25 | 2022-03-11 | 厦门大学 | 一种基于高自由度的水下doa估计方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109932680B (zh) * | 2019-04-04 | 2022-05-20 | 哈尔滨工程大学 | 一种基于平移互质阵列的非圆信号波达方向估计方法 |
CN110297214A (zh) * | 2019-07-17 | 2019-10-01 | 南京航空航天大学 | 多互质阵列协同室内辐射源定位装置和方法 |
CN110297213A (zh) * | 2019-07-17 | 2019-10-01 | 南京航空航天大学 | 基于装载互质线阵的无人机平台的辐射源定位装置及方法 |
CN111224246B (zh) * | 2020-03-02 | 2021-03-30 | 山东大学 | 一种实现均匀线阵宽带波束大角度扫描的方法 |
CN116224215A (zh) * | 2023-01-06 | 2023-06-06 | 南京航空航天大学 | 增广互质雷达中基于泰勒展开dft算法的doa估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388462A (zh) * | 2015-11-20 | 2016-03-09 | 成都理工大学 | 一种基于互质阵列的虚拟波束形成方法 |
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
CN106291473A (zh) * | 2016-08-19 | 2017-01-04 | 西安电子科技大学 | 嵌套式天线阵列设置方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN107092004A (zh) * | 2017-05-05 | 2017-08-25 | 浙江大学 | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 |
CN107505602A (zh) * | 2017-07-25 | 2017-12-22 | 南京航空航天大学 | 嵌套阵下基于dft的doa估计方法 |
-
2018
- 2018-07-27 CN CN201810840884.XA patent/CN108931758A/zh not_active Withdrawn
-
2019
- 2019-01-29 CN CN201910083523.XA patent/CN109597020A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388462A (zh) * | 2015-11-20 | 2016-03-09 | 成都理工大学 | 一种基于互质阵列的虚拟波束形成方法 |
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
CN106291473A (zh) * | 2016-08-19 | 2017-01-04 | 西安电子科技大学 | 嵌套式天线阵列设置方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN107092004A (zh) * | 2017-05-05 | 2017-08-25 | 浙江大学 | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 |
CN107505602A (zh) * | 2017-07-25 | 2017-12-22 | 南京航空航天大学 | 嵌套阵下基于dft的doa估计方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110927658A (zh) * | 2019-12-04 | 2020-03-27 | 南京理工大学实验小学 | 一种互质线阵中互质数的优选方法 |
CN111965591A (zh) * | 2020-07-21 | 2020-11-20 | 南京航空航天大学 | 一种基于四阶累积量矢量化dft的测向估计方法 |
CN111965591B (zh) * | 2020-07-21 | 2023-04-07 | 南京航空航天大学 | 一种基于四阶累积量矢量化dft的测向估计方法 |
CN112164898A (zh) * | 2020-09-09 | 2021-01-01 | 南京航空航天大学 | 一种双洞互质阵列天线结构 |
CN112164898B (zh) * | 2020-09-09 | 2021-07-20 | 南京航空航天大学 | 一种双洞互质阵列天线结构 |
CN112698264A (zh) * | 2020-12-10 | 2021-04-23 | 南京航空航天大学 | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 |
CN112698264B (zh) * | 2020-12-10 | 2023-12-05 | 南京航空航天大学 | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 |
CN112731280A (zh) * | 2020-12-24 | 2021-04-30 | 南京航空航天大学 | 互质阵列混合噪声环境下的esprit-doa估计方法 |
CN112731280B (zh) * | 2020-12-24 | 2023-11-07 | 南京航空航天大学 | 互质阵列混合噪声环境下的esprit-doa估计方法 |
CN113296049A (zh) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN113296049B (zh) * | 2021-04-28 | 2024-02-20 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN114167354A (zh) * | 2021-11-25 | 2022-03-11 | 厦门大学 | 一种基于高自由度的水下doa估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108931758A (zh) | 2018-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109597020A (zh) | 一种使用互质线阵进行低复杂度角度估计的方法 | |
Chen et al. | Introduction to Direction-of-arrival Estimation | |
Gao et al. | Decoupled frequency diverse array range–angle‐dependent beampattern synthesis using non‐linearly increasing frequency offsets | |
CN105912791B (zh) | 虚拟互质阵列中基于局部搜索的doa估计方法 | |
CN105510874B (zh) | 一种近场信源多参数联合估计的降维music方法 | |
CN106483493B (zh) | 一种稀疏双平行线阵及二维波达方向估计方法 | |
Zhang et al. | Direction finding in MIMO radar with unknown mutual coupling | |
CN108896954A (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN103017728B (zh) | 干扰环境下天线阵列方向向量的测定方法 | |
CN109633558A (zh) | 一种基于极化时频分布的波达方向估计算法 | |
CN109633522A (zh) | 基于改进的music算法的波达方向估计方法 | |
CN109375154A (zh) | 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法 | |
CN107907853A (zh) | 一种基于均匀圆阵差分相位的单分布源doa估计方法 | |
Chen et al. | A ULA-based MWC discrete compressed sampling structure for carrier frequency and AOA estimation | |
CN101252382B (zh) | 一种宽频段信号极化与doa估计方法及装置 | |
CN109254272A (zh) | 一种共点式极化mimo雷达的两维角度估计方法 | |
Zheng et al. | Polarisation smoothing for coherent source direction finding with multiple‐input and multiple‐output electromagnetic vector sensor array | |
CN106443571A (zh) | 一种保角变换旋转不变子空间波达方向估计方法 | |
Park et al. | The fast correlative interferometer direction finder using I/Q demodulator | |
Xu et al. | Low‐angle tracking algorithm using polarisation sensitive array for very‐high frequency radar | |
CN108089146A (zh) | 一种对预估角误差鲁棒的高分辨宽带波达方向估计方法 | |
Wu et al. | Joint doppler and channel estimation with nested arrays for millimeter wave communications | |
Naceur | Improved polynomial rooting of Capon's algorithm to estimate the direction-of-arrival in smart array antenna | |
Biswas et al. | New high resolution direction of arrival estimation using Compressive Sensing | |
Tian et al. | Localization of mixed coherently and incoherently distributed sources based on generalized array manifold |
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: 20190409 |