一种未知雷达目标脉冲序列的快速提取方法及系统与流程

本发明属于雷达脉冲流处理技术领域，具体涉及一种未知雷达目标脉冲序列的快速提取方法及系统。

背景技术：

在密集的雷达脉冲流中，快速筛选出雷达某种工作模式下的脉冲信号并获得尽量全面准确的信号参数，能赢取宝贵的时间，可以有针对性的快速对目标实施干扰，同时实时更新重点目标的参数、位置、活动情况、开关机规律等信息，便于迅速识别雷达目标的威胁等级和电磁环境态势，实现快速的威胁告警。

针对全脉冲数据中所含雷达信号种类未知的情况，从现有的雷达信号快速提取的相关方法，如基于自提取脉冲样本图的雷达信号快速提取方法，该方法是基于脉冲样本图，利用同种雷达信号脉冲之间的匹配相关性，通过脉冲的整体平移，在实现对脉冲样本图的自提取的同时，将可用此脉冲样本图描述的雷达辐射源信号筛选出来。其中，利用循环相关方法提取脉冲样本图，取出雷达信号单个脉冲的四维特征参数，即rf、pw、pri、mop。利用特征参数进行信号的匹配分选，参考图1，但该方法存在如下缺点：

1）在进行脉冲样本图自提取时，两脉冲的相似度测量要计算所有的特征参数，对于某一特征参数随机变化的雷达信号（如rf捷变类型），rf参数变化没有规律，对相似测度的准确性造成影响，降低了提取的准确度。

2）在进行两脉冲的匹配时，每次都计算所有的特征参数相似度，计算量大、耗时长。

技术实现要素：

本发明的目的之一，在于提供一种未知雷达目标脉冲序列的快速提取方法，该快速提取方法提取的脉冲样本图准确度高，对未知雷达目标信号参数的不同变化类型针对性强，且计算量小，运算效率高

本发明的目的之二，在于提供一种存储介质。

本发明的目的之三，在于提供一种未知雷达目标脉冲序列的快速提取系统。

为了达到上述目的之一，本发明采用如下技术方案实现：

一种未知雷达目标脉冲序列的快速提取方法，所述快速提取方法包括如下步骤：

步骤一、获取交错全脉冲序列；

步骤二、获取交错全脉冲序列中每个脉冲数据的信号到达方位角doa参数峰值并进行doa分块处理；

步骤三、获取每个doa分块内每个脉冲数据的脉冲宽度pw参数峰值并进行pw分块处理；

步骤四、获取每个pw分块内每个脉冲数据的载波频率rf参数峰值并进行rf分块处理，得到rf参数峰值对应的第一脉冲数据块和非rf参数峰值对应的第二脉冲数据块；

步骤五、判断未知雷达目标脉冲的rf参数值和rf参数峰值的差的绝对值是否小于阈值，如是，则对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理；如否，则对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理。

进一步的，步骤二的具体实现过程为：

步骤21、读取交错全脉冲序列中的脉冲数据并按照doa参数值进行排序，得到doa参数值的最大值doamax和最小值doamin；

步骤22、根据预先设定的doa值分辨率αgate以及doa参数值的最大值doamax和最小值doamin，计算doa峰值统计起始值doabegin、doa峰值统计结束值doaend以及doa峰值统计范围doawidth和doa峰值统计窗口数n_doawidth；

步骤23、设置doa峰值统计向量序列号i的初始值为1，第i个doa峰值统计向量peak_doa[i]的初始值为0；doa参数比较窗口αwindow的初始值为doabegin，doa峰值断点脉冲序号k的初始值为1；

步骤24、令从交错全脉冲序列中选取的脉冲序号j的初始值为1；

步骤25、判断i是否大于n_doawidth，如是，则进入步骤27；如否，则进入步骤26；

步骤26、当j<n且αwindow<x1(j)≤αwindow+αgate，则将peak_doa[i]加1后赋给peak_doa[i]，j加1后赋给j，返回步骤25；当j>n，则进入步骤27；当x1(j)>αwindow+αgate，则将k加1后赋给k，j加1后赋给j，αwindow加αgate后赋给αwindow，i加1后赋给i，返回步骤25；

其中，x1(j)为第j个脉冲对应的doa参数值；n为交错全脉冲序列中脉冲个数；

步骤27、统计doa峰值统计向量序列中doa参数值为波峰的脉冲个数m；并将交错全脉冲序列分成m个doa分块。

进一步的，步骤三的具体实现过程为：

步骤31、获取每个doa分块内脉冲个数n_doapid；并将每个doa分块内脉冲数据按照脉冲宽度pw参数值排序，得到对应doa分块内脉冲pw参数值的最大值pwmax和最小值pwmin；

步骤32、根据预先设定的pw值分辨率τgate以及每个doa分块内脉冲pw参数值的最大值pwmax和最小值pwmin，计算对应doa分块内pw峰值统计起始值pwbegin、pw峰值统计结束值pwend以及pw峰值统计范围pwwidth和pw峰值统计窗口数n_pwwidth；

步骤33、设置每个doa分块内pw峰值统计向量序列号i'的初始值为1，第i'个pw峰值统计向量peak_pw[i']的初始值为0；pw参数比较窗口τwindow的初始值为pwbegin，pw峰值断点脉冲序号k'的初始值为1；

步骤34、令对应doa分块内交错全脉冲序列中选取的脉冲序号j'的初始值为1；

步骤35、判断i'是否大于n_pwwidth，如是，则进入步骤37；如否，则进入步骤36；

步骤36、当j'<n_doapid且τwindow<x2(j')≤τwindow+τgate，则将peak_pw[i']加1后赋给peak_pw[i']，j'加1后赋给j'，返回步骤35；当j'>n_doapid，则进入步骤37；当x2(j')>τwindow+τgate，则将k'加1后赋给k'，j'加1后赋给j'，τwindow加上τgate后赋给τwindow，i'加1后赋给i'，返回步骤35；

其中，x2(j')为对应doa分块内第j'个脉冲对应的pw参数值；

步骤37、统计每个doa分块内pw峰值统计向量序列中pw参数值为波峰的脉冲个数l；并将对应doa分块内交错全脉冲序列中分成l个pw分块。

进一步的，步骤四的具体实现过程为：

步骤41、获取每个pw分块内脉冲个数n_pwpid；并将每个pw分块内脉冲数据按照脉冲宽度rf参数值排序，得到对应pw分块内脉冲rf参数值的最大值rfmax和最小值rfmin；

步骤42、根据预先设定的rf值分辨率ƒgate以及每个pw分块内脉冲rf参数值的最大值rfmax和最小值rfmin，计算对应pw分块内rf峰值统计起始值rfbegin、rf峰值统计结束值rfend以及rf峰值统计范围rfwidth和rf峰值统计窗口数n_rfwidth；

步骤43、设置每个pw分块内rf峰值统计向量序列号i''的初始值为1，第i''个rf峰值统计向量peak_pw[i'']的初始值为0；rf参数比较窗口ƒwindow的初始值为rfbegin，rf峰值断点脉冲序号k''的初始值为1；

步骤44、令对应pw分块内交错全脉冲序列中选取的脉冲序号j''的初始值为1；

步骤45、判断i''是否大于n_rfwidth，如是，则进入步骤47；如否，则进入步骤46；

步骤46、当j''<n_pwpid且ƒwindow<x3(j'')≤ƒwindow+ƒgate，则将peak_rf[i'']加1后赋给peak_rf[i'']，j''加1后赋给j''，返回步骤45；当j''>n_pwpid，则进入步骤37；当x3(j'')>ƒwindow+ƒgate，则将k''加1后赋给k''，j''加1后赋给j''，ƒwindow加上ƒgate后赋给ƒwindow，i''加1后赋给i''，返回步骤45；

其中，x3(j'')为对应pw分块内第j''个脉冲对应的rf参数值；

步骤47、将对应pw分块内rf参数值为波峰的脉冲数据归为第一脉冲数据块，剩余脉冲数据归为第二脉冲数据块。

进一步的，步骤五中，所述对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理的具体过程为：

步骤511、设置全脉冲移位位数i'''_shift的初始值为1，每一脉冲对应的脉冲匹配标志字的初始值均置零；

步骤512、令第一脉冲数据块中脉冲序列k'''的初始值为1，移位后的第一脉冲数据块中脉冲序列j'''的初始值为i'''_shift+1，设置移位对应的第一匹配脉冲个数number_matched的初始值为0；

步骤513、判断∣x3(k''')-x3(j''')∣是否小于等于rf值分辨率ƒgate，如是，进入步骤514；如否，则将j'''加1后赋给j'''，进入步骤515；

其中，x3(k''')和x3(j''')分别为第一脉冲数据块中第k'''个脉冲对应的rf参数值；x3(j''')为移位后的第一脉冲数据块中第j'''个脉冲对应的rf参数值；

步骤514、判断∣[x4(k''')-x4(1)]-[x4(j''')-x4(1+i'''_shift)]∣是否小于等于pri值分辨率tgate，如是，则将x4(k''')对应的脉冲匹配标示字置为1，并将number_matched加1后赋给number_matched，k'''加1后赋给k'''，j'''加1后赋给j'''，进入步骤515；如否，则将j'''加1后赋给j'''，进入步骤515；

其中，x4(k''')和x4(1)分别为第一脉冲数据块中第k'''和1个脉冲对应的pri参数值；x4(j''')和x4(1+i'''_shift)分别为移位后的第一脉冲数据块中第j'''和1+i'''_shift个脉冲对应的pri参数值；

步骤515、判断j'''是否小于n_number，如是，则返回步骤513；如否，则进入步骤516；

其中，n_number为第一脉冲数据块中脉冲总数目；

步骤516、判断number_matched是否大于number_threshold，如是，则此时对应的rf-pri投影样本图骨架周期为t_frame=x4(1+i'''_shift)-x4(1)，进入步骤517；如否，则将i'''_shift加1后赋给i'''_shift，返回步骤513；

其中，number_threshold为第一匹配脉冲个数门限值；

步骤517、提取脉冲匹配标志字置为1的脉冲，并将到达时间值最小的脉冲确定为首脉冲z(1)；

步骤518、以首脉冲z(1)为起始脉冲，以对应的rf-pri投影样本图骨架周期t_frame为时延进行脉冲匹配搜索。

进一步的，步骤五中，所述对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理的具体过程为：

步骤521、设置全脉冲移位位数i''''_shift的初始值为1，每一脉冲对应的脉冲匹配标志字的初始值均置零；

步骤522、令第二脉冲数据块中脉冲序列k''''的初始值为1，移位后的第二脉冲数据块中脉冲序列j''''的初始值为i''''_shift+1，设置移位对应的第二匹配脉冲个数number_matched'的初始值为0；

步骤523、判断∣[x4'(k'''')-x4'(1)]-[x4'(j'''')-x4'(1+i''''_shift)]∣是否小于等于pri值分辨率tgate，如是，则将x4'(k'''')对应的脉冲匹配标示字置为1，并将number_matched'加1后赋给number_matched'，k''''加1后赋给k''''，j''''加1后赋给j''''，进入步骤524；如否，则将j''''加1后赋给j''''，进入步骤524；

其中，x4'(k'''')和x4'(1)分别为第二脉冲数据块中第k''''和第1个脉冲对应的pri参数值；x4'(j'''')和x4'(1+i''''_shift)分别为移位后的第二脉冲数据块中第j''''和第1+i''''_shift个脉冲对应的pri参数值；

步骤524、判断j''''是否小于n_number'，如是，则返回步骤523；如否，则进入步骤525；

其中，n'_number为第一脉冲数据块中脉冲总数目；

步骤525、判断number_matched'是否大于number_threshold'，如是，则此时对应的pri投影样本图骨架周期为t_frame'=x4'(1+i''''_shift)-x4'(1)，进入步骤526；如否，则将i''''_shift加1后赋给i''''_shift，返回步骤523；

其中，number_threshold'为第二匹配脉冲个数门限值；

步骤526、提取脉冲匹配标志字置为1的脉冲，并将到达时间值最小的脉冲确定为首脉冲n(1)；

步骤527、以首脉冲n(1)为起始脉冲，以对应的pri投影样本图骨架周期t_frame'为时延进行未知雷达目标脉冲匹配搜索。

进一步的，步骤三中，采用多进程同时对多个doa分块进行pw分块处理。

进一步的，步骤四中，采用多进程同时对多个pw分块进行rf分块处理。

为了达到上述目的之二，本发明采用如下技术方案实现：

一种存储介质，所述存储介质存储有程序指令；所述程序指令被执行时，实现上述所述的快速提取方法。

为了达到上述目的之三，本发明采用如下技术方案实现：

一种未知雷达目标脉冲序列的快速提取系统，所述快速提取系统包括上述所述的存储介质。

本发明的有益效果：

本发明通过信号到达方位角doa参数峰值和脉冲宽度pw参数峰值实现对交错全脉冲序列doa分块、pw分块，从而剔除野值点（包括doa野值点和pw分块野值点）和杂声干扰；通过载波频率rf参数峰值，将块内数据中rf参数变化有规律的雷达目标脉冲（rf固定类型、rf跳变类型、rf脉组捷变类型，即第一脉冲数据块）与rf参数随机变化的雷达目标脉冲（rf脉间捷变类型，第二脉冲数据块）分开，一方面避免了将未统计出峰值的全脉冲序列剔除，从而导致rf脉间捷变类型目标信号丢失，造成漏警现象；二是将所有rf参数峰值对应脉冲分为一块，再进行后续处理，从而导致rf跳变、rf脉组捷变类型目标信号被分离成多个数据块，造成目标数据分批现象；并对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理，对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理，既减少了变化无规律的特征参数在匹配时的运算时间，降低了计算量，节省了处理时间，提高了处理速度，又避免了随机变化的特征参数在相似度测量时对匹配度值准确性的影响，提高了提取的准确度。

附图说明

图1为本发明的未知雷达目标脉冲序列的快速提取方法流程示意图；

图2为对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理流程示意图。

具体实施方式

以下结合附图对本发明的具体实施方式作出详细说明。

本实施例给出一种未知雷达目标脉冲序列的快速提取方法，参考图1，该快速提取方法包括如下步骤：

步骤一、获取交错全脉冲序列。

本实施例中的雷达脉冲描述字包括信号到达方位角（doa）参数、载波频率（rf）参数、脉冲宽度（pw）参数和到达时间（toa）参数。全脉冲序列数学模型y表示为：

；

其中，y为全脉冲序列数学模型；yn为第n个脉冲的特征参数向量，yn=[x1(n),x2(n),x3(n),x4(n)]，x1(n)~x4(n)分别为第n个脉冲对应的doa参数值、pw参数值、rf参数值和toa参数值；n=1,2,3,...,n，n为交错全脉冲序列的脉冲个数值。

步骤二、获取交错全脉冲序列中每个脉冲数据的信号到达方位角doa参数峰值并进行doa分块处理。

对doa参数峰值对应脉冲分块，假定峰值个数为m，则将全脉冲数据分为m块，对于未被分配到任一数据块的脉冲，则判定为野值点或干扰脉冲并剔除，剔除doa野值点和杂乱脉冲，可降低运算量，提高了运算效率，doa分块处理的具体实现过程为：

步骤21、读取交错全脉冲序列中的脉冲数据并按照doa参数值进行排序，得到doa参数值的最大值doamax和最小值doamin。

本实施例依据doa参数值从小到大的顺序，对全脉冲数据进行排序，得到doa参数值最大值和最小值分别为doamax和doamin。

当doa参数值在[-αlimit，αlimit]内时，则对应的全脉冲数据有效；反之，则剔除对应的全脉冲数据，αlimit为最大doa参数门限值，一般优选为30°。

步骤22、根据预先设定的doa值分辨率αgate以及doa参数值的最大值doamax和最小值doamin，计算doa峰值统计起始值doabegin、doa峰值统计结束值doaend以及doa峰值统计范围doawidth和doa峰值统计窗口数n_doawid。

本步骤中，doa峰值统计起始值doabegin和统计结束值doaend分别为：

doabegin=doamin-αgate-1（°）；

doaend=doamax+αgate+1（°）；

doa峰值统计范围为doawidth=doaend-doabegin；

doa峰值统计窗口数为n_doawidth=⌈doawidth/αgate⌉；

其中，doa值分辨率agate=0.5°，⌈doawidth/αgate⌉表示向上取整。

步骤23、步骤23、设置doa峰值统计向量序列号i的初始值为1，第i个doa峰值统计向量peak_doa[i]的初始值为0；doa参数比较窗口αwindow的初始值为doabegin，doa峰值断点脉冲序号k的初始值为1；

步骤24、令从交错全脉冲序列中选取的脉冲序号j的初始值为1；

步骤25、判断i是否大于n_doawidth，如是，则进入步骤27；如否，则进入步骤26；

步骤26、当j<n且αwindow<x1(j)≤αwindow+αgate，则将peak_doa[i]加1后赋给peak_doa[i]，j加1后赋给j，返回步骤25；当j>n，则进入步骤27；当x1(j)>αwindow+αgate，则将k加1后赋给k，j加1后赋给j，αwindow加αgate后赋给αwindow，i加1后赋给i，返回步骤25；

其中，x1(j)为第j个脉冲对应的doa参数值；n为交错全脉冲序列中脉冲个数。

步骤27、统计doa峰值统计向量序列中doa参数值为波峰的脉冲个数m；并将交错全脉冲序列分成m个doa分块。

本实施例的波峰条件为：peak_doa[i]大于peak_doa[i-1]，且peak_doa[i]大于等于peak_doa[i+1]；或当peak_doa[i]等于peak_doa[i-1]，且peak_doa[i]大于peak_doa[i+1]。

波谷的条件为：peak_doa[i]小于peak_doa[i-1]，且peak_doa[i]小于等于peak_doa[i+1]；或peak_doa[i]等于peak_doa[i-1]，且peak_doa[i]小于peak_doa[i+1]。

步骤三、获取每个doa分块内每个脉冲数据的脉冲宽度pw参数峰值并进行pw分块处理。

对pw参数峰值对应脉冲分块，假定峰值个数为l，则将全脉冲数据分为m块，对于未被分配到任一数据块的脉冲，则判定为野值点或干扰脉冲并剔除，剔除pw野值点和杂乱脉冲，可降低运算量，提高了运算效率，pw分块处理的具体实现过程为：

步骤31、获取每个doa分块内脉冲个数n_doapid；并将每个doa分块内脉冲数据按照脉冲宽度pw参数值排序，得到对应doa分块内脉冲pw参数值的最大值pwmax和最小值pwmin；

步骤32、根据预先设定的pw值分辨率τgate以及每个doa分块内脉冲pw参数值的最大值pwmax和最小值pwmin，计算对应doa分块内pw峰值统计起始值pwbegin、pw峰值统计结束值pwend以及pw峰值统计范围pwwidth和pw峰值统计窗口数n_pwwidth。

pwbegin=pwmin-τgate；

pwend=pwmax+τgate；

pwwidth=pwend-pwbegin；

n_pwwidth=⌈pwwidth/τgate⌉；

其中，pw值分辨率τgate=0.1µs，⌈pwwidth/τgate⌉表示向上取整。

步骤33、设置每个doa分块内pw峰值统计向量序列号i'的初始值为1，第i'个pw峰值统计向量peak_pw[i']的初始值为0；pw参数比较窗口τwindow的初始值为pwbegin，pw峰值断点脉冲序号k'的初始值为1；

步骤34、令对应doa块内交错全脉冲序列中选取的脉冲序号j'的初始值为1；

步骤35、判断i'是否大于n_pwwidth，如是，则进入步骤37；如否，则进入步骤36；

步骤36、当j'<n_doapid且τwindow<x2(j')≤τwindow+τgate，则将peak_pw[i']加1后赋给peak_pw[i']，j'加1后赋给j'，返回步骤35；当j'>n_doapid，则进入步骤37；当x2(j')>τwindow+τgate，则将k'加1后赋给k'，j'加1后赋给j'，τwindow加上τgate后赋给τwindow，i'加1后赋给i'，返回步骤35；其中，x2(j')为对应doa分块内第j'个脉冲对应的pw参数值；

步骤37、统计每个doa分块内pw峰值统计向量序列中pw参数值为波峰的脉冲个数l；并将对应doa分块内交错全脉冲序列中分成l个pw分块。

本实施例的波峰条件为：peak_pw[i']大于peak_pw[i'-1]，且peak_pw[i']大于等于peak_pw[i'+1]；或当peak_pw[i']等于peak_pw[i'-1]，且peak_pw[i']大于peak_pw[i'+1]。

波谷的条件为：peak_pw[i']小于peak_pw[i'-1]，且peak_pw[i']小于等于peak_pw[i'+1]；或peak_pw[i']等于peak_pw[i'-1]，且peak_pw[i']小于peak_pw[i'+1]。

步骤四、获取每个pw分块内每个脉冲数据的载波频率rf参数峰值并进行rf分块处理，得到rf参数峰值对应的第一脉冲数据块和非rf参数峰值对应的第二脉冲数据块。

本实施例通过rf参数峰值统计，将块内数据中rf参数变化有规律的雷达目标脉冲（rf固定类型、rf跳变类型、rf脉组捷变类型）与rf参数随机变化的雷达目标脉冲（rf脉间捷变类型）分开，一是可避免将未统计出峰值的全脉冲序列剔除，从而导致rf脉间捷变类型目标信号丢失，造成漏警现象；二是避免将每一个峰值对应脉冲分为一块，再进行后续处理，从而导致rf跳变、rf脉组捷变类型目标信号被分离成多个数据块，造成目标数据分批现象。rf分块处理的具体实现过程为：

步骤41、获取每个pw分块内脉冲个数n_pwpid；并将每个pw分块内脉冲数据按照脉冲宽度rf参数值排序，得到对应pw分块内脉冲rf参数值的最大值rfmax和最小值rfmin；

步骤42、根据预先设定的rf值分辨率ƒgate以及每个pw分块内脉冲rf参数值的最大值rfmax和最小值rfmin，计算对应pw分块内rf峰值统计起始值rfbegin、rf峰值统计结束值rfend以及rf峰值统计范围rfwidth和rf峰值统计窗口数n_rfwidth；

本实施例的rf值分辨率ƒgate=3mhz，为rf参数峰值统计时的分辨率，若两脉冲的rf值的差值小于ƒgate，则判定为rf值相同。

步骤43、设置每个pw分块内rf峰值统计向量序列号i''的初始值为1，第i''个rf峰值统计向量peak_pw[i'']的初始值为0；rf参数比较窗口ƒwindow的初始值为rfbegin，rf峰值断点脉冲序号k''的初始值为1；

步骤44、令对应pw块内交错全脉冲序列中选取的脉冲序号值j''的初始值为1；

步骤45、判断i''是否大于n_rfwidth，如是，则进入步骤47；如否，则进入步骤46；

步骤46、当j''<n_pwpid且ƒwindow<x3(j'')≤ƒwindow+ƒgate，则将peak_rf[i'']加1后赋给peak_rf[i'']，j''加1后赋给j''，返回步骤45；当j''>n_pwpid，则进入步骤37；当x3(j'')>ƒwindow+ƒgate，则将k''加1后赋给k''，j''加1后赋给j''，ƒwindow加上ƒgate后赋给ƒwindow，i''加1后赋给i''，返回步骤45；其中，x3(j'')为对应pw分块内第j''个脉冲对应的rf参数值；

步骤47、将对应pw分块内rf参数值为波峰的脉冲数据归为第一脉冲数据块，剩余脉冲数据归为第二脉冲数据块。

为达到了节约时间成本，加快处理速度，本实施例根据pw分块数l开辟子进程数为l，分别编号为子进程1，子进程2，…，子进程l后多个子进程同步进行rf分块处理。

步骤五、判断未知雷达目标脉冲的rf参数值和rf参数峰值的差的绝对值是否小于阈值，如是，则对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理；如否，则对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理。

参考图2，本实施例中对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理的具体过程为：

步骤511、设置全脉冲移位位数i'''_shift的初始值为1，每一脉冲对应的脉冲匹配标志字的初始值均置零；

步骤512、令第一脉冲数据块中脉冲序列k'''的初始值为1，移位后的第一脉冲数据块中脉冲序列j'''的初始值为i'''_shift+1，设置移位对应的第一匹配脉冲个数number_matched的初始值为0；

步骤513、判断∣x3(k''')-x3(j''')∣是否小于等于rf值分辨率ƒgate，如是，进入步骤514；如否，则将j'''加1后赋给j'''，进入步骤515；

其中，x3(k''')和x3(j''')分别为第一脉冲数据块中第k'''个脉冲对应的rf参数值；x3(j''')为移位后的第一脉冲数据块中第j'''个脉冲对应的rf参数值；

步骤514、判断∣[x4(k''')-x4(1)]-[x4(j''')-x4(1+i'''_shift)]∣是否小于等于pri值分辨率tgate，如是，则将x4(k''')对应的脉冲匹配标示字置为1，并将number_matched加1后赋给number_matched，k'''加1后赋给k'''，j'''加1后赋给j'''，进入步骤515；如否，则将j'''加1后赋给j'''，进入步骤515；

其中，x4(k''')和x4(1)分别为第一脉冲数据块中第k'''和1个脉冲对应的pri参数值；x4(j''')和x4(1+i'''_shift)分别为移位后的第一脉冲数据块中第j'''和1+i'''_shift个脉冲对应的pri参数值；

本实施例的pri值分辨率τgate=1µs。其中，雷达目标最大可能和最小可能的pri值分别为τlimit_max=20000µs和τlimit_min=5µs。

步骤515、判断j'''是否小于n_number，如是，则返回步骤513；如否，则进入步骤516；

其中，n_number为第一脉冲数据块中脉冲总数目；

步骤516、判断number_matched是否大于number_threshold，如是，则此时对应的rf-pri投影样本图骨架周期为t_frame=x4(1+i'''_shift)-x4(1)，进入步骤517；如否，则将i'''_shift加1后赋给i'''_shift，返回步骤513；

其中，number_threshold为第一匹配脉冲个数门限值；

步骤517、提取脉冲匹配标志字置为1的脉冲，并将到达时间值最小的脉冲确定为首脉冲z(1)；

步骤518、以首脉冲z(1)为起始脉冲，以对应的rf-pri投影样本图骨架周期t_frame为时延进行脉冲匹配搜索。

本实施例中，对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理的具体过程为：

步骤521、设置全脉冲移位位数i''''_shift的初始值为1，每一脉冲对应的脉冲匹配标志字的初始值均置零；

步骤522、令第二脉冲数据块中脉冲序列k''''的初始值为1，移位后的第二脉冲数据块中脉冲序列j''''的初始值为i''''_shift+1，设置移位对应的第二匹配脉冲个数number_matched'的初始值为0；

步骤523、判断∣[x4'(k'''')-x4'(1)]-[x4'(j'''')-x4'(1+i''''_shift)]∣是否小于等于pri值分辨率tgate，如是，则将x4'(k'''')对应的脉冲匹配标示字置为1，并将number_matched'加1后赋给number_matched'，k''''加1后赋给k''''，j''''加1后赋给j''''，进入步骤524；如否，则将j''''加1后赋给j''''，进入步骤524；

其中，x4'(k'''')和x4'(1)分别为第二脉冲数据块中第k''''和1个脉冲对应的pri参数值；x4'(j'''')和x4'(1+i''''_shift)分别为移位后的第二脉冲数据块中第j''''和1+i''''_shift个脉冲对应的pri参数值；

步骤524、判断j''''是否小于n_number'，如是，则返回步骤523；如否，则进入步骤525；

其中，n'_number为第一脉冲数据块中脉冲总数目；

步骤525、判断number_matched'是否大于number_threshold'，如是，则此时对应的pri投影样本图骨架周期为t_frame'=x4'(1+i''''_shift)-x4'(1)，进入步骤526；如否，则将i''''_shift加1后赋给i''''_shift，返回步骤523；

其中，number_threshold'为第二匹配脉冲个数门限值；

步骤526、提取脉冲匹配标志字置为1的脉冲，并将到达时间值最小的脉冲确定为首脉冲n(1)；

步骤527、以首脉冲n(1)为起始脉冲，以对应的pri投影样本图骨架周期t_frame'为时延进行未知雷达目标脉冲匹配搜索。

本实施例通过信号到达方位角doa参数峰值和脉冲宽度pw参数峰值实现对交错全脉冲序列doa分块、pw分块，从而剔除野值点（包括doa野值点和pw分块野值点）和杂声干扰；通过载波频率rf参数峰值，将块内数据中rf参数变化有规律的雷达目标脉冲（rf固定类型、rf跳变类型、rf脉组捷变类型，即第一脉冲数据块）与rf参数随机变化的雷达目标脉冲（rf脉间捷变类型，第二脉冲数据块）分开，一方面避免了将未统计出峰值的全脉冲序列剔除，从而导致rf脉间捷变类型目标信号丢失，造成漏警现象；二是将所有rf参数峰值对应脉冲分为一块，再进行后续处理，从而导致rf跳变、rf脉组捷变类型目标信号被分离成多个数据块，造成目标数据分批现象；并对第一脉冲数据块依次进行rf-pri投影样本图提取和未知雷达目标脉冲匹配处理，对第二脉冲数据块依次进行pri投影样本图提取和未知雷达目标脉冲匹配处理，既减少了变化无规律的特征参数在匹配时的运算时间，降低了计算量，节省了处理时间，提高了处理速度，又避免了随机变化的特征参数在相似度测量时对匹配度值准确性的影响，提高了提取的准确度。

另一实施例给出了一种存储介质，该存储介质存储有程序指令；所述程序指令被执行时，实现上述实施例给出的快速提取方法。

又一实施例给出了一种未知雷达目标脉冲序列的快速提取系统，所述快速提取系统包括上述实施例给出的存储介质。

以上实施方式仅用以说明本发明实施例的技术方案而非限制，尽管参照以上较佳实施方式对本发明实施例进行了详细说明，本领域的普通技术人员应当理解，可以对本发明实施例的技术方案进行修改或等同替换都不应脱离本发明实施例的技术方案的精神和范围。