本发明属于目标识别技术领域，具体的说是涉及一种多旋翼无人机识别中的块概率分布特征提取方法。本发明首先通过对无人机目标的雷达回波数据集进行短时付氏变换，获得时频谱图，然后，将时频谱图划分为多个具有交叠区域的子块，每一子块的元素按列顺序构成一个矢量，对该矢量进行量化，并统计整幅时频谱图的量化索引的直方图，作为块概率分布特征，完成对多旋翼无人机的识别。由于充分利用了子块的统计分布信息，能够提取到有关目标雷达回波数据中的局部结构分布特征，增大了目标之间的区分细节，从而改善了对目标的识别率。对4类多旋翼无人机的仿真实验结果验证了方法的有效性。
1.一种多旋翼无人机识别中的块概率分布特征提取方法，其特征在于，包括以下步骤：S1、设多旋翼无人机的雷达回波训练数据序列为n维列矢量xij，i＝1,2,…g，j＝1,2…Ni，其中，i表示无人机的类别，g表示类别数，Ni表示第i类无人机目标的训练序列数，总的训练序列数为S2、对第i类无人机目标的训练样本数据序列xij进行短时付氏变换，得时频谱图Sij：Sij＝[sij,km]K×M
其中，sij,km表示时频谱图Sij中的元素，k＝1,2,…K，m＝1,2,…M，k为行下标，表示频率变化方向，m为列下标，表示时间变化方向，K是Sij的行数，M是Sij的列数；
S3、以4x4的窗口在时频谱图Sij中按元素进行滑动，得到多个子块，组成子块集Pij：Pij＝[Qij,1 Qij,2 … Qij,L]Qij,l＝[qij,l,ef]4×4,e,f＝1,2,3,4l＝1,2,…L
其中，Qij,l表示子块集Pij中的第l个子块矩阵，qij,l,ef表示子块矩阵Qij,l中元素，L表示子块总数；将子块矩阵Qij,l中的元素按列构成如下子块列矢量qij,l：T
qij,l＝[qij,l,11 qij,l,21 … qij,l,44]S4、对时频谱图的子块列矢量qij,l进行量化，获得时频谱图的子块量化索引矢量，计算子块列矢量qij,l与码书中的码字的欧氏距离dij,l,u：dij,l,u＝||qij,l‑cu||2,u＝1,2,…64其中，||·||2表示矢量的二阶范数，cu是码书中的码字，该码书包含64个码字，由训练数据集进行聚类建立；取最小距离对应的码字下标作为子块Qij,l的量化索引将时频谱图中的所有子块对应的量化索引构成子块量化索引矢量uij：S5、提取块概率分布特征：
统计子块量化索引矢量uij中的各元素重复出现次数，并组成矢量：[λij,1 λij,2 … λij,64]
其中，λij,l为码字下标l的重复出现次数，l＝1,2,…,64，则时频谱图Sij对应的块概率分布特征hij：
一种多旋翼无人机识别中的块概率分布特征提取方法
技术领域
[0001] 本发明属于目标识别技术领域，具体的说是涉及一种多旋翼无人机识别中的块概率分布特征提取方法。
背景技术
[0002] 随着无人机在各领域的广泛应用，给国家的安全带来了巨大的威胁，有必要对无人机进行有效的监控、管理及打击。因此，准确识别出无人机的类型对反无人机作战中具有非常重要的实用意义。
[0003] 目前，子空间方法是一种识别无人机的有效方法，主要通过无人机目标的训练雷达回波数据的时频谱图建立投影子空间，提取目标特征。但是，常规子空间提取的特征属于全局结构特征，而忽视了局部结构特征，降低了目标之间的细节区分度，因此，基于常规子空间的多旋翼无人机识别方法的目标识别率有进一步改善的余地。
发明内容
[0004] 针对上述问题，本发明提出了一种块概率分布特征提取方法，通过将时频谱图划分为子块，并统计子块的分布信息，能够提取到有关目标细节的局部结构特征，增大了目标之间的区分度，从而改善了对目标的识别率。
[0005] 为实现上述目的，本发明的技术方案是：
[0006] 一种多旋翼无人机识别中的块概率分布特征提取方法，包括以下步骤：
[0007] S1、设多旋翼无人机的雷达回波训练数据序列为n维列矢量xij，i＝1,2,…g，j＝1,2…Ni，其中，i表示无人机的类别，g表示类别数，Ni表示第i类无人机目标的训练序列数，总的训练序列数为
[0008] S2、对第i类无人机目标的训练样本数据序列xij进行短时付氏变换，得时频谱图Sij：
[0009] Sij＝[sij,km]K×M
[0010] 其中，sij,km表示时频谱图Sij中的元素，k＝1,2,…K，m＝1,2,…M，k为行下标，表示频率变化方向，m为列下标，表示时间变化方向，K是Sij的行数，M是Sij的列数；
[0011] S3、以4x4的窗口在时频谱图Sij中按元素进行滑动，得到多个子块，组成子块集Pij：
[0012] Pij＝[Qij,1 Qij,2 … Qij,L]
[0013] Qij,l＝[qij,l,ef]4×4,e,f＝1,2,3,4
[0014] l＝1,2,…L
[0015] 其中，Qij,l表示子块集Pij中的第l个子块矩阵，qij,l,ef表示子块矩阵Qij,l中元素，L表示子块总数；将子块矩阵Qij,l中的元素按列构成如下子块列矢量qij,l：
[0016] qij,l＝[qij,l,11 qij,l,21 … qij,l,44]T
[0017] S4、对时频谱图的子块列矢量qij,l进行量化，获得时频谱图的子块量化索引矢量，计算子块列矢量qij,l与码书中的码字的欧氏距离dij,l,u：
[0018] dij,l,u＝||qij,l‑cu||2,u＝1,2,…64
[0019] 其中，||·||2表示矢量的二阶范数，cu是码书中的码字，该码书包含64个码字，由训练数据集进行聚类建立；取最小距离对应的码字下标作为子块Qij,l的量化索引[0020]
[0021] 将时频谱图中的所有子块对应的量化索引构成子块量化索引矢量uij：
[0022]
[0023] S5、提取块概率分布特征：
[0024] 统计子块量化索引矢量uij中的各元素重复出现次数，并组成矢量：
[0025] [λij,1 λij,2 … λij,64]
[0026] 其中，λij,l为码字下标l的重复出现次数，l＝1,2,…,64，则时频谱图Sij对应的块概率分布特征hij：
[0027]
[0028] 本发明的有益效果是：本发明充分利用了子块的统计分布信息，能够提取到有关目标雷达回波数据中的局部结构分布特征，增大了目标之间的区分细节，从而改善了对目标的识别率。
具体实施方式
[0029] 下面结合仿真实验说明本发明的实用性。
[0030] 仿真实验设计了4种类型的无人机，包括三旋翼无人机、四旋翼无人机、六旋翼无人机、八旋翼无人机，其仿真参数如表1所示。仿真雷达参数包括：雷达载频为24GHz；脉冲重复频率为100KHz；目标与雷达间的距离200m；无人机相对于雷达的俯仰角为10°、方位角为30°
[0031] 每类目标记录10s的雷达回波信号，并将其分为固定长度0.05s的段(至少包含一个旋转周期)，段之间的重叠为50％，每段包含0.05×100000＝5000个雷达回波采样数据点，每类共400段。在400段中随机选取200段作为训练数据集，其余200段作为测试数据集，则4类目标的训练数据集总共包括800个段，测试数据集包括800个段。对选取的训练数据集，利用本文方法提取块概率分布特征，建立模板库，同样提取测试样本的块概率分布特征，采用最小距离分类器进行分类，对4类多旋翼无人机的平均正确识别率达到97％。其中，信噪比为20dB。
[0032] 表1四种无人机的仿真参数
[0033]
[0034]
