[go: up one dir, main page]

CN110646818B - 一种高灵敏度的卫星导航精捕方法 - Google Patents

一种高灵敏度的卫星导航精捕方法 Download PDF

Info

Publication number
CN110646818B
CN110646818B CN201810671102.4A CN201810671102A CN110646818B CN 110646818 B CN110646818 B CN 110646818B CN 201810671102 A CN201810671102 A CN 201810671102A CN 110646818 B CN110646818 B CN 110646818B
Authority
CN
China
Prior art keywords
correlation
code
short
path
result
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
CN201810671102.4A
Other languages
English (en)
Other versions
CN110646818A (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.)
Beijing Automation Control Equipment Institute BACEI
Original Assignee
Beijing Automation Control Equipment Institute BACEI
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 Beijing Automation Control Equipment Institute BACEI filed Critical Beijing Automation Control Equipment Institute BACEI
Priority to CN201810671102.4A priority Critical patent/CN110646818B/zh
Publication of CN110646818A publication Critical patent/CN110646818A/zh
Application granted granted Critical
Publication of CN110646818B publication Critical patent/CN110646818B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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

本发明公开了一种高灵敏度的卫星导航精捕方法。普通的并行相关+FFT(离散傅里叶变换的快速算法)的捕获方法在卫星导航方面已经被广泛应用。本发明采用了一种延迟输出相关结果的并行相关+FFT的方法，该方法可以很大的节省多路相关完成后需要锁存相关值的寄存器，从而减少了占用FPGA的逻辑资源；同时还使用了一种特殊方法消除部分卫星信号的数据比特(还有如BD2的NH码等)对相关长度的影响，从而提高数据和伪码相关的长度，最终达到提高精捕灵敏度的目的。

Description

