发明内容
针对上述技术问题,本发明所要解决的技术问题是提供一种基于分段分裂基傅里叶变换,能够有效提高双频GPS信号捕获效率的基于分段分裂基傅里叶变换的双频GPS信号捕获方法。
本发明为了解决上述技术问题采用以下技术方案:本发明设计了一种基于分段分裂基傅里叶变换的双频GPS信号捕获方法,包括如下步骤:
步骤001.捕获GPS数字中频信号中L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A;
步骤002.根据L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A,获取GPS数字中频信号中L2C信号的载波多普勒频移fdr2、以及L2C信号中CM码的初始相位预估集合A,生成本地归零CM码信号;根据ωd(L2C/L1)=2πfdr2获得L2C信号的多普勒角频移ωd(L2C/L1);
步骤003.针对GPS数字中频信号和本地归零CM码信号的第m个相干积分块,分别进行相同的分段,并分别针对该两种信号的各分段信号进行分裂基傅里叶变换;分别针对经过分裂基傅里叶变换的各段本地归零CM码信号,取其共轭值;
步骤004.针对步骤003中,经过分裂基傅里叶变换的各段GPS数字中频信号和经过分裂基傅里叶变换、并取共轭值的各段本地归零CM码信号,按序列进行点对点相乘;并对按序列进行点对点相乘的结果做分裂基逆傅里叶变换,获得各段GPS数字中频信号与各段本地归零CM码信号之间按序列一一对应的CM码的相关积分ΔZCM(m,l),l=0,1…(l′-1),l'表示针对GPS数字中频信号和本地归零CM码信号的第m个相干积分块进行分段的段数;
步骤005.针对各个CM码的相关积分ΔZCM(m,l)进行移位相加,获得第m个相干积分块中CM码的相关积分ZCM(m);并对ZCM(m)取绝对值获得CM码的捕获判决函数
步骤006.根据CM码的捕获判决函数CM码的初始相位预估集合A、以及预设的检测门限η,获取CM码的初始相位
步骤007.根据CM码的初始相位CM码和CL码的相位关系,获得CL码的初始相位预估集合B,生成本地归零CL码信号,并按照步骤004至步骤006的方式获取CL码的初始相位。
作为本发明的一种优选技术方案:所述步骤001中,采用并行码相位搜索捕获算法捕获GPS数字中频信号中L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A。
作为本发明的一种优选技术方案:所述步骤001中,若采用并行码相位搜索捕获算法捕获GPS数字中频信号中L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A失败,则采用差分相干积分捕获算法进行获取。
作为本发明的一种优选技术方案:所述步骤002中,根据L1信号的载波多普勒频移fdr1和如下公式(1),获取GPS数字中频信号中L2C信号的载波多普勒频移fdr2,
其中,fL1、fL2分别为L1信号、L2C信号的载波频率;
根据C/A码的初始相位nC/A和如下公式(2),获取L2C信号中CM码的初始相位预估集合A,
A=nC/A+1023×(i-1),i=1、2…20(2)。
本发明所述一种基于分段分裂基傅里叶变换的双频GPS信号捕获方法采用以上技术方案与现有技术相比,具有以下技术效果:
本发明设计的基于分段分裂基傅里叶变换的双频GPS信号捕获方法,针对L2C信号捕获过程,将多普勒频移和码相位的二维搜索过程优化为一维搜索,并有效地减少了L2C信号的码相位搜索范围;同时采用分段分裂基FFT和分段分裂基IFFT进行相关运算,在很大程度上减少了捕获过程中的计算量和存储空间,提高了计算速度,为L1/L2C双频GPS信号接收机在嵌入式系统上的实现奠定了理论基础。
具体实施方式
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。
如图1和图2所示,本发明设计的基于分段分裂基傅里叶变换的双频GPS信号捕获方法在具体实际应用过程中,按如下步骤进行:
步骤001.采用并行码相位搜索捕获算法捕获GPS数字中频信号中L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A,若捕获失败,则采用差分相干积分捕获算法捕获GPS数字中频信号中L1信号的载波多普勒频移fdr1和C/A码的初始相位nC/A。
步骤002.根据L1信号的载波多普勒频移fdr1和如下公式(1),获取GPS数字中频信号中L2C信号的载波多普勒频移fdr2,
其中,fL1、fL2分别为L1信号、L2C信号的载波频率;
根据C/A码的初始相位nC/A和如下公式(2),获取L2C信号中CM码的初始相位预估集合A,
A=nC/A+1023×(i-1),i=1、2…20(2);
其中,CM码的初始相位存在20种预估值,然后生成本地归零CM码信号;
根据ωd(L2C/L1)=2πfdr2获得L2C信号的多普勒角频移ωd(L2C/L1)。
步骤003.针对GPS数字中频信号和本地归零CM码信号的第m个相干积分块,分别进行相同的分段,针对这两种信号分别分为l'段,且每段长度均相同,长度为L,分段完成后,通过补零的方式使各数据段的长度达到要求的长度;分别针对该两种信号的各分段信号进行分裂基傅里叶变换;分别针对经过分裂基傅里叶变换的各段本地归零CM码信号,取其共轭值;
针对GPS数字中频信号中第m个相干积分块的第l段信号yk进行分裂基傅里叶变换运算,如式(3)所示:
针对本地归零CM码信号中第m个相干积分块的第l段信号xk进行分裂基傅里叶变换运算,并取其共轭值,如式(4)所示:
其中,ncoh为第m个相干积分块信号总的数据长度,(为简化,假设ncoh能整除l′,即L=ncoh/l′)。
步骤004.针对步骤003中,经过分裂基傅里叶变换的各段GPS数字中频信号和经过分裂基傅里叶变换、并取共轭值的各段本地归零CM码信号,按序列进行点对点相乘;并对按序列进行点对点相乘的结果做分裂基逆傅里叶变换,获得各段GPS数字中频信号与各段本地归零CM码信号之间按序列一一对应的CM码的相关积分ΔZCM(m,l),l=0,1…(l′-1),l'表示针对GPS数字中频信号和本地归零CM码信号的第m个相干积分块进行分段的段数;
其中,第l段GPS数字中频信号和第l段本地归零CM码信号分别经过分裂基傅里叶变换后进行点对点相乘的结果做分裂基逆傅里叶变换,获得对应的CM码的相关积分ΔZCM(m,l),如下式(5)所示:
其中,式(3)-式(5)中N1=(2n)min≥(2L-1),n为正整数;l满足l=0,1…(l′-1);式(5)中的ΔZdiscard为要舍弃的相关积分结果。
分裂基逆傅里叶变换的一般计算方法通过下式进行运算:
其中,N为样本个数,
其中,相干积分块上彼此相对应的GPS数字中频信号和本地归零CM码信号分别经过分裂基傅里叶变换后进行点对点相乘的结果做分裂基逆傅里叶变换,获得对应的CM码的相关积分结果包括同向和正交分量,如第m个相干积分块上第l段GPS数字中频信号和第l段本地归零CM码信号分别经过分裂基傅里叶变换后进行点对点相乘的结果做分裂基逆傅里叶变换,获得对应的CM码的相关积分的结果包括同向和正交分量,其中,同向分量如下式(6)所示:
其中,正交分量如下式(7)所示:
其中,ωd(L2C/L1)=2πfdr2为L2C信号的多普勒角频移,和分别表示同向分量和正交分量,l=0,1…(l′-1),ncoh为第m个相干积分块信号总的数据长度,(为简化,假设ncoh能整除l′,即L=ncoh/l′),tk为采样时刻点,为CM码的初始相位,ωL2C_IF为L2C信号的中频。
针对第l段GPS数字中频信号和第l段本地归零CM码信号点对点相乘的结果做分裂基逆傅里叶变换,获得对应的CM码的相关积分ΔZCM(m,l)如下式(8)所示:
其中,j表示复数。
步骤005.针对各个CM码的相关积分ΔZCM(m,l)进行移位相加,获得第m个相干积分块中CM码的相关积分ZCM(m);并对ZCM(m)取绝对值获得CM码的捕获判决函数
其中,将ΔZCM(m,l)进行移位相加获得第m个相干积分块中CM码的相关积分ZCM(m),如下式(9)所示:
其中,Σ表示移位相加,ICM(m)和QCM(m)分别表示第m个相干积分块GPS数字中频信号和本地归零CM码信号点对点相乘的结果做分裂基逆傅里叶变换,获得对应的CM码的相关积分的结果的同向积分和正交积分。
最终获得的CM码的捕获判决函数如下式(10)所示:
步骤006.根据CM码的捕获判决函数CM码的初始相位预估集合A、以及预设的检测门限η,获取CM码的初始相位
其中,依次将CM码的初始相位预估集合A中的CM码初始相位预估值代入式(10)中,并根据预设的检测门限η,按照|Z|≥η规则,获得CM码的初始相位即实现CM码的捕获。
步骤007.根据CM码的初始相位CM码和CL码的相位关系,获得CL码的初始相位预估集合B,生成本地归零CL码信号,并按照步骤004至步骤006的方式获取CL码的初始相位,其中,由于每个周期的CL码正好包含75个周期的CM码,所以可以通过CM码和CL码的相位联系获得75种CL码可能的初始相位;并且由于L2C码由CM码和CL码依次相隔组成,因此CL码载波多普勒频移与CM信号的载波多普勒频移相同。
基于以上步骤001至步骤007,最终实现针对双频GPS信号的捕获,针对L2C信号捕获过程,将多普勒频移和码相位的二维搜索过程优化为一维搜索,并有效地减少了L2C信号的码相位搜索范围;同时采用分段分裂基FFT和分段分裂基IFFT进行相关运算,在很大程度上减少了捕获过程中的计算量和存储空间,提高了计算速度,为L1/L2C双频GPS信号接收机在嵌入式系统上的实现奠定了理论基础。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。