CN105372678B - 一种正弦boc调制信号的无模糊跟踪方法 - Google Patents
一种正弦boc调制信号的无模糊跟踪方法 Download PDFInfo
- Publication number
- CN105372678B CN105372678B CN201510727307.6A CN201510727307A CN105372678B CN 105372678 B CN105372678 B CN 105372678B CN 201510727307 A CN201510727307 A CN 201510727307A CN 105372678 B CN105372678 B CN 105372678B
- Authority
- CN
- China
- Prior art keywords
- signal
- boc
- local
- correlation function
- phase
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于全球导航卫星***(Global Navigation Satell ite System,GNSS)领域,具体涉及一种GNSS卫星导航信号,即二进制偏移载波(Binary Offset Carrier,BOC)调制信号的正弦BOC调制信号的无模糊跟踪方法。本发明的包括:接收BOC信号;将本地产生的脉冲信号分别与本地产生的PRN信号相乘;两路中频信号分别和超前、滞后本地辅助信号进行积分累加,即卷积运算,得到各个支路的积分累加运算结果;鉴相结果ε输入环路滤波器,经过转化,输出调节量,调节PRN。
Description
技术领域
本发明属于全球导航卫星***(Global Navigation Satellite System,GNSS)领域,具体涉及一种GNSS卫星导航信号,即二进制偏移载波(Binary Offset Carrier,BOC)调制信号的正弦BOC调制信号的无模糊跟踪方法。
背景技术
全球卫星导航***(GNSS)以其高精度、全球化、全天候的优质导航服务优势,在军事和民用领域有着广泛应用。为了实现频带共享和频谱分离,达到更高的跟踪精度和性能要求,现代化GPS、Galileo以及中国北斗都采用了BOC信号调制方式。相较于传统的BPSK调制方式,BOC调制具有诸多优点:可通过参数调整实现信号间的频谱分离,使信号能量集中分布在频带边缘;自相关函数主峰更窄,因此在信号跟踪过程中有更强的抗热噪声、抗窄带干扰以及抗多径能力,从而提高卫星导航信号接收机的测量精度。
然而,BOC信号的自相关函数存在多峰性,并且随着调制阶数的增加,副峰的数目以及峰值都会随之而增加。一方面增加了误捕概率,另一方面,在信号跟踪环路的鉴相输出曲线出现多个误锁点,使导航接收机易出现误捕和误锁,即模糊问题。
针对BOC信号的模糊问题,目前比较典型的方法主要包括Bump-Jumping方法、BPSK-like方法以及基于伪相关函数(Pseudo Correlation Function,PCF)的PUDLL(thePseudo-correlation-function based Unambiguous Delay Lock Loop,PUDLL)处理法。Bump-Jumping方法(即峰跳法)是通过增加远超前和远滞后两路相关器,来实时监测当前锁定的相关峰和其两侧相关峰的幅度差异来判断误锁是否发生。但是硬件结构复杂,适用于低阶BOC调制信号,对于信噪比和接收机带宽较为敏感,跟踪结果受多径影响较大。BPSK-like方法主要采用滤波来消除跟踪模糊性,基本思想是将BOC调制信号的两个边带当作两个BPSK信号分别处理,从而去除信号的多相关峰,也称作边带处理法。使用这种方法时接收机中需要滤波装置,增加了硬件复杂度。此外,滤波后BOC调制信号的相关峰被展宽,导致***的跟踪精度及抗多径能力都有所降低。伪相关函数(PCF)处理法则是通过使用不同的辅助信号与BOC调制信号进行组合运算,得到无边峰的伪相关函数。伪相关函数法存在能量损失,在调制阶数增加时尤其严重,容易造成误捕和误锁的发生。
发明内容
针对现有BOC信号模糊消除方法的不足,本发明的目的在于提出了一种完全移除了相关函数的边峰,消除了BOC信号跟踪时模糊问题的正弦BOC调制信号的无模糊跟踪方法。
本发明的目的是这样实现的:
(1)接收BOC信号r(t):
其中P代表信号能量;x(t)是基带sine-BOC信号;δ是以码片为单位的信号传输延时;ω和ψ分别为信号载波的角频率和相位;将接收信号r(t)与接收端的本地载波相乘,实现载波剥离,经过下变频得到同相I路和正交Q路两路中频信号:
(2)将本地产生的脉冲信号g1(t)和g2(t)分别与本地产生的PRN信号相乘,得到本地可调脉冲辅助信号x1(t)和x2(t);根据选择的延迟间隔Δ,对x1(t)和x2(t)经行超前、滞后处理,得到相应的超前滞后支路信号x1(t-Δ/2)、x1(t+Δ/2)、x2(t-Δ/2)和x2(t+Δ/2);
(3)两路中频信号分别和超前、滞后本地辅助信号进行积分累加,即卷积运算,得到各个支路的积分累加运算结果:
其中:i=1,2分别代表本地辅助信号x1(t)和x2(t);j=E,L对应超前、滞后相关器;δj代表超前、滞后支路对应的码片延迟,δE=Δ/2,δL=-Δ/2;εδ为以码片为单位的信号传输延时估计误差;
通过CCF生成器经行设定的绝对值运算:
输出结果进行超前减滞后能量鉴相运算,得到鉴相结果:
(4)鉴相结果ε输入环路滤波器,经过转化,输出调节量,调节PRN生成器和g1(t)、g2(t)生成器的相位,直至其与接收BOC信号相位一致,完成信号跟踪过程。
本发明的有益效果在于:
本发明运算复杂度适中,有很好的抗热噪声能力,并且对短时延多径信号有较好的抑制效果。
附图说明
图1是sine-BOC(1,1)的CCF示意图。
图2是不同λ的sine-BOC(1,1)、sine-BOC(2,1)和sine-BOC(5,2)的CCF。
图3是基于CCF的DLL结构图。
图4是本发明与常规方法的抗热噪声性能比较图。
图5是本发明与常规方法的抗多径性能比较图。
具体实施方式
下面结合附图对本发明的实施步骤作进一步的描述:
本发明基于接收机端的本地辅助信号提出了一种针对正弦BOC调制信号的无模糊跟踪方法。该方法首先在接收机端产生两个互为镜像的本地可调脉冲辅助信号x1和x2,然后利用这两种本地可调脉冲辅助信号分别与接收到的BOC信号进行互相关运算,再将得到的两个互相关运算结果进行指定的绝对值运算,从而得到无边峰的合成相关函数(ComposedCorrelation Function,CCF),使鉴相器输出的鉴相曲线中只存在唯一一个过零点,实现BOC信号的无模糊跟踪。本发明运算复杂度适中,有很好的抗热噪声能力,并且对短时延多径信号有较好的抑制效果。
本发明的关键在于本地可调脉冲辅助信号x1和x2的设计以及无边峰的合成相关函数的获得算法。
一般化的脉冲幅值调制信号可以表示为:
其中cl是一系列相互独立的随机变量(l=...,-1,0,1,...为随机变量编号),且当l≠m时,E{clcm}=0;g(t)是定义在区间0≤t≤Tc上的脉冲(Tc为单个随机变量的持续时间),在区间外,g(t)=0。
因此,sine-BOC(m,n)信号x(t)可以用以上形式表示,其中g(t)=sign{sin(2πfst)},0≤t≤Tc。fs为亚载波频率,且fs=m×1.023MHz;fc为伪码频率,且fc=n×1.023MHz,Tc=1/fc为码片宽度。BOC(m,n)信号调制阶数M为M=2m/n。
本地可调脉冲辅助信号x1(t)和x2(t)设计为以下形式(i=1,2):
其中,λ为脉冲宽度调节参数,M为BOC信号的阶数。
基带sine-BOC(m,n)信号x(t)与本地可调脉冲辅助信号的互相关函数,以x1(t)为例,可以表示为:
其中:τ为两个信号的相关间隔。经过进一步的计算化简,可以得到BOC信号x(t)与本地可调脉冲辅助信号x1(t)的互相关函数为:
其中,k=1,...,M-1。如上式所示,互相关函数是分段线性的。由于本地可调脉冲辅助信号x2(t)是x1(t)的镜像信号,因此有为了获得无边峰的合成相关函数R(τ),采用以下的绝对值运算:
为了完全消除边峰,并且保持区间线性性质,脉冲宽度调节参数的取值范围λ<1/2M。忽略接收机前端滤波时,在[-λ,λ]区间上的合成相关函数为:
可以看出,合成相关函数完全消除了边峰,使得鉴相输出曲线不存在误锁点,解决了模糊问题。
步骤1:接收BOC信号r(t)可以表示为
其中P代表信号能量;x(t)是基带sine-BOC信号;δ是以码片为单位的信号传输延时;ω和ψ分别为信号载波的角频率和相位。将接收信号r(t)与接收端的本地载波相乘,实现载波剥离,经过下变频得到同相(I路)和正交(Q路)两路中频信号:
步骤2:将本地产生的脉冲信号g1(t)和g2(t)分别与本地产生的PRN信号相乘,得到本地可调脉冲辅助信号x1(t)和x2(t)。根据选择的延迟间隔Δ,对x1(t)和x2(t)经行超前、滞后处理,得到相应的超前滞后支路信号x1(t-Δ/2)、x1(t+Δ/2)、x2(t-Δ/2)和x2(t+Δ/2)。
步骤3:两路中频信号分别和超前、滞后本地辅助信号进行积分累加,即卷积运算,得到各个支路的积分累加运算结果:
其中:i=1,2分别代表本地辅助信号x1(t)和x2(t);j=E,L对应超前、滞后相关器;δj代表超前、滞后支路对应的码片延迟,δE=Δ/2,δL=-Δ/2;εδ为以码片为单位的信号传输延时估计误差。
通过CCF生成器经行设定的绝对值运算:
输出结果进行超前减滞后能量鉴相运算,得到鉴相结果:
步骤4:鉴相结果ε输入环路滤波器,经过转化,输出调节量,调节PRN生成器和g1(t)、g2(t)生成器的相位,直至其与接收BOC信号相位一致,完成信号跟踪过程。
性能分析
结合附图对本发明的性能进行分析。
(1)抗热噪声性能分析
仅考虑热噪声干扰的情况下,在跟踪环路中采用非相干鉴相器,同时忽略数据调制。因此,接收BOC信号可以表示为:
其中,P代表信号能量;x(t)是基带sine-BOC信号;δ是以码片为单位的信号传输延时;ω和ψ分别为信号载波的角频率和相位;n(t)是单边功率谱密度为2N0的附加高斯白噪声。
为了分析简便,忽略多普勒频移,并且假设载波完全剥离。此时,积分累加器的输出:
i=1,2代表本地可调脉冲辅助信号x1(t)和x2(t),j=E,L则分别代表超前和滞后相关器。δj是相应的信号延迟(δE=Δ/2,δL=-Δ/2)。εδ是对信号传输延迟δ的估计误差。是对本地可调脉冲辅助信号与热噪声的相关计算输出结果:
由于产生于同一个噪声过程,和为相关高斯变量,服从以下联合正太分布:
其中
鉴相器输出为:
近似码跟踪误差标准差为:
其中,BL为环路滤波器带宽,T为想干积分时间,G为鉴相器增益,σε为鉴相器输出的标准差。
忽略前端滤波,通过蒙特卡洛仿真(105次),可以得到近似伪码跟踪误差标准差如图4所示。在图4中,为了便于进行性能对比,给出了sine-BOC(1,1)信号、ASPeCT算法(β=1)和PUDLL(k=0)算法的伪码跟踪误差标准差。其中,环路滤波器带宽BL=2Hz,相干积分时间T=1ms,超前滞后间隔Δ=0.1个码片。
从图4可以看出,采用本发明算法时的伪码跟踪误差标准差明显小于sine-BOC(1,1)信号、ASPeCT算法(β=1)和PUDLL(k=0)算法,尤其是在低于35dB的中低信噪比区域,即本发明的抗热噪声性能明显优于其他算法。
(2)抗多径性能分析
多径信号会导致鉴相器输出曲线的零点偏移,即多径误差。多径误差的大小主要取决于多径信号相当于直达信号的传输时延、相对幅值以及相位偏移。为了简化,仅考虑存在一路多径信号的情况。
图5给出了以sine-BOC(1,1)为例的多径误差包络对比图。其中多径信号幅值较直达信号小5dB,超前滞后间隔为0.25码片。为了便于进行性能对比,给出了sine-BOC(1,1)信号、ASPeCT算法(β=1)和PUDLL(k=0)算法的多径误差包络。通过对比可以看出,本发明对中短时延多径误差的抑制效果比较好,小于ASPeCT算法以及PUDLL算法的短时延多径误差,但是略大于sine-BOC(1,1)信号本身的短时延多径误差;本发明的长时延多径误差则明显大于sine-BOC(1,1)信号、ASPeCT算法以及PUDLL算法。另外,随着脉冲宽度λ的减小,本发明多径抑制效果变好。
Claims (1)
1.一种正弦BOC调制信号的无模糊跟踪方法,其特征在于包括如下步骤:
本地可调脉冲辅助信号x1(t)和x2(t)设计为以下形式i=1,2:
其中cl是一系列相互独立的随机变量,l=...,-1,0,1,...为随机变量编号,且当l≠m时,E{clcm}=0,Tc为单个随机变量的持续时间,g(t)是定义在区间0≤t≤Tc上的脉冲,λ为脉冲宽度调节参数,M为BOC信号的阶数;
基带sine-BOC(m,n)信号x(t)与本地可调脉冲辅助信号的互相关函数,表示为:
其中:τ为两个信号的相关间隔;经过进一步的计算化简,可以得到BOC信号x(t)与本地可调脉冲辅助信号x1(t)的互相关函数为:
其中,k=1,...,M-1;如上式所示,互相关函数是分段线性的;由于本地可调脉冲辅助信号x2(t)是x1(t)的镜像信号,有为了获得无边峰的合成相关函数R(τ),采用以下的绝对值运算:
为了完全消除边峰,并且保持区间线性性质,脉冲宽度调节参数的取值范围λ<1/2M;
忽略接收机前端滤波时,在[-λ,λ]区间上的合成相关函数为:
(1)接收BOC信号r(t):
其中P代表信号能量;x(t)是基带sine-BOC信号;δ是以码片为单位的信号传输延时;ω和ψ分别为信号载波的角频率和相位;将接收信号r(t)与接收端的本地载波相乘,实现载波剥离,经过下变频得到同相I路和正交Q路两路中频信号:
(2)将本地产生的脉冲信号g1(t)和g2(t)分别与本地产生的PRN信号相乘,得到本地可调脉冲辅助信号x1(t)和x2(t);根据选择的延迟间隔Δ,对x1(t)和x2(t)经行超前、滞后处理,得到相应的超前滞后支路信号x1(t-Δ/2)、x1(t+Δ/2)、x2(t-Δ/2)和x2(t+Δ/2);
(3)两路中频信号分别和超前、滞后本地辅助信号进行积分累加,即卷积运算,得到各个支路的积分累加运算结果:
其中:i=1,2分别代表本地辅助信号x1(t)和x2(t);j=E,L对应超前、滞后相关器;δj代表超前、滞后支路对应的码片延迟,δE=Δ/2,δL=-Δ/2;εδ为以码片为单位的信号传输延时估计误差;
通过CCF生成器经行设定的绝对值运算:
输出结果进行超前减滞后能量鉴相运算,得到鉴相结果:
(4)鉴相结果ε输入环路滤波器,经过转化,输出调节量,调节PRN生成器和g1(t)、g2(t)生成器的相位,直至其与接收BOC信号相位一致,完成信号跟踪过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510727307.6A CN105372678B (zh) | 2015-10-30 | 2015-10-30 | 一种正弦boc调制信号的无模糊跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510727307.6A CN105372678B (zh) | 2015-10-30 | 2015-10-30 | 一种正弦boc调制信号的无模糊跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105372678A CN105372678A (zh) | 2016-03-02 |
CN105372678B true CN105372678B (zh) | 2019-04-19 |
Family
ID=55375028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510727307.6A Active CN105372678B (zh) | 2015-10-30 | 2015-10-30 | 一种正弦boc调制信号的无模糊跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105372678B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106383083B (zh) * | 2016-11-01 | 2018-10-16 | 北京信息科技大学 | 流式细胞仪光谱重叠信号数字处理方法 |
CN108490462B (zh) * | 2018-03-21 | 2021-02-19 | 桂林电子科技大学 | 基于相关函数重构的boc无模糊跟踪方法 |
CN110109154B (zh) * | 2019-05-16 | 2021-02-19 | 桂林电子科技大学 | 一种BOC(n,n)移位相关无模糊捕获方法及装置 |
CN115657091B (zh) * | 2022-12-29 | 2023-03-31 | 成都奇芯微电子有限公司 | 一种高动态gnss跟踪的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424754A (zh) * | 2013-07-18 | 2013-12-04 | 哈尔滨工程大学 | 基于合成相关函数的mboc调制信号无模糊多径抑制方法 |
CN103645484A (zh) * | 2013-11-29 | 2014-03-19 | 哈尔滨工程大学 | 一种复合二进制偏移载波调制信号的处理方法 |
CN104898143A (zh) * | 2015-06-02 | 2015-09-09 | 刘哲 | 时分复用双参考波形的tmboc信号无模糊接收方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2832806B1 (fr) * | 2001-11-23 | 2006-07-21 | Thales Sa | Procede et dispositif de calcul de la fonction discriminante de signaux modules avec une ou plusieurs sous-porteuse |
KR101381104B1 (ko) * | 2013-05-22 | 2014-04-04 | 성균관대학교산학협력단 | Cboc 상관함수 생성 방법, cboc 신호 추적 방법 및 cboc 신호 추적 시스템 |
-
2015
- 2015-10-30 CN CN201510727307.6A patent/CN105372678B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424754A (zh) * | 2013-07-18 | 2013-12-04 | 哈尔滨工程大学 | 基于合成相关函数的mboc调制信号无模糊多径抑制方法 |
CN103645484A (zh) * | 2013-11-29 | 2014-03-19 | 哈尔滨工程大学 | 一种复合二进制偏移载波调制信号的处理方法 |
CN104898143A (zh) * | 2015-06-02 | 2015-09-09 | 刘哲 | 时分复用双参考波形的tmboc信号无模糊接收方法 |
Non-Patent Citations (2)
Title |
---|
BOC信号合成相关函数通用无模糊跟踪方法;任嘉伟等;《航空学报》;20140725;第35卷(第7期);第2031-2040页 |
Correlation Combination Ambiguity Removing Technology for Acquisition of Sine-Phased BOC Signals;LIU Wen et.al;《COMMUNICATIONS SYSTEM DESIGN》;20150430;第86-96页 |
Also Published As
Publication number | Publication date |
---|---|
CN105372678A (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101132191B (zh) | 一种gnss接收机基带信号处理方法 | |
CN104407323B (zh) | 一种高动态低信噪比扩频信号伪码时延测量方法 | |
CN105917622B (zh) | 用于接收复合信号的方法和接收器 | |
CN106772455B (zh) | 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法 | |
CN106291614B (zh) | 用于在多路径环境中跟踪卫星无线电导航信号的装置 | |
KR101467348B1 (ko) | 부분상관함수에 기초한 tmboc(6,1,4/33) 신호를 위한 비모호 상관함수 생성 방법, tmboc 신호 추적 장치 및 이를 이용한 위성 항법 신호 수신 시스템 | |
CN102944884B (zh) | Gnss接收机检测并消除窄带干扰的方法 | |
CN105372678B (zh) | 一种正弦boc调制信号的无模糊跟踪方法 | |
CN103424754B (zh) | 基于合成相关函数的mboc调制信号无模糊多径抑制方法 | |
WO2018103186A1 (zh) | 一种TD-AltBOC信号的接收方法和装置 | |
CN108897009B (zh) | 一种boc导航信号接收机及其码跟踪方法 | |
CN106918822B (zh) | 计算用于解析副载波跟踪模糊度的非模糊鉴别器的gnss接收器 | |
TW201445167A (zh) | Gps接收機及判斷gps接收機跟踪環路狀態的方法 | |
CN105445767B (zh) | 基于平均模糊函数的boc信号参数盲估计方法 | |
CN102023299B (zh) | 一种抑制多径的相关器方法 | |
KR101467320B1 (ko) | 균등 세분된 부분상관함수에 기초한 tmboc(6,1,4/33) 신호를 위한 비모호 상관함수 생성 방법, tmboc 신호 추적 장치 및 이를 이용한 위성 항법 신호 수신 시스템 | |
CN103760578B (zh) | 一种gnss卫星导航信号的无模糊跟踪方法 | |
KR101467234B1 (ko) | 부분상관함수들의 단계적 조합에 기초한 cboc(6,1,1/11) 신호를 위한 비모호 상관함수 생성 방법, cboc 신호 추적 장치 및 이를 이용한 위성 항법 신호 수신 시스템 | |
CN108469623A (zh) | 消除boc调制信号自相关函数边峰的方法 | |
KR101381104B1 (ko) | Cboc 상관함수 생성 방법, cboc 신호 추적 방법 및 cboc 신호 추적 시스템 | |
JP4700518B2 (ja) | 同期タイミング検出装置および受信機 | |
Chen et al. | Research and implementation of beidou-3 satellite multi-band signal acquisition and tracking method | |
KR101812323B1 (ko) | 사인 위상 boc 신호를 추적하는 방법 및 사인 위상 boc 신호 추적 장치 | |
JP5933559B2 (ja) | 信号捕捉装置及び方法 | |
Mao et al. | Unambiguous BPSK-like CSC method for Galileo acquisition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |