基于并行原型网络的小样本雷达辐射源信号识别方法与流程

1.本发明属于信号处理技术领域，特别涉及一种雷达信号识别方法，可用于电子情报侦察、电子支援和威胁告警系统。


背景技术：

2.随着电子信息领域的发展，电子对抗在电子情报侦察、电子支援和威胁告警系统中发挥着重要作用，雷达辐射源信号识别是电子对抗中的重要环节。
3.现阶段的雷达调制信号识别多是通过深度学习技术来实现，但是深度学习存在一个严重的问题，即需要大量的数据来支撑网络的训练。然而，新体制雷达的不断涌现和电磁环境的复杂化，导致大量高质量的辐射源信号标注样本难以获得，因此，如何在仅具有少量样本的情况下对同类型辐射源信号进行识别成为了目前学术界的热门研究方向。
4.当今小样本识别受到大量学者的关注，然而其主要集中在图像领域上，而雷达辐射源识别领域应用资料文献较少。其中，丁辰伟等人在2021年发表《基于小样本条件的雷达辐射源个体识别算法》，该方法对信号进行图像增强后利用对抗生成网络和卷积网络对其进行训练识别，该方法中数据的预处理部分包括短时傅里叶变换获得时频图，对时频图进行截取、灰度化、下采样等操作，存在信号预处理过程繁琐的问题。于浩洋等人在期刊《西安电子科技大学学报》2021年发表的《生成对抗网络小样本雷达调制信号识别算法》中，提出增强深度卷积生成对抗网络，通过改进的生成对抗网络生成与真实数据相近的数据，通过扩大训练数据的方式解决深度网络的训练问题，该类方法需要先后训练用于模拟数据的生成对抗网络和用于分类的卷积神经网络两个网络，需要耗费一定时间。


技术实现要素：

