CN101193087A - 一种信号捕获方法 - Google Patents

Abstract

本发明公开了一种信号捕获方法，用于在信号接收端的信号同步过程中捕获信号同步位置，在捕获过程的每次搜索中，从输入信号序列中抽取若干信号作为信号子序列进行自相关运算，寻找到本次搜索的自相关输出值的幅度值连续大于预先设定的幅度值门限T_AMP的最宽位置序列，并记录宽度值和终点位置，若所述宽度值大于预先设定的宽度门限T_WID，则表示本次搜索成功；若符合捕获条件捕获成功，根据所述宽度值、终点位置及信号子序列与信号序列的间隔值得到信号序列的初始同步位置。本发明采用对部分输入信号进行自相关运算，降低了捕获功能的复杂度，同时，可灵活设置捕获参数避免持续搜索或搜索时间过长等浪费终端电能的弊端，具有低复杂度、低功耗等特点。

Description

一种信号捕获方法

技术领域

本发明涉及通信技术领域的同步接收技术，尤其涉及同步接收过程中的一种信号捕获方法。

背景技术

所谓信号捕获，就是指接收系统信号的初始同步过程。接收系统只有在获得初始同步后，再进行精细同步，并且在接收信号的过程中不断进行同步跟踪，才可以确保接收系统获得真正的信号同步。因此获得信号的初始同步是接收机同步系统的第一个处理步骤，也是十分重要的处理步骤。信号捕获的目的是获得粗略的信号同步位置，而非确切的信号同步位置。并且由于无线信道的多样性，某种信号捕获方法很难保证在所有信道环境下都能够获得确切的信号同步位置；或者说，针对某些信道情况，该算法能够获得确切的信号同步位置，但对于其它的信道情况，该算法甚至存在失效的可能性。因此，信号捕获算法的一个最重要特征是应具有通用性，即可以在任何信道情况下，获得接收系统的粗略同步位置。进一步地，获得粗略同步后，再利用精细同步，就可以获得较确切的信号同步位置了。

欧洲的DVB(Digital Video Broadcasting，数字视频广播)系统采用了多载波OFDM(Orthogonal Frequency Division Multiplexing，正交频分复用)信道传输技术。DVB系统的信号捕获方法利用了OFDM符号的循环前缀(CP，Cyclic Prefix)部分。由于CP部分长度较短，如在DVB系统中2K模式下的最短CP长度仅为7微秒，因此无法保证信号捕获算法的鲁棒性；并且，由于CP部分无法准确标识每个信号帧的起始位置，因此即使可以获得某个OFDM符号的起始位置，接收系统也无法快速和确切地寻找到每个信号帧的起始位置。

为了避免DVB系统面临的同步困境，中国移动多媒体广播CMMB(China Mobile Multimedia Broadcasting)系统提出了一种新的帧结构，在每个时隙的起始处设置了一个信标，从而确保了接收机可以不利用OFDM符号的CP部分进行同步，仅仅利用信标部分就可以获得信号的同步位置。

如图1所示，给出了CMMB系统物理层的帧结构示意图。在CMMB系统物理层的帧结构中，将1秒等分为40个时隙，每个时隙由1个信标和53个OFDM符号组成。

如图2所示，给出了CMMB系统物理层的帧结构的每一时隙中信标的组成结构示意图。在每一时隙中，每个信标由一个TxID(发射机标识序列)和两个同样的同步信号(SYNC，Synchronous signal)组成。每个同步信号的时间长度为204.8微秒，远远大于OFDM符号的CP长度(51.2微秒)。

对于具有类似CMMB系统的两个同步信号的信号，现有信号捕获方法都是利用信标中的该两个同步信号即SYNC信号进行自相关操作来进行捕获。由于同步信号具有较大长度，因而必然地增大了信号捕获的实现复杂度。此外，由于网络覆盖尚未完善，在某些地段，可能存在移动终端无法接收到信号的情况，在这种情况下，目前的信号捕获方法由于持续搜索，可能将导致移动终端消耗大量的电能，不具有节电功能。

发明内容

本发明要解决的技术问题是提供一种信号捕获方法，用于对对于具有类似CMMB系统的两个同步信号的信号进行捕获，解决现有捕获方法中复杂度高、捕获过程耗电量高等问题，以低复杂度、低功耗、捕获方式灵活地实现信号捕获，确定粗同步位置。

本发明采用的技术方案是：

一种信号捕获方法，用于在信号接收端的信号同步过程中捕获信号同步位置，包括如下步骤：

A、在一设定时间段L_TIM内进行多次信号搜索，在每次搜索中从输入信号序列中抽取若干信号作为信号子序列进行自相关运算，获取本次信号搜索的自相关输出值；

B、在每次搜索中，寻找到本次搜索的自相关输出值的幅度值连续大于预先设定的幅度值门限T_AMP的最宽位置序列，并记录宽度值和终点位置，若所述宽度值大于预先设定的宽度门限T_WID，则表示本次搜索成功；

C、若在所述时间段L_TIM内连续搜索成功次数超过设定的次数N_SUC，则表示捕获成功，根据所述宽度值、终点位置及信号子序列与信号序列的间隔值得到信号序列的初始同步位置。

进一步地，所述步骤A中：

从输入信号序列中抽取信号子序列的步骤，是利用抽值因子F_D来进行的，即：从信号序列的初始信号开始每个抽值因子F_D个间隔抽取一个信号来组成所述信号子序列。

进一步地，所述步骤A中：

所述时间段L_TIM在信号捕获过程中是周期性时间段，即每隔时间段L_TIM进行一次信号捕获过程，在所述时间段L_TIM内设定有该时间段内最多允许进行的连续搜索次数N_SCH、每次搜索的时间长度T_ONE、判断捕获成功的连续成功搜索次数N_SUC。

进一步地，所述步骤A中：

所述对信号子序列进行自相关运算，是通过输入信号的共轭转置值与延迟后的输入信号相乘来得到，其中：

第一次搜索相应的第一个自相关输出值acorr_out(1)，通过下式获得：

acorr_out(1)＝in_sig(1:F_D:win_len)′*in_sig(delay+1:F_D:delay+win_len)；

其中，所述F_D表示抽值因子，所述in_sig(1:F_D:win_len)表示输入信号子序列；所述win_len表示窗口长度，用于表示为获得一个自相关输出值，对于参与运算的前后两段同步序列中的每段序列从序列的起点到序列的终点的距离，所述delay表示延时，所述右上标志“′”表示共轭转置操作；

对于第一个自相关输出值之后的自相关输出值acorr_out(k)，采用如下的滑动窗方法获得：

acorr_out(k)＝acorr_out(k-1)-

conj(in_sig((k-2)*F_D+1))*in_sig((k-2)*F_D+1+delay)+

conj(in_sig((k-2)*F_D+1+win_len))*

in_sig((k-2)*F_D+1+win_len+delay)；

其中，acorr_out(k)表示第k个自相关输出值，acorr_out(k-1)表示第k-1个自相关输出值，conj()函数表示取共轭函数，由上式可知，对于第K个自相关输出值，由前一次的自相关输出值以及k-2次输入信号、抽值因子利用上述的滑动窗方法得到。

进一步地，所述步骤B中通过如下步骤获取所述自相关输出值的幅度值：

分别获取自相关输出值的实部与虚部的绝对值，并从所述绝对值中选择出最大值以及最小值，再将该最大值加上最小值的二分之一的和，作为自相关输出值的幅度值acorr_amp。

进一步地，若步骤A中利用抽值因子获得信号子序列，相应地所述步骤C中，根据抽值因子F_D、终点位置end_idx及宽度值width确定信号序列的初始同步位置sync_idx为：

sync_idx＝F_D*(end_idx-floor(width/2))；

其中，所述floor()表示向下取整函数。

进一步地，所述方法进一步包括：

若捕获不成功，则依照继续进行捕获过程。

进一步地，所述方法进一步包括：

根据历史搜索结果调整所述步骤A、B、C中的搜索参数，所述搜索参数包括一次搜索的时间长度T_ONE、进行捕获的时间段的长度L_TIM、所述时间段L_TIM内的最多连续搜索次数N_SCH、以及连续成功搜索次数N_SUC。

进一步地，所述方法进一步包括：

若在所述时间段L_TIM内，最多连续搜索次数N_SCH与已进行的搜索次数之差小于或等于连续成功搜索次数N_SUC，则结束搜索。

进一步地，在单频干扰的信道环境下，所述步骤B中，通过直接搜索自相关输出幅度值的最大值来确定最宽的位置序列。

进一步地，所述方法进一步包括：

D、根据初始位置sync_idx、终点位置end_idx和宽度值width计算得到小数频偏。其中，计算小数频偏的步骤具体包括：

D1、判断初始位置sync_idx与本次搜索的起始位置之差是否大于所设定的每次搜索的时间长度的一半，

若是，则计算频率偏差索引foe_idx位置为：

foe_idx＝end_idx-floor(width/2)+62500；

若不是，则计算频率偏差索引foe_idx位置为：

foe_idx＝end_idx-floor(width/2)+125000；

D2、在确定foe_idx位置后，进一步估计小数频率偏差为：

$\mathrm{\Δf}=\frac{\mathrm{angle}(\mathrm{acorr}\_\mathrm{out}(\mathrm{foe}\_\mathrm{idx}))}{2\mathrm{\π}}*\frac{{10}^7}{2048}$

其中，angle()表示取角度函数函数，floor()表示向下取整函数。

本发明提出的CMMB系统信号捕获方法与目前采用的其它CMMB系统信号捕获方法同样都利用了信标中的两个SYNC信号，同样都执行了自相关操作，但针对现有捕获方法复杂度高以及不节电等问题，本发明采用对部分输入信号进行自相关运算，降低了捕获功能的复杂度，同时，在搜索过程中，通过设定搜索次数以及搜索成功条件，以及进一步可动态调整捕获参数，可以避免持续搜索或搜索时间过长等浪费终端电能的弊端，具有低复杂度、低功耗的显著优点。

附图说明

图1是背景技术中CMMB系统物理层的帧结构示意图；

图2是背景技术中CMMB系统物理层的帧结构的每一时隙中信标的组成结构示意图；

图3是本发明实施例中实现信号捕获的流程图；

图4是现有技术中对信号进行自相关运算后自相关输出值的波形图。

具体实施方式

下面将结合附图及实施例对本发明的技术方案进行更详细的说明。

本发明提出了一种低复杂度、低功耗的信号捕获方法，用于捕获信号，即实现信号的粗同步，可以应用于具有与CMMB系统类似帧结构的系统中，即在信标部分具有两个同样的同步信号。

如图3所示，给出了本发明信号捕获方法的一个实施例，信号捕获流程包括如下步骤：

步骤301、在一个设定的时间段L_TIM内进行连续最多N_SCH次信号搜索，每次搜索时间长度为T_ONE，在每次搜索中，从接收的信号序列中抽取若干信号作为信号子序列，对抽取出的信号子序列进行自相关运算，得到每次信号搜索的自相关输出值；

步骤302、每次搜索时，找到本次自相关输出幅度值连续大于幅度值门限T_AMP的最宽的位置序列，并记下宽度值和终点位置，若所述宽度值大于所设定的宽度门限T_WID，则表示本次搜索成功；

步骤303、若在所述时间段L_TIM内连续搜索成功次数超过设定的次数N_SUC，则表示捕获成功，则根据所述宽度值、终点位置及信号子序列与信号序列的间隔值得到信号序列的初始同步位置。

在上述捕获方法中，采用抽值因子F_D来从信号序列中获得一个信号子序列，例如in_sig()表示输入信号序列，在F_D＝4时，信号子序列in_sig(1:F_D:win_len)表示输入信号序列in_sig(1)、in_sig(5)、in_sig(9)......；

在判断捕获成功后，则根据抽值因子、终点位置及宽度值计算确定信号序列的初始同步位置sync_idx为：

sync_idx＝F_D*(end_idx-floor(width/2))；

先对二分之一宽度值向下取整，再利用终点位置end_idx减去取整得到的整数值，即可得到自相关输出值幅度波形中一个波峰的峰值位置，由于本发明采用抽值自相关运算，在计算初始同步位置时，需要进行校正，即还需要对峰值位置再乘以一个抽值因子F_D后，即可最终得到初始同步位置，即确定信号捕获所要的粗同步位置。

