一种基于查表法的北斗信号高精度捕获方法技术

本发明专利技术涉及一种基于查表法的北斗信号高精度捕获方法，属于卫星导航技术领域。本方法包括以下步骤：步骤一：剥离卫星信号伪码；步骤二：对信号做FFT；步骤三：设定门限值；步骤四：将FFT输出主峰值与门限值比较，若高于门限值说明捕获成功，捕获到的多普勒频移为fd；步骤五：将FFT输出次峰值与门限值比较，若高于门限值且次峰和主峰位置相邻，则可以对多普勒频移修正，否则将fd作为捕获结果；步骤六：以FFT输出主峰与次峰的比值作为地址查表，得到修正值fx；步骤七：判断主峰和次峰的位置关系，若次峰在主峰左侧，则修正后的多普勒频移为fd‑fx，否则为fd+fx。本方法应用在北斗信号的软件接收机中，能够在几乎没有增加计算负荷的情况下有效提高多普勒频移的捕获精度。

A High Accuracy Acquisition Method of Beidou Signal Based on Look-up Table Method

【技术实现步骤摘要】
一种基于查表法的北斗信号高精度捕获方法
本专利技术涉及一种基于查表法的北斗信号高精度捕获方法，属于卫星导航

技术介绍
随着北斗卫星导航系统的快速发展，使用北斗卫星进行高精度导航定位成为近几年的研究重点。信号捕获是北斗系统软件接收机中的关键技术，从PRN码(伪随机噪声码)、伪码相位和多普勒频移三个维度对信号进行捕获。由于捕获到的卫星信号中夹杂的噪声为高斯白噪声，单个复指数波形的最大似然频率估计值在接收信号的离散傅里叶变换(DFT)的峰值位置。为了更适合于计算机操作，快速傅里叶变换(FFT)通常用于单一的复指数波形处理。由于FFT输出的频谱是离散的，多普勒频移估计的精度由FFT的频率分辨率决定。卫星信号粗捕获后的载波频率精度无法保证输入信号进入跟踪环路的跟踪范围，必须对捕获后的粗略频率进行细化，使得载波频率落入跟踪环路的捕获带，才能使得跟踪环路正常工作。因此在卫星信号捕获后，还需要设计提高多普勒频移捕获精度的方法。现有的提高多普勒频移捕获精度的方法，无论是基于并行频率捕获法还是并行码相位捕获法，都是以增加额外的计算量为代价来提高FFT的分辨率，从而提高多普勒频移的捕获精度。现有的高精度捕获技术中，没有一种计算量小、有效的高精度捕获方法，造成现有高精度捕获方法消耗过多的运算资源。
技术实现思路
本专利技术提出一种基于查表法的北斗信号高精度捕获方法，将其应用在北斗信号的软件接收机中，能够在几乎没有增加计算负荷的情况下有效提高多普勒频移的捕获精度。本专利技术为解决其技术问题采用如下技术方案：一种基于查表法的北斗信号高精度捕获方法，包括以下步骤：步骤一：将输入的数字中频信号和本地伪码信号相关，剥离卫星信号中的伪随机码；步骤二：对剥离伪码后的信号进行快速傅里叶变换，通过FFT对多普勒频移值进行估计；步骤三：设定捕获门限值；步骤四：将FFT输出的主峰值结果与设定门限值进行比较，若高于门限值说明捕获成功，记录主峰位置，否则说明未捕获到卫星；步骤五：将FFT输出的次峰值结果与设定门限值进行比较，若高于门限值并且次峰值位置和主峰值位置相邻，说明对捕获到的多普勒频移值进行修正，否则直接将主峰值对应的多普勒频移值fd作为最终的捕获结果；步骤六：以FFT输出的主峰值与次峰值的比值作为地址进行查表，获得多普勒频移的修正值fx；步骤七：判断主峰值和次峰值的位置关系，若次峰值在主峰值左侧，则修正后的多普勒频移值为fd-fx，否则修正后的多普勒频移值为fd+fx。步骤一中所述输入的数字中频信号表示为：其中A是信号振幅；C(n)是伪随机码；D(n)是导航电文；fIF是载波中心频率；fd是多普勒频移；是初始载波相位；w(n)是高斯白噪声，n为输入的离散卫星信号的数据点中的第n个采样数据，exp表示以e为底的指数函数。步骤二中所述剥离伪码后的信号表示为：其中：i为虚数单位。步骤二中在对剥离伪码后的信号进行快速傅里叶变换前，首先要进行加窗处理，得到有限长信号S1(n)，过程表示为：其中，S(n)为伪码剥离后的卫星信号，TS为时域离散采样周期，n表示第n个采样数据,t为时间,L为FFT点数，δ(t-nTS)为冲激函数，为矩形函数。步骤三的具体过程如下：首先设置一个信号的捕获虚警率Pfa，再根据此虚警率的要求计算得到设定的门限值Yt，计算过程为：其中，σ为莱斯分布概率密度函数曲线的标准差。本专利技术的有益效果如下：本专利技术基于传统的并行频率捕获方法，在获得FFT输出结果后，记录主峰值和次峰值的大小和位置。若主峰值和次峰值均大于设定的门限值，并且主峰和次峰的位置相邻，就可以以主峰值和次峰值的比值为地址，通过查表的方式获得多普勒频移的修正值。基于查表法的北斗信号高精度捕获方法提高了载波频率的捕获精度，并且几乎没有增加计算负荷。将其应用在北斗信号的软件接收机中，当多普勒频移在表格的可查询范围内时，能够将多普勒频移的捕获精度提高到50Hz。附图说明图1是基于查表法的北斗信号高精度捕获方法的结构示意图。图2是本专利技术中基于查表法的多普勒频移修正方法流程图。图3(a)是并行频率捕获算法对多普勒频移为50Hz的北斗信号进行FFT得到的结果图；图3(b)是并行频率捕获算法对多普勒频移为100Hz的北斗信号进行FFT得到的结果图；图3(c)是并行频率捕获算法对多普勒频移为150Hz的北斗信号进行FFT得到的结果图；图3(d)是并行频率捕获算法对多普勒频移为200Hz的北斗信号进行FFT得到的结果图；图3(e)是并行频率捕获算法对多普勒频移为250Hz的北斗信号进行FFT得到的结果图；图3(f)是并行频率捕获算法对多普勒频移为300Hz的北斗信号进行FFT得到的结果图；图3(g)是并行频率捕获算法对多普勒频移为350Hz的北斗信号进行FFT得到的结果图；图3(h)是并行频率捕获算法对多普勒频移为400Hz的北斗信号进行FFT得到的结果图；图3(i)是并行频率捕获算法对多普勒频移为450Hz的北斗信号进行FFT得到的结果图；图3(j)是并行频率捕获算法对多普勒频移为500Hz的北斗信号进行FFT得到的结果图。图4是高精度多普勒频移捕获结果和多普勒频移真实值间的对比图。具体实施方式下面结合附图对本专利技术技术方案作进一步说明。基于查表法的北斗信号高精度捕获方法的结构示意图如图1所示。具体实施时，卫星信号经过射频前端处理后，输入捕获模块的卫星信号为复指数离散信号，其表达式为：其中A是信号振幅；C(n)是伪随机码；D(n)是导航电文；i为虚数单位,fIF是载波中心频率；fd是多普勒频移；是初始载波相位；w(n)是高斯白噪声，n为输入的离散卫星信号的数据点中的第n个采样数据，exp表示以e为底的指数函数。对该信号的处理过程为：第一步：将中频数字信号和本地伪码信号C(n-m)相乘，完成伪码的剥离，过程表示如下：S(n)＝SIF(n)C(n-m)其中，SIF(n)表示输入捕获模块的卫星信号,C(n-m)表示本地伪码信号，S(n)为伪码剥离后的卫星信号。m表示伪随机码的码片位置，因为北斗信号的码片长度为2046个码片，所以m＝0,1,2...2045。当m＝0时，本地伪码和输入信号的相位对齐，即可实现伪码剥离。只对1ms的卫星数据进行研究，由于用于做相干过程的信号长度远小于导航电文的周期长度，可以忽略导航电文对捕获过程的影响，所以剥离后的数字信号可以表示为：第二步：对信号进行FFT，首先对S(n)进行加矩形窗处理，得到信号S1(n)表示为：其中，TS为时域离散采样周期，t为时间，L为FFT点数，δ(t-nTS)为冲激函数，为矩形函数。将离散时域信号S1(n)转换为连续的频域信号S(f)，其表达式为：其中：w(f)为将时域上的噪声信号w(n)转换到频域的结果，f表示把时域信号变到频域信号之后的频率，k表示对频域信号离散加窗之后得到的信号的数据点中的第k个数据，sinc(LTSf)为采样函数。对信号S(f)进行离散和加窗处理，过程表示为：其中，为频域离散采样周期，S(k)为对信号S(f)进行离散和加窗后得到的有限长离散信号。第三步：设定捕获门限值。首先设置一个信号的捕获虚警率Pfa，再根据此虚警率的要求计算得到设定的门限值Yt，计算过程为：其中，σ为莱斯分布概率密度函数曲线的标准差。第四步：本文档来自技高网...

