基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法技术

本发明专利技术公开了一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，主要解决现有技术中虚拟阵列的非均匀性所导致的信息损失问题。其实现步骤是：接收端架构互质阵列；利用互质阵列接收入射信号并建模；计算互质阵列接收信号所对应的等价虚拟信号；构造内插虚拟阵列并建模；构造内插虚拟阵列的多采样快拍信号及其采样协方差矩阵；构造投影矩阵并定义与该投影矩阵相关的投影运算；根据原始虚拟阵列中的全部信息构建参考协方差矩阵，设计基于内插虚拟阵列协方差矩阵Toeplitz化重建的优化问题并求解；根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。本发明专利技术提高了信号波达方向估计的自由度及分辨率，可用于无源定位和目标探测。

【技术实现步骤摘要】
基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法
本专利技术属于信号处理
，尤其涉及对雷达信号、声学信号及电磁信号的波达方向估计，具体是一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，可用于无源定位和目标探测。
技术介绍
波达方向(Direction-of-Arrival,DOA)估计是阵列信号处理领域的一个重要分支，它是指利用阵列天线接收空域信号，并通过现代信号处理技术和各类优化方法实现对接收信号统计量的有效处理，从而实现信号的DOA估计，在雷达、声呐、语音、无线通信等领域有着重要的应用价值。DOA估计方法的自由度是指其能够估计的入射信号源的个数。现有的DOA估计方法通常采用均匀线性阵列进行信号的接收与建模，但是基于均匀线性阵列方法的自由度受限于实际天线阵元个数。具体而言，对于一个包含L个天线阵元的均匀线性阵列，其自由度为L-1。因此，当某个空域范围内入射信号源的个数大于阵列中天线阵元的个数时，现有采用均匀线性阵列的方法将无法进行有效的DOA估计。互质阵列能够在天线阵元个数一定的前提下增加DOA估计的自由度，因而受到了学术界的广泛关注。作为互质采样技术在空间域上的一个典型表现形式，互质阵列提供了一个系统化的稀疏阵列架构方案，并能够突破传统均匀线性阵列自由度受限的瓶颈，实现DOA估计方法自由度性能的提升。现有的基于互质阵列的DOA估计方法主要通过利用质数的性质将互质阵列推导到虚拟域，并形成等价虚拟均匀线性阵列接收信号以实现DOA估计。由于虚拟阵列中包含的虚拟阵元数大于实际的天线阵元数，自由度因此得到了有效的提升。但是由于从互质阵列推导而来的虚拟阵列属于非均匀阵列，因此很多现有基于均匀线性阵列的信号处理方法无法直接应用于虚拟阵列等价接收信号的DOA估计。当前采用互质阵列的DOA估计方法常用的一个解决方案是，仅利用虚拟阵列中连续的阵元部分形成一个虚拟均匀线阵以进行DOA估计，但是这造成了部分原始信息的丢失和相关估计性能的降低。
技术实现思路
本专利技术的目的在于针对上述现有技术存在的不足，提出一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，充分利用了非均匀虚拟阵列所提供的全部信息，从而提高了DOA估计的自由度与分辨率。本专利技术的目的是通过以下技术方案来实现的：一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，包含以下步骤：(1)接收端使用M+N-1个天线，并按照互质阵列结构进行架构；其中M与N为互质整数；(2)假设有K个来自θ1,θ2,…,θK方向的远场窄带非相干信号源，则(M+N-1)×1维互质阵列接收信号x(t)可建模为：其中，sk(t)为信号波形，n(t)为与各信号源相互独立的噪声分量，a(θk)为θk方向的导引矢量，表示为：其中，pid,i＝1,2,…,M+N-1表示互质阵列中第i个物理天线阵元的实际位置，且p1＝0；d为入射窄带信号波长λ的一半，即d＝λ/2，[·]T表示转置操作。共采集T个采样快拍，得到采样协方差矩阵这里，(·)H表示共轭转置；(3)计算互质阵列接收信号所对应的等价虚拟信号：矢量化互质阵列接收信号的采样协方差矩阵获得虚拟阵列等价接收信号v：其中，为(M+N-1)2×K维虚拟阵列导引矩阵，包含K个入射信号源的功率，为噪声功率，iv＝vec(IM+N-1)。这里，vec(·)表示矢量化操作，即把矩阵中的各列依次堆叠以形成一个新的矢量，(·)*表示共轭操作，表示克罗内克积，IM+N-1表示(M+N-1)×(M+N-1)维单位矩阵。矢量v对应的虚拟阵列中各虚拟阵元的位置为去除集合中各位置上重复的虚拟阵元，得到一个非均匀的虚拟阵列其对应的等价虚拟信号vc可通过选取矢量v中相对应的元素获得；(4)构造内插虚拟阵列及其接收信号并建模：首先对于非均匀的虚拟阵列在保留其原有虚拟阵元位置不变的前提下，向其中非连续的位置插入若干虚拟阵元，从而将非均匀虚拟阵列转化为间距为d、阵列孔径与互质阵列相同、且虚拟阵元数目增加的均匀虚拟阵列该内插均匀虚拟阵列共包含个虚拟阵元，其中|·|表示集合的势，其对应的等价虚拟信号vI可通过往矢量vc中插入0获得，插入0的位置与中插入的虚拟阵元的位置相对应；(5)构造内插虚拟阵列多采样快拍信号及其采样协方差矩阵：将切割为LI个长度为LI的连续子阵列，其中相应地，内插虚拟阵列的多采样快拍信号可通过截取矢量vI中对应的元素获得，即：vI,l,l＝1,2,…,LI由vI中第LI+1-l到第2LI-l个元素组成。接着，VI的采样协方差矩阵Rv可以由如下方式得到：其中，<vI>i表示位置为id的虚拟阵元所对应的等价接收信号；(6)构造投影矩阵并定义投影运算：投影矩阵P的维度与Rv相同，如果矩阵Rv中某个元素为0，则投影矩阵P中相同位置的元素值也为0；否则投影矩阵P中的元素值为1。定义为投影运算，其中括号内变量为与P维度相同的矩阵，投影运算通过变量矩阵中每一个元素与投影矩阵P中相应位置上的元素一一相乘实现，得到一个与矩阵P维度相同的矩阵；(7)设计基于内插虚拟阵列协方差矩阵Toeplitz化重建的优化问题并求解：根据接收信号理论协方差矩阵的Toeplitz性，利用(5)中得到的内插虚拟阵列协方差矩阵Rv作为参考值，寻找一个与其差异最小的低秩Toeplitz矩阵作为接收信号的协方差矩阵，可构建如下以矢量z为变量的优化问题：其中，是以矢量z为第一列的厄米特对称Toeplitz矩阵；∈为阈值常数，用于约束协方差矩阵的重建误差；保证了重建的协方差矩阵满足半正定的条件；‖·‖F表示Frobenius范数；rank(·)表示矩阵的秩。将上述非凸优化问题转化为凸优化问题，并求得最优化值相应地，重建的Toeplitz矩阵为内插虚拟阵列协方差矩阵；(8)根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。进一步地，步骤(1)所述的互质阵列结构可具体描述为：首先选取一对互质整数M、N；然后，构造一对稀疏均匀线性子阵列，其中第一个子阵列包含M个间距为Nd的天线阵元，其位置为0,Nd,…,(M-1)Nd，第二个子阵列包含N个间距为Md的天线阵元，其位置为0,Md,…,(N-1)Md；接着，将两个子阵列按照首个阵元重叠的方式进行子阵列组合，获得实际包含M+N-1个天线阵元的非均匀互质阵列架构。进一步地，步骤(7)中所构建的非凸优化问题可通过凸松弛技术，将优化问题目标函数中的矩阵的秩最小化操作替换为矩阵的迹最小化操作，得到以下以矢量z为变量的凸优化问题：其中，Tr(·)表示矩阵的迹。进一步地，步骤(7)中所构建的非凸优化问题可转化为如下以矢量z为变量的凸优化问题：其中μ为正则化参数，用于在最小化过程中权衡矩阵重建误差和矩阵的迹。进一步地，步骤(8)中的波达方向估计，可采用以下方法：多重信号分类方法、旋转不变子空间方法、求根多重信号分类方法、协方差矩阵稀疏重建方法等。进一步地，步骤8中，通过多重信号分类方法进行波达方向估计，具体为：画出虚拟域空间谱PMUSIC(θ)：其中d(θ)是LI×1维内插虚拟阵列导引矢量，对应于位置为由0到(LI-1)d的一段虚拟均匀阵列；En是LI×(LI本文档来自技高网...

【技术保护点】
一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，其特征在于，包含以下步骤：(1)接收端使用M+N‑1个天线，并按照互质阵列结构进行架构；其中M与N为互质整数；(2)假设有K个来自θ1,θ2,…,θK方向的远场窄带非相干信号源，则(M+N‑1)×1维互质阵列接收信号x(t)可建模为：

【技术特征摘要】
1.一种基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法，其特征在于，包含以下步骤：(1)接收端使用M+N-1个天线，并按照互质阵列结构进行架构；其中M与N为互质整数；(2)假设有K个来自θ1,θ2,…,θK方向的远场窄带非相干信号源，则(M+N-1)×1维互质阵列接收信号x(t)可建模为：其中，sk(t)为信号波形，n(t)为与各信号源相互独立的噪声分量，a(θk)为θk方向的导引矢量，表示为：其中，pid,i＝1,2,…,M+N-1表示互质阵列中第i个物理天线阵元的实际位置，且p1＝0；d为入射窄带信号波长λ的一半，即d＝λ/2，[·]T表示转置操作。共采集T个采样快拍，得到采样协方差矩阵这里，(·)H表示共轭转置；(3)计算互质阵列接收信号所对应的等价虚拟信号：矢量化互质阵列接收信号的采样协方差矩阵获得虚拟阵列等价接收信号v：其中，为(M+N-1)2×K维虚拟阵列导引矩阵，包含K个入射信号源的功率，为噪声功率，iv＝vec(IM+N-1)。这里，vec(·)表示矢量化操作，即把矩阵中的各列依次堆叠以形成一个新的矢量，(·)*表示共轭操作，表示克罗内克积，IM+N-1表示(M+N-1)×(M+N-1)维单位矩阵。矢量v对应的虚拟阵列中各虚拟阵元的位置为去除集合中各位置上重复的虚拟阵元，得到一个非均匀的虚拟阵列其对应的等价虚拟信号vc可通过选取矢量v中相对应的元素获得；(4)构造内插虚拟阵列及其接收信号并建模：首先对于非均匀的虚拟阵列在保留其原有虚拟阵元位置不变的前提下，向其中非连续的位置插入若干虚拟阵元，从而将非均匀虚拟阵列转化为间距为d、阵列孔径与互质阵列相同、且虚拟阵元数目增加的均匀虚拟阵列该内插均匀虚拟阵列共包含个虚拟阵元，其中|·|表示集合的势，其对应的等价虚拟信号vI可通过往矢量vc中插入0获得，插入0的位置与中插入的虚拟阵元的位置相对应；(5)构造内插虚拟阵列多采样快拍信号及其采样协方差矩阵：将切割为LI个长度为LI的连续子阵列，其中相应地，内插虚拟阵列的多采样快拍信号可通过截取矢量vI中对应的元素获得，即：vI,l,l＝1,2,…,LI由vI中第LI+1-l到第2LI-l个元素组成。接着，VI的采样协方差矩阵Rv可以由如下方式得到：其中，<vI>i表示位置为id的虚拟阵元所对应的等价接收信号；(6)构造投影矩阵并定义投影运算：投影矩阵P的维度与Rv相同，如果矩阵Rv中某个元素为0，则投影矩阵P中相同位置的元素值也为0；否则投影矩阵P中的元素值为1。定义为投影运算，其中括号内变量为与P维度相同的矩阵，投影...

【专利技术属性】
技术研发人员：周成伟，樊星，史治国，陈积明，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33