CN104007454A - 一种gps定位卫星l2c中频信号cl码捕获方法 - Google Patents

Abstract

本发明涉及一种GPS定位卫星L2C中频信号CL码捕获方法，基于折叠和平均采样思想，首先利用CL码信号不同段之间相互正交的特点，通过线性叠加形成折叠码后进行相关捕获，有效地减少了相关计算的次数，同时利用平均采样法减少每次傅里叶变换相关的采集点数目，并采用分裂基傅里叶变换代替传统的基2傅里叶变换进行相关运算，进一步减少了每次相关过程中所需要的计算量和存储空间，提高了捕获的实时性。

Description

一种GPS定位卫星L2C中频信号CL码捕获方法

技术领域

本发明涉及一种GPS定位卫星L2C中频信号CL码捕获方法。

背景技术

全球卫星定位系统(GPS)作为目前世界上应用范围最广泛、实用性最强的全球精密授时、测距、导航、定位系统，已广泛应用于汽车导航、目标监控、测绘等多个领域。但是随着时代的发展和用户对定位精度要求的不断提高，传统的L1民用单频GPS系统已经显示出一定的局限性。GPS现代化计划在L2波段上增加了由CM码和CL码经时分复用形成的L2C民用信号，减少了L1频段的C/A码在结构上的限制，而且与传统的民用信号C/A码相比，L2C信号具有更低的载波跟踪门限和数据解调门限，其伪码具有更优的相关性能，更适合在丛林、室内等信号较弱的场合应用。尤其是L2C信号中的CL码，由于该码没有调制导航数据，对积分长度没有限制，使其更适合对接收机灵敏度要求很高的情况。但是，CL码的长度远远大于C/A码以及CM码，如果采用传统的捕获算法对CL码直接进行捕获，现有的硬件和软件资源都无法满足捕获计算的要求。通常都是先进行CM码的捕获，然后利用CM码辅助CL码实现捕获，但是当CM码很难甚至无法被捕获到的时候，就需要一种可以对CL码进行快速直接捕获的算法。

发明内容

本发明所要解决的技术问题是提供一种无需借助其它码，能够直接快速的对L2C中频信号CL码实现捕获的GPS定位卫星L2C中频信号CL码捕获方法。

本发明为了解决上述技术问题采用以下技术方案：本发明设计了一种GPS定位卫星L2C中频信号CL码捕获方法，包括如下步骤：

步骤001.接收定位卫星发送的一个周期的L2C中频信号，并与预设多普勒频移值ω_D＝(ω_D)_min下、与L2C中频信号长度相同的本地载波信号按时间序列一一对应相乘，获取信号S，其中，L2C中频信号的周期为T，(ω_D)_min为多普勒频移搜索范围的最小值；

步骤002.将信号S按长度平均划分为K段，K≥2，各段信号的采样点数为N，任取其中一段信号，记为信号S1；然后按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d，记为平均采样信号S1_cy，并通过补零操作将平均采样信号S1_cy扩展为其中，0表示与平均采样信号S1_cy长度相同的零向量；

步骤003.本地生成一个周期为T、对应上述发送L2C中频信号的定位卫星的本地CL码信号，且长度与上述L2C中频信号长度相同，将其按长度平均划分为K段，CL₁、…CL_i、…CL_K，1≤i≤K；然后按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得K个新码YHC_i；按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s，1≤s≤n，n＝ceil(K/M)，ceil表示向上取整函数；按照步骤002中的平均采样法将各个折叠码HC_s的采样点数由2N降至2N_d，记为平均采样折叠码 ${\mathrm{HC}}_s^{\mathrm{cy}};$

步骤004.根据如下公式，将各个平均采样折叠码分别和S_zp按时间序列一一对应相关，获得相关结果记为YCorr_s，

$Y{\mathrm{Corr}}_s=\mathrm{SIFFT}\{\mathrm{SFFT}\left({\mathrm{HC}}_s^{\mathrm{cy}}\right).*\mathrm{conj}[\mathrm{SFFT}\left(S_{\mathrm{zp}}\right)\left]\right\}$

