CN104614740B - 一种导航信号数据导频联合跟踪方法及装置 - Google Patents

Abstract

本发明公开了一种导航信号数据导频联合跟踪方法及装置，包括：数据导频联合中频信号与载波数字控制振荡器(NCO)控制的本地载波相乘，完成载波剥离；载波剥离后的信号分别与码NCO控制的数据基带信号和导频基带信号相乘，完成码剥离；码剥离后的信号通过积分和清零处理，得到各支路相干积分结果；利用相干积分结果实现数据翻转检测和概率加权因子计算；数据翻转检测结果和概率加权因子辅助数据导频联合载波调整量估计，得到载波调整量；数据翻转检测结果和概率加权因子辅助数据导频联合码调整量估计，得到码调整量；载波调整量控制载波NCO，实现对数据导频联合载波信号的跟踪；码调整量控制码NCO，实现对数据导频联合基带信号的跟踪。

Description

一种导航信号数据导频联合跟踪方法及装置

技术领域

本发明属于导航信号跟踪技术领域，更具体地，涉及一种导航信号数据导频联合跟踪方法及装置。

背景技术

随着全球导航卫星系统(GNSS)的持续建设，导航信号体制引入了许多新的技术。与传统单数据通道的导航信号相比，导频通道的引入是现代化的导航信号体制的最大特点。对GPS和Galileo的ICD文档进行分析后可以发现，大多数新设计的导航信号均引入了导频分量，其中包括：GPS L1C、GPS L2C、GPS L5C、Galileo E1 OS、Galileo E5、Galileo E6。

导频通道的引入极大的提升了导航信号的跟踪性能。由于导频信号没有调制数据位，不存在数据位翻转的问题，所以在跟踪导频信号时，可以使用更长的相干积分时间，从而提高跟踪精度。除此之外，导频通道的载波跟踪可以使用对数据位翻转敏感的纯锁相环(Pure PLL)，而数据通道由于有数据位的翻转，其载波跟踪只能使用对数据位翻转不敏感的科斯塔斯环(Costas loop)。纯锁相环的动态牵引范围是360°，对应的相位噪声1-sigma经验门限为45°，而科斯塔斯环的动态牵引范围只有180°，对应的相位噪声1-sigma经验门限为15°，由此可见，纯锁相环的跟踪灵敏度要优于科斯塔斯环。但是这种单导频的跟踪方式却会带来另外一个问题：数据通道的信号没有参与跟踪，造成能量损失。所以在跟踪过程中如何将数据和导频联合成为现代化导航信号跟踪的研究热点。

对于导航信号数据导频联合跟踪，目前国内外已经有了许多研究成果。主要可以分为两类：一种是数据导频非相干累加，另一种是数据导频相干累加。在参考文献[1]“Trade-Off Between Data Rate and Signal PowerSplit in GNSS Signal Design”中提到了一种数据导频非相干联合跟踪的方法，该方法在数据通道采用科斯塔斯环，在导频通道采用纯锁相环，然后将两个环路的鉴别器输出按照一定权重相加，得到一个可供数据通道和导频通道共用的鉴别器输出。但是这种非相干的累加方法会带来平方损耗，使得跟踪精度并不是最优。除此之外，数据通道采用科斯塔斯环，其牵引范围小于导频通道的纯锁相环，这种联合方式相对单导频的跟踪灵敏度的提升并不明显。

在参考文献[2]“Dual Channel Optimizations of Tracking Schemesfor E1 CBOC Signal”中，提出了一种数据导频相干累加的方法，该方法通过对数据位的估计，遍历数据通道与导频通道的可能组合方式，然后通过一定的判决方式确定最终的组合方式，将组合后的相干积分结果用于跟踪。这种方法在高载噪比的时候，能达到最佳的跟踪性能，但是在低载噪比的时候性能不佳。因此，有必要为带导频通道的导航信号提供一种性能良好且结构简单的数据导频联合跟踪的方法。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供一种通用的、实现复杂度低、跟踪精度高、跟踪灵敏度高的导航信号数据导频联合跟踪方法及装置，用于带导频通道的导航信号的跟踪，并提供良好的跟踪性能。

为实现上述目的，按照本发明的一个方面，提供一种导航信号数据导频联合跟踪方法，包括以下步骤：

