CN106932795A - 一种gnss信号的矢量和标量混合跟踪方法及跟踪环路 - Google Patents

Abstract

本发明属于导航接收机设备研制领域，涉及一种GNSS信号的矢量和标量混合跟踪方法及跟踪环路。GNSS信号依次经过接收机中的天线、射频前端、AD转换器后转变为数字中频信号；本地载波生成装置NCO产生同相信号和正交信号；本地信号生成装置产生导频支路和数据支路本地伪码，分别与数字中频信号进行相乘混合；相关器分别对混合信号进行相关处理；鉴别器对相关值进行鉴别处理；对各个通道中鉴别器输出结果进行处理，获得精度更高的频率误差估计结果和精度值；接收机中的任意通道中的DUKF装置获得本地信号载波频率估计参数，输出至本地载波生成装置NCO，用于更新频率控制字。本发明具有较好的信号重捕获性能，降低了算法的计算复杂度。

Description

一种GNSS信号的矢量和标量混合跟踪方法及跟踪环路

技术领域

本发明属于导航接收机设备研制领域，具体涉及一种GNSS信号的矢量和标量混合跟踪方法及跟踪环路，其可运用在卫星导航***中接收终端类设备的研制中。

背景技术

随着卫星导航***发展，新一代导航信号将逐渐为用户提供服务，针对新一代导航信号中含有多路导航信号，如导频支路和数据支路，对不同支路的信号，其相干积分可以设计的不同，以提高信号的跟踪性能。目前GNSS信号的载波跟踪方法主要有两类，一类是基于标量跟踪环路(STL:Scalar Tracking Loop)的载波跟踪方法,另一类是基于矢量跟踪环路(VTL:Vector Tracking Loop)的载波跟踪方法。STL方法对每颗卫星进行独立跟踪，该方法计算复杂度低，易于接收机实现，VTL方法对所有可见卫星进行联合跟踪，对弱信号有较强的跟踪能力，可以快速实现信号的重捕获，提高可用性，但计算复杂度相对较高。

发明内容

针对新一代导航信号中导频和数据支路的特点，结合现有STL和VTL两种不同跟踪方法的特点，本发明提供了一种用于复合GNSS信号的矢量和标量混合跟踪环路(HTL：Hybrid Tracking Loop)及跟踪方法。通过利用双速率卡尔曼滤波器(DUKF：DualUpdate-rate Kalman Filter)将STL和VTL组合在一起构成HTL的信号载波跟踪方法。为实现上述目的，具体技术方案如下：

一种GNSS信号的矢量和标量混合跟踪方法，包括以下步骤：

步骤1，GNSS信号依次经过接收机中的天线，射频前端，AD转换器后变为数字中频信号r(t)；

步骤2，接收机有N个跟踪通道，每个跟踪通道中的处理方法相同，对于任意跟踪通道i中的本地载波生成装置NCO(Numerically Controlled Oscillator，缩写NCO)，其生成频率控制字为的两路信号，分别为同相载波信号和正交载波信号t表示时间，i＝1,2,…,N，N为正整数，具体为：

步骤3，接收机跟踪通道i中的本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，数据支路信号生成装置接收同相载波信号和正交载波信号分别与数据支路本地伪码c_d(t)相乘，产生数据支路的本地同相信号和正交信号导频支路信号生成装置接收同相载波信号和正交载波信号分别与导频支路的本地伪码c_p(t)相乘，产生导频支路的本地同相信号和正交信号信号称为本地复制信号，具体为：

步骤4，接收机跟踪通道i中的相关器进行相关处理，用于将本地复制信号 分别与数字中频信号r(t)相乘后的混合信号进行相干积累，设相干积分时间为T_c，对任意通道i，输出相关值为其中下标k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为T_c，故输出信号的积分区间为(k-1)·T_c到k·T_c，具体结果如下：

步骤5，接收机跟踪通道中的鉴别器装置对所述步骤4中输出的相关值进行处理，用于获得本地复制信号和数字中频信号间的误差估计参数，所述鉴别器包括数据支路鉴相器和导频支路鉴频器，经过鉴相器和鉴频器处理后，鉴相器输出误差估计参数为鉴频器输出误差估计参数为

其中，

其中atan表示反正切函数，atan2表示四象限反正切函数；N_p为导频支路相干积累次数；鉴相器的相干积分时间为T_c，鉴频器的相干积分时间为N_p·T_c，故鉴相器的结果每T_c时间有效一次，而鉴频器的结果每N_p·T_c时间输出一次结果。公式中I₁，I₂，Q₁，Q₂，m均表示计算过程中的中间量符号。

步骤6，接收机中的矢量频率跟踪环路(Vector Frequency Lock loop，简称VFLL)对各个通道中鉴频器输出结果进行处理，获得精度更高的频率误差估计结果其精度分别为

所述矢量频率跟踪环路处理过程包括步骤：

步骤61，根据各个通道中鉴频器的输出结果获取量测量Z_k及其噪声协方差矩阵R_z，其中N为接收到的卫星通道数；

VFLL的量测方程为：

其中为第k个跟踪历元接收机的运动状态，δv_x,δv_y,δv_z为ECEF坐标系下的三维速度误差，δa_x,δa_y,δa_z为ECEF坐标系下的三维加速度误差，δf为接收机上时钟的频率误差；H^V为测量矩阵，它由接收机到卫星的空间几何构型决定；为量测噪声，其协方差矩阵为其中diag()表示对角矩阵算子，为导频支路通道i中鉴频器输出结果的噪声方差，具体为：

其中Cⁱ/N₀表示通道i对应的信号载噪比，即信号功率与噪声的功率谱密度的比值；

步骤62，VFLL的迭代过程，具体描述如下：

***方程为

其中Φ^V为状态转移矩阵，具体表示为

其中

T_b＝N_p·T_c，为VFLL的更新间隔；

是统过程噪声，其协方差矩阵为Q^V，具体为：

其中

Q_f＝S_f·T_b

S_a为加速度噪声功率谱密度，S_f为时钟频率变化噪声功率谱密度。

根据步骤61中获取的量测信息，得到VFLL的滤波步骤如下：

Step 1,计算接收机状态矢量预测值及其协方差值

为对应的协方差矩阵；

Step 2，计算VFLL的增益矩阵

Step 3，更新接收机状态矢量及其协方差阵

其中I表示单位矩阵；

Step 4,计算各个通道频率估计误差

其估计精度满足

这样对于任意i通道而言，频率估计误差为估计精度为表示向量的第i个元素,表示矩阵的第i行第i列元素值

步骤7，接收机任意通道i中的DUKF装置，用于获得本地信号载波频率估计参数，将载波频率估计参数输入本地载波生成装置，更新频率控制字；

所述DUKF装置获得本地信号载波频率估计参数的步骤为：

步骤71，根据步骤5中鉴别器的输出结果和步骤6中VFLL的输出结果获取DUKF新息增量量测矩阵量测噪声矩阵当两种鉴别器的结果均有效时，计算公式如下：

其中为步骤5中鉴相器的输出结果,为步骤62输出的频率误差估计结果，H_d和H_p分别为两种不同量测量对应的量测矩阵，具体为

H_p＝[0 1 -(N_p -2)·T_c/2]

为数据支路鉴相器输出结果的噪声方差，具体为

为步骤62输出的频率误差估计精度。

当只有数据支路鉴相器有效时，只取鉴相器对应的项即可。

步骤72，对于任意通道i而言，其DUKF的迭代过程，具体描述如下：

DUKF的***方程为

其中为第k个跟踪历元通道i***状态向量，分别为信号的载波相位，多普勒频率和多普勒频率变化率，单位分别为周，Hz，Hz/s；w_k＝[ω_rf·w_b；ω_rf·w_d；(ω_rf/c)·w_a]^T为***噪声,w_b和w_d分别为由接收机中晶体振荡器的引起的相位噪声和频率噪声,其噪声谱密度分别为q_b和q_d；w_a是***频率变化率噪声，其功率谱密度为q_a。ω_rf表示载波频率，c为光速；Φ是***状态转移矩阵，具体为

w_k为***过程噪声，Q是w_k对应的过程噪声协方差矩阵，具体为

E[·]表示求均值符号；

结合步骤71中获取的新息信息，DUKF的滤波过程可以描述为

Step 1：计算***状态向量预测值

为第k-1个跟踪历元时刻的通道i***状态向量；

Step 2：计算***状态向量预测值的协方差矩阵

为的协方差矩阵；

Step 3：根据VFLL是否有结果输出获得量测信息若只有数据支路鉴相器有输出时，只取鉴相器对应项的数值；

Step 4：计算DUKF的增益矩阵

Step 4：根据新息更新状态估计结果：

Step 5：更新状态估计协方差矩阵：

步骤73，根据状态估计结果获得载波NCO的频率控制字即

其中，表示向量的第2个元素，至此，完成了一次DUKF的滤波处理过程。

本发明还提供了一种GNSS信号的矢量和标量混合跟踪环路，包括N个跟踪通道模块1和1个矢量频率跟踪环路2；所述N个跟踪通道模块具有相同的结构，包括本地载波生成装置11、本地信号生成装置12、第一乘法器13、第二乘法器14、第一相关器15、第二相关器16、鉴别器17，18和DUKF装置19；所述本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，用于产生导频支路和数据支路的伪码信号，并生成本地复制信号；所述鉴别器包括鉴相器17和鉴频器18，用于获取本地复制信号和接收信号之间的误差估计参数；

所述本地载波生成装置11根据输入的频率控制字，产生同相载波信号和正交载波信号；所述数据支路信号生成装置和导频支路信号生成装置的输入端分别与本地载波生成装置的输出端相连；

所述数据支路信号生成装置的输出端连接第一乘法器13的输入端，并将混合信号输出至第一相关器15的输入端；

所述导频支路信号生成装置的输出端连接第二乘法器14的输入端，并将混合信号输出至第二相关器16的输入端；

所述第一相关器15的输出端连接鉴相器17的输入端；鉴相器的输出端连接DUKF转置19的输入端；所述DUKF装置的输出端连接本地载波生成装置11；

所述第二相关器16的输出端连接鉴频器18的输入端；鉴频器的输出端连接所述矢量频率跟踪环路2的输入端；

所述矢量频率跟踪环路2的输出端分别输出至每个跟踪通道模块中的DUKF装置19的输入端。

采用本发明获得的有益技术效果：本发明通过利用DUKF滤波器，将不同更新速率下的标量跟踪环路和矢量跟踪环路了组合在一起，构成混合跟踪环路，对复合型GNSS信号进行联合跟踪。相比单独的标量跟踪环路而言，本发明混合跟踪环路具有较好的信号重捕获性能，相比单独的矢量跟踪环路而言，通过降低混合跟踪环路中矢量跟踪滤波器的更新频度，可以降低算法的计算复杂度。

附图说明

图1本发明方法流程示意图；

图2为VFLL装置的处理过程示意图；

图3为双速率卡尔曼滤波器(DUKF)的滤波过程示意图；

图4为本发明跟踪环路结构示意图；

图5为实施例中某场景下的GPS卫星星空图；

图6为实施例中某场景下本发明与现有技术对信号的跟踪结果对比图。

具体实施方式

下面结合附图和实施例对本发明作进一步说明。

如图1所示，为本发明流程图。本发明实施例给出了一种GNSS信号的矢量和标量混合跟踪方法，包括以下步骤：

步骤1，GNSS信号依次经过接收机中的天线，射频前端，AD转换器后变为数字中频信号r(t)；

步骤2，接收机有N个跟踪通道，每个跟踪通道中的处理方法相同，对于任意跟踪通道i中的本地载波生成装置NCO(Numerically Controlled Oscillator，缩写NCO)，其生成频率控制字为的两路信号，分别为同相载波信号和正交载波信号t表示时间，具体为：

步骤3，接收机跟踪通道i中的本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，数据支路信号生成装置接收同相载波信号和正交载波信号分别与数据支路本地伪码c_d(t)相乘，产生数据支路的本地同相信号和正交信号导频支路信号生成装置接收同相载波信号和正交载波信号分别与导频支路的本地伪码c_p(t)相乘，产生导频支路的本地同相信号和正交信号信号称为本地复制信号，具体为：

步骤4，接收机跟踪通道i中的相关器进行相关处理，用于将本地复制信号 和数字中频信号r(t)进行相干积累，设相干积分时间为T_c，对任意通道i，输出相关值为其中下标k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为T_c，故输出信号的积分区间为(k-1)·T_c到k·T_c，具体结果如下：

步骤5，接收机跟踪通道中的鉴别器装置对所述步骤4中输出的相关值进行处理，用于获得本地复制信号和数字中频信号间的误差估计参数，所述鉴别器包括数据支路鉴相器和导频支路鉴频器，经过鉴相器和鉴频器处理后，鉴相器输出误差估计参数为鉴频器输出误差估计参数为

其中

其中atan表示反正切函数，atan2表示四象限反正切函数；N_p为相干积累次数；鉴相器的相干积分时间为T_c，N_p为相干积累次数，鉴频器的相干积分时间为N_p·T_c，故鉴相器的结果每T_c时间有效一次，而鉴频器的结果每N_p·T_c时间输出一次结果。

步骤6，接收机中的矢量频率跟踪环路(Vector Frequency Lock loop，简称VFLL)对各个通道中鉴频器输出结果进行处理，获得精度更高的频率误差估计结果其精度分别为

图2给出了VFLL装置的处理过程示意图，其具体步骤包括：

步骤61，根据各个通道中鉴频器的输出结果获取量测量Z_k及其噪声协方差矩阵R_z，其中N为接收到的卫星通道数；

VFLL的量测方程为：

其中为第k跟踪历元接收机的运动状态，δv_x,δv_y,δv_z为ECEF坐标系(Earth-Centered,Earth-Fixed,缩写ECEF)下的三维速度误差，δa_x,δa_y,δa_z为ECEF坐标系下的三维加速度误差，δf为接收机上时钟的频率误差；H^V为测量矩阵，它由接收机到卫星的空间几何构型决定；为量测噪声，其协方差矩阵为 其中diag()表示对角矩阵算子，为导频支路通道i中鉴频器输出结果的噪声方差，具体为

其中Cⁱ/N₀表示通道i对应的信号载噪比，即信号功率与噪声的功率谱密度的比值；

步骤62，VFLL的迭代过程，具体描述如下：

***方程为

其中Φ^V为状态转移矩阵，具体表示为

其中

T_b＝N_p·T_c，为VFLL的更新间隔；

是***过程噪声，其协方差矩阵为Q^V，具体为：

其中

Q_f＝S_f·T_b

S_a为加速度噪声功率谱密度，S_f为时钟频率变化噪声功率谱密度。

根据步骤61中获取的量测信息，得到VFLL的滤波步骤如下：

Step 1,计算接收机状态矢量预测值及其协方差值

为对应的协方差矩阵；

Step 2，计算VFLL的增益矩阵

Step 3，更新接收机状态矢量及其协方差阵

其中I表示单位矩阵；

Step 4,计算各个通道频率估计误差

其估计精度满足

这样对于任意i通道而言，频率估计误差为估计精度为表示向量的第i个元素,表示矩阵的第i行第i列元素值

步骤7，接收机任意通道i中的DUKF装置，用于获得本地信号载波频率估计参数，将载波频率估计参数输入本地载波生成装置，更新频率控制字；

图3为DUKF滤波器的处理过程示意图，DUKF装置获得本地信号载波频率估计参数的步骤为：

步骤71，根据步骤5中鉴别器的输出结果和步骤6中VFLL的输出结果获取DUKF新息增量量测矩阵量测噪声矩阵当两种鉴别器的结果均有效时，计算公式为

其中为步骤5中鉴相器的输出结果,为步骤62输出的频率误差估计结果，H_d和H_p分别为两种不同量测量对应的量测矩阵，具体为

H_p＝[0 1 -(N_p -2)·T_c/2]

N_p为导频支路的相干积累次数；

为数据支路鉴相器输出结果的噪声方差，具体为

为步骤62输出的频率误差估计精度。

当只有数据支路鉴相器有效时，只取鉴相器对应的项即可。

步骤72，对于任意通道i而言，其DUKF的迭代过程，具体描述如下：

DUKF的***方程为

其中为第k跟踪历元通道i***状态向量，分别为信号的载波相位，多普勒频率和多普勒频率变化率，单位分别为周，Hz，Hz/s；w_k＝[ω_rf·w_b；ω_rf·w_d；(ω_rf/c)·w_a]^T为***噪声,w_b和w_d分别为由接收机中晶体振荡器的引起的相位噪声和频率噪声,其噪声谱密度分别为q_b和q_d；w_a是***频率变化率噪声，其功率谱密度为q_a。ω_rf表示载波频率；c≈3×10⁸m/s，为光速；Φ是***状态转移矩阵，具体为

w_k为DUKF***过程噪声，Q是w_k对应的过程噪声协方差矩阵，具体为

实施例中q_b和q_d通常取q_b＝2×10^-14，q_d＝2×10^-15；E[·]表示求均值符号；

结合步骤71中获取的新息信息，DUKF的滤波过程可以描述为

Step 1：计算***状态向量预测值

为第k-1个跟踪历元时刻的通道i***状态向量；

Step 2：计算***状态向量预测值的协方差矩阵

为的协方差矩阵；

Step 3：根据VFLL是否有结果输出获得量测信息若只有数据支路鉴相器有输出时，只取鉴相器对应项的数值；

Step 4：计算DUKF的增益矩阵

Step 4：根据新息更新状态估计结果：

Step 5：更新状态估计协方差矩阵：

步骤73，根据状态估计结果获得载波NCO的频率控制字即

其中，表示向量的第2个元素，至此，完成了一次DUKF的滤波处理过程。

