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

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

背景技术：

卫星导航系统已经成为当今世界国家综合国力及科学技术发展水平的重要标志，是国家经济基础的重要组成部分，它已经进入人们的日常生活中，与社会发展和经济建设息息相关。我国也已经投入大量的资金和人力，积极进行北斗卫星导航系统(beidounavigationsatellitesystem，bds)的研究。到2020年左右，我国将建成覆盖全球的北斗卫星导航系统，其将成为继美国gps、俄罗斯glonass和欧洲galileo之后的又一个全球全天候的卫星导航系统。随着卫星导航应用的进一步发展，要求接收机能工作于更为严苛的应用环境下，例如，城市峡谷、林荫大道等卫星信号严重遮挡环境。所谓的高灵敏度接收机(hs-gnss,highsensitivity-globalnavigationsatellitesystems),就是为满足此类对弱信号接收的需求。载波跟踪环路是弱信号接收的关键环节，其性能对高灵敏度接收机具有重要影响。通常，通过加长相干积分时间，来抵抗弱信号中噪声的影响，以提高载波跟踪环路性能。然而，由于gps以及bds等的电文数据比特调制速率为50bps，意味着相干积分时间不能超过一个电文比特长度(20ms)。为了突破对数据调制的限制，需要对接收信号中的电文比特进行剥离，而这要求电文比特是已知的。通过移动通信网络等辅助手段，来获得电文数据；但此类方案会增加接收机的复杂度以及成本，丧失接收机的独立性，并依赖移动通信网络等特定环境。

作为一种不依赖辅助网络的独立的高灵敏度接收机，asoloviev等人在“decodingnavigationdatamessagesfromweakgpssignals,”ieeetrans.aerospaceandelectronicsystems,vol.45,no.2,pp.660–665,april2009，提出了一种采用最大似然估计(maximumlikelihood，ml)的算法，来自预测电文比特。然而，ml估计算法需要2^n-1次相关运算，n为相干积分时间所对应的电文比特数。当相干积分时间较大时，例如n＝16，需要32768次相关运算，导致难以实现。yan,k.等人在"weakgpssignaltrackingusingfftdiscriminatorinopenloopreceiver",gpssolut(2016)20:225，提出了一种对电文数据进行平方的算法，来规避电文数据调制的影响。平方算法会带来平方损耗，导致环路性能下降。

技术实现要素：

本发明的目的是为了解决现有技术中基于上述情况，本发明提供了一种北斗导航弱信号的载波频率跟踪方法。通过差分解调，得到电文比特的相位翻转，从而预测电文比特，解除了相干积分时间的限制，可以利用长相干积分时间，改善载波跟踪环路的稳定性，提高灵敏度性能。

为解决上述技术问题，本发明提出一种北斗导航弱信号的载波频率跟踪方法，包括下列步骤：

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

s2、对累加运算结果，进行差分解调电文自预测，估计出发送的电文数据比特dn；

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

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

s5、频谱分析输出信号送入鉴频器，得到载波频率误差的估计值；

s6、频率误差的估计值送入环路滤波器进行滤波；

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

与现有技术相比，本发明的有益效果有：

根据本发明提供的方法，利用低复杂度的差分解调算法进行电文自预测，解除了电文数据调制对相干积分时间的限制，可以利用长相干积分时间，本发明实施例的载波频率跟踪环路更为稳定可靠，改善了弱信号接收性能。

附图说明

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

图2是本发明实施例的差分解调自预测电文的流程图。

图3是本发明实施例的鉴频器框图。

图4是载波跟踪环路的频率误差仿真结果示意图。

具体实施方式

本方法的创新原理是利用低复杂度的差分解调算法进行电文自预测，解除了电文数据调制对相干积分时间的限制，可以利用长相干积分时间。下面结合附图和实施方式对本发明作进一步说明。

具体实施方式1

图1所示为本发明实施例的的载波跟踪环路框图。接收信号经过伪码解扩以及载波剥离后，送入相关器，相关积分时间为一个电文比特。相关器输出可以表示为(参考公式(1))：

