基于特征提取的加权投票多分类器地面雷达目标自动识别方法与流程

1.本发明属于低分辨雷达数字信号处理技术领域，具体为一种基于特征提取的加权投票多分类器地面雷达目标自动识别方法。


背景技术：

2.随着现代战争信息集成的不断改进，雷达作为战场信息采集工具，其探测、测距等基本功能已经无法满足现代战场的需求。从雷达目标散射回波中提取其特征，并给出目标的性质、类型甚至具体模型成了迫切需要，雷达目标自动分类识别技术就由此诞生了。
3.根据分辨率性能，雷达可分为低分辨率窄带雷达和高分辨率宽带雷达。低分辨率窄带雷达由于成本低、实现方便等优点，应用仍然较为广泛，但由于其低分辨能力，回波中所包含的目标信息有限，故要从回波中提取相应的特征，必须采用各种信号处理方法对回波信息进行计算分析。雷达目标分类和识别过程有两个主要阶段：训练阶段和测试阶段。在训练阶段，从雷达回波数据中提取特征以进行分类器训练，而在测试阶段，将训练后的分类器用于目标分类测试。常用的目标分类器有朴素贝叶斯分类器、决策树分类器、lda分类器、k邻近分类器和支持向量机分类器等，大多分类设计都是基于单分类器进行的，但是基于单分类器，可能出现某些样本特征适合被特定分类器区分而不适合被另一些分类器区分的情况，导致一些特征利用率低、浪费计算量、识别准确率低。


技术实现要素：

