参差序列及脉组参差序列的PRI子周期提取方法与流程

本发明属于数据周期提取
技术领域：
，特别是涉及一种参差序列及脉组参差序列的pri子周期提取方法。
背景技术：
：近几年，随着各种周期信号的周期估计在各个行业领域的迅猛发展，其技术方法逐渐应用于军事和民用领域，特别是电子战领域的应用更为显著，多用于电子干扰与反干扰的脉冲重复间隔(pri)估计，经过几十年电子侦察的逐渐发展，脉冲重复间隔估计技术也越来越成熟，发明了序列搜索匹配和功率谱估计等一系列方法及改进型，但是随着电磁环境的复杂化，脉冲序列越来越复杂，脉冲重复间隔估计也面临着越来越大的挑战，如多目标辐射源造成脉冲序列混叠、人为因素干扰和测量装置不准确所造成的杂散脉冲、缺失脉冲和脉冲抖动等。辐射源覆盖范围不断变宽会导致nyquist采样频谱增加，从而增加了测量装置采样和数据存储、处理的难度，因此，传统估计方法受到了越来越大的挑战和限制。近十年，压缩感知理论的逐渐成熟和广泛应用，给我们带来了一种新的思路，即将压缩感知技术应用于脉冲重复间隔的估计。稀疏重构方法只能提取隐藏周期互不相等的脉冲序列周期，然而对于由多个相同隐藏周期脉冲序列组成的混叠脉冲序列，该混叠脉冲序列被称为参差脉冲序列。参差脉冲序列含有多个子周期，所有子周期的和称为框架周期。在出现较多的虚假脉冲和缺失脉冲的情况下，利用传统算法，如非脉冲分选序列搜索的方式，估计参差脉冲序列子周期的效果并不理想。这些传统方法往往只能得到参差脉冲序列的框架周期，而不能得到其子周期。同时，如果存在大量缺失脉冲，造成参差脉冲序列和脉组参差序列并没有完整的周期结构，这会给传统方法造成了极大的困难。传统的参差脉冲序列搜索方法主要是利用两个完整的周期结构，通过序列搜索来最终达到估计参差脉冲序列子周期的目的。如果脉冲序列中存在非完整的周期结构，即组成每个周期结构的框架周期是非完整的，参差脉冲序列搜索方法并不能准确地估计子周期。同时，在出现较多的虚假脉冲和缺失脉冲的情况下，此种方法的估计性能迅速下降，甚至失效。复杂电磁环境下的多目标混叠脉冲序列，其复杂性来源于多种样式的目标序列的混叠。例如各个周期不等的目标序列相混叠后，构成了一种具有子周期相乘结构的脉冲序列；各个周期相等的目标序列相混叠构成的参差脉冲序列或脉组参差序列；还有由其它样式构成的滑变脉冲序列等。对于这些各种类型的脉冲序列，传统算法很难提取出这些脉冲序列的子周期。现有技术中存在以下问题：1)对于一个由相同隐藏周期组成的混叠脉冲序列，即参差脉冲序列，传统方法很难从混叠脉冲序列中提取出序列的pri子周期；2)对于不完整的参差脉冲序列，连续的虚假脉冲和不完整的离散缺失都会导致传统估计方法的性能迅速下降；3)传统算法的适应范围较窄，能够估计的脉冲类型十分有限，不能提取加性结构和乘性结构的脉冲序列的pri子周期。因此，亟需一种参差脉冲序列及脉组参差脉冲序列的pri子周期提取方法。以解决现有技术中存在的难于从混叠脉冲序列中提取出序列的pri子周期、不完整的参差脉冲序列导致算法性能下降和适应范围较窄的问题，实现在连续的虚假脉冲和不完整的离散缺失情况下，准确提取各种混叠脉冲序列的pri子周期的目的。技术实现要素：本发明的目的在于提供一种参差序列及脉组参差序列的pri子周期提取方法，以解决现有技术中存在的难于从混叠脉冲序列中提取出序列的pri子周期、不完整的参差脉冲序列导致算法性能下降和适应范围较窄的问题，实现在连续的虚假脉冲和不完整的离散缺失情况下，准确提取各种混叠脉冲序列的pri子周期的目的。本发明所采用的技术方案是，参差序列及脉组参差序列的pri子周期提取方法，包括以下步骤：s1，构建固定重频脉冲序列的初始相位匹配库：设固定重频脉冲序列的周期为pg，则初始相位匹配库有pg行，每行有φ(pg)个元素，所述φ(pg)为欧拉函数；由于周期为pg的ramanujan子空间内有φ(pg)个列，则初始相位匹配库的第一至第φ(pg)行中，存在一个具有正单位幅值的正强度元素，而其余元素的幅值为零；初始相位匹配库的第φ(pg)+1行至第pg行的元素由下列公式确定：其中，dj是周期pg的第j个除数；si为初始相位匹配库中第i行的元素集，i＝1,2,…,pg；p为参差脉冲序列的框架周期；si中前φ(pg)行的元素集，即i＝1,2,…,φ(pg)，是已知的；利用公式(1)计算得出后pg-φ(pg)行的元素集，即i＝φ(pg)+1,…,pg，其过程为：首先计算出周期pg的所有除数d，即pg除以d得到一个整数；根据除数d的大小，依照从大到小的顺序循环计算公式(1)，直至算出初始相位匹配库的第φ(pg)+1行至第pg行的所有元素；s2，构建参差脉冲序列的初始相位匹配库，具体包括以下步骤：s21，构建含有单个初始相位的参差脉冲序列的初始相位匹配库k1：首先计算出周期为p的所有除数d，即p除以d得到一个整数；再将周期为p以及所有除数d的固定重频脉冲序列初始相位匹配库，根据从大到小的顺序整合到一个行数为p的参差脉冲序列初始相位匹配库k1中；s22，构建含有多个初始相位的参差脉冲序列的初始相位匹配库km：参差脉冲序列的框架周期为p由m个不相同初始相位的固定重频脉冲序列混叠组成，即含有m个不相同的初始相位；从所述s21中初始相位匹配库k1中任意提取出m个行，并且将其相应的元素相加作为初始相位匹配库km的一行，得到含有nu个行的初始相位匹配库km；所述nu由下列公式求得：其中，m为多个初始相位的参差脉冲序列中初始相位的个数；s3，利用多数匹配方法并结合ramanujan结构字典的稀疏恢复，确定参差脉冲序列pri子周期，即利用多数匹配方法对集合中的每个元素进行计算，集合中含有多少个元素，就进行多少次匹配计算。进一步的，所述s3中多数匹配方法包括以下步骤：s30，计算除数集合：设一个初始相位数小于mm的参差脉冲序列，其框架周期为p，这里mm为正整数，是初始相位数的上限，并且1＜mm≤p；首先计算由p的除数组成的集合ds，即：ds＝{d1,d2,…,dn}，p/di为整数；di≠1；i＝1,2,…,n；s31，计算周期稀疏解：从集合ds取出一个元素di，构成一个以此元素值为周期的ramanujan结构字典a，求解以下优化问题：min||dx||2s.t.y＝φx；φ＝ad(3)其中，y代表测量的脉冲序列，x代表周期稀疏解，d表示权重矩阵；min||dx||2表示计算一个矢量的2-范数的最小值；φ代表一个中间矩阵；优化问题的最优解，被称作第di个元素的周期稀疏解xi_p，xi_p是维数为φ(di)的行矢量，含有φ(di)个元素；照此原则，计算出集合ds中所有元素的周期稀疏解；s32，计算标准能量和标准能量比矢量：s33，计算匹配成功次数nz和能量标准偏差γ；s34，确定参差脉冲序列的初始相位：根据初始相位匹配库每一行的匹配成功次数nz和能量标准偏差γ，选择匹配成功行；在初始相位匹配库中，选出与匹配成功次数nz的最大值对应的那些行，作为匹配成功的候选行；若仅有一个最大值，则这个候选行就是匹配成功行；若有多个最大值，再比较与候选行对应的能量标准偏差γ，选择具有最小能量标准偏差的候选行，作为匹配成功行；然后，依据匹配成功行在初始相位匹配库的位置，可以确定参差脉冲序列初始相位的个数和数值；s35，确定参差脉冲序列pri子周期：若初始相位的个数为一个，则序列是一个固定重频序列；若初始相位的个数为m＜mm，则序列是一个参差脉冲序列；此参差脉冲序列的pri子周期由各个初始相位θ之间的差值来确定，即：pri子周期分别为θ1……θm是各个初始相位。进一步的，所述s32中计算标准能量和标准能量比矢量包括以下步骤：s321，剔除xi_p内幅值较小的元素：对于集合ds内的每个元素di，定义其周期稀疏解xi_p内第j个元素幅值的绝对值为此元素的强度，记作|xi_p|j，j＝1,2,…,φ(di)；对于xi_p内所有强度值，判断其最大值为|xi_p|max和最小值为|xi_p|min，是否满足如下公式：|xi_p|min×mm＜|xi_p|max(4)若满足公式(4)，则说明xi_p内存在幅值较小的元素，需要在xi_p内找出强度值小于所有强度平均值的元素，使其强度值等于0；若不满足公式(4)，则说明xi_p内不存在幅值较小的元素；s322，求出xi_p内强度处于容差范围内的元素个数：对于xi_p内强度不等于0的第j个元素，给定其容差范围，即：range(j)＝[|xi_p|j×(1-δ)，|xi_p|j×(1+δ)]，其中，δ是小于1的正数；在xi_p内，找出强度处于容差范围range(j)内的元素，并记录具有不同强度值的元素的个数，记作ng(j)；按此原则，依次给定xi_p内所有强度不等于0的元素的容差范围，求出强度处于此容差范围内的元素个数；s323，确定标准能量：对于xi_p内所有强度不等于0的元素，找出具有最大ng(j)值的元素，将所述具有最大ng(j)值的元素强度记为|xi_p|n；以|xi_p|n为中心，给定其容差范围，即：range＝[|xi_p|n×(1-δ)，|xi_p|n×(1+δ)]；在xi_p内，找出强度处于此容差范围内的所有元素，并对具有不同强度值的所有元素求出其强度平均值，记作为元素di的标准能量；确定为每个元素di的标准能量；s324，确定标准能量比矢量：在xi_p内，找出所有不处于由下式给出的容差范围内的元素，使其幅值等于0；确定每个元素的标准能量比矢量即：其中，符号表示四舍五入取整；最终，计算确定出集合d内的所有元素的标准能量和标准能量比矢量。进一步的，所述s33中计算匹配成功次数nz和能量标准偏差γ的过程包括以下步骤：s331，从所述s22中的初始相位匹配库km中，取出一行数据构成一个匹配矢量rz；按照标准能量比矢量的长度，将匹配矢量rz相应地分成n个子矢量所述每个子矢量的长度与标准能量比矢量的长度相同，i＝1,2,…,n；逐一循环对子矢量与标准能量比矢量进行匹配，即：若标准能量比矢量中的每个元素值与子矢量中的每个元素值相同，则认为此次匹配成功，匹配成功次数nz增加1次，所述匹配成功次数nz的初始值为0；若标准能量比矢量中的每个元素值与子矢量中的每个元素值不相同，则认为此次匹配不成功，匹配成功次数nz保持不变；直至循环结束后，检查是否所有匹配都不成功；如果所有匹配都不成功，则匹配成功个数nz＝0，能量标准偏差γi＝0；否则，有一次匹配成功就进入下面计算：s332，计算总能量ed，即：其中，表示矢量中所有元素绝对值之和，ds表示所有匹配成功的标准能量比矢量的集合；s333，对于匹配不成功的标准能量比矢量重新计算标准能量比矢量将矢量与初始相位组合库的子矢量进行匹配；若两个矢量的每个元素值都相同，则认为此次匹配成功，匹配成功次数nz增加1次，并把这个矢量加入到集合ds中；同时按公式(6)重新计算ed；若两个矢量的每个元素值都不完全相同，则认为此次匹配不成功，计算能量标准偏差γi，即：其中，abs()表示绝对值操作；此次匹配结束后，重新计算下一个匹配不成功的标准能量比矢量，并且进行匹配；直至所有匹配不成功的标准能量比矢量都循环一遍为止；s334，计算总能量标准偏差γ，即：本发明的有益效果是，在已知参差脉冲序列及脉组参差序列的框架周期前提下，利用初始相位匹配库和多数匹配算法，提取出固定重频脉冲序列的初始相位，进而利用提取的初始相位与pri子周期的关系，实现了参差脉冲序列及脉组参差脉冲序列的pri子周期的提取。在存在大量的缺失和虚假脉冲的情况下，特别是对于具有不完整周期结构的脉冲序列，此提取方法仍然具有良好的提取性能。另外，此提取方法的适应性强，能够估计的脉冲类型较多，对加性结构和乘性结构的脉冲序列，诸如混叠抖动脉冲序列等，也能进行子周期的提取。附图说明为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1：本发明方法流程图。图2：不完整参差脉冲序列的形成过程。图3：在不同虚假率情况下，随着缺失率的变化，成功提取出参差脉冲序列三个pri子周期的概率曲线图。图4：在不同虚假率情况下，随着缺失率的变化，成功提取出参差脉冲序列三个pri子周期的概率曲线图。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明方法的流程如图1所示。在此实施例中，假设参差脉冲序列由三个固定重频脉冲序列叠加而成。三个固定重频脉冲序列的周期相同，即t1＝t2＝t3＝18；初始相位分别为θ1＝6，θ2＝10和θ3＝17。参差脉冲序列的形成过程如图2所示。从图2看出，参差脉冲序列的框架周期p＝18；三个pri子周期分别为和因各种原因，会造成参差脉冲序列中的个别脉冲缺失；因随机抖动误差，会使参差脉冲序列中的个别脉冲偏离原始位置，因而形成带有脉冲缺失和随机位置误差的不完整参差脉冲序列。下面介绍本发明如何从不完整参差脉冲序列中，提取参差脉冲序列pri子周期tirr的方法。可知，参差脉冲序列的pri子周期可以利用参差脉冲序列的框架周期和固定重频脉冲序列的初始相位来计算。假设参差脉冲序列的框架周期是已知的，或者通过估计算法已经得到。而固定重频脉冲序列的初始相位，则需要利用初始相位匹配库和匹配算法来获得。下面介绍初始相位匹配库的构建步骤和匹配算法的计算步骤。1、初始相位匹配库的构建1.1固定重频脉冲序列的初始相位匹配库构建假设固定重频脉冲序列的周期为pg，则初始相位匹配库有pg行，每行有φ(pg)个元素，这里φ(pg)被称为欧拉函数。由于周期为pg的ramanujan子空间内有φ(pg)个列，则初始相位匹配库的第一至第φ(pg)行中，将有一个具有正单位幅值的元素(简称正强度元素)，而其余元素的幅值为零。初始相位匹配库的第φ(pg)+1行至第pg行的元素由公式(1)确定：式中，dj是周期pg的第j个除数；si(i＝1,2,…,pg)表示初始相位匹配库中第i行的元素集，它是一个矢量，其中前φ(pg)行的元素集，即si(i＝1,2,…,φ(pg))，已由上面给出，它们是已知的。而后pg-φ(pg)行的元素集，即si(i＝φ(pg)+1,…,pg)，需按下述方法利用公式(1)计算得到：首先计算出周期pg的所有除数d，即pg除以d得到一个整数；根据除数d的大小，依照从大到小的顺序循环计算公式(1)，直至算出初始相位匹配库的第φ(pg)+1行至第pg行的所有元素为止。以周期pg＝18为例，具体说明计算过程。当pg＝18时，φ(pg)＝6。因此，初始相位匹配库的第1至第6行的元素集为s1＝[r1+，0，0，0，0，0]，s2＝[0，r2+，0，0，0，0]，s3＝[0，0，r3+，0，0，0]，s4＝[0，0，0，r4+，0，0]，s5＝[0，0，0，0，r5+，0]，s6＝[0，0，0，0，0，r6+]，其中上标“+”代表正强度元素，即元素幅值为正；而上标“-”代表负强度元素，即元素幅值为负。计算可得周期为18的所有除数为9，6，3，2。取第一个除数d1＝9，由于pg/d1＝18/9＝2，因此公式(1)左侧应有两项相加，即这里i可取1，2，3，4，5，6。例如当i取1时，s1+s10＝0。于是，s10＝-s1＝-[r1+，0，0，0，0，0]。因为r1+是正强度元素，所以初始相位匹配库的第10行元素集含有一个负强度元素，即s10＝[r1-，0，0，0，0，0]。依此计算方法，当i可取2，3，4，5，6时，可得s11＝[0，r2-，0，0，0，0]，s12＝[0，0，r3-，0，0，0]，s13＝[0，0，0，r4-，0，0]，s14＝[0，0，0，0，r5-，0]，s15＝[0，0，0，0，0，r6-]。取第二个除数d2＝6，公式(1)左侧应有三项相加，即这里i可取1，2，3，4，5，6。例如当i取1时，s1+s7+s13＝0。于是，s7＝-s1-s13＝-[r1+，0，0，0，0，0]-[0，0，0，r4-，0，0]。因而初始相位匹配库的第7行元素集含有一个正强度元素和一个负强度元素，即s7＝[r1-，0，0，r4+，0，0]。依此计算方法，当i可取2，3，4，5，6时，可得s8＝[0，r2-，0，0，r5+，0]，s9＝[0，0，r3-，0，0，r6+]，s16＝[r1+，0，0，r4-，0，0]，s17＝[0，r2+，0，0，r5-，0]，s18＝[0，0，r3+，0，0，r6-]。至此，初始相位匹配库的第1行至第18行所有元素全部得到，初始相位匹配库构建完成，如表1所示。表1周期为18的固定重频脉冲序列初始相位匹配库1.2参差脉冲序列的初始相位匹配库构建参差脉冲序列是由多个固定重频脉冲序列叠加而成的。参差脉冲序列的特点是组成的多个固定重频脉冲序列的周期相等，但初始相位不同。因此，构建初始相位匹配库需要分两步进行。a、含有单个初始相位的参差脉冲序列假设参差脉冲序列的框架周期为p。首先计算出周期为p的所有除数d，即p除以d得到一个整数；然后，将周期为p以及所有除数d的固定重频脉冲序列初始相位匹配库，根据从大到小的顺序整合到一个行数为p的参差脉冲序列初始相位匹配库中。由于除数d小于p，且p除以d是一个整数，因而周期为d的固定重频脉冲序列初始相位匹配库将在参差脉冲序列初始相位匹配库中，循环出现p除以d次。例如参差脉冲序列的框架周期p＝18，则周期为p的所有除数d是9、6、3、2。通过整合，共有18行的参差脉冲序列初始相位匹配库如表2所示。表2框架周期为18、含有单个初始相位的参差脉冲序列初始相位匹配库行数p＝18d＝9d＝6d＝3d＝21r1+r1+r1+r1+r1+2r2+r2+r2+r2+r1-3r3+r3+r1-r2+r1-r2-r1+4r4+r4+r1-r1+r1-5r5+r5+r2-r2+r1+6r6+r6+r1+r2-r1-r2-r1-7r1-r4+r1-r4-r1+r1+r1+8r2-r5+r2-r5-r2+r2+r1-9r3-r6+r3-r6-r1-r2+r1-r2-r1+10r1-r1+r1-r1+r1-11r2-r2+r2-r2+r1+12r3-r3+r1+r2-r1-r2-r1-13r4-r4+r1+r1+r1+14r5-r5+r2+r2+r1-15r6-r6+r1-r2+r1-r2-r1+16r1+r4-r1-r4-r1-r1+r1-17r2+r5-r2-r5-r2-r2+r1+18r3+r6-r3-r6-r1+r2-r1-r2-r1-注：在表2中，为了简单明了，省去所有元素集的零元素。另外，下脚标的数字表示此非零元素处于元素集的位置。b、含有多个初始相位的参差脉冲序列当已知含有单个初始相位参差脉冲序列的初始相位匹配库(简称为k1)时，通过对此初始相位匹配库的各行进行组合，可构建含有m个不同初始相位参差脉冲序列的初始相位匹配库(简称为km)。构建方法如下：假设参差脉冲序列的框架周期为p，此参差脉冲序列由m个不相同初始相位的固定重频脉冲序列混叠组成，因此，此参差脉冲序列含有m个不相同的初始相位。从参差脉冲序列框架周期为p的初始相位匹配库k1中，任意提取出m个行，并且将其相应的元素相加作为初始相位匹配库km的一行。由于在p行中，任取m行进行组合，共可得到nu个行，因此，初始相位匹配库km共有nu个行。这里，nu由公式(2)求得：例如初始相位个数m＝2，周期p＝18时，初始相位匹配库k2共有153行。其初始相位匹配库如表3所示。表3框架周期为18、含有两个不同初始相位的参差脉冲序列初始相位匹配库注意：在表2中，由于行数较多，只列举其中一部分行。标注“0”表示此元素集的所有元素全为零。对于每一行，标注出了它对应的两个初始相位值θ。例如，第18行对应的一个初始相位是2，另一个初始相位是4。2、匹配算法假设一个参差脉冲序列的框架周期为p，具有小于mm个初始相位，这里mm为正整数，且1＜mm≤p。计算由p的除数(包括p)组成的集合ds，即：ds＝{d1,d2,…,dn}，这里：p/di＝整数；di≠1；i＝1,2,…,n。s1，计算周期稀疏解从集合ds取出一个元素di，构成一个以此元素值为周期的ramanujan结构字典a，求解以下优化问题：min||dx||2s.t.y＝φx；φ＝ad(3)其中y代表测量的脉冲序列，x代表周期稀疏解，d表示权重矩阵，起到了均衡各个周期列权重的作用。此问题的最优解，被称作第di个元素的周期稀疏解xi_p，这里xi_p是维数为φ(di)的行矢量，即它含有φ(di)个元素。照此原则，计算出集合ds中所有元素的周期稀疏解。s2，计算标准能量和标准能量比矢量对于集合ds内的每个元素di，定义其周期稀疏解xi_p内第j个元素幅值的绝对值为此元素的强度，记作|xi_p|j，j＝1,2,…,φ(di)。s2-1，剔除xi_p内幅值较小的元素对于xi_p内所有强度值，其最大值为|xi_p|max和最小值为|xi_p|min，若下式成立|xi_p|min×mm＜|xi_p|max(4)则说明xi_p内存在幅值较小的元素。因此，需要在xi_p内找出强度值小于所有强度平均值的元素，使其强度值等于0。若(3)式不成立，则说明xi_p内不存在幅值较小的元素。s2-2，求出xi_p内强度处于容差范围内的元素个数对于xi_p内强度不等于0的第j个元素，给定其容差范围，即range(j)＝[|xi_p|j×(1-δ)，|xi_p|j×(1+δ)]，其中δ是小于1的正数并事先给定的。在xi_p内，找出强度处于容差范围range(j)内的元素。在这些元素中，记录具有不同强度值的元素的个数，记作ng(j)。按此原则，依次给定xi_p内所有强度不等于0的元素的容差范围，求出强度处于此容差范围内的元素个数。s2-3，确定标准能量对于xi_p内所有强度不等于0的元素，找出具有最大ng(j)值的那个元素，将这个元素的强度记为|xi_p|n。以|xi_p|n为中心，给定其容差范围，即range＝[|xi_p|n×(1-δ)，|xi_p|n×(1+δ)]。在xi_p内，找出强度处于此容差范围内的所有元素。在这些元素中，对具有不同强度值的所有元素求出其强度平均值，记作确定为每个元素di的标准能量。s2-4，确定标准能量比矢量在xi_p内，首先找出所有不处于由下式给出的容差范围内的元素，使其幅值等于0。然后，确定每个元素的标准能量比矢量即：这里符号表示四舍五入取整。按以上4步，计算集合ds内的所有元素的标准能量和标准能量比矢量。s3，计算匹配成功次数nz和能量标准偏差γ从事先计算好的初始相位匹配库(见1.2中的表3)中，按从小到大的顺序取出每一行与标准能量比矢量匹配，匹配结果由匹配成功次数nz和能量标准偏差γ体现。具体匹配步骤如下：s3-1，从初始相位匹配库中，取出一行数据构成一个匹配矢量rz。按照标准能量比矢量的长度，将匹配矢量rz相应地分成n个子矢量这里每个子矢量的长度与标准能量比矢量的长度相同(i＝1,2,…,n)。然后，逐一地循环对子矢量与标准能量比矢量进行匹配，即：如果标准能量比矢量中的每个元素值与子矢量中的每个元素值相同，则认为此次匹配成功，匹配成功次数nz增加1次(匹配成功次数nz的初始值为0)。否则认为此次匹配不成功，匹配成功个数nz保持不变。直至循环结束后，检查是否所有匹配都不成功，如果所有匹配都不成功，则匹配成功个数nz＝0，能量标准偏差γi＝0。否则，进入下面各步计算。s3-2，计算总能量ed，即式中表示矢量中所有元素绝对值之和，ds表示所有匹配成功的标准能量比矢量集合。s3-3，对于匹配不成功的标准能量比矢量重新计算标准能量比矢量然后，将矢量与初始相位组合库的子矢量进行匹配。如果两个矢量的每个元素值都相同，则认为此次匹配成功，匹配成功次数nz增加1次，并把这个矢量加入到集合ds中。同时按式(6)重新计算ed。如果两个矢量的每个元素值都不完全相同，则认为此次匹配不成功，需计算能量标准偏差γi，即式中abs()表示绝对值操作。此次匹配结束后，转入s3-3，重新计算下一个匹配不成功的标准能量比矢量，并且进行匹配。直至所有匹配不成功的标准能量比矢量都循环一遍为止。s3-4，计算总能量标准偏差γ，即s4，确定参差脉冲序列的初始相位首先，根据初始相位匹配库每一行的匹配成功次数nz和能量标准偏差γ，选择匹配成功行。在初始相位匹配库中，选出与匹配成功次数nz的最大值对应的那些行，作为匹配成功的候选行；若仅有一个最大值，则这个候选行就是匹配成功行；若有多个最大值，再比较与候选行对应的总能量标准偏差γ，选择具有最小能量标准偏差的候选行，作为匹配成功行。然后，依据匹配成功行在初始相位匹配库的位置，可以确定参差脉冲序列初始相位的个数和数值。s5，确定参差脉冲序列pri子周期若初始相位的个数为一个，则此序列是一个固定重频序列；若初始相位的个数为m＜mm，则此序列是一个参差脉冲序列。此参差脉冲序列的pri子周期由各个初始相位θ之间的差值来确定，即：pri子周期分别为这里p为参差脉冲序列的框架周期。3、具体实施效果在不同缺失率和虚假率条件下，针对此实施例中列举的、具有三个pri子周期的参差脉冲序列，利用上述匹配方法，成功提取三个pri子周期的概率变化曲线如图3；而利用传统的谱估计算法，成功提取三个pri子周期的概率变化曲线如图4。在图中，缺失率被定义为相同长度脉冲序列中缺失的脉冲个数与完整无缺失的脉冲个数的比值，虚假率(即：falserate)被定义为相同长度脉冲序列中杂散脉冲的个数与完整无缺失脉冲序列个数的比值。缺失表现方式包括了连续缺失和离散缺失，首脉冲和末端脉冲的缺失，产生虚假脉冲方式为在脉冲序列中人为地均匀随机添加虚假脉冲。脉冲序列长度为300。图中曲线是100次montecarlo仿真实验结果的平均值。从图3中看出，随着缺失率不断升高，虚假率越小，成功提取概率受到的影响越小；在具有相同缺失率的情况下，虚假率变化越大，成功提取概率也随之严重下降。这说明当缺失率和虚假率为相同值时，虚假率比缺失率对成功提取概率具有更强的影响。比较图3和图4可知，相对于传统的谱估计方法，本发明提出方法都具有更好的提取性能。特别在虚假率达到0.3以上，且缺失率高达0.6时，本发明提出方法依然能够具有0.4的成功提取概率，而传统的谱估计方法的成功提取概率已接近于0。这说明在高缺失和高虚假的情况下，本发明提出方法仍具有更高的估计成功提取概率。以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。当前第1页12