一种北斗导航弱信号接收的频率跟踪方法与流程

本发明属于卫星导航领域，特别是涉及一种北斗导航弱信号接收的载波频率跟踪方法。

背景技术：

卫星导航系统已经成为当今世界国家综合国力及科学技术发展水平的重要标志，是国家经济基础的重要组成部分，它已经进入人们的日常生活中，与社会发展和经济建设息息相关。我国也已经投入大量的资金和人力，积极进行北斗卫星导航系统(beidounavigationsatellitesystem，bds)的研究。到2020年左右，我国将建成覆盖全球的北斗卫星导航系统，其将成为继美国gps、俄罗斯glonass和欧洲galileo之后的又一个全球全天候的卫星导航系统。

随着卫星导航应用的进一步发展，要求接收机能工作于更为严苛的应用环境下，例如，城市峡谷、林荫大道等卫星信号严重遮挡环境。所谓的高灵敏度接收机(hs-gnss,highsensitivity-globalnavigationsatellitesystems),就是为满足此类对弱信号接收的需求。载波跟踪环路是弱信号接收的关键环节，其性能对高灵敏度接收机具有重要影响。通常，通过加长相干积分时间，来抵抗弱信号中噪声的影响，以提高载波跟踪环路性能。基于fft鉴频的载波频率跟踪，利用了fft频谱分析，得到频率误差信号。同时，fft变换也是一个相干积分过程。然而，受限于fft实现的复杂度以及频率跟踪环路迟延，fft鉴频算法的相干积分时间不能太长。另一方面，较短的相干积分时间会带来fft频谱分析分辨率低的问题。

技术实现要素：

本发明提供一种北斗导航弱信号接收的频率跟踪方法，可用于弥补fft鉴频器频率分辨率低的不足，改善频率跟踪环路的稳定性。

本发明提出的一种北斗导航弱信号接收的频率跟踪方法，包括以下步骤：

s1、接收信号经过伪码解扩以及载波剥离后，进入相关器进行累加运算，取得相关值；

s2、根据累加运算结果，进行电文自预测，估计得出发送的电文数据比特，取得电文剥离所必需的电文比特；

s3、使用估计得出的电文数据比特，对接收信号进行电文剥离，以去除电文数据调制对相干积分的影响；

s4、对电文剥离后的数据进行缓存，缓存到预定的n个数据；

s5、对所述n个数据进行fft鉴频，得到粗鉴频的频率误差值df1；

s6、对所述n个数据进行超前滞后非相干鉴频，得到细鉴频的频率误差值df2；

s7、将粗鉴频的频率误差值df1和细鉴频的频率误差值df2进行相加得到总频率误差值；

s8、将所述总频率误差值送入环路滤波器进行滤波；

s9、环路滤波器输出信号驱动数控振荡器，进行频率跟踪。

优选地，s5包括以下步骤：s5.1、对输入的所述n个数据进行fft变换；s5.2、对fft变换输出进行模平方运算；s5.3、从模平方后的n个数据中搜索最大值，并得到最大值所对应的索引idx，其中，当索引idx>n/2时，idx＝idx–n；s5.4、计算粗鉴频的频率误差df1＝idx/(tcoh*n)，其中，tcoh表示相关积分时间。

优选地，s6包括以下步骤：s6.1、利用上述的频率误差信号df1和超前滞后非相干鉴频器的频率鉴别间隔fsh，进行频率搬移；s6.2、对频率搬移后的信号进行相干累加得到xl＝∑zl(i)，xr＝∑zr(i)，i＝1，2，…n；其中，z1表示超前支路进行频率搬移后的信号值，zr表示滞后支路进行频率搬移后的信号值；s6.3、累加后的信号再进行模平方el＝|xl|²，er＝|xr|²；s6.4、对模平方输出信号进行频率误差计算，得到细鉴频的频率误差df2。

优选地，s6.1中超前支路进行频率搬移后的公式为：zl＝y*exp(-j2pi(df1–fsh))；滞后支路进行频率搬移后的公式为：zr＝y*exp(-j2pi(df1+fsh))，其中，y＝yi+jyq为相关器输出的复数信号，pi为圆周率，fsh为超前滞后非相干鉴频器的频率鉴别间隔。

优选地，s6.1中频率鉴别间隔fsh满足条件fsh≤1/(tcoh*n)。

优选地，s6.4中细鉴频的超前滞后非相干鉴频的频率误差计算公式为：

一种计算机可读存储介质，其存储有与计算设备结合使用的计算机程序，所述计算机程序可被处理器执行以上述任一项所述方法。