步骤1将数字中频信号与本地载波相乘，完成载波剥离，所述本地载波采用同向和正交两个支路；

步骤2载波剥离后的信号分别与数据基带信号和导频基带信号相乘，完成码剥离，基带信号采用超前、即时和滞后三个支路，其中所述基带信号包括所述数据基带信号和所述导频基带信号；

步骤3码剥离后的信号通过积分和清零处理，得到各个支路的相干积分结果，其中，导频通道的相干积分结果包括：同相超前I_PE、同相即时I_PP、同相滞后I_PL、正交超前Q_PE、正交即时Q_PP以及正交滞后Q_PL，数据通道的相干积分结果包括：同相超前I_DE、同相即时I_DP、同相滞后I_DL、正交超前Q_DE、正交即时Q_DP以及正交滞后Q_DL；

步骤4利用各即时支路的相干积分结果实现数据翻转检测，得到数据翻转检测输出Flip＝sign(I_PP×I_DP+Q_PP×Q_DP)，其中，sign表示取符号位函数，sign(x)表示当x大于等于0时输出为+1，当x小于零时输出为-1；

步骤5利用所述数据通道和所述导频通道的功率对比计算最优功率加权因子其中，P₁和P₂分别表示所述数据通道和所述导频通道信号功率；

步骤6利用各即时支路的相干积分结果计算概率加权因子α₂，所述概率加权因子α₂为与所述数据翻转检测正确概率单调性相同的权重因子；

步骤7计算得到同相即时I_P和正交即时Q_P相干积分结果如下所示：

I_P＝I_PP+Flip×I_DP×α₁ Q_P＝Q_PP+Flip×Q_DP×α₁

然后将所述同相即时I_P和所述正交即时Q_P相干积分结果进行鉴相，得到鉴相输出并将所述鉴相输出和所述概率加权因子α₂相乘，然后经环路滤波得到数据导频联合载波调整量，所述数据导频联合载波调整量反馈到控制所述本地载波的载波数字控制振荡器，对所述载波数字控制振荡器进行调整，实现对数据导频联合信号载波的跟踪；

步骤8计算得到数据导频联合同相超前I_E、同相滞后I_L、正交超前Q_E和正交滞后Q_L相干积分结果如下所示：

I_E＝I_PE+Flip×I_DE×α₁ I_L＝I_PL+Flip×I_DL×α₁

Q_E＝Q_PE+Flip×Q_DE×α₁ Q_L＝Q_PL+Flip×Q_DL×α₁

然后将所述同相超前I_E、所述同相滞后I_L、所述正交超前Q_E和所述正交滞后Q_L相干积分结果进行鉴相，得到鉴相结果Δτ，所述鉴相结果Δτ和概率加权因子α₂相乘，然后经环路滤波得到数据导频联合码调整量，所述码调整量反馈到控制所述基带信号的码数字控制振荡器，对所述码数字控制振荡器进行调整，实现对所述基带信号的跟踪。

按照本发明的另一方面，提供一种导航信号数据导频联合跟踪装置，包括：载波剥离器、码剥离器、本地载波生成器、基带信号生成器、载波数字控制振荡器、码数字控制振荡器、积分和清零器、数据翻转检测器、概率加权因子计算器、功率加权因子计算器、数据导频联合载波调整量估计器以及数据导频联合码调整量估计器，其中：

所述载波数字控制振荡器用于控制本地载波生成器生成本地载波；所述码数字控制振荡器用于控制基带信号生成器生成基带信号，所述基带信号包括数据基带信号和导频基带信号；所述载波剥离器用于实现所述本地载波与数字中频信号的载波剥离；所述码剥离器用于实现载波剥离后的信号与所述基带信号的码剥离；所述积分和清零器用于得到码剥离后的信号的相干积分结果；所述相干积分结果中的即时支路的相干积分结果送入所述数据翻转检测器得到数据翻转检测输出；所述即时支路的相干积分结果还送入所述概率加权因子计算器得到概率加权因子；所述功率加权因子计算器用于对比数据通道和导频通道的功率，得到功率加权因子；所述即时支路的相干积分结果还送入数据导频联合载波调整量估计器得到载波调整量，通过所述载波调整量控制所述载波数字控制振荡器，实现对数据导频联合信号载波的跟踪；所述相干积分结果中的超前支路和滞后支路的相干积分结果送入所述数据导频联合码调整量估计器得到码调整量，通过所述码调整量控制所述码数字控制振荡器，实现对所述基带信号的跟踪。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：

由于本发明数据翻转检测器的引入，解决了数据位翻转的问题，使得数据通道的相干积分结果和导频通道的相干结果可以进行相干累加，从而使得跟踪精度在高载噪比环境下能够达到最优。概率加权因子的引入，使得鉴别器输出结果按照最大似然概率加权，从而在低载噪比的环境下能够保证良好的跟踪性能，提高跟踪灵敏度。而且这两种创新性的结构引入都是在相干积分之后进行，属于低速处理，没有增加高速处理部分的复杂度，使得整个跟踪设备比较简单，实现复杂度低。

附图说明

图1为本发明导航信号数据导频联合跟踪方法的流程图；

图2为本发明数据导频联合跟踪方法的一实施例的整体实现框图；

图3为本发明数据翻转检测实现框图；

图4为本发明数据导频联合载波调整量估计实现框图；

图5为本发明数据导频联合码调整量估计实现框图；

图6为本发明数据导频联合跟踪装置的整体框图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

图1所示为本发明导航信号数据导频联合跟踪方法的流程图，数字中频信号与载波数字控制振荡器(numerically controlled oscillator，以下简称NCO)控制的本地载波相乘，完成载波剥离；载波剥离后的信号分别与码NCO控制的数据基带信号和导频基带信号相乘，完成码剥离；码剥离后的信号通过积分和清零处理，得到各个支路相干积分结果；利用该相干积分结果实现数据翻转检测和概率加权因子计算；利用数据通道和导频通道的功率对比计算功率加权因子；数据翻转检测结果、功率加权因子和概率加权因子辅助数据导频联合载波调整量估计，得到载波调整量；数据翻转检测结果、功率加权因子和概率加权因子辅助数据导频联合码调整量估计，得到码调整量；载波调整量控制载波NCO，实现对数据导频联合信号载波信号的跟踪；码调整量控制码NCO，实现对基带信号的跟踪。

为了描述方便，下文给出了一种典型实施例，但是该实施例仅仅是为了说明的目的而举的示例，而不是用来进行限制。本领域技术人员应当理解，凡在本申请的教导和权利要求保护范围下所作的任何修改、等同替换等，均应包含在本申请要求保护的范围内。

图2所示为本发明数据导频联合跟踪方法的一实施例的整体实现框图。本发明实施例的本地载波采用同向(I)和正交(Q)两个支路，基带信号采用超前(E)、即时(P)、滞后(L)三个支路。本发明实施例的跟踪起点为数字中频信号，对射频模拟信号的AD转换和下变频不予考虑。针对的信号体制是带导频通道的导航信号，接收的数字中频信号数学表达式如下公式(1)所示：

其中：P₁表示数据通道功率；d(t)表示当前时刻t的数据比特符号；c_d表示数据通道扩频码；τ表示信号时延；ω_IF表示信号角频率；表示载波初始相位；表示数据基带信号和导频基带信号之间的相位差；P₂表示导频通道功率；c_p表示导频通道扩频码；n(t)表示均值为0功率谱密度为N₀的白噪声。这里给出接收信号的数学表达式并不是为了限定本发明，凡是含有数据通道和导频通道并且二者之间相位差恒定的信号都不脱离本发明的申请范围。

如图2所示的典型实施例的详细步骤如下：

步骤1：由于是数据导频联合跟踪，数据通道和导频通道信号应该同频并且相位差恒定，所以数据通道和导频通道共用载波NCO 1和码NCO 6。如图2所示，载波NCO 1分别控制余弦映射表2和正弦映射表3。余弦映射表2产生同相载波信号，正弦映射表3产生正交载波信号。同相载波信号和数字中频信号在采样脉冲控制下通过乘法器4相乘得到同相支路信号I，正交载波信号和数字中频信号在采样脉冲控制下通过乘法器5相乘得到正交支路信号Q。

步骤2：导频通道和数据通道共用的码NCO 6分别驱动导频基带信号发生器7和数据基带信号发生器8生成超前的导频基带信号和超前的数据基带信号。超前的导频基带信号在延时器9的作用下产生即时导频基带信号和滞后导频基带信号。超前的数据基带信号在延时器10的作用下产生即时的数据基带信号和滞后的数据基带信号。