5.本发明提供了基于并行原型网络的小样本雷达辐射源信号识别方法，其通过以时域和频域输出结果的平均作为最终识别结果，有效利用信号的多个特征对信号进行分辨识别的同时避免直接将时频特征进行结合可能造成的特征冗余，并且以任务训练方式训练并行原型网络，通过不断适应每个具体任务，使并行原型网络在该训练方式下获取如何解决每一个任务的抽象概念，解决小样本问题。
6.为了实现上述技术效果，本发明通过以下技术方案予以实现。
7.基于并行原型网络的小样本雷达辐射源信号识别方法，包括以下步骤：
8.数据集的获得：利用仿生软件，仿生形成时频数据形式的已知雷达辐射信号数据集、样本雷达辐射信号数据集和待识别雷达射频信号数据集；所述已知雷达辐射信号数据集包括cw、lfm、nlfm和bpsk四类信号数据集，所述样本雷达辐射信号数据集和待识别雷达射频信号数据集均包括bfsk、qpsk和qfsk三类信号数据集；
9.元训练集、元测试集以及测试集的生成：分别将获取的已知雷达辐射信号数据集、样本雷达辐射信号数据集和待识别雷达射频信号数据集进行傅里叶变换，得到时域数据集和频域数据集共同构成的元训练集da、元测试集db以及测试集dc；
10.初始并行原型网络的生成：分别依次构建时域和频域的特征提取模块和度量模块，然后将时域和频域的度量模块组合，形成识别模块，以识别模块、2个所述度量模块以及2个所述特征提取模块组合，形成初始并行原型网络；
11.并行原型网络的训练和测试：将元训练集在初始并行原型网络中进行元训练任务的训练，得到第一权重，然后将元测试集在载有第一权重的训练并行原型网络上进行元测试集的继续训练，得到第二权重，用第二权重覆盖第一权重，得到训练好的并行原型网络；
12.信号识别：将测试集中获取的测试任务，输入到载有第二权重的并行原型网络中，获取得到雷达信号的识别测试结果。
13.本技术方案中，首先，通过元训练集(即包含雷达辐射信号数据集及其频域数据)训练网络，使网络获得如何解决每一个任务的抽象概念；再通过元测试集(即样本雷达辐射信号数据集及其频域数据)训练网络，使网络在上面已经学会如何解决每一个任务的基础上，学会如何解决当前设置的小样本问题(即识别qpsk,qfsk,bfsk三类信号)，实现对小样本的识别。
14.本技术方案中，其次，由于cw,lfm,nlfm,bpsk,bfsk,qpsk,qfsk这7类信号属于比较常用的信号类型，故本实验选取这7类信号作为示例说明原理。
15.这7类信号中可以任意选取几类信号作为“具有大量已知信号”的数据集，几类信号作为“小样本待识别”数据集，两个数据集不能有类别交集，为了模拟大量已知信号样本和小样本。
16.本技术方案选取bfsk,qpsk,qfsk三类信号作为小样本信号，生成少量的含有这三类信号的数据集作为现在获取的少量标记样本信号(记为b)，同时也生成含有大量这三类信号的数据集作为待识别数据集用以验证本方法的测试识别率(记为c)，并选取剩下的四类作为大量已知信号(记为a)。
17.通过本技术方案，用这四类大量已知信号集a和三类小样本信号集b训练网络，使网络最终能够识别这三类小样本信号，最后利用信号集c验证本文的小样本识别方法识别率。
18.本技术方案中，最后，时域和频域信号分别输入各自的特征提取模块，不仅可以更好的提取样本的时域特征和频域特征，并能充分利用小样本的时域和频域信息，为识别提供依据。进而，通过时域和频域度量模块分别计算时域和频域上待识别信号特征与已知信号特征原型的距离作为判断，并通过识别模块平均时域和频域的判断做出最终识别，由于综合了时域和频域的判断，并避免了在前期对时、频特征直接结合进行判断可能造成的特征冗余问题。
19.作为本发明的进一步改进，所述数据集的获取中，统一信噪比下，每类已知雷达辐射信号数据集、每类样本雷达辐射信号数据集和每类待识别雷达射频信号数据集的信号数量分别是1500-2500个、8-12个以及700-1300个。
20.本技术方案中，信噪比选择-4db模拟复杂的信号环境。小样本识别的目的和方法就是在仅有少量标注样本(即b)的情况下，借助大量其他类型的有标注数据(即a)学会如何识别b这些类型的样本，以便在下次采集到这些类型的未标记样本(即c)时可以识别并对其标注。因此，a数量大，b数量小，a每类信号选取2000个模拟大量已知，实际应用时并不规定必须2000个，只要是大量就可以，此处选择2000做实验仿真；b每类信号选取10个模拟少量
标注，实际应用时并不规定必须是10个，只是本实验选取10个作为少量，进行试验；c每类选取1000模拟与b同类的未标注样本，实际应用时并不是必须1000个，此处选择1000个做实验仿真，验证本方法的识别率。
21.作为本发明的进一步改进，所述元训练集、元测试集以及测试集的生成步骤中，还包括标准化处理，所述标准化处理具体为采用z-score标准化方法对元训练集、元测试集及测试集的数据进行预处理，获得均值为0、标准差为1的预处理后的元训练集、元测试集及测试集。
22.本技术方案中，通过标准化的处理，用以保证不同维度的数据具有相同的分布规模，更利于网络优化，缩短训练时间；由于本技术方案中，有多种类型的数据，进而通过这样的处理，使其更加统一化。
23.作为本发明的进一步改进，所述并行原型网络的生成中，所述时域特征提取模块包括四个卷积单元和平坦层，四个所述卷积单元中，前3个卷积层均包括卷积层、激活层、批量归一化层和最大池化层，第4个卷积层包括卷积层、激活层和批量归一化层。
24.本技术方案中，对卷积单元进行规定，使得卷积神经网络通过局部感受野、权重共享和降采样3中策略降低了网络模型的复杂度，减少了网络参数，使网络的适应性更强。
25.作为本发明的进一步改进，所述初始并行原型网络的生成中，时域和频域的特征提取模块的工作过程为：所述元训练集或元测试集输入任务的时域/频域支持集和时域/频域查询集，输出时域/频域支持集和时域/频域查询集的特征信号。
26.本技术方案中，时域和频域度量模块搭建方法完全相同，只是输入输出不同：具体地，时域度量模块输入时域支持集和查询集信号特征，输出时域查询集每个信号所属时域支持集每个信号类别的概率；频域度量模块输入频域支持集和查询集的信号特征，输出频域查询集的每个信号所属频域支持集每个信号类别的概率。
27.作为本发明的进一步改进，所述初始并行原型网络的生成中，所述度量模块的构建具体为：选择欧几里得距离作为度量方法，具有一个激活函数为softmax的激活层，输入时域/频域支持集和时域/频域查询集的信号特征，输出时域查询集的每个信号所属时域支持集每个信号类别的概率。
28.本技术方案中，时域度量模块和频域度量模块搭建方法完全相同，而通过激活层以及欧几里得等的处理，使得数据规范化，进而最后输出的信号更容易识别。
29.作为本发明的进一步改进，所述识别模块的形成具体为：通过平均时域与频域识别类别概率，并选取查询集每条信号所述支持集各类别概率中的最大概率作为查询集该条信号的所属类别进行输出。
30.本技术方案中，通过平均值函数对查询集各信号所属的类别概率进行时、频域判断结果综合，得到综合判断的查询集各信号所属支持集各类别的概率，具体如下：
[0031][0032]
通过最大值函数选取查询集每条信号所述支持集各类别概率中的最大概率作为查询集该条信号的所属类别。
[0033]
作为本发明的进一步改进，所述初始并行原型网络的生成之后，并行原型网络的训练和测试之前，还包括元训练任务的获取，具体为：
[0034]
1)在预处理后的元训练集中，在每类信号类型中，均随机抽取八条信号a1,a2,
a3,
···
,a8，将八条信号a1,a2,a3,
···
,a8中的前三条a1,a2,a3作为时域支持集t
sa
，后五条a4,a5,
···
a8作为时域查询集t
qa
；
[0035]
2)在预处理后的元训练集中，分别选择时域a1,a2,a3,
···
,a8八条信号对应的频域信号b1,b2,b3,
···
,b8，将b1,b2,b3作为频域支持集f
sa
，b4,b5,
···
b8作为频域查询集f
qa
；
[0036]
3)将所述t
sa
、t
qa
、f
sa
、f
qa
共同组成元训练任务u；
[0037]
4)重复步骤1)至步骤3)d次，完成元训练任务u的获取。
[0038]
本技术方案中，通过多次的元训练任务，使得并行原型网络通过元训练任务学会如何处理不同的任务，获得任务处理的抽象概念，为后续的进一步训练提供好的基础。
[0039]
作为本发明的进一步改进，所述元训练任务的获取之后，还包括元测试任务的获取，具体为：
[0040]
a)在预处理后的元测试集db的三类信号中，每类信号中，均随机抽取e1,e2,e3,
···
,e8八条信号，将e1,e2,e3作为时域支持集t
sb
，e4,e5,
···
e8作为时域查询集t
qb
；
[0041]
b)在预处理后的元测试集db频域信号的bfsk，qpsk，qfsk这三类信号中，分别选择6.1.1)中时域e1,e2,e3,
···
,e8八条信号对应的频域信号c1,c2,c3,
···
,c8，将c1,c2,c3作为频域支持集f
sb
，c4,c5,
···
c8作为频域查询集f
qb
；
[0042]
c)将a)和b)中的t
sb
、t
qb
、f
sb
、f
qb
共同组成元测试任务w；
[0043]
d)重复步骤a)至c)p次，获取p个元测试任务。
[0044]
本技术方案中，通过多次的元测试任务，使得使得并行网络在已有上面任务处理抽象概念的基础上通过元训练任务学会如何解决当前的小样本问题，为后续的测试提供好的基础。
[0045]
作为本发明的进一步改进，所述信号识别前，还包括测试任务的获取，具体为：
[0046]
a)在预处理后的元测试db时域信号的三类信号中，每类信号中，均随机抽取q1,q2,q3三条信号作为时域支持集t
sc
，将预处理后的测试集中三类信号的所有时域信号作为时域查询集t
qc
；
[0047]
b)在预处理后的元测试集db频域信号的bfsk，qpsk，qfsk这三类信号中，分别选择a)中时域q1,q2,q3三条信号对应的频域信号u1,u2,u3，将c1,c2,c3作为频域支持集f
sc
，将预处理后的测试集dc中bfsk，qpsk，qfsk这三类信号的所有频域信号作为频域查询集f
qc
；
[0048]
c)将a)和b)中的t
sc
、t
qc
、f
sc
、f
qc
共同组成测试任务r。
[0049]
本技术方案中，对于测试任务也进行多册的预处理等，确保方案最优化。
附图说明
[0050]
图1为本发明提供的基于并行原型网络的小样本雷达辐射源信号识别方法的流程图；
[0051]
图2为本发明提供的实施例3中的流程图；
[0052]
图3为本发明提供的并行原型网络的结构图。
具体实施方式
[0053]
下面结合附图所示的各实施方式对本发明进行详细说明，但应当说明的是，这些
实施方式并非对本发明的限制，本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代，均属于本发明的保护范围之内。
[0054]
在本实施例的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明创造和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明创造的限制
[0055]
此外，术语“第一”、“第二”、“第三”等仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明创造的描述中，除非另有说明，“多个”的含义是两个或两个以上。术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以通过具体情况理解上述术语在本发明创造中的具体含义。
[0056]
实施例1
[0057]
参照附图1所示，本实施例中，基于并行原型网络的小样本雷达辐射源信号识别方法，包括以下步骤：
[0058]
数据集的获得：利用仿生软件，仿生形成时频数据形式的已知雷达辐射信号数据集、样本雷达辐射信号数据集和待识别雷达射频信号数据集；所述已知雷达辐射信号数据集包括cw、lfm、nlfm和bpsk四类信号数据集，所述样本雷达辐射信号数据集和待识别雷达射频信号数据集均包括bfsk、qpsk和qfsk三类信号数据集；
[0059]
元训练集、元测试集以及测试集的生成：分别将获取的已知雷达辐射信号数据集、样本雷达辐射信号数据集和待识别雷达射频信号数据集进行傅里叶变换，得到时域数据集和频域数据集共同构成的元训练集da、元测试集db以及测试集dc；
[0060]
初始并行原型网络的生成：分别依次构建时域和频域的特征提取模块和度量模块，然后将时域和频域的度量模块组合，形成识别模块，以识别模块、2个所述度量模块以及2个所述特征提取模块组合，形成初始并行原型网络；
[0061]
并行原型网络的训练和测试：将元训练集在初始并行原型网络中进行元训练任务的训练，得到第一权重，然后将元测试集在载有第一权重的训练并行原型网络上进行元测试集的继续训练，得到第二权重，用第二权重覆盖第一权重，得到训练好的并行原型网络；
[0062]
信号识别：将测试集中获取的测试任务，输入到载有第二权重的并行原型网络中，获取得到雷达信号的识别测试结果。
[0063]
本实施例中，首先，通过元训练集(即包含雷达辐射信号数据集及其频域数据)训练网络，使网络获得如何解决每一个任务的抽象概念；再通过元测试集(即样本雷达辐射信号数据集及其频域数据)训练网络，使网络在上面已经学会如何解决每一个任务的基础上，学会如何解决当前设置的小样本问题(即识别qpsk,qfsk,bfsk三类信号)，实现对小样本的识别。
[0064]
本实施例中，其次，由于cw,lfm,nlfm,bpsk,bfsk,qpsk,qfsk这7类信号属于比较常用的信号类型，故本实验选取这7类信号作为示例说明原理。
[0065]
这7类信号中可以任意选取几类信号作为“具有大量已知信号”的数据集，几类信号作为“小样本待识别”数据集，两个数据集不能有类别交集，为了模拟大量已知信号样本和小样本。
[0066]
本实施例选取bfsk,qpsk,qfsk三类信号作为小样本信号，生成少量的含有这三类信号的数据集作为现在获取的少量标记样本信号(记为b)，同时也生成含有大量这三类信号的数据集作为待识别数据集用以验证本方法的测试识别率(记为c)，并选取剩下的四类作为大量已知信号(记为a)。
[0067]
通过本技术方案，用这四类大量已知信号集a和三类小样本信号集b训练网络，使网络最终能够识别这三类小样本信号，最后利用信号集c验证本文的小样本识别方法识别率。
[0068]
进一步地，本实施例中，已知雷达辐射信号数据集(记为a)是为了模拟已有的有标注的大量数据，样本雷达辐射信号数据集(记为b)模拟已有的有标注的少量数据，a和b类别完全不同，待识别雷达信号数据集(记为c)模拟与b同类的未标注数据。
[0069]
小样本识别的目的和方法就是在仅有少量标注样本(即b)的情况下，借助大量其他类型的有标注数据(即a)学会如何识别b这些类型的样本，以便在下次采集到这些类型的未标记样本(即c)时可以识别并对其标注。因此，a数量大，b数量小。
[0070]
本实施例中，最后，时域和频域信号分别输入各自的特征提取模块，不仅可以更好的提取样本的时域特征和频域特征，并能充分利用小样本的时域和频域信息，为识别提供依据。进而，通过时域和频域度量模块分别计算时域和频域上待识别信号特征与已知信号特征原型的距离作为判断，并通过识别模块平均时域和频域的判断做出最终识别，由于综合了时域和频域的判断，并避免了在前期对时、频特征直接结合进行判断可能造成的特征冗余问题。
[0071]
本发明的技术思路是：对获取的雷达辐射源时域信号进行快速傅里叶变换，并将其时域和频域信号共同送入并行原型网络进行训练，最终以时域和频域输出结果的平均作为最终识别结果，有效利用信号的多个特征对信号进行分辨识别的同时避免直接将时频特征进行结合可能造成的特征冗余，并且，原型网络通过已知信号样本的特征原型与未知样本的特征进行距离计算，可以有效地识别未知样本的信号类型。此外，以任务训练方式训练并行原型网络，通过不断适应每个具体任务，使并行原型网络在该训练方式下获取如何解决每一个任务的抽象概念，解决小样本问题。
[0072]
实施例2
[0073]
参照附图1-2所示，本实施例中，针对实施例1的方案，进行细化描述。
[0074]
本发明的实现方案包括如下：
[0075]
1)用matlab软件仿真生成雷达信号的数据集：
[0076]
1a)生成元训练集da，模拟拥有的大量已知雷达辐射源信号(此处，matlab仿真大量信号，作为已知标签的大量信号)，该数据集包括cw，lfm，nlfm，bpsk共4类信号，每一类信号在snr＝-4db下产生2000个信号作为时域数据集，并对该数据集进行快速傅里叶变换，获得信号的频域数据集，时域数据集和频域数据集共同构成本实验的元训练集da；
[0077]
1b)生成元测试集db，模拟仅有的少量已知雷达辐射源信号，该数据集包括bfsk，qpsk，qfsk共3类信号，每一类信号在snr＝-4db下产生10个信号作为时域数据集，并对该数
据集进行快速傅里叶变换，获得信号的频域数据集，时域数据集和频域数据集共同构成本实验的元测试集db；
[0078]
1c)生成测试集dc，模拟待识别的雷达辐射源信号，该数据集包括bfsk，qpsk，qfsk共3类信号，每一类信号在snr＝-4db下产生1000个信号作为时域数据集，并对该数据集进行快速傅里叶变换，获得信号的频域数据集，时域数据集和频域数据集共同构成本实验的测试集dc；
[0079]
2)通过z-score标准化方法对元训练集da、元测试集db及测试集dc的数据进行预处理，获得均值为0、标准差为1的预处理后的元训练集da、元测试集db及测试集dc，用以保证不同维度的数据具有相同的分布规模，更利于网络优化，缩短训练时间；
[0080]
3)在预处理后的元训练集da中获取d个元训练任务u；
[0081]
4)构建并行原型网络：
[0082]
4a)搭建特征提取模块：特征提取模块包含时域特征提取模块和频域特征提取模块，构建方法如下：
[0083]
4a1)时域特征提取模块具有四个卷积单元和一个平坦层，其中第一、二、三个卷积单元构造相同(第四个卷积单元通过减少滤波器个数浓缩信号特征，使其更简洁)，卷积单元组成均为卷积层、激活层、批量归一化层和最大池化层，其中，卷积层的卷积核大小为16，滤波器个数为64，步长为1；池化层池化大小为3，步长为3；激活层选择激活函数为relu；第四个卷积单元由卷积层、激活层和批量归一化层组成，其中，卷积层的卷积核大小为16，滤波器个数为16，步长为1；激活层选择激活函数为relu；
[0084]
4a2)频域特征提取模块搭建方法同4a1)；
[0085]
此处，卷积单元利用卷积层充分提取信号的局部特征；利用激活层将特征映射到高维的非线性区间，解决线性模型所不能解决的问题；利用批量归一化层加快网络的训练和收敛速度，降低梯度爆炸和梯度消失情况的发生；利用池化层在保留主要特征的同时减少参数，即对特征进行降维，提高模型泛化能力。卷积单元这样设置，既保证了输入信号特征能被有效提取，也能防止网络的过拟合，并且提高训练速度。
[0086]
4b)搭建度量模块：度量模块包含时域度量模块和频域度量模块，构建方法如下：
[0087]
4b1)时域度量模块选择欧几里得距离作为度量方法，并具有一个激活函数为softmax的激活层；
[0088]
4b2)频域度量模块搭建方法同4b1)；
[0089]
4c)搭建识别模块，该模块由平均值函数和最大值函数构成，通过对时域和频域度量模块的输出进行平均计算，综合时域和频域量模块的输出，并取对时域和频域度量模块的输出进行平均计算之后的最大值索引作为最后的识别结果；
[0090]
4d)对4a)、4b)、4c)共同组成并行原型网络；
[0091]
5)对4)搭建好的并行原型网络进行训练，训练过程如下：
[0092]
将3)中的任务输入4)搭建好的并行原型网络，对并行原型网络进行训练，直到识别正确率达到预期，停止训练，保存并行原型网络权重m；
[0093]
6)对5)训练好的网络模型继续训练，训练过程如下：
[0094]
6a)在预处理后的元测试集db中获取p个元测试任务w；
[0095]
6b)在4)的并行原型网络中载入5)获得的并行原型网络权重m，将6a)中的任务输
入并行原型网络，对并行原型网络进行训练，直到识别正确率达到预期，停止训练，保存训练好的并行原型网络权重k；
[0096]
7)载入6)获得的并行原型网络权重k，进行测试，测试过程如下：
[0097]
7a)在预处理后的元测试集db和测试集dc中获取测试任务r；
[0098]
7b)将7a)中的任务r输入载入并行原型网络权重k的并行原型网络，获得辐射源信号识别测试结果。
[0099]
本发明具有如下优点：
[0100]
1)本发明将原型网络应用于小样本雷达辐射源识别，解决在雷达辐射元信号识别领域的小样本识别应用较少的问题；
[0101]
2)本发明改进原型网络，构建并行原型网络，分别输入时域信号和频域信号，综合高效利用小样本特征，实现雷达辐射源小样本识别并提高识别率；
[0102]
3)本发明通过平均时域和频域识别结果作为最终识别结果，而非进行时频特征结合，避免特征融合可能造成的冗余问题。
[0103]
实施例3
[0104]
本实施例中，以具体参数为主进行介绍。
[0105]
参照图1-3，本实施例的实现包括如下步骤：
[0106]
步骤1：用matlab软件仿真生成雷达信号数据集。
[0107]
1.1)生成元训练集da，模拟拥有的大量已知雷达辐射源信号：
[0108]
1.1.1)生成元训练集da的时域数据集，该数据集包括cw，lfm，nlfm，bpsk共4种类信号，每一类在-4db条件下产生2000个时域信号，则共8000个信号作为时域数据集进行实验，信号设置如表1所示：
[0109]
表1 da信号参数设置
[0110][0111]
[0112]
统一采样频率fs为1ghz，脉宽pw设置为1.025us，采样点数为1025个点，nlfm信号的调制信号频率f设置为10mhz。为使仿真信号更接近原始雷达辐射源信号，载频fc、信号带宽bw、编码方式等信号调制参数不取固定参数，其中，四类信号的载频fc均在180～220mhz间随机取值，lfm信号带宽bw在9～11mhz间随机取值，bpsk信号的巴克码位数在5，7，11，13四种选择中随机确定；参照表1中，载频fc、信号带宽bw、编码方式等信号调制参数不取固定参数，以使仿真信号更接近真实雷达辐射源信号。
[0113]
本实施例中，信噪比选择-4db模拟复杂的信号环境。小样本识别的目的和方法就是在仅有少量标注样本(即b)的情况下，借助大量其他类型的有标注数据(即a)学会如何识别b这些类型的样本，以便在下次采集到这些类型的未标记样本(即c)时可以识别并对其标注。因此，a数量大，b数量小，a每类信号选取2000个模拟大量已知，实际应用时并不规定必须2000个，只要是大量就可以，此处选择2000做实验仿真；b每类信号选取10个模拟少量标注，实际应用时并不规定必须是10个，只是本实验选取10个作为少量，进行试验；c每类选取1000模拟与b同类的未标注样本，实际应用时并不是必须1000个，此处选择1000个做实验仿真，验证本方法的识别率。
[0114]
1.1.2)对1.1.1)产生的元训练集da的时域数据集进行快速傅里叶变换，得到具有8000个信号的频域数据集；
[0115]
1.1.3)1.1.1)的时域数据集和1.1.2)的频域数据集组成本实验的元训练集da；
[0116]
1.2)生成元测试集db，模拟仅有的少量已知雷达辐射源信号：
[0117]
1.2.1)生成元测试集db的时域数据集，该数据集包括bfsk，qpsk，qfsk共3类信号，每一类在-4db条件下产生10个时域信号，则共30个信号作为时域数据集进行实验，信号设置如表2所示：
[0118]
表2 db信号参数设置
[0119]
[0120][0121]
统一采样频率fs为1ghz，脉宽pw设置为1.025us，采样点数为1025个点，bpsk信号的两个载频f
c1
、f
c2
分别在80～120mhz、130～170mhz间随机取值，巴克码位数在5，7，11，13四种选择中随机确定；qpsk信号载频fc在180～220mhz间随机取值，采用16位frank码；qfsk四个载频f
c1
、f
c2
、f
c3
、f
c4
分别在80～120mhz、130～170mhz、180～220mhz、230～270mhz间随机取值，采用16位frank码；
[0122]
1.2.2)对1.2.1)产生的元测试集db的时域数据集进行快速傅里叶变换，得到具有30个信号的频域数据集；
[0123]
1.2.3)1.2.1)的时域数据集和1.2.2)的频域数据集组成本实验的元测试集db；
[0124]
1.3)生成测试集dc，模拟待识别的雷达辐射源信号：
[0125]
数据集包括bfsk，qpsk，qfsk共3类信号，每一类信号在snr＝-4db下产生1000个信号作为时域数据集，则共3000个信号作为时域数据集进行实验，信号设置如表2所示；
[0126]
1.3.2)对1.3.1)产生的测试集dc的时域数据集进行快速傅里叶变换，得到具有3000个信号的频域数据集；
[0127]
1.3.3)1.3.1)的时域数据集和1.3.2)的频域数据集组成本实验的测试集dc；
[0128]
步骤2：对元训练集da、元测试集db及测试集dc进行预处理。
[0129]
通过z-score标准化方法对元训练集da、元测试集db及测试集dc进行预处理，获得均值为0、标准差为1的预处理后的元训练集da、元测试集db及测试集dc，用以保证不同维度的数据具有相同的分布规模，更利于网络优化，缩短训练时间；
[0130]
其中，z-score标准化方法如下：
[0131]
设某条信号的序列为x＝[x1,x2,x3,
···
,xn]，其中n为采样点数。以第i个采样点为例，标准化之后的采样值为：
[0132][0133]
其中，xi是信号序列x的第i个采样点，是该条信号序列的平均值，s是该条信号序列的标准差，其计算表达式如下：
[0134][0135][0136]
本实施例中，通过标准化的处理，用以保证不同维度的数据具有相同的分布规模，更利于网络优化，缩短训练时间；由于本实施例中，有多种类型的数据，进而通过这样的处
理，使其更加统一化。
[0137]
步骤3：在预处理后的元训练集da中获取d个元训练任务u。
[0138]
3.1)在预处理后的元训练集da中的cw，lfm，nlfm，bpsk四类信号中随机抽取a1,a2,a3三个类型；
[0139]
3.2)在预处理后的元训练集da时域信号的a1,a2,a3这三类信号中，分别随机抽取a1,a2,a3,
···
,a8八条信号，将a1,a2,a3作为时域支持集t
sa
，a4,a5,
···
a8作为时域查询集t
qa
；
[0140]
3.3)在预处理后的元训练集da频域信号的a1,a2,a3这三类信号中，分别选择3.2)中时域a1,a2,a3,
···
,a8八条信号对应的频域信号b1,b2,b3,
···
,b8，将b1,b2,b3作为频域支持集f
sa
，b4,b5,
···
b8作为频域查询集f
qa
；
[0141]
3.4)3.2)和3.3)中的t
sa
、t
qa
、f
sa
、f
qa
共同组成元训练任务u；
[0142]
3.5)重复步骤3.1)至3.4)d次，；
[0143]
步骤4：构建并行原型网络。
[0144]
4.1)搭建特征提取模块：特征提取模块包含时域特征提取模块和频域特征提取模块，构建方法如下：
[0145]
4.1.1)时域特征提取模块具有四个卷积单元和一个平坦层，输入第i个训练任务u的时域支持集t
sa
和时域查询集t
qa
，输出时域支持集t
sa
和时域查询集t
qa
的信号特征；
[0146]
其中第一、二、三个卷积单元构造相同，卷积单元组成均为卷积层、激活层、批量归一化层和最大池化层，其中，卷积层的卷积核大小为16，滤波器个数为64，步伐为1；池化层池化大小为3，步伐为3；激活层选择激活函数为relu；第四个卷积单元由卷积层、激活层和批量归一化层组成，其中，卷积层的卷积核大小为16，滤波器个数为16，步伐为1；激活层选择激活函数为relu；
[0147]
4.1.2)频域特征提取模块搭建方法同4.1.1)，输入第i个训练任务u的频域支持集f
sa
和频域查询集f
qa
，输出频域支持集f
sa
和频域查询集f
qa
的信号特征；
[0148]
此处，时域和频域度量模块搭建方法完全相同，只是输入输出不同：具体地，时域度量模块输入时域支持集和查询集信号特征，输出时域查询集每个信号所属时域支持集每个信号类别的概率；频域度量模块输入频域支持集和查询集的信号特征，输出频域查询集的每个信号所属频域支持集每个信号类别的概率。
[0149]
4.2)搭建度量模块：度量模块包含时域度量模块和频域度量模块，构建方法如下：
[0150]
4.2.1)时域度量模块选择欧几里得距离作为度量方法，具有一个激活函数为softmax的激活层，输入时域支持集t
sa
和时域查询集t
qa
的信号特征，输出时域查询集t
qa
的每个信号所属时域支持集t
sa
每个信号类别的概率；
[0151]
4.2.2)频域度量模块搭建方法同4.2.1)，输入频域支持集f
sa
和频域查询集f
qa
的信号特征，输出频域查询集f
qa
的每个信号所属频域支持集f
sa
每个信号类别的概率；
[0152]
此处具体是，通过平均时域与频域识别类别概率，并选取查询集每条信号所述支持集各类别概率中的最大概率作为查询集该条信号的所属类别进行输出。
[0153]
本实施例中，通过平均值函数对查询集各信号所属的类别概率进行时、频域判断结果综合，得到综合判断的查询集各信号所属支持集各类别的概率，具体如下：
[0154]
[0155]
通过最大值函数选取查询集每条信号所述支持集各类别概率中的最大概率作为查询集该条信号的所属类别。
[0156]
4.3)搭建识别模块，该模块由平均值函数和最大值函数构成，通过对时域和频域度量模块的输出进行平均计算，平均时域和频域的识别结果，综合时域和频域的判断，选择平均结果中概率最大的类别作为查询集的最终类别，而非进行时频特征结合，避免特征融合可能造成的冗余问题，计算识别准确率；
[0157]
4.4)4.1)、4.2)、4.3)共同组成并行原型网络；
[0158]
步骤5：训练并行原型网络。
[0159]
将步骤3中的任务输入步骤4搭建好的并行原型网络，训练并行原型网络，直到识别正确率达到预期，停止训练，保存并行原型网络权重m；
[0160]
步骤6：对步骤4的并行原型网络载入并行原型网络权重m，继续训练并行原型网络。
[0161]
6.1)在预处理后的元测试集中获取p个元测试任务w，获取方法如下：
[0162]
6.1.1)在预处理后的元测试集db时域信号的bfsk，qpsk，qfsk这三类信号中，分别随机抽取e1,e2,e3,
···
,e8八条信号，将e1,e2,e3作为时域支持集t
sb
，e4,e5,
···
e8作为时域查询集t
qb
；
[0163]
6.1.2)在预处理后的元测试集db频域信号的bfsk，qpsk，qfsk这三类信号中，分别选择6.1.1)中时域e1,e2,e3,
···
,e8八条信号对应的频域信号c1,c2,c3,
···
,c8，将c1,c2,c3作为频域支持集f
sb
，c4,c5,
···
c8作为频域查询集f
qb
；
[0164]
6.1.3)6.1.1)和6.1.2)中的t
sb
、t
qb
、f
sb
、f
qb
共同组成元测试任务w；
[0165]
6.1.4)重复步骤6.1.1)至6.1.3)p次，获取p个元测试任务；
[0166]
6.2)在步骤4的并行原型网络中载入步骤5获得的并行原型网络权重m，将6.1)中的任务输入并行原型网络，训练并行原型网络，直到识别正确率达到预期，停止训练，保存训练好的并行原型网络权重k；
[0167]
步骤7：对步骤4的并行原型网络载入并行原型网络权重k，进行测试。
[0168]
7.1)获取测试任务r，获取方式如下：
[0169]
7.1.1)在预处理后的元测试集db时域信号的bfsk，qpsk，qfsk这三类信号中，分别随机抽取q1,q2,q3三条信号作为时域支持集t
sc
，将预处理后的测试集dc中bfsk，qpsk，qfsk这三类信号的所有时域信号作为时域查询集t
qc
；
[0170]
7.1.2)在预处理后的元测试集db频域信号的bfsk，qpsk，qfsk这三类信号中，分别选择7.1.1)中时域q1,q2,q3三条信号对应的频域信号u1,u2,u3，将c1,c2,c3作为频域支持集f
sc
，将预处理后的测试集dc中bfsk，qpsk，qfsk这三类信号的所有频域信号作为频域查询集f
qc
；
[0171]
7.1.3)7.1.1)和7.1.2)中的t
sc
、t
qc
、f
sc
、f
qc
共同组成测试任务r；
[0172]
7.2)将7.1)中的任务r输入载入并行原型网络权重k的并行原型网络，获得辐射源信号识别测试结果。
[0173]
本发明能充分利用并有效提取雷达辐射源小样本信号特征，并有较好的识别率，可用于雷达辐射源小样本信号的识别。
[0174]
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说
明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
[0175]
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将权利要求落在等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
[0176]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。