一种稀疏双平行线阵及二维波达方向估计方法技术

本发明专利技术公开了一种稀疏双平行线阵及二维波达方向估计方法。本发明专利技术的线阵包括两个SULA，且阵元数相差1，阵元数多的子阵还包括一个距原点的间距小于或等于0.5λ的辅助阵元，SULA阵元间距为λ，且子阵间的间距为0.5λ。对应的DOA估计方法为：基于接收数据的自相关矩阵得到线性变换矩阵，再基于维度与信号个数相同的单位阵构造传播矩阵Pc；通过对Pc进行上下划分和相应的矩阵变换，求得包含入射信号与y轴的夹角的余弦值；基于两个SULA所对应的总的阵列流型矩阵，求得入射信号与x轴的夹角的两组余弦值；求结合基于辅助阵元的接收数据进行去模糊；利用入射信号与x轴和y轴的夹角的余弦值求俯仰角和方位角的估计值。本发明专利技术用于无线通信，其结构简单，低成本，高测向精度。

【技术实现步骤摘要】

本专利技术属于无线移动通信
，特别是涉及一种基于稀疏均匀线阵的平行线性阵列及其二维波达方向(DOA)估计方法。
技术介绍
随着空分多址技术和智能天线技术的发展，利用信号的波达方向(DOA)完成信号的空域捕获和跟踪吸引了大量国内外学者的研究。近二十年来，基于无源阵列接收信号的二维(2D)波达方向(DOA)估计问题越来越得到广泛的关注，尤其是其在军事雷达、声呐、导航、通信等领域的应用。现有的二维波达方向估计方法大多是基于由阵元间距等于半倍波长的简化面阵构成的，包括L形阵列、十字形阵列和双平行线阵等。其中，双平行线阵由于结构简单、易于实现、具有较强的方法适用性等优点得到了广泛的关注和应用。近几十年来，人们已经做了很多的利用双平行线阵估计2-DDOA的研究并提出了大量的算法，比较常见的主要有以下三种方案：(1)利用两个阵元间距为半倍波长的均匀线阵(ULA)构造的双平行线阵，其中，第一子阵比第二子阵多一个阵元。信号的俯仰角和方位角是基于传播算子法求得的，但是由于所采用的是常规的ULA，其测向精度并不是很好，此外，该方案估计边缘角度时性能很差，具体参考文献1：J.Li,X.Zhang,H.Chen,ImprovedtwodimensionalDOAestimationalgorithmfortwoparalleluniformlineararraysusingpropagatormethod,SignalProcess.vol.92,no.12,pp.3032–3038,2012.；(2)利用两个阵元数相等的双平行线阵和一系列的线性变换操作，将复数域求解问题转换到实数域求解，一定程度上降低了计算复杂度，但是估计精度的提高也是有限的，具体参考文献2：YANGTao,YUXiaohong,GANLu.New2-DUnitaryESPRITalgorithmforDOAestimation.ComputerEngineeringandApplications,2012,48(15)：123-128；(3)利用两个阵元数相同的线性子阵，将二维谱峰搜索问题简化为两个一维的处理问题，虽然其计算复杂度比传统的二维谱峰搜索低一些，但是仍然有一次谱峰搜索，因此计算量比方案(1)高，此外，该方案的估计精度有限，具体参考文献3：ieqiXia，YiZheng,QunWan,andXuegangWang，“DecoupledEstimationof2-DAnglesofArrivalUsingTwoParallelUniformLinearArrays.”IEEETRANSACTIONSONANTENNASANDPROPAGATION,VOL.55,NO.9,SEPTEMBER2007。
技术实现思路
本专利技术的专利技术目的在于：针对现有的双平行线阵在估计算法复杂度，估计精度和系统成本方面的矛盾，提出一种结构简单的基于稀疏线阵的稀疏双平行线阵及其相应的二维波达角(DOA)估计算法，以达到降低计算复杂度、系统成本、简化处理程序、有效提高估计精度。本专利技术的稀疏双平行线阵，包括相互平行的第一子阵和第二子阵，两个子阵的间距为0.5λ，两个子阵的一个端点阵元位置对齐，将对齐的第一子阵的端点阵元定义为原点，第一子阵包括M+1个以阵元间距λ均匀线性分布的阵元，以及一个距原点的间距小于或等于0.5λ的辅助阵元(优选0.5λ)，即第一子阵的阵元总数为M+2；第二子阵包括M个以阵元间距λ均匀线性分布的阵元，其中λ表示载波波长。同时，本专利技术还公开了一种用于本专利技术的稀疏双平行线阵的二维波达方向估计方法，包括下列步骤：步骤1：稀疏双平行线阵接收K个不相关信源的入射信号，得到两个子阵中各阵元的接收数据，其中K小于或等于M；步骤2：获取传播矩阵Pc：步骤201：将第一子阵中除辅助阵元外的所有阵元的接收数据记为x1(t)，将第二子阵的所有阵元的接收数据记为x2(t)，堆叠x1(t)和x2(t)得到接收数据z(t)，即计算接收数据z(t)在N次采样下的互相关矩阵并将的前K列记为后2M+1-K列记为其中N为预设采样次数；步骤202：根据得到矩阵其中符号(·)+表示M-P广义逆；步骤203：堆叠矩阵和K×K维的单位矩阵I，得到(2M+1)×K维传播矩阵Pc；步骤3：用α、β表示入射信号的入射方向在水平方向、竖直方向的夹角，通过下标k(k＝1,2…,K)区分不同的入射信号，即αk、βk，分别获取cosβk、cosαk的估计值301：将传播矩阵Pc的前M行记为P1，最后M行记为P3，然后利用P1和P3之间的线性关系求得包含β全部信息的矩阵对矩阵F1进行特征分解，得到K个特征值λ1,…,λK和K个特征向量，将特征向量按列排放构成矩阵用于的计算；基于矩阵F1的K个特征值得到其中子阵间距d＝0.5λ，即302：基于矩阵得到阵列流型矩阵的估计值将的前M行记为的第M+2至2M行记为根据得到矩阵B1，其中符号(·)H表示矩阵的共轭转置将的第2至M+1行记为的第M+3至2M+1行记为根据得到矩阵B2；根据得到矩阵而的对角线元素与F1的特征值是一一对应，即：由的第k个对角线元素求得的和F1的第k个特征值求得的是自动配对的。步骤303：获取cosαk的两个备选估计值：根据公式计算第第一备选估计值cosα′k，其中阵元间距d1＝λ，即根据公式计算空间相位差若则第第二备选估计值否则其中阵元间距d1＝λ；步骤304：对cosαk的两个备选估计值去模糊，得到将第一子阵的所有阵元的接收数据记为x(t)，并计算x(t)在N次采样下的自相关矩阵对进行特征值分解，得到M+2个特征值，将前M+2-K个最小特征值对应的特征向量记为噪声子空间Un；分别将K个α′k、α″k代入第一子阵的阵列流型矩阵中，得到和其中a(αk′)和a(α″k)表示导向矢量；若则否则步骤4：根据计算K个入射信号的俯仰角和方位角估计值俯仰角估计值方位角估计值其中综上所述，由于采用了上述技术方案，本专利技术的有益效果是：阵列结构简单，易于实现，可以实现低成本，低运算量，高测向精度的二维DOA估计，该方法可以应用于雷达，声呐及无线通信等领域，解决现有常规阵列所不能解决的问题。附图说明图1是本专利技术提出的基于SULA和SLA构造的稀疏L阵的阵列结构；图2是本专利技术所提稀疏双平行线阵及算法与现有算法的俯仰角和方位角的估计值的均方根误差随信噪比变化的轨迹图；图3是本专利技术所提稀疏双平行线阵及算法与现有算法的俯仰角和方位角的估计值的均方根误差随采样快拍数变化的轨迹图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面结合实施方式和附图，对本专利技术作进一步地详细描述。实施例参见图1，本专利技术的稀疏双平行线阵包括子阵间距为d的子阵1和子阵2，其中子阵1位于x轴包含M1＝M+2个阵元，其中稀疏均匀线阵(SULA)的阵元数为M+1，SULA的阵元间距d1，子阵1还包括一个距原点的间距为d的辅助阵元；子阵2平行于子阵1，包括M2＝M个阵元，其中M个阵元皆为SULA阵元，阵元间距为d1，且左端点阵元与子阵1的左端点对齐。其中d＝0.5λ，d1＝λ，λ为载波波长。本实施例中M＝5。将图1所示的稀疏双平行线阵用于对K＝2个不相关信源的入射信号的二维DOA估计，用α、β表本文档来自技高网...