其中，SFFT和SIFFT分别表示分裂基傅里叶变换和分裂基逆傅里叶变换，conj表示共轭运算；然后分别获取各个相关结果YCorr_s中的前一半数据，即为各个平均采样折叠码分别与S_zp的最终相关结果Corr_s，对n个最终相关结果Corr_s进行比较，判断是否存在相关峰值，是则进入步骤006；否则进入下一步骤；

步骤005.依次在{(ω_D)_min+Δω_D、(ω_D)_min+2Δω_D、…(ω_D)_max}范围中，改变多普勒频移值ω_D的值，跳转至步骤001，其中，Δω_D为多普勒频移的搜索长度，(ω_D)_max为多普勒频移搜索范围的最大值；

步骤006.根据步骤004中获得的相关峰值，确定L2C中频信号CL码信号的捕获范围，实现CL码信号的捕获，获得CL码信号的初始相位。

作为本发明的一种优选技术方案：所述步骤002中，按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d的具体过程如下：

设计将信号S1的采样点数由N降至N_d需要进行N_p次p个采样点的平均和N_q次q个采样点的平均，即针对信号S1中的采样点依序先将(N_p*p)个采样点，每p个采样点累加为1个采样点，重新获得N_p个采样点；再针对剩余的(N_q*q)个采样点，每q个采样点累加为1个采样点，重新获得N_q个采样点；其中，N_d＝N_p+N_q，p＝ceil(N/N_d)，q＝floor(N/N_d)，floor表示向下取整函数。

作为本发明的一种优选技术方案：所述步骤003中，针对本地CL码信号划分的K段信号，按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得K个新码YHC_i的具体过程如下：

当i＜K时，YHC_i＝[CL_i CL_i+1]；当i＝K时，YHC_K＝[CL_K CL₁]。

作为本发明的一种优选技术方案：所述步骤003中，按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s的具体过程如下：

当s＜n时， ${\mathrm{HC}}_s=\underset{j=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)};$ 当s＝n时， ${\mathrm{HC}}_n=\underset{j=1}{\overset{(\mathrm{Mn}+K)}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)}.$

作为本发明的一种优选技术方案：所述步骤006具体包括如下步骤：

步骤00601.根据所述最终相关结果Corr_s比较中相关峰值所对应的第crw个平均采样折叠码1≤crw≤n，获取该平均采样折叠码所对应线性叠加的新码YHC_i，且该平均采样折叠码所对应线性叠加的新码YHC_i的个数为L个，然后根据如下公式，将该L个新码YHC_i分别与[S1 0]按时间序列一一对应相关，获得相关结果记为YRCorr_i，

针对L个相关结果YRCorr_i，分别获取各个相关结果YRCorr_i中的前一半数据，即为第crw个平均采样折叠码所对应线性叠加的各个新码YHC_i分别与S_zp的最终相关结果RCorr_i，对L个最终相关结果RCorr_i进行比较，捕获相关峰值，获取该相关峰值所对应的第cr个最终相关结果RCorr_cr中第cn个采样点，1≤cr≤L，1≤cn≤N；

步骤00602.根据ccn＝((crw-1)×M+cr-1)×N+cn，获得CL码信号的初始相位ccn。

本发明所述一种GPS定位卫星L2C中频信号CL码捕获方法采用以上技术方案与现有技术相比，具有以下技术效果：本发明设计的GPS定位卫星L2C中频信号CL码捕获方法，基于折叠和平均采样思想，首先利用CL码信号不同段之间相互正交的特点，通过线性叠加形成折叠码后进行相关捕获，有效地减少了相关计算的次数，同时利用平均采样法减少每次傅里叶变换相关的采集点数目，并采用分裂基傅里叶变换代替传统的基2傅里叶变换进行相关运算，进一步减少了每次相关过程中所需要的计算量和存储空间，提高了捕获的实时性，也为双频GPS信号接收机在资源有限的嵌入式系统上的实现奠定了一定的理论基础。

附图说明

图1是本发明设计的GPS定位卫星L2C中频信号CL码捕获方法的流程示意图；

图2是本发明设计中平均采样示意图。

具体实施方式

下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。

如图1所示，本发明设计了一种GPS定位卫星L2C中频信号CL码捕获方法，包括如下步骤：

步骤001.接收定位卫星发送的一个周期的L2C中频信号，并与预设多普勒频移值ω_D＝(ω_D)_min下、与L2C中频信号长度相同的本地载波信号按时间序列一一对应相乘，获取信号S，其中，L2C中频信号的周期为T，(ω_D)_min为多普勒频移搜索范围的最小值；

步骤002.将信号S按长度平均划分为K段，K≥2，各段信号的采样点数为N，任取其中一段信号，记为信号S1；然后按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d，记为平均采样信号S1_cy，其中，如图2所示，平均采样的具体过程如下：

设计将信号S1的采样点数由N降至N_d需要进行N_p次p个采样点的平均和N_q次q个采样点的平均，即针对信号S1中的采样点依序先将(N_p*_p)个采样点，每p个采样点累加为1个采样点，重新获得N_p个采样点；再针对剩余的(N_q*_q)个采样点，每q个采样点累加为1个采样点，重新获得N_q个采样点；其中，N_d＝N_p+N_q，p＝ceil(N/N_d)，q＝floor(N/N_d)，floor表示向下取整函数；

并通过补零操作将平均采样信号S1_cy扩展为其中，0表示与平均采样信号S1_cy长度相同的零向量；

步骤003.本地生成一个周期为T、对应上述发送L2C中频信号的定位卫星的本地CL码信号，且长度与上述L2C中频信号长度相同，将其按长度平均划分为K段，CL₁、…CL_i、…CL_K，1≤i≤K；然后按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得K个新码YHC_i，具体过程如下：

当i＜K时，YHC_i＝[CL_i CL_i+1]；当i＝K时，YHC_K＝[CL_K CL₁]；

按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s，1≤s≤n，n＝ceil(K/M)，ceil表示向上取整函数，具体过程如下：

当s＜n时， ${\mathrm{HC}}_s=\underset{j=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)};$ 当s＝n时， ${\mathrm{HC}}_n=\underset{j=1}{\overset{(\mathrm{Mn}+K)}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)}.$

按照步骤002中的平均采样法将各个折叠码HC_s的采样点数由2N降至2N_d，记为平均采样折叠码

步骤004.根据如下公式，将各个平均采样折叠码分别和S_zp按时间序列一一对应相关，获得相关结果记为YCorr_s，

$Y{\mathrm{Corr}}_s=\mathrm{SIFFT}\{\mathrm{SFFT}\left({\mathrm{HC}}_s^{\mathrm{cy}}\right).*\mathrm{conj}[\mathrm{SFFT}\left(S_{\mathrm{zp}}\right)\left]\right\}$

其中，SFFT和SIFFT分别表示分裂基傅里叶变换和分裂基逆傅里叶变换，conj表示共轭运算；然后分别获取各个相关结果YCorr_s中的前一半数据，即为各个平均采样折叠码分别与S_zp的最终相关结果Corr_s，对n个最终相关结果Corr_s进行比较，判断是否存在相关峰值，是则进入步骤006；否则进入下一步骤；

步骤005.依次在{(ω_D)_min+Δω_D、(ω_D)_min+2Δω_D、…(ω_D)_max}范围中，改变多普勒频移值ω_D的值，跳转至步骤001，其中，Δω_D为多普勒频移的搜索长度，(ω_D)_max为多普勒频移搜索范围的最大值；

步骤006.根据步骤004中获得的相关峰值，确定L2C中频信号CL码信号的捕获范围，实现CL码信号的捕获，获得CL码信号的初始相位，具体包括如下步骤：

步骤00601.根据所述最终相关结果Corr_s比较中相关峰值所对应的第crw个平均采样折叠码1≤crw≤n，获取该平均采样折叠码所对应线性叠加的新码YHC_i，且该平均采样折叠码所对应线性叠加的新码YHC_i的个数为L个，然后根据如下公式，将该L个新码YHC_i分别与[S1 0]按时间序列一一对应相关，获得相关结果记为YRCorr_i，