一种高灵敏度的卫星导航精捕方法
技术领域
本发明属于一种卫星导航接收机信号处理方法，具体涉及一种高灵敏度的卫星导航精捕方法。
背景技术
卫星导航系统(北斗、GPS、GLONASS等)是一种高精度的导航定位系统。卫星接收机通过接收天线接收卫星信号对卫星信号进行下变频处理，然后对下变频信号进行基带数字信号处理，经过导航解算处理后，可以给出精确的经度、纬度、高度、速度及时间等信息。目前卫星导航广泛的运用在军事及民用领域。
基带数字信号处理主要包括对卫星信号的捕获、跟踪等功能。由于算法限制，捕获灵敏度一般要远低于跟踪灵敏度，使得提高对卫星导航信号的捕获能力变得至关重要。在捕获卫星信号时，特别是长码捕获时，具备精捕功能更是必不可少的。一个好的精捕方法占用FPGA资源量极少，但能很大的减小虚警概率，提高牵引成功率。
目前广泛使用的精捕方法均为通过采用1ms(最小单元的相关运算)相干积分及多次非相干积分来提高精捕灵敏度，这种方法带来的问题如下：
(1)需要用长时间的非相干积分才能达到比较高的灵敏度；
(2)长时间的非相干积分无法满足高动态场景的卫星信号捕获；
(3)长时间的非相干积分捕获时间较长。
发明内容
本发明的目的是提供一种高灵敏度的卫星导航精捕方法，通过提高相干积分时间，在提高精捕灵敏度的同时，不增加捕获时间，且满足高动态需求。
本发明是这样实现的，一种高灵敏度的卫星导航精捕方法，包括以下步骤：
步骤1：产生通道本地码、设置载波频率；
步骤2：产生I、Q两路混频结果；
步骤3：I、Q两路信号与本地码进行多路短时相关运算；
步骤4：储存多路相关运算结果；
步骤5：短时相关结果与数据跳变系数进行相关运算；
步骤6：相关运算结果进行FFT；
步骤7：FFT结果进行求模运算并输出捕获结果。
所述的步骤1包括以下操作：
根据粗捕得到的码相位置入跟踪通道，通过停止通道码发生器来延时通道码相位(长码为置时间到相关芯片等待芯片产生码)，将本地码调整到超前粗捕结果L个码片；
根据粗捕得到多普勒频率f0，并将f0置入跟踪通道，将对应跟踪通道的载波频率调整为f(本地载波初始频率)+f0，同时产生相应频率的I、Q两路正交载波信号。
所述的步骤2包括以下操作：
将I、Q两路正交载波信号与中频信号进行相关运算，产生I、Q两路混频信号。
所述的步骤3包括以下操作：
对跟踪通道输出的本地码进行N路码片延迟，产生N路并行码输出，将N路并行码与输入的I、Q两路信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后将N路短时累加结果延迟输出。
所述的步骤4包括以下操作：
每次acc_complete(该信号为第一路短时相关完成标志)为高电平后输出N组短时相关结果到RAM中，N个时钟后完成该次写操作，一次捕获进行M次写操作，RAM深度为M*N。
所述的步骤5包括以下操作：
将RAM中缓存的短时相关结果与事先预存在ROM中的数据跳变系数进行对应相乘，产生新的短时相关结果。
所述的步骤6包括以下操作：
将步骤5相关运算完的数据进行FFT运算，每次捕获共进行P*N次Q点FFT运算(其中Q大于等于M)。
所述的步骤7包括以下操作：
对FFT结果进行求模运算，并比较出最大值并保持最大值对应的捕获码相位值及多普勒频率值。
所述的步骤3中的对跟踪通道输出的本地码进行N路码片延迟，产生N路并行码输出，将N路并行码与输入的I、Q两路信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后通过N路短时相关完成标记控制N路短时累加结果延迟输出。
所述的步骤5中的将RAM中缓存的短时相关结果与事先预存在ROM中的数据跳变系数进行对应相乘，产生新的短时相关结果。
本发明的优点是，(1)采用一种延迟输出相关结果的并行相关+FFT的方法(见步骤3)，很大的节省多路相关完成后需要锁存相关值的寄存器，从而减少了占用FPGA的逻辑资源，从而减小了系统功耗；(2)使用了一种从ROM读取数据跳变系数与相关结果相乘的方法(见步骤5)，消除了卫星信号数据跳变或NH码跳变等对相关长度的限制影响，从而提高数据和伪码相。
具体实施方式
下面结合具体实施例对本发明进行详细介绍：
本发明的原理是：对跟踪通道输出的本地码进行N路码片延迟，产生N路并行码输出，将N路并行码与输入的基带信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后将N路短时累加结果延迟输出，并与ROM中读出的数据跳变系数相乘，进行FFT变换并将变换结果取模比较最大值，最后将最大值及最大值对应的码相位和多普勒频率锁存。本发明的方法具体包括以下步骤：
步骤1：产生通道本地码、设置载波频率
根据粗捕得到的码相位调整跟踪通道，调整跟踪通道是通过暂停通道码发生器的运行来延时通道码相位(长码为置时间到相关芯片等待芯片产生码)，最终将本地码调整到超前粗捕结果L个码片。
通过粗捕得到了载波多普勒频率f0，将f0换算成本地压控振荡器增量置入跟踪通道的压控震荡器来调整跟踪通道的载波频率，经过通道调整后的载波频率调整为f(本地载波初始频率)+f0，同时产生相应频率的I、Q两路正交本地载波信号。
步骤2：产生I、Q两路混频结果：
将步骤1输出的I、Q两路正交本地载波信号分别与射频前端经过数字信号处理的中频信号进行混频，产生I、Q两路混频结果。
步骤3：I、Q两路信号与本地码进行多路短时相关运算：
对步骤1跟踪通道输出的本地码进行N路码片延迟(一般N＝2*L*TCODE/TCLK，其中TCODE为一个码片的时长、TCLK为工作时钟的周期)，产生N路并行码输出，然后将N路并行码与输入的I、Q两路信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后通过N路短时相关完成标记控制N路短时累加结果延迟输出。
具体如下：
(1)当即时路(第1路未延迟)的相关结果每次累加到T0(ms)时长的最后1个工作时钟周期时，产生1个只有1个工作时钟周期高电平的累加完成信号acc_complete，当acc_complete为1时输出该即时路的累加结果，同时将累加值置为当前即时的相关值，开始进行下一轮累加运算；
(2)在acc_complete为高电平的后1个工作时钟周期，输出延迟1路的累加结果，同时将累加值置为当前延迟1路的相关值，开始进行下一轮累加运算；
(3)在acc_complete为高电平的后2拍，输出延迟2路的累加结果，同时将累加值置为当前延迟2路的相关值，开始进行下一轮累加；
(4)依此直到在acc_complete为高电平的后N-1拍，输出延迟N-1路的累加结果，同时将累加值置为当前延迟N-1路的相关值，开始进行下一轮累加；即完成了1次N路短时相关运算；
(5)在下1个acc_complete到来时进行上述同样操作，直到完成T1(ms)的相关运算(该操作进行M次，M×T0＝T1,T0必须能整除1)；
(6)通过(1)～(4)的方法，使N路的相关结果不是在同一个时刻计算完成，而是在acc_complete为高电平后的N个工作时钟依次每个时钟计算完一路的短时相关结果并将短时相关值到缓存到RAM中(该短时相关结果只保存1个时钟后即开始参与下一次相关累加运算中，因此该值一旦输出即将其输入到RAM中)；acc_complete为高电平N个工作时钟后完成一次N路相关结果的输出；
(7)在下1个acc_complete到来时进行上述同样操作，直到完成T1(ms)的相关运算(该操作进行M次，M×T0＝T1,T0必须能整除1)。
步骤4：储存多路相关运算结果
每次acc_complete为高电平的N个工作时钟内依次将步骤3输出的N组短时相关结果输入到RAM中，N个时钟后完成该次写操作；一次捕获进行M次写操作，RAM深度为M*N。
步骤5：短时相关结果与数据跳变系数进行相关运算
将RAM中缓存的短时相关结果与事先预存在ROM中的数据跳变系数进行对应相乘，具体如下：
(1)从RAM中读取短时相关结果，同时从ROM中循环读取预存的第一组数据跳变系数，RAM和ROM的读使能信号一致，每读M个数后停止读Q-M个工作时钟。RAM的读地址初始值为0，在读使能为高时一直加1，每次加到M*N-1即停止累加，地址置0；ROM的读地址初始值为0，读使能为高时每1/T0个工作时钟地址一直加1，每次加到T1-1时地址置0，再重新开始加，按此操作循环N次。依此进行后，RAM中N路短时相关结果根据从ROM中循环读取的第1组对应的数据跳变系数进行正负判断，当为1时对相关结果取负输出，当为0时对相关结果取正输出；
(2)依次如(1)操作直到第P次，从RAM中读取短时相关结果，同时从ROM中循环读取预存的第P组数据跳变系数，RAM和ROM的读使能信号一致，每读M个数后停止读Q-M个工作时钟。RAM的读地址初始值为0，在读使能为高时一直加1每次加到M*N-1即停止累加，地址置0；ROM的读地址初始值为(P-1)*T1，读使能为高时每1/T0个工作时钟地址一直加1，每次加到P*T1-1时地址置0，再重新开始加，按此操作循环N次。依此进行后，RAM中N路短时相关结果根据从ROM中循环读取的第P组对应的数据跳变系数进行正负判断，当为1时对相关结果取负输出，当为0时对相关结果取正输出；
(3)ROM缓存深度为T1*P，共产生P*N组相关结果，每组M个相关结果，其中以20ms长度的GPS信号为例ROM数据跳变系数设置方法如下(其它卫星信号依次拓展即可)：将GPS数据跳变的可能均考虑到，因为GPS数据bit长为20ms，因此共20种可能，所以这里的T1等于20(1ms个数)、P等于20(数据跳变可能性)，具体见表1。
表1 20ms长度的GPS信号
Figure GDA0003801882870000061
Figure GDA0003801882870000071
步骤6：相关运算结果进行FFT
对步骤5运算完的相关结果进行FFT运算，每次捕获共进行P*N次Q点FFT运算(其中Q大于等于M)。
步骤7：FFT结果进行求模运算并输出捕获结果
对FFT结果进行求模运算，并比较出P*N*Q个FFT输出结果的最大值，并保持最大值及最大值对应的码相位值及多普勒频率值。最后比较最大值与门限的大小，若最大值大于门限值，则进牵引，若最大值小于门限，则判断捕获失败。
在FPGA上实现该算法，其中M＝40，N＝20，T1＝10，T0＝0.25，Q＝64。设置输入信号与码片相位差为1个时钟周期，输入信号与本地载波频率差为2000Hz。通过modsim10.1a仿真可得最后的捕获峰值为13894，捕获地址对应的用2进制表示为11’b00001100000，其中低6位代表多普勒频率对应频率为2000Hz，高5位对应码片偏差为1。仿真结果与预置testbench的输入信号一致。

Claims (2)

1.一种高灵敏度的卫星导航精捕方法，其特征在于：包括以下步骤：
步骤1：产生通道本地码、设置载波频率；
步骤2：产生I、Q两路混频结果；
步骤3：I、Q两路信号与本地码进行多路短时相关运算；
步骤4：储存多路相关运算结果；
步骤5：短时相关结果与数据跳变系数进行相关运算；
步骤6：相关运算结果进行FFT；
步骤7：FFT结果进行求模运算并输出捕获结果；
所述的步骤1包括以下操作：
根据粗捕得到的码相位置入跟踪通道，通过停止通道码发生器来延时通道码相位，其中，长码为置时间到相关芯片等待芯片产生码，将本地码调整到超前粗捕结果L个码片；
根据粗捕得到多普勒频率f0，并将f0置入跟踪通道，将对应跟踪通道的载波频率调整为f+f0，其中，f为本地载波初始频率，同时产生相应频率的I、Q两路正交载波信号；
所述的步骤2包括以下操作：
将I、Q两路正交载波信号与中频信号进行相关运算，产生I、Q两路混频信号；
所述的步骤3包括以下操作：
对跟踪通道输出的本地码进行N路码片延迟，产生N路并行码输出，将N路并行码与输入的I、Q两路信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后将N路短时累加结果延迟输出；
所述的步骤4包括以下操作：
每次acc_complete为高电平后输出N组短时相关结果到RAM中，其中acc_complete信号为第一路短时相关完成标志，N个时钟后完成该次写操作，一次捕获进行M次写操作，RAM深度为M*N；
所述的步骤5包括以下操作：
将RAM中缓存的短时相关结果与事先预存在ROM中的数据跳变系数进行对应相乘，产生新的短时相关结果；
所述的步骤6包括以下操作：
将步骤5相关运算完的数据进行FFT运算，每次捕获共进行P*N次Q点FFT运算，其中Q大于等于M；
所述的步骤7包括以下操作：
对FFT结果进行求模运算，并比较出最大值并保持最大值对应的捕获码相位值及多普勒频率值。
2.如权利要求1所述的一种高灵敏度的卫星导航精捕方法，其特征在于：所述的步骤3中的对跟踪通道输出的本地码进行N路码片延迟，产生N路并行码输出，将N路并行码与输入的I、Q两路信号进行相关运算，并将N路并行相关结果分别进行N路短时累加，最后通过N路短时相关完成标记控制N路短时累加结果延迟输出。
CN201810671102.4A 2018-06-26 2018-06-26 一种高灵敏度的卫星导航精捕方法 Active CN110646818B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810671102.4A CN110646818B (zh) 2018-06-26 2018-06-26 一种高灵敏度的卫星导航精捕方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810671102.4A CN110646818B (zh) 2018-06-26 2018-06-26 一种高灵敏度的卫星导航精捕方法

