技术领域
本发明属于跳频分选技术领域，具体涉及一种基于跳频聚类分选结果的错误检测和修正方法。
背景技术
在现代信息技术高速发展的条件下，为了取得战争的胜利，各国军队都在大量使用常规或非常规的通信设备。由于无线通信所具有的开放性和无界性的特点，战场通信能够为己方带来即时的信息交互。1995年，美国国防部提出信息战的概念，随着信号处理技术的迅速发展，信息战已经作为一种作战方式融入到现代战争中，成为影响战场局势的重要因素。通信对抗是战场上有效打击对手的一种信息化手段，研究的内容是在非协作通信条件下对电磁环境中的隐蔽信号检测和估计，进而获取信息，对敌方的协作通信实施干扰进而完成对电子通信载体进行精准打击。因此，为了保障己方通信的安全性与隐蔽性，低截获、抗干扰等技术手段被广泛应用于通信电子对抗之中，而跳频通信便是最主要的常用技术之一。跳频分选技术是现代通信对抗中的一项关键技术，具有重大的战略意义。
在传统跳频分选中，会对侦收到的跳频信号进行参数估计生成每一跳信号的Hop描述字(Hop Descirption Word，HDW)，包括跳周期、跳时和到达角等，并依据相似性使用聚类算法进行盲分选，将属于同一跳信号的Hop筛选分类，分选成单独完整的跳频信号。
当不同跳频信号的HDW数据集在样本空间发生混叠时，聚类分选会产生错误，因此依据跳频信号在时间上的连续性和单一性，提出了基于序列差分的跳频聚类分选错误检测和修正方法，能够一定程度地提升分选准确率。
发明内容
本发明的目的是为了能够对跳频聚类分选结果进行检测和修正，提高准确率，传统跳频分选会依据特征参数的相似性使用聚类算法进行分选，并且以聚类的结果作为分选结果，当特征参数集在样本空间中产生重叠时，可能会产生违背跳频信号时间特性的分选错误问题，提出了一种基于序列差分的跳频聚类分选结果的错误检测和修正方法，能够提高跳频信号分选的准确率。
本发明技术方案为：基于序列差分的跳频聚类分选结果的错误检测和修正方法，包括以下步骤：
S1、搭建数据采集平台，并采集待分选的多个跳频信号测量参数组成的数据集；
S2、使用聚类算法进行分选，输出分选结果；
S3、根据分选后单个跳频信号的时频图，提取其中每一跳的中心时刻，组成中心时刻序列；
S4、根据中心时刻序列作出曲线，将序列依次差分并与设定的检测门限进行比较，并判断错误类型；
S5、将不同跳频信号之间发生错误的时间范围和类型进行匹配，若互相吻合则进行错误修正。
进一步地，步骤S2中当跳频信号的参数数据集在样本空间中产生交叠时，会发生分选错误，但由于跳频信号的时间连续性选择和单一性原则，分选结果可能会产生违背以上两种原则的错误，因而需要对其检测和修正。
进一步地，步骤S3中提取中心时刻包括以下分步骤：
S31、对于时频变换后的时频矩阵TF
S32、对经过细化处理的时频矩阵数据中TF′
得每个频率点对应的峰值位置，得到峰值矩阵Peak
S33、对于每一个时间点沿着频率轴利用二次差分法获得每个时间点对应的峰值位置，获得峰值矩阵Peak
S34、将Peak
S35、对峰值进行修正其对应的时间轴坐标即为跳中心时刻(Center Time，CT)。
进一步地，步骤S4中对中心序列的拟合采用更合理的曲线拟合，理想状况下，隶属于同一跳频信号的CT值序列可以在直角坐标系中连成一条直线以CT值的序号数作为x轴，以CT值的大小作为y轴，可以将直线表示出来：
y＝ax+b
式中，a为斜率对应的该序列的跳周期，b为截距，对应第一跳的持续时间。由于噪声以及测量误差的影响，各跳CT值构建的序列并非是一条完美的直线，CT测量值会分布在直线的周围，虽然对前后的CT序列值以此进行差分操作。
进一步地，步骤S4中序列差分、门限检测包括以下分步骤：
S41、根据分选后的某一个整体的跳频信号的CT序列，CT＝[CT
S42、对同一个跳频信号的CT序列进行差分操作获得ΔCT，其中ΔCT
S43、通过设定检测两个门限，判断是否产生I类错误和II类错误。若
Threshold
进一步地，步骤S32中的二次差分法为：
对于一个序列x(n)，n＝1，2，...，N，将x(n)相邻的两个数据的值进行比较得到一次差分的结果，设为y
进一步地，步骤S35中的峰值修正包括以下分步骤：
S351、搜寻峰值矩阵Peak(n，k)数值为1的坐标即峰值所在位置，得到(n1，k1)；
S352、以(n1，k1)为中心设置一个矩形的邻域Ω
S253、搜寻Peak(n，k)以(n1，k1)为中心邻域Ω
S354、令Peak(n′，k′)＝1，而原先的峰值全部都置0；
S355、对峰值矩阵Peak(n，k)所有峰值进行修正，得到修正后的峰值矩阵Peak′(n，k)；
S356、取每一个峰值位置的时间轴坐标，作为每跳信号的CT值。
本发明的有益效果是：
对于传统分选跳频聚类分选在数据集混叠情况下会发生分选错误，因此本发明提出了对于聚类分选后的跳频信号进行中心时刻序列提取并差分与门限进行检测，能够对违背跳频时间特性原则的错误进行检测和修正，提高分选准确率。
附图说明
图1所示为本发明实例提供的基于序列差分的跳频聚类分选结果的错误检测和修正方法流程图。
图2所示为本发明实例提供的两个跳频信号的参数测量数据集合和聚类分选结果图。
图3所示为本发明实例提供的分选错误结果类型示意图。
图4所示为本发明实例提供的二次差分法后提取峰值示意图。
图5所示为本发明实例提供的峰值修正后示意图。
图6所示为本发明实例提供的理想中心时刻曲线示意图。
图7所示为本发明实例提供的I类错误曲线和门限检测示意图。
图8所示为本发明实例提供的II类错误曲线和门限检测示意图。
具体实施方法
现在将参考附图来详细描述本发明的示例性实施方式。应当理解，附图中示出和描述的实施方式仅仅是示例性的，意在阐释本发明的原理和精神，而并非限制本发明的范围。
本发明实力提供一种一种基于跳频聚类分选结果的错误检测和修正方法，如图1所示，包括以下步骤：
S1、搭建数据采集平台，并采集生成待分选的多个跳频信号，如下表所示，测量参数组成的数据集；
S2、使用聚类算法进行分选，输出分选结果，如图2所示；
本发明实例中，步骤S2中跳频信号的参数数据集在样本空间中产生交叠时，会发生分选错误，但由于跳频信号的时间连续性选择和单一性原则，分选结果可能会产生违背以上两种原则的错误，如图3所示，因而需要对其检测和修正。
S3、根据分选后单个跳频信号的时频图，提取其中每一跳的中心时刻，组成中心时刻序列；
本发明实例中，步骤S3中提取中心时刻包括以下分步骤：
S31、对于时频变换后的时频矩阵TF
S32、对经过细化处理的时频矩阵数据中TF′
得每个频率点对应的峰值位置，得到峰值矩阵Peak
本发明实例中，步骤S32中的二次差分法为：
对于一个序列x(n)，n＝1，2，...，N，将x(n)相邻的两个数据的值进行比较得到一次差分的结果，设为y
S33、对于每一个时间点沿着频率轴利用二次差分法获得每个时间点对应的峰值位置，获得峰值矩阵Peak
S34、将Peak
S35、对峰值进行修正其对应的时间轴坐标即为跳中心时刻(Center Time，CT)，如图5所示。
本发明实例中，步骤S35中的峰值修正包括以下分步骤：
S351、搜寻峰值矩阵Peak(n，k)数值为1的坐标即峰值所在位置，得到(n1，k1)；
S352、以(n1，k1)为中心设置一个矩形的邻域Ω
S253、搜寻Peak(n，k)以(n1，k1)为中心邻域Ω
S354、令Peak(n′，k′)＝1，而原先的峰值全部都置0；
S355、对峰值矩阵Peak(n，k)所有峰值进行修正，得到修正后的峰值矩阵Peak′(n，k)；
S356、取每一个峰值位置的时间轴坐标，作为每跳信号的CT值。
S4、根据中心时刻序列作出曲线；将序列依次差分并与设定的检测门限进行比较，并判断错误类型；
本发明实例中，步骤S4中理想状况下，隶属于同一跳频信号的CT值序列可以在直角坐标系中连成一条直线以CT值的序号数作为x轴，以CT值的大小作为y轴，可以将直线表示出来，如图6所示：
y＝ax+b
式中，a为斜率对应的该序列的跳周期，b为截距，对应第一跳的持续时间。由于噪声以及测量误差的影响，各跳CT值构建的序列并非是一条完美的直线，CT测量值会分布在直线的周围，虽然对前后的CT序列值以此进行差分操作。
进一步地，本发明实例中，步骤S4中序列差分、门限检测包括以下分步骤：
S41、根据分选后的某一个整体的跳频信号的CT序列，CT＝[CT
S42、对同一个跳频信号的CT序列进行差分操作获得ΔCT，其中ΔCT
S43、通过设定检测两个门限，判断是否产生I类错误和II类错误。若
Threshold
S5、将不同跳频信号之间发生错误的时间范围和类型进行匹配，若互相吻合则进行错误修正，能够提升分选准确率，如下表所示。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