步骤3：导频通道超前、即时、滞后三路基带信号分别与相位移动后的同相支路信号I通过乘法器11、12、13相乘，再通过积分和清零器23、24、25累加，得到导频同相超前、即时和滞后相干积分结果(即I_PE、I_PP、I_PL)；

导频通道超前、即时、滞后三路基带信号分别与相位移动后的正交支路信号Q通过乘法器14、15、16相乘，再通过积分和清零器26、27、28累加，得到导频正交超前、即时和滞后相干积分结果(即Q_PE、Q_PP、Q_PL)；

数据通道超前、即时、滞后三路基带信号分别与同相支路信号I通过乘法器20、21、22相乘，再通过积分和清零器34、33、32累加，得到数据同相超前、即时和滞后相干积分结果(即I_DE、I_DP、I_DL)；

数据通道超前、即时、滞后三路基带信号分别与正交支路信号Q通过乘法器17、18、19相乘，再通过积分和清零器31、30、29累加，得到数据正交超前、即时和滞后相干积分结果(即Q_DE、Q_DP、Q_DL)。

步骤4：利用上述即时支路相干积分结果I_PP、Q_PP、I_DP、Q_DP进行数据翻转检测，数据翻转检测的实现框图如图3所示。I_PP和I_DP通过乘法器35相乘，Q_PP和Q_DP通过乘法器36相乘，两个相乘结果通过加法器37相加，再由符号提取器38对相加的结果取符号，得到数据翻转检测输出Flip：

Flip＝sign(I_PP×I_DP+Q_PP×Q_DP)

其中，sign表示取符号位函数，sign(x)表示当x大于等于0时输出为+1，当x小于零时输出为-1，即数据翻转检测输出为正负1，输出Flip为正1时表示数据通道与导频通道同相，为负1时表示数据通道与导频通道相差180°。

步骤5：功率加权因子的计算。由于实际的信号体制中，数据导频功率比并不是1:1，所以将数据导频相干积分结果进行相干累加时需要考虑功率的影响。功率加权因子α₁会直接影响数据导频联合跟踪精度，本发明实施例给出最优化的功率加权因子，如下公式(2)所示：

${\mathrm{\α}}_1=\sqrt{\frac{P_1}{P_2}}---\left(2\right)$

其中，P₁和P₂分别表示数据通道和导频通道信号功率。选取该加权因子，数据导频联合的跟踪精度能够达到最优。

步骤6：概率加权因子计算。对于数据翻转的检测，由于有噪声的存在，所以检测结果并不是百分之百正确，而是以一定的概率正确。对于正确概率高的结果，可以认为本次鉴别器的输出可靠的概率高，所以应该以较大的权重输出；相反，对于正确概率低的结果，可以认为本次鉴别器的输出可靠的概率低，所以应该以较小的权重输出。

根据上述公式(1)，数据通道和导频通道即时支路的相干积结果的数学表达式如下所示：

其中，T_p表示相干积分时间；d表示当前数据位；n_DIP、n_DQP、n_PIP、n_PQP分别表示各支路相干积分后归一化的噪声。

定义θ为矢量1(I_DP，Q_DP)和矢量2(I_PP，Q_PP)的夹角。没有信号存在时，由于n_DIP、n_DQP、n_PIP、n_PQP是零均值的白噪声，所以θ应该在0～180°内均匀分布。有信号存在时，假设d＝1，则两个矢量信号部分应该是同相，则θ角度越小的概率越高。根据前面数据翻转检测器的原理，当θ<90°时，判断d＝1，此时θ角度越小，判断正确的概率越高，而θ角度越接近90°，判断正确的概率越低；假设d＝-1，则两个矢量信号部分应该是反相，则θ角度越大的概率越高。根据前面数据翻转检测器的原理，当θ>90°时，判断d＝-1，此时θ角度越接近180°，判断正确的概率越高，而θ角度越接近90°，判断正确的概率越低。实际上在θ＝90°时，可以认为由于有噪声的存在，此时已经完全无法判断d的正负，这种结果完全不可靠，应以权重0输出。但是实际的计算这种最大似然概率权重的过程非常复杂，不利于工程实现，凡是与数据翻转检测正确概率单调性相同的权重因子都可以作为近似的概率加权因子，近似度越高，跟踪结果越好。在本发明实施例中用一个和最大似然概率权重单调性相同且近似度较高的权重因子来替换。定义概率加权因子α₂如下公式(7)所示：

${\mathrm{\&alpha;}}_2=\frac{Q_{\mathrm{PP}}\&times;Q_{\mathrm{DP}}+I_{\mathrm{PP}}\&times;I_{\mathrm{DP}}}{\sqrt{({I_{\mathrm{DP}}}^2+{Q_{\mathrm{DP}}}^2)({I_{\mathrm{PP}}}^2+{Q_{\mathrm{PP}}}^2)}}---\left(7\right)$

这种基于最大似然概率的加权因子并不唯一，还有许多近似的方法，任何与数据翻转检测正确概率单调性相同的权重因子均不脱离本申请的范围。

步骤7：数据导频联合载波调整量估计实现框图如图4所示。即时支路相干积分结果(即I_PP、Q_PP、I_DP、Q_DP)用来进行数据导频联合载波调整量估计。I_DP和Q_DP分别通过乘法器39、40与数据翻转检测器输出Flip相乘，消除数据翻转的影响。相乘的结果再分别通过乘法器41、42与功率加权因子α₁相乘，其结果再分别通过加法器43、44与I_PP和I_DP进行相干累加得到I_P和Q_P。I_P和Q_P再通过鉴相器45进行鉴相，得到鉴相输出由于消除了数据位翻转，可以使用纯锁相环鉴相器，例如：四象限反正切arctan(I_P，Q_P)。鉴相输出通过乘法器46与概率加权因子α₂相乘，再通过环路滤波器47进行滤波处理，得到最终的载波调整量。载波调整量反馈到图2中的载波NCO 1，对载波NCO 1进行调整，实现对数据导频联合信号载波的跟踪。

步骤8：数据导频联合码调整量估计实现框图如图5所示。超前支路和滞后支路的相干积分结果(即I_PE、Q_PE、I_PL、Q_PL、I_DE、Q_DE、I_DL、Q_DL)用来进行数据导频联合码调整量估计。I_DE、Q_DE、I_DL、Q_DL分别通过乘法器49、50、51、52与数据翻转检测器输出Flip相乘，消除数据翻转的影响。相乘的结果再分别通过乘法器53、54、55、56与功率加权因子α₁相乘，其结果再分别通过加法器57、58、59、60与I_PE、Q_PE、I_PL、Q_PL行相干累加得到联合相干积分结果I_E、Q_E、I_L、Q_L。码相位鉴别器有许多不同的实现方法，为了保持发明的完整性，给出一种非相干超前减滞后鉴别器，其它的实现方式不脱离本申请的范围。联合相干积分结果I_E、Q_E和I_L、Q_L分别通过非相干累加器61、62进行非相干累加得到超前非相干累加结果E_S和滞后非相干累加结果L_S，其数学表达式如下公式(8)所示：

$\begin{array}{cc}{E}_S=\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}(I_E^2\left(k\right)+Q_E^2\left(k\right))& L_S=\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}(I_L^2\left(k\right)+Q_L^2\left(k\right))---\left(8\right)\end{array}$

其中，N表示非相干累加次数。超前非相干累加结果E_S和滞后非相干累加结果L_S再通过鉴相器63进行鉴相，得到鉴相结果Δτ。这里选取非相干超期减滞后平方鉴相器，并做归一化处理，采用该鉴相器只是为了保持本发明的完整性，其它的实现方式不脱离本申请的范围。其数学表达式如下公式(9)所示：

$\mathrm{\&Delta;\&tau;}=\frac{E_S-L_S}{E_S+L_S}---\left(9\right)$

鉴相结果Δτ通过乘法器64与概率加权因子α₂相乘，再通过环路滤波器65进行滤波处理，得到最终的码调整量。码调整量反馈到图2中的码NCO 6，对码NCO 6进行调整，实现对基带信号的跟踪。

图6所示为本发明数据导频联合跟踪装置的整体框图，载波NCO 104控制本地载波生成器102生成本地载波，码NCO 105控制基带信号生成器103生成基带信号，需要注意的是，基带信号生成器103生成的基带信号既包括数据基带信号，也包括导频基带信号。本地载波与数字中频信号通过载波剥离器100实现载波剥离，载波剥离后的信号与基带信号通过码剥离器101实现码剥离。码剥离后的信号再通过积分和清零器106得到相干积分结果，其中即时支路的相干积分结果送入数据翻转检测器107，得到数据翻转检测输出Flip，数据翻转检测器107的实现原理如图3所示。即时支路的相干积分结果还送入概率加权因子计算器108，得到概率加权因子α₂，其实现原理由公式(7)定义。利用数据通道和导频通道的功率对比，通过功率加权因子计算器109，得到功率加权因子α₁。即时支路相干积分结果还送入数据导频联合载波调整量估计器110，其一种典型的实施例如图4所示，数据导频联合载波调整量估计器110得到载波调整量，通过该载波调整量控制载波NCO 104，实现对数据导频联合信号载波的跟踪。积分和清零器106得到相干积分结果中的超前支路和滞后支路的相干积分结果送入数据导频联合码调整量估计器111，其一种典型的实施例如图5所示，数据导频联合码调整量估计器111得到码调整量，通过该码调整量控制码NCO105，实现对基带信号的跟踪。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种导航信号数据导频联合跟踪方法，其特征在于，包括：

步骤1将数字中频信号与本地载波相乘，完成载波剥离，所述本地载波采用同相和正交两个支路；

步骤2载波剥离后的信号分别与数据基带信号和导频基带信号相乘，完成码剥离，基带信号采用超前、即时和滞后三个支路，其中所述基带信号包括所述数据基带信号和所述导频基带信号；

步骤3码剥离后的信号通过积分和清零处理，得到各个支路的相干积分结果，其中，导频通道的相干积分结果包括：同相超前I_PE、同相即时I_PP、同相滞后I_PL、正交超前Q_PE、正交即时Q_PP以及正交滞后Q_PL，数据通道的相干积分结果包括：同相超前I_DE、同相即时I_DP、同相滞后I_DL、正交超前Q_DE、正交即时Q_DP以及正交滞后Q_DL；

步骤4利用各即时支路的相干积分结果实现数据翻转检测，得到数据翻转检测输出Flip＝sign(I_PP×I_DP+Q_PP×Q_DP)，其中，sign表示取符号位函数，sign(x)表示当x大于等于0时输出为+1，当x小于零时输出为-1；

步骤5利用所述数据通道和所述导频通道的功率对比计算最优功率加权因子其中，P₁和P₂分别表示所述数据通道和所述导频通道信号功率；

步骤6利用各即时支路的相干积分结果计算概率加权因子α₂，所述概率加权因子α₂为与所述数据翻转检测正确概率单调性相同的权重因子；

步骤7计算得到同相即时I_P和正交即时Q_P相干积分结果如下所示：

I_P＝I_PP+Flip×I_DP×α₁ Q_P＝Q_PP+Flip×Q_DP×α₁

然后将所述同相即时I_P和所述正交即时Q_P相干积分结果进行鉴相，得到鉴相输出并将所述鉴相输出和所述概率加权因子α₂相乘，然后经环路滤波得到数据导频联合载波调整量，所述数据导频联合载波调整量反馈到控制所述本地载波的载波数字控制振荡器，对所述载波数字控制振荡器进行调整，实现对数据导频联合信号载波的跟踪；

步骤8计算得到数据导频联合同相超前I_E、同相滞后I_L、正交超前Q_E和正交滞后Q_L相干积分结果如下所示：

I_E＝I_PE+Flip×I_DE×α₁ I_L＝I_PL+Flip×I_DL×α₁

Q_E＝Q_PE+Flip×Q_DE×α₁ Q_L＝Q_PL+Flip×Q_DL×α₁

然后将所述同相超前I_E、所述同相滞后I_L、所述正交超前Q_E和所述正交滞后Q_L相干积分结果进行鉴相，得到鉴相结果Δτ，所述鉴相结果Δτ和概率加权因子α₂相乘，然后经环路滤波得到数据导频联合码调整量，所述码调整量反馈到控制所述基带信号的码数字控制振荡器，对所述码数字控制振荡器进行调整，实现对所述基带信号的跟踪。

2.如权利要求1所述的方法，其特征在于，所述步骤1中接收的所述数字中频信号含有所述数据通道和所述导频通道并且二者之间相位差恒定。

3.如权利要求1所述的方法，其特征在于，所述步骤1中接收的所述数字中频信号的表达式如下所示：

其中，P₁表示数据通道功率；d(t)表示当前时刻t的数据比特符号；c_d表示数据通道扩频码；τ表示信号时延；ω_IF表示信号角频率；表示载波初始相位；表示所述数据基带信号和所述导频基带信号之间的相位差；P₂表示导频通道功率；c_p表示导频通道扩频码；n(t)表示均值为0功率谱密度为N₀的白噪声。

4.如权利要求1-3中任一项所述的方法，其特征在于，所述步骤3中所述导频通道超前、即时、滞后三路基带信号分别与同相支路信号相乘，再分别进行累加，得到导频同相超前、即时和滞后相干积分结果；所述导频通道超前、即时、滞后三路基带信号还分别与正交支路信号相乘，再分别进行累加，得到导频正交超前、即时和滞后相干积分结果；

所述数据通道超前、即时、滞后三路基带信号分别与同相支路信号相乘，再分别进行累加，得到数据同相超前、即时和滞后相干积分结果；所述数据通道超前、即时、滞后三路基带信号还分别与正交支路信号相乘，再分别进行累加，得到数据正交超前、即时和滞后相干积分结果。

5.如权利要求1-3中任一项所述的方法，其特征在于，所述步骤4中所述数据翻转检测输出Flip为正1时表示所述数据通道与所述导频通道同相，为负1时表示所述数据通道与所述导频通道相差180°。

6.如权利要求1-3中任一项所述的方法，其特征在于，所述步骤6中所述概率加权因子为

7.如权利要求1所述的方法，其特征在于，所述步骤7中使用纯锁相环鉴相器对所述同相即时I_P和所述正交即时Q_P相干积分结果进行鉴相。

8.如权利要求1所述的方法，其特征在于，所述步骤8中联合相干积分结果I_E、Q_E和I_L、Q_L分别进行非相干累加得到超前非相干累加结果E_S和滞后非相干累加结果L_S，其数学表达式如下所示：

$\begin{array}{cc}{E}_S=\underset{k=1}{\overset{N}{\&Sigma;}}\left(I_E^2\right(k)+Q_E^2(k\left)\right)& L_S=\underset{k=1}{\overset{N}{\&Sigma;}}\left(I_L^2\right(k)+Q_L^2(k\left)\right)\end{array}$

其中，N表示非相干累加次数；所述超前非相干累加结果E_S和所述滞后非相干累加结果L_S再进行鉴相，得到所述鉴相结果Δτ。

9.一种导航信号数据导频联合跟踪装置，其特征在于，包括：载波剥离器、码剥离器、本地载波生成器、基带信号生成器、载波数字控制振荡器、码数字控制振荡器、积分和清零器、数据翻转检测器、概率加权因子计算器、功率加权因子计算器、数据导频联合载波调整量估计器以及数据导频联合码调整量估计器，其中：

所述载波数字控制振荡器用于控制本地载波生成器生成本地载波；所述码数字控制振荡器用于控制基带信号生成器生成基带信号，所述基带信号包括数据基带信号和导频基带信号；所述载波剥离器用于实现所述本地载波与数字中频信号的载波剥离；所述码剥离器用于实现载波剥离后的信号与所述基带信号的码剥离；所述积分和清零器用于得到码剥离后的信号的相干积分结果；所述相干积分结果中的即时支路的相干积分结果送入所述数据翻转检测器得到数据翻转检测输出；所述即时支路的相干积分结果还送入所述概率加权因子计算器得到概率加权因子；所述功率加权因子计算器用于对比数据通道和导频通道的功率，得到功率加权因子；所述即时支路的相干积分结果还送入数据导频联合载波调整量估计器得到载波调整量，通过所述载波调整量控制所述载波数字控制振荡器，实现对数据导频联合信号载波的跟踪；所述相干积分结果中的超前支路和滞后支路的相干积分结果送入所述数据导频联合码调整量估计器得到码调整量，通过所述码调整量控制所述码数字控制振荡器，实现对所述基带信号的跟踪。