【技术保护点】
1.一种基于查表法的北斗信号高精度捕获方法，其特征在于，包括以下步骤：步骤一：将输入的数字中频信号和本地伪码信号相关，剥离卫星信号中的伪随机码；步骤二：对剥离伪码后的信号进行快速傅里叶变换，通过FFT对多普勒频移值进行估计；步骤三：设定捕获门限值；步骤四：将FFT输出的主峰值结果与设定门限值进行比较，若高于门限值说明捕获成功，记录主峰位置，否则说明未捕获到卫星；步骤五：将FFT输出的次峰值结果与设定门限值进行比较，若高于门限值并且次峰值位置和主峰值位置相邻，说明对捕获到的多普勒频移值进行修正，否则直接将主峰值对应的多普勒频移值fd作为最终的捕获结果；步骤六：以FFT输出的主峰值与次峰值的比值作为地址进行查表，获得多普勒频移的修正值fx；步骤七：判断主峰值和次峰值的位置关系，若次峰值在主峰值左侧，则修正后的多普勒频移值为fd‑fx，否则修正后的多普勒频移值为fd+fx。

【技术特征摘要】
1.一种基于查表法的北斗信号高精度捕获方法，其特征在于，包括以下步骤：步骤一：将输入的数字中频信号和本地伪码信号相关，剥离卫星信号中的伪随机码；步骤二：对剥离伪码后的信号进行快速傅里叶变换，通过FFT对多普勒频移值进行估计；步骤三：设定捕获门限值；步骤四：将FFT输出的主峰值结果与设定门限值进行比较，若高于门限值说明捕获成功，记录主峰位置，否则说明未捕获到卫星；步骤五：将FFT输出的次峰值结果与设定门限值进行比较，若高于门限值并且次峰值位置和主峰值位置相邻，说明对捕获到的多普勒频移值进行修正，否则直接将主峰值对应的多普勒频移值fd作为最终的捕获结果；步骤六：以FFT输出的主峰值与次峰值的比值作为地址进行查表，获得多普勒频移的修正值fx；步骤七：判断主峰值和次峰值的位置关系，若次峰值在主峰值左侧，则修正后的多普勒频移值为fd-fx，否则修正后的多普勒频移值为fd+fx。2.根据权利要求1所述的一种基于查表法的北斗信号高精度捕获方法，其特征在于，步骤一中所述输入的数字中频信号表示为：其中A是信号...

【专利技术属性】
技术研发人员：曾庆喜，陈则王，高唱，邱文旗，吕查德，刘德辉，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32