若捕获不成功，则依照捕获规则继续捕获过程。

上述步骤中所用参数需要预先设定，涉及如下参数：

设定“抽值因子”F_D，由于根据抽值因子抽取部分信号进行自相关运算，减少了运算量，可以实现低复杂度处理；

设定“幅度门限”值T_AMP用于根据该门限值判断幅度是否满足信号要求，即大于该门限可认为是信号幅度；

设定“宽度门限”值T_WID用于在幅度波形中判断波峰宽度大于该门限时认为是一个信号峰；

一次搜索的时间长度T_ONE，用于表示完成一次搜索设定的时间值；

一个时间段的长度L_TIM，表示进行捕获搜索的时间段长度；

一段时间内的最多(连续)搜索次数N_SCH、表示例如L_TIM内连续最多搜索的次数；

以及判断搜索成功时的(连续)成功搜索次数N_SUC，若连续成功搜索次数超过N_SUC则认为本次捕获成功。

若上述参数设计合理，则还可以达到节电的效果(低功耗)。设定上述参数后，若捕获不成功，则依照捕获规则继续捕获。必要时，可以更新上述参数来调整捕获规则的细节。由于本发明本身减小了计算量即可节省运算耗能，同时，本发明在捕获过程中还可根据历史捕获结果来动态调整上述参数，例如，若多次捕获不成功时，可以增大时间段长度，减少连续搜索次数。

同时，根据具体参数，还可以进行计算上的简化，如N_SCH为4和N_SUC＝3时，若前2次搜索不成功，则不必进行随后2次的搜索，因为即使搜索成功也不会满足捕获成功的条件，这就减少了搜索次数，可以低功耗；若不考虑低功耗，则可以继续搜索。

在获得第一个同步(SYNC)信号的起始位置后，可以随周期(时隙)调整。

进一步地，为了使得同步位置更准确还可进行小数频偏计算进行频率校正，那么在捕获成功后，继续做一个时隙的自相关，并根据周期调整后的起始位置、抽值因子和自相关输出值计算小数频偏，然后终止捕获过程。接收机的整体状态机转移到下一个状态。

以下是本发明实施例中具体算法的详细描述。

步骤1、在信号捕获过程中，设定抽值因子F_D＝4，根据抽值因子F_D从输入的信号序列中选择等间隔数目的信号进行自相关运算。所述抽值因子F_D表示在输入信号序列中每隔F_D个数目选择一个信号，对选出的信号子序列进行自相关运算。

A1、对于第一个自相关输出值，通过输入信号的共轭转置值与延迟后的输入信号进行相乘来得到。在具体实现时，第一个自相关输出值acorr_out(1)如下所示：

acorr_out(1)＝in_sig(1:F_D:win_len)′*in_sig(delay+1:F_D:delay+win_len)；

其中，in_sig()表示输入信号序列，为列向量，在F_D＝4时，in_sig(1:F_D:win_len)表示输入信号序列in_sig(1)、in_sig(5)、in_sig(9)......；

win_len表示窗口长度，此处设win_len＝2048，该窗口长度用于表示为获得一个自相关输出值，对于参与运算的前后两段同步序列中的每段序列，从序列的起点到序列的终点的距离，以出现的数据样本而非实际参与运算的数据样本的数目表示，因为可能有的数据样本不参与运算。

delay表示延迟，对于信标中的两个连续同步信号，必然存在一定的延时，此处设delay＝2048；

右上标′，表示共轭转置操作。

A2、对于第一个自相关输出值之后的自相关输出值，可以采用如下的滑动窗方法获得：

acorr_out(k)＝acorr_out(k-1)-

conj(in_sig((k-2)*F_D+1))*in_sig((k-2)*F_D+1+delay)+

conj(in_sig((k-2)*F_D+1+win_len))*

in_sig((k-2)*F_D+1+win_len+delay)；

其中，acorr_out(k)表示第k个自相关输出值，acorr_out(k-1)表示第k-1个自相关输出值，conj()函数表示取共轭函数，由上式可知，对于第K个自相关输出值(对应第k次搜索)，可以由前一次的自相关输出值以及k-2次输入信号、抽值因子利用上述的滑动窗方法得到。