zn(δfd)＝bnexp{j2πδfdntcoh+jδθ0)}sinc(πδfdtcoh)

+ηn,(n＝1,2,…n)(1)

式中bn为发送电文比特，δfd为跟踪环路的残留载波频率误差，tcoh为相关积分时间(20ms)，δθ0为载波相位误差，ηn为背景噪声。sinc函数表示载波剥离后的频谱。

相关器输出的相关值，送入差分解调电文自预测，估计出发送的电文数据比特dn。使用自预测电文比特dn对接收信号进行电文剥离，以去除电文数据调制对后续的相干积分的影响。

电文剥离后的数据送入缓存器，缓存到预定的n个数据后，送入fft进行频谱分析。通过fft变换，卫星信号频率点的功率，比噪声功率增强了n倍，即fft频谱分析过程亦为一个相干累加过程。因此，跟踪环路总的相干积分时间t＝tcoh*n。优选地，n＝32，相干积分时间为640ms。此时，频谱分析的分辨率为1/t＝1.5625hz。

fft输出信号送入鉴频器，得到载波频率误差的估计值。频率误差信号送入环路滤波器进行滤波后，驱动数控振荡器(nco，numericallycontrolledoscillator)，完成了载波频率跟踪。

图2是本发明实施例的差分解调电文自预测的流程图。首先，对相关器输出的20ms积分的相关值，进行相邻数据的差分检测，得到(参考公式(2))：

更进一步，利用差分检测结果，求取相邻相关值(即相邻电文数据)的相位差(参考公式(3))：

θn(δfd)＝tan^-1(yn(δfd))(3)

式中tan^-1为反正切函数。

gps/bds的电文调制为二进制相位调制(bpsk)，任何两个发送调制信号之间的相位差必定为“0”，表示无相位差；或者“π”，表示相位反相。通过门限π/2的比较，得到相邻比特是否相位反相的判决(参考公式(4)):

当trn＝1时，表示第n比特与第n-1比特相比，相位反相。假设固定第一个比特d1＝“+1”或者“-1”，剩余的其他n-1个比特可以估计为(参考公式(5)):

dn＝dn-1·(-trn),n＝2,3,…,n(5)

由于环路中采用了非相干解调，对于dn,n＝1,2,3,…,n的绝对相位并不敏感，而只关心它们的相对一致性。因此，d1可以任意地固定为“+1”或者“-1”，不会影响环路性能。

图3是本发明实施例的鉴频器框图。首先求取fft输出信号的模平方(功率值)，再通过最大值搜索得到功率峰值。功率峰值即代表卫星信号所在频点的强度，而峰值频点的索引值表示了频率误差。当环路中的无载波频率误差时，峰值频点的索引落于零点。fft输出功率值进行平均，得到功率平均值。用功率平均值对峰值功率进行归一化，得到归一化的峰值功率，并与预定门限进行比较。当归一化峰值功率大于预定门限，优选地，门限取值2.0，认为该鉴频输出可靠，输出到环路滤波器；否则的话，鉴频器输出置“0”，表示当前的鉴频值无效。通过上述的鉴频输出控制，提高了鉴频器的可靠性，从而改善了频率跟踪环路的可靠性和稳定性。

图4是载波跟踪环路的频率误差仿真结果示意图。仿真中，卫星信号初始的解调载噪比为48db-hz，在10秒后加入31dbc的噪声，解调载噪比约为17db-hz；相干积分时间t＝640ms；卫星信号的多普勒加速度为0.5hz/s。图4中，“差分解调”为本发明实施例的仿真结果，”平方法”为平方算法的仿真结果。本发明的载波频率跟踪环路的残留频率误差波动小，取得了良好效果。

根据本发明提供的方法，利用低复杂度的差分解调算法进行电文自预测，解除了电文数据调制对相干积分时间的限制，可以利用长相干积分时间，本发明实施例的载波频率跟踪环路更为稳定可靠，改善了弱信号接收性能。