CN113567922A - 一种基于深度展开网络的近场源定位方法及系统 - Google Patents

一种基于深度展开网络的近场源定位方法及系统 Download PDF

Info

Publication number
CN113567922A
CN113567922A CN202110827588.8A CN202110827588A CN113567922A CN 113567922 A CN113567922 A CN 113567922A CN 202110827588 A CN202110827588 A CN 202110827588A CN 113567922 A CN113567922 A CN 113567922A
Authority
CN
China
Prior art keywords
field source
network
arrival
distance
deep
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
CN202110827588.8A
Other languages
English (en)
Other versions
CN113567922B (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202110827588.8A priority Critical patent/CN113567922B/zh
Publication of CN113567922A publication Critical patent/CN113567922A/zh
Application granted granted Critical
Publication of CN113567922B publication Critical patent/CN113567922B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于深度展开网络的近场源定位方法及系统，方法包括：S1、利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；S2、将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；S3、利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；S4、将距离相位向量的信息输入距离深度展开网络，计算近场源的距离。本发明本发明相较于传统神经网络，深度展开网络的参数具有可解释性，具有泛化能力，而相较于ISTA算法，本发明能够提高近场源波达方向的分辨力，加快收敛速度，提高参数估计精度。

Description

一种基于深度展开网络的近场源定位方法及系统
技术领域
本发明涉及阵列信号处理和深度学习技术领域，更具体地说，特别涉及一种基于深度展开网络的近场源定位方法及系统。
背景技术
近场源定位在近炸引信和无源雷达中发挥着重要作用。空间中的辐射源按照与接收阵列的距离可以分为远场源和近场源，远场源可以看做平面波入射，需要对波达方向参数进行估计；近场源需要由球面波进行描述，需要对波达方向和距离参数进行估计。
现有技术“Deep Learning Based Localization of Near-Field Sources withExact Spherical Wavefront Model[J]”(27th European Signal ProcessingConference(EUSIPCO),2019:第1页～第5页)采用卷积神经网络对近场源进行定位，该方法收敛速度快，但是泛化能力较差，即对没有经过训练的近场源参数进行估计时会出现较大的误差。
现有技术“Convolution Neural Networks for Localization of Near-FieldSources via Symmetric Double-Nested Array[J]”(Wireless Communications andMobile Computing，2021，DOI:10.1155/2021/9996780)采用嵌套对称阵列，能够在阵元数相同的情况下能够增大阵列孔径，从而提高近场源参数估计精度，但是该方法为“黑箱模型”，网络的参数不具备可解释性。
现有技术“Direction-of-arrival estimation with circular array usingcompressed sensing in 20GHz band[J]”(IEEEAntennas andWireless PropagationLetters,2021，20(5):第703页～第707页)将迭代收缩阈值算法(Iterative ShrinkageThresholdingAlgorithm，ISTA)展开为网络级联形式，实现对远场源波达方向估计。相较于传统的基于深度学习的数据驱动类方法，该方法将模型驱动的方法融入到网络结构设计中，使得网络的参数具备可解释性，泛化能力得到提高。但是该方法不能够实现近场源定位，并且只能对实数处理，不能对复数信号进行处理。
为此，有必要设计一种基于深度展开网络的近场源定位方法及系统来克服上述技术问题。
发明内容
本发明的目的在于提供一种基于深度展开网络的近场源定位方法及系统，以克服现有技术所存在的缺陷。
为了达到上述目的，本发明采用的技术方案如下：
一种基于深度展开网络的近场源定位方法，包括：
S1、利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；
S2、将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；
S3、利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；
S4、将距离相位向量的信息输入距离深度展开网络，计算近场源的距离。
进一步地，所述步骤S1具体包括：
嵌套对称阵列包括2M+1个阵元，设xm(n)表示第m个阵元的第n个采样点的数值，m＝-M,...,-2,-1,0,1,2,...,M，n＝1,2,...,N，N为快拍数，利用下式计算第m1个阵元、第-m1个阵元、第m2个阵元和第-m2个阵元的四阶累积量c4(m1,m2)：
Figure BDA0003174148910000021
其中(·)*表示取共轭复数运算，m1,m2＝-M,...,-2,-1,0,1,2,…,M，由四阶累积量构成M×M维的近场源波达方向的四阶累积量矩阵C，其中第(m1,m2)个元素为c4(m1,m2)。
进一步地，所述步骤S2具体包括：
首先、将四阶累积量矩阵C的实部和虚部向量化处理，得到向量
Figure BDA0003174148910000022
Figure BDA0003174148910000023
其中
Figure BDA0003174148910000024
(·)T表示向量转置，real(·)表示取实部运算，imag(·)表示取虚部运算，
Figure BDA0003174148910000025
表示四阶累积量矩阵C的第m3列，m3＝-M,...,-2，-1,0,1,2,...,M；
然后、将
Figure BDA0003174148910000026
输入波达方向深度展开网络，得到网络的输出
Figure BDA0003174148910000027
即为近场源波达方向空间谱，Lθ表示网络的层数，对近场源波达方向空间谱
Figure BDA0003174148910000031
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的波达方向估计值
Figure BDA0003174148910000032
k＝1,2,…,K，K表示近场源的数量。
进一步地，所述波达方向深度展开网络的第1层的输出为
Figure BDA0003174148910000033
第lθ层的输出为
Figure BDA0003174148910000034
lθ＝2,3,...,Lθ，hst(·,ε)表示非线性变换函数，即
hst(β,ε)＝sgn(β)⊙max(|β|-ε,0)
sgn(·)表示符号函数，⊙表示Hadamard积，波达方向深度展开网络的初始化参数为：
Figure BDA0003174148910000035
Figure BDA0003174148910000036
ε(θ)＝0.05
其中α(θ)＝0.9/δ(θ)，δ(θ)表示
Figure BDA0003174148910000037
的最大特征值，I(θ)表示单位矩阵，
Figure BDA0003174148910000038
表示波达方向冗余字典矩阵，即
Figure BDA0003174148910000039
Figure BDA00031741489100000314
Figure BDA00031741489100000310
其中，
Figure BDA00031741489100000315
表示Kronecker积，
Figure BDA00031741489100000311
p＝1,2，…，P，P表示波达方向采样点数，λ表示近场源的波长，
Figure BDA00031741489100000312
表示第m3个阵元的位置。
进一步地，在使用波达方向深度展开网络之前需对其进行训练，训练中使用随机梯度下降对网络参数Ψ(θ)、Φ(θ)、ε(θ)进行更新，优化的目标函数为：
Figure BDA00031741489100000313
其中，
Figure BDA0003174148910000041
表示2-范数的平方，||·||1表示1-范数，μ表示正则化参数。
进一步地，所述步骤S3具体包括：
利用公式计算第k个近场源的距离相位向量gk：
Figure BDA0003174148910000042
其中arg(·)表示取相位运算，(·)-1表示求逆运算，
Figure BDA0003174148910000043
Figure BDA0003174148910000044
U表示噪声子空间，
Figure BDA0003174148910000045
可由下式计算
Figure BDA0003174148910000046
进一步地，所述步骤S4具体包括：
首先、由距离相位向量gk的实部和虚部，得到向量
Figure BDA0003174148910000047
Figure BDA0003174148910000048
然后、将向量
Figure BDA0003174148910000049
输入距离深度展开网络，得到网络的输出
Figure BDA00031741489100000410
即为第k个近场源距离空间谱，Lr表示网络的层数；
最后、对近场源距离空间谱
Figure BDA00031741489100000411
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的距离估计值
Figure BDA00031741489100000412
k＝1,2,…,K，K表示近场源的数量。
进一步地，所述距离深度展开网络的第1层的输出为
Figure BDA0003174148910000051
第lr层的输出为
Figure BDA0003174148910000052
其中，lr＝2,3,...,Lr；
网络的初始化参数如下：
Figure BDA0003174148910000053
Figure BDA0003174148910000054
ε(r)＝0.05
其中，α(r)＝0.9/δ(r)，δ(r)表示
Figure BDA0003174148910000055
的最大特征值，I(r)表示单位矩阵，
Figure BDA0003174148910000056
表示距离冗余字典矩阵：
Figure BDA0003174148910000057
B(r)＝[b(r1) b(r2) ... b(rQ)]
其中，
Figure BDA0003174148910000058
Figure BDA0003174148910000059
m4＝0，1，...,M，rq表示空间中距离采样，q＝1，2,…,Q，Q表示距离采样点数。
进一步地，在使用距离深度展开网络之前需对其进行训练，训练中使用SGD对网络参数Ψ(r)、Φ(r)、ε(r)进行更新，优化的目标函数为：
Figure BDA00031741489100000510
本发明还提供一种根据上述的基于深度展开网络的近场源定位方法的系统，包括：
四阶累积量矩阵模块，用于利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；
第一计算模块，将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；
近场源距离相位向量模块，用于利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；
第二计算模块，用于将距离相位向量的信息输入距离深度展开网络，计算近场源的距离；
所述四阶累积量矩阵模块、第一计算模块、近场源距离相位向量模块和第二计算模块依次连接。
与现有技术相比，本发明的优点在于：本发明相较于传统神经网络，深度展开网络的参数具有可解释性，具有泛化能力，而相较于ISTA算法，本发明能够提高近场源波达方向的分辨力，加快收敛速度，提高参数估计精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1本发明提供的基于深度展开网络的近场源定位方法流程图。
图2本发明涉及的波达方向深度展开网络的结构示意图。
图3本发明涉及的距离深度展开网络的结构示意图。
图4利用嵌套对称阵列进行近场源定位的场景示意图。
图5利用本发明得到的近场源在不同波达方向下的估计结果。
图6利用本发明得到的近场源在不同波达方向下的估计误差。
图7利用本发明得到的近场源在不同距离下的估计结果。
图8利用本发明得到的近场源在不同距离下的估计误差。
图9使用不同方法得到的近场源波达方向空间谱。
图10使用不同方法得到的近场源波达方向NMSE。
图11使用不同方法得到的近场源距离NMSE。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。
参阅图1所示，本实施例公开了一种基于深度展开网络的近场源定位方法，包括：
步骤S1、利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵，具体为：
由于已知天线阵列为嵌套对称排布，包括2M+1个阵元。设xm(n)表示第m个阵元的第n个采样点的数值，m＝-M,...,-2,-1，0,1,2,...,M，n＝1,2,...,N，N为快拍数，利用下式计算第m1个阵元、第-m1个阵元、第m2个阵元和第-m2个阵元的四阶累积量c4(m1,m2)：
Figure BDA0003174148910000071
其中，(·)*表示取共轭复数运算，m1,m2＝-M,...,-2,-1,0,1,2,...,M，由四阶累积量构成M×M维的近场源波达方向的四阶累积量矩阵C，其中第(m1,m2)个元素为c4(m1,m2)。
在本实施例中，天线阵列优选嵌套对称线阵，这种阵列能够利用阵元对称的特点得到仅包含近场源波达方向的四阶累积量矩阵，并且在阵元数相同的情况下增大阵列孔径，提高参数估计精度。
步骤S2、将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向，具体为：
首先、将四阶累积量矩阵C的实部和虚部向量化处理，得到向量
Figure BDA0003174148910000072
Figure BDA0003174148910000073
其中，
Figure BDA0003174148910000074
(·)T表示向量转置，real(·)表示取实部运算，imag(·)表示取虚部运算，
Figure BDA0003174148910000075
表示四阶累积量矩阵C的第m3列，m3＝-M,...,-2,-1,0,1,2,...,M。
然后、将向量
Figure BDA0003174148910000076
输入波达方向深度展开网络，得到网络的输出
Figure BDA0003174148910000077
即为近场源波达方向空间谱，Lθ表示网络的层数。
最后、对近场源波达方向空间谱
Figure BDA0003174148910000081
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的波达方向估计值
Figure BDA0003174148910000082
k＝1,2,…,K，K表示近场源的数量。
其中，波达方向深度展开网络的结构如图2所示，其中第1层的输出为：
Figure BDA0003174148910000083
第lθ层的输出为：
Figure BDA0003174148910000084
lθ＝2,3,...,Lθ，hst(·,ε)表示非线性变换函数，即：
hst(β,ε)＝sgn(β)⊙max(|β|-ε,0)
sgn(·)表示符号函数，⊙表示Hadamard积,波达方向深度展开网络的初始化参数如下：
Figure BDA0003174148910000085
Figure BDA0003174148910000086
ε(θ)＝0.05
其中，α(θ)＝0.9/δ(θ)，δ(θ)表示
Figure BDA0003174148910000087
的最大特征值，I(θ)表示单位矩阵，
Figure BDA0003174148910000088
表示波达方向冗余字典矩阵，即：
Figure BDA0003174148910000089
Figure BDA00031741489100000813
Figure BDA00031741489100000810
其中，
Figure BDA00031741489100000814
表示Kronecker积，
Figure BDA00031741489100000811
p＝1,2,…，P，P表示波达方向采样点数，λ表示近场源的波长，
Figure BDA00031741489100000812
表示第m3个阵元的位置。
值得说明的是，在使用波达方向深度展开网络是必须先对它进行训练，训练中使用随机梯度下降(Stochasticgradientdescent，SGD)对网络参数Ψ(θ)、Φ(θ)、ε(θ)进行更新，优化的目标函数为：
Figure BDA0003174148910000091
其中、
Figure BDA0003174148910000092
表示2-范数的平方，||·||1表示1-范数，μ表示正则化参数。
步骤S3、利用波达方向深度展开网络的输出信息，得到近场源距离相位向量，具体为：利用下式计算第k个近场源的距离相位向量gk：
Figure BDA0003174148910000093
其中，arg(·)表示取相位运算，(·)-1表示求逆运算，
Figure BDA0003174148910000094
Figure BDA0003174148910000095
U表示噪声子空间，
Figure BDA0003174148910000096
可由下式计算：
Figure BDA0003174148910000097
步骤S4、将距离相位向量的信息输入距离深度展开网络，计算近场源的距离，具体为：
首先、由距离相位向量gk的实部和虚部，得到向量
Figure BDA0003174148910000098
Figure BDA0003174148910000099
然后、将向量
Figure BDA00031741489100000910
输入距离深度展开网络，得到网络的输出
Figure BDA00031741489100000911
即为第k个近场源距离空间谱，Lr表示网络的层数。
最后、对近场源距离空间谱
Figure BDA00031741489100000912
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的距离估计值
Figure BDA0003174148910000101
k＝1,2,…,K，K表示近场源的数量。
本实施例中的距离深度展开网络的结构如图3所示，其中第1层的输出为：
Figure BDA0003174148910000102
第lr层的输出为：
Figure BDA0003174148910000103
其中lr＝2,3,...,Lr。
距离深度展开网络的初始化参数如下：
Figure BDA0003174148910000104
Figure BDA0003174148910000105
ε(r)＝0.05
其中α(r)＝0.9/δ(r)，δ(r)表示
Figure BDA0003174148910000106
的最大特征值，I(r)表示单位矩阵，
Figure BDA0003174148910000107
表示距离冗余字典矩阵，即
Figure BDA0003174148910000108
B(r)＝[b(r1) b(r2) … b(rQ)]
其中，
Figure BDA0003174148910000109
Figure BDA00031741489100001010
m4＝0，1，...，M，rq表示空间中距离采样，q＝1，2，…，Q，Q表示距离采样点数。
在本实施例中，在使用距离深度展开网络之前必须先对其进行训练，训练中使用SGD对网络参数Ψ(r)、Φ(r)、ε(r)进行更新，优化的目标函数为：
Figure BDA00031741489100001011
本发明还提供一种根据上述的基于深度展开网络的近场源定位方法的系统，包括：四阶累积量矩阵模块，用于利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；第一计算模块，将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；近场源距离相位向量模块，用于利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；第二计算模块，用于将距离相位向量的信息输入距离深度展开网络，计算近场源的距离；其中，四阶累积量矩阵模块、第一计算模块、近场源距离相位向量模块和第二计算模块依次连接。
为验证本发明对近场源的定位性能，通过四个仿真实验进行说明。
仿真实验一
在该仿真实验中，两级嵌套对称线阵如图4所示，阵元个数为17，M＝8，实心圆圈表示第一级子阵，间距为λ/8，空心圆圈表示第二级子阵，间距为5λ/8。以1°和1λ为间隔分别对波达方向[-60°,60°]和距离[5λ,30λ]进行均匀采样产生近场源的训练样本，快拍数N为10000。波达方向深度展开网络的层数Lθ和距离深度展开网络的层数Lr均为50，在训练两个深度展开网络的过程中，epoch均设置为300，mini-batch均设置为32。
仿真实验一用于验证本发明对近场源波达方向估计的泛化能力，实验中近场源的距离设置为10λ，在空域中产生3组测试样本，第1组测试样本波达方向设置为-59.99°，-58.99°，…，-0.99°，0.01°，1.01°，…，59.01°，第2组测试样本波达方向设置为-59.90°，-58.90°，…，-0.90°，0.10°，1.10°，…，59.10°，第3组测试样本波达方向设置为-59.70°，-58.70°，…，-0.70°，0.30°，1.30°，…，59.30°，共计360个测试样本。图5为利用本发明获得的波达方向估计结果，横坐标为测试样本编号，纵坐标为波达方向的估计结果，图6为利用本发明获得的波达方向估计误差，横坐标为测试样本编号，纵坐标为波达方向的估计误差。可以看出第一组和第二组测试样本的波达方向估计误差分别为0.01°和0.1°，第三组测试样本的波达方向估计误差在1°以内，表明本发明对近场源波达方向估计具有泛化能力。
仿真实验二
该仿真实验用于验证本发明对近场源距离估计的泛化能力，实验中近场源的波达方向设置为30°，在空域中产生3组测试样本，第1组测试样本距离设置为5.01λ，6.01λ，…，29.01λ，第2组测试样本距离设置为5.10λ，6.10λ，…，29.10λ，第3组测试样本距离设置为5.30λ，6.30λ，…，29.30λ，共计75个测试样本。图7为利用本发明获得的距离估计结果，横坐标为测试样本编号，纵坐标为距离估计结果，图8为利用本发明获得的距离估计误差，横坐标为测试样本编号，纵坐标为距离估计误差。可以看出第一组和第二组测试样本的距离估计误差为0.01λ和0.1λ，第三组测试样本的距离估计误差在1λ以内，表明本发明对近场源距离估计具有泛化能力。
仿真实验三
该仿真实验用于验证本发明对近场源波达方向的分辨力，实验中两个近场源的位置分别设置为(15°，25λ)、(20°，10λ)。图9表示近场源空间谱，其中横轴表示波达方向，纵轴表示空间谱强度，菱形点表示真实波达方向，实线表示利用本发明得到的近场源的波达方向空间谱，虚线表示利用ISTA算法经过50次迭代得到的近场源的波达方向空间谱，可以看出本发明得到的空间谱出现两个峰值，两个峰值所在位置分别对应近场源的波达方向，而ISTA算法无法对两个近场源的波达方向进行分辨，表明本发明可以提高波达方向的分辨力。
仿真实验四
该仿真实验用于验证本发明收敛速度。图10和图11分别为近场源的波达方向和距离参数的归一化均方根误差(NormalizedMean SquaredError，NMSE)，其中横坐标表示深度展开网络层数或者ISTA算法的迭代次数，纵坐标表示训练样本的NMSE，加“○”的实线表示本发明的NMSE，加“◇”的虚线表示ISTA算法的NMSE。可以看出随着网络层数/迭代次数增加，本发明和ISTA算法NMSE逐渐减小，在相同的网络层数或迭代次数下，本发明具有更快的收敛速度和更高的参数估计精度。
虽然结合附图描述了本发明的实施方式，但是专利所有者可以在所附权利要求的范围之内做出各种变形或修改，只要不超过本发明的权利要求所描述的保护范围，都应当在本发明的保护范围之内。

Claims (10)

1.一种基于深度展开网络的近场源定位方法，其特征在于，包括：
S1、利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；
S2、将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；
S3、利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；
S4、将距离相位向量的信息输入距离深度展开网络，计算近场源的距离。
2.根据权利要求1所述的基于深度展开网络的近场源定位方法，其特征在于，所述步骤S1具体包括：
嵌套对称阵列包括2M+1个阵元，设xm(n)表示第m个阵元的第n个采样点的数值，m＝-M，...，-2，-1，0，1，2，...，M，n＝1，2，...，N，N为快拍数，利用下式计算第m1个阵元、第-m1个阵元、第m2个阵元和第-m2个阵元的四阶累积量c4(m1，m2)：
Figure FDA0003174148900000011
其中(·)*表示取共轭复数运算，m1，m2＝-M，...，-2，-1，0，1，2，...，M，由四阶累积量构成M×M维的近场源波达方向的四阶累积量矩阵C，其中第(m1，m2)个元素为c4(m1，m2)。
3.根据权利要求2所述的基于深度展开网络的近场源定位方法，其特征在于，所述步骤S2具体包括：
首先、将四阶累积量矩阵C的实部和虚部向量化处理，得到向量
Figure FDA0003174148900000012
Figure FDA0003174148900000013
其中
Figure FDA0003174148900000014
(·)T表示向量转置，real(·)表示取实部运算，imag(·)表示取虚部运算，
Figure FDA0003174148900000015
表示四阶累积量矩阵C的第m3列，m3＝-M，…，-2，-1，0，1，2，...，M；
然后、将
Figure FDA0003174148900000021
输入波达方向深度展开网络，得到网络的输出
Figure FDA0003174148900000022
即为近场源波达方向空间谱，Lθ表示网络的层数，对近场源波达方向空间谱
Figure FDA0003174148900000023
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的波达方向估计值
Figure FDA0003174148900000024
K表示近场源的数量。
4.根据权利要求3所述的基于深度展开网络的近场源定位方法，其特征在于，所述波达方向深度展开网络的第1层的输出为
Figure FDA0003174148900000025
第lθ层的输出为
Figure FDA0003174148900000026
lθ＝2，3，...，Lθ，hst(·，ε)表示非线性变换函数，即
hst(β，ε)＝sgn(β)⊙max(|β|-ε，0)
sgn(·)表示符号函数，⊙表示Hadamard积，波达方向深度展开网络的初始化参数为：
Figure FDA0003174148900000027
Figure FDA0003174148900000028
ε(θ)＝0.05
其中α(θ)＝0.9/δ(θ)，δ(θ)表示
Figure FDA0003174148900000029
的最大特征值，I(θ)表示单位矩阵，
Figure FDA00031741489000000210
表示波达方向冗余字典矩阵，即
Figure FDA00031741489000000211
Figure FDA00031741489000000212
Figure FDA00031741489000000213
其中，
Figure FDA00031741489000000214
表示Kronecker积，
Figure FDA00031741489000000215
P表示波达方向采样点数，λ表示近场源的波长，
Figure FDA00031741489000000216
表示第m3个阵元的位置。
5.根据权利要求4所述的基于深度展开网络的近场源定位方法，其特征在于，在使用波达方向深度展开网络之前需对其进行训练，训练中使用随机梯度下降对网络参数Ψ(θ)、Φ(θ)、ε(θ)进行更新，优化的目标函数为：
Figure FDA0003174148900000031
其中，
Figure FDA0003174148900000032
表示2-范数的平方，||·||1表示1-范数，μ表示正则化参数。
6.根据权利要求5所述的基于深度展开网络的近场源定位方法，其特征在于，所述步骤S3具体包括：
利用公式计算第k个近场源的距离相位向量gk：
Figure FDA0003174148900000033
其中arg(·)表示取相位运算，(·)-1表示求逆运算，
Figure FDA0003174148900000034
Figure FDA0003174148900000035
U表示噪声子空间，
Figure FDA0003174148900000036
可由下式计算
Figure FDA0003174148900000037
7.根据权利要求6所述的基于深度展开网络的近场源定位方法，其特征在于，所述步骤S4具体包括：
首先、由距离相位向量gk的实部和虚部，得到向量
Figure FDA0003174148900000038
Figure FDA0003174148900000039
然后、将向量
Figure FDA0003174148900000041
输入距离深度展开网络，得到网络的输出
Figure FDA0003174148900000042
即为第k个近场源距离空间谱，Lr表示网络的层数；
最后、对近场源距离空间谱
Figure FDA0003174148900000043
进行谱峰搜索，由谱峰所对应的位置即可确定近场源的距离估计值
Figure FDA0003174148900000044
K表示近场源的数量。
8.根据权利要求7所述的基于深度展开网络的近场源定位方法，其特征在于，所述距离深度展开网络的第1层的输出为
Figure FDA0003174148900000045
第lr层的输出为
Figure FDA0003174148900000046
其中，lr＝2，3，...，Lr；
网络的初始化参数如下：
Figure FDA0003174148900000047
Figure FDA0003174148900000048
ε(r)＝0.05
其中，α(r)＝0.9/δ(r)，δ(r)表示
Figure FDA0003174148900000049
的最大特征值，I(r)表示单位矩阵，
Figure FDA00031741489000000410
表示距离冗余字典矩阵：
Figure FDA00031741489000000411
B(r)＝[b(r1) b(r2) … b(rQ)]
其中，
Figure FDA00031741489000000412
Figure FDA00031741489000000413
rq表示空间中距离采样，q＝1，2，…，Q，Q表示距离采样点数。
9.根据权利要求8所述的基于深度展开网络的近场源定位方法，其特征在于，在使用距离深度展开网络之前需对其进行训练，训练中使用SGD对网络参数Ψ(r)、Φ(r)、ε(r)进行更新，优化的目标函数为：
Figure FDA0003174148900000051
10.根据权利要求1-9任意一项所述的基于深度展开网络的近场源定位方法的系统，其特征在于，包括：
四阶累积量矩阵模块，用于利用嵌套对称阵列得到近场源波达方向的四阶累积量矩阵；
第一计算模块，将四阶累积量矩阵的信息输入波达方向深度展开网络，计算近场源的波达方向；
近场源距离相位向量模块，用于利用波达方向深度展开网络的输出信息，得到近场源距离相位向量；
第二计算模块，用于将距离相位向量的信息输入距离深度展开网络，计算近场源的距离；
所述四阶累积量矩阵模块、第一计算模块、近场源距离相位向量模块和第二计算模块依次连接。
CN202110827588.8A 2021-07-21 2021-07-21 一种基于深度展开网络的近场源定位方法及系统 Active CN113567922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110827588.8A CN113567922B (zh) 2021-07-21 2021-07-21 一种基于深度展开网络的近场源定位方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110827588.8A CN113567922B (zh) 2021-07-21 2021-07-21 一种基于深度展开网络的近场源定位方法及系统

Publications (2)

Publication Number Publication Date
CN113567922A true CN113567922A (zh) 2021-10-29
CN113567922B CN113567922B (zh) 2023-08-01

Family

ID=78166107

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110827588.8A Active CN113567922B (zh) 2021-07-21 2021-07-21 一种基于深度展开网络的近场源定位方法及系统

Country Status (1)

Country Link
CN (1) CN113567922B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114187A (zh) * 2021-11-18 2022-03-01 中国人民解放军国防科技大学 网格失配情况下基于深度展开admm网络测向方法
CN116054909A (zh) * 2022-12-16 2023-05-02 广东工业大学 一种基于改进型传播算子的卫星通信抗干扰方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160034810A1 (en) * 2014-07-30 2016-02-04 Mitsubishi Electric Research Laboratories, Inc. Neural Networks for Transforming Signals
US20180240219A1 (en) * 2017-02-22 2018-08-23 Siemens Healthcare Gmbh Denoising medical images by learning sparse image representations with a deep unfolding approach
US20180285695A1 (en) * 2017-03-28 2018-10-04 Siemens Healthcare Gmbh Magnetic Resonance Image Reconstruction System and Method
CN111308416A (zh) * 2019-11-29 2020-06-19 北京理工大学 基于四阶累积量的近场非圆信源参数估计方法
CN111965598A (zh) * 2020-08-17 2020-11-20 中国人民解放军国防科技大学 一种利用非均匀阵列进行辐射源定位方法
CN112699735A (zh) * 2020-12-07 2021-04-23 中国人民解放军国防科技大学 基于深度学习的近场源定位方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160034810A1 (en) * 2014-07-30 2016-02-04 Mitsubishi Electric Research Laboratories, Inc. Neural Networks for Transforming Signals
US20180240219A1 (en) * 2017-02-22 2018-08-23 Siemens Healthcare Gmbh Denoising medical images by learning sparse image representations with a deep unfolding approach
US20180285695A1 (en) * 2017-03-28 2018-10-04 Siemens Healthcare Gmbh Magnetic Resonance Image Reconstruction System and Method
CN111308416A (zh) * 2019-11-29 2020-06-19 北京理工大学 基于四阶累积量的近场非圆信源参数估计方法
CN111965598A (zh) * 2020-08-17 2020-11-20 中国人民解放军国防科技大学 一种利用非均匀阵列进行辐射源定位方法
CN112699735A (zh) * 2020-12-07 2021-04-23 中国人民解放军国防科技大学 基于深度学习的近场源定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RONG FU ET AL.: "Compressed LISTA Exploiting Toeplitz Structure", IEEE *
刘振 等: "基于矩阵差分的远场和近场混合源定位方法", 雷达学报 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114187A (zh) * 2021-11-18 2022-03-01 中国人民解放军国防科技大学 网格失配情况下基于深度展开admm网络测向方法
CN114114187B (zh) * 2021-11-18 2022-05-17 中国人民解放军国防科技大学 网格失配情况下基于深度展开admm网络测向方法
CN116054909A (zh) * 2022-12-16 2023-05-02 广东工业大学 一种基于改进型传播算子的卫星通信抗干扰方法
CN116054909B (zh) * 2022-12-16 2023-08-08 广东工业大学 一种基于改进型传播算子的卫星通信抗干扰方法

Also Published As

Publication number Publication date
CN113567922B (zh) 2023-08-01

Similar Documents

Publication Publication Date Title
WO2021139208A1 (zh) 一种基于特定频率组合信号的一维doa估计方法
CN109655799B (zh) 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法
Alinezhad et al. DOA estimation in conformal arrays based on the nested array principles
CN113567922A (zh) 一种基于深度展开网络的近场源定位方法及系统
CN111929637B (zh) 基于互质阵列差和虚拟扩充的一维波达方向估计方法
Wang et al. Sparsity-aware DOA estimation of quasi-stationary signals using nested arrays
CN112666513B (zh) 一种改进式music波达方向估计方法
Xiang et al. Angle separation learning for coherent DOA estimation with deep sparse prior
Li et al. DOA estimation of non-circular source for large uniform linear array with a single snapshot: Extended DFT method
Guo et al. Robust DOA estimation for burst impulsive noise
Pavel et al. Machine learning-based direction-of-arrival estimation exploiting distributed sparse arrays
Chen et al. A robust sparse Bayesian learning-based DOA estimation method with phase calibration
CN113567921B (zh) 一种基于深度展开网络的混合源定位方法及系统
Yu et al. A novel 2D off-grid DOA estimation method based on compressive sensing and least square optimization
Wang et al. Direction finding method via acoustic vector sensor array with fluctuating misorientation
Chen et al. Robust direction-of-arrival estimation via sparse representation and deep residual convolutional network for co-prime arrays
Zhang et al. Two-dimensional underdetermined doa estimation of quasi-stationary signals via sparse bayesian learning
Santori et al. Sensor self-localization for antenna arrays subject to bending and vibrations
Shi et al. A fourth-order cumulant orthonormal propagator rooting method based on Toeplitz approximation
Su et al. Deep unfolding network for near-field source localization via symmetric nested array
Ganti Calibration and direction of arrival performance of sonar arrays composed of multiple sub-arrays
Guo et al. Off-Grid Sparse Bayesian Learning Algorithm for Compressed Sparse Array
Yang et al. Joint calibration of array shape and sensor gain/phase for highly deformed arrays using wideband signals
Shi et al. Sparsity-based DOA estimation of coherent and uncorrelated targets with co-prime MIMO radar
Ganguly et al. Compressive sensing based 2-D DOA estimation by a sparse L-shaped co-prime array

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