基于深度卷积神经网络的雷达辐射源识别方法与流程

本发明属于计算机
技术领域：
，具体是雷达辐射源信号识别领域，具体涉及一种基于深度卷积神经网络的雷达辐射源识别方法。
背景技术：
：雷达辐射源信号识别当今电子战的重要功能之一，主要就是通过各个辐射源在发射信号时所表现出来的差异，将被测雷达辐射源参数与预先累积的参数进行比较，来对各个辐射源进行识别，从而实现对雷达信号进行截获、定位、分析和识别，进而可以对接收的信号进行分辨，判断出究竟是来自哪个雷达辐射源，最终完成对雷达辐射源信号地识别。所以雷达辐射源识别算法是非常具有挑战性的问题。近年来，深度学习受到广大研究者的欢迎，进而基于深度学习的方法在许多研究领域取得了显著的成功，例如语音识别，图像分类、图像识别、物体检测、计算机视觉和大数据特征提取等方面都具有广泛的应用。在雷达辐射源识别算法发展中，一系列方法相继提出。将bp神经网络应用到辐射源识别算法中，但是由于多个神经元组成复杂多变量的非线性组合，增加了神经网络在雷达识别算法上的难度，促使研究者另辟新路。进而提出支持向量机(svm)运用到雷达辐射源的信号识别中，svm在处理小样本数据时，具有更加出色的识别能力，但是在处理语音或者图像这些问题，svm会有明显的不足，它的泛化能力受到约束。为了解决上述中存在的问题，许多研究者将深度学习识别算法引入到雷达辐射源，这在深度学习应用领域中也是非常新颖，因此基于深度学习的雷达辐射源识别算法越来越受到人们的重视。研究者进行大量的探索，提出了一些方法，其中一种称为深度卷积神经网络，它可以用原始图像作为输入，直接从图像中提特征。将特征提取与分类器相结合，用卷积神经网络进行数学建模，可以实现端到端学习，具有较强的泛化能力。技术实现要素：本发明针对现有技术的不足，提供一种基于深度卷积神经网络的雷达辐射源识别方法。本发明将雷达辐射源模型产生的一维波形信号，经过短时傅里叶变换(stft)处理之后转化为频谱图，再针对一维波形信号和频谱图设计不同的网络结构。本发明方法具体包括以下步骤：步骤(1).数据产生：选取八种雷达信号，分别为lfm信号、lfmcw、lfm-bc信号、frank-lfm信号、s型nlfm信号、costas编码信号、fsk/psk信号、p3码编码信号；其中，对于lfm信号、lfmcw、lfm-bc信号、frank-lfm信号、s型nlfm信号改变带宽b，从而改变调频斜率μ＝b/t的大小，分别产生五类样本数据；对于costas编码、fsk/psk信号，将其编码进行全排列，产生两类数据样本；对于p3码编码信号，令ncτc＝25，产生一类数据样本；分别从8类样本数据中取出信号，将其波形提取出来，将这些雷达波形信号记为：s＝{(si(t),yi)|i∈λ}；其中，si(t)是第i个雷达信号样本，si(t)＝[si(0),si(1),...,si(n-1)]，i为雷达信号样本个数，n是采样点的个数，yi∈{0,..,c-1}表示第i个雷达信号样本si(t)的类别，一共有c类雷达信号，λ是样本的索引集合。步骤(2).数据集预处理：将步骤(1)产生的8类样本数据进行短时傅里叶变换，在变换过程中对数据进行翻转，实现数据增强，获得翻两倍的数据集；经过短时傅里叶变换之后，得到频谱图。短时傅里叶变换的数学描述如下：其中n为短时傅里叶变换过程中需要时间，k为窗滑动的次数，w(·)为窗函数，j为虚数，stfti(n,k)为第i个雷达信号波形经过短时傅里叶变换的结果，具有2维结构。步骤(3).建立深度卷积神经网络模型获取雷达辐射源信号特征：将步骤(1)产生的波形信号和步骤(2)进行翻转之后的频谱图分别输入到深度卷积神经网络，进行深度卷积处理，得到经过短时傅里叶变换后频谱图的特征。步骤(4).雷达辐射源信号分类：由步骤(3)得到特征图中的雷达信号特征信息，用softmax分类器将由步骤(1)产生不同类型的雷达辐射源信号进行分类识别；其中softmax分类器的模型公式为：其中η为雷达辐射源信号特征，y为分类中的某一类，k为分类的全部总数，p为分类器的输出概率，按照输出概率值进行最终的结果分类。进一步，步骤(1)中五种信号改变带宽b具体如下：①lfm表示为：式中，其中，s(t)为产生的lfm信号，t为时间变量，a为幅度、t为线性调频信号的持续时间、b为lfm信号带宽、f0为起始频率、为初始相位，μ＝b/t为调频斜率，j为虚数；lfm的瞬时频率f(t)＝f0+μt(2)；选取雷达信号仿真参数为信号幅度a＝1，起始频率f0＝3ghz，采样频率fs＝1024mhz，脉冲宽度t＝25μs，带宽b＝20mhz，初始相位②lfmcw表示为：令一次观测时间为tobs，t1为一个lfm调制周期的时间长度，则观测时间内的lfm调制脉冲数q＝tobs/t1，h(t)为产生的lfmcw信号，μ＝b/t1为调频斜率；mod(·)表示取模运算，mod(a,b)表示a除以b所得的余数，τbias为信号的时间偏移；选取a＝1，f0＝3ghz，fs＝1024mhz，b＝20mhz，t1＝1μs，时间偏移τbias＝0。③lfm-bc表示为：式中，l(t)为产生的lfm-bc信号，nt为码元数，f0为起始频率，为二相编码相位参数取1或-1，π；选取a＝1，f0＝3ghz，fs＝1024mhz，t＝5μs，b＝20mhz。④frank-lfm表示为：其相位矩阵中的元素将相位矩阵中的元素按行依次串行排列，可得长度为nm＝m2的相位序列式中，f(t)为产生的frank-lfm信号，nm为码元数，τm为码元宽度；t＝nmτm，t为线性调频信号的持续时间；为编码相位参数，取1或-1；f0为起始频率，μ＝b/t为调频斜率，b为lfm信号带宽；选取a＝1，fs＝1024mhz，b＝20mhz，τm＝1μs，m＝5。⑤s型nlfm表示为：w(f)＝0.54+0.46cos(2πf/b)(6)；得到信号的群延时：t(f)＝0.54kf+(0.27kb/π)sin(2πf/b)，-b/2≤f≤b/2(7)；令t(f)|f＝b/2＝t2/2，则常数k＝(t2/b)/0.54，带入式(7)得到：其中，k为常数，t2为信号持续脉宽，b为lfm信号带宽；w(f)为产生的s型nlfm；选取a＝1，f0＝3ghz，fs＝1024mhz，b＝20mhz，t2＝25μs。根据上述给出的信号模型，对lfm、lfmcw、lfm-bc、frank-lfm、s型nlfm这五类信号改变b的值，从而改变调频斜率μ＝b/t的大小，分别产生五类样本数据。进一步，步骤(1)中对于costas编码、fsk/psk信号，将其编码进行全排列，具体是：①设线性调频信号的持续时间为t，将t分为p个相同宽度的码元，码元宽度τp＝t/p，频率编码长度为p，编码序列为a＝{a1,a2,…,ap}；costas编码表示为：式中，c(t)为产生的costas编码信号，fp＝ap/τp，ap是编码序列中的元素，f0为起始频率，fp为跳频频率；选取a＝1，f0＝3ghz，fs＝1024mhz，τp＝5μs，costas码序列[21534]。②若已知fsk信号的跳频序列为{f0,f1,…,fk-1}，bpsk信号的码元为fsk信号的持续周期为t4，则信号产生如下：将t4等分为k个宽度为tf的子脉冲，在各个子脉冲间进行频率跳变，再将每个子脉冲等分为np个宽度为tp的跳频子脉冲，对各个跳频子脉冲进行相位编码，得到fsk/psk信号，表示为：式中，k(t)为产生的fsk/psk信号，调制相位θi＝πci，ci为二相编码，为初始相位；选取a＝1，f0＝3ghz，fs＝1024mhz，tp＝1μs，fsk序列采用costas序列即[21534]，psk采用5位barker码[111-11]。对于costas编码、fsk/psk信号，将其编码进行全排列，产生两类数据样本。进一步，步骤(1)中对于p3码编码信号，具体是：p3码编码信号表示为：式中，p(t)为产生的p3码编码信号nc为码元数，τc为码元宽度，t＝ncτc为线性调频信号的持续时间，f0为起始频率，为初始相位，p3码的第k个码元相位表示为其中k＝0,1,2,…,nc-1；选取a＝1，f0＝3ghz，fs＝1024mhz，τc＝1μs，nc＝25；对于p3编码信号，保证nctc＝25，产生一类数据样本。经过stft处理得到的频谱图作为深度卷积神经网络的输入，可使雷达辐射源信号分类识别的准确率比波形信号作为网络输入高出15％。对比目前雷达辐射源识别方法，本次发明方法更能有效的对雷达信号进行分类。本发明方法对stft处理后的频谱图的泛化能力更强，深度卷积神经网络确实学到了具有区分性的特征表示，从而得到更高的分类识别准确率。也说明频谱图给网络的信息更多，能提取出的特征表示也更多，所以深度卷积神经网络对频谱图的识别效果更好。本方法运用卷积神经网络进行网络结构建模，极大地减小了模型的复杂程度，具有旋转、位移和缩放不变性等特点，可以将提取的雷达辐射源信号特征精确地区分出来，最终完成雷达辐射源信号的分类识别。附图说明图1.本发明中8类雷达信号的实部波形图；图2.对应图1中的雷达辐射源信号实部波形经过stft后相对应的频谱图；图3.stft之后频谱图所设计的网络架构图；图4.网络模型混淆矩阵示意图具体实施方式下面结合附图对本发明进行进一步的说明。本发明提出基于深度卷积神经网络的雷达辐射源识别方法，下面根据该方法的流程，并结合附图说明，其具体实施方式如下：1.数据产生：如图1，本发明选取8种雷达信号，分别为线性调频信号(lfm)、线性调频连续波(lfmcw)、lfm-bc信号、frank-lfm信号、s型nlfm信号、costas编码信号、p3码编码信号、fsk/psk信号。如步骤(1)所述，对lfm、lfmcw、lfm-bc、frank-lfm、s型nlfm这五类信号改变带宽b的值，从而改变调频斜率μ＝b/t的大小，产生五类数据样本；对costas编码、sk/psk将其编码进行全排列，产生两类数据样本；对p3码编码信号，保证ncτc＝25，产生一类数据，据此可得八类数据样本，并将八类雷达信号的波形提取出来。图1为八类雷达信号的实部波形图(a-h依次为：costas编码、frank-lfm、fsk/psk、lfm-bc、lfm、lfmcw、s型nlfm、p3码编码信号)。2.数据集预处理：由步骤(1)可知产生八类样本数据，但是每类样本个数对深度卷积神经网络来说，数据量较小，在处理雷达辐射源信号分类识别问题上很难，所以需要通过数据增强来扩大数据量从而有效地防止过拟合问题。为解决过拟合问题，我们将步骤(1)产生的样本进行短时傅里叶变换(stft)，在变换过程中对数据进行翻转，实现数据增强，从而获得翻两倍的数据集。另外经过stft之后，会得到频谱图。经过波形信号转换为频谱图，还有数据的增强，这对之后网络结构的训练和雷达辐射源的分类识别会更有效。其中stft的数学描述为：其中w(·)为窗函数，stfti(n,k)为短时傅里叶变换结果，具有2维结构。3.参数的设置：在步骤2的stft过程中，采用‘hamming’窗，窗长大小为256，信号长度为25600，重叠点数为128个采样点，窗滑动次数为199。每个样本经过stft后会产生成为32*199的二维图像。图2是对应图1中8类雷达辐射源信号实部波形经过stft后得到相对应的频谱图。cnn网络结构在训练过程中选择迭代次数为20000，批操作大小为16，选择学习率为0.01，每隔2000次降低为原来的0.5倍，总共训练1000次。动量和权重衰减系数分别为0.9和0.000005。4.cnn网络结构设计：虽然在经过stft得到的频谱图能够描述雷达信号的时频特性，但是对于不同雷达信号的分类识别问题仍然还是需要提取具有区分性的特征表示。因此，通过建立深度卷积神经网络模型来解决该问题。设计不同的卷积神经网络模型来不断提高识别率，表1、表2分别是针对一维实部波形和stft后频谱图设计的两个不同的网络结构。表1一维波形实验的网络结构表2二维频谱实验的网络结构(表1中m1d1、m1d2是针对一维实部波形所选出的网络模型，表2中m2d1、m2d2是针对stft之后得到的实部和双通道频谱图所设计的网络结构。)其中在表2中m2d2的网络架构图如图3。5.雷达辐射源信号分类由步骤(3)得到特征图中的雷达信号特征信息，用softmax分类器将由步骤(1)产生不同类型的雷达辐射源信号进行分类识别；其中softmax分类器的模型公式为：其中η为雷达辐射源信号特征，y为分类中的某一类，k为分类的全部总数，p为分类器的输出概率，按照输出概率值进行最终的结果分类。本实施例对表1、表2中网络结构进行训练和测试，其本算法的仿真实验结果如表3所示：表3不同模型测试结果模型训练集准确率测试集准确率m1d179.17％77.34％m1d272.15％71.49％m2d189.0％89.35％m2d288.05％88.89％m3d191.28％90.9％m3d292.8％92.4％(注：m1d1、m1d2是一维实部波形设计的网络模型；m2d1、m2d2是stft后得到实部频谱图设计的网络结构；m3d1与m3d2是stft之后的双通道(实部和虚部)频谱图所设计的网络结构)。为了更清晰地看清每一类的分类情况，我们将分类结果转化为可视化工具——混淆矩阵。而且也可以从混淆矩阵中看雷达辐射源每类错分的个数。表3中m3d2的混淆矩阵如图4所示。当前第1页12