【技术保护点】
一种稀疏双平行线阵，包括相互平行的第一子阵和第二子阵，两个子阵的间距为0.5λ，两个子阵的一个端点阵元位置对齐，将对齐的第一子阵的端点阵元定义为原点，其特征在于，第一子阵包括M+1个以阵元间距λ均匀线性分布的阵元；第一子阵还包括一个距原点的间距小于或等于0.5λ的辅助阵元；第二子阵包括M个以阵元间距λ均匀线性分布的阵元，其中λ表示载波波长。

【技术特征摘要】
1.一种稀疏双平行线阵，包括相互平行的第一子阵和第二子阵，两个子阵的间距为0.5λ，两个子阵的一个端点阵元位置对齐，将对齐的第一子阵的端点阵元定义为原点，其特征在于，第一子阵包括M+1个以阵元间距λ均匀线性分布的阵元；第一子阵还包括一个距原点的间距小于或等于0.5λ的辅助阵元；第二子阵包括M个以阵元间距λ均匀线性分布的阵元，其中λ表示载波波长。2.一种用于权利要求1所述的稀疏双平行线阵的二维波达方向估计方法，其特征在于，包括下列步骤：步骤1：稀疏双平行线阵接收K个不相关信源的入射信号，得到两个子阵中各阵元的接收数据，其中K小于或等于M；步骤2：获取传播矩阵Pc：步骤201：将第一子阵中除辅助阵元外的所有阵元的接收数据记为x1(t)，将第二子阵的所有阵元的接收数据记为x2(t)，堆叠x1(t)和x2(t)得到接收数据z(t)；计算接收数据z(t)在N次采样下的互相关矩阵并将的前K列记为后2M+1-K列记为其中N为预设采样次数；步骤202：根据得到矩阵其中符号(·)+表示M-P广义逆；步骤203：堆叠矩阵和K×K维的单位矩阵I，得到(2M+1)×K维传播矩阵Pc；步骤3：用αk、βk表示K个入射信号的入射方向在水平方...

【专利技术属性】
技术研发人员：杨雨轩，郑植，孟会鹏，葛琰，刘柯宏，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51