CN109633713B - 一种北斗导航弱信号接收的频率跟踪方法 - Google Patents
一种北斗导航弱信号接收的频率跟踪方法 Download PDFInfo
- Publication number
- CN109633713B CN109633713B CN201811626011.5A CN201811626011A CN109633713B CN 109633713 B CN109633713 B CN 109633713B CN 201811626011 A CN201811626011 A CN 201811626011A CN 109633713 B CN109633713 B CN 109633713B
- Authority
- CN
- China
- Prior art keywords
- frequency
- discrimination
- frequency error
- signal
- data
- 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
Images
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
-
- 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/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/254—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to Doppler shift of satellite signals
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种北斗导航弱信号接收的频率跟踪方法,包括以下步骤:S1、信号经伪码解扩进入相关器;S2、根据累加运算结果,进行电文自预测;S3、使用预测值对接收信号进行电文剥离;S4、对电文剥离后的数据进行缓存;S5、得到数据的粗鉴频的频率误差值df1;S6、得到数据的细鉴频的频率误差df2;S7、将粗鉴频的频率误差值df1和细鉴频的频率误差值df2进行相加得到总频率误差值;S8、将总频率误差信号送入环路滤波器进行滤波;S9、环路滤波器输出信号驱动数控振荡器,进行频率跟踪。本发明通过采用进行过频率误差补偿的超前滞后非相干鉴频,频率跟踪环路更能细致地实施频率跟踪,提高了环路的稳定性。
Description
技术领域
本发明属于卫星导航领域,特别是涉及一种北斗导航弱信号接收的载波频率跟踪方法。
背景技术
卫星导航系统已经成为当今世界国家综合国力及科学技术发展水平的重要标志,是国家经济基础的重要组成部分,它已经进入人们的日常生活中,与社会发展和经济建设息息相关。我国也已经投入大量的资金和人力,积极进行北斗卫星导航系统(BeiDouNavigation Satellite System,BDS)的研究。到2020年左右,我国将建成覆盖全球的北斗卫星导航系统,其将成为继美国GPS、俄罗斯GLONASS和欧洲Galileo之后的又一个全球全天候的卫星导航系统。
随着卫星导航应用的进一步发展,要求接收机能工作于更为严苛的应用环境下,例如,城市峡谷、林荫大道等卫星信号严重遮挡环境。所谓的高灵敏度接收机(HS-GNSS,High Sensitivity-Global Navigation Satellite Systems),就是为满足此类对弱信号接收的需求。载波跟踪环路是弱信号接收的关键环节,其性能对高灵敏度接收机具有重要影响。通常,通过加长相干积分时间,来抵抗弱信号中噪声的影响,以提高载波跟踪环路性能。基于FFT鉴频的载波频率跟踪,利用了FFT频谱分析,得到频率误差信号。同时,FFT变换也是一个相干积分过程。然而,受限于FFT实现的复杂度以及频率跟踪环路迟延,FFT鉴频算法的相干积分时间不能太长。另一方面,较短的相干积分时间会带来FFT频谱分析分辨率低的问题。
发明内容
本发明提供一种北斗导航弱信号接收的频率跟踪方法,可用于弥补FFT鉴频器频率分辨率低的不足,改善频率跟踪环路的稳定性。
本发明提出的一种北斗导航弱信号接收的频率跟踪方法,包括以下步骤:
S1、接收信号经过伪码解扩以及载波剥离后,进入相关器进行累加运算,取得相关值;
S2、根据累加运算结果,进行电文自预测,通过电文自预测算法估计得出卫星发送信号中的电文数据比特,取得电文剥离所必需的电文比特;
S3、使用步骤S2中自预测算法估计得出的电文数据比特,对接收信号进行电文剥离,以去除电文数据调制对相干积分的影响;
S4、对电文剥离后的数据进行缓存,缓存到预定的N个数据;
S5、对所述N个数据进行FFT鉴频,得到粗鉴频的频率误差值df1;
S6、对所述N个数据进行超前滞后非相干鉴频,得到细鉴频的频率误差值df2;
S7、将粗鉴频的频率误差值df1和细鉴频的频率误差值df2进行相加得到总频率误差值;
S8、将所述总频率误差值送入环路滤波器进行滤波;
S9、环路滤波器输出信号驱动数控振荡器,进行频率跟踪。
优选地,S5包括以下步骤:S5.1、对输入的所述N个数据进行FFT变换;S5.2、对FFT变换输出进行模平方运算;S5.3、从模平方后的N个数据中搜索最大值,并得到最大值所对应的索引idx,其中,当索引idx>N/2时,idx=idx–N;S5.4、计算粗鉴频的频率误差df1=idx/(Tcoh*N),其中,Tcoh表示相关积分时间。
优选地,S6包括以下步骤:S6.1、利用上述的频率误差信号df1和超前滞后非相干鉴频器的频率鉴别间隔fsh,进行频率搬移;S6.2、对频率搬移后的信号进行相干累加得到Xl=∑Zl(i),Xr=∑Zr(i),i=1,2,…N;其中,Z1表示超前支路进行频率搬移后的信号值,Zr表示滞后支路进行频率搬移后的信号值;S6.3、累加后的信号再进行模平方El=|Xl|2,Er=|Xr|2;S6.4、对模平方输出信号进行频率误差计算,得到细鉴频的频率误差df2。
优选地,S6.1中超前支路进行频率搬移后的公式为:Zl=Y*exp(-j2pi(df1–fsh));滞后支路进行频率搬移后的公式为:Zr=Y*exp(-j2pi(df1+fsh)),其中,Y=Yi+jYq为相关器输出的复数信号,pi为圆周率,fsh为超前滞后非相干鉴频器的频率鉴别间隔。
优选地,S6.1中频率鉴别间隔fsh满足条件fsh≤1/(Tcoh*N)。
优选地,S6.4中细鉴频的超前滞后非相干鉴频的频率误差计算公式为:
一种计算机可读存储介质,其存储有与计算设备结合使用的计算机程序,所述计算机程序可被处理器执行以上述任一项所述方法。
本发明的有益效果包括:通过采用进行过频率误差补偿的超前滞后非相干鉴频,频率跟踪环路更能细致地实施频率跟踪,提高了环路的稳定性。
附图说明
图1是本发明实施例的载波频率跟踪环路框图。
图2是本发明实施例的FFT鉴频器的处理流程图。
图3是本发明实施例的超前滞后非相干鉴频器的处理流程图。
图4是本发明实施例的频率跟踪环路的多普勒频率仿真结果示意图。
具体实施方式
下面结合具体实施方式并对照附图对本发明作进一步详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
如图1所示,为本发明实施例的的载波频率跟踪环路框图。接收信号经过伪码解扩以及载波剥离后,送入相关器,相关积分时间(Tcoh)为一个电文比特,即相关积分时间Tcoh=20ms。
相关器输出的相关值,送入电文自预测,估计出卫星发送信号中的电文数据比特Dn。使用自预测算法估计得出的电文数据比特Dn对接收信号进行电文剥离,以去除电文数据调制对后续的相干积分的影响。电文自预测算法,可以采用已知的最大似然估计(Maximum likelihood,ML)的算法,或者差分解调估计算法,这里不作限定。
电文剥离后的数据送入缓存器,缓存到预定的N个数据后,送入FFT鉴频器和超前滞后非相干鉴频器。
如图2所示为FFT鉴频器的处理流程图,包括下列步骤(下述步骤中N代表上述预定的数据的个数):
对输入的N个采样数据进行FFT变换;
对FFT变换输出进行模平方运算;
从模平方后的N个数据中搜索最大值,并得到最大值所对应的索引idx;
当且仅当索引idx>N/2时,idx=idx–N;
计算粗鉴频的频率误差df1=idx/(Tcoh*N)。
粗鉴频的频率误差df1送入超前滞后非相干鉴频器。利用df1,非相干鉴频器可以对输入信号中的频率误差进行初略的校正,以便进行细鉴频。
超前滞后非相干鉴频器的处理流程如图3所示,包括下列步骤:
超前支路进行频率搬移,得到搬移后的信号值Z1,如公式(1):
Zl=Y*exp(-j2pi(df1–fsh)) (1)
式中,pi为圆周率,Y为相关器输出的复数信号,Y=Yi+jYq,fsh为超前滞后非相干鉴频器的频率鉴别间隔,其必须满足条件fsh≤1/(Tcoh*N)。间隔越小,鉴频输出的线性度越好,优选地,fsh可选择为1/(Tcoh*N*8)。
滞后支路进行频率搬移,得到搬移后的信号值Zr,如公式(2),表示滞后支路进行频率搬移后的信号值:
Zr=Y*exp(-j2pi(df1+fsh)) (2)
对频率搬移后的信号进行相干累加,如公式(3):
Xl=∑Zl(i),i=1,2,…N
Xr=∑Zr(i),i=1,2,…N (3)
公式(3)中,X1表示总的超前支路进行频率搬移后的信号值之和,Xr表示总的滞后支路进展频率搬移后的信号值之和。
对上述进行相干累加后的信号再进行模平方,如公式(4):
El=|Xl|2
Er=|Xr|2 (4)
最后,进行频率误差计算,得到细鉴频的频率误差df2,频率误差计算公式为:
FFT鉴频得到的粗鉴频信号df1和非相干鉴频器得到的细鉴频信号df2相加,得到总的载波频率误差信号。频率误差信号送入环路滤波器进行滤波后,驱动数控振荡器(NCO,numerically controlled oscillator),完成载波频率跟踪。
如图4所示的频率跟踪环路的多普勒频率仿真结果示意图,图中的1代表仿真参考,该仿真曲线在10s后没有加入噪声,故多普勒频率非常稳定;2代表本发明实施例的仿真结果;3代表采用FFT鉴频算法的仿真结果。从图中可以看出,本发明的频率误差补偿的超前滞后非相干鉴频比传统的FFT鉴频更为稳定。
本发明还提供一种计算机可读存储介质,其存储有与计算设备结合使用的计算机程序,该计算机程序可被处理器执行以实现上述任一种所述方法。
本领域技术人员将认识到,对以上描述做出众多变通是可能的,所以实施例和附图仅是用来描述一个或多个特定实施方式。
尽管已经描述和叙述了被看作本发明的示范实施例,本领域技术人员将会明白,可以对其做出各种改变和替换,而不会脱离本发明的精神。另外,可以做出许多修改以将特定情况适配到本发明的教义,而不会脱离在此描述的本发明中心概念。所以,本发明不受限于在此披露的特定实施例,但本发明可能还包括属于本发明范围的所有实施例及其等同物。
Claims (7)
1.一种北斗导航弱信号接收的频率跟踪方法,其特征在于,包括以下步骤:
S1、接收信号经过伪码解扩以及载波剥离后,进入相关器进行累加运算,取得相关值;
S2、根据累加运算结果,进行电文自预测,通过电文自预测算法估计得出卫星发送信号中的电文数据比特;
S3、使用步骤S2中电文自预测算法估计得出的电文数据比特,对接收信号进行电文剥离,以去除电文数据调制对相干积分的影响;
S4、对电文剥离后的数据进行缓存,缓存到预定的N个数据;
S5、对所述N个数据进行FFT鉴频,得到粗鉴频的频率误差值df1;
S6、对所述N个数据进行超前和滞后非相干鉴频,得到细鉴频的频率误差值df2;
S7、将粗鉴频的频率误差值df1和细鉴频的频率误差值df2进行相加得到总频率误差值;
S8、将所述总频率误差值送入环路滤波器进行滤波;
S9、环路滤波器输出信号驱动数控振荡器,进行频率跟踪。
2.如权利要求1所述的北斗导航弱信号接收的频率跟踪方法,其特征在于,S5包括以下步骤:
S5.1、对输入的所述N个数据进行FFT变换;
S5.2、对FFT变换输出进行模平方运算;
S5.3、从模平方后的N个数据中搜索最大值,并得到最大值所对应的索引idx,其中,当索引idx>N/2时,idx=idx–N;
S5.4、计算粗鉴频的频率误差值df1=idx/(Tcoh*N),其中,Tcoh表示相关积分时间。
3.如权利要求1所述的北斗导航弱信号接收的频率跟踪方法,其特征在于,S6包括以下步骤:
S6.1、利用粗鉴频的频率误差值df1和超前滞后非相干鉴频器的频率鉴别间隔fsh,进行频率搬移;
S6.2、对频率搬移后的信号进行相干累加得到Xl=∑Zl(i),Xr=∑Zr(i),i=1,2,…N;其中,Z1表示超前支路进行频率搬移后的信号值,Zr表示滞后支路进行频率搬移后的信号值;
S6.3、对S6.2中进行相干累加后的信号值X1和Xr再分别进行模平方,即El=|Xl|2,Er=|Xr|2;
S6.4、对模平方输出信号进行频率误差计算,得到细鉴频的频率误差值df2。
4.如权利要求3所述的北斗导航弱信号接收的频率跟踪方法,其特征在于:S6.1中超前支路进行频率搬移后的公式为:Zl=Y*exp(-j2pi(df1–fsh)),滞后支路进行频率搬移的公式为:Zr=Y*exp(-j2pi(df1+fsh)),其中,Y=Yi+j Yq为相关器输出的复数信号,pi为圆周率,fsh为超前滞后非相干鉴频器的频率鉴别间隔。
5.如权利要求3所述的北斗导航弱信号接收的频率跟踪方法,其特征在于:S6.1中频率鉴别间隔fsh满足条件fsh≤1/(Tcoh*N)。
7.一种计算机可读存储介质,其存储有与计算设备结合使用的计算机程序,所述计算机程序可被处理器执行以实现权利要求1至6任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811626011.5A CN109633713B (zh) | 2018-12-28 | 2018-12-28 | 一种北斗导航弱信号接收的频率跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811626011.5A CN109633713B (zh) | 2018-12-28 | 2018-12-28 | 一种北斗导航弱信号接收的频率跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633713A CN109633713A (zh) | 2019-04-16 |
CN109633713B true CN109633713B (zh) | 2020-06-23 |
Family
ID=66078765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811626011.5A Active CN109633713B (zh) | 2018-12-28 | 2018-12-28 | 一种北斗导航弱信号接收的频率跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633713B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112596082B (zh) * | 2020-11-16 | 2023-12-26 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 北斗b3i信号的载波频率牵引方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495417A (zh) * | 2011-12-07 | 2012-06-13 | 无锡德思普科技有限公司 | 全球导航卫星系统接收机的高鲁棒性载波跟踪系统及方法 |
EP2806609A1 (en) * | 2013-05-24 | 2014-11-26 | O2 Micro, Inc. | Receivers and methods for detecting false lock of a carrier tracking loop |
CN104765049A (zh) * | 2015-04-20 | 2015-07-08 | 和芯星通科技(北京)有限公司 | 一种导航卫星信号跟踪方法及装置 |
CN105425257A (zh) * | 2015-11-03 | 2016-03-23 | 上海交通大学 | 一种高动态gnss载波信号的跟踪方法及系统 |
CN105572701A (zh) * | 2015-12-17 | 2016-05-11 | 深圳市力合微电子股份有限公司 | 一种用于北斗二代接收机弱信号跟踪的码环鉴相方法 |
CN106019333A (zh) * | 2016-05-16 | 2016-10-12 | 中国人民解放军国防科学技术大学 | 一种基于非相干鉴别器的北斗导航信号矢量跟踪方法 |
-
2018
- 2018-12-28 CN CN201811626011.5A patent/CN109633713B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495417A (zh) * | 2011-12-07 | 2012-06-13 | 无锡德思普科技有限公司 | 全球导航卫星系统接收机的高鲁棒性载波跟踪系统及方法 |
EP2806609A1 (en) * | 2013-05-24 | 2014-11-26 | O2 Micro, Inc. | Receivers and methods for detecting false lock of a carrier tracking loop |
CN104765049A (zh) * | 2015-04-20 | 2015-07-08 | 和芯星通科技(北京)有限公司 | 一种导航卫星信号跟踪方法及装置 |
CN105425257A (zh) * | 2015-11-03 | 2016-03-23 | 上海交通大学 | 一种高动态gnss载波信号的跟踪方法及系统 |
CN105572701A (zh) * | 2015-12-17 | 2016-05-11 | 深圳市力合微电子股份有限公司 | 一种用于北斗二代接收机弱信号跟踪的码环鉴相方法 |
CN106019333A (zh) * | 2016-05-16 | 2016-10-12 | 中国人民解放军国防科学技术大学 | 一种基于非相干鉴别器的北斗导航信号矢量跟踪方法 |
Non-Patent Citations (2)
Title |
---|
一种改善GNSS弱信号动态跟踪性能的FFT鉴频方法;祁发瑞 等;《航空学报》;20180522;第39卷(第8期);321932-1至7 * |
卫星导航接收机载波捕获与跟踪环路关键技术研究;韩吉南;《中国博士学位论文全文数据库 信息科技辑》;20180615;I136-115 * |
Also Published As
Publication number | Publication date |
---|---|
CN109633713A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20120326926A1 (en) | High sensitivity gps/gnss receiver | |
CN101807940B (zh) | 一种gnss接收系统的抗多径干扰装置及其方法 | |
Yang et al. | Generalized multifrequency GPS carrier tracking architecture: Design and performance analysis | |
CN103116169A (zh) | 一种基于矢量跟踪环路的抗干扰方法 | |
CN103592662A (zh) | 一种gps信号接收机的载波跟踪方法及环路 | |
CN102565825B (zh) | 接收信号可靠度判定装置、方法及码相位误差算出方法 | |
Wang et al. | A scheme for weak GPS signal acquisition aided by SINS information | |
Tamazin et al. | Robust fine acquisition algorithm for GPS receiver with limited resources | |
CN109633713B (zh) | 一种北斗导航弱信号接收的频率跟踪方法 | |
JP4777353B2 (ja) | Gps測位方法及びgps測位装置 | |
US8494094B2 (en) | Demodulation of data collected prior to bit edge detection | |
Fan et al. | A two-step stochastic hybrid estimation for GNSS carrier phase tracking in urban environments | |
CN104459747A (zh) | 一种gnss/ins组合导航方法及系统 | |
Vilà-Valls et al. | Multi-frequency GNSS robust carrier tracking for ionospheric scintillation mitigation | |
Amani et al. | Correlator‐based multipath detection technique for a global positioning system/GNSS receiver | |
CN110780320B (zh) | 一种软硬件一体化的卫星导航信号处理方法 | |
CN111458730B (zh) | 基于多普勒残差估计的gnss载波跟踪方法及接收机 | |
CN102565822A (zh) | 对gps l5信号的捕获方法及装置 | |
JP5519223B2 (ja) | 衛星信号受信装置 | |
Tang et al. | The explicit tuning investigation and validation of a full Kalman filter-based tracking loop in GNSS receivers | |
CN103983989B (zh) | 一种用于卫星定位导航接收系统的数字非整数跟踪环路 | |
Wu et al. | A BDS-3 B1C/B2a dual-frequency joint tracking architecture based on adaptive Kalman filter and extended integration time | |
Yao et al. | Spreading code phase measurement technique for snapshot GNSS receiver | |
CN107843905B (zh) | 一种快速高动态gnss频率搜索方法 | |
Jianfeng et al. | Low C/N0 carrier tracking loop based on optimal estimation algorithm in GPS software receivers |
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 |