4.本发明的目的在于提出一种基于特征提取的加权投票多分类器地面雷达目标自动识别方法。
5.实现本发明目的的技术方案为：一种基于特征提取的加权投票多分类器地面雷达目标自动识别方法，包含以下步骤：
6.步骤1，对时域回波信号进行时域特征提取，对下变频处理后的基带回波时域信号进行fft变换，进行频域特征提取；
7.步骤2，对提取后的特征进行5折交叉验证朴素贝叶斯分类器、决策树分类器、lda分类器、k邻近分类器和支持向量机分类器训练，并通过贝叶斯优化对支持向量机和k邻近分类器的超参数进行优化，得到最终的五个子分类器；
8.步骤3，根据混淆矩阵计算五个子分类器的精度权值，通过加权投票的方式得到最终的分类器；
9.步骤4，将待测回波数据进行时域、频域特征提取，输入到子分类器并通过加权投票得到最终分类结果。
10.优选地，时域特征包括：时域波形标准差、时域波形不规则度、时域波形平滑度，时域波形幅值变异系数，时域波形相位变异系数，时域波形等效散射中心维度，时域波形强点功率占比；频域特征包括：地物杂波特性、主峰谱宽、展宽谱宽、功率谱熵、功率谱一阶原点
矩、二阶中心距、相对有效反射面积、峰值能量比、高频成分含量、奇异值波形衰减速度、奇异值波形稳定程度、奇异值波形能量集中程度、emd分解熵、emd分解峰值比，共计21个特征值。
11.本发明与现有技术相比，其显著优点为：本发明实现简单，在低分辨地面雷达系统中的目标分类识别准确性较高，能够高效而又稳定的进行目标分类识别；本发明利用雷达回波数据，提取融合的时频特征，建立融合的多分类器，通过五折交叉验证进行分类器训练，采用加权投票法得到最终结果，提高了分类识别系统的性能。
12.下面结合附图对本发明作进一步详细描述。
附图说明
13.图1是本发明实施例的流程图。
14.图2是基于特征提取的加权投票多分类器地面雷达目标识别系统组成框图。
具体实施方式
15.如图1所示，一种基于特征提取的加权投票多分类器地面雷达目标自动识别方法包括：利用训练集数据对时域回波信号进行时域特征提取，对下变频处理后的基带回波时域信号进行fft变换，进行频域特征提取；通过提取的21个时域和频域特征进行5种子分类器的5折交叉验证训练；根据混淆矩阵计算出5个子分类器的精度并进行归一化作为各个分类器的权值；利用训练好的分类器对测试数据集进行目标分类识别，通过加权投票得到的分类结果进行判定目标是人还是车。具体步骤为：
16.步骤1，对时域回波信号进行时域特征提取，对下变频处理后的基带回波时域信号进行fft变换，进行频域特征提取；
17.具体地，时域特征包括：时域波形标准差、时域波形不规则度、时域波形平滑度，时域波形幅值变异系数，时域波形相位变异系数，时域波形等效散射中心维度，时域波形强点功率占比；
18.频域特征包括：地物杂波特性、主峰谱宽、展宽谱宽、功率谱熵、功率谱一阶原点矩、二阶中心距、相对有效反射面积、峰值能量比、高频成分含量、奇异值波形衰减速度、奇异值波形稳定程度、奇异值波形能量集中程度、emd分解熵、emd分解峰值比。
19.进一步的，设雷达时域回波序列为x(i),i＝1,2,...,n，则时域波形标准差的计算方法为：
[0020][0021]
式中，为目标样本时域波形内n个采样点的均值。
[0022]
进一步的，时域波形幅值变异系数计算方法为：
[0023]
首先对时域序列进行归一化，有：
[0024][0025]
则归一化方差为：
[0026][0027]
其中，x
mean
为归一化值的均值。
[0028]
则幅值变异系数为：
[0029]camplitude
＝x
std
/x
mean
ꢀꢀ
(4)
[0030]
进一步的，时域波形相位变异系数计算方法为：
[0031]
首先，根据iq数据求出序列各点的相位值：
[0032][0033]
则序列相邻两点间的相位差为：
[0034][0035]
有
[0036][0037]
则时域波形相位变异系数为：
[0038][0039]
进一步的，时域波形等效散射中心维度的计算方法为：
[0040][0041]
其中，为单位阶跃函数。
[0042]
进一步的，时域波形强点功率占比的计算方法为：
[0043][0044]
其中，y(i),i＝1,2,...,n为时域回波由大到小降序排列形成的序列。
[0045]
进一步的，对时域波形求fft，得到功率谱x(n),n＝1,2,...,n，并求其对数值得到功率值p(n),n＝1,2,...,n，则地物杂波特性的计算方法为：
[0046]
设杂波区的谱线为nc根，目标所在谱线为n
t
，目标的功率为pr，去除目标所在谱线及其左右各一根谱线后对杂波区外的谱线功率求和平均得到检测门限功率，对杂波区内的谱线求和平均得到杂波功率，两者的比值即为地物杂波特性。
[0047]
进一步的，相对有效波形反射截面积的计算方法为：
[0048]
σ(db)＝10lg(pr) 40lg(r
t
)-40lg(r
max
)-10lg(gr)
ꢀꢀ
(11)
[0049]
进一步的，emd分解熵、emd分解能量比的计算方法为：
[0050]
(1)先对时域信号进行emd分解。
[0051]
emd分解的步骤为：
[0052]
(1.1)初始化r(t)＝x(t)，对r(t)的上包络和下包络分别进行三次样条插值法
[0053]
拟合出采样波形的包络曲线max_r(t)和min_r(t)；
[0054]
(1.2)对两条包络曲线求均值得到m(t)，计算出中间信号
[0055]
d(t)＝r(t)-m(t)
ꢀꢀ
(12)
[0056]
(1.3)判断中间信号d(t)是否满足imf的两个条件，如果满足，该信号就是一个imf分量；如果不满足，令r(t)＝d(t)，迭代(1.1)～(1.3)。
[0057]
(1.4)r(t)＝x(t)-d(t)，重复(1.1)～(1.4)，得到k个imf分量。
[0058]
(2)求得本质模态函数(imf1)分量后进行傅里叶变换；
[0059]
(3)对(2)的计算结果求波形熵作为emd的分解熵；
[0060]
(4)将imf1频带的能量总和与原始频谱的对应谱线范围的能量总和之比作为emd分解能量比。
[0061]
步骤2，利用提取的特征进行5折交叉验证朴素贝叶斯分类器、决策树分类器、lda分类器、k邻近分类器和支持向量机分类器训练，并通过贝叶斯优化对支持向量机和k邻近分类器的超参数进行优化，得到最终的五个子分类器；
[0062]
步骤3，根据混淆矩阵计算五个子分类器的总体精度，通过归一化计算将结果作为各个分类器的权值w(i)，即为
[0063][0064]
步骤4，将待测回波数据进行时域、频域特征提取得到21个特征值后，分别输入到五个子分类器得到初步的分类结果，并通过加权投票得到最终分类结果
[0065][0066]
其中，
[0067]
如果计算值r为正，则判为车；如果计算值为0，表示未识别；如果计算值r为负，则判为人。
[0068]
下面结合附图2及具体实施例对本发明做进一步详细说明。
[0069]
实施例
[0070]
如图2所示，本发明为基于特征提取的加权投票多分类器地面雷达目标识别系统，具体步骤为：
[0071]
步骤1，对时域回波信号进行时域特征提取，对下变频处理后的基带回波时域信号进行fft变换，进行频域特征提取；
[0072]
时域特征包括：(1)时域波形标准差、(2)时域波形不规则度、(3)时域波形平滑度，(4)时域波形幅值变异系数，(5)时域波形相位变异系数，(6)时域波形等效散射中心维度，(7)时域波形强点功率占比；频域特征包括：(8)地物杂波特性、(9)主峰谱宽、(10)展宽谱宽、(11)功率谱熵、(12)一阶原点矩、(13)二阶中心距、(14)相对有效反射面积、(15)峰值能量比、(16)高频成分含量、(17)奇异值波形衰减速度、(18)奇异值波形稳定程度、(19)奇异值波形能量集中程度、(20)emd分解熵、(21)emd分解峰值比。具体计算方法为：
[0073]
步骤1-1，设雷达时域回波序列为x(i),i＝1,2,...,n，则时域波形标准差的计算方法为：
[0074][0075]
式中，为目标样本时域波形内n个采样点的均值。
[0076]
步骤1-2，时域波形不规则度的计算方法为：
[0077][0078]
步骤1-3，时域波形平滑度根据一阶差分编码值求得，计算方法为：
[0079]
首先计算时域波形序列的一阶差分值：
[0080][0081]
接着求当前点与上一点差分值符号的异同，可通过乘积的正负判断：
[0082]
sign(i)＝firstorder_diff(i)*firstorder_diff(i-1)
ꢀꢀ
(4)
[0083]
接着，计算一阶差分编码值，设初始一阶差分编码值为1，则
[0084][0085]
最后，对一阶差分编码值求和得到时域波形平滑度，有：
[0086][0087]
步骤1-4，时域波形幅值变异系数计算方法为：
[0088]
首先对时域序列进行归一化，有：
[0089][0090]
则归一化方差为：
[0091][0092]
其中，x
mean
为归一化值的均值。
[0093]
则幅值变异系数为：
[0094]camplitude
＝x
std
/x
mean
ꢀꢀ
(9)
[0095]
步骤1-5，时域波形相位变异系数计算方法为：
[0096]
首先，根据iq数据求出序列各点的相位值：
[0097][0098]
则序列相邻两点间的相位差为：
[0099][0100]
有
[0101][0102]
则时域波形相位变异系数为：
[0103][0104]
步骤1-6，时域波形等效散射中心维度的计算方法为：
[0105][0106]
步骤1-7，时域波形强点功率占比的计算方法为：
[0107][0108]
其中，y(i),i＝1,2,...,n为时域回波由大到小降序排列形成的序列。
[0109]
步骤1-8，对时域波形求fft，得到功率谱x(n),n＝1,2,...,n，并求其对数值得到功率值p(n),n＝1,2,...,n，则地物杂波特性的计算方法为：
[0110]
设杂波区的谱线为nc根，目标所在谱线为n
t
，目标的功率为pr，去除目标所在谱线及其左右各一根谱线后对杂波区外的谱线功率求和平均得到检测门限功率，对杂波区内的谱线求和平均得到杂波功率，两者的比值即为地物杂波特性。
[0111]
步骤1-9，频域主峰谱宽的计算方法为：
[0112]
去除杂波区的谱线后，对主频域间的谱线找最大值，得到目标所在谱线n
t
及其功率值pr，以目标谱线为中心记录左、右功率下降6db的谱线位置，左侧记为n
low
，右侧记为n
high
，则频域主峰谱宽的计算方法为：
[0113]
width_6db＝n
high-n
low
1
[0114]
步骤1-10，展宽谱宽的计算方法为：
[0115]
取阈值为所有谱线功率值的均值，以目标谱线为中心，遍历目标谱线左右各十根，记录谱线值与阈值的差值大于目标谱线与阈值差值的0.4倍的谱线数量，将该数量记作展宽谱宽。
[0116]
步骤1-11，功率谱熵的计算方法为：
[0117]
对杂波区外的谱线进行规范化，再求其信息熵值，有
[0118][0119][0120]
步骤1-12，功率谱一阶原点矩的计算方法为：
[0121]
将谱线进行归一化
[0122]
[0123]
可得到一阶原点矩为
[0124][0125]
步骤1-13，二阶中心矩为
[0126][0127]
步骤1-14，相对有效波形反射截面积的计算方法为：
[0128]
根据雷达接收功率公式可推得
[0129][0130]
其中，k雷达固定参数，可当做常量，则可推得相对有效反射截面积(db)的计算方法为：
[0131]
σ(db)＝10lg(pr) 40lg(r
t
)-40lg(r
max
)-10lg(gr)
ꢀꢀ
(22)
[0132]
其中，r
t
为目标的距离单元号，pr为目标接收功率，gr接收通道增益。
[0133]
步骤1-15，峰值能量比的计算方法为：
[0134]
将主峰谱宽内的谱线值进行求和，其与整个频带内的功率和之比即为峰值能量比。
[0135]
步骤1-16，高频成分含量的计算方法为：
[0136]
对功率谱进行归一化得到
[0137][0138]
将归一化的功率谱按照间隔p划分为个频带，则高频成分序列为
[0139][0140]
对其再进行归一化，选取阈值，判断序列中超过阈值的个数，记为高频成分含量。
[0141]
步骤1-17～19，奇异值波形衰减速度、奇异值波形稳定程度、奇异值波形能量集中程度的计算方法为：
[0142]
对功率谱求自相关矩阵，再进行奇异值分解，将奇异值中超过阈值的作为主成分含量，将主成分个数作为波形衰减速度、对奇异值求方差作为波形稳定程度、对奇异值求波形熵为能量集中程度。
[0143]
步骤1-20～21，emd分解熵、emd分解能量比的计算方法为：
[0144]
(1)先对时域信号进行emd分解。
[0145]
emd分解的步骤为：
[0146]
(1.1)初始化r(t)＝x(t)，对r(t)波形的上包络和下包络分别进行三次样条插值法拟合出包络曲线max_r(t)和min_r(t)；
[0147]
(1.2)对两条包络曲线求均值得到m(t)，计算出中间信号
[0148]
d(t)＝r(t)-m(t)
ꢀꢀ
(25)
[0149]
(1.3)判断中间信号d(t)是否满足imf的两个条件，如果满足，该信号就是一个imf分量；如果不满足，令r(t)＝d(t)，迭代(1.1)～(1.3)。
[0150]
(1.4)r(t)＝x(t)-d(t)，重复(1.1)～(1.4)，得到k个imf分量。
[0151]
(2)求得本质模态函数(imf1)分量后进行傅里叶变换；
[0152]
(3)对(2)的计算结果求波形熵作为emd的分解熵；
[0153]
(4)将imf1频带的能量总和与原始频谱的对应谱线范围的能量总和之比作为emd分解能量比。
[0154]
步骤2，用提取后的特征进行5折交叉验证朴素贝叶斯分类器、决策树分类器、lda分类器、k邻近分类器和支持向量机分类器训练，并通过贝叶斯优化对支持向量机和k邻近分类器的超参数进行优化，得到最终的五个子分类器。
[0155]
将上述的7个时域特征和14个频域特征值作为特征量，选取m个样本，记为(fi,li),i∈[1,m]，将m个样本分为5组，每次选择其中4组作为训练样本集,剩余1组作为测试样本集，以此类推遍历5组数据得出5次实验结果，并将5次结果均值作为识别结果，进行子分类器5折交叉验证训练。
[0156]
步骤3，根据混淆矩阵计算五个子分类器的精度权值，通过加权投票的方式得到最终的分类器；
[0157]
将训练得到的五个子分类器根据混淆矩阵计算分类器权值w(i)，即为
[0158][0159]
其中，ai为第i个分类器的分类精度。
[0160]
则加权投票多分类器的输出结果即为
[0161][0162]
其中，
[0163]
步骤4，将待测回波数据进行时域、频域特征提取，输入到子分类器并通过加权投票得到最终分类结果。
[0164]
将待测回波数据进行时域、频域特征提取得到21个特征值，将特征值分别输入到五个子分类器得到初步的分类结果并通过步骤3的权值进行加权投票计算得到最终分类结果，如果计算值为正，则判为车；如果计算值为0，表示未识别；如果计算值为负，则判为人。
[0165]
本发明实现简单，先根据雷达时域回波提取特征，进行fft后在频域提取波形特征，共计21种特征，特征量丰富可靠；采用了五种常规分类器，根据人、车各5000个样本进行五折交叉验证训练分类器，再通过人、车各1000个测试数据进行测试，通过投票法集成多分类器的结果，准确率达到95.4％，展现了良好的分类识别能力。