一种导航接收机中检测观测量有效性的方法
【技术领域】
[0001] 本发明涉及导航接收机的基带信号处理算法，尤其涉及在导航接收机中用于检测 跟踪环路输出的观测量的有效性的方法。
【背景技术】
[0002] 导航接收机的内部结构按照工作流程的先后顺序，通常分为射频（RF)前端处理、 基带数字信号处理（DSP)和定位导航运算三大功能模块。其中，射频前端处理模块通过天 线接收所有可见卫星的信号，并将L波段的信号下变频到中频（或者接近基带），然后进行 采样与量化，将数字中频信号输入到基带处理模块；基带数字信号处理部分主要进行捕获、 跟踪可见卫星信号，对导航电文进行位同步、帧同步、解调电文信息，同时产生伪距与载波 相位两个基本距离测量值；将测量值与电文信息输入到定位导航运算模块，计算出可见卫 星位置，并由最小二乘法或者Kalman滤波估计出接收机位置。
[0003] 接收机对每颗卫星的伪距与载波相位测量值直接决定了定位解算的性能。测量值 的误差不仅仅来源于卫星导航系统的时钟误差、星历误差、电离层与对流层延时等基本的 误差源，同时随着接收机在室内、城市峡谷、树阴遮挡、高架桥下、干扰环境等中的应用，弱 信号、多路径以及干扰对测量值噪声带来了更重要的误差影响，所以在定位解算前，有必要 对测量值的准确度以及精度进行评估和检测。其主要内容包括两个方面：一、对于误差过 大，一旦参与定位便会产生不可靠结果的观测量，应该给予剔除，不参与定位解算；二、对于 误差在合理范围内的每颗卫星信号的观测量，应该评估各个观测量的好坏，进而在定位算 法中正确地区别对待不同质量的测量值，以提高定位性能。
[0004] 目前，导航接收机中检测观测量有效性的方法有两种途径：一、在跟踪环路中，加 入锁定判别器，对载波跟踪环和码跟踪环的性能进行检测，如果发现跟踪的载波频率或者 码相位失锁，则认为信号丢失并且观测量不可靠；二、在定位导航运算模块中，加入自主正 直性检测算法，利用伪距残差等信息，判断定位解算是否可靠。
[0005] 但是以上方法中，第一种方法只能检测跟踪环路性能，即码相位的准确性，而伪距 是由码相位，码片数，导航电文等信息构成，所以不能检测出整体伪距的正确性；第二种方 法基于定位解算中的冗余信息，要求至少5颗卫星才能判定定位是否有效，6颗卫星才能剔 除错误观测量，并且算法复杂，计算量大。
【发明内容】
[0006] 本发明要解决的技术问题为：提供一种新的导航接收机中检测观测量有效性的方 法，克服现有技术中存在的锁定判别器只能检测跟踪环路性能，不能检测出整体伪距的正 确性，基于定位解算中的冗余信息的检测要求至少5颗卫星才能判定定位是否有效，6颗卫 星才能剔除错误观测量，并且算法复杂，计算量大的问题。
[0007] 本发明的技术方案为：
[0008] -种导航接收机中检测观测量有效性的方法，包括：S1、信号能量检测模块根据接 收到的信号的能量，计算出信号的载噪比；S2、载波跟踪锁定检测模块判断载波跟踪的稳定 性，如果判断结果为"稳定"，则进入步骤S3,如果判断结果为"不稳定"，则进入步骤S8 ;S3、 观测量收集模块在每一历元提取载波跟踪环的载波多普勒频率观测量以及码跟踪环的伪 距观测量；S4、利用载波多普勒和伪距两种观测量，分别计算出对应的前后两个历元的伪距 差，即Δ伪距；S5、伪距观测量检测模块计算两种Δ伪距的差值，并与门限比较，判断当前 历元伪距观测量是否正确，如果正确则进入步骤S6、否则进入步骤S8 ;S6、将当前历元的Δ 伪距之差与之前η次历元的Δ伪距之差求和，并与门限比较，判断伪距观测量是否正确， 如果正确则进入步骤S7、否则进入步骤S8 ;S7、对每一颗保持跟踪的卫星均进行步骤Sl至 S6,收集所有进行步骤S6并判断观测量正确的卫星信号，输入到定位解算模块，进行接收 机位置、速度、时间的解算；S8,当S2判断结果为不稳定、S5判断结果为不正确、S6判断结 果为不正确时，进行对变量的复位，舍弃观测量。
[0009] 进一步地，包括：步骤S2中载波跟踪锁定检测模块判断载波跟踪的稳定性包 括，将一个数据比特20ms长度分成前后两个10ms，I支路的前一个IOms的累加值为 I1GO，后一个IOms的累加值为I2GO，Q支路的前一个IOms的累加值为Q1GO，后一个 IOms 的累加值为 Q2 (k);求出叉积 Cross = I1 (k) *Q2 (k) -I2 (k) =I=Q1 (k)，以及点积 Dot = I1 (k) *I2 (k) +Q1 (k) *Q2 (k);在Is内对Cross与Dot进行50次累加，然后对两个累加值Crossli 与0的1(求四象限反正切ATAN2(Cr〇SS K/DotK)，并且与门限T1比较，得出判断结果。
[0010] 进一步地，步骤S4中包括以下步骤：在历元n，利用卫星电文中的时间信息和码 跟踪环路得到的码相位观测量，计算出卫星到接收机的伪距Pr raffi(η)，前后两个历元的伪 距差即为由码观测量计算出的Δ伪距，g卩ΔΡι^ζΡΓα^ιΟ-Ρι^Οι-Ι);在历元η，利用 载波跟踪环路得到的载波积分多普勒观测量Φ(η)，计算出接收机的A伪距，g卩ΔΡ Γωκκγ =_(Φ (η) - IF)* λ，其中IF为载波中频频率，λ为光的波长。
[0011] 进一步地，步骤S5中伪距观测量检测模块包括以下步骤：计算两种Δ伪距 的差再求绝对值，即I Δ Prra ω (n) I = I Δ PrraDE (η) - Δ PrMRRY (η) I，再与门限T2比较；如 果大于门限，就判定伪距观测量不正确；如果小于门限，若n〈 = N，则将△ Pr03 ω (η) 与前 η-1 次的 Δ Prra CA? 行累加，求出 MUL_A Pr c。CA(n)，即 MUL_A Prra CA(n) = MUL_ Δ Pr03 ω (n-l) + Δ Pr03 ω (n)，再取绝对值与门限MUL_Tn进行比较，如果大于，就判断伪距观 测量不正确；若n>N，则将η置0,同时将MUL_APr ra ω(η)置0 ;并在下一个历元n+1，重复上 述步骤。
[0012] 进一步地，噪声检测通道为一个指定的跟踪通道，设置本地码相位为0,载波频率 为中频频率IF，使噪声检测通道与接收信号没有相关性。
[0013] 本发明与现有技术相比的优点在于：
[0014] 根据本发明的方法，可以监测载波跟踪环跟踪性能，依据载波跟踪环精度远高于 码跟踪环精度的条件，可以进一步利用载波积分多普勒监测码观测量的准确性，评估观测 量误差的大小，剔除错误的观测量。在不需要冗余信息的条件下，利用较小的计算量解决了 接收机在复杂环境下定位稳定性的问题。
【附图说明】
[0015] 图1是根据本发明的一种导航接收机中检测观测量有效性的方法的示意流程图；
[0016] 图2是本发明中使用的估计接收信号载噪比的方法示意图；
[0017] 图3是本发明中的载波跟踪锁定检测模块的方法示意图。
【具体实施方式】
[0018] 本发明提供一种导航接收机中检测观测量有效性的方法，其特征在于包括以下步 骤：
[0019] (1)信号能量检测模块估计接收到的信号的能量大小，计算出信号的载噪比；
[0020] (2)在载波跟踪环中，载波跟踪锁定检测模块判断载波跟踪的稳定性；
[0021] (3)观测量收集模块在每一历元提取载波跟踪环的载波多普勒频率观测量以及码 跟踪环的伪距观测量；
[0022] (4) △伪距计算模块利用载波多普勒和伪距两种观测量分别计算出对应的前后两 个历元的伪距差，即A伪距；
[0023] (5)伪距观测量检测模块计算这两种Δ伪距的差值，并与门限比较，判断当前历 元伪距观测量是否正确；
[0024] (6)再将当前历元的Δ伪距之差与之前η次历元的Δ伪距之差求和，并与门限比 较，判断伪距观测量是否正确。
[0025] 图1是按照本发明的一实施方式在导航接收机中检测观测量有效性的方法的示 意流程图。该流程开始步骤S101。然后在步骤S102,信号能量检测模块计算出接收信号的 载噪比，用于为本方法实施方式的以下步骤中的门限值设置提供依据。在本技术领域中已 存在现有技术，可以使用多种方法计算接收信号的载噪比，例如噪声通道法，窄带比宽带法 等。
[0026] 在步骤S103中，载波跟踪锁定检测模块用于判断载波跟踪的稳定性。导航接收机 的载波跟踪环路一般采用频率锁定环（FLL)跟踪信号的载波频率，或者采用COSTAS相位锁 定环（PLL)跟踪信号的载波相位，用于实现精确地跟踪卫星信号的多普勒频率以及载波相 位变化，并且产生的载波积分多普勒频率观测量精确的体现了伪距率（即A伪距）的大 小。
[0027] 载波跟踪锁定检测模块利用同相支路（I路）与正交支路（Q路）的相关值判断载 波相位以及频率的锁定状态。在相位锁定情况下，I路相关值包含了有用信号能量，而Q路 相关值只包含噪声，ATAN(Q/I)代表了鉴别相位的误差，应该在0值附近摆动；叉积Cross =IdQ2-Ii^Q1，以及点积Dot = I1Wi^QdQ2，下标1与2分别代表前一时刻与后一时刻的相 关值，则ATAN2 (Cross/Dot)代表了鉴别频率的误差，在频率锁定情况下，应该在0值附近摆 动。
[0028] 如果步骤S103的结果为"否"，则流程进入步骤SllO ;否则流程进入步骤S104。
[0029] 在步骤S104中，在历元n，观测量收集模块提取当前历元的载波跟踪环的载波积 分多普勒频率观测量Φ (η)以及码跟踪环的伪距观测量PrOT1E(n)。
[0030] 在步骤S105中，利用步骤S104得到的两个观测量分别求出伪距率。
[0031] 第一种方法利用前后历元时刻的伪距的差求出Δ伪距，即ΔΡΓ_ = Pr CODE (η) _ΡΓ CODE (η_ I)；
[0032] 第二种方法利用当前历元的载波积分多普勒观测量Φ (η)，计算Δ伪距，即 △PrCARRY= _(Φ (η) - IF)* λ，其中IF为没有多普勒效应的载波中频频率，λ为光的波长。
[0033] 在步骤S106中，计算两种Δ伪距的差再求绝对值，即I APrra ω(η) I = A PrCdDE (n)_ A PrCARRY (η) I 〇
[