也就是说我们得到第一个自相关输出值后，可以依次得到后续的若干个自相关输出值。

例如，对于上式，k＝2时，acorr_out(2)＝acorr_out(1)-

conj(in_sig(1))*in_sig(1+delay)+

conj(in_sig(1+win_len))*in_sig(1+win_len+delay)；

例如，对于上式，k＝3时，acorr_out(3)＝acorr_out(2)-

conj(in_sig(F_D+1))*in_sig(F_D+1+delay)+

conj(in_sig(F_D+1+win_len))*in_sig(F_D+1+win_len+delay)；

步骤2、获取自相关输出值的幅度值，并根据捕获判断规则以及所述自相关输出值的幅度值，判断是否捕获成功。

B1、对于搜索获得的各个自相关输出值acorr_out取绝对值则可得到其幅度值acorr_amp，本实施例中可以采用如下近似计算方法来确定acorr_amp：

acorr_amp＝max(abs(real(acorr_out))，abs(imag(acorr_out)))+0.5*min(abs(real(acorr_out))，abs(imag(acorr_out)))；

即：分别获取自相关输出值的实部与虚部的绝对值，并从中选择出最大值以及最小值，再将该最大值加上最小值的二分之一的和，作为自相关输出值的绝对值(即幅度值)acorr_amp。

B2、设定用于捕获判断的各种参数。

本发明实施例中进行捕获判断时，需要在一个时间段内进行若干次连续搜索，对与每一次搜索具有一个搜索时间长度，判断连续搜索成功次数为若干次后即可认为搜索成功。针对上述判断流程，可设定如下参数：

一个时间段长度L_TIM，用于表示需要进行搜索的时间段长度，此处可以设为15个时隙(time slots)；

一次搜索的时间长度T_ONE，用于表示每次搜索所需要花费的时间，也即该次搜索持续的时间，此处设为1.5个时隙；

一个时间段内的最多的连续搜索次数N_SCH，即表示在设定的时间段L_TIM内至多可进行搜索的次数，此处设为4次；

连续成功搜索次数N_SUC，用于判断是否捕获成功，若在设定时间段L_TIM内连续搜索成功达到N_SUC次，则可判定本次搜索成功，此处设N_SUC为2次。

同时，为了判断信号，还需要设定自相关输出波形的幅度门限值T_AMP以及输出波形的宽度门限值T_WID，该两个门限值的设定可以根据系统性能要求及信号的功率、延迟等特性进行确定，此处设置幅度门限值T_AMP为127，宽度门限值T_WID为63。

B3、在每次搜索中，寻找到自相关输出幅度值连续大于T_AMP的最宽的位置序列，并记下宽度和终点位置。

如图4所示，显示了现有技术中对信号进行自相关运算后，得到自相关输出值的波形图。进行搜索时，在每次搜索中根据计算得到的本次自相关输出的幅度值，判断是否大于所设定的幅度值门限T_AMP，若大于则表明此处为信号对应区域，若连续大于T_AMP则标明此处为信号对应的波峰，此时，找到最宽的位置序列，记下波形的宽度值(即大于T_AMP的终点值与起点值之差)，以及终点位置。

若出现多个最宽序列，则取最后一个最宽序列。若宽度大于波形宽度值门限T_WID，则表明本次搜索成功；否则表示本次搜索不成功，依据相应的参数进行下一步操作。若连续搜索成功2次，则表明捕获成功。

由于自相关输出值是基于前次搜索的自相关输出值进行的，因此，在搜索过程中可以动态调整上述设定的参数，例如，若多次判断搜索失败，一方面可以增加搜索时间段的长度，或者可以减少搜索时间段内连续搜索的次数来达到节省电能的目的。

需要指出的是，在某些信道环境下，如存在单频干扰的信道环境下，也可以直接搜索自相关输出幅度值的最大值。事实上，在这些信道环境下，最大值往往处于最宽序列里。

步骤3、在捕获成功后，根据最后获得的自相关的宽度值width和终点位置end_idx，以及抽值因子F_D，可以计算出接收到的信号中的第一个同步(SYNC)信号的起始位置sync_idx，公式如下：

sync_idx＝F_D*(end_idx-floor(width/2))；

上式中floor()表示向下取整，终点位置end_idx减去波形宽度的一半即可得到波峰的位置，再与抽值因子F_D相乘，即可得到第一个同步信号的起始位置sync_idx。

由于本发明实施例中，利用抽值因子F_D，对信号序列并不是全部进行自相关运算，而是从中间隔F_D个信号进行自相关运算，因此，需要确定原输入信号序列中的初始位置，利用F_D进行校正后即可得到相对准确的粗同步位置。

步骤4、在信号同步过程中，为了进一步优化同步结果，往往还需要进行小数频偏处理，对搜索确定的粗同步位置进行校正。

D1，判断sync_idx与本次搜索的起始位置之差是否大于所设定的每次搜索的长度的一半，此处即判断是否大于0.75个时隙长度：

若是，则根据下述公式计算频率偏差索引foe_idx位置，继续进行自相关操作至该位置，然后终止自相关操作，然后计算小数频偏，然后捕获过程终止；其中，foe_idx＝end_idx-floor(width/2)+62500；

若不是，则根据下述公式计算foe_idx，继续进行自相关操作至该位置，然后终止自相关操作，然后计算小数频偏，然后捕获过程终止，其中，

foe_idx＝end_idx-floor(width/2)+125000；

在上述计算foe_idx的步骤中，所述参数62500、125000仅为示例性的，根据具体情况可以选择不同数值，本示例中将一个时隙分为250000个样本点，62500/250000＝0.25＝1-0.75个时隙(time slot)，同理，125000＝0.5个时隙(time slot)。具体数值的确定也可以依照仿真验证来得到。

D2、在确定foe_idx位置后，进一步估计小数频率偏差。小数频偏的计算公式如下(单位Hz)：

$\mathrm{\Δf}=\frac{\mathrm{angle}(\mathrm{acorr}\_\mathrm{out}(\mathrm{foe}\_\mathrm{idx}))}{2\mathrm{\π}}*\frac{10^7}{2048}$

尽管本发明提出的CMMB系统信号捕获方法与目前采用的其它CMMB系统信号捕获方法同样都利用了信标中的两个SYNC信号，同样都执行了自相关操作，但目前采用的其它CMMB信号捕获方法表现出了较大的实现复杂度。

针对现有捕获方法复杂度高以及不节电等问题，本发明的信号捕获方法采用对部分输入信号进行自相关运算，降低了捕获功能的复杂度，同时，在搜索过程中，通过设定搜索次数以及搜索成功条件，可以避免持续搜索或搜索时间过长等浪费终端电能的弊端，具有低复杂度、低功耗的显著优点。

当然，本发明还可有其他多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (12)

1.一种信号捕获方法，用于在信号接收端的信号同步过程中捕获信号同步位置，其特征在于，包括如下步骤：

A、在一设定时间段L_TIM内进行多次信号搜索，在每次搜索中从输入信号序列中抽取若干信号作为信号子序列进行自相关运算，获取本次信号搜索的自相关输出值；

B、在每次搜索中，寻找到本次搜索的自相关输出值的幅度值连续大于预先设定的幅度值门限T_AMP的最宽位置序列，并记录宽度值和终点位置，若所述宽度值大于预先设定的宽度门限T_WID，则表示本次搜索成功；

C、若在所述时间段L_TIM内连续搜索成功次数超过设定的次数N_SUC，则表示捕获成功，根据所述宽度值、终点位置及信号子序列与信号序列的间隔值得到信号序列的初始同步位置。

2.如权利要求1所述的信号捕获方法，其特征在于，所述步骤A中：

从输入信号序列中抽取信号子序列的步骤，是利用抽值因子F_D来进行的，即：从信号序列的初始信号开始每个抽值因子F_D个间隔抽取一个信号来组成所述信号子序列。

3.如权利要求1所述的信号捕获方法，其特征在于，所述步骤A中：

所述时间段L_TIM在信号捕获过程中是周期性时间段，即每隔时间段L_TIM进行一次信号捕获过程，在所述时间段L_TIM内设定有该时间段内最多允许进行的连续搜索次数N_SCH、每次搜索的时间长度T_ONE、判断捕获成功的连续成功搜索次数N_SUC。

4.如权利要求2所述的信号捕获方法，其特征在于，所述步骤A中：

