本发明属于无线通信和雷达信号处理技术领域,具体涉及一种l型稀疏阵doa估计的精确配对方法及装置。
背景技术:
波达方向(doa)估计是天线阵列通过接收信号以特定的方法辨别空间信源的入射方向。该技术最早应用于军事领域,主要针对敌方目标进行定位并实施检测和精准打击。近年来,在雷达、声呐、导航、地震、生物医学、射电天文学等诸多领域有着广泛应用。
其中,二维doa估计在这些领域当中又扮演者至关重要的角色,因为在实际应用中仅有一维doa信息是远远不够的,例如:移动通信等数据传输的过程中往往需要知道入射信号方位角和俯仰角的二维doa信息。现有的二维doa估计方法大多是基于阵元间距等于半波长的简化面阵,如l形阵列、双平行线阵、十字形阵列等。其中,l型阵列由于结构简单、更低的克拉美罗界以及更好的估计性能,得到了广泛关注和应用。
l型阵列列通常被划分为两个线性子阵,每个子阵分别位于一个轴上,并通过两个子阵各自的一维doa估计得到方位角和俯仰角的估计值。这种方式避免了二维谱峰搜索,但是却造成了方位角和俯仰角不能正确配对的问题。针对l型阵列,已经有很多专家学者研究了关于二维doa估计的配对方法。
另一方面,针对传统均匀线阵自由度较低、估计精度和分辨率较低的缺点,人们又进行了稀疏阵列的相关研究,所研究的稀疏阵列主要指的是[p.palandp.p.vaidyanathan,nestedarrays:anovelapproachtoarrayprocessingwithenhanceddegreesoffreedom,ieeetrans.signalprocess.,58(8)(2010),4167–4181.]中的嵌套阵和[p.pal,p.p.vaidyanathan.coprimesamplingandthemusicalgorithm.ieeedsp/speworkshop2011,sedona,2011,289–294]中的互质阵,c.niu等人还将稀疏阵与l型阵列相结合提出了基于l型嵌套阵的idp算法[c.niu,y.zhang,j.guo.interlaceddouble-precision2-dangleestimationalgorithmusingl-shapednestedarrays.ieeesignalprocessingletters.23(4)(2016)522-526.]。由于l型稀疏阵与传统l型阵列的阵列结构具有较大差异,因此原先l型阵列的相关配对算法对于l型稀疏阵不再适用。但是,目前l型稀疏阵的相关研究甚少,基于l型稀疏阵配对方法的研究更少,亟需研究基于l型稀疏阵的配对方法来充分发挥稀疏阵列自由度大、分辨率高等优势,因此,该项研究具有很好的发展前景与实用价值,也是当前研究的热点与难点。
技术实现要素:
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种l型稀疏阵doa估计的精确配对方法及装置,该发明实现了二维角度的精确配对、降低了doa估计的系统成本并显著提高了估计性能、估计精度和分辨率。
为实现上述目的及其他相关目的,本发明提供一种l型稀疏阵doa估计的精确配对方法,所述l型稀疏阵包括具有n1个阵元的第一子阵和具有n2个阵元的第二子阵,所述第一子阵位于z轴上,所述第二子阵位于x轴上,所述l型稀疏阵接收k个不相关的远场窄带信号,信号入射方向和z轴、x轴的夹角分别为φ和β,入射信号在xoy平面内的投影与x轴的夹角为θ;第一子阵和第二子阵的接收数据矢量分别为x1(t)和x2(t);该配对方法包括:
分别求出第一子阵、第二子阵物理阵列的自相关矩阵和互相关矩阵;
计算第一子阵虚拟优化阵的自相关矩阵和第二子阵虚拟优化阵的自相关矩阵;
根据所述第一子阵虚拟优化阵的自相关矩阵计算俯仰角的估计值;
根据所述第二子阵虚拟优化阵的自相关矩阵和所述俯仰角的估计值求得方位角的模糊估计值;
估计信号的功率矢量,并利用该矢量解模糊,实现二维doa估计以及方位角和俯仰角的精确配对。
可选地,所述分别求出第一子阵、第二子阵物理阵列的自相关矩阵和互相关矩阵,包括:
计算第一子阵接收数据x1(t)在n次快拍下的自相关矩阵
计算第二子阵接收数据x2(t)在n次快拍下的自相关矩阵
计算第一子阵和第二子阵接收数据x1(t)和x2(t)的在n次快拍下的互相关矩阵
可选地,所述计算第一子阵虚拟优化阵的自相关矩阵和第二子阵虚拟优化阵的自相关矩阵包括:
将自相关矩阵
对矢量z1去冗余得到观测矢量
记录观测矢量
基于观测矢量
将自相关矩阵
对矢量z2去冗余得到观测矢量
基于观测矢量
可选地,所述根据所述第一子阵虚拟优化阵的自相关矩阵计算俯仰角的估计值;
将第一子阵虚拟优化阵的自相关矩阵
利用esprit算法求得俯仰角的估计值
可选地,根据所述第二子阵虚拟优化阵的自相关矩阵和所述俯仰角的估计值对应的方位角的模糊估计值;
根据第二子阵虚拟优化阵的自相关矩阵
可选地,所述估计信号的功率矢量,并利用该矢量解模糊,实现二维doa估计以及方位角和俯仰角的精确配对,包括:
估计第一子阵的方向矩阵
利用俯仰角的估计值φk估计第一子阵的方向矩阵
估计第一子阵虚拟阵列的去冗余方向矩阵
根据第一子阵的方向矩阵
根据所述第一子阵虚拟阵列的去冗余方向矩阵
基于所述俯仰角的估计值
根据所述俯仰角的估计值
为实现上述目的及其他相关目的,本发明还提供一种l型稀疏阵doa估计的精确配对装置,所述l型稀疏阵包括具有n1个阵元的第一子阵和具有n2个阵元的第二子阵,所述第一子阵位于z轴上,所述第二子阵位于x轴上,所述l型稀疏阵接收k个不相关的远场窄带信号,信号入射方向和z轴、x轴的夹角分别为φ和β,入射信号在xoy平面内的投影与x轴的夹角为θ;第一子阵和第二子阵的接收数据矢量分别为x1(t)和x2(t);该配对装置包括:
自相关与互相关模块,用于分别求出第一子阵、第二子阵物理阵列的自相关矩阵和互相关矩阵;
自相关模块,用于计算第一子阵虚拟优化阵的自相关矩阵和第二子阵虚拟优化阵的自相关矩阵;
第一估计模块,用于根据所述第一子阵虚拟优化阵的自相关矩阵计算俯仰角的估计值;
第二估计模块,根据所述第二子阵虚拟优化阵的自相关矩阵和所述俯仰角的估计值求得方位角的模糊估计值;
配对模块,用于估计信号的功率矢量,并利用该矢量解模糊,实现二维doa估计以及方位角和俯仰角的精确配对。
可选地,所述自相关与互相关模块,包括:
第一自相关单元,用于计算第一子阵接收数据x1(t)在n次快拍下的自相关矩阵
第二自相关单元,用于计算第二子阵接收数据x2(t)在n次快拍下的自相关矩阵
互相关单元,用于计算第一子阵和第二子阵接收数据x1(t)和x2(t)的在n次快拍下的互相关矩阵
可选地,所述自相关模块包括:
第一向量化模块,用于将自相关矩阵
第一去冗余模块,用于对矢量z1去冗余得到观测矢量
第一记录模块,用于记录观测矢量
第一构建模块,用于基于观测矢量
第二向量化模块,用于将自相关矩阵
第二去冗余模块,用于对矢量z2去冗余得到观测矢量
第二构建模块,用于基于观测矢量
可选地,所述配对模块包括:
第一估计单元,用于估计第一子阵的方向矩阵
第二估计单元,用于利用俯仰角的估计值φk估计第一子阵的方向矩阵
第三估计单元,用于估计第一子阵虚拟阵列的去冗余方向矩阵
第一方向矩计算单元,用于根据第一子阵的方向矩阵
功率矢量估计单元,用于根据所述第一子阵虚拟阵列的去冗余方向矩阵
模糊估计单元,基于所述俯仰角的估计值
配对单元,用于根据所述俯仰角的估计值
如上所述,本发明的一种l型稀疏阵doa估计的精确配对方法及装置,具有以下有益效果:
本发明使用稀疏阵的虚拟阵元来进行估计,相较于传统的基于l型均匀线阵(ula)的估计算法,自由度得到了极大提升,并显著提高了doa的估计性能;不仅可以估计二维doa信息,还可以估计信号的功率;利用信号的功率矢量实现了l型稀疏阵doa估计的精确配对,配对成功率高。
附图说明
为了进一步阐述本发明所描述的内容,下面结合附图对本发明的具体实施方式作进一步详细的说明。应当理解,这些附图仅作为典型示例,而不应看作是对本发明的范围的限定。
图1为本发明一实施例一种l型稀疏阵doa估计的精确配对方法中阵列设置示意图;
图2为本发明一实施例中一种l型稀疏阵doa估计的精确配对方法中方位角和俯仰角的求根均方误差随snr变化关系示意图;
图3为本发明一实施例中一种l型稀疏阵doa估计的精确配对方法中方位角和俯仰角的求根均方误差随快拍数变化关系示意图;
图4为本发明一实施例中一种l型稀疏阵doa估计的精确配对方法中方位角和俯仰角的检测概率随snr变化关系示意图;
图5为本发明一实施例中一种l型稀疏阵doa估计的精确配对方法中方位角和俯仰角的检测概率随快拍数变化关系示意图;
图6为本发明一实施例中一种l型稀疏阵doa估计的精确配对方法的流程图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
如图6所示,一种l型稀疏阵doa估计的精确配对方法,该方法包括以下步骤:
步骤1:设置天线阵列。
采用l型稀疏阵,该阵列由两个完全相同的稀疏非均匀子阵组成,子阵可以为嵌套阵或者互质阵,若子阵为嵌套阵则l型稀疏阵即是l型嵌套阵,同理,若子阵为互质阵则l型稀疏阵即是l型互质阵。为方便描述,本文采用l型嵌套阵来进行讨论。不妨设第一子阵位于z轴上,则第二子阵位于x轴上。显然,两个稀疏子阵相互垂直。同理,第一子阵也可以位于x轴上,第二子阵则位于z轴上。如图1所示,每个子阵都有n=n1+n2个阵元,n1为每个子阵密集ula的阵元数,n2为每个子阵稀疏ula的阵元数。λ表示信号波长,并且第一子阵和第二子阵共用原点处的阵元。
不妨用dz,i表示第一子阵第i个阵元的阵元位置,同样的,用dx,i表示第二子阵第i个阵元的阵元位置,其中i=1,2,...,n。
假设有k个非相关远场窄带信号sk(t)从方向(θk,φk)入射到阵列,其中,k=1,2,…,k,θk和φk分别表示第k个信号的方位角和俯仰角。噪声为独立同分布的加性高斯白噪声,且与信号独立。则l型嵌套阵中两个子阵的接收信号矢量可分别表示为:
其中,a1=[a1(φ1),a1(φ2),…,a1(φk)]表示第一子阵的阵列流型矩阵,a2=[a2(θ1,φ1),a2(θ2,φ2),…,a2(θk,φk)]表示第二子阵的阵列流型矩阵,
步骤2:分别求出第一子阵、第二子阵的自相关矩阵和互相关矩阵:
利用x1(t)求解第一子阵的自相关矩阵为:
其中,
其中,t为快拍数。同理,可以估计第二子阵的自相关矩阵为:
第一子阵和第二子阵的互相关矩阵为:
步骤3:求得第一子阵虚拟优化阵的自相关矩阵以及位置矢量;
向量化矩阵
其中,
其中,
接下来,基于矢量
则构建的
步骤4:根据所述第一子阵虚拟优化阵的自相关矩阵计算俯仰角的估计值,具体可以采用一维doa估计算法估计俯仰角:
将步骤3中第一子阵虚拟优化阵对应的自相关矩阵
其中,λs是k×k维对角矩阵,包含
然后,基于信号子空间us采用一维doa估计方法估计俯仰角,这里不妨采用esprit算法来进行估计。
将信号子空间us的前γ-1行记为us1,将us的后γ-1行记为us2,得到第一子阵的相位矩阵
其中,angle(·)为取相位算子,d为阵列的单位间距,λ为入射信号的波长。
步骤5:根据所述第二子阵虚拟优化阵的自相关矩阵和所述俯仰角的估计值对应的方位角的模糊估计值;
与步骤3类似,基于步骤2得到的自相关矩阵
然后,同步骤4,将第二子阵虚拟优化阵对应的自相关矩阵
接着,采用一维doa估计方法得到方位角所有的模糊估计值,不妨采用esprit算法来进行估计。
将信号子空间u′s的前γ-1行记为u′s1,将u′s的后γ-1行记为u′s2,得到第一子阵的相位矩阵
步骤6:估计信号的功率矢量,并利用该矢量解模糊,实现二维doa估计以及方位角和俯仰角的精确配对:
(6a)估计第一子阵的方向矩阵:
利用步骤4中俯仰角的估计值
其中,
(6b)估计第一子阵虚拟阵列的去冗余方向矩阵:
结合步骤(6a)中第一子阵的方向矩阵
其中,⊙表示求khatri-rao积。与矢量z1类似,矢量z1有冗余的元素,同样的,方向矩阵
其中,
(6c)估计信号功率矢量:
由于
其中,功率矢量
(6d)求得第二子阵方向矩阵的模糊估计值:
基于步骤4中俯仰角的估计值
其中,
(6e)通过最大化相关度量解模糊,实现二维角度的精确配对:
给定俯仰角估计值
其中,⊙表示求khtri-rao积,||·||f表示求矩阵的f范数。
则使得相关度量ηi最大的i即是最优的i值,与俯仰角
这样就完成了基于l型稀疏阵的二维doa估计,并实现了方位角
为了分析本发明所提算法与文献[j.-f.gu,p.wei.jointsvdoftwocross-correlationmatricestoachieveautomaticpairingin2-dangleestimationproblems.ieeeantennasandwirelesspropagationletters.6(2007)553-556.]的jsvd算法以及文献[c.niu,y.zhang,j.guo.interlaceddouble-precision2-dangleestimationalgorithmusingl-shapednestedarrays.ieeesignalprocessingletters.23(4)(2016)522-526.]的idp算法的估计性能,设计了两组仿真实验来进行比较。其中,本发明与文献[c.niu,y.zhang,j.guo.interlaceddouble-precision2-dangleestimationalgorithmusingl-shapednestedarrays.ieeesignalprocessingletters.23(4)(2016)522-526.]均采用l型嵌套阵,阵列参数为n1=n2=3,文献[j.-f.gu,p.wei.jointsvdoftwocross-correlationmatricestoachieveautomaticpairingin2-dangleestimationproblems.ieeeantennasandwirelesspropagationletters.6(2007)553-556.]采用的l型阵列,阵列参数为m=6,显然l型嵌套阵和l型阵列的阵元总数相同,均有11个阵元(两个轴的子阵共用零位置阵元)。信号数为2,入射方向分别为(θ1,φ1)=(55°,70°)和(θ2,φ2)=(60°,60°)。定义检测概率为:如果方位角和俯仰角的估计值与真实值的偏差均不超过标准值(本实验设为1.5°),则检测成功。否则检测失败。
第一组试验快拍数为5000,并进行1000次独立试验,方位角和俯仰角估计的求根均方误差(rmse)随信噪比(snr)变化的关系如图2所示,方位角和俯仰角的检测概率随snr变化关如图4所示。
另一组试验信噪比为5db,同样进行1000次独立试验,方位角和俯仰角的求根均方误差(rmse)随快拍数变化的关系如图3所示,方位角和俯仰角的检测概率随快拍数的变化关如图5所示。
从图中可以看出,本发明所提的基于l型稀疏阵doa估计的精确配对方法能够很好的提高二维doa估计性能,降低系统成本,并且无需谱搜索,计算复杂度较低。同时,方位角和俯仰角的检测概率高、配对成功率高。
本实施例还公开一种l型稀疏阵doa估计的精确配对装置,所述l型稀疏阵包括具有n1个阵元的第一子阵和具有n2个阵元的第二子阵,所述第一子阵位于z轴上,所述第二子阵位于x轴上,所述l型稀疏阵接收k个不相关的远场窄带信号,信号入射方向和z轴、x轴的夹角分别为φ和β,入射信号在xoy平面内的投影与x轴的夹角为θ;第一子阵和第二子阵的接收数据矢量分别为x1(t)和x2(t);该配对装置包括:
自相关与互相关模块,用于分别求出第一子阵、第二子阵物理阵列的自相关矩阵和互相关矩阵;
自相关模块,用于计算第一子阵虚拟优化阵的自相关矩阵和第二子阵虚拟优化阵的自相关矩阵;
第一估计模块,用于根据所述第一子阵虚拟优化阵的自相关矩阵计算俯仰角的估计值;
第二估计模块,根据所述第二子阵虚拟优化阵的自相关矩阵和所述俯仰角的估计值求得方位角的模糊估计值;
配对模块,用于估计信号的功率矢量,并利用该矢量解模糊,实现二维doa估计以及方位角和俯仰角的精确配对。
在一些实施例中,所述自相关与互相关模块,包括:
第一自相关单元,用于计算第一子阵接收数据x1(t)在n次快拍下的自相关矩阵
第二自相关单元,用于计算第二子阵接收数据x2(t)在n次快拍下的自相关矩阵
互相关单元,用于计算第一子阵和第二子阵接收数据x1(t)和x2(t)的在n次快拍下的互相关矩阵
在一些实施例中,所述自相关模块包括:
第一向量化模块,用于将自相关矩阵
第一去冗余模块,用于对矢量z1去冗余得到观测矢量
第一记录模块,用于记录观测矢量
第一构建模块,用于基于观测矢量
第二向量化模块,用于将自相关矩阵
第二去冗余模块,用于对矢量z2去冗余得到观测矢量
第二构建模块,用于基于观测矢量
在一些实施例中,所述配对模块包括:
第一估计单元,用于估计第一子阵的方向矩阵
第二估计单元,用于利用俯仰角的估计值φk估计第一子阵的方向矩阵
第三估计单元,用于估计第一子阵虚拟阵列的去冗余方向矩阵
第一方向矩计算单元,用于根据第一子阵的方向矩阵
功率矢量估计单元,用于根据所述第一子阵虚拟阵列的去冗余方向矩阵
模糊估计单元,基于所述俯仰角的估计值
配对单元,用于根据所述俯仰角的估计值
需要说明的是,由于装置部分的实施例与方法部分的实施例相互对应,因此装置部分的实施例的内容请参见方法部分的实施例的描述,这里暂不赘述。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本发明所提供的实施例中,应该理解到,所揭露的装置/终端设备和方法,可以通过其它的方式实现。例如,以上所描述的装置/终端设备实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。