一种利用迭代相位补偿技术估计任意阵列波达方向方法技术

本发明专利技术涉及一种利用迭代相位补偿技术估计任意阵列波达方向方法，属于波达方向估计技术领域。该方法先根据实际阵列阵元位置构造一个虚拟均匀线阵作为参考阵列，然后引入迭代相位补偿操作，即将上一次经典ESPRIT算法输出的估计角度作为相位补偿角，根据实际阵列和参考阵列的位置差计算出相位补偿矩阵，并用该相位补偿矩阵对原始信号子空间进行相位补偿，从而对更新后的信号子空间迭代施行经典ESPRIT算法，直到获得稳定收敛的角度估计，并将其作为信号波达方向估计。本发明专利技术能快速获得基于任意非规则线阵波达方向估计，并且估计结果具有统计无偏性。

An iterative phase compensation technique for estimating the direction of arrival of arbitrary array

The invention relates to an iterative phase compensation technique for estimating the direction of the arbitrary array of arrival (DOA), which belongs to the field of the direction of arrival estimation. Firstly, according to the actual array position to construct a virtual uniform linear array as a reference array, then introduce the iterative phase compensation operation, to estimate angle a classic ESPRIT algorithm as the phase angle compensation, according to the actual array and reference array position error matrix to calculate phase compensation, phase compensation of the original signal the subspace and the phase compensation matrix, thus the implementation of the classical ESPRIT algorithm for signal sub space iteration after the update, until the convergence angle estimation, and as signal DOA estimation. The invention can quickly obtain the estimation of the direction of arrival based on arbitrary irregular linear array, and the estimated results have statistical unbiased.

【技术实现步骤摘要】
一种利用迭代相位补偿技术估计任意阵列波达方向方法
本专利技术属于波达方向估计
，涉及一种利用迭代相位补偿技术估计任意阵列波达方向方法。
技术介绍
空间信号波达方向(DoA，Direction-of-Arrival)估计，或者称为空间谱(SpatialSpectrum)估计，指的是利用天线系统接收空间来波信号并通过接收数据来确定信号方向的一种技术。DoA估计是现代阵列信号处理领域两大分支之一，具有较大的理论和实用价值。利用相控阵列，结合先进的DoA估计技术来对空间目标进行定位，已经在雷达、声纳、气象预测、天文观测等领域有着大量应用。新一代宽带无线通信系统，如4G/5G，利用“智能天线”(即阵列天线)来准确估计用户的方向，从而利用波束成形技术生成窄的无线波束与之通信，能有效降低系统能耗，从而提高系统容量。此外在地震源的确定，射线断层成像等都需要用到DoA估计技术。传统的DoA估计在雷达技术中又称为测角，即利用雷达天线发射较窄的无线电波对目标空间进行扫描，然后通过接收回波信号的有无来判断目标的存在与否，同时也可估计出该回波的DOA参数，即目标的方向。利用窄波束扫描再进行估计的方法称为Bartlett波束形成法，这种经典的波束形成技术获得的DoA估计的空间分辨率受瑞利限限制，即不能分辨出同一个波束宽度内的多个空间目标；而后基于最小方差无畸变响应(MVDR，MinimumVarianceDistortionlessResponse)原理的Capon算法或者最大熵法(MEM，MaximumEntropyMethod)等自适应波束形成技术取得了高于瑞利限的空间分辨能力，被称为高分辨率空间谱估计技术；而由多重信号分类(MUSIC，MultipleSignalClassification)算法，利用旋转不变技术估计信号参数(ESPRIT，EstimationofSignalParametersViaRotationalInvarianceTechniques)算法以及最大似然(ML，MaximumLikelihood)算法为代表的超分辨率DoA估计算法在DoA估计领域内具有里程碑的意义。通过这些算法，可以获得针状的空间谱估计，空间分辨能力获得数量级的提升，因此超分辨率DoA估计算法获得了大量的关注和研究。一般来讲，ML类DoA估计算法估计精度高，均方误差小，但要涉及高维参数搜索，运算量较大；虽然有牛顿法，准牛顿法，交替迭代法等优化方法，但容易获得局部最优解，并且运算量仍然较高；MUSIC类算法可以普遍适用于各种阵列结构，它将ML算法的多维参数搜索转变为1维参数搜索(1D,1-DimentionDoA估计)或者2维参数搜索(2DDoA估计)。对于非相干信号源，已经证明MUSIC算法是ML算法的大样本实现，但MUSIC方法基于穷尽搜索策略，运算量也相当大，不利于实时应用；经典ESPRIT算法基于代数形式的解，运算效率最高，但要求阵列具有移不变结构，因此限制了该算法的应用范围。除了规则的阵列结构，实际中还存在着大量的非规则阵列，即任意结构的阵列，比如代表着高精技术的共形相控阵，一般就具有非规则的阵列结构。针对这种任意结构的阵列的DoA估计，要想得到高效率的DoA估计结果，目前主要有两种方法可供利用：一是采用阵列内插技术(AIT，ArrayInterpolationTechnique)，即将实际阵列流形在一定的空间角域内映射为均匀阵列流形，然后再使用ESPRIT进行估计，这种方法存在着映射误差，并且在事先确定的空间角域外误差更为明显；另一种方法称为流形分离技术(MST，ManifoldSeperationTechnique)，即将实际阵列的流形矢量进行Jacobi-Anger展开，再进行模式截断，从而将实际的流形矩阵分解为一个采样矩阵与一个范德蒙结构矩阵的乘积，从而利用求根MUSIC(Root-MUSIC)得到DoA估计的代数形式解。MST存在有截断误差，这将导致DoA估计存在偏差，同时对于稀疏阵列，为获得小的估计偏差需要截取一个较大模式数，最终会求解一个高阶代数方程，导致运算量较大。因此，设计一种能适用于非规则任意结构阵列，并且能获得统计无偏结果的DoA快速估计算法，就具有较大的理论和工程实际价值。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种利用迭代相位补偿技术估计任意阵列波达方向方法，基于迭代相位补偿技术，输出代数形式DoA估计式，能高效获得统计无偏的DoA估计结果。为达到上述目的，本专利技术提供如下技术方案：一种利用迭代相位补偿技术估计任意阵列波达方向方法，该方法包括以下步骤：S1：由实际阵列阵元位置坐标构造一个虚拟均匀线阵作为参考阵列，并确定实际阵列与虚拟参考阵列的位置差向量；S2：由接收到的阵列快拍数据获得信号子空间的估计；S3：初始化相位补偿角，初始化相位补偿矩阵，设定阈值a，初始化相位补偿角集合，初始化输出集合；S4：将相位补偿矩阵与信号子空间矩阵相乘，完成相位补偿操作，并对经过相位补偿后的信号子空间施行经典的ESPRIT算法得到角度估计值，计算每个输出估计角度与输入的相位补偿角的绝对差，并比较该绝对差和设定阈值的大小，小于阈值则将该估计角度放入输出集合中，否则将该估计角度放入相位补偿角集合中并顺序编号；S5：判断输出集合中的元素个数是否与信号个数相等，不相等则继续执行S6，否则转向执行S7；S6：从相位补偿角集合中按顺序取出一个相位补偿角，根据该相位补偿角结合位置差向量计算出相位补偿矩阵并转向S4；S7：结束迭代验证并将输出集合中的角度估计值作为波达方向的最终估计值。进一步，在步骤S1中，所述由实际阵列阵元位置坐标构造一个虚拟均匀线阵作为参考阵列具体为：按实际阵列阵元x坐标从小到大顺序编号，得到阵元位置向量，令虚拟参考均匀线阵的两个端阵元位置与实际阵列的两个端阵元位置重合，虚拟参考阵列的其余阵元位置在端阵元之间均匀布置，虚拟参考阵元也按x坐标从小到大顺序编号，得到虚拟阵列阵元位置向量，确定实际阵列与虚拟参考阵列的位置差向量，最后计算出虚拟参考阵相邻阵元之间的距离。进一步，在步骤S4中，所述对经过相位补偿后的信号子空间施行经典的ESPRIT算法得到角度估计具体为：从经过相位补偿后的信号子空间选择出两个子阵列的信号子空间，然后再利用最小二乘算法计算这两个子阵列信号空间的拟合矩阵，并对拟合矩阵进行特征分解，获得拟合矩阵的特征值，最后利用每个特征值的相角参数获得角度值估计。进一步，在步骤S6中，所述根据取得的该相位补偿角结合位置差向量计算出相位补偿矩阵并转向S4具体为：从相位补偿角集合中顺序取得一个相位补偿角，并结合S1中的位置差向量根据相位补偿矩阵计算公式计算出基于该相位补偿角的相位补偿矩阵并转向执行步骤S4。本专利技术的有益效果在于：(1)本专利技术针对不规则阵列引入一个虚拟均匀线阵作为参考阵列，该虚拟阵列由实际阵列的阵元位置决定，并由实际阵列与虚拟参考阵列对应阵元之间的位置差以及一个相位补偿角度来定义相位补偿矩阵。相位补偿矩阵是一个酉对角矩阵，利用该酉矩阵对原始的信号子空间进行相位补偿，从而对相控阵列阵元空间位置的不规则空间采样进行补偿，使得经过补偿后的信号子空间逐步具有经典ESPRIT算法所要求的移不变性本文档来自技高网...