所述对信号子序列进行自相关运算，是通过输入信号的共轭转置值与延迟后的输入信号相乘来得到，其中：

第一次搜索相应的第一个自相关输出值acorr_out(1)，通过下式获得：

acorr_out(1)＝in_sig(1:F_D:win_len)′*in_sig(delay+1:F_D:delay+win_len)；

其中，所述F_D表示抽值因子，所述in_sig(1:F_D:win_len)表示输入信号子序列；所述win_len表示窗口长度，用于表示为获得一个自相关输出值，对于参与运算的前后两段同步序列中的每段序列从序列的起点到序列的终点的距离，所述delay表示延时，所述右上标志“′”表示共轭转置操作；

对于第一个自相关输出值之后的自相关输出值acorr_out(k)，采用如下的滑动窗方法获得：

acorr_out(k)＝acorr_out(k-1)-

conj(in_sig((k-2)*F_D+1))*in_sig((k-2)*F_D+1+delay)+

conj(in_sig((k-2)*F_D+1+win_len))*

in_sig((k-2)*F_D+1+win_len+delay)；

其中，acorr_out(k)表示第k个自相关输出值，acorr_out(k-1)表示第k-1个自相关输出值，conj()函数表示取共轭函数，由上式可知，对于第K个自相关输出值，由前一次的自相关输出值以及k-2次输入信号、抽值因子利用上述的滑动窗方法得到。

5.如权利要求1所述的信号捕获方法，其特征在于，所述步骤B中通过如下步骤获取所述自相关输出值的幅度值：

分别获取自相关输出值的实部与虚部的绝对值，并从所述绝对值中选择出最大值以及最小值，再将该最大值加上最小值的二分之一的和，作为自相关输出值的幅度值acorr_amp。

6.如权利要求2所述的信号捕获方法，其特征在于，若步骤A中利用抽值因子获得信号子序列，相应地所述步骤C中，根据抽值因子F_D、终点位置end_idx及宽度值width确定信号序列的初始同步位置sync_idx为：

sync_idx＝F_D*(end_idx-floor(width/2))；

其中，所述floor()表示向下取整函数。

7.如权利要求1所述的信号捕获方法，其特征在于，所述方法进一步包括：

若捕获不成功，则依照继续进行捕获过程。

8.如权利要求1所述的信号捕获方法，其特征在于，所述方法进一步包括：

根据历史搜索结果调整所述步骤A、B、C中的搜索参数，所述搜索参数包括一次搜索的时间长度T_ONE、进行捕获的时间段的长度L_TIM、所述时间段L_TIM内的最多连续搜索次数N_SCH、以及连续成功搜索次数N_SUC。

9.如权利要求7或8所述的信号捕获方法，其特征在于，所述方法进一步包括：

若在所述时间段L_TIM内，最多连续搜索次数N_SCH与已进行的搜索次数之差小于或等于连续成功搜索次数N_SUC，则结束搜索。

10.如权利要求1所述的信号捕获方法，其特征在于，在单频干扰的信道环境下，所述步骤B中，通过直接搜索自相关输出幅度值的最大值来确定最宽的位置序列。

11.如权利要求6所述的信号捕获方法，其特征在于，所述方法进一步包括：

D、根据初始位置sync_idx、终点位置end_idx和宽度值width计算得到小数频偏。

12.如权利要求11所述的信号捕获方法，其特征在于，计算小数频偏的步骤具体包括：

D1、判断初始位置sync_idx与本次搜索的起始位置之差是否大于所设定的每次搜索的时间长度的一半，

若是，则计算频率偏差索引foe_idx位置为：

foe_idx＝end_idx-floor(width/2)+62500；

若不是，则计算频率偏差索引foe_idx位置为：

foe_idx＝end_idx-floor(width/2)+125000；

D2、在确定foe_idx位置后，进一步估计小数频率偏差为：

$\mathrm{\Δf}=\frac{\mathrm{angle}(\mathrm{acorr}\_\mathrm{out}(\mathrm{foe}\_\mathrm{idx}))}{2\mathrm{\π}}*\frac{{10}^7}{2048}$

其中，angle()表示取角度函数函数，floor()表示向下取整函数。

Images