CN102621561A - 一种卫星导航接收机环路自调节方法 - Google Patents

Abstract

本发明公开一种卫星导航接收机环路自调节方法，在保留接收机原有跟踪环路结构的基础上增加判断以及对环路参数进行优化的步骤，能够根据载体动态信息自主调节环路阶数以及噪声带宽。该方法不需要额外增加完整的惯导系统。

Description

一种卫星导航接收机环路自调节方法

技术领域

本发明属于卫星信号接收机跟踪环路优化方法的技术领域，涉及一种利用三轴加速度计芯片输出信息作为环路参数优化调整的判断依据，使跟踪环路对自身阶数以及带宽进行自主调节的跟踪环路性能优化方法。

背景技术

卫星导航技术自引入我国以来，其应用领域不断扩展，在民用、军事、航天等方面均存在极大的运用前景。当搭载软件接收机的载体具有较高动态时，接收机与卫星之间的相对速度、加速度以及加速度变化率会导致其接收到的卫星信号在载波上存在很大且随时间变化的多普勒频移，当多普勒频移足够大时，接收机的跟踪环路将不能保持锁定，致使跟踪环路失锁。故提高环路的动态适应性能是实现高性能卫星导航接收机的关键技术之一，要提高跟踪环路的动态适应性能，究其根本就是要使高动态情况下的环路跟踪误差尽量小。

目前，针对高动态卫星信号的跟踪，主要处理方法有两种：一种是利用惯性器件进行辅助，一种是采用一些复杂的处理算法。利用惯性器件对跟踪环路进行辅助是补偿由于载体动态导致的卫星信号动态的最直接也是较为有效的方法，能够提高接收机对动态的适应能力和抗干扰能力，不过该方法实现复杂，且需要改变接收机原有跟踪环路的结构，实现起来较为不便；在无法利用惯性器件提供多普勒频移先验信息时，往往只能采用最大似然估计、扩展卡尔曼滤波自适应最小平均等一些复杂算法，这些算法将导致运算量的明显增加。

由于跟踪环路误差多由环路参数与载体动态不匹配所导致，故调节环路参数使其能够较好地对卫星信号进行跟踪是一种简单且易实现的跟踪环路性能优化方法，本案即是基于前述思想提出。

发明内容

本发明的目的，在于提供一种卫星导航接收机环路自调节方法，其不需要额外增加完整的惯导系统，在保留接收机原有跟踪环路结构的基础上仅需增加判断以及对环路参数进行优化的步骤，能够根据载体动态信息自主调节环路阶数以及噪声带宽。

为了达成上述目的，本发明的解决方案是：

一种卫星导航接收机环路自调节方法，包括如下步骤：

(1)设置环路初始阶数和初始带宽，采集并存储加速度信息，计算加速度变化率；

(2)读取经过射频前端处理的1ms卫星数字中频信号；

(3)分别通过码生成器和载波发生器产生与步骤(2)中卫星信号的输入码和输入载波尽可能一致的码信号和载波信号，并生成对应的超前、滞后码；

(4)通过混频机制利用步骤(3)中的本地载波信号剥离步骤(2)中卫星信号中的载波，利用得到的两路混频结果分别与步骤(3)中生成的超前、即时和滞后码做相关运算，得到六路相关积分结果；

(5)根据步骤(4)中得到的相关积分结果估算载噪比，并通过二象限反正切函数鉴相器进行载波鉴相，得到鉴相器输出；

(6)载波跟踪环路参数自调节步骤，该步骤分为以下子步骤：

(61)环路锁定检测：利用步骤(4)中跟踪环路相关积分结果以及步骤(5)得到的载噪比估计值信息进行环路跟踪状态的综合判别，若该通道信号该时刻发生失锁现象，则转步骤(8)，若判断锁定信号，则继续执行步骤(62)；

(62)环路阶数调节：通过步骤(1)得到的加速度、加速度变化率判断所需的环路阶数，并计算平均加速度以及平均加速度变化率，当平均加速度变化率小于门限值且平均加速度不为0时，若当前环路阶数高于载体动态对应的理想阶数，则逐次向下调节环路阶数至理想阶数，若当前环路阶数即为理想阶数则保持不变；当平均加速度变化率大于门限值时，若当前环路阶数低于载体动态对应的理想阶数，则逐次向上调节环路阶数至理想阶数，若当前环路结束即为理想阶数则保持不变；若非以上两种情况，则跳出该步骤，直接进入步骤(8)；

(63)环路带宽调节：利用N阶载波跟踪环路最优带宽计算公式计算当前载波跟踪环路的最优带宽，调节环路带宽为最优带宽；

(7)利用步骤(6)调节过的环路滤波器对步骤(5)得到的载波鉴相结果进行滤波，滤波结果通过载波NCO对本地载波频率进行修正；

(8)码鉴相、鉴相结果滤波及码频率修正步骤：利用对步骤(4)中得到的相关积分结果通过非相干超前减滞后幅值法对码环进行鉴相，得到鉴相结果，将该鉴相结果经过一阶低通滤波器滤波后，控制码环NCO对本地码频率进行修正。

上述步骤(1)中，以最小采样周期采集并存储三轴加速度计芯片输出的加速度信息a，利用差分方法求得加速度变化率

其中a_k和

分别为k时刻载体加速度以及加速度变化率，ΔT为k、(k-1)时刻间的时间间隔。

上述步骤(4)中，六路相关积分结果的表达式分别如下：

$I_E=\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})T$

$I_P=\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)T$

$I_L=\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})T$

$Q_E=-\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})\mathrm{\Δ\θT}$

$Q_P=-\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)\mathrm{\Δ\θT}$

$Q_L=-\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})\mathrm{\Δ\θT}$

其中，I_E、I_P、I_L分别为同相支路超前、即时、滞后相干积分值，Q_E、Q_P、Q_L分别为正交支路超前、即时、滞后相干积分值，A表示信号强度，Δt表示本地码与输入码之间的相位差，C(Δt)表示相位差Δt引起的相关峰值，Δω和Δθ分别表示本地载波与输入载波之间的频率差以及相位差，τ为超前滞后的码片移动大小，T为相关积分时间。

上述步骤(5)中，估算载噪比的步骤为：

(51)构造I、Q平方和序列步骤：对每个时刻的相关积分结果I_P、Q_P计算其平方和，即每1ms进行一次，构造序列Z，其中

i为表征当前时刻的下标；

(52)载波功率均值计算步骤：计算步骤(51)中序列Z的均值

以及方差

其中 $\stackrel{\‾}{Z}=(1/i){\mathrm{\Σ}}_i{Z}_i,$ ${\mathrm{\σ}}_Z^2=\left(\frac{1}{k-1}\right){\mathrm{\Σ}}_i{(Z_i-\stackrel{\‾}{Z})}^2,$ 在此基础上利用计算载波功率均值，N表示序列Z中样本总数，A_f表示前端输出信号幅度；

(53)I、Q噪声方差

估算步骤：利用步骤(52)中得到的序列Z的均值

以及方差

计算I、Q噪声方差 ${\mathrm{\σ}}_{\mathrm{IQ}}^2=\frac{1}{2}(\stackrel{\‾}{Z}-\sqrt{{\stackrel{\‾}{Z}}^2-{\mathrm{\σ}}_Z^2});$

(54)载噪比C/N₀估算步骤：在由步骤(52)得到

由步骤(63)得到已知积分时间T的前提下，估算载噪比

上述步骤(5)中，利用步骤(4)得到的相关积分结果I_P和Q_P通过二象限反正切函数鉴相器进行载波鉴相，得到鉴相器输出

采用上述方案后，本发明提出一种在不改变卫星信号接收机原有跟踪环路结构的基础上，综合三轴加速度计芯片输出以及对信号载噪比估算得到的结果，对环路阶数、带宽进行实时调节，以保持信号的跟踪同时使跟踪误差最小的方法。该方法不需要完整的惯性导航模块，只需要一块加速度计芯片，实现灵活、成本低，较传统跟踪环路实现没有明显运算量增加，能够提高接收机跟踪环路性能。

附图说明

图1是本发明的原理流程示意图；

图2是本发明中测试信号载波跟踪环路鉴相器输出示意图；

图3是本发明中测试信号载波跟踪环路频率误差示意图。

具体实施方式

以下将结合附图，对本发明的技术方案进行详细说明。

本发明提供一种卫星导航接收机环路自调节方法，跟踪环路运行流程如图1所示。首先设置环路的初始跟踪带宽以及阶数，使能三轴加速度计芯片的工作。加速度计芯片每个采样周期输出载体的加速度值，该方法对加速度计芯片输出进行处理，并存储加速度、加速度变化率值，作为环路阶数判断的依据。每个环路处理周期内，对每个通道进行如下循环处理：读取1ms数据，建立本地码信息以及本地载波，将本地载波与输入中频信号混频，随后进行相关积分处理，再利用相关积分结果进行鉴相，如果正在处理的是第一个通道，则需要进行环路参数自调节的操作，然后再对鉴相结果进行滤波并修正本地信号，若非第一个通道则可直接对鉴相结果滤波并修正本地信号。

以下将结合具体实施例对本发明作详细介绍。

本发明为一种卫星导航接收机环路自调节方法，包括如下步骤：

(1)环路初始设置及使能步骤：设置环路初始阶数为二阶、初始带宽为20Hz，环路开始工作的同时，三轴加速度计芯片也开始工作；

(2)加速度计芯片输出数据处理步骤：以最小采样周期采集并存储三轴加速度计芯片输出的加速度信息a，利用差分方法求得加速度变化率

其中a_k和

分别为k时刻载体加速度以及加速度变化率，ΔT为k、(k-1)时刻间的时间间隔(即加速度计芯片采样周期)；

(3)卫星信号读取：读取1ms的经过射频前端处理的卫星数字中频信号(跟踪环路的循环周期为1ms)；

(4)本地信号产生步骤：通过码生成器产生一个与步骤(3)读取的卫星信号中的输入码尽可能一致的码信号，并生成对应的超前、滞后码，通过载波发生器生成与步骤(3)读取的卫星信号中的输入载波尽可能一致的载波信号；

(5)混频及相关累加步骤：通过混频机制利用步骤(4)产生的本地载波尽可能地剥离步骤(3)读取的卫星信号中的载波，利用得到的两路混频结果分别与步骤(4)中生成的超前、即时以及滞后码作相关运算，得到六路相关积分结果：

$I_E=\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})T$

$I_P=\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)T$

$I_L=\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})T$

$Q_E=-\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})\mathrm{\Δ\θT}$

$Q_P=-\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)\mathrm{\Δ\θT}$

$Q_L=-\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})\mathrm{\Δ\θT}$

其中，I_E、I_P、I_L分别为同相支路超前、即时、滞后相干积分值，Q_E、Q_P、Q_L分别为正交支路超前、即时、滞后相干积分值，A表示信号强度，Δt表示本地码与输入码之间的相位差，C(Δt)表示相位差Δt引起的相关峰值，Δω和Δθ分别表示本地载波与输入载波之间的频率差以及相位差，τ为超前滞后的码片移动大小，T为相关积分时间；

(6)载波鉴相步骤：利用步骤(5)得到的相关积分结果I_P和Q_P通过二象限反正切函数鉴相器进行载波鉴相，得到鉴相器输出

(7)载噪比估算步骤：利用步骤(5)中得到的相关积分结果I_P、Q_P进行载噪比估计，该步骤分为以下子步骤：

(71)构造I、Q平方和序列步骤：对每个时刻的相关积分结果I_P、Q_P计算其平方和，即每1ms进行一次，构造序列Z，其中

i为表征当前时刻的下标；

(72)载波功率均值计算步骤：计算步骤(71)中序列Z的均值以及方差其中 $\stackrel{\‾}{Z}=(1/i){\mathrm{\Σ}}_i{Z}_i,$ ${\mathrm{\σ}}_Z^2=\left(\frac{1}{k-1}\right){\mathrm{\Σ}}_i{(Z_i-\stackrel{\‾}{Z})}^2,$ 在此基础上利用

计算载波功率均值，N表示序列Z中样本总数，本实施例中取N＝1000，即每1s能够得到一个载噪比估计值，A_f表示前端输出信号幅度；

(73)I、Q噪声方差

估算步骤：利用步骤(62)中得到的序列Z的均值

以及方差

计算I、Q噪声方差 ${\mathrm{\σ}}_{\mathrm{IQ}}^2=\frac{1}{2}(\stackrel{\‾}{Z}-\sqrt{{\stackrel{\‾}{Z}}^2-{\mathrm{\σ}}_Z^2});$

(74)载噪比C/N₀估算步骤：在由步骤(72)得到

由步骤(73)得到

已知积分时间T的前提下，估算载噪比

(8)载波跟踪环路参数自调节步骤，该步骤分为以下子步骤：

(81)环路锁定检测：利用步骤(5)中跟踪环路相关积分结果以及步骤(7)得到的载噪比估计值信息进行环路跟踪状态的综合判别，判别方法为：取当前时刻前14个相邻积分周期的相关积分结果IP值，顺序取连续的5个积分周期的相关积分结果IP作为一组(可构成10组)，判断每组中5个IP值的符号，如出现两次或以上符号变化的组数大于5，且载噪比下降大于一定阈值(阈值一般可设为20)，则判定该通道信号该时刻发生失锁现象，若判断出信号失锁，则跳出该步骤，直接进行到步骤(10)，若判断锁定信号，则继续执行步骤(82)；

(82)环路阶数调节：通过步骤(2)得到的加速度a、加速度变化率

判断所需的环路阶数(即鉴相滤波器阶数+1)，取当前时刻往前的10组加速度、加速度变化率值做平均，得到平均加速度

以及平均加速度变化率 $\stackrel{\‾}{\stackrel{\·}{a}}=\frac{({\stackrel{\·}{a}}_{k-9}+{\stackrel{\·}{a}}_{k-8}+...+{\stackrel{\·}{a}}_{k-1}+{\stackrel{\·}{a}}_k)}{10},$ 若

(

为用于判断载体是否存在加速度变化率的门限值)且

则调节环路阶数为二阶(使用一阶滤波器)，若

则调节环路阶数为三阶(使用二阶滤波器)，若非以上两种情况则跳出该步骤，直接进行步骤(10)；

(83)环路带宽调节：利用N阶载波跟踪环路最优带宽计算公式

计算当前载波跟踪环路的最优带宽，调节环路带宽为最优带宽，其中：N为步骤(82)中确定的环路阶数，α为对应不同环路阶数的系数：一阶环路α＝1/4，二阶环路三阶环路α＝5/6，R^(N)表示卫星与接收机间连线距离对时间的N次导数其一、二、三阶导分别与载体的速度、加速度、加速度变化率对应，c/n₀为信号载噪比，T为相关积分时间，λ_L为载波波长；

(9)对载波鉴相结果滤波及本地载波频率修正步骤：利用步骤(8)调节过的环路滤波器对步骤(6)得到的载波鉴相结果进行滤波，滤波结果通过载波NCO对本地载波频率进行修正；

(10)码鉴相、鉴相结果滤波及码频率修正步骤：利用对步骤(5)中得到的I_E、I_L、Q_E、Q_L通过非相干超前减滞后幅值法对码环进行鉴相，得到鉴相结果

鉴相结果δ_cp经过一阶低通滤波器滤波后，控制码环NCO对本地码频率进行修正。

跟踪环路自调节方法仿真实验效果：

设计测试场景：载体以1g的初始加速度做匀加速直线运动，10s后加载2g/s的加速度变化率做变加速直线运动。从载波鉴相器输出图2以及载波跟踪环路频率误差图3可以看出，调节阶数、带宽后，环路能够保持稳定跟踪且总体性能有所提高。

综上，本发明针对载体高动态环境下由接收机卫星间相对速度、加速度以及加速度变化率导致的较大的且随时间变化的多普勒频移极易导致接收机跟踪环路不能稳定工作的问题，在基本不改变接收机跟踪环路原有结构的基础上，充分利用环路内部信息以及加速度计芯片输出进行环路参数的优化调节，有效提高了跟踪环路的动态性能，为实现适应高动态环境的接收机跟踪过程提供了一种新的、经济且易于实现的方法，具有很强的工程应用价值。

以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (5)

1.一种卫星导航接收机环路自调节方法，其特征在于包括如下步骤：

(1)设置环路初始阶数和初始带宽，采集并存储加速度信息，计算加速度变化率；

(2)读取经过射频前端处理的1ms卫星数字中频信号；

(3)分别通过码生成器和载波发生器产生与步骤(2)中卫星信号的输入码和输入载波尽可能一致的码信号和载波信号，并生成对应的超前、滞后码；

(4)通过混频机制利用步骤(3)中的本地载波信号剥离步骤(2)中卫星信号中的载波，利用得到的两路混频结果分别与步骤(3)中生成的超前、即时和滞后码做相关运算，得到六路相关积分结果；

(5)根据步骤(4)中得到的相关积分结果估算载噪比，并通过二象限反正切函数鉴相器进行载波鉴相，得到鉴相器输出；

(6)载波跟踪环路参数自调节步骤，该步骤分为以下子步骤：

(61)环路锁定检测：利用步骤(4)中跟踪环路相关积分结果以及步骤(5)得到的载噪比估计值信息进行环路跟踪状态的综合判别，若该通道信号该时刻发生失锁现象，则转步骤(8)，若判断锁定信号，则继续执行步骤(62)；

(62)环路阶数调节：通过步骤(1)得到的加速度、加速度变化率判断所需的环路阶数，并计算平均加速度以及平均加速度变化率，当平均加速度变化率小于门限值且平均加速度不为0时，若当前环路阶数高于载体动态对应的理想阶数，则逐次向下调节环路阶数至理想阶数，若当前环路阶数即为理想阶数则保持不变；当平均加速度变化率大于门限值时，若当前环路阶数低于载体动态对应的理想阶数，则逐次向上调节环路阶数至理想阶数，若当前环路结束即为理想阶数则保持不变；若非以上两种情况，则跳出该步骤，直接进入步骤(8)；

(63)环路带宽调节：利用N阶载波跟踪环路最优带宽计算公式计算当前载波跟踪环路的最优带宽，调节环路带宽为最优带宽；

(7)利用步骤(6)调节过的环路滤波器对步骤(5)得到的载波鉴相结果进行滤波，滤波结果通过载波NCO对本地载波频率进行修正；

(8)码鉴相、鉴相结果滤波及码频率修正步骤：利用对步骤(4)中得到的相关积分结果通过非相干超前减滞后幅值法对码环进行鉴相，得到鉴相结果，将该鉴相结果经过一阶低通滤波器滤波后，控制码环NCO对本地码频率进行修正。

2.如权利要求1所述的一种卫星导航接收机环路自调节方法，其特征在于：所述步骤(1)中，以最小采样周期采集并存储三轴加速度计芯片输出的加速度信息a，利用差分方法求得加速度变化率

其中a_k和

分别为k时刻载体加速度以及加速度变化率，ΔT为k、(k-1)时刻间的时间间隔。

3.如权利要求2所述的一种卫星导航接收机环路自调节方法，其特征在于：所述步骤(4)中，六路相关积分结果的表达式分别如下：

$I_E=\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})T$

$I_P=\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)T$

$I_L=\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\sin (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\sin \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})T$

$Q_E=-\frac{\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}-\mathrm{\τ})\mathrm{\Δ\θT}$

$Q_P=-\frac{\mathrm{AC}\left(\mathrm{\Δt}\right)}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}\left(\mathrm{\Δt}\right)\mathrm{\Δ\θT}$

$Q_L=-\frac{\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})}{2\mathrm{\Δ\ω}}[\cos (\mathrm{\Δ\ωT}+\mathrm{\Δ\θ})-\cos \left(\mathrm{\Δ\θ}\right)]\≈\frac{1}{2}\mathrm{AC}(\mathrm{\Δt}+\mathrm{\τ})\mathrm{\Δ\θT}$

其中，I_E、I_P、I_L分别为同相支路超前、即时、滞后相干积分值，Q_E、Q_P、Q_L分别为正交支路超前、即时、滞后相干积分值，A表示信号强度，Δt表示本地码与输入码之间的相位差，C(Δt)表示相位差Δt引起的相关峰值，Δω和Δθ分别表示本地载波与输入载波之间的频率差以及相位差，τ为超前滞后的码片移动大小，T为相关积分时间。

4.如权利要求3所述的一种卫星导航接收机环路自调节方法，其特征在于：所述步骤(5)中，估算载噪比的步骤为：

(51)构造I、Q平方和序列步骤：对每个时刻的相关积分结果I_P、Q_P计算其平方和，即每1ms进行一次，构造序列Z，其中

i为表征当前时刻的下标；

(52)载波功率均值计算步骤：计算步骤(51)中序列Z的均值

以及方差

其中 $\stackrel{\‾}{Z}=(1/i){\mathrm{\Σ}}_i{Z}_i,$ ${\mathrm{\σ}}_Z^2=\left(\frac{1}{k-1}\right){\mathrm{\Σ}}_i{(Z_i-\stackrel{\‾}{Z})}^2,$ 在此基础上利用

计算载波功率均值，N表示序列Z中样本总数，A_f表示前端输出信号幅度；

(53)I、Q噪声方差

估算步骤：利用步骤(52)中得到的序列Z的均值

以及方差

计算I、Q噪声方差 ${\mathrm{\σ}}_{\mathrm{IQ}}^2=\frac{1}{2}(\stackrel{\‾}{Z}-\sqrt{{\stackrel{\‾}{Z}}^2-{\mathrm{\σ}}_Z^2});$

(54)载噪比C/N₀估算步骤：在由步骤(52)得到

由步骤(63)得到

已知积分时间T的前提下，估算载噪比

5.如权利要求3所述的一种卫星导航接收机环路自调节方法，其特征在于：所述步骤(5)中，利用步骤(4)得到的相关积分结果I_P和Q_P通过二象限反正切函数鉴相器进行载波鉴相，得到鉴相器输出

Images