CN107102291B - 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 - Google Patents
基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 Download PDFInfo
- Publication number
- CN107102291B CN107102291B CN201710302902.4A CN201710302902A CN107102291B CN 107102291 B CN107102291 B CN 107102291B CN 201710302902 A CN201710302902 A CN 201710302902A CN 107102291 B CN107102291 B CN 107102291B
- Authority
- CN
- China
- Prior art keywords
- array
- virtual
- interpolation
- virtual array
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/78—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using electromagnetic waves other than radio waves
- G01S3/782—Systems for determining direction or deviation from predetermined direction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,主要解决现有技术中虚拟阵列的非均匀性所导致的信息损失问题。其实现步骤是:接收端架构互质阵列;利用互质阵列接收入射信号并建模;计算互质阵列接收信号所对应的等价虚拟信号;构造内插虚拟阵列并建模;构造内插虚拟阵列的多采样快拍信号及其采样协方差矩阵;构造投影矩阵并定义与该投影矩阵相关的投影运算;设计基于内插虚拟阵列信号协方差矩阵核范数最小化的优化问题并求解;根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。本发明提高了波达方向估计的自由度及分辨率,可用于无源定位和目标探测。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及对雷达信号、声学信号及电磁信号的波达方向估计,具体是一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,可用于无源定位和目标探测。
背景技术
波达方向(Direction-of-Arrival,DOA)估计是阵列信号处理领域的一个重要分支,它是指利用阵列天线接收空域信号,并通过现代信号处理技术和各类优化方法实现对接收信号统计量的有效处理,从而实现信号的DOA估计,在雷达、声呐、语音、无线通信等领域有着重要的应用价值。
DOA估计方法的自由度是指其能够估计的入射信号源的个数。现有的DOA估计方法通常采用均匀线性阵列进行信号的接收与建模,但是基于均匀线性阵列方法的自由度受限于实际天线阵元个数。具体而言,对于一个包含L个天线阵元的均匀线性阵列,其自由度为L-1。因此,当某个空域范围内入射信号源的个数大于或等于阵列中天线阵元的个数时,现有采用均匀线性阵列的方法将无法进行有效的DOA估计。
互质阵列能够在天线阵元个数一定的前提下增加DOA估计的自由度,因而受到了学术界的广泛关注。作为互质采样技术在空间域上的一个典型表现形式,互质阵列提供了一个***化的稀疏阵列架构方案,并能够突破传统均匀线性阵列自由度受限的瓶颈,实现DOA估计方法自由度性能的提升。现有的基于互质阵列的DOA估计方法主要通过利用质数的性质将互质阵列推导到虚拟域,并形成等价虚拟均匀线性阵列接收信号以实现DOA估计。由于虚拟阵列中包含的虚拟阵元数大于实际的天线阵元数,自由度因此得到了有效的提升。但是由于从互质阵列推导而来的虚拟阵列属于非均匀阵列,因此很多现有基于均匀线性阵列的信号处理方法无法直接应用于虚拟阵列等价接收信号以实现有效的DOA估计。当前采用互质阵列的DOA估计方法常用的一个解决方案是,仅利用虚拟阵列中连续的阵元部分形成一个虚拟均匀线阵以进行DOA估计,但是这造成了部分原始信息的丢失和相关估计性能的降低。
同时,目前众多DOA估计方法在优化问题的设计过程中,需要预先设置信号假定波达方向的空间网格点。随着对波达方向估计结果精度要求的提高,这些DOA估计方法需要预先设置的空间网格点将变得越来越密集,这导致了计算复杂度的急剧增加。不仅如此,在实际情况中,难免会有一些信号的波达方向无法完全落在预先设置的网格点上,从而造成了固有的模型失配误差。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,充分利用了非均匀虚拟阵列所提供的全部信息,并保证了无网格化的波达方向估计,从而提高了DOA估计的自由度与分辨率,并在一定程度上降低了DOA估计的计算复杂度。
本发明的目的是通过以下技术方案来实现的:一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,包含以下步骤:
(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)设计基于内插虚拟阵列信号协方差矩阵核范数最小化的优化问题并求解。利用(5)得到的内插虚拟阵列协方差矩阵Rv作为参考值,寻找一个核范数最小的Toeplitz矩阵作为内插虚拟阵列信号的协方差矩阵,且要求其与Rv的差异小于某一阈值,可构建如下以矢量z为变量的优化问题:
其中,表示的核范数,表示以矢量z为第一列的厄米特对称Toeplitz矩阵;∈为阈值常数,用于约束协方差矩阵的重建误差;保证了重建的协方差矩阵满足半正定的条件;‖·‖F表示Frobenius范数。求解上述凸优化问题可得到最优化值相应地,重建的Toeplitz矩阵为内插虚拟阵列协方差矩阵;
(8)根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。
进一步地,步骤(1)所述的互质阵列结构可具体描述为:首先选取一对互质整数M、N;然后,构造一对稀疏均匀线性子阵列,其中第一个子阵列包含M个间距为Nd的天线阵元,其位置为0,Nd,…,(M-1)Nd,第二个子阵列包含N个间距为Md的天线阵元,其位置为0,Md,…,(N-1)Md;接着,将两个子阵列按照首个阵元重叠的方式进行子阵列组合,获得实际包含M+N-1个天线阵元的非均匀互质阵列架构。
进一步地,步骤(5)所构建的VI的采样协方差矩阵Rv也可以由下述方法等价得到:
进一步地,步骤(7)中的凸优化问题可转化为如下以矢量z为变量的优化问题:
其中μ为正则化参数,用于在最小化过程中权衡矩阵重建误差和z的核范数。
进一步地,步骤(8)中的波达方向估计,可采用以下方法:多重信号分类方法、旋转不变子空间方法、求根多重信号分类方法、协方差矩阵稀疏重建方法等。
进一步地,步骤8中,通过多重信号分类方法进行波达方向估计,具体为:画出虚拟域空间谱PMUSIC(θ):
其中d(θ)是LI×1维内插虚拟阵列导引矢量,对应于位置为由0到(LI-1)d的一段虚拟均匀阵列;En是LI×(LI-K)维矩阵,表示内插虚拟阵列协方差矩阵的噪声子空间;θ是假定的信号波达方向;通过谱峰搜索寻找空间谱PMUSIC(θ)上的峰值,并将这些峰值所对应的响应值从大到小排列,取前K个峰值所对应的角度方向,即为波达方向估计结果。
本发明与现有技术相比具有以下优点:
(1)本发明在互质阵列等价虚拟域上引入了阵列内插的思想,充分利用了虚拟阵列提供的全部信息。通过在非均匀虚拟阵列中内插虚拟阵元的方式构建出均匀线性虚拟阵列,在保留了由原始非均匀虚拟阵列接收到的全部信息的同时,使得构建的虚拟域信号模型满足奈奎斯特采样定律;
(2)本发明基于内插虚拟阵列信号协方差矩阵核范数最小化的思想设计优化问题,在优化问题设计的过程中无需预先定义空间网格点,实现了无网格化的波达方向估计,同时保证了波达方向估计的分辨率以及计算效率;
(3)本发明所提出的基于内插虚拟阵列协方差矩阵重建的优化问题保证了优化求解结果为厄米特对称的Toeplitz矩阵,使得最优解与理论协方差矩阵之间的误差更小。由于均匀线性阵列非相干接收信号的理论协方差矩阵满足Toeplitz结构,因此利用它的Toeplitz特性作为先验约束条件进行协方差矩阵的重建,可以使得重建结果与真实值差异更小,从而提高DOA估计的性能。
附图说明
图1是本发明的方法总体流程框图。
图2是本发明中组成互质阵列的一对稀疏均匀子阵列结构示意图。
图3是本发明中互质阵列的结构示意图。
图4是本发明中内插虚拟阵列的结构示意图。
图5是本发明中内插虚拟阵列分割方法的示意图。
图6是用于体现本发明所提方法自由度性能的空间功率谱示意图。
图7是用于体现本发明所提方法分辨率性能的归一化空间谱示意图。
具体实施方式
以下参照附图,对本发明的技术方案和效果作进一步的详细说明。
对于DOA估计在实际***中的应用,互质阵列由于其可以通过等价虚拟阵列信号的计算和统计信号处理,突破物理阵元数量对自由度的限制而备受关注。但是受限于虚拟阵列的非均匀性,目前很多方法都会选择利用其中连续的虚拟阵元部分进行DOA估计,从而造成了信息损失。同时,很多方法在进行DOA估计之前会预先设置假定波达信号方向的空间网格点,这造成了固有失配误差以及计算复杂度与估计精度之间的矛盾。为了充分利用非均匀虚拟阵列中所包含的所有信息,并避免由于预定义空间网格点所造成的估计分辨率受限问题,本发明提供了一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,参照图1,本发明的实现步骤如下:
步骤一:在接收端使用M+N-1个天线阵元架构互质阵列;首先,选取一组互质整数M、N;然后,参照图2,构造一对稀疏均匀线性子阵列,其中第一个子阵列包含M个间距为Nd的天线阵元,其位置为0,Nd,…,(M-1)Nd;第二个子阵列包含N个间距为Md的天线阵元,其位置为0,Md,…,(N-1)Md;单位间距d取为入射窄带信号波长的一半,即d=λ/2;接着,将两个子阵列的首个天线阵元视为参考阵元,参照图3,将两个子阵的参考阵元重叠以实现子阵列组合,获得实际包含M+N-1个天线阵元的非均匀互质阵列架构。
步骤二:采用互质阵列接收信号并建模。假设有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表示共轭转置操作。
步骤三:计算互质阵列接收信号所对应的等价虚拟信号。矢量化互质阵列接收信号的采样协方差矩阵获得虚拟阵列等价接收信号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,对于非均匀的虚拟阵列在保留其原有虚拟阵元位置不变的前提下,向其中存在孔洞的位置***若干虚拟阵元(如图4中的空心圆所示),从而将非均匀虚拟阵列转化为间距为d、阵列孔径与互质阵列相同、且虚拟阵元数目增加的均匀虚拟阵列内插虚拟阵列共包含个虚拟阵元,其中|·|表示集合的势。内插虚拟阵列对应的等价虚拟信号vI可通过往矢量vc中孔洞的相应位置填入0获得。
步骤五:构造内插虚拟阵列多采样快拍信号及其采样协方差矩阵。参照图5,将内插虚拟阵列切割为LI个长度为LI的连续子阵列,其中
由于中的虚拟阵元以零位对称分布,始终为奇数,故LI为整数。相应地,内插虚拟阵列的多采样快拍信号可通过截取矢量vI中对应的元素获得,即:VI=[vI,1,vI,2,…,vI,LI],其中vI,l,l=1,2,…,LI由vI中第LI+1-l到第2LI-l个元素组成。接着,VI的采样协方差矩阵Rv可以由如下方式得到:
其中,〈vI〉i表示位置为id的虚拟阵元所对应的等价接收信号。由于内插虚拟阵列中虚拟阵元关于零位对称分布,因此其上的等价虚拟接收信号对应于零位呈共轭关系,所以上述采样协方差矩阵也可以通过如下方式等价得到:
步骤六:构造投影矩阵并定义投影运算。由于步骤五所得的协方差矩阵Rv中包含有在步骤四中***的0,因此其相应位置对角线上的元素全部为0。根据这样的结构定义一个与Rv维度相同的投影矩阵P,如果Rv中某一位置上的元素是0,则投影矩阵P中相同位置的元素值也为0;反之则投影矩阵P中相应位置的元素值为1。定义为投影运算,其中括号内变量为与P维度相同的矩阵,投影运算通过变量矩阵的每一个元素与投影矩阵P中相应位置上的元素一一相乘实现,得到一个与矩阵P维度相同的矩阵。
步骤七:设计基于内插虚拟阵列信号协方差矩阵核范数最小化的优化问题并求解。利用步骤五得到的内插虚拟阵列协方差矩阵Rv作为参考值,寻找一个核范数最小的Toeplitz矩阵作为内插虚拟阵列信号的协方差矩阵,且要求其与Rv的差异小于某一阈值,可构建如下以矢量z为变量的优化问题:
其中,表示的核范数,表示以矢量z为第一列的厄米特对称Toeplitz矩阵;∈为阈值常数,用于约束协方差矩阵的重建误差;保证了重建的协方差矩阵满足半正定的条件;‖·‖F表示Frobenius范数。求解上述凸优化问题可得到最优化值上述凸优化问题可转化为以下以矢量z为变量的优化问题:
其中μ为正则化参数,用于在最小化过程中权衡矩阵重建误差和z的核范数。求解上述凸优化问题可得到最优化值相应地,重建的Toeplitz矩阵为内插虚拟阵列协方差矩阵。
步骤八:根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。通过引入经典的方法,如多重信号分类方法、旋转不变子空间方法、求根多重信号分类方法、协方差矩阵稀疏重建方法等对重建的内插虚拟阵列协方差矩阵进行操作,可以求得波达方向估计结果。以多重信号分类方法为例,画出虚拟域空间谱PMUSIC(θ):
其中d(θ)是LI×1维内插虚拟阵列导引矢量,对应于位置为由0到(LI-1)d的一段虚拟均匀阵列;En是LI×(LI-K)维矩阵,表示内插虚拟阵列协方差矩阵的噪声子空间;θ是假定的信号波达方向;通过谱峰搜索寻找空间谱PMUSIC上的峰值,并将这些峰值所对应的响应值从大到小排列,取前K个峰值所对应的角度方向,即为波达方向估计结果。
本发明一方面引入虚拟阵列内插的思想,在推导的原始虚拟阵列基础上内***虚拟阵元,从而将原始的非均匀虚拟阵列转化为虚拟均匀阵列,同时保留了原始非均匀虚拟阵列上的所有信息,避免了因原始虚拟阵列的非均匀性所导致的统计信号处理模型失配及传统方法截取虚拟均匀子阵所导致的信息损失问题;另一方面,引入了基于虚拟阵列信号协方差矩阵核范数最小化的思想来设计优化问题,以重建内插虚拟阵列的协方差矩阵,实现了虚拟域上的无网格化波达方向估计。
下面结合仿真实例对本发明的效果做进一步的描述。
仿真实例1:采用互质阵列接收入射信号,其参数选取为M=3,N=5,即架构的互质阵列共包含M+N-1=7个物理阵元。假定入射窄带信号个数为9,且入射方向均匀分布于-50°至50°这一空间角度域范围内;信噪比设置为30dB,采样快拍数T=500;正则化参数μ设置为0.25。
本发明所提出的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法空间功率谱如图6所示,其中垂直虚线代表入射信号源的实际方向。可以看出,本发明所提方法能够有效分辨这9个入射信号源。而对于传统采用均匀线性阵列的方法,利用7个物理天线阵元最多只能分辨6个入射信号,以上结果体现了本发明所提方法实现了自由度的增加。
仿真实例2:采用互质阵列接收入射信号,其参数同样选取为M=3,N=5,即架构的互质阵列共包含M+N-1=7个物理天线阵元;假定入射窄带信号个数为2,且入射方向为-0.5°至0.5°,其余参数设置与仿真实例1保持一致。由图7所示的归一化空间谱可以看出,本发明所提方法可以有效地分辨出这两个近距离信号的波达方向,说明了本方法良好的分辨率性能。
综上所述,本发明所提方法充分利用了非均匀虚拟阵列上的全部信息,能够在信号源个数大于等于物理天线个数的情况下实现入射信号的无网格化估计,增加了DOA估计的自由度和分辨率。此外,与传统采用均匀线性阵列的方法相比,本发明所提方法在实际应用中所需的物理天线阵元及射频模块也能够相应减少,体现了经济性和高效性。
Claims (6)
1.一种基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于,包含以下步骤:
(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的连续子阵列,其中
相应地,内插均匀虚拟阵列SI的多采样快拍信号可通过截取矢量vI中对应的元素获得,即:vI,l,l=1,2,...,LI由vI中第LI+1-l到第2LI-l个元素组成;接着,VI的采样协方差矩阵Rv可以由如下方式得到:
其中,表示位置为id的虚拟阵元所对应的等价接收信号;
(6)构造投影矩阵并定义投影运算:投影矩阵P的维度与Rv相同,如果矩阵Rv中某个元素为0,则投影矩阵P中相同位置的元素值也为0;反之投影矩阵P中相应位置的元素值为1;定义为投影运算,其中括号内变量为与P维度相同的矩阵,投影运算通过变量矩阵中每一个元素与投影矩阵P中相应位置上的元素一一相乘实现,得到一个与矩阵P维度相同的矩阵;
(7)设计基于内插虚拟阵列信号协方差矩阵核范数最小化的优化问题并求解;利用(5)得到的内插虚拟阵列协方差矩阵Rv作为参考值,寻找一个核范数最小的Toeplitz矩阵作为内插虚拟阵列信号的协方差矩阵,且要求其与Rv的差异小于某一阈值,可构建如下以矢量z为变量的优化问题:
其中,表示的核范数,表示以矢量z为第一列的厄米特对称Toeplitz矩阵;∈为阈值常数,用于约束协方差矩阵的重建误差;保证了重建的协方差矩阵满足半正定的条件;||·||F表示Frobenius范数;求解上述优化问题可得到最优化值相应地,重建的Toeplitz矩阵为内插虚拟阵列协方差矩阵;
(8)根据重建的内插虚拟阵列协方差矩阵进行波达方向估计。
2.根据权利要求1所述的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于:步骤(1)所述的互质阵列结构可具体描述为:首先选取一对互质整数M、N;然后,构造一对稀疏均匀线性子阵列,其中第一个子阵列包含M个间距为Nd的天线阵元,其位置为0,Nd,...,(M-1)Nd,第二个子阵列包含N个间距为Md的天线阵元,其位置为0,Md,...,(N-1)Md;接着,将两个子阵列按照首个阵元重叠的方式进行子阵列组合,获得实际包含M+N-1个天线阵元的非均匀互质阵列架构。
3.根据权利要求1所述的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于:步骤(5)所构建的VI的采样协方差矩阵Rv也可以由下述方法等价得到:
4.根据权利要求1所述的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于:步骤(7)中的凸优化问题可转化为如下以矢量z为变量的优化问题:
其中μ为正则化参数,用于在最小化过程中权衡矩阵重建误差和z的核范数。
5.根据权利要求1所述的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于:步骤(8)中的波达方向估计,可采用以下方法之一:多重信号分类方法、旋转不变子空间方法、求根多重信号分类方法、协方差矩阵稀疏重建方法。
6.根据权利要求1所述的基于虚拟阵列内插的无网格化互质阵列波达方向估计方法,其特征在于:步骤(8)中,通过多重信号分类方法进行波达方向估计,具体为:画出虚拟域空间谱PMUSIC(θ):
其中d(θ)是LI×1维内插虚拟阵列导引矢量,对应于位置为由0到(LI-1)d的一段虚拟均匀阵列;En是LI×(LI-K)维矩阵,表示重建的内插虚拟阵列协方差矩阵的噪声子空间;θ是假定的信号波达方向;通过谱峰搜索寻找空间谱PMUSIC(θ)上的峰值,并将这些峰值所对应的响应值从大到小排列,取前K个峰值所对应的角度方向,即为波达方向估计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710302902.4A CN107102291B (zh) | 2017-05-03 | 2017-05-03 | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710302902.4A CN107102291B (zh) | 2017-05-03 | 2017-05-03 | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107102291A CN107102291A (zh) | 2017-08-29 |
CN107102291B true CN107102291B (zh) | 2019-07-23 |
Family
ID=59657497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710302902.4A Active CN107102291B (zh) | 2017-05-03 | 2017-05-03 | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107102291B (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7023565B2 (ja) * | 2017-10-06 | 2022-02-22 | 日本無線株式会社 | アレーアンテナ装置 |
CN107907852B (zh) * | 2017-10-27 | 2021-08-03 | 大连大学 | 基于空间平滑的协方差矩阵秩最小化doa估计方法 |
CN107870315B (zh) * | 2017-11-06 | 2021-07-30 | 重庆邮电大学 | 一种利用迭代相位补偿技术估计任意阵列波达方向方法 |
CN109239649B (zh) * | 2018-04-04 | 2023-02-10 | 中国人民解放军空军预警学院 | 一种阵列误差条件下的互质阵列doa估计新方法 |
CN108872929B (zh) * | 2018-04-12 | 2021-03-23 | 浙江大学 | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 |
CN108922553B (zh) * | 2018-07-19 | 2020-10-09 | 苏州思必驰信息科技有限公司 | 用于音箱设备的波达方向估计方法及*** |
CN109061555B (zh) * | 2018-08-27 | 2022-10-11 | 电子科技大学 | 嵌套阵列下混合相干doa估计方法 |
CN109375152B (zh) * | 2018-09-05 | 2020-08-07 | 南京航空航天大学 | 电磁矢量嵌套l阵下低复杂度的doa与极化联合估计方法 |
CN109471086B (zh) * | 2018-10-18 | 2020-11-24 | 浙江大学 | 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法 |
CN109507636B (zh) * | 2018-11-16 | 2022-08-16 | 南京邮电大学 | 基于虚拟域信号重构的波达方向估计方法 |
CN109557503B (zh) * | 2018-12-19 | 2023-03-14 | 成都理工大学 | 基于相关矩阵重建解相干的mimo互质阵列doa估计方法 |
CN109444810B (zh) * | 2018-12-24 | 2022-11-01 | 哈尔滨工程大学 | 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法 |
CN109901101A (zh) * | 2019-02-25 | 2019-06-18 | 西安电子科技大学 | 基于电磁矢量传感器互质阵列相干信号到达角估计方法 |
CN110412535B (zh) * | 2019-08-10 | 2021-08-03 | 浙江大学 | 一种序贯的空时自适应处理参数估计方法 |
JP7044290B2 (ja) | 2020-05-03 | 2022-03-30 | 浙江大学 | 構造化仮想ドメインのテンソル信号の処理に基づく互いに素なエリアアレイの二次元到来方向の推定方法 |
CN111665468B (zh) * | 2020-06-08 | 2022-12-02 | 浙江大学 | 基于单比特量化信号虚拟域统计量重建的互质阵列波达方向估计方法 |
CN111983553B (zh) * | 2020-08-20 | 2024-02-20 | 上海无线电设备研究所 | 一种基于互质多载频稀疏阵列的无网格doa估计方法 |
CN112285642B (zh) * | 2020-09-22 | 2023-09-29 | 华南理工大学 | 一种无重叠优化互质阵列的信号波达方向估计方法 |
CN112305495B (zh) * | 2020-10-22 | 2023-10-13 | 南昌工程学院 | 一种基于原子范数最小的互质阵列协方差矩阵重构方法 |
CN112505675B (zh) * | 2021-02-08 | 2021-06-08 | 网络通信与安全紫金山实验室 | 目标角度和距离定位方法、装置、雷达和存储介质 |
CN113093093B (zh) * | 2021-04-07 | 2023-08-18 | 南京邮电大学 | 基于线性阵列波达方向估计的车辆定位方法 |
CN113820655A (zh) * | 2021-09-18 | 2021-12-21 | 宜宾电子科技大学研究院 | 一种基于Toeplitz矩阵重构和矩阵填充的互质阵相干信号DOA估计方法 |
CN114019446B (zh) * | 2021-10-19 | 2024-04-12 | 南京航空航天大学 | 一种基于去噪核范数最小化的互质相干信源估计方法 |
CN114371440A (zh) * | 2022-01-14 | 2022-04-19 | 天津大学 | 基于信息几何的互质阵doa估计方法 |
CN114444298A (zh) * | 2022-01-21 | 2022-05-06 | 浙江大学 | 基于虚拟域张量填充的互质面阵二维波达方向估计方法 |
CN115236586B (zh) * | 2022-06-30 | 2023-04-18 | 哈尔滨工程大学 | 一种基于数据预处理的极地冰下doa估计方法 |
CN115438604B (zh) * | 2022-11-08 | 2023-03-24 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于质数体系的网格标识方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于互质阵列mimo雷达doa与dod估计方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN106569171A (zh) * | 2016-11-08 | 2017-04-19 | 西安电子科技大学 | 基于双层混合阵列的波达方向角估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10054666B2 (en) * | 2014-09-26 | 2018-08-21 | The United States Of America, As Represented By The Secretary Of The Navy | Sparse space-time adaptive array architecture |
-
2017
- 2017-05-03 CN CN201710302902.4A patent/CN107102291B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于互质阵列mimo雷达doa与dod估计方法 |
CN106324558A (zh) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | 基于互质阵列的宽带信号doa估计方法 |
CN106569171A (zh) * | 2016-11-08 | 2017-04-19 | 西安电子科技大学 | 基于双层混合阵列的波达方向角估计方法 |
Non-Patent Citations (2)
Title |
---|
A Grid-Less Approach to Underdetermined Direction of Arrival Estimation Via Low Rank Matrix Denoising;Piya Pal,et al;《IEEE SIGNAL PROCESSING LETTERS》;20140630;第21卷(第6期);p737-741 |
Coprime Coarray Interpolation for DOA Estimation via Nuclear Norm Minimization;Chun-Lin Liu,et al;《IEEE》;20161231;p2639-2642 |
Also Published As
Publication number | Publication date |
---|---|
CN107102291A (zh) | 2017-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN107315160B (zh) | 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法 | |
CN107329108B (zh) | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 | |
CN107561484B (zh) | 基于内插互质阵列协方差矩阵重建的波达方向估计方法 | |
CN108872929A (zh) | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 | |
CN107329110B (zh) | 基于稀疏阵列直接内插的波达方向估计方法 | |
Guo et al. | Millimeter-wave channel estimation based on 2-D beamspace MUSIC method | |
CN107290709B (zh) | 基于范德蒙分解的互质阵列波达方向估计方法 | |
CN107422295B (zh) | 基于互质阵列虚拟域等价信号原子范数表示的波达方向估计方法 | |
CN107037392B (zh) | 一种基于压缩感知的自由度增加型互质阵列波达方向估计方法 | |
CN107589399B (zh) | 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法 | |
CN106972882B (zh) | 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法 | |
CN108710102B (zh) | 基于互质阵列二阶等价虚拟信号离散傅里叶逆变换的波达方向估计方法 | |
CN107104720B (zh) | 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法 | |
CN111665468B (zh) | 基于单比特量化信号虚拟域统计量重建的互质阵列波达方向估计方法 | |
CN107092004B (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN107015190A (zh) | 基于虚拟阵列协方差矩阵稀疏重建的互质阵列波达方向估计方法 | |
Zheng et al. | Two-dimensional DOA estimation for coprime planar array: A coarray tensor-based solution | |
CN110244273A (zh) | 一种基于均匀分布式阵列的目标角度估计方法 | |
CN108614234B (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
Qi et al. | Time-frequency DOA estimation of chirp signals based on multi-subarray | |
CN108680892B (zh) | 基于角度-空间频率域快速傅里叶变换的互质阵列波达方向估计方法 | |
Tian et al. | 2D-DOA estimation in arc-array with a DNN based covariance matrix completion strategy | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
CN109471087B (zh) | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |