基于深度学习的高速机动目标检测方法与流程

本发明属于雷达技术领域，涉及目标检测方法，可用于对高速机动目标进行检测。

背景技术：

航空航天技术的不断发展使得越来越多的高速、机动目标出现在雷达探测领域，对于此类目标，通过传统的脉冲压缩和动目标检测mtd方法提高雷达探测性能将面临如下问题：一方面是目标高速运动带来的回波脉内多普勒调制，导致脉冲压缩时有较高的信噪比损失；另一方面是在相参积累时需要校正目标速度引起的一阶距离走动和加速度引起的距离弯曲rc，补偿目标加速度引起的多普勒走动dfm。针对距离走动的问题，r.tao研究了做加速运动的高速目标相参积累过程中的距离走动与多普勒走动效应，分析了目标发生距离走动与多普勒走动的临界条件，并提出了基于缩放处理和分数阶傅里叶变换spfrft的距离走动校正和多普勒走动补偿算法。j.su针对高速目标相参积累检测问题，提出了基于kt和去调频处理的校正补偿算法实现了目标能量的相参积累。m.xing提出了基于kt与最小熵准则的高速机动目标相参积累与运动参数估计算法，并分析了运动参数估计误差对相参积累性能的影响。c.pang研究了高信噪比和低信噪比两种情况下的距离走动校正与多普勒走动补偿算法，指出高信噪比时可以先采用二维中值滤波校正距离走动，然后通过瞬时互相关方法估计出目标的加速度并补偿多普勒走动；低信噪比下则还需要结合分数阶傅里叶变换与离散多项式变换。在ar算法的基础上，x.rao提出了基于改进ar与分数阶傅里叶变换iar-frft的匀加速目标相参积累算法。

上述算法由于对回波信号的输入信噪比要求较高，不适用于低信噪比下的目标相参积累；并且由于具有较高的计算复杂度，不利于工程实现。

技术实现要素：

本发明的目的在于克服上述已有技术的不足，提出一种基于深度学习的高速机动目标检测方法，以减小由目标速度引起的脉冲压缩信噪比损失，克服距离走动导致的相参积累困难，提高在低信噪比条件下的高速机动目标检测性能；并降低计算复杂度，以利于工程实现。

实现本发明目的的技术方案，包括如下：

1)构造单脉冲检测网络训练数据集，对单脉冲检测网络进行训练：

1a)对所有的雷达回波数据按照快时间沿行排列，慢时间沿列排列得到m个行数为g、列数为p的复矩阵{um}m＝1,...,m，g≥1，p≥1,并以滑窗的形式对每个复矩阵进行分解得到n个长度为l的复数向量{en}n＝1,...,n；

1b)将每个复数向量en转换为归一化特征累积矢量qn，再将qn的实部和虚部分别作为矩阵的第一行和第二行，并将该矩阵记为xn；假设已知回波数据中所有目标所在的距离门，根据en的第一个距离单元内是否存在目标产生对应的标号yn，即yn取1表示en的第一个距离单元内有目标，取0则表示无目标；将集合{xn,yn}n＝1,...,n作为单脉冲检测网络的训练数据集；

1c)使用一维残差网络作为单脉冲检测网络，基于训练数据集{xn,yn}n＝1,...,n，使用交叉熵代价函数和反向传播算法对单脉冲检测网络进行训练，直到代价函数收敛时，得到训练好的单脉冲检测网络；

2)构造融合网络训练数据集,对融合网络进行训练：

2a)将雷达回波数据um输入训练好的单脉冲检测网络，输出为一个g行r列的矩阵ym，并确定与ym对应的标号zm；

2b)用ym及zm组成集合{ym,zm}m＝1,...,m，并将该集合作为融合网络的训练数据集；

2c)使用深度卷积网络作为融合网络，基于训练数据集{ym,zm}m＝1,...,m，采用交叉熵代价函数和反向传播算法对融合网络进行训练，直到代价函数收敛，得到训练好的融合网络；

3)使用训练好的单脉冲检测网络和训练好的融合网络对待检测的回波数据进行目标检测：

3a)将待检测的雷达回波数据按照快时间沿行排列，慢时间沿列排列的方法表示为一个行数为g、列数为p的矩阵utst；将矩阵utst输入训练好的单脉冲检测网络做滑窗检测，得到单脉冲检测结果矩阵ytst；

3b)将单脉冲检测结果矩阵ytst输入训练好的融合网络，得到网络输出btst，该btst是一个长度为r的实数向量，其第r个元素btst(r)表示第r个距离门含有目标的概率，1≤r≤r，若btst(r)≥0.5，则表示第r个距离门有目标,否则，该距离门没有目标，至此完成对回波数据的目标检测。

本发明与现有技术相比，具有如下优点：

1)本发明使用单脉冲检测网络对回波数据进行滑窗检测，信噪比损失小于传统的脉冲压缩方法；

2)本发明使用融合网络对单脉冲检测结果做融合处理，目标所有的距离走动引起的能量扩散都在融合网络感受野内，由于综合考虑所有距离单元上的目标信息，相比于传统的动目标检测方法，具有更好的检测性能；

3)本发明的整个检测流程基于神经网络方法实现，通过gpu平台进行同质运算，有更高的计算效率，利于工程实现。

附图说明

图1是本发明的实现流程图；

图2是本发明中融合网络的结构图；

图3是本发明中融合网络的降维生成层示意图；

图4是用本发明对单脉冲的检测曲线图；

图5是用本发明对多脉冲的检测曲线图。

具体实施方式

下面结合附图，对本发明的实施例及效果做详细说明：

参照图1，本实例的实现步骤如下：

步骤1，构造单脉冲检测网络训练数据集。

本实例针对单通道雷达，假设该雷达发射相位编码波形，发射波形为s，在发射脉宽内的采样点数为l，则s是一个长度为l的恒模复数向量，该雷达检测目标的运动模型都是匀加速运动，其速度范围是[vmin,vmax]，加速度的范围是[amin,amax]，其训练数据集的构造如下：

1a)对所有的雷达回波数据按照快时间沿行排列，慢时间沿列排列，得到m个行数为g、列数为p的复矩阵{um}m＝1,...,m，g≥1，p≥1；

1b)以滑窗的形式对每个复矩阵进行分解得到n个长度为l的复数向量{en}n＝1,...,n，其中，每个向量en通过如下公式得到：

en＝um(g,r:r+l-1)，1≤g≤g，1≤r≤r

式中，um(g,r:r+l-1)表示回波数据矩阵中第g行，第r列到第r+l-1列元素；l表示滑窗长度；r表示对每一行进行滑窗分解的次数；en表示滑窗分解得到的数据，是长度为l的复向量；

1c)将每个复数向量en按如下公式转换为归一化特征累积矢量qn：

其中s^*表示发射波形s的共轭；表示hadamard积；cumsum(·)表示累加操作；norm(·)表示向量2范数；max{·}表示在实数向量中取最大值；

1d)将qn的实部和虚部分别作为矩阵的第一行和第二行，并将该矩阵记为xn；

xn(1,:)＝real(qn)

xn(2,:)＝imag(qn)

real(·)表示对复数取实部；imag(·)表示对复数取虚部；xn表示第n个单脉冲检测网络的输入数据，是一个行数为2，列数为l的实数矩阵；xn(1,:)表示xn的第一行元素；xn(2,:)表示xn的第二行元素；

1e)假设已知回波数据中所有目标所在的距离门，根据en的第一个距离单元内是否存在目标产生对应的标号yn：

yn＝indicator(en的第一个距离门内包含目标)，

其中，en表示滑窗分解得到的数据，是一长度为l的复向量；indicator(·)表示指示函数，当en的第一个距离单元内有目标时，则yn取值为1；当en的第一个距离单元内无目标时，则yn取值为0；

1f)将矩阵xn和与之对应的标号yn组成集合{xn,yn}n＝1,...,n，用该集合{xn,yn}n＝1,...,n作为单脉冲检测网络的训练数据集。

步骤2，对单脉冲检测网络进行训练。

2a)使用一维残差网络作为单脉冲检测网络；

2b)基于训练数据集{xn,yn}n＝1,...,n，使用交叉熵代价函数和反向传播算法对单脉冲检测网络进行训练，直到代价函数收敛时，得到训练好的单脉冲检测网络；

所述交叉熵代价函数，其表达式如下：

其中，c表示交叉熵，n表示样本总数，x表示样本输入，y表示样本输出，a表示神经元实际输出；

步骤3，构造融合网络训练数据集。

3a)将每个雷达回波数据um输入训练好的单脉冲检测网络，输出一个g行r列的矩阵ym，该矩阵ym中的元素ym(g,r)表示回波数据第g个脉冲第r个距离门含有目标的概率，1≤g≤g，1≤r≤r；

3b)生成与ym对应的标号zm，zm是一个长度为l的向量，其第r个元素zm(r)按如下公式产生：

zm(r)＝indicator(um(1,r)中包含目标)，

其中，indicator(·)表示指示函数，当雷达回波数据um第一个脉冲回波的第r个距离门含有目标时，zm(r)取值为1；当雷达回波数据um第一个脉冲回波的第r个距离门不含目标时，zm(r)取值为0；

3c)用ym及zm组成集合{ym,zm}m＝1,...,m，并将该集合作为融合网络的训练数据集。

步骤4,对融合网络进行训练。

4a)使用深度卷积网络作为融合网络，该融合网络结构如图2所示：其输入层和输出层均为卷积层，隐藏层共6层，其中前5个隐藏层为降维生成层，如图3所示，第6个隐藏层为卷积层；

4b)基于训练数据集{ym,zm}m＝1,...,m，采用交叉熵代价函数和反向传播算法对融合网络进行训练，直到代价函数收敛，得到训练好的融合网络。

步骤5,使用训练好的单脉冲检测网络和训练好的融合网络对待检测的回波数据进行目标检测。

5a)将待检测的雷达回波数据按照快时间沿行排列，慢时间沿列排列的方法表示为一个行数为g、列数为p的矩阵utst；将矩阵utst输入训练好的单脉冲检测网络做滑窗检测，得到单脉冲检测结果矩阵ytst；

5b)将单脉冲检测结果矩阵ytst输入训练好的融合网络，得到网络输出btst，该btst是一个长度为r的实数向量，其第r个元素btst(r)表示第r个距离门含有目标的概率，1≤r≤r；

5c)根据虚警率设置检测门限λ＝0.5，将第r个距离门含有目标的概率btst(r)与检测门限λ比较：若btst(r)≥λ，则表示第r个距离门有目标,否则，该距离门没有目标，至此完成对回波数据的目标检测。

本发明的效果可以通过以下仿真进一步说明

1.仿真参数

使用一维残差网络作为单脉冲检测网络，其参数如表1所示；使用图2所示的深度卷积网络作为融合网络，其参数如表2所示；设定虚警率为10^-4。

表1.单脉冲检测网络结构

表2.融合网络结构

2.仿真内容

仿真1，使用交叉熵代价函数和反向传播算法对单脉冲检测网络进行训练，并对其性能进行测试，绘制单脉冲检测网络检测概率与信噪比的关系图，结果如图4所示。从图4可以发现，当参考单元数为20，使用相位编码信号或线性调频信号的单元平均-恒虚警ca-cfar检测效果都劣于单脉冲检测网络的检测效果。

仿真2，对比本发明和ca-cfar检测方法的检测概率与信噪比曲线，ca-cfar检测方法中使用相位编码信号，参考单元数为20，相参积累方法分别使用无keystone变换和keystone变换的积累方法，其结果如图5所示。从图5可以发现，虽然本发明在信噪比过低时检测效果不如ca-cfar的检测效果，但是随着信噪比增加，本发明的检测性能快速提升，总体检测效果优于ca-cfar的检测效果。