无人机与飞鸟目标雷达识别方法与流程

1.本发明涉及雷达信号处理、数据处理、目标识别等问题，可广泛用于探鸟雷达、反无人机雷达等领域。


背景技术：

2.随着无人机领域技术的快速发展，因其具有成本低、体积小、操作简单等优势，已经广泛应用于军事领域和民用领域，但随着非合作无人机黑飞、干扰等事件的频发，已经成为一大重要隐患，而在航空器起飞或者飞行过程与鸟类相撞的事件也层出不穷，飞鸟的防范压力也越来越大，这两类目标都对机场、安保等重要场所均产生了严重威胁，而一般根据不同的作战环境要求，针对这两者类型目标，在发现后需要采取不同的反制手段和措施，如机场一般在发现飞鸟后，会采取驱鸟设备将其驱赶，而对无人机目标则会采取预警、协调等措施，因此对无人机和飞鸟的识别是极其必要的。
3.目前，对无人机和飞鸟的监视与识别方法主要分为两类，一类从目标回波特性角度出发，由于无人机多旋翼转动与飞鸟翅膀扇动产生的回波信号存在特性差异，所以一般是采用时频分析算法提取两者的微多普勒特征。法国航天航空研究中心采用dvb-t信号进行了多旋翼无人机微多普勒测量实验，荷兰robin公司雷达则采用发射调频连续波波形，同时具备无人机探测和识别的功能，在随后多年发展中，国内外单位通过设计优化雷达工作系统、工作模式、发射接收波形、延长驻留时间等方法来进一步提取精细化特征，但在实际应用到复杂多变的环境中，目前还是较为困难，尤其当目标处于远距离时，两者的回波特征微弱导致特征提取困难、时频分辨力低等原因，限制了大范围的应用与发展。另一类方法则是从无人机与飞鸟的飞行轨迹特征方面出发，通过记录大量历史数据，熟悉了解飞鸟的活动规律和飞行轨迹，从数据层提取无人机与飞鸟的运动轨迹特征，目前有采用多运动模型并通过计算模型间转换频率的方法来进行识别，也有利用航迹轨迹信息建立特征量并利用机器学习或深度学习的方法，这几种方法都进一步验证了该类方法的可行性。


技术实现要素：

4.本发明提出了一种无人机与飞鸟目标雷达识别方法，首先基于大量历史雷达系统航迹数据分析和研究无人机与飞鸟在运动轨迹上的特征差异，提出了以航迹层信息为基础，建立航向角标准差、航向震荡频率、速度均值、速度标准差，速度震荡频率五个特征量，并采用支持向量机算法，将历史样本分为训练样本和测试样本，对训练样本采用交叉训练方法获得最优超参数，并利用测试样本验证了识别分类的准确率。
5.本发明提出的无人机与飞鸟目标雷达识别方法，与现有识别方法的不同点在于：第一，未采用复杂的时频分析方法，而是采用了轻量性的航迹数据集，大大降低了雷达系统的识别成本，且普适性强；第二，为了更好体现两者飞行轨迹上的特征差异，采用了支持向量机算法模型，提出了时间相关的航向震荡频率、速度震荡频率特征向量描述方法，并结合其他特征量，实现在线动态修正并预测，能够更好提高识别准确率。
附图说明
6.图1是本发明一实施例的无人机与飞鸟目标识别框图；
7.图2是本发明一实施例的典型无人机飞行轨迹示意图；
8.图3是本发明一实施例的典型鸟类飞行轨迹示意图；
9.图4是本发明一实施例的无人机目标跟踪识别ppi图；
10.图5是本发明一实施例的飞鸟目标跟踪识别ppi图。
具体实施方式
11.本发明提出了一种无人机与飞鸟目标雷达识别方法，首先在离线条件下，通过雷达系统记录的大量历史实验航迹数据进行分析，再对比无人机和飞鸟在运动轨迹上的差异性，进行特征分析和特征计算，确定了航向角标准差、航向震荡频率、速度均值、速度标准差、速度震荡频率五个特征量，然后将数据集分为训练集、测试集，将训练集数据作为输入，利用svm算法计算出模型和最优超参数，并利用测试集验证了分类的准确率，具体流程如图1。
12.1.1.特征分析
13.从无人机与鸟类的活动轨迹来描述，一般飞鸟的飞行轨迹相对灵活且时间较短，而典型无人机作业一般是按照预先设定的航线任务来飞行的，常见地分为连续直线飞行、航点悬停、折返飞行、多方位多点飞行等路线，因此该类无人机目标的飞行轨迹时间较长，且一般相对稳定，但如果是人工干预控制操作飞行，其在转向、折返、悬停上则具有较高自由度，将连续人工干预操作无人机机动飞行称为强机动模式，按照航点任务作业称为弱机动模式。
14.利用某相控阵雷达系统平台，其具备探测跟踪无人机和飞鸟功能，通过大量实验工作与历史航迹数据记录分析，典型的无人机和鸟类在航迹层级别的飞行轨迹如图2、图3，从距离、方位角、航迹长度可以直观体现出两者的区别，为了完善数据库，其中无人机航迹数据层中包含了大量无人机不同作业类型的航迹，例如手动操作飞行、直线飞行、8字飞行、折返飞行、悬停飞行等不同轨迹。
15.下面从运动轨迹的局部和总体特征来分别描述无人机与飞鸟的特征量的差异，这些特征量分别是航向角标准差、航向震荡频率、平均速度、速度标准差、速度震荡频率。
16.航向角标准差：对于一般飞鸟的轨迹，其方向性是相对比较随机的，标准差变化范围较大，而无人机在执行弱机动模式飞行时，其标准差值较小，对于这类可以通过对比航向角标准差的差异能较明显体现出两者飞行轨迹的区别，而在强机动模式下时则值域变换范围较大，较难区分。
17.航向震荡频率：上述航向角标准差只能体现两者的总体特征，并不能够很好描述局部航向变化特征，为此引入航向震荡频率来表征两者的区别，无人机在弱机动模式下其航向震荡次频率值较小，而在飞行过程中如果只是部分时间人工干预控制操作让其处于强机动模式飞行，随着跟踪次数增加，航向震荡频率值也会逐渐变小，而如果是一直处于强机动模式飞行，则航向震荡频率值则一直处于较高的值，因此对于无人机其航向震荡频率值总是处于较大或者较小值，而飞鸟的飞行轨迹较为灵活且相对平滑，加上一般雷达对其跟踪时间较短，因此其航向震荡频率数值处于较小与较大之间范围，所以航向震荡频率具备
区分无人机与飞鸟的潜力。
18.平均速度：对于区分飞鸟和一般的无人机不太明显，但是仍然可以作为区分一些大型无人机和飞鸟目标的特征，例如一般固定翼无人机的飞行速度远超于飞鸟。
19.速度标准差：对于飞鸟而言，其速度相对比较稳定，其标准差的数值一般较小，而无人机由于其飞行相对自由，例如悬停、加速、减速等等机动都会影响飞行速度，因此在其机动过程中标准差的数值相对较大，能够作为区分无人机与飞鸟的一个特征。
20.速度震荡频率：对于飞鸟而言，其速度不会产生突变的情况，而在无人机转向飞行或者悬停过程中，需要加速或者减速，而上述速度标准差只能体现两者的总体特征，并不能够很好描述局部速度变化特征，为此引入速度震荡频率来体现速度局部变化特征，一般飞鸟的速度震荡频率数值很小，对于无人机而言，在加速或者减速过程虽然只需要较短时间，但此时的速度震荡频率值相对较大，因此可以通过长时间积累，也可以作为区分飞鸟与强机动模式下无人机的飞行特征。
21.1.2.特征计算
22.定义某条长度为n的航迹目标集合符号记为s＝{s(1),s(2),...,s(n)}，每个数据采样点记为符号s，其中第i个航迹采样点分别代表目标在该点的距离、方位角、俯仰角和速度。
23.1.航向角标准差
24.为方便计算航向角，将航迹数据集合下每个采样点s(i)的距离、方位、俯仰转化为直角坐标系下(x(i),y(i),z(i))，转换公式如下：
[0025][0026]
根据上一时刻s(i-1)和当前时刻s(i)算出直角坐标系下的位置信息x(i-1)、y(i-1)、x(i)、y(i)，计算当前i时刻的航向角，符号记为h(i)，计算公式如下：
[0027][0028]
其中计算航向角时需考虑x(i)-x(i-1)、y(i)-y(i-1)值所处象限区域，从而将-90
°
～90
°
的值域转化为0
°
～360
°
，那么航迹集合s航向角集合h为
[0029]
h＝{h(2),h(3),...,h(n)}
ꢀꢀꢀꢀꢀꢀ
(3)
[0030]
要计算航向角的标准差，首先要计算航向角的均值，但是存在某些目标轨迹的航向角跨0
°
或360
°
或者由于测量角度误差的影响而导致产生的跨越情况，因此对于航向角不能采用常规均值方法计算，本发明采用向量化的方法来描述这类周期性角度的数据，计算周期性数据均值与标准差的方法是首先将第i个航向角h(i)转化为二维平面上的单位向量该航向角在x轴坐标值对应于cosh(i)，y轴坐标值对应于sinh(i)，公式如下：
[0031][0032]
把h集合中n-1个航向角转换成矢量并且相加得到合成向量计算公式如下，航
向角均值hm和标准差h
std
分别定义为：
[0033][0034][0035][0036]
其中航向角均值需要考虑所处象限区域，将航向角均值转化为0
°
～360
°
。
[0037]
2.航向震荡频率
[0038]
为了反映航向变化局部特征，引入航向震荡频率，首先对航迹目标集合s之间的相邻航向角计算差值，记为符号dh(i)，公式如下：
[0039]dh
(i)＝h(i)-h(i-1),i≥3
ꢀꢀꢀꢀꢀ
(8)
[0040]
并判断角度差dh(i)与设定阈值γh之间的大小，计算出航向角之间的差异符号值o(i)，定义如下：
[0041][0042]
再根据航向角符号值集合o，定义两种航向震荡模式：
[0043]
o(i-1)+o(i)＝0且o(i-1)≠o(i)
ꢀꢀꢀꢀꢀ
(10)
[0044]
o(i-1)+o(i+1)＝0且o(i-1)≠o(i+1),o(i)＝0
ꢀꢀ
(11)
[0045]
当满足上述公式(10)或者(11)时，判决处于震荡模式并统计震荡次数，符号记为uh，
[0046][0047]
其中γh为航向角差异判定的阈值，可以根据实际雷达测角误差设定，在得出航向震荡次数后，定义航向震荡频率为
[0048][0049]
3.平均速度
[0050]
航迹目标集合s的平均速度计算公式如下：
[0051][0052]
4.速度标准差
[0053]
航迹目标集合s的速度标准差计算公式如下：
[0054][0055]
5.速度震荡频率
[0056]
为了反映速度变化局部特征，引入速度震荡频率，在当前时刻i取历史l长度段样
本点的速度值，计算该长度下的速度均值v
ml
，定义如下：
[0057][0058]
计算当前点速度v(i)和历史均值速度v
ml
之间的差值，符号记为τv(i)，定义如下
[0059]
τv(i)＝v(i)-v
ml
,i≥l+1
ꢀꢀꢀꢀꢀ
(17)
[0060]
设定速度变化阈值γv，当差值大于γv，则该点定义为速度震荡，符号记为uv(i)，定义如下：
[0061][0062]
其中可根据实际雷达跟踪数据率和跟踪精度调整阈值γv和l长度，速度震荡频率定义为
[0063][0064]
最后得出航迹计算特征量为航向角标准差、航向震荡频率、速度均值、速度标准差、速度震荡频率，符号为g＝[h
std
,fh,vm,v
std
,fv]，利用这五个特征量可能在跟踪开始阶段误判目标的类型，但随着目标稳定跟踪时间和次数的增加，时间相关的特征量值可以通过在线动态修正并重新对目标类型进行预测，识别的准确率也会逐步提高。
[0065]
1.3.离线训练与测试
[0066]
1.模型选择
[0067]
在上述得到5个特征量后，考虑样本特征量是非线性的且为了防止过拟合，采用软间隔的支持向量机算法，其模型公式如下：
[0068][0069]
其中w为超平面法向量，c为常数，且大于0，ξ为松弛变量，φ(x)为映射函数，b是超平面的截距，xi为样本矢量，上式中计算会引入核函数κ(xi,xj)，本发明中采用高斯核函数，公式如下：
[0070][0071]
式中σ为高斯核函数的带宽，后续采用smo(sequential minimal optimization)算法求解模型。
[0072]
2.训练模型与验证
[0073]
在确定支持向量机模型后，不同大小的超参数对识别的准确性有较为明显影响，在高斯核函数中需要调整的超参数是惩罚因子c和核函数的带宽σ，首先利用随机采样法将航迹数据集分为训练和测试样本集，然后将训练样本数据集采用5折交叉训练法，寻得最优超参数。获得训练模型和最优超参数后，再利用测试样本集对模型进行验证评估，其中测试
样本集数中无人机和鸟类航迹数量为114、76条，识别混淆矩阵如表1，其中准确识别率达到了87％，验证了该方法的有效性。
[0074]
表1
[0075][0076]
1.4.在线实验结果
[0077]
实验过程中以某相控阵雷达系统为平台，部署在飞鸟较为活跃的环境，无人机型号为大疆精灵4，考虑到实验的完整性，无人机模拟了多次不同作业和人工操作飞行场景的任务，其中在某次跟飞实验中，雷达目标跟踪与识别的平面显示器(ppi)画面如图4、图5，由图可见，无人机目标的批号为223，方位角在356
°
方向，运动轨迹为首先经过一段背离雷达方向的飞行，在达到指定航向任务点时拐弯调转方向，并按照进入雷达的方向飞行，而飞鸟目标的批号为337，方位角在33
°
方向，飞行轨迹较为无序，雷达在稳定持续跟踪情况下，均正确识别了无人机与飞鸟目标，实验结果验证了该方法在工程角度实现的轻量性与可行性。

