一种使用互质线阵进行低复杂度角度估计的方法技术

本发明专利技术公开了一种使用互质线阵进行低复杂度角度估计的方法，首先，对阵列接收信号的协方差矩阵进行向量化处理，删除冗余行并截取连续阵元部分，得到阵元间距为λ/2的虚拟均匀线阵的接收信号；然后，对虚拟接收信号进行离散傅里叶变换(Discrete Fourier Transform,DFT)并找出K个最大峰值(K为信源数)的位置，由此可得波达方向的初始估计；最后，对经过相位旋转的虚拟接收信号进行DFT变换，在小范围内对旋转相位η进行搜索，分别找出使得K个峰值收缩至最大的旋转相位，并由此得到波达方向的精确估计。本发明专利技术将DFT变换与互质阵列的DOA估计问题相结合，不需要对虚拟接收信号进行解相关处理，能够获得相对较高的空间自由度，同时具有较低的计算复杂度和较高的估计精度。

【技术实现步骤摘要】
一种使用互质线阵进行低复杂度角度估计的方法
本专利技术涉及空间谱估计领域，特别是涉及一种使用互质线阵进行低复杂度角度估计的方法。
技术介绍
当前，一种被称为互质阵的阵列天线布局得到了关注，该阵列结构突破了传统天线阵阵元间距半波长的限制，使得天线孔径得到极大的扩展，能够获得角度估计性能的提升。使用互质线阵进行角度估计的方法有解模糊和虚拟化两类。其中，解模糊方法能够获得比均匀线阵更好的估计性能，但使用两个子阵分别进行角度估计的方法使得空间自由度大大降低。而虚拟化方法则能够获得比实际阵元数更高的空间自由度，极大地增加了可探测的信源数。由于虚拟化后等效的虚拟信号是单快拍的相干信号，基于虚拟化的角度估计算法通常需要使用空间平滑技术对接收信号进行解相关处理，这一过程降低了阵列的空间自由度。
技术实现思路
本专利技术所要解决的技术问题是提供一种使用互质线阵进行角度估计的方法，该方法能够获得较大的空间自由度，同时具有较低的计算复杂度和较高的估计精度。可以应用于无线通信、声呐、定位等领域中。本专利技术为解决上述技术问题采用以下技术方案：提出一种结合虚拟化和离散傅里叶变换(DiscreteFourierTransform,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{Tsinθ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、空间自由度(Degreeoffreedom,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+G本文档来自技高网...

【技术保护点】
1.一种使用互质线阵进行低复杂度角度估计的方法，其特征在于，包括如下步骤：步骤1：使用互质线阵结构的阵列天线接收信号，得到阵列接收信号矩阵；步骤2：计算所述接收信号的协方差矩阵；步骤3：将所述协方差矩阵向量化处理，将获得的向量按相位进行排序并删除冗余行，得到虚拟阵列接收信号；步骤4：截取所述虚拟阵列接收信号的连续阵元部分，得到连续虚拟阵列接收信号z；步骤5：对所述连续虚拟阵列接收信号z进行DFT运算，得到接收信号谱yini，接收信号谱yini的K个最大峰值的位置为

【技术特征摘要】
2018.07.27 CN 201810840884X1.一种使用互质线阵进行低复杂度角度估计的方法，其特征在于，包括如下步骤：步骤1：使用互质线阵结构的阵列天线接收信号，得到阵列接收信号矩阵；步骤2：计算所述接收信号的协方差矩阵；步骤3：将所述协方差矩阵向量化处理，将获得的向量按相位进行排序并删除冗余行，得到虚拟阵列接收信号；步骤4：截取所述虚拟阵列接收信号的连续阵元部分，得到连续虚拟阵列接收信号z；步骤5：对所述连续虚拟阵列接收信号z进行DFT运算，得到接收信号谱yini，接收信号谱yini的K个最大峰值的位置为K为信号源数目；步骤6：对所述连续虚拟阵列接收信号z进行相位旋转，之后进行DFT变换，得到...

【专利技术属性】
技术研发人员：张小飞，汪云飞，沈金清，郑旺，施展，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32