Publications (2)

Publication Number Publication Date
CN110646818A CN110646818A (zh) 2020-01-03
CN110646818B true CN110646818B (zh) 2023-01-13

Family

ID=69008781

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810671102.4A Active CN110646818B (zh) 2018-06-26 2018-06-26 一种高灵敏度的卫星导航精捕方法

Country Status (1)

Country Link
CN (1) CN110646818B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111796306B (zh) * 2020-07-31 2023-05-02 北京中捷时代航空科技有限公司 一种导航卫星信号接收方法及接收机
CN114741133B (zh) * 2022-04-21 2023-10-27 中国航空无线电电子研究所 一种基于模型的综合模块化航电系统资源配置与评估方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101718874A (zh) * 2009-12-14 2010-06-02 浙江大学 高灵敏度导航卫星接收机中的信号快速捕获方法及装置
KR20110068128A (ko) * 2009-12-15 2011-06-22 한국전자통신연구원 신호 획득 장치와 방법 및 신호 획득 장치의 파라미터 설정 방법
CN102608626A (zh) * 2012-03-09 2012-07-25 暨南大学 一种高灵敏度卫星导航信号捕获方法及系统
CN103698780A (zh) * 2013-09-05 2014-04-02 北京捷星广达科技有限责任公司 一种针对二次编码调制的卫星导航信号的捕获引擎
CN106707309A (zh) * 2016-11-23 2017-05-24 北京自动化控制设备研究所 一种载波伪码二维搜索的精捕方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101718874A (zh) * 2009-12-14 2010-06-02 浙江大学 高灵敏度导航卫星接收机中的信号快速捕获方法及装置
KR20110068128A (ko) * 2009-12-15 2011-06-22 한국전자통신연구원 신호 획득 장치와 방법 및 신호 획득 장치의 파라미터 설정 방법
CN102608626A (zh) * 2012-03-09 2012-07-25 暨南大学 一种高灵敏度卫星导航信号捕获方法及系统
CN103698780A (zh) * 2013-09-05 2014-04-02 北京捷星广达科技有限责任公司 一种针对二次编码调制的卫星导航信号的捕获引擎
CN106707309A (zh) * 2016-11-23 2017-05-24 北京自动化控制设备研究所 一种载波伪码二维搜索的精捕方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
全球卫星导航系统信号捕获方法研究综述;李灯熬 等;《计算机工程与设计》;20160131;第5-11页 *