技术特征：
1.一种无人机与飞鸟目标雷达识别方法，其特征在于，包括：步骤s1，对离线状态下采集到的无人机和飞鸟数据进行分析，根据无人机任务作业的特点进行强机动模式与弱机动模式描述，基于这两种模式下无人机与飞鸟特征上的差异性，提出了时间相关的航向震荡频率与速度震荡频率特征量描述方法，并结合速度均值、速度标准差、航向角标准差三个特征量，从而提高至多维特征来增加区分无人机与飞鸟的能力。2.如权利要求1所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s1，包括：步骤s11，离线状态下，基于航向角标准差、航向震荡频率、平均速度、速度标准差、速度震荡频率共五个特征量，采用支持向量机算法，并将历史数据分为训练集和测试集，对训练集采用交叉训练方法获得最优超参数和模型，并利用测试集进行验证，在线状态下，利用最优模型和超参数，实时计算稳定跟踪目标的五个相关特征量，并且随着跟踪次数的增加可以重复修正特征量值，从而逐步提升识别的准确率。3.如权利要求2所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s11，包括：计算航向角标准差采用向量化的方式，将第i个航向角h(i)转化为二维平面上的单位向量该航向角在x轴坐标值对应于cosh(i)，y轴坐标值对应于sinh(i)，矢量把航向角h集合中n-1个航向角转换成矢量并且相加得到合成向量则航向角标准差h
std
定义为避免了由于测量误差带来航向角跨0
°
或360
°
的影响。4.如权利要求3所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s11，包括：对航迹目标集合s之间的相邻航向角计算差值，记为符号d
h
(i)，公式如下：d
h
(i)＝h(i)-h(i-1)，i≥3
ꢀꢀꢀꢀꢀꢀꢀꢀ
(8)并判断角度差d
h
(i)与设定阈值γ
h
之间的大小，计算出航向角之间的差异符号值o(i)，定义如下：再根据航向角符号值集合o，定义两种航向震荡模式：o(i-1)+o(i)＝0且o(i-1)≠o(i)
ꢀꢀꢀꢀꢀꢀꢀ
(10)o(i-1)+o(i+1)＝0且o(i-1)≠o(i+1),o(i)＝0
ꢀꢀꢀꢀꢀꢀ
(11)当满足上述公式(10)或者(11)时，判决处于震荡模式并统计震荡次数，符号记为u
h
，
其中，γ
h
为航向角差异判定的阈值，可以根据实际雷达测角误差设定，在得出航向震荡次数后，定义航向震荡频率为5.如权利要求4所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s11，包括：航迹目标集合s的平均速度计算公式如下：6.如权利要求5所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s11，包括：航迹目标集合s的速度标准差计算公式如下：7.如权利要求6所述的无人机与飞鸟目标雷达识别方法，其特征在于，所述步骤s11，包括：在当前时刻i取历史l长度段样本点的速度值，计算该长度下的速度均值v
ml
，定义如下：计算当前点速度v(i)和历史均值速度v
ml
之间的差值，符号记为τ
v
(i)，定义如下τ
v
(i)＝v(i)-v
ml
，i≥l+1
ꢀꢀꢀꢀꢀꢀ
(17)设定速度变化阈值γ
v
，当差值大于γ
v
，则该点定义为速度震荡，符号记为u
v
(i)，定义如下：其中，根据实际雷达跟踪数据率和跟踪精度调整阈值γ
v
和l长度，速度震荡频率定义为

技术总结
本发明提出了一种无人机与飞鸟目标雷达识别方法，首先对比两者目标在运动轨迹上差异性，进行特征分析，提出了时间相关的航向震荡频率与速度震荡频率特征量描述方法，并在离线状态下，利用实测雷达系统记录的航迹数据，提取两者的有效特征量，然后利用支持向量机算法对样本进行训练并在获得最优模型参数后，通过测试样本进行测试，测试分类结果显示能够达到较高的准确识别率，最后在线状态下通过大量的跟飞实验，其结果既表明该方法的正确性，也体现了该发明在工程实现角度上的轻量性、实用性、适用性，具有较高价值。具有较高价值。具有较高价值。


技术研发人员：冯正康 张良俊 曾正新 马小艳 黄盼盼 叶舟 赵婵娟 秦琨
受保护的技术使用者：上海航天电子通讯设备研究所
技术研发日：2023.06.29
技术公布日：2023/10/6