如图4所示，为本发明提供的GNSS信号的矢量和标量混合跟踪环路结构示意图，包括N个跟踪通道模块1和1个VFLL(2)；所述N个跟踪通道模块具有相同的结构，包括本地载波生成装置11、本地信号生成装置12、第一乘法器13、第二乘法器14、第一相关器15、第二相关器16、鉴别器17，18和DUKF装置19；所述本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，用于产生导频支路和数据支路的伪码信号，并生成本地复制信号；所述鉴别器包括鉴相器17和鉴频器18，用于获取本地复制信号和接收信号之间的误差估计参数；所述本地载波生成装置11根据输入的频率控制字，产生同相载波信号和正交载波信号；所述数据支路信号生成装置和导频支路信号生成装置的输入端分别与本地载波生成装置的输出端相连；所述数据支路信号生成装置的输出端连接第一乘法器13的输入端，并将混合信号输出至第一相关器15的输入端；所述导频支路信号生成装置的输出端连接第二乘法器14的输入端，并将混合信号输出至第二相关器16的输入端；所述第一相关器15的输出端连接鉴相器17的输入端；鉴相器的输出端连接DUKF转置19的输入端；所述DUKF装置的输出端连接本地载波生成装置11；所述第二相关器16的输出端连接鉴频器18的输入端；鉴频器的输出端连接所述VFLL(2)的输入端；所述VFLL(2)的输出端分别输出至每个跟踪通道模块中的DUKF装置19的输入端。

如图5为某仿真场景下的GPS卫星星空图，其中共有8颗可见卫星，卫星PRN号分别为4，9，14，18，19，21，22，24。

如图6为本实施例在图4星空图下，利用标量跟踪环路和混合跟踪环路对4号卫星的跟踪结果，其中所有可见卫星在在前20s内信号强度为35dBHz，从20s到60s内，4号星和9号星的信号强度降至5dBHz，60s后恢复正常，从40s到80s内14号和18号星的信号强度降低至5dBHz，80s后恢复正常。图中DU-STL(10,20)表示采用双速率标量跟踪环路，数据支路环路更新间隔为10ms，导频支路环路更新间隔为20ms，DU-HTL(10,20)表示采用双速率混合跟踪环路(本发明方法)，数据支路环路更新间隔为10ms，导频支路更新间隔为20ms，同样的DU-HTL(10,50)和DU-HTL(10,100)分别对应导频支路更新间隔为50ms或100ms的双速率混合跟踪环路。从图中的跟踪结果可以看出在20s到60s，DU-STL和DU-HTL两种方法均无法保持对载波相位的正常锁定，但是从图中信号载波频率的跟踪误差结果可以看出，对于DU-HTL方法而言，通过利用其中的VFLL环路，可以使得信号的载波频率的跟踪误差保持在一定范围内，使得信号处于频率锁定状态，但是对于DU-STL方法而言，依然无法保证信号的载波频率处于锁定状态，当信号在60s的时候恢复至35dBHz时，DU-HTL方法可以快速重新锁定信号的载波相位，但是此时DU-STL方法却依然误差正常跟踪该卫星信号，从而验证了DU-HTL比DU-STL方法具有更好的信号重捕获和跟踪连续性能。对比不同参数下的DU-HTL跟踪结果可以看出，当VFLL的更新间隔越大，其频率跟踪误差也越大。

综上所述，虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明，任何本领域普通技术人员，在不脱离本发明的精神和范围内，当可作各种更动与润饰，因此本发明的保护范围当视权利要求书界定的范围为准。

Claims (4)

1.一种GNSS信号的矢量和标量混合跟踪方法，其特征在于，包括以下步骤：

步骤1，GNSS信号依次经过接收机中的天线，射频前端，AD转换器后变为数字中频信号r(t)；

步骤2，接收机有N个跟踪通道，每个跟踪通道中的处理方法相同，对于任意跟踪通道i中的本地载波生成装置NCO(Numerically Controlled Oscillator，缩写NCO)，生成频率控制字为的两路信号，分别为同相载波信号和正交载波信号t表示时间，具体为：

$s_I^i\left(t\right)=cos(2{\mathrm{\πf}}_{NCO}^i\·t)$

$s_Q^i\left(t\right)=-\sin (2{\mathrm{\πf}}_{NCO}^i\·t)$

步骤3，接收机跟踪通道i中的本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，数据支路信号生成装置接收同相载波信号和正交载波信号分别与数据支路本地伪码c_d(t)相乘，产生数据支路的本地同相信号和正交信号导频支路信号生成装置接收同相载波信号和正交载波信号分别与导频支路的本地伪码c_p(t)相乘，产生导频支路的本地同相信号和正交信号信号称为本地复制信号，具体为：

$s_{Id}^i\left(t\right)=s_I^i\left(t\right)\·c_d\left(t\right)$

$s_{Qd}^i\left(t\right)=s_Q^i\left(t\right)\·c_d\left(t\right)$

$s_{Ip}^i\left(t\right)=s_I^i\left(t\right)\·c_p\left(t\right)$

$s_{Qp}^i\left(t\right)=s_Q^i\left(t\right)\·c_p\left(t\right)$

步骤4，接收机跟踪通道i中的相关器进行相关处理，用于将本地复制信号分别与数字中频信号r(t)相乘后的混合信号进行相干积累，设相干积分时间为T_c，对任意通道i，输出相关值为其中下标k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为T_c，故输出信号的积分区间为(k-1)·T_c到k·T_c，具体结果如下：

$I_{p,k}^i={\∫}_{(k-1)\·T_c}^{k\·T_c}{s}_{Ip}^i\left(t\right)\·r\left(t\right)dt$

$Q_{p,k}^i={\∫}_{(k-1)\·T_c}^{k\·T_c}{s}_{Qp}^i\left(t\right)\·r\left(t\right)dt$

$I_{d,k}^i={\∫}_{(k-1)\·T_c}^{k\·T_c}{s}_{Id}^i\left(t\right)\·r\left(t\right)dt$

$Q_{d,k}^i={\∫}_{(k-1)\·T_c}^{k\·T_c}{s}_{Qd}^i\left(t\right)\·r\left(t\right)dt$

步骤5，接收机跟踪通道中的鉴别器装置对所述步骤4中输出的相关值进行处理，用于获得本地复制信号和数字中频信号间的误差估计参数，所述鉴别器包括数据支路鉴相器和导频支路鉴频器，经过鉴相器和鉴频器处理后，输出误差估计参数分别为和

${\mathrm{\ϵ}}_{d,k}^i=atan(Q_{d,k}^i/I_{d,k}^i)$

${\mathrm{\ϵ}}_{p,k}^i=\frac{atan2(I_1{Q}_2-Q_1{I}_2,I_1{I}_2+Q_1{Q}_2)}{N_p\·T_c}$

其中

$I_1=\underset{m=k-2\·N_p+1}{\overset{k-N_p}{\Σ}}{I}_{p,m}^i$

$I_2=\underset{m=k-N_p+1}{\overset{k}{\mathrm{\Σ}}}{I}_{p,m}^i$

$Q_1=\underset{m=k-2\·N_p+1}{\overset{k-N_p}{\Σ}}{Q}_{p,m}^i$

$Q_2=\underset{m=k-N_p+1}{\overset{k}{\Σ}}{Q}_{p,m}^i$

其中atan表示反正切函数，atan2表示四象限反正切函数；N_p为导频支路相干积累次数；

步骤6，接收机中的矢量频率跟踪环路(Vector Frequency Lock loop，缩写：VFLL)对各个通道中鉴频器输出结果进行处理，获得精度更高的频率误差估计结果其精度值分别为

步骤7，接收机中的任意通道i中的DUKF装置(Dual Update-rate Kalman Filter，缩写DUKF)根据步骤6中输出的频率误差估计结果和精度值，获得本地信号载波频率估计参数，并将载波频率估计参数输出至本地载波生成装置NCO，用于更新频率控制字。

2.如权利要求1所述的一种GNSS信号的矢量和标量混合跟踪方法，其特征在于，所述步骤6中的矢量频率跟踪环路处理过程如下：

步骤61，根据各个通道中鉴频器的输出结果获取量测量Z_k及其噪声协方差矩阵R_z，其中

VFLL的量测方程为：

$Z_k=H^V\·X_k^V+w_k^V$

其中为第k跟踪历元接收机的运动状态，δv_x,δv_y,δv_z为ECEF坐标系下的三维速度误差，δa_x,δa_y,δa_z为ECEF坐标系下的三维加速度误差，δf为接收机上时钟的频率误差；H^V为测量矩阵；为量测噪声，其协方差矩阵为其中diag()表示对角矩阵算子，为导频支路通道i中鉴频器输出结果ε_p ⁱ的噪声方差，具体为

$R_p^i={\left(\frac{1}{2{\mathrm{\πT}}_c}\right)}^2\·\frac{1}{C^i/N_0\·N_p\·T_c}\·(1+\frac{1}{C^i/N_0\·N_p\·T_c})$

其中Cⁱ/N₀表示通道i对应的信号载噪比，即信号功率与噪声的功率谱密度的比值；

步骤62，VFLL的迭代过程，具体描述如下：

***方程为

$X_{k+1}^V={\mathrm{\Φ}}^V\·X_k^V+{\mathrm{\ω}}_k^V$

其中Φ^V为状态转移矩阵，具体表示为

${\mathrm{\Φ}}^V={\left[\begin{array}{cccc}{\mathrm{\Φ}}_x& 0& 0& 0\\ 0& {\mathrm{\Φ}}_y& 0& 0\\ 0& 0& {\mathrm{\Φ}}_z& 0\\ 0& 0& 0& 1\end{array}\right]}_{7\×7}$

其中

${\mathrm{\Φ}}_x={\mathrm{\Φ}}_y={\mathrm{\Φ}}_z=\left[\begin{array}{cc}1& T_b\\ 0& 1\end{array}\right]$

T_b＝N_p·T_c,为VFLL的更新间隔；

是***过程噪声，其协方差矩阵为Q^V，具体为：

$Q^V={\left[\begin{array}{cccc}{Q}_x& 0& 0& 0\\ 0& Q_y& 0& 0\\ 0& 0& Q_z& 0\\ 0& 0& 0& Q_f\end{array}\right]}_{7\×7}$

其中

$Q_x=Q_y=Q_z=S_a\left[\begin{array}{cc}{T}_b^3/3& T_b^2/2\\ T_b^2/2& T_b\end{array}\right]$

Q_f＝S_f·T_b

S_a为加速度噪声功率谱密度，S_f为时钟频率变化噪声功率谱密度；

根据步骤61中获取的量测信息，得到VFLL的滤波步骤如下：

Step 1,计算接收机状态矢量预测值及其协方差值

$X_k^{V-}={\mathrm{\Φ}}^V\·X_{k-1}^V$

$P_k^{V-}={\mathrm{\Φ}}^V\·P_{k-1}^V\·{\left({\mathrm{\Φ}}^V\right)}^T+Q^V$

为对应的协方差矩阵；

Step 2,计算VFLL的增益矩阵

$G_k^V=P_k^{V-}\·{\left(H^V\right)}^T{(H^V\·P_k^{V-}\·{\left(H^V\right)}^T+R_z)}^{-1}$

Step 3,更新接收机状态矢量及其协方差阵

$X_k^V=X_k^{V-}+G_k^V\·(Z_k-H^V\·X_k^{V-})$

$P_k^V=(I-G_k^V\·H^V)\·P_k^{V-}$

其中I表示单位矩阵；

Step 4,计算各个通道频率估计误差

${\hat{Z}}_k=H^V\·X_k^V$

其估计精度满足

$P_k^z=H^V\·P_k^V\·{\left(H^V\right)}^T$

这样对于任意i通道而言，频率估计误差为估计精度为 表示向量的第i个元素,表示矩阵的第i行第i列元素值。

3.如权利要求1所述的一种GNSS信号的矢量和标量混合跟踪方法，其特征在于：所述步骤7中的DUKF装置获得本地信号载波频率估计参数具体过程如下：

步骤71，根据步骤5中鉴别器的输出结果和步骤6中VFLL的输出结果获取DUKF新息增量量测矩阵量测噪声矩阵当两种鉴别器的结果均有效时，计算公式为

$\left\{\begin{array}{c}{Y}_k^i=\[{\mathrm{\ϵ}}_{d,k}^i;{\widehat{\mathrm{\ϵ}}}_{p,k}^i\]\\ H_k^i=\[H_d;H_p\]\\ R_k^i=diag(\[R_d^i;{\hat{R}}_{p,k}^i\])\end{array}\right.$

其中为步骤5中鉴相器的输出结果,为步骤62输出的频率误差估计结果，H_d和H_p分别为两种不同量测量对应的量测矩阵，具体为

H_d＝[1 -T_c/2 T_c ²/6]

H_p＝[0 1 -(N_p-2)·T_c/2]

为数据支路鉴相器输出结果的噪声方差，具体为

$R_d^i={\left(\frac{1}{2\mathrm{\π}}\right)}^2\·\frac{1}{2\·C^i/N_0\·T_c}\·\left(1+\frac{1}{2\·C^i/N_0\·T_c}\right)$

为步骤62输出的频率误差估计精度；

当只有数据支路鉴相器有效时，只取鉴相器对应项的数值；

步骤72，对于任意通道i而言，其DUKF的迭代过程，具体描述如下：

DUKF的***方程为

$X_{k+1}^i=\mathrm{\Φ}\·X_k^i+w_k$

其中为第k跟踪历元通道i***状态向量，分别为信号的载波相位，多普勒频率和多普勒频率变化率，单位分别为周，Hz，Hz/s；w_k＝[ω_rf·w_b；ω_rf·w_d；(ω_rf/c)·w_a]^T为***噪声,w_b和w_d分别为由接收机中晶体振荡器的引起的相位噪声和频率噪声,其噪声谱密度分别为q_b和q_d；w_a是***频率变化率噪声，其功率谱密度为q_a；ω_rf表示载波频率，c为光速；Φ是***状态转移矩阵，具体为

$\mathrm{\Φ}=\left[\begin{array}{ccc}1& T_c& T_c^2/2\\ 0& 1& T_c\\ 0& 0& 1\end{array}\right]$

w_k为***过程噪声，Q是w_k对应的过程噪声协方差矩阵，具体为

$\begin{array}{c}Q=E\[w_k\·w_k^T\]\\ ={\mathrm{\ω}}_{rf}^2{q}_b\left[\begin{array}{ccc}{T}_c& 0& 0\\ 0& 0& 0\\ 0& 0& 0\end{array}\right]+{\mathrm{\ω}}_{rf}^2{q}_d\left[\begin{array}{ccc}{T_c}^3/3& {T_c}^2/2& 0\\ {T_c}^2/2& T_c& 0\\ 0& 0& 0\end{array}\right]\\ +{\left(\frac{{\mathrm{\ω}}_{rf}}{c}\right)}^2{q}_a\left[\begin{array}{ccc}{T_c}^5/20& {T_c}^4/8& {T_c}^3/6\\ {T_c}^4/8& {T_c}^3/6& {T_c}^2/2\\ {T_c}^3/6& {T_c}^2/2& T_c\end{array}\right]\end{array}$

结合步骤71中获取的新息信息，DUKF的滤波过程描述为

Step 1：计算***状态向量预测值

$X_k^{i-}=\mathrm{\Φ}\·X_{k-1}^i$

Step 2：计算***状态向量预测值的协方差矩阵

$P_k^{i-}=\mathrm{\Φ}\·P_{k-1}^i\·{\mathrm{\Φ}}^T+Q$

Step 3：根据VFLL是否有结果输出获得量测信息若只有数据支路鉴相器有输出时，只取鉴相器对应项的数值；

Step 4：计算DUKF的增益矩阵

$G_k^i=P_k^{i-}\·{\left(H_k^i\right)}^T\·{(R_k^i+H_k^i\·P_k^{i-}\·{\left(H_k^i\right)}^T)}^{-1}$

Step 4：根据新息更新状态估计结果

$X_k^i=X_k^{i-}+G_k^i\·Y_k^i$

Step 5：更新状态估计协方差矩阵

$P_k^i=P_k^{i-}-G_k^i\·H_k^i\·P_k^{i-}$

步骤73，根据状态估计结果获得载波NCO的频率控制字即

$f_{NCO}^i=X_k^i\left(2\right)$

其中，表示向量的第2个元素，至此，完成了一次DUKF的滤波处理过程。

4.一种GNSS信号的矢量和标量混合跟踪环路，其特征在于，包括N个跟踪通道模块(1)和1个VFLL(2)；所述N个跟踪通道模块具有相同的结构，包括本地载波生成装置NCO(11)、本地信号生成装置(12)、第一乘法器(13)、第二乘法器(14)、第一相关器(15)、第二相关器(16)、鉴别器(17，18)和DUKF装置(19)；所述本地信号生成装置包括数据支路信号生成装置和导频支路信号生成装置，用于产生导频支路和数据支路的伪码信号，并生成本地复制信号；所述鉴别器包括鉴相器(17)和鉴频器(18)，用于获取本地复制信号和接收信号之间的误差估计参数；

所述本地载波生成装置NCO(11)根据输入的频率控制字，产生同相载波信号和正交载波信号；所述数据支路信号生成装置和导频支路信号生成装置的输入端分别与本地载波生成装置NCO的输出端相连；

所述数据支路信号生成装置的输出端连接第一乘法器(13)的输入端，并将混合信号输出至第一相关器(15)的输入端；

所述导频支路信号生成装置的输出端连接第二乘法器(14)的输入端，并将混合信号输出至第二相关器(16)的输入端；

所述第一相关器(15)的输出端连接鉴相器(17)的输入端；鉴相器(17)的输出端连接DUKF转置(19)的输入端；所述DUKF装置的输出端连接本地载波生成装置NCO(11)；

所述第二相关器(16)的输出端连接鉴频器(18)的输入端；鉴频器(18)的输出端连接所述VFLL(2)的输入端；

所述VFLL(2)的输出端分别输出至每个跟踪通道模块中的DUKF装置(19)的输入端。