针对L个相关结果YRCorr_i，分别获取各个相关结果YRCorr_i中的前一半数据，即为第crw个平均采样折叠码所对应线性叠加的各个新码YHC_i分别与S_zp的最终相关结果RCorr_i，对L个最终相关结果RCorr_i进行比较，捕获相关峰值，获取该相关峰值所对应的第cr个最终相关结果RCorr_cr中第cn个采样点，1≤cr≤L，1≤cn≤N；

步骤00602.根据ccn＝((crw-1)×M+cr-1)×N+cn，获得CL码信号的初始相位ccn。

本发明设计了一种GPS定位卫星L2C中频信号CL码捕获方法在实际应用过程当中，具体包括如下步骤：

步骤001.接收定位卫星发送的一个周期的L2C中频信号，并与预设多普勒频移值ω_D＝(ω_D)_min下、与L2C中频信号长度相同的本地载波信号按时间序列一一对应相乘，获取信号S，其中，L2C中频信号的周期为T＝1.5s，(ω_D)_min为多普勒频移搜索范围的最小值；

步骤002.将信号S按长度平均划分为75段，每段长20ms，各段长20ms信号的采样点数为N，任取其中一段信号，记为信号S1；然后按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d，记为平均采样信号S1_cy，此时，信号S1_cy的长度为20ms，其中，平均采样的具体过程如下：

设计将信号S1的采样点数由N降至N_d需要进行N_p次p个采样点的平均和N_q次q个采样点的平均，即针对信号S1中的采样点依序先将(N_p*_p)个采样点，每p个采样点累加为1个采样点，重新获得N_p个采样点；再针对剩余的(N_q*q)个采样点，每q个采样点累加为1个采样点，重新获得N_q个采样点；其中，N_d＝N_p+N_q，p＝ceil(N/N_d)，q＝floor(N/N_d)，floor表示向下取整函数；

并通过补零操作将平均采样信号S1_cy扩展为其中，0表示与平均采样信号S1_cy长度相同的零向量，则S_zp的长度为40ms；

步骤003.本地生成一个周期为T＝1.5s、对应上述发送L2C中频信号的定位卫星的本地CL码信号，且长度与上述L2C中频信号长度相同，将其按长度平均划分为75段，CL₁、…CL_i、…CL₇₅，1≤i≤75；然后按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得75个新码YHC_i，各个新码YHC_i的长度为40ms，具体过程如下：

当i＜75时，YHC_i＝[CL_i CL_i+1]；当i＝75时，YHC₇₅＝[CL₇₅CL₁]；

按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s，各个折叠码HC_s的长度为40ms，1≤s≤n，n＝ceil(75/M)，ceil表示向上取整函数，例如ceil(3.14)＝4，具体过程如下：

当s＜n时， ${\mathrm{HC}}_s=\underset{j=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)};$ 当s＝n时， ${\mathrm{HC}}_n=\underset{j=1}{\overset{(\mathrm{Mn}+75)}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)};$

按照步骤002中的平均采样法将各个折叠码HC_s的采样点数由2N降至2N_d，记为平均采样折叠码各个平均采样折叠码的长度为40ms；

步骤004.根据如下公式，将各个平均采样折叠码分别和S_zp按时间序列一一对应相关，获得相关结果记为YCorr_s，

$Y{\mathrm{Corr}}_s=\mathrm{SIFFT}\{\mathrm{SFFT}\left({\mathrm{HC}}_s^{\mathrm{cy}}\right).*\mathrm{conj}[\mathrm{SFFT}\left(S_{\mathrm{zp}}\right)\left]\right\}$

其中，SFFT和SIFFT分别表示分裂基傅里叶变换和分裂基逆傅里叶变换，conj表示共轭运算；然后分别获取各个相关结果YCorr_s中的前一半数据，即为各个平均采样折叠码分别与S_zp的最终相关结果Corr_s，对n个最终相关结果Corr_s进行比较，判断是否存在相关峰值，是则进入步骤006；否则进入下一步骤；

