CN113514808A - 一种用于判定小型无人机目标个数的智能辨识方法 - Google Patents

一种用于判定小型无人机目标个数的智能辨识方法 Download PDF

Info

Publication number
CN113514808A
CN113514808A CN202110400484.9A CN202110400484A CN113514808A CN 113514808 A CN113514808 A CN 113514808A CN 202110400484 A CN202110400484 A CN 202110400484A CN 113514808 A CN113514808 A CN 113514808A
Authority
CN
China
Prior art keywords
circle
matrix
unmanned aerial
neural network
channel
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
CN202110400484.9A
Other languages
English (en)
Other versions
CN113514808B (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.)
Tong Fang Electronic Science & Technology Co ltd
Civil Aviation Flight University of China
Original Assignee
Tong Fang Electronic Science & Technology Co ltd
Civil Aviation Flight University of China
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 Tong Fang Electronic Science & Technology Co ltd, Civil Aviation Flight University of China filed Critical Tong Fang Electronic Science & Technology Co ltd
Priority to CN202110400484.9A priority Critical patent/CN113514808B/zh
Publication of CN113514808A publication Critical patent/CN113514808A/zh
Application granted granted Critical
Publication of CN113514808B publication Critical patent/CN113514808B/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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411—Identification of targets based on measurements of radar reflectivity
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

一种用于判定小型无人机目标个数的智能辨识方法，其特征在于，包括以下步骤：S1：将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号；S2：对各个通道接收信号的数据集进行通道校正，获得各通道校正后观测信号数据集；S3：各通道校正后观测信号数据集求协方差矩阵，再经酉变换计算得到盖尔圆矩阵；S4：通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络；S5：改变预先设定的无人机个数，多次重复步骤S1～S4，生成一定数量的训练样本用于训练步骤S4中所述的分类神经网络。

Description

一种用于判定小型无人机目标个数的智能辨识方法
技术领域
本发明属于小型无人机目标智能识别与控制技术领域，具体涉及一种用于判定小型无人机目标个数的智能辨识方法。
背景技术
近年来，随着小型无人机应用技术迅猛发展，小型无人机在农业、电力、测绘、航拍、安全等领域发挥了重大作用，小型无人机应用模式也已从单机应用向多无人机系统方向发展。作为小型无人机识别与控制技术的重要技术组成部分，无人机识别与控制技术也由单无人机识别向多无人机识别与控制技术方向发展，在无人机目标识别与控制任务处理过程中，能否准确识别无人机目标数目直接关系到能否准确锁定无人机目标方位。如果对无人机目标数目判别结果与实际目标数目不符，则会直接影响到后续无人机目标方位判定，进而无法为无人机后续控制提供关键参数。因此无人机目标智能识别对实现无人机目标协议控制任务、制定无人机控制链路信号定向辐射方位策略、以及评估无人机目标所在空域的飞行安全风险等具有十分重要的意义。
在技术层面，无人机目标个数识别属于阵列信号处理领域中的信源个数估计问题。当前，基于天线阵列输出数据进行无人机目标个数判别的方法总体上分为三类。一类是基于对天线阵列输出快拍数据的协方差矩阵进行特征分解，根据特征分解后的特征值大小判定无人机目标个数，但是这类判别方法至今为止缺乏理论依据。另一类则是基于信息论准则的判别方法，主要是基于H.Akaike提出的最小信息准则(AIC)(H Akaike.A new lookat the statistical model identification[J].IEEE Transactions on AutomaticControl,1974,19(6):716-723.)和基于M.Wax提出的最小描述长度准则(MDL)(Wax M,Ziskind I.Detection ofthe number ofcoherent signals by the MDL principle[J].IEEE Transactions on Acoustics Speech and Signal Processing,1989,37(8):1190-1196.)，以及在二者基础上提出的各种改进方法。然而，基于AIC准则和基于MDL准则的各种方法均存在仅适用于白噪声场景下目标个数估计的技术缺陷，对于色噪声环境下的目标个数估计准确性具有不确定性。第三类目标个数判定方法主要针对色噪声场景下提出的，主要思想是基于盖尔圆盘准则(GDE)的目标个数判定方法,虽然能够弥补信息论准则在色噪声下信源个数估计失效的不足，但是GDE准则在低信噪比，小快拍数条件下的性能会急剧下降甚至失效。为改善GDE准则存在的不足，2010年，Yin-zhen提出了一种新的酉变换方法，并结合GDE准则进行目标个数估计。实验表明，该方法在快拍数为90时的信源个数检测准确率达到90％，但是要求的信噪比在10dB以上(Yin-zhen Mei,Peng Wang.Estimation ofsource number using Gerschgorin disks based on acoustic vector sensor array[C]//2010International Conference on ComputerApplication and System Modeling(ICCASM 2010),Taiyuan,2010:V7-291-V7-295.)。2011年，Z.-Liu等人利用接收信号协方差矩阵的对角线平均值来构建新的协方差矩阵，并结合GDE准则估计信号源数，该算法在非平稳的色噪声中信噪比为-8dB﹑快拍数为100时的检测准确率就已经达到80％以上，该方法虽能同时保证在低信噪比小快拍数下检测准确率高，但仅仅是在10个阵元估计2个信源的条件下(Z-Liu,Z-Lu,Z-Huang,et al.Improved Gerschgorin disk estimator forsource enumeration with robustness against spatially non-uniform noise[J].IETRadar,Sonar&Navigation,2011,5(9):952-957.)。2013年，董姝敏则对盖尔圆半径进行压缩，并结合盖尔圆心值提出了一种基于自适应调整因子的GDE准则，该方法在信噪比为-4dB、快拍数大于2000的条件下检测准确率为80％；而在信噪比为15dB、快拍数为50时的检测准确率达到90％，虽然避免了对阵元和信源个数的要求，但该方法无法同时保证在低信噪比和小快拍数下的检测性能(董姝敏,梁国龙.改进的盖尔圆源数目估计方法[J].哈尔滨工程大学学报,2013,34(04):440-444.)。2016年，罗海坤提出了加权盖尔圆(WGDE)准则，该准则中的特征加权矩阵对增广盖尔圆矩阵中的半径值做特征加权变换从而进一步增大了信号盖尔圆和噪声盖尔圆之间半径的差异，使得在快拍数为64时的信源估计准确率就能达到90％，但是要求的信噪比为13dB，同样无法同时保证在低信噪比和小快拍数下的检测性能(一种基于改进盖氏圆盘法的信源数估计新算法[J].电子信息对抗技术,2016,31(06):24-28.)。2018年，褚鼎立则将GDE准则和AIC准则结合，提出了一种基于盖尔圆准则的信源数目判定改进方法，该方法克服了因特征值的无序性所导致估计错误的问题，使得该准则在信噪比为-5dB时的检测准确率为93％，但该方法依然要求较大的快拍数(褚鼎立,陈红,蔡晓霞.基于盖尔圆准则的信源数目估计改进算法[J].探测与控制学报,2018,40(04):109-115.)。最近，张冰玉在加权盖尔圆准则的基础上，提出了一种在增广加权盖尔圆矩阵中获取多重特征并融合的信源个数估计方法，该方法不仅在信源数只比阵元数少一个的情况下能够有效的估计信源个数，而且对于在低信噪比﹑小快拍数环境下的信源个数估计也有良好的性能，但前提是接收阵列为均匀圆阵(张冰玉,潘晴,田妮莉,Everett X.Wang.一种基于多重特征融合的信源个数估计方法[J/OL].计算机工程:1-5[2021-02-02].)。
基于盖尔圆盘理论的无人机目标个数判别方法虽然弥补了现有信息论准判别方法在不适用于色噪声条件下目标个数判定的问题，但基于盖尔圆盘准则的目标个数判定方法中存在人为选取的调节参数。对于不同信噪比、不同快拍数条件下，对应的最优调节参数不同，调节参数的选择至今尚没有较好的理论支撑，存在选取上的主观随意性。实际使用过程中，对于固定的基于经验的调节参数，无人机目标个数判别准确性急剧下降甚至完全失效。
发明内容
本发明的目的在于提供一种用于判定小型无人机目标个数的智能辨识方法，通过构建观测信号数据集的盖尔圆矩阵并绘制盖尔圆图，训练分类神经网络判别小型无人机个数。避免了GDE准则中最佳调节参数选取问题，同时提升了对无人机目标个数判别的环境适应能力。
为了实现上述目的，本发明采取以下技术方案：
一种用于判定小型无人机目标个数的智能辨识方法，其特征在于，包括以下步骤：
S1：将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号；
S2：对各个通道接收信号的数据集进行通道校正，获得各通道校正后观测信号数据集；
S3：各通道校正后观测信号数据集求协方差矩阵，再经酉变换计算得到盖尔圆矩阵；
S4：通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络(具体包括全连接类神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络)；
S5：改变预先设定的无人机个数，多次重复步骤S1-S4，生成一定数量的训练样本用于训练步骤S4中所述的分类神经网络，直至分类神经网络的准确率达到要求。
进一步地，步骤S1中所述的将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号，包括：
所述阵列包含M个阵元，间距为d，采样快拍数L，无人机个数为F，假设F个远场窄带信号入射到阵列上，其中F＝1,2，...，K，且K＜M，在阵列信号处理领域中将接收信号表示为:x(t)＝As(t)+n(t),式中x(t)＝[x1(t),x2(t),...,xM(t)]T为t时刻的接收信号矢量,s(t)＝[s1(t)，...，sk(t)]T为t时刻信号源矢量，n(t)＝[n1(t)，...nM(t)]T为t时刻色噪声矢量，a＝[α(Θ1)，...，α(Θk)]为阵列流型矢量，Θk＝(φk，θk)为第k个信号源的入射角，其中φk和θk分别为第k个信号源的方位角和俯仰角，且φk∈(0,2π)、θk∈(0,π/2)，
Figure BDA0003020049880000031
为第k个信号源的导向矢量，βm＝2π(m-1)/M为逆时针第m个阵元与X轴的夹角，
Figure BDA0003020049880000041
m＝1,2，...，M，λ为载波波长。
进一步地，步骤S2中所述的对各个通道接收信号的数据集进行通道校正，获得各通道观测信号数据集，包括：
提取通道校正系数：
令校正信号数据集为J＝[a1+b1i，...，aM+bMi]T，i为虚数单位，则通道校正系数矩阵为：
Figure BDA0003020049880000042
式中，T为通道校正系数矩阵。
通道校正：
设各通道接收信号的数据集为：
Figure BDA0003020049880000043
式中，Y表示各通道接收信号的数据集，cmn+emni为矩阵Y中第m行，第n列的元素，N为每一通道数据集的数据点最大个数，m＝1,2，…，M；n＝1,2，…N。
则X＝Diag(T)·Y所示的矩阵为经通道校正后各通道观测信号数据集，式中，Diag(T)为以矩阵T中的元素为对角元素的对角矩阵。
进一步地，步骤S3中所述各通道观测信号数据集求得协方差矩阵，再经酉变换计算得到盖尔圆矩阵，包括：
观测信号数据集的协方差矩阵为
Figure BDA0003020049880000044
其中，(.)H为共轭转置，对Rxx作酉变换，得到盖尔圆矩阵，即：
Figure BDA0003020049880000045
上式中，
Figure BDA0003020049880000051
UM-1为RM-1的特征向量构成的酉矩阵，RM-1为Rxx去除最后一行最后一列后构成的子矩阵。盖尔圆的圆心和半径分别为ci＝λi，ri＝|ρi|。信号盖尔圆半径远大于零，噪声盖尔圆半径则接近于零。
进一步地，步骤S4中所述通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络(具体包括全连接神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络)，包括：
通过盖尔圆矩阵获得绘制盖尔圆图所需的圆心和半径：
圆心矢量c＝diag(real(RT))/sum(diag(real(RT)))，式中，c＝[c1，c2，...，cm，...，cM]T为盖尔圆圆心矢量，cm为第m个盖尔圆圆心，real(.)为括号内矩阵的所有元素取实部后构成的矩阵，diag(.)为取括号内的矩阵的对角线上的元素并构成一列向量。sum(.)表示对括号内的矩阵的所有元素求和。
半径矢量ru＝abs(RT(M))/sum(abs(RT(M)))，式中，ru＝(r1，...，rm，...，rM)T盖尔圆半径矢量，rm为第m个盖尔圆半径，RT(M)为矩阵RT的最后一列构成的列向量，abs(.)为括号内的矩阵中的每一个元素取复数幅值后构成的矩阵。
在二维直角坐标系中，以cm·cos(π/4)为圆心点横坐标，cm·sin(π/4)为圆心点纵坐标，rm为半径在二维直角坐标系下绘制盖尔圆图并标记后依次送入分类神经网络中(具体包括全连接神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络)。
进一步地，步骤e中所述改变预先设定的无人机个数，多次重复步骤a-d，生成一定数量的训练样本用于训练步骤d中所述的分类神经网络，直至分类神经网络的准确率达到要求，包括：
假设目标无人机个数为F，依次设定F＝1,2，...，K，K为目标无人机个数最大值,且K＜M。重复步骤a-d，每设定一个无人机个数F采集H次经无人机反射过来的回波信号，取前G次采集的信号为训练集，则训练样本个数为Q＝K·G，其余为测试集，测试样本个数为W＝K·(H-G)。当无人机个数为F时，将对应的H个样本标记为F，神经网络目标输出为F。当将测试集送入神经网络中分类的准确率达到95％以上时停止训练，同种训练模式下，准确率最先达到95％以上的神经网络为最优网络。
本发明与现有技术相比，具有以下有益效果：
本申请先将各通道接收信号数据集进行通道校正，构建观测信号数据集的盖尔圆矩阵并在坐标系下绘制盖尔圆图，通过分类神经网络判别集群无人机目标个数，避免了GDE准则中调整因子的选择至今尚没有较好的理论支持，有较大的不确定性的缺陷，同时进一步提高GDE准则在低信噪比，小快拍数条件下的性能。
本发明提出了一种用于判定小型无人机目标数目的智能辨识方法。利用外辐射源发射一定频率的信号，由阵列天线接收经集群无人机反射后的回波信号，构建经过通道校正后的观测信号数据集的协方差矩阵并对协方差矩阵作酉变换获得盖尔圆矩阵，从中获得能将信号盖尔圆和噪声盖尔圆可视化的盖尔圆心值和盖尔圆半径并在坐标系中绘制盖尔圆可视化图，最后将图片送入分类神经网络中进行训练，直到神经网络的准确率达到要求。
附图说明
图1为本发明提出的判别方法的流程图；
图2为盖尔圆可视化制作流程图；
图3为实施例中无人机目标数目为1(F＝1)时的盖尔圆可视化训练样本图；
图4为实施例中无人机目标数目为2(F＝2)时的盖尔圆可视化训练样本图；
图5为实施例中无人机目标数目为3(F＝3)时的盖尔圆可视化训练样本图；
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
实施例
参见图1，本申请公开了一种用于判定小型无人机目标个数的智能辨识方法，包括以下步骤：
步骤A、将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号：
A1.所述阵列包含M个阵元，间距为d，采样快拍数L，无人机个数为F，假设F个远场窄带信号入射到阵列上，其中F＝1,2，...，K，且K＜M，在阵列信号处理领域中将接收信号表示为:x(t)＝As(t)+n(t),式中x(t)＝[x1(t),x2(t),...,xM(t)]T为t时刻的接收信号矢量,s(t)＝[s1(t)，...，sk(t)]T为t时刻信号源矢量，n(t)＝[n1(t)，...nM(t)]T为t时刻色噪声矢量，a＝[α(Θ1)，...，α(Θk)]为阵列流型矢量，Θk＝(φk，θk)为第k个信号源的入射角，其中φk和θk分别为第k个信号源的方位角和俯仰角，且φk∈(0,2π)、θk∈(0,π/2)，
Figure BDA0003020049880000071
为第k个信号源的导向矢量，βm＝2π(m-1)/M为逆时针第m个阵元与X轴的夹角，
Figure BDA0003020049880000072
m＝1,2，...，M，λ为载波波长。
步骤B、对各个通道接收信号的数据集进行通道校正，获得各通道观测信号数据集，包括：
B1.提取通道校正系数：
B11.令校正信号数据集为J＝[a1+b1i，...，aM+bMi]T，i为虚数单位，则通道校正系数矩阵为：
Figure BDA0003020049880000073
式中，T为通道校正系数矩阵。
B2.通道校正：
B21.设各通道接收信号的数据集为：
Figure BDA0003020049880000074
式中，Y表示各通道接收信号的数据集，cmn+emni为矩阵Y中第m行，第n列的元素，N为每一通道数据集的数据点最大个数，m＝1,2，…，M，n＝1,2，…N。
B22.则X＝Diag(T)·Y所示的矩阵为经通道校正后各通道观测信号数据集，式中，Diag(T)为以矩阵T中的元素为对角元素的对角矩阵。
步骤C、各通道观测信号数据集求得协方差矩阵，再经酉变换计算得到盖尔圆矩阵：
C1.观测信号数据集的协方差矩阵为
Figure BDA0003020049880000081
其中，(.)H为共轭转置，对Rxx作酉变换，得到盖尔圆矩阵，即：
Figure BDA0003020049880000082
上式中，
Figure BDA0003020049880000083
UM-1为RM-1的特征向量构成的酉矩阵，RM-1为Rxx去除最后一行最后一列后构成的子矩阵。盖尔圆的圆心和半径分别为ci＝λi，ri＝|ρi|。信号盖尔圆半径远大于零，噪声盖尔圆半径则接近于零。
步骤D、参见图2，通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络(具体包括全连接神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络)：
D1.通过盖尔圆矩阵获得绘制盖尔圆图所需的圆心和半径：
D11.圆心矢量c＝diag(real(RT))/sum(diag(real(RT)))，式中，c＝[c1，c2，...，cm，...，cM]T为盖尔圆圆心矢量，cm为第m个盖尔圆圆心，real(.)为括号内矩阵的所有元素取实部后构成的矩阵，diag(.)为取括号内的矩阵的对角线上的元素并构成一列向量。sum(.)表示对括号内的矩阵的所有元素求和。
D12.半径矢量ru＝abs(RT(M))/sum(abs(RT(M)))，式中，ru＝(r1，...，rm，...，rM)T盖尔圆半径矢量，rm为第m个盖尔圆半径，RT(M)为矩阵RT的最后一列构成的列向量，abs(.)为括号内的矩阵中的每一个元素取复数幅值后构成的矩阵。
D2.在二维直角坐标系中，以cm·cos(π/4)为圆心点横坐标，cm·sin(π/4)为圆心点纵坐标，rm为半径在二维直角坐标系下绘制盖尔圆图并标记后依次送入分类神经网络中(具体包括全连接神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络)。
步骤E、改变预先设定的无人机个数，多次重复步骤A-D，生成一定数量的训练样本用于训练步骤D中所述的分类神经网络，直至分类神经网络的准确率达到要求：
E1.假设目标无人机个数为F，依次设定F＝1,2，...，K，K为目标无人机个数最大值,且K＜M。重复步骤A-D，每设定一个无人机个数F采集H次经无人机反射过来的回波信号，取前G次采集的信号为训练集，则训练样本个数为Q＝K·G，其余为测试集，测试样本个数为W＝K·(H-G)。当无人机个数为F时，将对应的H个样本标记为F，神经网络目标输出为F。当将测试集送入神经网络中分类的准确率达到95％以上时停止训练，同种训练模式下，准确率最先达到95％以上的神经网络为最优网络。
图3，图4，图5分别展示了F＝1,2,3时通过本实施例获得的盖尔圆可视化训练样本图，通过本实例，卷积类神经网络的准确率最先达到95％以上。
按照上述实施例，便可很好地实现本发明。值得说明的是，基于上述结构设计的前提下，为解决同样的技术问题，即使在本发明上做出的一些无实质性的改动或润色，所采用的技术方案的实质仍然与本发明一样，故其也应当在本发明的保护范围内。

Claims (7)

1.一种用于判定小型无人机目标个数的智能辨识方法，其特征在于，包括以下步骤：
S1：将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号；
S2：对各个通道接收信号的数据集进行通道校正，获得各通道校正后观测信号数据集；
S3：各通道校正后观测信号数据集求协方差矩阵，再经酉变换计算得到盖尔圆矩阵；
S4：通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络；
S5：改变预先设定的无人机个数，多次重复步骤S1～S4，生成一定数量的训练样本用于训练步骤S4中所述的分类神经网络，直至分类神经网络的准确率达到要求。
2.根据权利要求1所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，步骤S1中所述的将包含有多个阵元的阵列在色噪声环境的远场中任意放置，在预先设定无人机个数和采样快拍数下，采样经无人机反射过来的回波信号，具体方法为：
所述阵列包含M个阵元，间距为d，采样快拍数L，无人机个数为F，F个远场窄带信号入射到阵列上，其中F＝1,2，...，K，且K＜M，在阵列信号处理领域中将接收信号表示为:x(t)＝As(t)+n(t),式中x(t)＝[x1(t),x2(t),...,xM(t)]T为t时刻的接收信号矢量,s(t)＝[s1(t)，...，sk(t)]T为t时刻信号源矢量，n(t)＝[n1(t)，...nM(t)]T为t时刻色噪声矢量，a＝[α(Θ1)，...，α(Θk)]为阵列流型矢量，Θk＝(φk，θk)为第k个信号源的入射角，其中φk和θk分别为第k个信号源的方位角和俯仰角，且φk∈(0,2π)、θk∈(0,π/2)，
Figure FDA0003020049870000011
为第k个信号源的导向矢量，βm＝2π(m-1)/M为逆时针第m个阵元与X轴的夹角，
Figure FDA0003020049870000012
λ为载波波长。
3.根据权利要求1所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，步骤S2中所述的对各个通道接收信号的数据集进行通道校正，获得各通道观测信号数据集，具体方法为：
提取通道校正系数：
令校正信号数据集为J＝[a1+b1i，...，aM+bMi]T，i为虚数单位，则通道校正系数矩阵为：
Figure FDA0003020049870000021
式中，T为通道校正系数矩阵。
通道校正：
设各通道接收信号的数据集为：
Figure FDA0003020049870000022
式中，Y表示各通道接收信号的数据集，cmn+emni为矩阵Y中第m行，第n列的元素，N为每一通道数据集的数据点最大个数，m＝1,2，…，M，n＝1,2，…N。
则X＝Diag(T)·Y所示的矩阵为经通道校正后各通道观测信号数据集，式中，Diag(T)为以矩阵T中的元素为对角元素的对角矩阵。
4.根据权利要求1所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，步骤S3中所述各通道观测信号数据集求得协方差矩阵，再经酉变换计算得到盖尔圆矩阵，具体方法为：
观测信号数据集的协方差矩阵为
Figure FDA0003020049870000023
其中，(.)H为共轭转置，对Rxx作酉变换，得到盖尔圆矩阵，即：
Figure FDA0003020049870000024
上式中，
Figure FDA0003020049870000025
UM-1为RM-1的特征向量构成的酉矩阵，RM-1为Rxx去除最后一行最后一列后构成的子矩阵。盖尔圆的圆心和半径分别为ci＝λi，ri＝|ρi|。信号盖尔圆半径远大于零，噪声盖尔圆半径则接近于零。
5.根据权利要求1所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，步骤S4中所述通过盖尔圆矩阵确定用于集群无人机个数估计的盖尔圆圆心和盖尔圆半径，在二维直角坐标系下绘制盖尔圆图，并将二维盖尔圆图标记后送入分类神经网络，具体方法为：
通过盖尔圆矩阵获得绘制盖尔圆图所需的圆心和半径：
圆心矢量c＝diag(real(RT))/sum(diag(real(RT)))，式中，c＝[c1，c2，...，cm，...，cM]T为盖尔圆圆心矢量，cm为第m个盖尔圆圆心，real(.)为括号内矩阵的所有元素取实部后构成的矩阵，diag(.)为取括号内的矩阵的对角线上的元素并构成一列向量。sum(.)表示对括号内的矩阵的所有元素求和；
半径矢量ru＝abs(RT(M))/sum(abs(RT(M)))，式中，ru＝(r1，...，rm，...，rM)T盖尔圆半径矢量，rm为第m个盖尔圆半径，RT(M)为矩阵RT的最后一列构成的列向量，abs(.)为括号内的矩阵中的每一个元素取复数幅值后构成的矩阵；
在二维直角坐标系中，以cm·cos(π/4)为圆心点横坐标，cm·sin(π/4)为圆心点纵坐标，rm为半径在二维直角坐标系下绘制盖尔圆图并标记后依次送入分类神经网络中。
6.根据权利要求1所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，步骤S5中所述改变预先设定的无人机个数，多次重复步骤a-d，生成一定数量的训练样本用于训练步骤d中所述的分类神经网络，直至分类神经网络的准确率达到要求，具体方法为：
目标无人机个数为F，依次设定F＝1,2，...，K，K为目标无人机个数最大值,且K＜M。重复步骤a-d，每设定一个无人机个数F采集H次经无人机反射过来的回波信号，取前G次采集的信号为训练集，则训练样本个数为Q＝K·G，其余为测试集，测试样本个数为W＝K·(H-G)。当无人机个数为F时，将对应的H个样本标记为F，神经网络目标输出为F。当将测试集送入神经网络中分类的准确率达到95％以上时停止训练，同种训练模式下，准确率最先达到95％以上的神经网络为最优网络。
7.根据权利要求5所述的用于判定小型无人机目标个数的智能辨识方法，其特征在于，
分类神经网络包括全连接神经网络、卷积类神经网络、残差神经网络、自编码器神经网络、以及由上述四类网络组合构成的混合神经网络。
CN202110400484.9A 2021-04-14 2021-04-14 一种用于判定小型无人机目标个数的智能辨识方法 Active CN113514808B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110400484.9A CN113514808B (zh) 2021-04-14 2021-04-14 一种用于判定小型无人机目标个数的智能辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110400484.9A CN113514808B (zh) 2021-04-14 2021-04-14 一种用于判定小型无人机目标个数的智能辨识方法

Publications (2)

Publication Number Publication Date
CN113514808A true CN113514808A (zh) 2021-10-19
CN113514808B CN113514808B (zh) 2023-04-25

Family

ID=78061361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110400484.9A Active CN113514808B (zh) 2021-04-14 2021-04-14 一种用于判定小型无人机目标个数的智能辨识方法

Country Status (1)

Country Link
CN (1) CN113514808B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1096077A (en) * 1964-04-09 1967-12-20 Thomson Houston Comp Francaise Method and apparatus for adding a plurality of primary sources
WO1993010500A1 (en) * 1991-11-13 1993-05-27 Iowa State University Research Foundation, Inc. Improved memory capacity neural network
US20080231505A1 (en) * 2007-03-23 2008-09-25 Weiqing Zhu Method of Source Number Estimation and Its Application in Method of Direction of Arrival Estimation
CN104880723A (zh) * 2015-04-29 2015-09-02 中国电子科技集团公司第三十六研究所 一种单星测向定位地面辐射源的方法和装置
CN106019212A (zh) * 2016-05-16 2016-10-12 中国电子科技集团公司第十研究所 一种用于多通道间非理想因素的校正方法
CN106059639A (zh) * 2016-06-21 2016-10-26 西安电子科技大学 基于矩阵盖尔圆的发射天线数目盲估计方法
CN106707228A (zh) * 2016-11-30 2017-05-24 广东工业大学 一种结合四阶累积量与盖尔圆改进的信号源个数估计方法
CN110346793A (zh) * 2019-07-24 2019-10-18 电子科技大学 一种分布式阵列雷达迭代自适应高分辨成像方法
CN111007488A (zh) * 2019-11-21 2020-04-14 中国人民解放军63892部队 基于盖尔圆变换和修正Rao得分检验的信源数估计方法
CN111738291A (zh) * 2020-05-18 2020-10-02 广东工业大学 一种基于增广加权盖尔圆矩阵的信源个数估计建模方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1096077A (en) * 1964-04-09 1967-12-20 Thomson Houston Comp Francaise Method and apparatus for adding a plurality of primary sources
WO1993010500A1 (en) * 1991-11-13 1993-05-27 Iowa State University Research Foundation, Inc. Improved memory capacity neural network
US20080231505A1 (en) * 2007-03-23 2008-09-25 Weiqing Zhu Method of Source Number Estimation and Its Application in Method of Direction of Arrival Estimation
CN104880723A (zh) * 2015-04-29 2015-09-02 中国电子科技集团公司第三十六研究所 一种单星测向定位地面辐射源的方法和装置
CN106019212A (zh) * 2016-05-16 2016-10-12 中国电子科技集团公司第十研究所 一种用于多通道间非理想因素的校正方法
CN106059639A (zh) * 2016-06-21 2016-10-26 西安电子科技大学 基于矩阵盖尔圆的发射天线数目盲估计方法
CN106707228A (zh) * 2016-11-30 2017-05-24 广东工业大学 一种结合四阶累积量与盖尔圆改进的信号源个数估计方法
CN110346793A (zh) * 2019-07-24 2019-10-18 电子科技大学 一种分布式阵列雷达迭代自适应高分辨成像方法
CN111007488A (zh) * 2019-11-21 2020-04-14 中国人民解放军63892部队 基于盖尔圆变换和修正Rao得分检验的信源数估计方法
CN111738291A (zh) * 2020-05-18 2020-10-02 广东工业大学 一种基于增广加权盖尔圆矩阵的信源个数估计建模方法

Also Published As

Publication number Publication date
CN113514808B (zh) 2023-04-25

Similar Documents

Publication Publication Date Title
CN110045322B (zh) 一种基于高分辨测向语图智能识别的短波自动测向方法
CN112346030A (zh) 无人机群的超分辨波达方向估计方法
CN112034418A (zh) 基于频域Bark子带的波束扫描方法及声源定向装置
CN110806557A (zh) 一种冲击噪声背景下单基地mimo雷达doa估计方法
CN111896913B (zh) 高频雷达单极子/交叉环天线通道增益校准方法及装置
CN112630784A (zh) 基于凸优化和神经网络的平面阵列幅相误差校正方法
CN111913155A (zh) 基于阵列雷达的二维doa估计方法
CN113884980A (zh) 基于加权特征空间投影的鲁棒盖氏圆信源数估计方法
CN111738291A (zh) 一种基于增广加权盖尔圆矩阵的信源个数估计建模方法
CN111812581A (zh) 基于原子范数的球面阵列声源波达方向估计方法
CN111817803A (zh) 一种基于相关系数和k-均值聚类算法的频谱感知方法、系统及计算机可读存储介质
CN111142062A (zh) 一种利用Toeplitz特性的无栅格目标波达方向估计方法
CN111551892B (zh) 一种稳健自适应波束形成方法及装置
CN113514808A (zh) 一种用于判定小型无人机目标个数的智能辨识方法
CN112699992A (zh) 基于卷积神经网络的混合源定位方法
CN110261814B (zh) 基于空间谱重构和导向矢量直接估计的波束形成方法
CN112630760A (zh) 多目标条件下的强干扰抑制波束形成器设计方法
WO2022166477A1 (zh) 定位方法、装置、基站、计算机设备和存储介质
CN112415469B (zh) 一种两维数字阵列雷达快速干扰测向方法
CN111323750B (zh) 一种基于声矢量阵列网络的直接定位方法
CN114966523A (zh) 大规模mimo测向中基于幂迭代算法的doa估计方法
Polikarovskykh et al. Neural Network Method of Directing Finder Signals Processing in Perimeter Protection Systems
CN113075613A (zh) 一种卫星移动通信系统中的干扰测向方法
CN112737644A (zh) 无人机蜂群的自定位方法及装置
CN111337879B (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