【技术保护点】
一种利用迭代相位补偿技术估计任意阵列波达方向方法，其特征在于：该方法包括以下步骤：S1：由实际阵列阵元位置坐标构造一个虚拟均匀线阵作为参考阵列，并确定实际阵列与虚拟参考阵列的位置差向量；S2：由接收到的阵列快拍数据获得信号子空间的估计；S3：初始化相位补偿角，初始化相位补偿矩阵，设定阈值a，初始化相位补偿角集合，初始化角度估计输出集合；S4：将相位补偿矩阵与信号子空间矩阵相乘，完成相位补偿操作，并对经过相位补偿后的信号子空间施行经典的ESPRIT算法得到角度估计值，计算每个输出估计角度与输入相位补偿角的绝对差，并比较该绝对差和设定阈值的大小，小于阈值则将该估计角度放入输出集合中，否则将该估计角度放入相位补偿角集合中并顺序编号；S5：判断输出集合中的元素个数是否与信号个数相等，不相等则继续执行S6，否则转向执行S7；S6：从相位补偿角集合中按顺序取出一个相位补偿角，根据该相位补偿角结合位置差向量计算出相位补偿矩阵并转向S4；S7：结束迭代验证并将输出集合中的角度估计值作为波达方向的最终估计值。

【技术特征摘要】
1.一种利用迭代相位补偿技术估计任意阵列波达方向方法，其特征在于：该方法包括以下步骤：S1：由实际阵列阵元位置坐标构造一个虚拟均匀线阵作为参考阵列，并确定实际阵列与虚拟参考阵列的位置差向量；S2：由接收到的阵列快拍数据获得信号子空间的估计；S3：初始化相位补偿角，初始化相位补偿矩阵，设定阈值a，初始化相位补偿角集合，初始化角度估计输出集合；S4：将相位补偿矩阵与信号子空间矩阵相乘，完成相位补偿操作，并对经过相位补偿后的信号子空间施行经典的ESPRIT算法得到角度估计值，计算每个输出估计角度与输入相位补偿角的绝对差，并比较该绝对差和设定阈值的大小，小于阈值则将该估计角度放入输出集合中，否则将该估计角度放入相位补偿角集合中并顺序编号；S5：判断输出集合中的元素个数是否与信号个数相等，不相等则继续执行S6，否则转向执行S7；S6：从相位补偿角集合中按顺序取出一个相位补偿角，根据该相位补偿角结合位置差向量计算出相位补偿矩阵并转向S4；S7：结束迭代验证并将输出集合中的角度估计值作为波达方向的最终估计值。2.根据权利要求1所述的一种利用迭代相位补偿技术估计任意阵列波达方向方法，其特征在于：在步骤S1中，所述由实际阵列阵元位置坐标构造一个虚拟...

【专利技术属性】
技术研发人员：刘松，庞育才，翁明江，余翔，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：重庆,50