CN109143276B - 一种基于cusum算法的组合抗多径鉴相器的设计方法 - Google Patents

一种基于cusum算法的组合抗多径鉴相器的设计方法 Download PDF

Info

Publication number
CN109143276B
CN109143276B CN201810156977.0A CN201810156977A CN109143276B CN 109143276 B CN109143276 B CN 109143276B CN 201810156977 A CN201810156977 A CN 201810156977A CN 109143276 B CN109143276 B CN 109143276B
Authority
CN
China
Prior art keywords
environment
phase
multipath
phase discriminator
cusum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810156977.0A
Other languages
English (en)
Other versions
CN109143276A (zh
Inventor
朱魁华
王超
邢超
张沛尧
刘强
董胜飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Gallop Sky Electronic Technology Co ltd
Shanghai Huace Navigation Technology Ltd
Original Assignee
Nanjing Gallop Sky Electronic Technology Co ltd
Shanghai Huace Navigation Technology Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Gallop Sky Electronic Technology Co ltd, Shanghai Huace Navigation Technology Ltd filed Critical Nanjing Gallop Sky Electronic Technology Co ltd
Priority to CN201810156977.0A priority Critical patent/CN109143276B/zh
Publication of CN109143276A publication Critical patent/CN109143276A/zh
Application granted granted Critical
Publication of CN109143276B publication Critical patent/CN109143276B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/22Multipath-related issues

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Noise Elimination (AREA)

Abstract

本发明提供了一种基于CUSUM算法的组合抗多径鉴相器的设计方法,包括以下步骤:步骤(1):利用CUSUM控制图算法的变点理论进行良好环境和多径环境的识别与监测;步骤(2):利用strobe抗多径鉴相器和常规鉴相器的综合使用,来减小strobe鉴相环路在常规环境中错锁的可能性,本发明利用了CUSUM算法对于微小异常变化的敏感性以及常规鉴相器与strobe鉴相器部分优缺点互补的特性,设计了一种实用性更高的组合鉴相器,能够在不中断跟踪环路的情况下,有效地抑制strobe鉴相环路在抗多径的同时又带来的可能的相位错锁情况的发生。

Description

一种基于CUSUM算法的组合抗多径鉴相器的设计方法
技术领域
本发明涉及GNSS接收机或通信设备接收机的基带信号处理领域,具体涉及到一种基于CUSUM算法的组合抗多径鉴相器的设计方法。
背景技术
多径信号是影响GNSS接收机定位精度的主要影响因素之一。综合考虑抗多径算法实现复杂度以及抗多径效果,类似于strobe鉴相器一类的多相关器鉴相器是目前GNSS接收机较为常用的抗多径手段之一,strobe鉴相器公式如下表示:
Figure BDA0001581707760000011
其中,dc为鉴相结果,IE1、QE1为超前支路1的IQ相关值,IE2、QE2为超前支路2的IQ相关值,IL1、QL1为滞后支路1的IQ相关值,IL2、QL2为滞后支路2的IQ相关值,IP、QP为即时支路的IQ相关值。
但是在实际的使用当中,由于strobe鉴相器本身有效鉴相范围较窄,且范围外的鉴相输出为0,这导致当信号载噪比较小,鉴相噪声较大时,strobe鉴相环路可能锁定在错误的码相位上,结果恶化了定位精度。缓解上述strobe鉴相器使用的不利情况会进一步增强该鉴相器的实用性。
发明内容
针对现有技术存在上述的不足,本发明提供了一种基于CUSUM算法的组合抗多径鉴相器的设计方法,本发明利用了CUSUM算法对于微小异常变化的敏感性以及常规鉴相器与strobe鉴相器部分优缺点互补的特性,设计了一种实用性更高的组合鉴相器,能够在不中断跟踪环路的情况下,有效地抑制strobe鉴相环路在抗多径的同时又带来的可能的相位错锁情况的发生。
本发明的技术方案为,一种基于CUSUM算法的组合抗多径鉴相器的设计方法,包括以下步骤:
步骤(1):利用CUSUM控制图算法的变点理论进行良好环境和多径环境的识别与监测;
步骤(2):利用strobe抗多径鉴相器和常规鉴相器的综合使用,来减小strobe鉴相环路在常规环境中错锁的可能性。
上述的设计方法,其中,所述步骤(1)中包括:以相关值的载噪比CN0以及常规鉴相器和strobe鉴相器的差值为统计量,利用CUSUM算法实时监测这两类统计量指标的变化情况,综合两类指标的变化情况识别当前环境是否处于良好环境。
上述的设计方法,其中,利用多径环境识别模块给出的环境判别结果,并依据两类鉴相器的鉴相差值与阈值的比较情况,进行鉴相器的切换,使得在良好环境中使用常规鉴相器,而在多径恶劣的环境中使用strobe抗多径鉴相器,减少因strobe鉴相器有效鉴相范围较窄而导致的环路错锁的情况发生。
上述的设计方法,其中,所述的多径环境识别模块,使用了累积和(CUSUM)算法,对信号的信噪比和两类鉴相器输出差值的统计分布的变化情况进行监测。
上述的设计方法,其中包括:假设x(n)为n时刻接收机估计接收信号的载噪比,该变量近似服从高斯分布,nv为良好环境变成多径环境的临界时刻,H0为认为当前环境为良好环境的假设,H1为认为当前环境为多径环境的假设,即:
Figure BDA0001581707760000021
Figure BDA0001581707760000022
同时,利用似然比定义变化指标为:
Figure BDA0001581707760000023
则CUSUM根据如下累积和量进行变点判决:
dc(n)=max(0,dc(n-1)+LLR(n))
当某一时刻nv,累积和量大于判决阈值时:
Figure BDA0001581707760000031
则认为统计量分布发生变化,即GNSS接收机运行环境发生改变,
同样,在不同环境下两类鉴相器差值的统计分布特性也会不一样,利用CUSUM算法也可统计相应的累积和变量df(n)进行监测。在本发明中,依据不同的策略,可综合使用这两类累积和变量监测接收机运行环境的变化:
保守判决:
可认为当两项累积和变量均大于相应判决阈值时,才认为环境发生改变。
激进判决:
只要有其中一项累积和变量超过判决阈值,则认为环境发生变化。
其中,所述的组合鉴相器模块,利用了多径环境识别模块以及两类鉴相器输出差值,进行鉴相器的切换。当两类鉴相器输出差值超过某一阈值,则依据多径环境识别模块的输出判断当前运行环境是否为多径环境:
若是多径环境,则不作处理。
若是良好环境,则将锁相环路鉴相器从strobe抗多径鉴相器切换为常规鉴相器。以此避免strobe鉴相环路在常规环境下对与弱信号的码相位错锁的发生。
本发明提供的一种基于CUSUM算法的组合抗多径鉴相器的设计方法具有以下有益效果:1、利用了CUSUM算法对于微小异常变化的敏感性以及常规鉴相器与strobe鉴相器部分优缺点互补的特性,设计了一种实用性更高的组合鉴相器,能够在不中断跟踪环路的情况下,有效地抑制strobe鉴相环路在抗多径的同时又带来的可能的相位错锁情况的发生。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明及其特征、外形和优点将会变得更明显。在全部附图中相同的标记指示相同的部分。并未刻意按照比例绘制附图,重点在于示出本发明的主旨。
图1为本发明提供的一种基于CUSUM算法的组合抗多径鉴相器的设计方法中的跟踪环路结构图;
图2为本发明提供的一种基于CUSUM算法的组合抗多径鉴相器的设计方法中的运行流程图。
具体实施方式
在下文的描述中,给出了大量具体的细节以便提供对本发明更为彻底的理解。然而,对于本领域技术人员而言显而易见的是,本发明可以无需一个或多个这些细节而得以实施。在其他的例子中,为了避免与本发明发生混淆,对于本领域公知的一些技术特征未进行描述。
为了彻底理解本发明,将在下列的描述中提出详细的步骤以及详细的结构,以便阐释本发明的技术方案。本发明的较佳实施例详细描述如下,然而除了这些详细描述外,本发明还可以具有其他实施方式。
参照图1-图2所示,本发明提供了一种基于CUSUM算法的组合抗多径鉴相器的设计方法,包括以下步骤:
步骤(1):利用CUSUM控制图算法的变点理论进行良好环境和多径环境的识别与监测,具体为,以相关值的载噪比CN0以及常规鉴相器和strobe鉴相器的差值为统计量,利用CUSUM算法实时监测这两类统计量指标的变化情况,综合两类指标的变化情况识别当前环境是否处于良好环境;
步骤(2):利用strobe抗多径鉴相器和常规鉴相器的综合使用,来减小strobe鉴相环路在常规环境中错锁的可能性,通过利用多径环境识别模块给出的环境判别结果,并依据两类鉴相器的鉴相差值与阈值的比较情况,进行鉴相器的切换,使得在良好环境中使用常规鉴相器,而在多径恶劣的环境中使用strobe抗多径鉴相器,减少因strobe鉴相器有效鉴相范围较窄而导致的环路错锁的情况发生。
在本发明一优选的实施例中,多径环境识别模块,使用了累积和(CUSUM)算法,对信号的信噪比和两类鉴相器输出差值的统计分布的变化情况进行监测,其中包括:假设x(n)为n时刻接收机估计接收信号的载噪比,该变量近似服从高斯分布,nv为良好环境变成多径环境的临界时刻,H0为认为当前环境为良好环境的假设,H1为认为当前环境为多径环境的假设,即:
Figure BDA0001581707760000051
Figure BDA0001581707760000052
同时,利用似然比定义变化指标为:
Figure BDA0001581707760000053
则CUSUM根据如下累积和量进行变点判决:
dc(n)=max(0,dc(n-1)+LLR(n))
当某一时刻nv,累积和量大于判决阈值时:
Figure BDA0001581707760000054
则认为统计量分布发生变化,即GNSS接收机运行环境发生改变,
同样,在不同环境下两类鉴相器差值的统计分布特性也会不一样,利用CUSUM算法也可统计相应的累积和变量df(n)进行监测。
在本发明的工作流程为:相关值生成器对接收码信号和本地码信号进行相关运算,将生成的相关值输出至常规码鉴相器和strobe鉴相器以及多径环境识别模块;常规码鉴相器和strobe鉴相器利用相关值分别计算得到相应的码相差鉴相结果,并依据鉴相器切换模块选择的结果输出至环路滤波器,同时将和的鉴相器差值输出至鉴相器切换模块以及多径环境识别模块;多径环境识别模块根据相关值以及鉴相差值进行统计计算,将环境识别的结果输出至鉴相器切换模块;鉴相器切换模块依据环境识别结果以及鉴相器差值滤波结果进行鉴相器的选择。后续的环路滤波器和本地码生成器如常规跟踪环路中的部分一样。Strobe鉴相器为主要鉴相器,鉴相器切换模块首先判断当前环境是否为良好环境。若当前运行环境为多径环境,则设置鉴相器为strobe鉴相器。若为良好环境,则进一步判断两类鉴相器差值滤波值是否超过阈值,若超过,则认为strobe鉴相环路发生错锁,将鉴相器切换到常规鉴相器,直到鉴相差值回归到正常水平后再切换到strobe鉴相器。若未超过阈值,依旧保持为strobe鉴相器。
以上对本发明的较佳实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,其中未尽详细描述的设备和结构应该理解为用本领域中的普通方式予以实施;任何熟悉本领域的技术人员,在不脱离本发明技术方案范围情况下,都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰,或修改为等同变化的等效实施例,这并不影响本发明的实质内容。因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。

Claims (3)

1.一种基于CUSUM算法的组合抗多径鉴相器的设计方法,其特征在于,包括以下步骤:
步骤(1):利用CUSUM控制图算法的变点理论进行良好环境和多径环境的识别与监测;
步骤(2):利用strobe抗多径鉴相器和常规鉴相器的综合使用,来减小strobe鉴相环路在常规环境中错锁的可能性,所述步骤(1)中包括:以相关值的载噪比CN0以及常规鉴相器和strobe鉴相器的差值为统计量,利用CUSUM算法实时监测这两类统计量指标的变化情况,综合两类指标的变化情况识别当前环境是否处于良好环境,利用多径环境识别模块给出的环境判别结果,并依据两类鉴相器的鉴相差值与阈值的比较情况,进行鉴相器的切换,使得在良好环境中使用常规鉴相器,而在多径恶劣的环境中使用strobe抗多径鉴相器,减少因strobe鉴相器有效鉴相范围较窄而导致的环路错锁的情况发生。
2.如权利要求1所述的一种基于CUSUM算法的组合抗多径鉴相器的设计方法,其特征在于,所述的多径环境识别模块,使用了累积和CUSUM算法,对信号的信噪比和两类鉴相器输出差值的统计分布的变化情况进行监测。
3.如权利要求2所述的一种基于CUSUM算法的组合抗多径鉴相器的设计方法,其特征在于,其中包括:假设x(n)为n时刻接收机估计接收信号的载噪比,该变量近似服从高斯分布,nv为良好环境变成多径环境的临界时刻,H0为认为当前环境为良好环境的假设,H1为认为当前环境为多径环境的假设,即:
H0
Figure FDA0003810943030000011
H1
Figure FDA0003810943030000012
同时,利用似然比定义变化指标为:
Figure FDA0003810943030000013
则CUSUM根据如下累积和量进行变点判决:
dc(n)=max(0,dc(n-1)+LLR(n))
当某一时刻nv,累积和量大于判决阈值时:
Figure FDA0003810943030000021
则认为统计量分布发生变化,即GNSS接收机运行环境发生改变,
同样,在不同环境下两类鉴相器差值的统计分布特性也会不一样,利用CUSUM算法也可统计相应的累积和变量df(n)进行监测。
CN201810156977.0A 2018-02-24 2018-02-24 一种基于cusum算法的组合抗多径鉴相器的设计方法 Active CN109143276B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810156977.0A CN109143276B (zh) 2018-02-24 2018-02-24 一种基于cusum算法的组合抗多径鉴相器的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810156977.0A CN109143276B (zh) 2018-02-24 2018-02-24 一种基于cusum算法的组合抗多径鉴相器的设计方法

Publications (2)

Publication Number Publication Date
CN109143276A CN109143276A (zh) 2019-01-04
CN109143276B true CN109143276B (zh) 2022-09-30

Family

ID=64801680

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810156977.0A Active CN109143276B (zh) 2018-02-24 2018-02-24 一种基于cusum算法的组合抗多径鉴相器的设计方法

Country Status (1)

Country Link
CN (1) CN109143276B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5808582A (en) * 1996-09-13 1998-09-15 Litton Consulting Group, Inc. Global positioning system receiver with improved multipath signal rejection
US8908744B1 (en) * 2010-02-10 2014-12-09 Marvell International Ltd. Discriminator system for timing error detection in presence and absence of multipath conditions
CN102023299B (zh) * 2010-11-03 2012-11-28 东南大学 一种抑制多径的相关器方法
CN102183770A (zh) * 2011-05-20 2011-09-14 哈尔滨工程大学 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法
CN104849730B (zh) * 2015-06-02 2023-07-04 中国人民解放军国防科技大学 伽利略卫星导航系统e1频点信号无模糊抗多径方法
CN106291606B (zh) * 2016-07-22 2018-12-21 北京航天自动控制研究所 一种抗多径高精度码环鉴相方法

Also Published As

Publication number Publication date
CN109143276A (zh) 2019-01-04

Similar Documents

Publication Publication Date Title
US7957496B2 (en) Suppressing interference for wireless reception and improvements relating to processing a frequency shift keyed signal
EP0945994B1 (en) Adaptive acquisition system for spread spectrum transmissions compensating for frequency offset and noise
CA2113710C (en) Carrier detector
EP1389860A2 (en) Frequency offset estimator
CN109787197B (zh) 基于Bhattacharyya距离算法的线路纵联保护方法
CN109143276B (zh) 一种基于cusum算法的组合抗多径鉴相器的设计方法
US6795452B2 (en) Method of tracking time intervals for a communication signal
US6229862B1 (en) Selective clock recovery from portions of digital data signal best suited therefor
US8416907B2 (en) Method and apparatus for integral state initialization and quality of lock monitoring in a clock and data recovery system
US20050220228A1 (en) Accurate signal detection in a wireless environment
US6529566B1 (en) Method and apparatus for facilitating a classification of a baseband signal
CN109039511A (zh) 一种接收机抗干扰及防错锁的信号处理方法
Ahmad et al. A SARS method for reliable spectrum sensing in multiband communication systems
CN104242928B (zh) 锁频环中一种基于叉积算法的频率锁定检测器处理方法
MXPA05005618A (es) Un detector de bloqueo de circuito de rastreo de portador.
CN101854189B (zh) 一种带有频差估计的软信号识别和扩频码捕捉方法
CN117008160B (zh) Gnss接收机卫星失锁快速判决方法、装置和gnss接收机
Nair et al. An energy detector for cognitive radios in channels at low SNR using adaptive threshold
CN116430415B (zh) 一种自适应多峰值多门限检测的捕获验证策略方法
EP1507377A1 (en) Automatic frequency control with symbol error detection
KR100379394B1 (ko) 디지털 방송 수신기의 락 검출 장치 및 방법
NL2024309B1 (en) Frequency Shift Keying signal demodulator system and method
CN110224961B (zh) 一种检测qpsk信号频率锁定状态的方法及装置
JP4595211B2 (ja) ディジタル復調装置
Ke Adaptive-Threshold Detection Technology

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant