CN110161453A - 基于多输入残差网络的互质阵列宽窄带目标doa估计算法 - Google Patents

基于多输入残差网络的互质阵列宽窄带目标doa估计算法 Download PDF

Info

Publication number
CN110161453A
CN110161453A CN201910436028.2A CN201910436028A CN110161453A CN 110161453 A CN110161453 A CN 110161453A CN 201910436028 A CN201910436028 A CN 201910436028A CN 110161453 A CN110161453 A CN 110161453A
Authority
CN
China
Prior art keywords
frequency
array
matrix
correlation matrix
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.)
Granted
Application number
CN201910436028.2A
Other languages
English (en)
Other versions
CN110161453B (zh
Inventor
贾勇
郭勇
钟晓玲
肖钧友
晏超
王刚
宋瑞源
陈胜亿
胡月杨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201910436028.2A priority Critical patent/CN110161453B/zh
Publication of CN110161453A publication Critical patent/CN110161453A/zh
Application granted granted Critical
Publication of CN110161453B publication Critical patent/CN110161453B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Biomedical Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Discrete Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，首先通过两个均匀阵列构造一个互质阵列，确定每个阵列中阵元的位置坐标，获取每个阵元接收到的目标信号；对空间目标信号进行快拍采样，获取目标信号矩阵C，然后进行离散傅里叶变换，得到频谱集G；并处理频谱集G，分别得到单频相关矩阵和多频相关矩阵R；再得到纯实数多频相关矩阵Rr和实数化的纯虚数多频相关矩阵RI；输入多输入残差神经网络；完成DOA估计并输出信号的空间谱。本方法通过残差神经网络有效避免外界干扰的影响从而提高估计的精确度。

Description

基于多输入残差网络的互质阵列宽窄带目标DOA估计算法
技术领域
本发明涉及阵列的DOA估计技术领域，特别是一种基于多输入残差网络的互质阵列宽窄带目标DOA估计算法。
背景技术
波达方向DOA估计可以确定多个空间目标的方位角位置信息，具有高分辨率，广泛应用于通信、雷达、声呐、地震传感等领域。近年提出的互质阵列在阵元位置确定、相邻阵元耦合互扰等方面具有突出的优势，逐渐成为关注的热点。
目前，利用均匀阵列估计目标方位角，最大可分辨的目标数目受限于物理阵元的数目，为突破此限制，考虑到方位估计自由度取决于物理阵元位置差值，最小冗余阵列、(超级)嵌套阵列、互质阵列等非均匀稀疏阵列被引入目标方位参数估计中，利用协同阵的概念，提升自由度，突破物理阵元数目的限制，同时，在同等阵元数目下，非均匀稀疏阵列扩大了阵列孔径，提高了目标方位估计的分辨率。
针对非均匀稀疏阵列方位估计算法的研究，主要集中在传统信号处理类算法，如Capon算法、MUSIC算法、ESPRIT算法、子空间拟合算法、稀疏重构算法等，既打破了阵元数目对最大可分辨目标数目的限制，又提升了目标方位分辨率，但在时效性和鲁棒性(低信噪比、低快拍、大方位角、相干目标)方面存在明显缺陷。此类算法或多或少存在计算复杂时效性不足、低信噪比低快拍估计精度差、大方位角目标和相干目标估计性能下降、环境适应能力不足等问题。
基于互质接收阵列的非相干目标DOA估计方法，构建适用于非均匀互质阵列的深度学习方位识别网络，在保留传统信号处理类算法具有的方位分辨率高、最大可分辨目标数目突破物理阵元数目限制的优势下，进一步增强非均匀互质阵列对目标方位估计的时效性和环境适应性，以及在低信噪比、低快拍、大方位角、相干目标条件下的鲁棒性。由于考虑到经典的深度神经网络在训练过程中可能会出现梯度爆炸和消失的缘故，本发明选择了残差神经网络。
该算法主要解决面向非均匀互质接收阵列的、适用于窄带/宽带和多个相干/非相干信号源的特征表达域的形成问题和对应的深度学习网络的构建问题。其实质是通过将时域多快拍数据转换到频域多频率单快拍数据再针对各个频率进行相关计算，实现对窄带和宽带的兼容；同时利用“差协同阵”的概念，对不同频率的相关矩阵进行矢量化处理，相关矢量合并形成多频相关矩阵特征表达域，既解除了相干源的影响，实现了相干与非相干信号源兼容，又扩充了特征表达维度，提升了多信号源DOA估计能力。对此多频相关矩阵特征表达域，针对性设计残差神经网络从复杂的多频相关矩阵提取与信号源DOA相关的空域特征映射形成空间谱。
发明内容
有鉴于此，本发明的目的在于提供一种基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，该方法结合了阵列DOA估计与深度学习。
为达到上述目的，本发明提供如下技术方案：
本发明提供的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，包括以下步骤：
(1)使用两个均匀阵列构造一个互质阵列，确定每个阵列中阵元的位置坐标，获取每个阵元接收到的目标信号；
(2)目标向整个空域辐射信号波，对空间目标信号进行快拍采样，获取目标信号矩阵C，所述目标信号矩阵C的每一排为同一阵元接收的信号，每一列为同一快拍下整个阵列接收的信号；
(3)将目标信号矩阵C中每一排数据进行离散傅里叶变换，得到频谱集G；
(4)将频谱集G按频率划分为N个窄带信号，对每个窄带进行单频率相关处理得到单频相关矩阵
(5)将单频相关矩阵中的元素一一取出并按照差值大小进行重新排序，生成单频率相关向量：将宽带中所有的相关向量按频率大小从低到高进行组合，生成多频相关矩阵R；
(7)将多频相关矩阵R分解为两个对应的纯实数多频相关矩阵Rr和纯虚数多频相关矩阵Ri，提取目标信号表征域矩阵，得到纯实数多频相关矩阵Rr和实数化的纯虚数多频相关矩阵RI，至此完成阵列信号特征表达域处理与提取过程；
(8)建立提取表征域矩阵中角度信息的残差神经网络，用于鉴别目标波达方向，输入网络的数据为纯实数多频相关矩阵Rr和实数化的纯虚数多频相关矩阵RI，构造多输入残差神经网络；完成DOA估计并输出信号的空间谱。
进一步，所述互质矩阵是使用两个均匀阵列按照以下方式进行构造的：
均匀阵列P的阵元坐标集合为A1＝{A1 i|i＝0,Qd,...(P-1)Qd}；
均匀矩阵Q的阵元坐标集合为A2＝{A2 i|i＝Pd,2Pd,...(2Q-1)Pd}；
其中，P和Q为互质关系，Q<P，λ为波长，对应信号频率为fmax且为宽带信号中最大频率，基础间距为 表示均匀阵列P的第i个阵元的坐标；表示均匀阵列Q的第i个阵元的坐标；i表示阵元的序数；p表示均匀阵列P的阵元数目；Q表示均匀阵列Q的阵元数目；
其中阵元接收到的目标信号按照以下公式表示：
xm(t)＝sm(t)e-i2πd(m-1)sin(θ+nm(t)；
其中，xm(t)为坐标m的阵元接收的信号，θ为信号源入射角度，nm(t)为坐标m的阵元接收的高斯白噪声，Sm(t)表示信号源的信号。
进一步，所述目标信号矩阵C具体如下表示：
其中，每一排为同一阵元接收的信号，每一列为同一快拍下整个阵列接收的信号；信号波的辐射频率为0-fmax；其中，fmax表示最大辐射频率；K表示对空间目标信号进行快拍采样的次数。
进一步，所述频谱集G的表示具体如下：
其中，为坐标i的阵元接收的第n个频点的频谱，N表示离散傅里叶变换点数。
进一步，所述单频相关矩阵是按照按照以下步骤进行处理的：
将频谱集G按频率划分为N个窄带信号，每个窄带表示为：
Xn＝gn；其中，gn表示所有阵元离散傅里叶变换后第n个频点的频谱；n表示第n个频点；
对每个窄带Xn按照如下公式进行单频率相关处理；
其中，xn(0)表示坐标为原点的阵元接到到的空间中的信号，xn(Q)表示坐标为Q的阵元接到到的空间中的信号，rn(0)表示表示坐标为原点的阵元接收到的信息，rn(Q)表示表示坐标为Q的阵元接收到的信息，H表示矩阵的共轭转置，右上角的*表示共轭操作，单频相关矩阵中每个元素由两个物理阵元决定，对两个物理阵元求差值得到以单位间距d为增加量的阵元信息，从(1-2Q)Pd延伸到(2Q-1)Pd的2P(2Q-1)+1个差值。
进一步，所述多频相关矩阵R按照以下方式生成：
将单频相关矩阵中的元素一一取出并按照差值大小进行重新排序，生成频点为n的单频相关向量：
r＝[rn((1-2Q)P … rn(-1) rn(0) rn(1) … rn((2Q-1)P)]；
其中，r表示相关向量；n表示频点n；
将宽带中所有的相关向量按频率大小从低到高进行组合，生成多频相关矩阵R；所述多频相关矩阵R的表示如下：
其中，r1表示频点为1的单频相关矩阵矢量化后的单频相关向量；rN表示频点为N的单频相关矩阵矢量化后的单频相关向量；1和N分别对应前面的N点DFT的第1点和第N点。
进一步，所述实数化的纯虚数多频相关矩阵RI是按照以下步骤形成：
将多频相关矩阵R分解为两个对应的纯实数矩阵Rr和纯虚数矩阵Ri，所述纯实数矩阵和纯虚数矩阵的大小与多频相关矩阵R相同；
所述纯实数矩阵Rr如下所示：
其中，表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯实数部分；rr N表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯实数部分；
所述纯虚数矩阵Ri如下所示：
其中，表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；
按照以下公式计算得到实数化的纯虚数多频相关矩阵RI：
RI＝Ri.*(-i)；
其中，纯虚数矩阵Ri与-i之间采用点乘计算，得到与纯虚数矩阵Ri大小相同的实数化的纯虚数多频相关矩阵RI。
进一步，所述多输入残差神经网络具体结构为：
并行网络第一层是卷积神经网络；二、三并行层为两个残差块，之后在第三个残差块汇合，合并后的网络有6个残差块，再连接一层平均池化层，网络在最后连接一层有1800个神经元的全连接层，完成基于角度的分类问题，最终完成DOA估计并输出信号的空间谱，分辨力为0.1°。
进一步，还包括生成对应角度标签向量T，所述对应角度标签向量为后期神经网络训练提供标签，配置向量过程为：对应角度标签向量中在相应角度处置1，其余角度处置0。
本发明的有益效果在于：
本发明提出了一种对多目标宽窄带信号源DOA估计的方法，包括对目标数量的估计，对多目标方位角的估计。该方法对传统DOA估计中数据处理部分进行了改进，传统的信号处理算法存在计算复杂时效性不足、低信噪比低快拍估计精度差、大方位角目标和相干目标估计性能下降、环境适应能力不足等问题。为了改进传统方法的问题，首先对阵列采集的目标信号数据进行表征域提取，从原始数据中提前提取有关于目标角度信息的数据，可以为后期训练减小负担。再送入深度神经网络进行训练，让网络学习到有关目标角度的特征并构建网络参数，最后送入数据得到多目标的空间谱。通过残差神经网络可以有效的避免外界干扰的影响从而提高估计的精确度。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
图1为阵列结构。
图2为多输入残差神经网络结构。
图3为多频特征表达域的提取流程图。
图4为整体流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
实施例1
如图1所示，图1为阵列结构，图中黑色实心三角代表接收阵元；本实施例提供的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，包括以下步骤：
(1)使用两个均匀阵列构造一个互质阵列，
一个子阵的阵元坐标集合为A1＝{A1 i|i＝0,Qd,…(P-1)Qd}；
另一个子阵的阵元坐标集合为A2＝{A2 i|i＝Pd,2Pd,…(2Q-1)Pd}；
其中，P和Q为互质关系，Q<P，λ为波长，对应信号频率为fmax且为宽带信号中最大频率，基础间距为
阵元接收到的目标信号按照以下公式表示：
xm(t)＝sm(t)e-i2πd(m-1)sin(θ+nm(t)；
其中，xm(t)为坐标m的阵元接收的信号，θ为信号源入射角度，nm(t)为坐标m的阵元接收的高斯白噪声。
(2)目标向整个空域辐射频率0-fmax的信号波，对空间目标信号进行K次快拍采样，获得(2Q+P-1)×K的目标信号矩阵C，所述目标信号矩阵C采用如下矩阵表示：
每一排为同一阵元接收的信号，每一列为同一快拍下整个阵列接收的信号；
同时生成对应角度标签向量T，所述对应角度标签向量大小为1×1800，为后期神经网络训练提供标签，配置向量过程为：对应角度标签向量中在相应角度处置1，其余角度处置0；
(3)将目标信号矩阵C中每一排数据进行N点的离散傅里叶变换，得到如下频谱集：
其中，为坐标i的阵元接收的第n个频点的频谱。
(4)将频谱集G按频率划分为N个窄带信号，每个窄带表示为：
Xn＝gn；
对每个窄带Xn按照如下公式进行单频率相关处理；
其中，xn(0)表示坐标为原点的阵元接到到的空间中的信号，xn(Q)表示坐标为Q的阵元接到到的空间中的信号，rn(0)表示表示坐标为原点的阵元接收到的信息，rn(Q)表示表示坐标为Q的阵元接收到的信息，其中会有物理阵元没有的阵元，称之为虚拟阵元；H表示矩阵的共轭转置，右上角的*表示共轭操作，单频相关矩阵中每个元素由两个物理阵元决定，对两个物理阵元求差值得到以单位间距d为增加量的阵元信息，从(1-2Q)Pd延伸到(2Q-1)Pd的2P(2Q-1)+1个差值。
(5)将中的元素一一取出并按照差值大小进行重新排序，生成频点为n的相关向量：
r＝[rn((1-2Q)P … rn(-1) rn(0) rn(1) … rn((2Q-1)P)]
r表示相关向量；
将宽带中所有的相关向量按频率大小从低到高进行组合，生成矩阵R；
其中，r1表示频点为1的单频相关矩阵矢量化后的单频相关向量；rN表示频点为N的单频相关矩阵矢量化后的单频相关向量；1和N分别对应前面的N点DFT的第1点和第N点。
(7)将矩阵R分解为两个对应的纯实数矩阵和纯虚数矩阵，矩阵大小与R相同；
所述纯实数矩阵Rr如下所示：
其中，表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯实数部分；表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯实数部分；
所述纯虚数矩阵Ri如下所示：
其中，表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；
将Ri与-i进行点乘，得到与Ri矩阵大小相同的矩阵RI；
RI＝Ri.*(-i)；
至此完成目标信号表征域提取，得到矩阵Rr和RI。至此完成阵列信号特征表达域处理与提取过程，如图3所示，图3为多频特征表达域的提取流程图。
(8)建立提取表征域矩阵中角度信息的残差神经网络，用于鉴别目标波达方向，输入网络的数据为Rr和RI，构造多输入残差神经网络，具体如下：
如图2所示，图2为多输入残差神经网络结构图，其中，3×3代表卷积核的大小，64、128、256、512分别代表对应残差块中卷积核数量，Conv为卷积层。/2是对卷积核的数量调整，以匹配该残差块输出，Avg pool为平均池化层，FC1800层是一层有1800个神经元的全连接层。
本方法中使用的是多输入残差网络，具体结构为：并行网络第一层是卷积神经网络；二、三并行层为两个残差块，之后在第三个残差块汇合，合并后的网络有6个残差块，再连接一层平均池化层，网络在最后连接一层有1800个神经元的全连接层，完成基于角度的分类问题，最终完成DOA估计并输出信号的空间谱，分辨力为0.1°。
整个网络中并行的第一个残差块的步长为1，其余的步长均为2，激活函数均采用ReLu激活函数。
(9)本方法采用的数据集通过matlab仿真生成，数据集具体构成如下：
1、单信号源：目标入射角度从0°到180°均匀的1800个角度遍历生成仿真数据，目标频率随机分布在fmax以内，并生成对应角度标签一起以txt文件形式存放。一共遍历15次，单一信号源数据集中包含27000个数据。
2、多信号源：通过matlab随机生成目标数小于20的多目标信号，目标入射角度随机分布在从0°到180°均匀的1800个角度，信号频率随机分布在fmax以内，并生成对应角度标签一起以txt文件形式存放。多信号源数据集包含40000个数据。
将上述数据集中单信号源数据和多信号源数据各自按照7：3的比例随机分出训练集和测试集。将各自的训练集和测试集合并在一起，就准备好了训练集和测试集。
将训练集数据与对应的角度标签T送入深度神经网络中完成训练。输出层的角度判决准确率达到97.6％，超过预设值96％，最后得到训练完成的深度神经网络。
将测试集数据送入深度神经网络中进行测试，实验表明训练好的深度神经网络准确度能够达到98.1％,实时性和抗干扰能力很好，达到本发明的预设效果。
(10)本实施例整体流程图如图4所示。其中，多频特征表达域的整体提取流程图如图4所示；本实施例使用matlab，Python3.7.1，Pycharm以及pytorch作为实验平台，通过编程实现本方法，实验证明，使用新的方法进行DOA估计，可以很好的进行目标数目的估计和空间谱的生成，其中，互质阵列可以弥补实际情况中物理阵元数量不足的问题，以提高阵列同时鉴别更多目标的能力，深度神经网络可以从大量的训练集数据中学习到角度信息的特征，并且从大量的数据中去除干扰因素，比如背景的噪声信号、不同的温度、湿度以及阵元微小的间距误差等等。传统的方法无法考虑到非常多的因素所以通过深度学习方法就可以规避误差，这是一个优势。当在同时估计大量目标的时候就需要很深的神经网络，经典的网络可能会出现梯度爆炸和梯度消失的情况，残差神经网络中残差块的首尾之间有一条传递捷径，可以增强梯度的传导能力，就能避免以上情况的发生，所以残差神经网络是一个很好的选择。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (9)

1.基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：包括以下步骤：
(1)使用两个均匀阵列构造一个互质阵列，确定每个阵列中阵元的位置坐标，获取每个阵元接收到的目标信号；
(2)目标向整个空域辐射信号波，对空间目标信号进行快拍采样，获取目标信号矩阵C，所述目标信号矩阵C的每一排为同一阵元接收的信号，每一列为同一快拍下整个阵列接收的信号；
(3)将目标信号矩阵C中每一排数据进行离散傅里叶变换，得到频谱集G；
(4)将频谱集G按频率划分为N个窄带信号，对每个窄带进行单频率相关处理得到单频相关矩阵
(5)将单频相关矩阵中的元素一一取出并按照差值大小进行重新排序，生成单频率相关向量：将宽带中所有的相关向量按频率大小从低到高进行组合，生成多频相关矩阵R；
(7)将多频相关矩阵R分解为两个对应的纯实数多频相关矩阵Rr和纯虚数多频相关矩阵Ri，提取目标信号表征域矩阵，得到纯实数多频相关矩阵Rr和实数化的纯虚数多频相关矩阵RI，至此完成阵列信号特征表达域处理与提取过程；
(8)建立提取表征域矩阵中角度信息的残差神经网络，用于鉴别目标波达方向，输入网络的数据为纯实数多频相关矩阵Rr和实数化的纯虚数多频相关矩阵RI，构造多输入残差神经网络；完成DOA估计并输出信号的空间谱。
2.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述互质矩阵是使用两个均匀阵列按照以下方式进行构造的：
均匀阵列P的阵元坐标集合为A1＝{A1 i|i＝0,Qd,…(P-1)Qd}；
均匀矩阵Q的阵元坐标集合为A2＝{A2 i|i＝Pd,2Pd,…(2Q-1)Pd}；
其中，P和Q为互质关系，Q<P，λ为波长，对应信号频率为fmax且为宽带信号中最大频率，基础间距为 表示均匀阵列P的第i个阵元的坐标；表示均匀阵列Q的第i个阵元的坐标；i表示阵元的序数；p表示均匀阵列P的阵元数目；Q表示均匀阵列Q的阵元数目；
其中阵元接收到的目标信号按照以下公式表示：
xm(t)＝sm(t)e-i2πd(m-1)sin(θ)+nm(t)；
其中，xm(t)为坐标m的阵元接收的信号，θ为信号源入射角度，nm(t)为坐标m的阵元接收的高斯白噪声，Sm(t)表示信号源的信号。
3.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述目标信号矩阵C具体如下表示：
其中，每一排为同一阵元接收的信号，每一列为同一快拍下整个阵列接收的信号；信号波的辐射频率为0-fmax；其中，fmax表示最大辐射频率；K表示对空间目标信号进行快拍采样的次数。
4.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述频谱集G的表示具体如下：
其中，为坐标i的阵元接收的第n个频点的频谱，N表示离散傅里叶变换点数。
5.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述单频相关矩阵是按照按照以下步骤进行处理的：
将频谱集G按频率划分为N个窄带信号，每个窄带表示为：
Xn＝gn；其中，gn表示所有阵元离散傅里叶变换后第n个频点的频谱；n表示第n个频点；
对每个窄带Xn按照如下公式进行单频率相关处理；
其中，xn(0)表示坐标为原点的阵元接到到的空间中的信号，xn(Q)表示坐标为Q的阵元接到到的空间中的信号，rn(0)表示表示坐标为原点的阵元接收到的信息，rn(Q)表示表示坐标为Q的阵元接收到的信息，H表示矩阵的共轭转置，右上角的*表示共轭操作，单频相关矩阵中每个元素由两个物理阵元决定，对两个物理阵元求差值得到以单位间距d为增加量的阵元信息，从(1-2Q)Pd延伸到(2Q-1)Pd的2P(2Q-1)+1个差值。
6.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述多频相关矩阵R按照以下方式生成：
将单频相关矩阵中的元素一一取出并按照差值大小进行重新排序，生成频点为n的单频相关向量：
r＝[rn((1-2Q)P…rn(-1) rn(0) rn(1)…rn((2Q-1)P)]；
其中，r表示相关向量；n表示第n个频点；
将宽带中所有的相关向量按频率大小从低到高进行组合，生成多频相关矩阵R；所述多频相关矩阵R的表示如下：
其中，r1表示频点为1的单频相关矩阵矢量化后的单频相关向量；rN表示频点为N的单频相关矩阵矢量化后的单频相关向量；1和N分别对应前面的N点DFT的第1点和第N点。
7.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述实数化的纯虚数多频相关矩阵RI是按照以下步骤形成：
将多频相关矩阵R分解为两个对应的纯实数矩阵Rr和纯虚数矩阵Ri，所述纯实数矩阵和纯虚数矩阵的大小与多频相关矩阵R相同；
所述纯实数矩阵Rr如下所示：
其中，rr 1表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯实数部分；rr N表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯实数部分；
所述纯虚数矩阵Ri如下所示：
其中，ri 1表示频点为1的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；ri N表示频点为N的单频相关矩阵矢量化后的单频相关向量的纯虚数部分；
按照以下公式计算得到实数化的纯虚数多频相关矩阵RI：
RI＝Ri.*(-i)；
其中，纯虚数矩阵Ri与-i之间采用点乘计算，得到与纯虚数矩阵Ri大小相同的实数化的纯虚数多频相关矩阵RI。
8.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：所述多输入残差神经网络具体结构为：
并行网络第一层是卷积神经网络；二、三并行层为两个残差块，之后在第三个残差块汇合，合并后的网络有6个残差块，再连接一层平均池化层，网络在最后连接一层有1800个神经元的全连接层，完成基于角度的分类问题，最终完成DOA估计并输出信号的空间谱，分辨力为0.1°。
9.如权利要求1所述的基于多输入残差网络的互质阵列宽窄带目标DOA估计算法，其特征在于：还包括生成对应角度标签向量T，所述对应角度标签向量为后期神经网络训练提供标签，配置向量过程为：对应角度标签向量中在相应角度处置1，其余角度处置0。
CN201910436028.2A 2019-05-23 2019-05-23 基于多输入残差网络的互质阵列宽窄带目标doa估计算法 Expired - Fee Related CN110161453B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910436028.2A CN110161453B (zh) 2019-05-23 2019-05-23 基于多输入残差网络的互质阵列宽窄带目标doa估计算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910436028.2A CN110161453B (zh) 2019-05-23 2019-05-23 基于多输入残差网络的互质阵列宽窄带目标doa估计算法

Publications (2)

Publication Number Publication Date
CN110161453A true CN110161453A (zh) 2019-08-23
CN110161453B CN110161453B (zh) 2021-07-02

Family

ID=67632468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910436028.2A Expired - Fee Related CN110161453B (zh) 2019-05-23 2019-05-23 基于多输入残差网络的互质阵列宽窄带目标doa估计算法

Country Status (1)

Country Link
CN (1) CN110161453B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110554352A (zh) * 2019-09-11 2019-12-10 哈尔滨工业大学 基于vgg16神经网络的航天测控系统干扰源波达方向估计方法
CN111142066A (zh) * 2019-12-25 2020-05-12 达闼科技成都有限公司 波达方向估计方法、服务器以及计算机可读存储介质
CN111983553A (zh) * 2020-08-20 2020-11-24 上海无线电设备研究所 一种基于互质多载频稀疏阵列的无网格doa估计方法
US20210150286A1 (en) * 2019-11-20 2021-05-20 Rohde & Schwarz Gmbh & Co. Kg Method and system for detecting and/or classifying a wanted signal
CN115792849A (zh) * 2022-11-23 2023-03-14 哈尔滨工程大学 一种基于sac算法的一维非均匀阵列设计方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324558A (zh) * 2016-08-30 2017-01-11 东北大学秦皇岛分校 基于互质阵列的宽带信号doa估计方法
CN106646344A (zh) * 2016-12-16 2017-05-10 西北工业大学 一种利用互质阵的波达方向估计方法
CN106896340A (zh) * 2017-01-20 2017-06-27 浙江大学 一种基于压缩感知的互质阵列高精度波达方向估计方法
CN107545245A (zh) * 2017-08-14 2018-01-05 中国科学院半导体研究所 一种年龄估计方法及设备
CN107589399A (zh) * 2017-08-24 2018-01-16 浙江大学 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法
CN109085531A (zh) * 2018-08-27 2018-12-25 西安电子科技大学 基于神经网络的近场源到达角估计方法
CN109255308A (zh) * 2018-11-02 2019-01-22 陕西理工大学 存在阵列误差的神经网络到达角估计方法
CN109471086A (zh) * 2018-10-18 2019-03-15 浙江大学 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324558A (zh) * 2016-08-30 2017-01-11 东北大学秦皇岛分校 基于互质阵列的宽带信号doa估计方法
CN106646344A (zh) * 2016-12-16 2017-05-10 西北工业大学 一种利用互质阵的波达方向估计方法
CN106896340A (zh) * 2017-01-20 2017-06-27 浙江大学 一种基于压缩感知的互质阵列高精度波达方向估计方法
CN107545245A (zh) * 2017-08-14 2018-01-05 中国科学院半导体研究所 一种年龄估计方法及设备
CN107589399A (zh) * 2017-08-24 2018-01-16 浙江大学 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法
CN109085531A (zh) * 2018-08-27 2018-12-25 西安电子科技大学 基于神经网络的近场源到达角估计方法
CN109471086A (zh) * 2018-10-18 2019-03-15 浙江大学 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法
CN109255308A (zh) * 2018-11-02 2019-01-22 陕西理工大学 存在阵列误差的神经网络到达角估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王华奎等: "基于互质阵列的宽带信号DOA估计", 《声学技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110554352A (zh) * 2019-09-11 2019-12-10 哈尔滨工业大学 基于vgg16神经网络的航天测控系统干扰源波达方向估计方法
US20210150286A1 (en) * 2019-11-20 2021-05-20 Rohde & Schwarz Gmbh & Co. Kg Method and system for detecting and/or classifying a wanted signal
CN111142066A (zh) * 2019-12-25 2020-05-12 达闼科技成都有限公司 波达方向估计方法、服务器以及计算机可读存储介质
CN111983553A (zh) * 2020-08-20 2020-11-24 上海无线电设备研究所 一种基于互质多载频稀疏阵列的无网格doa估计方法
CN111983553B (zh) * 2020-08-20 2024-02-20 上海无线电设备研究所 一种基于互质多载频稀疏阵列的无网格doa估计方法
CN115792849A (zh) * 2022-11-23 2023-03-14 哈尔滨工程大学 一种基于sac算法的一维非均匀阵列设计方法及系统

Also Published As

Publication number Publication date
CN110161453B (zh) 2021-07-02

Similar Documents

Publication Publication Date Title
CN110161453A (zh) 基于多输入残差网络的互质阵列宽窄带目标doa估计算法
Shen et al. Wind-speed inversion from HF radar first-order backscatter signal
CN110133576A (zh) 基于级联残差网络的双基互质mimo阵列方位估计算法
CN110275166A (zh) 基于admm的快速稀疏孔径isar自聚焦与成像方法
CN109471082A (zh) 基于信号子空间重构的阵元缺损mimo雷达角度估计方法
CN105699945A (zh) 频控阵mimo雷达系统中的波形优化设计方法
CN108549059A (zh) 一种复杂地形条件下的低空目标仰角估计方法
CN106483516A (zh) 基于先验知识的雷达杂波空时自适应处理方法
CN107167783A (zh) 一种共形阵列杂波协方差矩阵的稀疏重构方法
CN106772273B (zh) 一种基于动态孔径的sar虚假目标干扰抑制方法及系统
CN111950438B (zh) 基于深度学习的天宫二号成像高度计有效波高反演方法
CN109597046A (zh) 基于一维卷积神经网络的米波雷达doa估计方法
CN107576961A (zh) 一种互质降采样间歇合成孔径雷达稀疏成像方法
CN106680815A (zh) 基于张量稀疏表示的mimo雷达成像方法
CN106526563B (zh) 一种基于互相关虚拟阵的五元体积阵多目标方位估计方法
CN107783078A (zh) 一种波束‑多普勒酉esprit多目标角度估计方法
CN105388462A (zh) 一种基于互质阵列的虚拟波束形成方法
Li et al. Extraction of wind direction spreading factor from broad-beam high-frequency surface wave radar data
CN109283497A (zh) 双基地fda-mimo雷达距离欺骗干扰识别方法
CN113868583A (zh) 一种子阵波束聚焦的声源距离计算方法及系统
CN101470195B (zh) 一种非参数型高分辨波束形成方法与装置
CN110346752A (zh) 基于互质稀疏阵的无模糊测向方法
CN109298382A (zh) 一种基于期望极大算法的非均匀直线阵波达方向角估计方法
CN109100711A (zh) 一种深海环境下单基地主动声纳低运算量三维定位方法
CN108398659A (zh) 一种矩阵束与求根music结合的波达方向估计方法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210702