发明名称 一种GPS定位卫星L2C中频信号CL码捕获方法
摘要 本发明涉及一种GPS定位卫星L2C中频信号CL码捕获方法,基于折叠和平均采样思想,首先利用CL码信号不同段之间相互正交的特点,通过线性叠加形成折叠码后进行相关捕获,有效地减少了相关计算的次数,同时利用平均采样法减少每次傅里叶变换相关的采集点数目,并采用分裂基傅里叶变换代替传统的基2傅里叶变换进行相关运算,进一步减少了每次相关过程中所需要的计算量和存储空间,提高了捕获的实时性。
申请公布号 CN104007454B 申请公布日期 2016.04.27
申请号 CN201410191908.5 申请日期 2014.05.07
申请人 南京航空航天大学;江苏新安电器有限公司 发明人 曾庆喜;唐琳琳;徐亮;张鹏娜;杨英杰;唐靖
分类号 G01S19/30(2010.01)I 主分类号 G01S19/30(2010.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 杨楠
主权项 一种GPS定位卫星L2C中频信号CL码捕获方法,其特征在于,包括如下步骤:步骤001.接收定位卫星发送的一个周期的L2C中频信号,并与多普勒频移值ω<sub>D</sub>=(ω<sub>D</sub>)<sub>min</sub>下、与L2C中频信号长度相同的本地载波信号按时间序列一一对应相乘,获取信号S,其中,L2C中频信号的周期为T,(ω<sub>D</sub>)<sub>min</sub>为多普勒频移搜索范围的最小值;步骤002.将信号S按长度平均划分为K段,K≥2,各段信号的采样点数为N,任取其中一段信号,记为信号S1;然后按预设采样点数N<sub>d</sub>,通过平均采样法将信号S1的采样点数由N降至N<sub>d</sub>,记为平均采样信号S1<sub>cy</sub>,并通过补零操作将平均采样信号S1<sub>cy</sub>扩展为<img file="FDA0000901441340000015.GIF" wi="320" he="87" />其中,0表示与平均采样信号S1<sub>cy</sub>长度相同的零向量;步骤003.本地生成一个周期为T、对应上述发送L2C中频信号的定位卫星的本地CL码信号,且长度与上述L2C中频信号长度相同,将其按长度平均划分为K段,CL<sub>1</sub>、…CL<sub>i</sub>、…CL<sub>K</sub>,1≤i≤K;然后按首尾相连、依次将相邻两段信号组合成一个新码信号的方式,获得K个新码YHC<sub>i</sub>;按预设叠加段数M,通过线性叠加法依次将K个新码YHC<sub>i</sub>叠加构成n个折叠码HC<sub>s</sub>,1≤s≤n,n=ceil(K/M),ceil表示向上取整函数;按照步骤002中的平均采样法将各个折叠码HC<sub>s</sub>的采样点数由2N降至2N<sub>d</sub>,记为平均采样折叠码<img file="FDA0000901441340000011.GIF" wi="136" he="71" />步骤004.根据如下公式,将各个平均采样折叠码<img file="FDA0000901441340000012.GIF" wi="116" he="71" />分别和S<sub>zp</sub>按时间序列一一对应相关,获得相关结果记为YCorr<sub>s</sub>,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>YCorr</mi><mi>s</mi></msub><mo>=</mo><mi>S</mi><mi>I</mi><mi>F</mi><mi>F</mi><mi>T</mi><mo>{</mo><mi>S</mi><mi>F</mi><mi>F</mi><mi>T</mi><mrow><mo>(</mo><msubsup><mi>HC</mi><mi>s</mi><mrow><mi>c</mi><mi>y</mi></mrow></msubsup><mo>)</mo></mrow><mo>.</mo><mo>*</mo><mi>c</mi><mi>o</mi><mi>n</mi><mi>j</mi><mo>&lsqb;</mo><mi>S</mi><mi>F</mi><mi>F</mi><mi>T</mi><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>z</mi><mi>p</mi></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FDA0000901441340000013.GIF" wi="883" he="94" /></maths>其中,SFFT和SIFFT分别表示分裂基傅里叶变换和分裂基逆傅里叶变换,conj表示共轭运算;然后分别获取各个相关结果YCorr<sub>s</sub>中的前一半数据,即为各个平均采样折叠码<img file="FDA0000901441340000014.GIF" wi="119" he="70" />分别与S<sub>zp</sub>的最终相关结果Corr<sub>s</sub>,对n个最终相关结果Corr<sub>s</sub>进行比较,判断是否存在相关峰值,是则进入步骤006;否则进入下一步骤;步骤005.依次在{(ω<sub>D</sub>)<sub>min</sub>+Δω<sub>D</sub>、(ω<sub>D</sub>)<sub>min</sub>+2Δω<sub>D</sub>、…(ω<sub>D</sub>)<sub>max</sub>}范围中,改变多普勒频移值ω<sub>D</sub>的值,跳转至步骤001,其中,Δω<sub>D</sub>为多普勒频移的搜索长度,(ω<sub>D</sub>)<sub>max</sub>为多普勒频移搜索范围的最大值;步骤006.根据步骤004中获得的相关峰值,确定L2C中频信号CL码信号的捕获范围,实现CL码信号的捕获,获得CL码信号的初始相位。
地址 210016 江苏省南京市秦淮区御道街29号