本发明的有益效果包括：通过采用进行过频率误差补偿的超前滞后非相干鉴频，频率跟踪环路更能细致地实施频率跟踪，提高了环路的稳定性。

附图说明

图1是本发明实施例的载波频率跟踪环路框图。

图2是本发明实施例的fft鉴频器的处理流程图。

图3是本发明实施例的超前滞后非相干鉴频器的处理流程图。

图4是本发明实施例的频率跟踪环路的多普勒频率仿真结果示意图。

具体实施方式

下面结合具体实施方式并对照附图对本发明作进一步详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。

如图1所示，为本发明实施例的的载波频率跟踪环路框图。接收信号经过伪码解扩以及载波剥离后，送入相关器，相关积分时间(tcoh)为一个电文比特，即相关积分时间tcoh＝20ms。

相关器输出的相关值，送入电文自预测，估计出发送的电文数据比特dn。使用自预测电文比特dn对接收信号进行电文剥离，以去除电文数据调制对后续的相干积分的影响。电文自预测算法，可以采用已知的最大似然估计(maximumlikelihood，ml)的算法，或者差分解调估计算法，这里不作限定。

电文剥离后的数据送入缓存器，缓存到预定的n个数据后，送入fft鉴频器和超前滞后非相干鉴频器。

如图2所示为fft鉴频器的处理流程图，包括下列步骤(下述步骤中n代表上述预定的数据的个数)：

对输入的n个采样数据进行fft变换；

对fft变换输出进行模平方运算；

从模平方后的n个数据中搜索最大值，并得到最大值所对应的索引idx；

当且仅当索引idx>n/2时，idx＝idx–n；

计算粗鉴频的频率误差df1＝idx/(tcoh*n)。

粗鉴频的频率误差df1送入超前滞后非相干鉴频器。利用df1，非相干鉴频器可以对输入信号中的频率误差进行初略的校正，以便进行细鉴频。

超前滞后非相干鉴频器的处理流程如图3所示，包括下列步骤：

超前支路进行频率搬移，得到搬移后的信号值z1，如公式(1)：

zl＝y*exp(-j2pi(df1–fsh))(1)

式中，pi为圆周率，y为相关器输出的复数信号，y＝yi+jyq，fsh为超前滞后非相干鉴频器的频率鉴别间隔，其必须满足条件fsh≤1/(tcoh*n)。间隔越小，鉴频输出的线性度越好，优选地，fsh可选择为1/(tcoh*n*8)。

滞后支路进行频率搬移，得到搬移后的信号值zr，如公式(2)，表示滞后支路进行频率搬移后的信号值：

zr＝y*exp(-j2pi(df1+fsh))(2)

对频率搬移后的信号进行相干累加，如公式(3)：

xl＝∑zl(i),i＝1,2,…n

xr＝∑zr(i),i＝1,2,…n(3)

公式(3)中，x1表示总的超前支路进行频率搬移后的信号值之和，xr表示总的滞后支路进展频率搬移后的信号值之和。

对上述进行相干累加后的信号再进行模平方，如公式(4)：

el＝|xl|²

er＝|xr|²(4)

最后，进行频率误差计算，得到细鉴频的频率误差df2，频率误差计算公式为：

fft鉴频得到的粗鉴频信号df1和非相干鉴频器得到的细鉴频信号df2相加，得到总的载波频率误差信号。频率误差信号送入环路滤波器进行滤波后，驱动数控振荡器(nco，numericallycontrolledoscillator)，完成载波频率跟踪。

如图4所示的频率跟踪环路的多普勒频率仿真结果示意图，图中的1代表仿真参考，该仿真曲线在10s后没有加入噪声，故多普勒频率非常稳定；2代表本发明实施例的仿真结果；3代表采用fft鉴频算法的仿真结果。从图中可以看出，本发明的频率误差补偿的超前滞后非相干鉴频比传统的fft鉴频更为稳定。

本发明还提供一种计算机可读存储介质，其存储有与计算设备结合使用的计算机程序，该计算机程序可被处理器执行以实现上述任一种所述方法。

本领域技术人员将认识到，对以上描述做出众多变通是可能的，所以实施例和附图仅是用来描述一个或多个特定实施方式。

尽管已经描述和叙述了被看作本发明的示范实施例，本领域技术人员将会明白，可以对其做出各种改变和替换，而不会脱离本发明的精神。另外，可以做出许多修改以将特定情况适配到本发明的教义，而不会脱离在此描述的本发明中心概念。所以，本发明不受限于在此披露的特定实施例，但本发明可能还包括属于本发明范围的所有实施例及其等同物。