步骤005.依次在{(ω_D)_min+Δω_D、(ω_D)_min+2Δω_D、…(ω_D)_max}范围中，改变多普勒频移值ω_D的值，跳转至步骤001，其中，Δω_D为多普勒频移的搜索长度，(ω_D)_max为多普勒频移搜索范围的最大值；

步骤006.根据步骤004中获得的相关峰值，确定L2C中频信号CL码信号的捕获范围，实现CL码信号的捕获，获得CL码信号的初始相位，具体包括如下步骤：

步骤00601.根据所述最终相关结果Corr_s比较中相关峰值所对应的第crw个平均采样折叠码1≤crw≤n，获取该平均采样折叠码所对应线性叠加的新码YHC_i，且该平均采样折叠码所对应线性叠加的新码YHC_i的个数为L个，然后根据如下公式，将该L个新码YHC_i分别与[S1 0]按时间序列一一对应相关，获得相关结果记为YRCorr_i，其中[S1 0]为通过补零操作将信号S1扩展为[S1 0]，其中，0表示与平均采样信号S1长度相同的零向量；

针对L个相关结果YRCorr_i，分别获取各个相关结果YRCorr_i中的前一半数据，即为第crw个平均采样折叠码所对应线性叠加的各个新码YHC_i分别与S_zp的最终相关结果RCorr_i，对L个最终相关结果RCorr_i进行比较，捕获相关峰值，获取该相关峰值所对应的第cr个最终相关结果RCorr_cr中第cn个采样点，1≤cr≤L，，最终相关结果RCorr_i的长度为20ms，对应采样点数目为N，1≤cn≤N；

步骤00602.根据ccn＝((crw-1)×M+cr-1)×N+cn，获得CL码信号的初始相位ccn。

本发明设计的GPS定位卫星L2C中频信号CL码捕获方法，基于折叠和平均采样思想，首先利用CL码信号不同段之间相互正交的特点，通过线性叠加形成折叠码后进行相关捕获，有效地减少了相关计算的次数，同时利用平均采样法减少每次傅里叶变换相关的采集点数目，并采用分裂基傅里叶变换代替传统的基2傅里叶变换进行相关运算，进一步减少了每次相关过程中所需要的计算量和存储空间，提高了捕获的实时性，也为双频GPS信号接收机在资源有限的嵌入式系统上的实现奠定了一定的理论基础。

上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (5)

1.一种GPS定位卫星L2C中频信号CL码捕获方法，其特征在于，包括如下步骤：

步骤001.接收定位卫星发送的一个周期的L2C中频信号，并与预设多普勒频移值ω_D＝(ω_D)_min下、与L2C中频信号长度相同的本地载波信号按时间序列一一对应相乘，获取信号S，其中，L2C中频信号的周期为T，(ω_D)_min为多普勒频移搜索范围的最小值；

步骤002.将信号S按长度平均划分为K段，K≥2，各段信号的采样点数为N，任取其中一段信号，记为信号S1；然后按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d，记为平均采样信号S1_cy，并通过补零操作将平均采样信号S1_cy扩展为其中，0表示与平均采样信号S1_cy长度相同的零向量；

步骤003.本地生成一个周期为T、对应上述发送L2C中频信号的定位卫星的本地CL码信号，且长度与上述L2C中频信号长度相同，将其按长度平均划分为K段，CL₁、…CL_i、…CL_K，1≤i≤K；然后按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得K个新码YHC_i；按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s，1≤s≤n，n＝ceil(K/M)，ceil表示向上取整函数；按照步骤002中的平均采样法将各个折叠码HC_s的采样点数由2N降至2N_d，记为平均采样折叠码

步骤004.根据如下公式，将各个平均采样折叠码分别和S_zp按时间序列一一对应相关，获得相关结果记为YCorr_s，

$Y{\mathrm{Corr}}_s=\mathrm{SIFFT}\{\mathrm{SFFT}\left({\mathrm{HC}}_s^{\mathrm{cy}}\right).*\mathrm{conj}[\mathrm{SFFT}\left(S_{\mathrm{zp}}\right)\left]\right\}$

其中，SFFT和SIFFT分别表示分裂基傅里叶变换和分裂基逆傅里叶变换，conj表示共轭运算；然后分别获取各个相关结果YCorr_s中的前一半数据，即为各个平均采样折叠码分别与S_zp的最终相关结果Corr_s，对n个最终相关结果Corr_s进行比较，判断是否存在相关峰值，是则进入步骤006；否则进入下一步骤；

步骤005.依次在{(ω_D)_min+Δω_D、(ω_D)_min+2Δω_D、…(ω_D)_max}范围中，改变多普勒频移值ω_D的值，跳转至步骤001，其中，Δω_D为多普勒频移的搜索长度，(ω_D)_max为多普勒频移搜索范围的最大值；

步骤006.根据步骤004中获得的相关峰值，确定L2C中频信号CL码信号的捕获范围，实现CL码信号的捕获，获得CL码信号的初始相位。

2.根据权利要求1所述一种GPS定位卫星L2C中频信号CL码捕获方法，其特征在于：所述步骤002中，按预设采样点数N_d，通过平均采样法将信号S1的采样点数由N降至N_d的具体过程如下：

设计将信号S1的采样点数由N降至N_d需要进行N_p次p个采样点的平均和N_q次q个采样点的平均，即针对信号S1中的采样点依序先将(N_p*p)个采样点，每p个采样点累加为1个采样点，重新获得N_p个采样点；再针对剩余的(N_q*q)个采样点，每q个采样点累加为1个采样点，重新获得N_q个采样点；其中，N_d＝N_p+N_q，p＝ceil(N/N_d)，q＝floor(N/N_d)，floor表示向下取整函数。

3.根据权利要求1所述一种GPS定位卫星L2C中频信号CL码捕获方法，其特征在于：所述步骤003中，针对本地CL码信号划分的K段信号，按首尾相连、依次将相邻两段信号组合成一个新码信号的方式，获得K个新码YHC_i的具体过程如下：

当i＜K时，YHC_i＝[CL_i CL_i+1]；当i＝K时，YHC_K＝[CL_K CL₁]。

4.根据权利要求1所述一种GPS定位卫星L2C中频信号CL码捕获方法，其特征在于：所述步骤003中，按预设叠加段数M，通过线性叠加法依次将K个新码YHC_i叠加构成n个折叠码HC_s的具体过程如下：

当s＜n时， ${\mathrm{HC}}_s=\underset{j=1}{\overset{M}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)};$ 当s＝n时， ${\mathrm{HC}}_n=\underset{j=1}{\overset{(\mathrm{Mn}+K)}{\mathrm{\Σ}}}\mathrm{YH}{C}_{((s-1)\×M+j)}.$

5.根据权利要求1所述一种GPS定位卫星L2C中频信号CL码捕获方法，其特征在于：所述步骤006具体包括如下步骤：

步骤00601.根据所述最终相关结果Corr_s比较中相关峰值所对应的第crw个平均采样折叠码1≤crw≤n，获取该平均采样折叠码所对应线性叠加的新码YHC_i，且该平均采样折叠码所对应线性叠加的新码YHC_i的个数为L个，然后根据如下公式，将该L个新码YHC_i分别与[S1 0]按时间序列一一对应相关，获得相关结果记为YRCorr_i，

针对L个相关结果YRCorr_i，分别获取各个相关结果YRCorr_i中的前一半数据，即为第crw个平均采样折叠码所对应线性叠加的各个新码YHC_i分别与S_zp的最终相关结果RCorr_i，对L个最终相关结果RCorr_i进行比较，捕获相关峰值，获取该相关峰值所对应的第cr个最终相关结果RCorr_cr中第cn个采样点，1≤cr≤L，1≤cn≤N；

步骤00602.根据ccn＝((crw-1)×M+cr-1)×N+cn，获得CL码信号的初始相位ccn。