CN104765050A - 一种新的北斗信号二次捕获算法 - Google Patents
一种新的北斗信号二次捕获算法 Download PDFInfo
- Publication number
- CN104765050A CN104765050A CN201510191930.4A CN201510191930A CN104765050A CN 104765050 A CN104765050 A CN 104765050A CN 201510191930 A CN201510191930 A CN 201510191930A CN 104765050 A CN104765050 A CN 104765050A
- Authority
- CN
- China
- Prior art keywords
- signal
- local
- secondary coding
- algorithm
- acquisition
- 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.)
- Pending
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
-
- 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/30—Acquisition or tracking or demodulation of signals transmitted by the system code 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)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及北斗卫星信号捕获方法,具体为一种新的北斗信号二次捕获算法,本发明首先在主码层面进行串行捕获;而后在整个二次编码周期上进行并行二次捕获的方法,以去除二次编码对积分时间的限制,达到提高接收机灵敏度的目的。本发明提出一种基于串并结合的北斗分层编码信号的捕获算法。本充分利用二次编码的特性搜索PRN码和二次编码的码相位和多普勒频偏。积分时间不再局限在一个主码周期内而是扩展到几个主码周期,在算法复杂度合理上升的同时提高了捕获灵敏度。而且为了在信号足够强的情况下减少捕获时间,算法采用大频偏搜索步长先捕获主码相位而后进行二次编码的精细搜索,因此能够实现快速捕获强信号。
Description
技术领域
本发明涉及北斗卫星信号捕获方法,具体为一种新的北斗信号二次捕获算法。
背景技术
北斗卫星导航系统(BeiDou Navigation Satellite System,BDS)中需要长时间的积分来获得较高的灵敏度,而BDS二代中新型扩频码的采用给接收机的设计带来了新的挑战,新型扩频码采用二次编码与主码的分层结构,这种结构大大限制了导频通道相干积分时间的增加。
BDS接收机的信号处理包括三个阶段:捕获、跟踪和位置,速度,时间(Position Velocity Time,PVT)解算。接收机通过捕获可以粗略地估算信号多普勒频偏和扩频码相位,跟踪阶段更加精准地同步多普勒频偏与码相位,进而提取导航电文和伪距信息,最终解算出PVT信息。随着BDS新技术的不断应用,未来接收机将面临复杂环境(如深山、峡谷、城市、室内)的严峻挑战。因此,接收机的捕获部分必须精心设计,确保在微弱信号环境下正常工作。事实上,捕获的目的是搜索所有可见卫星,并确定不同卫星扩频码(PRN码)相位。
接收信号与相应的本地载波和扩频码(PRN码)副本做相关运算,如果相关值超过提前设定的阈值,则捕获卫星成功,同时得到多普勒频偏和PRN码相位。通过在不同的载频和码相位上进行重复搜索,直到搜索到所有可见卫星。一般相关时间为单PRN码周期,如果需要提高接收机的灵敏度,单周期的相关时间是不够的。因此,可以叠加不同码周期的相关结果来提高灵敏度,分为相干积分(即平方之前求和)和非相干累加(平方之后求和)。在低信噪比的环境中,相干积分比非相干累加效果要好。关于这两种方法已经有文章见诸于报,本文不再赘述。对于传统GPS的L1信号,可以通过简单的方法克服导航电文数据位翻转,实现相干积分。然而,BDS二代中扩频码特殊的编码结构需要在捕获时考虑所有可能的二次编码NH码相位组合,致使相干积分的实现难度变大。
随着BDS二代的广泛使用,如何增加相干积分时间已成为研究热点。在[2]中作者讨论了一种树形结构的长时间相干积分方法,主要原理是:在指定长度的二次编码码片上尝试所有可能的相位组合来增加相干积分时间,这种方法虽然一定程度上提高了灵敏度,但相位组合的数量是按照二次编码码片数呈指数增长,因而只适用于较短的二次编码捕获,而且没有获取码相位。因此本文提出了一种新的算法,在一个PRN码周期上并行搜索PRN码相位,在整个NH码周期上进行并行二次搜索NH码相位,提高灵敏度的同时也捕获了NH码相位。
发明内容
本发明为了解决BDS二代中扩频码特殊的编码结构大大限制了导频通道相干积分时间的增加的问题,提供了一种新的北斗信号二次捕获算法。
本发明是采用如下的技术方案实现的:一种新的北斗信号二次捕获算法,包括以下步骤:
步骤一:BDS接收机接收到的RF信号经过射频前端放大、下变频和数字化后得到中频信号;
步骤二:对中频信号每隔1ms选取信号,一共选取m组作为输入信号,m可任意选定;
步骤三:设定多普勒频偏,生成1ms扩频码作为本地扩频码;
步骤四:将m组输入信号分别与经过本地副载波调制的本地扩频码进行一个周期的循环相乘,而后将乘积进行NSP点的累加得到相关输出,Nsp是一周期本地扩频码的样点数;
步骤五:相关输出经过Ns=20次的串行捕获,得到串行捕获结果;
步骤六:选取20ms本地二次编码补20ms零构成本地备用二次编码,本地备用二次编码经过FFT变换后取复共轭得到本地二次编码复共轭信号;
步骤七:将串行捕获结果经FFT变换后与本地二次编码复共轭信号相乘,相乘结果经IFFT变换到时域得到时域信号;
步骤八:将时域信号与设定的门限阈值进行比较,检测是否捕获到卫星信号,若时域信号没有超过阈值则将在下一个多普勒频偏上重复步骤三到步骤八,直到捕获到卫星信号或者搜索完所有的可能多普勒频偏。
本发明提出一种基于串并结合的北斗分层编码信号的捕获算法。本充分利用二次编码的特性搜索PRN码和二次编码的码相位和多普勒频偏。积分时间不再局限在一个主码周期内而是扩展到几个主码周期,在算法复杂度合理上升的同时提高了捕获灵敏度。而且为了在信号足够强的情况下减少捕获时间,算法采用大频偏搜索步长先捕获主码相位而后进行二次编码的精细搜索,因此能够实现快速捕获强信号。
附图说明
图1是本发明的串并结合北斗卫星信号捕获流程图。
图2是本发明北斗信号捕获装置的结构图。
具体实施方式
一种新的北斗信号二次捕获算法,包括以下步骤:
步骤一:BDS接收机接收到的RF信号经过射频前端放大、下变频和数字化后得到中频信号IF,X(nTs);
步骤二:对中频信号IF,X(nTs)每隔1ms选取信号,一共选取m组作为输入信号;
步骤三:设定多普勒频偏,生成1ms扩频码作为本地扩频码;
步骤四:将m组输入信号分别与经过本地副载波调制的本地扩频码进行一个周期的循环相乘,而后将乘积进行NSP点的累加得到相关输出yp(k),Nsp是一周期本地扩频码的样点数;
步骤五:相关输出yp(k)经过Ns=20次的串行捕获,得到串行捕获结果ysc,Ns是本地二次编码的码片数;
步骤六:选取20ms本地二次编码补20ms零构成本地备用二次编码,本地备用二次编码经过FFT变换后取复共轭得到本地二次编码复共轭信号;
步骤七:将串行捕获结果ysc经FFT变换后与本地二次编码复共轭信号相乘,相乘结果经IFFT变换到时域得到时域信号
步骤八:将时域信号与设定的门限阈值进行比较,检测是否捕获到卫星信号,若时域信号没有超过阈值则将在下一个多普勒频偏上重复步骤三到步骤八,直到捕获到卫星信号或者搜索完所有的可能多普勒频偏。
具体实施时,接收机接收到的卫星信号经过射频前端放大、下变频、数字化后得到中频信号IF,X(nTs),中频信号表达式如下:
其中Ts是采样周期;fIF是中心频率;τ是PRN码相位;fd是多普勒频偏;ηIF是高斯白噪声;X是任意北斗卫星信号(i.e.B1,B2);SX(t)是低通信号:sX(t)=sX,I(t)+jsX,Q(t)(2),现只考虑导频通道,公式(1)可以改写为:
串行搜索就是将中频信号IF,X(nTs)与经过本地副载波调制的本地扩频码进行一个PRN码周期的循环相乘,而后将乘积进行Nsp点的累加,其中Nsp是一周期PRN码的样点数,同相和正交分量共同组成FFT部分的输入,上述过程的累加结果yp(k)可以表示为:
yp(k)=yI(k)+jyQ(k),k∈1,2,...,Nsc (4),其中
其中,是第X路卫星信号的自相关函数;和fd是中频信号相对于本地扩频码的相位延迟与多普勒频偏;Tint是积分时间;eφ是载波频率误差;sc(t)是二次编码序列,τ0是二次编码码相位;ηI和ηQ分别为能量为的高斯噪声;k是一个主码周期上的样点数,k∈{1,2,...,Nsc}。
需要指出,这些相关计算的结果都与本地副载波有关。一旦本地副载波与中频信号同步时,所有的功率将被集中在同相部分。然而,由于本地副载波和中频信号之间的相位漂移,能量将分配在同相和正交分量。Ns是二次编码的码片数,将存储Ns次的串行捕获结果ysc作为FFT部分的输入,串行捕获结果ysc表达式如下:
ysc=[yp(1),yp(2),...,yp(Ns)] (7),为了捕获二次编码采用FFT单元进行搜索,FFT单元的输出ysc随后与本地二次编码FFT变换的复共轭相乘,乘积经过IFFT变换到时域,时域信号表达式如下:
Claims (1)
1.一种新的北斗信号二次捕获算法,其特征在于包括以下步骤:
步骤一:BDS接收机接收到的卫星信号经过射频前端放大、下变频和数字化后得到中频信号;
步骤二:对中频信号每隔1ms选取信号,一共选取m组作为输入信号,m可任意选定;
步骤三:设定多普勒频偏,生成1ms扩频码作为本地扩频码;
步骤四:将m组输入信号分别与经过本地副载波调制的本地扩频码进行一个周期的循环相乘,而后将乘积进行NSP点的累加得到相关输出,Nsp是一周期本地扩频码的样点数;
步骤五:相关输出经过Ns=20次的串行捕获,得到串行捕获结果;
步骤六:选取20ms本地二次编码补20ms零构成本地备用二次编码,本地备用二次编码经过FFT变换后取复共轭得到本地二次编码复共轭信号;
步骤七:将串行捕获结果经FFT变换后与本地二次编码复共轭信号相乘,相乘结果经IFFT变换到时域得到时域信号;
步骤八:将时域信号与设定的门限阈值进行比较,检测是否捕获到卫星信号,若时域信号没有超过阈值则将在下一个多普勒频偏上重复步骤三到步骤八,直到捕获到卫星信号或者搜索完所有的可能多普勒频偏。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510191930.4A CN104765050A (zh) | 2015-04-21 | 2015-04-21 | 一种新的北斗信号二次捕获算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510191930.4A CN104765050A (zh) | 2015-04-21 | 2015-04-21 | 一种新的北斗信号二次捕获算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104765050A true CN104765050A (zh) | 2015-07-08 |
Family
ID=53646998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510191930.4A Pending CN104765050A (zh) | 2015-04-21 | 2015-04-21 | 一种新的北斗信号二次捕获算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104765050A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105717522A (zh) * | 2016-02-23 | 2016-06-29 | 电子科技大学 | “北斗二代”b1频段弱信号捕获方法 |
CN106338749A (zh) * | 2015-12-30 | 2017-01-18 | 深圳艾科创新微电子有限公司 | 北斗卫星导航接收机及其nh码剥离、自相关方法及装置 |
CN106772356A (zh) * | 2017-01-05 | 2017-05-31 | 西安电子科技大学 | 单通道单脉冲系统的扩频角跟踪信号捕获方法 |
CN108011652A (zh) * | 2016-10-28 | 2018-05-08 | 上海复控华龙微系统技术有限公司 | 一种扩频信号捕获的方法和装置 |
CN109143285A (zh) * | 2017-06-27 | 2019-01-04 | 航天恒星科技有限公司 | 应用于姿态多变高动态目标的定位通报系统 |
CN109239743A (zh) * | 2018-09-17 | 2019-01-18 | 西安开阳微电子有限公司 | 一种卫星信号捕获方法及装置 |
CN112578411A (zh) * | 2020-11-06 | 2021-03-30 | 中国科学院国家空间科学中心 | 一种微弱bds-3 b1c基带信号的捕获方法及系统 |
CN112803968A (zh) * | 2020-12-30 | 2021-05-14 | 南京天际易达通信技术有限公司 | 一种无人机机载测控方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2110682A2 (en) * | 2008-04-18 | 2009-10-21 | Honeywell International Inc. | Systems and methods for fast GNSS signals acquisition |
CN103645483A (zh) * | 2013-12-09 | 2014-03-19 | 西安电子科技大学昆山创新研究院 | 一种弱信号环境下北斗信号捕获方法 |
CN104133225A (zh) * | 2014-07-29 | 2014-11-05 | 豪芯微电子科技(上海)有限公司 | 北斗信号处理方法和装置 |
-
2015
- 2015-04-21 CN CN201510191930.4A patent/CN104765050A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2110682A2 (en) * | 2008-04-18 | 2009-10-21 | Honeywell International Inc. | Systems and methods for fast GNSS signals acquisition |
CN103645483A (zh) * | 2013-12-09 | 2014-03-19 | 西安电子科技大学昆山创新研究院 | 一种弱信号环境下北斗信号捕获方法 |
CN104133225A (zh) * | 2014-07-29 | 2014-11-05 | 豪芯微电子科技(上海)有限公司 | 北斗信号处理方法和装置 |
Non-Patent Citations (1)
Title |
---|
DENG-AO LI ET AL.: "Analysis and Application of Secondary Code in BDS", 《2014 ASIA-PACIFIC SERVICES COMPUTING CONFERENCE》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106338749A (zh) * | 2015-12-30 | 2017-01-18 | 深圳艾科创新微电子有限公司 | 北斗卫星导航接收机及其nh码剥离、自相关方法及装置 |
CN105717522B (zh) * | 2016-02-23 | 2019-01-01 | 电子科技大学 | “北斗二代”b1频段弱信号捕获方法 |
CN105717522A (zh) * | 2016-02-23 | 2016-06-29 | 电子科技大学 | “北斗二代”b1频段弱信号捕获方法 |
CN108011652B (zh) * | 2016-10-28 | 2020-05-12 | 上海复控华龙微系统技术有限公司 | 一种扩频信号捕获的方法和装置 |
CN108011652A (zh) * | 2016-10-28 | 2018-05-08 | 上海复控华龙微系统技术有限公司 | 一种扩频信号捕获的方法和装置 |
CN106772356A (zh) * | 2017-01-05 | 2017-05-31 | 西安电子科技大学 | 单通道单脉冲系统的扩频角跟踪信号捕获方法 |
CN109143285A (zh) * | 2017-06-27 | 2019-01-04 | 航天恒星科技有限公司 | 应用于姿态多变高动态目标的定位通报系统 |
CN109239743A (zh) * | 2018-09-17 | 2019-01-18 | 西安开阳微电子有限公司 | 一种卫星信号捕获方法及装置 |
CN109239743B (zh) * | 2018-09-17 | 2023-04-21 | 西安开阳微电子有限公司 | 一种卫星信号捕获方法及装置 |
CN112578411A (zh) * | 2020-11-06 | 2021-03-30 | 中国科学院国家空间科学中心 | 一种微弱bds-3 b1c基带信号的捕获方法及系统 |
CN112578411B (zh) * | 2020-11-06 | 2023-10-13 | 中国科学院国家空间科学中心 | 一种微弱bds-3 b1c基带信号的捕获方法及系统 |
CN112803968A (zh) * | 2020-12-30 | 2021-05-14 | 南京天际易达通信技术有限公司 | 一种无人机机载测控方法 |
CN112803968B (zh) * | 2020-12-30 | 2021-07-30 | 南京天际易达通信技术有限公司 | 一种无人机机载测控方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104765050A (zh) | 一种新的北斗信号二次捕获算法 | |
CN101132191B (zh) | 一种gnss接收机基带信号处理方法 | |
CN104280750B (zh) | 基于部分匹配滤波fft算法的长码捕获系统及方法 | |
CN106772455B (zh) | 一种基于惯性信息辅助与参数估计的gnss反欺骗环路跟踪方法 | |
CN102928853B (zh) | 一种捕获北斗d1卫星导航系统弱信号的方法 | |
CN105277955A (zh) | 一种用于gps卫星导航的基带信号处理方法 | |
CN103926601B (zh) | 基于合成相关函数的boc(15,2.5)调制方式捕获方法 | |
CN108254767B (zh) | 一种boc信号捕获方法和基带同步接收机 | |
CN103499825B (zh) | 一种高轨航天器bds弱信号快速捕获方法 | |
CN104570017B (zh) | 一种针对北斗导航卫星系统D1导航电文中的Neumann‑Hoffman码规避与剥离方法 | |
CN105759288A (zh) | 基于特征序列的北斗b1i弱信号捕获方法 | |
CN104833988A (zh) | 大多普勒动态下高数据速率直扩信号的捕获系统及方法 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN105553506B (zh) | 一种长码扩频信号的快速捕获方法及装置 | |
CN103091687A (zh) | 北斗接收机多卫星信号联合捕获方法 | |
CN103543456A (zh) | 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法 | |
CN104181556A (zh) | 一种基于重叠差分循环相干积分的boc调制信号捕获方法 | |
CN104199057A (zh) | 一种基于压缩感知技术的boc信号无模糊捕获算法 | |
CN101520505A (zh) | 全球卫星导航系统弱信号下的自适应相干累积联合捕获方法 | |
CN103926604A (zh) | 基于重叠差分循环相干积分的弱信号捕获方法 | |
CN103941269A (zh) | 用于卫星导航系统的pn码捕获方法 | |
CN103439718A (zh) | 一种高阶boc调制信号的无模糊跟踪单元 | |
CN103197328A (zh) | 一种基于2ms数据存储的高动态卫星接收机快速捕获方法 | |
CN102680988B (zh) | 高灵敏度导航卫星信号非线性捕获方法及装置 | |
CN104035109A (zh) | 基于重叠1/5比特差分循环相干积分的弱信号捕获方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150708 |