Also Published As

Publication number Publication date
CN110646818A (zh) 2020-01-03

Similar Documents

Publication Publication Date Title
CN103033828B (zh) 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法
CN104422939B (zh) 用于关联从导航卫星系统接收到的信号的方法和装置
CN104833988B (zh) 大多普勒动态下高数据速率直扩信号的捕获系统及方法
CN101118281B (zh) 一种用于gps接收机对本地晶体振荡器频率进行自校正的方法
CN107153206A (zh) 一种基于fft的北斗卫星信号快速捕获方法
CN102353968A (zh) 一种基于fpga的gps信号捕获方法及其系统
CN100338475C (zh) 具有用于改进信号的接收动态的装置的射频信号接收机
CN110646818B (zh) 一种高灵敏度的卫星导航精捕方法
CN100581068C (zh) 一种多星座卫星信号的混合串并行相关捕获方法
CN106707309A (zh) 一种载波伪码二维搜索的精捕方法
CN113031019A (zh) 一种基于相关峰动态拖曳的卫星导航再生式欺骗产生方法
CN102854516A (zh) 一种gnss接收机中的载噪比估计方法及其系统
CN101881818B (zh) 对微弱信号进行快速检测的装置及方法
CN113466898B (zh) 用于gnss的相关器电路、gnss接收机及其接收信号的方法
CN106526633A (zh) 一种gnss基带信号的捕获方法及装置
CN104536020A (zh) 一种gnss信号捕获方法
CN104931980A (zh) 载波相位测量半周模糊解除方法
CN102508265B (zh) 基于信号分离估计理论的卫星导航信号多径干扰抑制方法
CN105527636A (zh) 导航信号捕获转跟踪方法及系统
CN103954978B (zh) 基于pxi架构动态导航信号源的实现方法
CN104614742B (zh) 北斗星基高精度差分信息接收实现方法
Zeng et al. Realization of baseband signal processing for Beidou/GPS multi-mode receiver by FPGA
CN1954233B (zh) 一种全球卫星定位装置
CN110149197A (zh) 一种用于时钟同步系统的高精度同步方法及同步系统
CN108152843B (zh) 一种附加拟稳基准的钟差数据融合计算方法

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