一种基于同构序列的脉冲分选方法
【专利摘要】本发明提出一种基于同构序列搜索的脉冲分选方法，解决实际中淹没于大量脉冲流里的脉组检测问题。步骤S1，由脉冲流的到达时间序列的一阶后向差分得到脉冲重复时间间隔序列，并对脉冲重复时间间隔序列的每个元素的值作量化处理；步骤S2，利用后缀数组和最大公共前缀的方法，筛选脉冲重复时间间隔序列的重复子串：步骤S3，对于有包含关系的子串，删除长度较短的子串，对于有交叠关系的子串，进行合并衔接操作；步骤S4，对剩余脉冲流，构造脉冲流到达时间差值矩阵；步骤S5，提取差值矩阵每一行的正实数序列组成一维数组，进行排序得到若干子集；步骤S6，搜索每个子集的和的最大公共子序列，确定目标脉冲的位置；步骤S7，谐波校验，脉冲丢失校验。
【专利说明】—种基于同构序列的脉冲分选方法
【技术领域】
[0001]本发明涉及一种基于同构序列的脉冲分选方法，属于信号处理领域。
【背景技术】
[0002]随着电磁环境的日益复杂化，如何从高度密集的脉冲信号环境中分选出感兴趣的信号，挖掘其中的规律，一直都是雷达、声呐等领域的研究热点。提取有规律的脉组序列是脉冲分选的重要方法。早期脉冲分选的经典方法是对脉冲参数(PRI)进行统计直方图分析，直方图法主要用于分析常规脉冲信号，无法消除谐波的影响，更无法适应现代复杂多变的脉冲类型。PRI变换算法是利用了相位因子达到抑制谐波的作用，但该方法对参数初值设置较为敏感，且计算量大。状态转移矩阵法从研究相邻两个PRI值的关联性的角度，分析整个脉冲序列的PRI特性，由此得到的分选结果有局限性。
[0003]直方图、PRI变换、状态转移矩阵法的作用都是分析PRI的整体规律，还需要借助别的方法把符合该规律的脉冲提取出来。而且它们的分选效果在一定程度上依赖于目标脉冲的累加效果，若目标脉冲占整个脉冲流的比例很小，则目标脉冲的PRI特性将不能被检测出。而实际上，由于接收系统采用空域、频域二维方式搜索目标，接收机在一定时间内以一定的概率对目标脉冲的录取，或者由于辐射源自身扫描和灵活的波束调度，波束驻留时间很短，这些原因都导致了接收机只能获取部分脉冲流。辐射源发出的脉冲多成组工作，因此接收到的大量脉冲流里面实际上只存在少量的目标辐射源的脉冲组。这些少量脉组的差值直方图、PRI谱峰都很低，且状态转移概率也很低，无法明显地区分于干扰脉冲。若这些脉组还发生了脉冲丢失，那上述算法的检测性能会更低。
[0004]因此，不仅为了更快地处理普通情形下的脉冲分选，还要解决这种淹没于大量脉冲流里的脉冲组检测问题，同时要求对脉冲组发生脉冲丢失或受到干扰的情形有良好的鲁棒性，我们提出了基于同构序列的脉冲分选方法，所述同构脉冲序列，即脉冲序列的PRI值相同(或近似相同)，如果这些脉冲序列没有发生脉冲丢失或没有受到脉冲干扰，那么称之为显同构序列；如果同一辐射源发出的脉冲序列存在脉冲丢失或脉冲干扰，或多个辐射源发出的脉冲序列交织在一起，那么称这种脉冲序列为隐同构序列。
【发明内容】
[0005]本发明是为了克服现有技术的缺陷，除了能够更快地处理普通情况下的脉冲分选之外，还要解决实际中淹没于大量脉冲流里的脉组检测问题，同时为了进一步解决当脉组受到干扰或发生脉冲丢失情形下的检测问题，本发明提出了一种基于同构序列的脉冲分选方法。
[0006]本发明提出的一种基于同构序列的脉冲分选方法，该方法包括以下步骤:
[0007]步骤SI，由脉冲流的到达时间序列的一阶后向差分得到脉冲重复时间间隔序列，并对脉冲重复时间间隔序列的每个元素的值作量化处理；
[0008]步骤S2，利用后缀数组和最大公共前缀的方法，筛选上述脉冲重复时间间隔序列的重复子串:
[0009]步骤S3，对于有包含关系的子串，删除长度较短的子串，对于有交叠关系的子串，进行合并衔接操作；
[0010]步骤S4，对剩余脉冲流，构造脉冲流到达时间差值矩阵；
[0011]步骤S5，提取上述差值矩阵每一行的正实数序列组成一维数组Φ，进行排序得到若干子集{Φρ Φ2，...，Φ,}，每个子集内元素相同，对应一个脉冲重复时间间隔值，记录每个子集Φ?的元素在差值矩阵中对应的位置，子集每个元素行序号构成一维数组Rowi,列序号构成一维数组Columni ；
[0012]步骤S6,搜索每个子集Φ i对应的Rowi和Columni的最大公共子序列,确定目标脉冲的位置；
[0013]步骤S7，谐波校验:如果子集对应的脉冲重复时间间隔是子集Φ?对应的脉冲重复时间间隔的整数倍，且子集I对应的脉冲包含子集Φ]对应的脉冲，则认为子集对应的脉冲序列是谐波影响构成的；脉冲丢失校验:对于某子集Φ?中搜索到的目标脉冲，如果每个脉冲的列序号等于相邻后一个脉冲的行序号，则不存在脉冲丢失现象；若存在某个脉冲的列序号不等于相邻后一个脉冲的行序号，则存在脉冲丢失现象；至此完成脉冲分选。
[0014]有益效果:本发明提出一种基于同构序列的脉冲分选方法，把在大量脉冲流中提取目标辐射源脉组的问题转化成了提取PRI序列的重复子串问题，把后缀数组和最大公共前缀、最长公共子序列的概念引入脉冲分选，实现了 PRI重复子串的快速筛选，若脉组发生了脉冲丢失或受到干扰，把提取目标脉冲序列转化为了在所构造的TOA差值矩阵中搜索真实PRI值的问题，可有效处理脉冲干扰和丢失情形，借助容差的设置，该方法对抖动、参差重频类型的脉冲有一定的检测能力。
【专利附图】
【附图说明】
[0015]图1是本发明的流程示意图；
[0016]图2是本发明的脉冲序列的到达时间和载频二维分布图；
[0017]图3是本发明的搜索到辐射源目标脉冲的到达时间和载频的二维分布图；
[0018]图4是本发明的脉冲序列的到达时间和载频二维分布图；
[0019]图5是本发明的搜索到的辐射源目标脉冲的到达时间和载频的二维分布图。
【具体实施方式】
[0020]为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
[0021]本发明提供的方法既可以在个人计算机、工控机及服务器上以软件的形式安装并执行，也可将方法做成嵌入式芯片以硬件的形式来体现。
[0022]本发明所提出的基于同构序列的脉冲分选方法包括以下步骤:
[0023]步骤SI，由脉冲流的到达时间(TOA)序列的一阶后向差分得到脉冲重复时间间隔(PRI)序列，并对此PRI序列的每个元素的值作量化处理；
[0024]所述脉冲流的到达时间序列表示为:
[0025]TOA = It1, t2,..., tn},其中&表示η个脉冲中第i个脉冲的前沿到达时间，i =1，…，η ;所述脉冲流的重复时间间隔序列表示为:PRI = (PRI1, PRI2,PRRn_J，其中PRIj=tJ+1-tj表示第j个脉冲与第j+Ι个脉冲的时间间隔，j = I,…，n-1。
[0026]所述的PRI量化处理可以减小后续处理的压力，把要考量的PRI的范围划分成若干个子区间，落在同一子区间的PRI值用同一个值代替，区间的长度取决于实际测量的精度。
[0027]步骤S2，利用后缀数组和最大公共前缀的方法，筛选上述PRI序列的重复子串:
[0028]所述步骤S2进一步包括以下步骤:
[0029]步骤S21，构造PRI序列的后缀，并进行排序得到后缀数组SA:
[0030]后缀:从长度为η-1的PRI序列的某个位置i开始到序列末尾的一个特殊子串，记为 Sufix (i), I ^ i ^ η-1 ；
[0031]后缀数组:PRI序列的所有后缀都按字典顺序进行排序后，将后缀在PRI序列中的位置用一个一维数组记录，这个数组称为后缀数组，记为SA,且Sufix(SA[i])〈Sufix(SA[i+l])，I < i〈n-l,
[0032]长度为n-1的PRI序列有n-1个后缀,可采用倍增算法、Difference Cove mod3等算法对这η-1个后缀进行排序，构造后缀数组SA。
[0033]步骤S22:比较相邻后缀，提取它们的最长公共前缀和对应的SA值:
[0034]字典比较:比较两个串的大小，先比较第一个元素的大小，如果相等再比较第二个元素大小，依次类推；
[0035]前缀:从PRI序列的开头到某个位置i的一个特殊子串，记为Prefix (i), I ^ i ^ n-1,
[0036]提取相邻后缀的最大公共前缀，由于量化误差、测量误差的存在，要求对应元素之差在一定的阈值内，阈值的大小与脉冲重频类型有关。
[0037]步骤S3，对于有包含关系的子串，删除长度较短的子串，对于有交叠关系的子串，对它们进行合并衔接操作:
[0038]所述步骤S3进一步包括以下步骤:
[0039]S31:按照SA顺序，对搜索到的子串进行排序；
[0040]排序后的子串可表示为Substringa)，长度记为length(i)，该子串的首个元素在整个PRI序列里的位置index (i)；
[0041]S32:若两个相邻的子串存在包含关系，则删除其中长度小的子串；若两个相邻的子串存在交置关系，则合并衔接这两个子串:
[0042]包含、交叠关系可以描述为:
[0043]包含:若两个相邻子串分别为Substring (i)，Substring (i+Ι),长度分另丨J 为 length (i), length (i+1),索引值分别为 index (i),index (i+1),不妨设length (i)〈length (i+1),若满足以下条件,贝U删除 Substring (i)，
「 ? \SA(j)>SA{/+\)
[0044]i
ISA{f) + lcngth^f) — I S SA{i +1) + length^ + 1)-1
[0045]交叠:若两个相邻子串满足下面的条件，
【权利要求】
1.一种基于同构序列的脉冲分选方法，其特征在于，该方法包括以下步骤: 步骤Si，由脉冲流的到达时间序列的一阶后向差分得到脉冲重复时间间隔序列，并对脉冲重复时间间隔序列的每个元素的值作量化处理； 步骤S2，利用后缀数组和最大公共前缀的方法，筛选上述脉冲重复时间间隔序列的重复子串: 步骤S3，对于有包含关系的子串，删除长度较短的子串，对于有交叠关系的子串，进行合并衔接操作； 步骤S4，对剩余脉冲流，构造脉冲流到达时间差值矩阵； 步骤S5，提取上述差值矩阵每一行的正实数序列组成一维数组Φ，进行排序得到若干子集{Φρ Φ2，...，Φ,}，每个子集内元素相同，对应一个脉冲重复时间间隔值，记录每个子集Φ?的元素在差值矩阵中对应的位置，子集每个元素行序号构成一维数组Rowi,列序号构成一维数组Columni ； 步骤S6,搜索每个子集Φ i对应的Rowi和Columni的最大公共子序列,确定目标脉冲的位置； 步骤S7，谐波校验:如果子集Φ]对应的脉冲重复时间间隔是子集(^对应的脉冲重复时间间隔的整数倍，且子集Φ i对应的脉冲包含子集Φ j对应的脉冲，则认为子集Φ j对应的脉冲序列是谐波影响构成的；脉冲丢失校验:对于某子集Φ?中搜索到的目标脉冲，如果每个脉冲的列序号等于相邻后一个脉冲的行序号，则不存在脉冲丢失现象；若存在某个脉冲的列序号不等于相邻后一个脉冲的行序号，则存在脉冲丢失现象；至此完成脉冲分选。
2.如权利要求1所述的一种基于同构序列的脉冲分选方法，其特征在于，所述步骤SI中的量化处理为将要考量的PRI的范围划分成若干个子区间，落在同一子区间的PRI值用同一个值代替，区间的长度取决于实际测量的精度。
3.如权利要求1所述的一种基于同构序列的脉冲分选方法，其特征在于，所述步骤S2的筛选脉冲重复时间间隔序列的重复子串进一步包括以下步骤: 步骤S21:构造PRI序列的后缀，并进行排序得到后缀数组SA ； 步骤S22:比较相邻后缀，提取它们的最长公共前缀和对应的SA值。
4.如权利要求3所述的一种基于同构序列的脉冲分选方法，其特征在于:所述步骤S21的后缀、后缀数组描述为: 后缀:从长度为η的序列S的某个位置i开始到序列末尾的一个特殊子串，记为Sufix (i), I ^ i ^ η ； 后缀数组:序列S的所有后缀都按字典顺序进行排序后，将后缀在序列S中的位置用一个一维数组记录，这个数组称为后缀数组，记为SA,且Sufix(SA[i])〈Sufix(SA[i+l])，I ( i<η。
【文档编号】G01S7/292GK104198999SQ201410386459
【公开日】2014年12月10日 申请日期:2014年8月7日 优先权日:2014年8月7日
【发明者】任丽香, 李枫, 曹军, 徐嘉岑, 龙腾 申请人:北京理工大学