申请/专利权人：山东农业大学

申请日：2017-11-30

公开（公告）日：2020-09-15

公开（公告）号：CN108594164B

主分类号：G01S3/00(20060101)

分类号：G01S3/00(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.09.15#授权;2018.10.26#实质审查的生效;2018.09.28#公开

摘要：本发明公开了一种平面阵列DOA估计方法及设备。所述平面阵列DOA估计方法，包括：步骤1：使用平面阵列接收来自至少一个信源的信号，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；步骤2：基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；步骤3：基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。能够适用于阵列尺寸受限及实时性要求高的场合，可有效地以较少的阵元提供低复杂度、高准确度的波达方向估计。

主权项：1.一种平面阵列DOA估计方法，其特征在于，包括：步骤1：使用平面阵列接收来自至少一个信源的信号，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；步骤2：基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；步骤3：基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值；所述相互平行的互质线性子阵列对包括：第一线性子阵列和第二线性子阵列，所述第一线性子阵列包含M1个阵元，所述M1个阵元沿y轴方向以M2λ2为间隔进行布置，所述第二线性子阵列包含M2个阵元，所述M2个阵元沿y轴方向以M1λ2为间隔进行布置，所述第一线性子阵列的第一个阵元和所述第二线性子阵列的第一个阵元沿x轴方向对齐且间隔小于或等于λ2，其中，M1和M2是互质的正整数，λ为信号波长；其中，所述x轴与y轴垂直；所述步骤2包括：步骤2-1：获取所述互质线性子阵列对的实际接收信号x1t和x2t，其中， x1t和x2t分别表示第一线性子阵列和第二线性子阵列在t时刻的实际接收信号；矩阵A1＝[a1α1,…,a1αK]和A2＝[a2α1,…,a2αK]分别表示第一线性子阵列和第二线性子阵列沿y轴的流型矩阵，和分别表示第一线性子阵列和第二线性子阵列对于第k个信源的导向矢量；矩阵为对角矩阵；st＝[s1t,s2t,…,sKt]T为由入射角度分别为αk,βk,k＝1,2,…,K的K个信源信号组成的信号矢量，上标T表示转置运算，αk表示入射方向与y轴之间的夹角，βk表示入射方向与x轴之间的夹角，Ak为信号振幅，ωk为信号频率；矢量n1t和n2t分别表示第一线性子阵列和第二线性子阵列所实际接收到的均值为零、方差为的加性高斯白噪声矢量，矢量n1t和n2t与信号矢量st不相关；步骤2-2：基于第一线性子阵列的实际接收信号x1t和第二线性子阵列的实际接收信号x2t，计算和构建所述互质线性子阵列对的虚拟接收信号R1τ和R2τ： 其中， 其中，上标*表示共轭运算，Rsτ表示K个信源信号在不同时刻下的自相关矢量，A1、A2分别表示沿y轴的不同虚拟流型矩阵，R1τ和R2τ分别表示以第一线性子阵列的第一个阵元an,bn＝0,0为中心，经共轭增广处理后的虚拟接收信号矢量，R1′τ和R2′τ分别表示以第二线性子阵列的第一个阵元为中心，经共轭增广处理后的虚拟接收信号矢量，和R1-τ分别表示矩阵和R1-τ*的倒数M1-1行子矩阵，R1′-τ表示矩阵R1′-τ*的倒数M1-1行子矩阵，和R2′-τ分别表示和R2′-τ*的倒数M2-1行子矩阵，其中，表示通过第一线性子阵列和第二线性子阵列的两个阵元处的实际接收信号之间的相关运算所直接得到的虚拟接收信号，xmt和xnt分别表示由所述x轴和y轴所确定的平面坐标系中的am,bm和an,bn位置处的两个阵元的实际接收信号，所述对应中以第一线性子阵列的第一个阵元an,bn＝0,0为中心的情况，对应中以第二线性子阵列的第一个阵元为中心的情况。

全文数据：一种平面阵列DOA估计方法及设备技术领域[0001]本发明涉及通信信号处理领域，尤其涉及一种平面阵列DOA估计方法及设备。背景技术[0002]波达方向（DOA估计是阵列信号处理中的重要研究内容，在雷达、声呐等领域应用广泛。按照一定规律将多个天线排列构成天线阵列，可用来测定辐射源的来波方向，从而实现辐射源的测向。在民用领域，快速准确的测向定位是实现无线电频谱监测、非法用频设备如伪基站、黑广播等查找与定位的迫切要求。在军事领域，快速、准确、隐蔽地对目标辐射源进行测向定位，既能最大限度保护己方，又能精准打击敌方军事目标，是关乎战争结果的重要因素。[0003]传统波达方向估计方法，如多重信号分类（MUSIC、旋转不变技术估计算法ESPRIT，利用N天线均匀线性阵列时，最多可分辨N-I个信号源。为提高分辨能力，非均匀的阵列结构如互质阵列逐渐引起了研究者的重视。互质阵列由呈互质关系的两个均匀子阵构成，可检测多于天线数目的辐射源。由于DOA估计性能受到阵列孔径的限制，在尺寸受限的场合下布设大孔径阵列天线非常困难，导致估计性能不高。同时二维空间谱搜索和二维角度配对将导致计算复杂度过高，因此，难以应用至实时性要求高的场合。[0004]因此，至少需要提出新的技术方案来对现有技术方案的不足之处进行改进。发明内容[0005]本发明的目的是通过以下技术方案实现的。[0006]根据本发明的平面阵列DOA估计方法，包括：[0007]步骤1:使用平面阵列接收来自至少一个信源的信号，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；[0008]步骤2:基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；[0009]步骤3:基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。[0010]根据本发明的平面阵列DOA估计方法，所述相互平行的互质线性子阵列对包括：[0011]第一线性子阵列和第二线性子阵列，所述第一线性子阵列包含M1个阵元，所述此个阵元沿y轴方向Μ2λ2为间隔进行布置，所述第二线性子阵列包含M2个阵元，所述M2个阵元沿y轴方向以此12为间隔进行布置，所述第一线性子阵列的第一个阵元和所述第二线性子阵列的第一个阵元沿X轴方向对齐且间隔小于或等于λ2,其中，MdPM2是互质的正整数，λ为信号波长。[0012]根据本发明的平面阵列DOA估计方法，所述步骤2包括：[0013]步骤2-1:获取所述互质线性子阵列对的实际接收信号^⑴和X2t，其中，[0014][0015]X1⑴和X2⑴分别表示第一线性子阵列和第二线性子阵列在t时刻的实际接收信号，矩阵Ai=[aiCX1，…，aiακ]和A2=[a2αι，…，a2ακ]分别表不第一线性子阵列和第二线性子阵列沿y轴的流型矩阵，和分别表示第一线性子阵列和第二线性子阵列对于第k个信源的导向矢量，矩阵为对角矩阵，s⑴=[S1t，s2t，···，8κα]τ为由入射角度分别为akJk，k=l，2，"_，K，（K彡1的K个信源信号组成的信号矢量，上标T表示转置运算，ak表示入射方向与y轴之间的夹角，0k表示入射方向与X轴之间的夹角Ak为信号振幅，cok为信号频率，矢量m⑴和n2⑴分别表示第一线性子阵列和第二线性子阵列所实际接收到的均值为零、方差为K进行迭代运算，直至获得与一维字典中各栅格点一一对应的能量值矢量P=[PiP1，…，Pd]T的估计值矢量L对于第i次迭代，包括以下步骤：[0046]固定Θ，按照下式对能量值矢量p进行更新，得到p1，[0047][0048]固定P1，按照公式更新©⑷，或者，按照公式更新α1，其中，为变量α的梯度变化最大的方向，〜为步长，[0049]步骤3-2-2:将估计值矢量中的第k个非零项所对应的角度值作为第k个夹角ak的角度估计值，并且按照下式获取相应的夹角的角度估计值，[0050][0051]其中表示中的第k个非零元素。[0052]根据本发明的平面阵列DOA估计设备，所述设备包括平面阵列、处理器和存储有可执行指令的存储器，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对，所述处理器执行所述可执行指令来完成根据上文所述的方法中的步骤。[0053]根据本发明的平面阵列DOA估计设备，包括：[0054]平面阵列模块，用于接收来自至少一个信源的信号，所述平面阵列模块包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；[0055]虚拟接收信号计算和构建模块，其与平面阵列模块连接，用于基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；[0056]DOA估计模块，其与虚拟接收信号计算和构建模块连接，用于基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。[0057]本发明的优点在于:能够利用阵列接收信号的空时差异性进行共辄增广处理，能够在不增加实际阵元的条件下，构建具有更多阵元和更大观测孔径的虚拟阵列，改善估计性能。能够利用平行子阵列之间的互相关特性，将二维测向问题转为一维问题，来降低实现复杂度。能够适用于阵列尺寸受限及实时性要求高的场合，可有效地以较少的阵元提供低复杂度、高准确度的波达方向估计。附图说明[0058]通过阅读下文具体实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出具体实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：[0059]图1示出了根据本发明实施方式的平面阵列DOA估计方法的示意流程图。[0060]图2示出了根据本发明实施方式的相互平行的互质线性子阵列的结构示意图。[0061]图3示出了根据本发明实施方式的虚拟阵列对的结构示意图。[0062]图4示出了根据本发明实施方式的平面阵列DOA估计方法所得到的估计角度与真实角度之间的关系示意图。[0063]图5示出了根据本发明实施方式的平面阵列DOA估计方法的均方根误差与信噪比之间的关系不意图。[0064]图6示出了根据本发明实施方式的平面阵列DOA估计方法的均方根误差与快拍数之间的关系不意图。具体实施方式[0065]下面将参照附图更详细地描述本公开的示例性实施方式。虽然附图中显示了本公开的示例性实施方式，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。[0066]图1示出了根据本发明实施方式的平面阵列DOA估计方法100的示意流程图。[0067]如图1所示，平面阵列DOA估计方法100包括以下步骤：[0068]步骤S102:使用平面阵列接收来自至少一个信源的信号，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对。[0069]步骤S104:基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号。[0070]步骤S106:基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。[0071]图2示出了根据本发明实施方式的相互平行的互质线性子阵列200的结构示意图。[0072]如图2所示，平面阵列DOA估计方法100所使用的相互平行的互质线性子阵列对200包括第一线性子阵列（即，图2所示的子阵列1和第二线性子阵列（S卩，图2所示的子阵列2。[0073]第一线性子阵列包含M1个阵元，所述施个阵元沿y轴方向Μ2λ2为间隔进行布置，所述第二线性子阵列包含M2个阵元，所述此个阵元沿y轴方向以施12为间隔进行布置，所述第一线性子阵列的第一个阵元和所述第二线性子阵列的第一个阵元沿X轴方向对齐且间隔小于或等于λ2,其中，MjPM2是互质的正整数，λ为信号波长。[0074]尽管在图2中仅仅示出了一个相互平行的互质线性子阵列对200,然而，平面阵列DOA估计方法100所使用的平面阵列可以包括多个相互平行的互质线性子阵列对200。[0075]尽管在图1中未示出，然而，可选地，上述步骤S104可以包括以下步骤：[0076]步骤2-1:获取所述互质线性子阵列对的实际接收信号^⑴和X2t。[0077]S卩，⑴构建二维接收阵列模型。[0078]S卩，利用如图2所示的二维接收阵列模型接收来自至少一个信源的信号。[0079]其中，[0080][0081]X1t和X2t分别表示第一线性子阵列和第二线性子阵列在t时刻的实际接收信号，矩阵Ai=[aiCX1，…，aiακ]和A2=[a2αι，…，a2ακ]分别表不第一线性子阵列和第二线性子阵列沿y轴的流型矩阵和分别表示第一线性子阵列和第二线性子阵列对于第k个信源的导向矢量，矩阵为对角矩阵，s⑴=[S1t，s2t，···，8κα]τ为由入射角度分别为akJk，k=l，2，"_，K，（K彡1的K个信源信号组成的信号矢量，上标T表示转置运算，ak表示入射方向与y轴之间的夹角，0k表示入射方向与X轴之间的夹角.Ak为信号振幅，COk为信号频率，矢量m⑴和n2⑴分别表示第一线性子阵列和第二线性子阵列所实际接收到的均值为零、方差为σ,丨的加性高斯白噪声矢量，矢量m⑴和n2t与信号矢量st不相关。[0082]步骤2-2:基于第一线性子阵列的实际接收信号^t和第二线性子阵列的实际接收信号X2⑴计算和构建所述互质线性子阵列对的虚拟接收信号R1⑴和R2τ。[0083]即，（2共辄增广空时处理扩展孔径[0084]令am，bm和an，bn分别表示阵列中两个不同阵元的位置，其对应的接收信号分别表示为^⑴和xnt。则不同时间标签下的互相关函数定义为：[0086]其中，上标*表示共辄运算，和分别表示信号源和噪声在不同时刻下的自相关函数，可分别表示为：[0088]由此可知，互相关函数可进一步简化为：[0089][0090]由此可知，互相关函数中出现了不同位置阵元的位置差分，此时可形成更多的虚拟阵元。下面分别以第一线性子阵列（即，图2所示的子阵列1和第二线性子阵列（S卩，图2所示的子阵列2的首个阵元为中心进行处理。[0091]1以子阵列1的首个阵元S卩，第一线性子阵列的第一个阵元）（an，bn=0,0为中心，可得：[0092][0093]定义矩阵Rsτ:[0094][0095]表示K个信源信号在不同时刻下的自相关矢量。[0096]定义如下两个矩阵：[0099]可得：[0100][0101]根据共辄对称性，有，则成立。分别令和R1_τ表示矩阵和R1_τΓ的最后M1-I行子矩阵，则有：[0102][0103]可以看出，R1⑴和R2⑴分别表示以子阵列1的首个阵元为中心，经共辄增广处理后的虚拟接收信号矢量，心和知表示相应的阵列流型矩阵。[0104]2以子阵列2的首个阵元（S卩，第二线性子阵列的第一个阵元）为中心，可得：[0106]定义如下两个矩阵：[0109]则有：[0110][0111]由于成立，则：[0112][0113]令RlMτ为矩阵Rn_τ*的后Mi-I行子矩阵，和R2M0分别为和Rc^_τΓ的后M2-I行矩阵，则：[0114][0115]可以看出，RlMτ和R2M⑴分别以子阵列2首个阵元为中心进行共辄增广处理的虚拟接收信号矢量，和分别表示相应的阵列流型矩阵。[0116]即，（3虚拟阵列搭建。[0117]经共辄增广处理后，对矩阵R1_、RW⑴、R2^t和Rcn⑴进行整合，有：[0118][0119]对R⑵⑴和R1’-⑴进行整合，则有：[0120][0121]可以看出，R1τ和此⑴可视为虚拟阵列的等效接收信号。该虚拟阵列从子阵列间的互相关矩阵的共辄增广处理中得来，由两个阵列的接收信号在不同时刻和不同位置的差分形成。与原始阵列相比，该虚拟阵列包含更多阵元，扩展了阵列孔径，从而进一步改善了阵列的自由度和检测能力。[0122]图3示出了根据本发明实施方式的虚拟阵列对300的结构示意图。[0123]如图3所示，虚拟阵列对300从子阵列（S卩，第一线性子阵列和第二线性子阵列）间的互相关矩阵的共辄增广处理中得来，由两个子阵列的接收信号在不同时刻和不同位置的差分形成。即，R1τ和R2τ可视为虚拟阵列的等效接收信号。与图2所示的原始子阵列对200相比，该虚拟阵列对300包含更多阵元，扩展了阵列孔径，从而进一步改善了阵列的自由度和检测能力。[0124]尽管在图1中未示出，然而，可选地，上述步骤S106可以包括以下步骤：[0125]步骤3-1:对所述互质线性子阵列对的虚拟接收信号R1τ和办⑴进行相关运算得到虚拟协方差矩阵RC，对虚拟协方差矩阵Rc的矩阵表达式进行矢量化处理。[0126]可选地，所述步骤3-1包括:通过以下公式对虚拟协方差矩阵Rc的矩阵表达式进行矢量化处理。[0127]S卩，⑷降维处理。[0128]该虚拟阵列中两子阵的等效接收信号可分别表示为：[0131]其中等效信号矢量Seqv=[RsTs，Rs2Ts，···，RsNpTs]，TS为等效抽样周期，Np是等效快拍数。[0132]两虚拟子阵的互协方差矩阵为：[0133][0134]其中，Reqv=E[seqvSeqvη]为对角矩阵，上标H表示共辄转置运算，其第k个对角元素为β矩阵ΦReqv亦为对角矩阵。[0135]对互协方差矩阵Rc矢量化，可得：[[0137]其中等效阵列流型矩阵，符号θ表示Khatri-Rao乘积，矢量u中包含着对角矩阵（ORe3qv中的对角元素。由于Re3qv为实值对角矩阵而矩阵Φ中对角元素的相位信息与β相关，这意味着矢量u中元素的相位信息与真实相位β相关，从而估计出u之后，其对应的角度β也可轻易计算得出。[0138]步骤3-2:基于经过矢量化处理的虚拟协方差矩阵Rc的表达式，使用一维字典来估计所述至少一个信源的DOA的数值。[0139]可选地，所述步骤3-2包括以下步骤：[0140]即，⑸一维字典学习。[0141]步骤3-2-1:基于一维字典{θιθ^—,θο}DK进行迭代运算，直至获得与一维字典中各栅格点对应的能量值矢量P=[Pi，Pi，…，PD]τ的估计值矢量[0142]对于第i次迭代，包括以下步骤：[0M3]1固定Θ，按照下式对能量值矢量p进行更新，得到p1，[0144][0145]例如，该优化问题可利用凸优化工具包CVX求解，得到稀疏解表示为p1=CVXr，Θ㈣）。[0146]2固定P1，按照公式更新©⑷，或者，按照公式更新α1，其中，为变量α的梯度变化最大的方向，〜为步长。[0147]步骤3-2-2:将估计值矢量中的第k个非零项所对应的角度值作为第k个夹角ak的角度估计值，并且按照下式获取对应的夹角的角度估计值，[0148][0149]其中，表示中的第k个非零元素。[0150]通过上述步骤3-2-1和步骤3-2-2,在有限次迭代更新中，就可以得到稀疏信号的估计值及更新后字典矩阵^,矩阵中的非零项所对应的字典位置表示角度α的估计值，非零元素本身的相位即为邱勺估计值。[0151]S卩，通过上述步骤3-2-1和步骤3-2-2实现了，对一维角度域进行离散栅格化：{θ^θ2，···，θ〇}DK，建立如下稀疏重构优化问题：[0154]其中Θ为从栅格点{0^02，···，θ〇}所构建的字典矩阵，P=[Pl，Pl，…，pD]T为对应栅格点的能量。角度值可通过寻找P中的非零项估计，具体来说根据P中的非零项位置查找Θ即可确定角度α的估计值，P中的非零项本身所对应的相位为角度β的估计值。注意到上述优化问题只涉及一维字典矩阵学习来重建稀疏信号，因此其复杂度得到极大降低。[0155]上述优化问题可转换为无约束优化问题，[0156][0157]其中η为规则化参数用来平衡稀疏性和准确度。[0158]针对有限栅格所造成的栅格失配问题，设计了基于迭代字典学习的快速测向方法。在稀疏重构过程中，角度a4Pi3k实现了自动配对，降低了实现复杂度。[0159]为了使本领域技术人员更直观地了解上文提出的平面阵列DOA估计方法100的技术效果，下文中给出了采用平面阵列DOA估计方法100所得到的部分仿真结果。[0160]图4示出了根据本发明实施方式的平面阵列DOA估计方法所得到的估计角度与真实角度之间的关系示意图。所使用的信源个数为16,快拍数为500,信噪比为5分贝，平行互质阵列第一线性子阵列和第二线性子阵列的阵元数分别为3和4。[0161]从图4可以看出，根据本发明的平面阵列DOA估计方法100能够成功检测多于阵元个数的信源。因此，当利用接收信号的空时差异性对阵列进行共辄增广，所构建的虚拟阵列可有效扩展阵列孔径，显著提高阵列的检测能力和自由度。[0162]图5示出了根据本发明实施方式的平面阵列DOA估计方法的均方根误差与信噪比之间的关系示意图。所使用的快拍数为500,信噪比为-5至15分贝。[0163]如图5所示，根据本发明的平面阵列DOA估计方法100在相同快拍数500不同信噪比下的均方根误差明显小于现有技术的稀疏重构方法和多项式求根方法的均方根误差。[0164]图6示出了根据本发明实施方式的平面阵列DOA估计方法的均方根误差与快拍数之间的关系示意图。所使用的信噪比为30分贝，快拍数为100至1000。[0165]如图6所示，根据本发明的平面阵列DOA估计方法100在相同信噪比（30分贝）不同快拍数下的均方根误差也明显小于现有技术的稀疏重构方法和多项式求根方法的均方根误差。[0166]以上实验结果表明，根据本发明的平面阵列DOA估计方法100在不增加阵元个数及成本的条件下，可有效实现阵列孔径的扩展，提高检测能力。降维处理将二维问题转为一维问题，有效降低实现复杂度。基于字典学习的稀疏优化问题将有助于提高测向精度。因此，根据本发明的平面阵列DOA估计方法100有助于提高估计精度和降低实现复杂度，对改善尺寸受限及实时性要求高的场合下的性能有着重要的应用价值。[0167]综上所述，根据本发明的上述技术方案至少包括⑴构建二维接收阵列模型；（2共辄增广空时处理扩展孔径；（3虚拟阵列搭建；（4利用虚拟阵列两个平行子阵之间的互相关矩阵，对其进行矢量化处理，将原始二维测向问题转化为一维问题；（5—维字典学习，实现角度的自动配对等处理步骤，至少具有以下优点：[0168]1能够利用平行互质阵列接收信号的空时差异性进行共辄增广处理，以此实现孔径扩展，增强阵列自由度及检测能力。[0169]2能够利用平行子阵列之间的互相关特性，将二维测向问题转为一维问题，来降低实现复杂度。[0170]3能够在不增加实际阵元的条件下，构建具有更多阵元和更大观测孔径的虚拟阵列，改善估计性能，适用于阵列尺寸受限及实时性要求高的场合，可有效地以较少的阵元提供低复杂度、高准确度的波达方向估计。[0171]⑷能够实现角度的自动配对，降低实现复杂度。[0172]结合上文提出的平面阵列DOA估计方法100,还提出了一种平面阵列DOA估计设备，所述设备包括平面阵列、处理器和存储有可执行指令的存储器，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对，所述处理器执行所述可执行指令来完成根据上文所述的平面阵列DOA估计方法100中的步骤。[0173]结合上文提出的平面阵列DOA估计方法100,还提出了另一种平面阵列DOA估计设备，包括：[0174]平面阵列模块，用于接收来自至少一个信源的信号，所述平面阵列模块包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；[0175]虚拟接收信号计算和构建模块，其与平面阵列模块连接，用于基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；[0176]DOA估计模块，其与虚拟接收信号计算和构建模块连接，用于基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。[0177]以上所述，仅为本发明示例性的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

权利要求：1.一种平面阵列DOA估计方法，其特征在于，包括：步骤1:使用平面阵列接收来自至少一个信源的信号，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；步骤2:基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；步骤3:基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。2.根据权利要求1所述的平面阵列DOA估计方法，其特征在于，所述相互平行的互质线性子阵列对包括：第一线性子阵列和第二线性子阵列，所述第一线性子阵列包含M1个阵元，所述M1个阵元沿y轴方向以跑12为间隔进行布置，所述第二线性子阵列包含M2个阵元，所述跑个阵元沿y轴方向以此12为间隔进行布置，所述第一线性子阵列的第一个阵元和所述第二线性子阵列的第一个阵元沿X轴方向对齐且间隔小于或等于λ2,其中，MjPM2是互质的正整数，λ为信号波长。3.根据权利要求2所述的平面阵列DOA估计方法，其特征在于，所述步骤2包括：步骤2-1:获取所述互质线性子阵列对的实际接收信号^⑴和Χ2⑴，其中，X1t和X2⑴分别表示第一线性子阵列和第二线性子阵列在t时刻的实际接收信号，矩阵Ai=[aiαι，…，aiα£]和A2=[a2αι，…，a2ακ]分别表不第一线性子阵列和第二线性子阵列沿y轴的流型矩阵，和分别表示第一线性子阵列和第二线性子阵列对于第k个信源的导向矢量，矩阵为对角矩阵，S⑴=[sit，S2⑴，…，SK⑴]T为由入射角度分别为ak，，k=l，2，-_，K，（Κ彡1的K个信源信号组成的信号矢量，上标T表示转置运算，ak表示入射方向与y轴之间的夹角，知表示入射方向与X轴之间的夹角，，Ak为信号振幅，Cok为信号频率，矢量m⑴和n2⑴分别表示第一线性子阵列和第二线性子阵列所实际接收到的均值为零、方差为的加性高斯白噪声矢量，矢量m⑴和Π2⑴与信号矢量St不相关；步骤2-2:基于第一线性子阵列的实际接收信号^t和第二线性子阵列的实际接收信号X2⑴，计算和构建所述互质线性子阵列对的虚拟接收信号R1⑴和R2τ:其中，其中，上标*表示共辄运算，Rsτ表示K个信源信号在不同时刻下的自相关矢量，、、如分别表示沿y轴的不同虚拟流型矩阵，R1τ和R2τ分别表示以第一线性子阵列的第一个阵元an，bn=0,0为中心，经共辄增广处理后的虚拟接收信号矢量，R1~⑴和R2~⑴分别表示以第二线性子阵列的第一个阵元为中心，经共辄增广处理后的虚拟接收信号矢量，和RH⑴分别表示矩阵和R1_τΓ的倒数M1-I行子矩阵，R1M⑴表示矩阵Rn-τ广的倒数._1行子矩阵，和R2M⑴分别表示和Rcn_τΓ的倒数M2-I行子矩阵，其中，表示通过第一线性子阵列和第二线性子阵列的两个阵元处的实际接收信号之间的相关运算所直接得到的虚拟接收信号，Xmt*xnt分别表示由所述X轴和y轴所确定的平面坐标系中的am，bm和an，bn位置处的两个阵元的实际接收信号，所述对应中以第一线性子阵列的第一个阵元an，bn=0,0为中心的情况I对应中以第二线性子阵列的第一个阵元为中心的情况。4.根据权利要求3所述的平面阵列DOA估计方法，其特征在于，所述步骤3包括：步骤3-1:对所述互质线性子阵列对的虚拟接收信号R1τ和办τ进行相关运算得到虚拟协方差矩阵Rc，对虚拟协方差矩阵Rc的矩阵表达式进行矢量化处理；步骤3-2:基于经过矢量化处理的虚拟协方差矩阵Rc的表达式，使用一维字典来估计所述至少一个信源的DOA的数值。5.根据权利要求4所述的平面阵列DOA估计方法，其特征在于，所述步骤3-1包括:通过以下公式对虚拟协方差矩阵Re的矩阵表达式进行矢量化处理，Seqv=[RsTs，Rs2TS，…，RsNpTs]，Reqv=E[seqvseqvH]，其中，上标H表示共辄转置运算，r为等效接收矢量，为等效阵列流型矩阵，符号©表示Khatri-Rao乘积，u为与真实相位β相关的等效信源矢量，Seqv为等效信号矢量，Ts为等效抽样周期，Np是等效快拍数，ITv为对角矩阵，其第k个对角元素为u，矩阵ΦΙ^ν亦为对角矩阵。6.根据权利要求5所述的平面阵列DOA估计方法，其特征在于，所述步骤3-2包括：步骤3-2-1:基于一维字典，θ2，…，0D}DΚ进行迭代运算，直至获得与一维字典中各栅格点--对应的能量值矢量P=[Ρι,Ρι,…，PD]T的估计值矢量，对于第i次迭代，包括以下步骤：固定Θ1_1，按照下式对能量值矢量P进行更新，得到P1，固定P1，按照公式1更新©⑷，或者，按照公式更新α1，其中，▽,为变量α的梯度变化最大的方向，μα为步长，步骤3-2-2:将估计值矢量中的第k个非零项所对应的角度值作为第k个夹角ak的角度估计值，并且按照下式获取对应的的角度估计值，其中，4表示中的第k个非零元素。7.—种平面阵列DOA估计设备，所述设备包括平面阵列、处理器和存储有可执行指令的存储器，其特征在于，所述平面阵列包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对，所述处理器执行所述可执行指令来完成根据权利要求1至6中任一项所述的方法中的步骤。8.—种平面阵列DOA估计设备，其特征在于，包括：平面阵列模块，用于接收来自至少一个信源的信号，所述平面阵列模块包括至少一个平行互质阵列，所述平行互质阵列包括相互平行的互质线性子阵列对；虚拟接收信号计算和构建模块，其与平面阵列模块连接，用于基于所述互质线性子阵列对的实际接收信号计算和构建所述互质线性子阵列对的虚拟接收信号；DOA估计模块，其与虚拟接收信号计算和构建模块连接，用于基于所述互质线性子阵列对的虚拟接收信号估计所述至少一个信源的DOA的数值。

百度查询： 山东农业大学 一种平面阵列DOA估计方法及设备

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD