一种基于贝塞尔曲线拟合雷达辐射源个体识别方法及系统技术方案

本发明专利技术涉及一种基于贝塞尔曲线拟合雷达辐射源个体识别方法及系统，技术要点为：首先提取雷达脉冲信号UPMOP特征，再用贝塞尔曲线对该特征进行拟合，最后取拟合后曲线的控制点作为个体特征输入神经网络分类器。从雷达辐射源个体识别面临的特征维度高、识别效率低等实际出发，选取了稳定性较好的UPMOP作为个体特征，利用贝塞尔曲线实现对UPMOP的拟合，提取拟合后的控制点作为新的个体识别特征向量，实现了特征维度的大幅降低。解决了现有技术中存在的建模困难、指纹特征不易测量、运算复杂度高的问题。

【技术实现步骤摘要】
一种基于贝塞尔曲线拟合雷达辐射源个体识别方法及系统
本专利技术属于信号处理
，具体涉及一种基于贝塞尔曲线拟合雷达辐射源个体识别方法及系统。
技术介绍
辐射源个体识别的核心是个体特征提取，通过对接收信号进行测量，将辐射源唯一电磁特征与辐射源个体关联起来(指纹识别技术)，对采集到的信号进行特征(指纹)提取，完成特征选择和融合后送入分类器。对信号的分析包括信号暂态特征分析和信号稳态特征分析，信号稳态特征分析是指辐射源在稳定工作时，不仅存在为有效可靠传递信息而人为添加的调幅、调频以及编码灯有意调制特征，同时还有辐射源对信号的无意调制特征。无意调制特征的出现主要是由于辐射源内部频率源不稳定、辐射源内部诸多器件的非线性特性和器件加工过程中工艺缺陷造成的器件特性和技术指标的离散型引起的。早期由Dr.Lawrence等人提出的利用脉内无意调制(UnintentionalModulationonPulse,UMOP)作为识别特征的方法取得了一定的效果。后来，基于域变换的很多优秀算法(包络分析、双谱、循环谱和小波变换等)被用来提取个体特征，这些算法虽对冗余特征进行了相应的优化，但存在个体特征维度高以及以仿真实验为主，未针对实际信号进行测试等缺陷。许丹、黄渊凌等人另辟蹊径，从指纹产生的机理出发，对发射机频率源畸变、功放的非线性等特性进行分析、建模，所提取的个体特征达到了较好的识别效果，但是，运算过程中存在较多的个体特征，运算量大且效率较低，存在运算复杂度高以及缺乏实测数据验证等问题。
技术实现思路
本专利技术的目的在于提供一种基于贝塞尔曲线拟合雷达辐射源个体识别方法及系统，用于解决现有技术中存在的运算复杂度高的问题。为解决上述技术问题，本专利技术的技术方案为：本专利技术提供了一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，包括以下步骤：1)提取雷达脉冲信号的脉内无意调制特征D，特征D的表达式为：D＝{y(n)|n∈[0,N-1]}其中N为脉冲持续时间内的采样点数，y(n)为脉内无意调相特征；2)利用贝塞尔曲线对特征D所在的曲线进行拟合；3)提取步骤2)中拟合后曲线的控制点，将控制点作为雷达脉冲信号的个体特征，对其进行识别。有益效果：本专利技术提供了一种基于贝塞尔曲线拟合的雷达辐射源个体特征提取方法。首先提取雷达脉冲信号UPMOP特征，再用贝塞尔曲线对该特征进行拟合，最后取拟合后曲线的控制点作为个体特征输入神经网络分类器。从雷达辐射源个体识别面临的特征维度高、识别效率低等实际出发，选取了稳定性较好的UPMOP作为个体特征，利用贝塞尔曲线实现对UPMOP的拟合，提取拟合后的控制点作为新的个体识别特征向量，实现了特征维度的大幅降低。进一步的，步骤2)中目标函数为：其中F(x)为残差的F范数，称为伯恩斯坦多项式，Pii∈[0,d]为控制点，为二项式系数。进一步的，步骤3)中提取控制点的方法为高斯牛顿法或LM算法。进一步的，通过神经网络分类器对所述控制点进行识别。本专利技术还提供了一种基于贝塞尔曲线拟合雷达辐射源个体识别系统，包括处理器和储存器，所述处理器执行存储在储存器中的如下方法指令：1)提取雷达脉冲信号的脉内无意调制特征D，特征D的表达式为：D＝{y(n)|n∈[0,N-1]}其中N为脉冲持续时间内的采样点数，y(n)为脉内无意调相特征；2)利用贝塞尔曲线对特征D所在的曲线进行拟合；3)提取步骤2)中拟合后曲线的控制点，将控制点作为雷达脉冲信号的个体特征，对其进行识别。有益效果：该系统包括处理器和储存器，处理器发出指令首先提取雷达脉冲信号UPMOP特征，再用贝塞尔曲线对该特征进行拟合，最后取拟合后曲线的控制点作为个体特征输入神经网络分类器。从雷达辐射源个体识别面临的特征维度高、识别效率低等实际出发，选取了稳定性较好的UPMOP作为个体特征，利用贝塞尔曲线实现对UPMOP的拟合，提取拟合后的控制点作为新的个体识别特征向量，实现了特征维度的大幅降低。进一步的，步骤2)中目标函数为：其中F(x)为残差的F范数，称为伯恩斯坦多项式，Pii∈[0,d]为控制点，为二项式系数。进一步的，步骤3)中提取控制点的方法为高斯牛顿法或LM算法。进一步的，通过神经网络分类器对所述控制点进行识别。附图说明图1是3阶贝塞尔曲线示意图；图2是基于贝塞尔曲线拟合的特征提取流程图；图3是拟合平均相对误差随阶数变化示意图；图4是3阶贝塞尔曲线拟合示意图；图5a是6部发射机中前3部发射机UPMOP拟合结果示意图；图5b是6部发射机中后3部发射机UPMOP拟合结果示意图；图6是5部发射机控制点聚类示意图。具体实施方式为使本专利技术的目的、技术方案及优点更加清楚，下面结合附图及实施例，对本专利技术作进一步的详细说明。接收到的辐射源信号一般模型为：其中，A为脉冲幅度，是初始相位，fc为信号载频，w(n)是均值为零方差为σ2的高斯白噪声，N脉冲持续时间内的采样点数，是相位调制函数，常用的有以下几种：(1)单频脉冲，(2)线性调频脉冲，K为调制率；(3)相位编码，如二相编码雷达，在子脉冲持续时间内，c(n)＝0或1。ΔA(n)和分别代表包络和相位的无意调制，对于实际的雷达发射机而言，模拟电路的不理想和环境的差异性总是不可避免的，导致ΔA(n)和总是存在，通过对这两个参数的估计可以反映出发射机的个体差异性。贝塞尔曲线表达式为其中称为伯恩斯坦多项式，Pii∈[0,d]为控制点，为二项式系数。贝塞尔曲线由线段和控制点组成，给定二维平面内线段的端点和控制点，则可唯一确定一条光滑曲线，反之，指定平面内光滑曲线及阶数则可以计算出控制点的唯一解。如图1所示，为3阶贝塞尔曲线，其中P1、P2为曲线的控制点，用于描述曲线的弯曲程度。UPMOP(脉内无意调相)特征可以看做是长度与脉冲宽度相当的光滑曲线，显然将该曲线作为特征向量时维度随着脉冲宽度增加而增大；与此同时，贝塞尔曲线在低阶的情况下就能够实现对曲线的良好拟合，以拟合后曲线的控制点作为个体特征向量可以实现大幅度降维，因此，本专利技术选取贝塞尔曲线来拟合UPMOP。下面为具体实施例如图2所示，首先对脉冲信号的描述字进行测量，其次对每个脉冲信号进行时相估计，提取其UPMOP特征。然后构造残差的F范数，与初始参数μ、增长因子β、贝塞尔曲线迭代停止条件ε和最大迭代次数kmax输入LM算法计算过程，当达到停止条件或最大迭代次数时输出控制点Pi。具体为：定义脉内无意调相UPMOP为D，表达式为D＝{y(n)|n∈[0,N-1]}其中N为脉冲持续时间内的采样点数，y(n)为脉内无意调相特征，即用y(n)表示上述x(n)中的然后构造目标函数，定义F(x)为残差的F范数，即由上式可以看出，贝塞尔曲线拟合UPMOP，可以归结为非线性最小二乘问题。针对上式，若采用高斯牛顿法，可得hgn＝-(JTJ)-1JTFx＝x+hgn＝x-(JTJ)-1JTF其中，J为F(x)的雅克比矩阵，hgn为高斯牛顿的方向。由于高斯牛顿法在JTJ奇异或接近奇异的情况下，求(JTJ)-1会遇到不稳定或偏差较大的问题，有时甚至不可逆。为此，LM算法对其进行了修正，在高斯牛顿的基础上，把一个正定对角矩阵μI加到JTJ上去，改变原矩阵的特征值结构，使其变成条件数较好的对称正定矩阵，从而使得运算过程得以继续进行本文档来自技高网...

【技术保护点】
1.一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，其特征在于，包括以下步骤：1)提取雷达脉冲信号的脉内无意调制特征D，特征D的表达式为：D＝{y(n)|n∈[0,N‑1]}其中N为脉冲持续时间内的采样点数，y(n)为脉内无意调相特征；2)利用贝塞尔曲线对特征D所在的曲线进行拟合；3)提取步骤2)中拟合后曲线的控制点，将控制点作为雷达脉冲信号的个体特征，对其进行识别。

【技术特征摘要】
1.一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，其特征在于，包括以下步骤：1)提取雷达脉冲信号的脉内无意调制特征D，特征D的表达式为：D＝{y(n)|n∈[0,N-1]}其中N为脉冲持续时间内的采样点数，y(n)为脉内无意调相特征；2)利用贝塞尔曲线对特征D所在的曲线进行拟合；3)提取步骤2)中拟合后曲线的控制点，将控制点作为雷达脉冲信号的个体特征，对其进行识别。2.根据权利要求1所述的一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，其特征在于，步骤2)中拟合通过构造目标函数来实现的，所构造的目标函数为：其中F(x)为残差的F范数，称为伯恩斯坦多项式，Pii∈[0,d]为控制点，为二项式系数。3.根据权利要求1所述的一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，其特征在于，步骤3)中提取控制点的方法为高斯牛顿法或LM算法。4.根据权利要求1或2或3所述的一种基于贝塞尔曲线拟合雷达辐射源个体识别方法，其特征在于，通过神经网络分类器对所述控制点进行识别。5.一...

【专利技术属性】
技术研发人员：邢小鹏，黄洁，陈世文，赵拥军，
申请(专利权)人：中国人民解放军战略支援部队信息工程大学，
类型：发明
国别省市：河南,41