CN109581434B - 一种B2a信号捕获方法及装置 - Google Patents

一种B2a信号捕获方法及装置 Download PDF

Info

Publication number
CN109581434B
CN109581434B CN201811469341.8A CN201811469341A CN109581434B CN 109581434 B CN109581434 B CN 109581434B CN 201811469341 A CN201811469341 A CN 201811469341A CN 109581434 B CN109581434 B CN 109581434B
Authority
CN
China
Prior art keywords
signal
code phase
current satellite
satellite
signal transmitted
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
CN201811469341.8A
Other languages
English (en)
Other versions
CN109581434A (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.)
Xi'an Kaiyang Microelectronic Co ltd
Original Assignee
Xi'an Kaiyang Microelectronic Co 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 Xi'an Kaiyang Microelectronic Co ltd filed Critical Xi'an Kaiyang Microelectronic Co ltd
Priority to CN201811469341.8A priority Critical patent/CN109581434B/zh
Publication of CN109581434A publication Critical patent/CN109581434A/zh
Application granted granted Critical
Publication of CN109581434B publication Critical patent/CN109581434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related
    • 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/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/29Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明属于卫星导航技术领域领域,公开了一种B2a信号捕获方法及装置。该方法包括:捕获当前卫星发送的B1C信号的多普勒频率和码相位,根据B1C信号对应的多普勒频率和码相位确定B2a信号的码相位搜索范围,进而可以在较小的码相位搜索范围内对B2a码相位进行捕获搜索。本发明能够解决由于B2a信号码周期短、积分时长不够所导致的现有捕获方法捕获灵敏度不足的问题,能够快速捕获B2a信号,且简单易行,灵敏度较高。

Description

一种B2a信号捕获方法及装置
技术领域
本发明涉及卫星导航技术领域,尤其涉及一种B2a信号捕获方法及装置。
背景技术
北斗卫星导航系统是由中国自行研制的全球卫星导航系统,是全球四大卫星导航系统之一。北斗卫星导航系统简称北斗系统,分为三个建设发展阶段,第一阶段为北斗一代,第二阶段为北斗二代,第三阶段为北斗三代。
目前,北斗三代系统空间星座与用户终端之间公开服务信号为BIC信号和B2a信号。B1C信号和B2a信号在北斗三代中圆地球轨道卫星和倾斜地球同步轨道卫星上播发。其中,B1C信号的主码周期为10毫秒,使用传统的码相位并行捕获方法可以得到10毫秒的相干积分结果;B2a信号的主码周期仅为1毫秒,每毫秒调制1个子码,子码周期为5,所以每毫秒的子码调制符号可能不同。因此,如果使用现有的码相位并行捕获方法对B2a信号直接进行捕获,最多只能得到1毫秒的积分结果,这会使得B2a信号的捕获灵敏度非常低,甚至出现捕获不到B2a信号的情况。
综上所述,现有的码相位并行捕获方法不再适用于新的B2a信号,目前急需一种适用于B2a信号的捕获方法。
发明内容
有鉴于此,本发明提供一种B2a信号捕获方法及装置,以至少解决由于B2a信号码周期短、积分时长不够所导致的现有捕获方法捕获灵敏度不足的问题,能够快速捕获B2a信号,且简单易行,灵敏度较高。
为达到上述目的,本发明采用如下技术方案:
第一方面,提供一种B2a信号捕获方法,包括:
步骤S1,获取北斗天线接收到的北斗卫星信号,并对所述北斗卫星信号依次进行下变频处理和模拟/数字A/D转换,得到数字中频信号;其中,所述北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号;
步骤S2,根据所述数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位;
步骤S3,利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获。第二方面,提供一种B2a信号捕获装置,包括:数字中频信号获取单元、B1C信号捕获单元以及B2a信号捕获单元;
其中,所述数字中频信号获取单元,用于获取北斗天线接收到的北斗卫星信号,并对所述北斗卫星信号依次进行下变频处理和A/D转换,得到数字中频信号;其中,所述北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号;
所述B1C信号捕获单元,用于根据所述数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位;根据当前卫星发送的B1C信号的多普勒频率fdop_B1C,计算当前卫星发送的B2a信号的多普勒频率;根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围;
所述B2a信号捕获单元,用于利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获。
上述本发明提供的B2a信号捕获方法及装置,基于B1C信号和B2a信号之间的特性——二者频率存在一定的比例关系且比特边界沿基本对齐,通过利用B1C信号的捕获结果辅助实现B2a信号的捕获。具体来说,即是先捕获得到当前卫星发送的B1C信号对应的多普勒频率和码相位,进而根据B1C信号对应的多普勒频率和码相位确定B2a信号的码相位起始沿,进而可以在较小范围内对B2a码相位进行捕获搜索,因此本发明提供的B2a信号捕获方法及装置能够减少计算量,提高效率。此外,现有的码相位并行捕获方法最多只能得到1毫秒的积分结果,而基于本发明提供的B2a信号捕获方法及装置能够获得1毫秒的积分结果,因此相比现有技术,本发明提供的B2a信号捕获方法及装置能够提高捕获灵敏度。经验证,本发明提供的B2a信号捕获方法及装置将捕获灵敏度提高了7db左右。
综上,本发明提供的B2a信号捕获方法及装置,能够解决由于B2a信号码周期短、积分时长不够所导致的现有捕获方法捕获灵敏度不足的问题,快速捕获B2a信号,且简单易行,灵敏度较高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种B2a信号捕获方法的流程示意图;
图2为现有的并行码相位捕获方法的原理框图;
图3为本发明实施例提供的另一种B2a信号捕获方法的流程示意图;
图4为本发明实施例提供的一种B2a信号捕获装置的组成示意图;
图5为本发明实施例提供的另一种B2a信号捕获装置的组成示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1所示为本发明实施例提供的一种B2a信号捕获方法的流程示意图。
如图1所示,本发明实施例提供的B2a信号捕获方法,包括以下步骤:
步骤S1,获取北斗天线接收到的北斗卫星信号,并对北斗卫星信号依次进行下变频处理和模拟/数字A/D转换,得到数字中频信号。
其中,所述的北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号。本领域技术人员可以理解,捕获的目的即在于确定该北斗卫星信号中是否包含指定卫星发送的信号,并估算信号的载波多普勒频率和码相位。
步骤S2,根据数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位。
其中,当前卫星为当前指定的捕获卫星。
需要说明的是,这一步骤中,可利用现有的并行码相位捕获方法对当前卫星发送的B1C信号进行捕获,以获得对应的载波多普勒频率和码相位。具体的,请参考图2,将数字中频信号分别与I支路和Q支路上某一频率的复制正弦和复制余弦载波信号混频后,对复数形式的混频结果i+jq进行快速傅里叶变换,得到傅里叶变换结果,同时对伪码发生器产生的伪码进行快速傅里叶变换,并求复数共轭,将伪码傅里叶变换的复数共轭与混频结果的傅里叶变换相乘,接着将所得的乘积经傅里叶反变换并取模,得到在时域内的相关结果,最后将相关幅值与捕获门限进行比较,若相关幅值大于捕获门限,则捕获成功,并进入对多普勒频率估计的精细化过程,否则,载波数控振荡器(Numerically Controlled Oscillator,NCO)复制另一个频率值的正弦和余弦载波,然后类似地完成对下一个频带的搜索与检测。
步骤S3,利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获。
优选的,如图3所示,本发明实施例提供的B2a信号捕获方法中,步骤S3具体可以包括步骤S3a至步骤S3e:
步骤S3a,根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率,以及根据当前卫星发送的B1C信号的码相位,确定B2a信号的码相位搜索范围。具体的,步骤S3a中,根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率,可以包括:
根据当前卫星发送的B1C信号的多普勒频率,利用第一预设公式,计算当前卫星发送的B2a信号的多普勒频率。
其中,第一预设公式如公式1所示:
式中,fdop_B2a表示当前卫星发送的B2a信号的多普勒频率,fB1C表示B1C信号的载波频率,fB2a表示B2a信号的载波频率,fdop_B1C表示当前卫星发送的B1C信号的多普勒频率。
具体的,步骤S3a中,根据当前卫星发送的B1C信号的码相位,确定B2a信号的码相位搜索范围,具体可以包括:
根据当前卫星发送的B1C信号的码相位,利用第二预设公式估算当前卫星发送的B2a信号的比特边界值;
根据B2a信号的码相位边界值,确定B2a信号的码相位搜索范围。
其中,第二预设公式如公式2所示:
式中,n表示当前卫星发送的B2a信号的比特边界值,L表示主码长度,pB1C表示当前卫星发送的B1C信号的码相位,fs表示采样频率,fCA表示码速率。
需要说明的是,B2a信号的码相位搜索范围具体可以是以估算的B2a信号的码相位边界值为中心的一定范围内的码相位。优选的,可以估算的B2a信号的码相位边界值为中心的左右各32个半码片构成的码相位范围作为码相位搜索范围。
步骤S3b,根据B2a信号的多普勒频率和中频频率,生成本地载波。
具体的,可根据如下所示的公式3生成本地载波:
式中,fc表示本地载波,fi表示中频频率,中频频率的具体取值取决于设备的采样频率。
步骤S3c,将数字中频信号的起始位置依次移动至B2a信号的码相位搜索范围内的各半码片处,得到各半码片对应的卫星中频信号;再将本地载波分别与各半码片对应的卫星中频信号进行混频以剥离载波,进而对剥离载波后的信号进行重采样,得到各半码片对应的重采样数据。
步骤S3d,生成当前卫星对应的B2a信号的伪随机码,利用伪随机码对各半码片对应的重采样数据分别进行循环相关,得到各半码片对应的1ms相干积分结果;根据子码符号,分别对各半码片对应的1ms相干积分结果进行5ms的相干累加,得到各半码片对应的5ms相干累加数据。
其中,需要说明的是,本领域技术人员可以理解,之所以获得5ms相干累加数据,是因为B2a信号的子码周期为5。
步骤S3e,根据各半码片对应的5ms相干累加数据,检测是否成功捕获当前卫星发送的B2a信号,并在成功捕获当前卫星发送的B2a信号后,确定当前卫星发送的B2a信号的码相位。
具体的,步骤S3e具体包括:
对N组相干累加数据进行取模操作,得到对应的N个模值M1,M2,…,Mi,…,MN,确定N个模值M1,M2,…,Mi,…,MN中的最大模值Mmax
判断最大模值Mmax是否超过预设门限值:若是,则确定成功捕获当前卫星发送的B2a信号,且最大模值Mmax对应的半码片位置即为当前卫星发送的B2a信号的码相位;若否,则确定未成功捕获当前卫星发送的B2a信号。
其中,N表示码相位搜索范围内的半码片个数,i表示码相位搜索范围内的第i个半码片。
需要说明的是,在确定成功捕获当前卫星发送的B2a信号后,步骤S3计算得到的B2a信号的多普勒频率即为当前卫星发送的B2a信号的多普勒频率。
进一步的,完成对当前卫星的捕获后,即可重新指定卫星,重复执行步骤S2-S7对其进行捕获,直至完成对全部卫星的捕获。
基于B1C信号和B2a信号之间的特性——二者频率存在一定的比例关系且比特边界沿基本对齐,本发明实施例提供的B2a信号捕获方法通过利用B1C信号的捕获结果辅助实现B2a信号的捕获。具体来说,即是先捕获得到当前卫星发送的B1C信号对应的多普勒频率和码相位,进而根据B1C信号对应的多普勒频率和码相位确定B2a信号的码相位起始沿,进而可以在较小范围内对B2a码相位进行捕获搜索,因此本发明实施例提供的B2a信号捕获方法能够减少计算量,提高效率。此外,现有的码相位并行捕获方法最多只能得到1毫秒的积分结果,而基于本发明实施例提供的B2a信号捕获方法能够获得1毫秒的积分结果,因此相比现有技术,本发明实施例提供的B2a信号捕获方法能够提高捕获灵敏度。经验证,本发明实施例提供的B2a信号捕获方法将捕获灵敏度提高了7db左右。综上所述,本发明实施例提供的B2a信号捕获方法能够解决由于B2a信号码周期短、积分时长不够所导致的现有捕获方法捕获灵敏度不足的问题,快速捕获B2a信号,且简单易行,灵敏度较高。
基于上述B2a信号捕获方法,本发明实施例还提供了一种B2a信号捕获装置30,如图4所示。
参见图4,本发明实施例提供的B2a信号捕获装置包括:数字中频信号获取单元10、B1C信号捕获单元20以及B2a信号捕获单元30。
其中,数字中频信号获取单元10,用于获取北斗天线接收到的北斗卫星信号,并对北斗卫星信号依次进行下变频处理和A/D转换,得到数字中频信号;其中,北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号。
B1C信号捕获单元20,用于根据数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位;根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率;根据当前卫星发送的B1C信号的码相位,确定B2a信号的码相位搜索范围。
B2a信号捕获单元30,用于利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获。
优选的,如图5所示,本发明实施例提供的B2a信号捕获装置中,B2a信号捕获单元具体包括:本地载波发生器301、混频器302、重采样模块303、伪随机码发生器304、相关器305、累加器306以及捕获检测模块307。
其中,本地载波发生器301,用于根据B2a信号的多普勒频率和中频频率,生成本地载波。
数字中频信号获取单元10,还用于将数字中频信号的起始位置依次移动至B2a信号的码相位搜索范围内的各半码片处,得到各半码片对应的卫星中频信号。
混频器302,用于将本地载波分别与各半码片对应的卫星中频信号进行混频以剥离载波。
重采样模块303,用于对剥离载波后的信号进行重采样,得到各半码片对应的重采样数据。
伪随机码发生器304,用于生成当前卫星对应的B2a信号的伪随机码。
相关器305,用于利用伪随机码对各半码片对应的重采样数据分别进行循环相关,得到各半码片对应的1ms相干积分结果。
累加器306,用于根据子码符号,分别对各半码片对应的1ms相干积分结果进行5ms的相干累加,得到各半码片对应的5ms相干累加数据。
捕获检测模块307,用于根据各半码片对应的5ms相干累加数据,检测是否成功捕获当前卫星发送的B2a信号,并在成功捕获当前卫星发送的B2a信号后,确定当前卫星发送的B2a信号的码相位。
进一步的,本发明实施例提供的B2a信号捕获装置中,B1C捕获单元具体用于:
根据当前卫星发送的B1C信号的多普勒频率,利用第一预设公式,计算当前卫星发送的B2a信号的多普勒频率。
其中,第一预设公式为前述的公式1。
更进一步的,本发明实施例提供的B2a信号捕获装置中,B1C捕获单元具体还用于:
根据当前卫星发送的B1C信号的码相位,利用第二预设公式估算当前卫星发送的B2a信号的比特边界值;
根据B2a信号的码相位边界值,确定B2a信号的码相位搜索范围。
其中,第二预设公式为前述的公式2。
优选的,本发明实施例提供的B2a信号捕获装置中,B2a信号的码相位搜索范围为以B2a信号的码相位边界值为中心的左右各32个半码片构成的码相位范围。
一种具体实现方式中,在图2所示的B2a信号捕获装置中,捕获检测模块307具体用于:
对N组相干累加数据进行取模操作,得到对应的N个模值M1,M2,…,Mi,…,MN,确定N个模值M1,M2,…,Mi,…,MN中的最大模值Mmax
判断最大模值Mmax是否超过预设门限值:若是,则确定成功捕获当前卫星发送的B2a信号,且最大模值Mmax对应的半码片位置即为当前卫星发送的B2a信号的码相位;若否,则确定未成功捕获当前卫星发送的B2a信号。
其中,N表示码相位搜索范围内的半码片个数,i表示码相位搜索范围内的第i个半码片。
基于B1C信号和B2a信号之间的特性一一二者频率存在一定的比例关系且比特边界沿基本对齐,本发明实施例提供的B2a信号捕获装置通过利用B1C信号的捕获结果辅助实现B2a信号的捕获。具体来说,即是先捕获得到当前卫星发送的B1C信号对应的多普勒频率和码相位,进而根据B1C信号对应的多普勒频率和码相位确定B2a信号的码相位起始沿,进而可以在较小范围内对B2a码相位进行捕获搜索,因此本发明实施例提供的B2a信号捕获装置能够减少计算量,提高效率。此外,现有捕获装置最多只能得到1毫秒的积分结果,而基于本发明实施例提供的B2a信号捕获装置能够获得1毫秒的积分结果,因此相比现有技术,本发明实施例提供的B2a信号捕获装置能够提高捕获灵敏度。经验证,本发明实施例提供的B2a信号捕获装置将捕获灵敏度提高了7db左右。综上所述,本发明实施例提供的B2a信号捕获装置能够解决由于B2a信号码周期短、积分时长不够所导致的现有技术捕获灵敏度不足的问题,快速捕获B2a信号,且简单易行,灵敏度较高。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完戍前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (10)

1.一种B2a信号捕获方法,其特征在于,包括:
步骤S1,获取北斗天线接收到的北斗卫星信号,并对所述北斗卫星信号依次进行下变频处理和模拟/数字A/D转换,得到数字中频信号;其中,所述北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号;
步骤S2,根据所述数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位;
步骤S3,利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获;
具体步骤如下:
步骤S3a,根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率;根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围;
所述根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围,包括:
根据当前卫星发送的B1C信号的码相位,利用第二预设公式估算当前卫星发送的B2a信号的比特边界值;其中,所述第二预设公式包括:式中,n表示当前卫星发送的B2a信号的比特边界值,L表示主码长度,pB1C表示当前卫星发送的B1C信号的码相位,fs表示采样频率,fCA表示码速率;
根据所述B2a信号的码相位边界值,确定所述B2a信号的码相位搜索范围;
步骤S3b,根据所述B2a信号的多普勒频率和中频频率,生成本地载波;
步骤S3c,将所述数字中频信号的起始位置依次移动至所述B2a信号的码相位搜索范围内的各半码片处,得到各半码片对应的卫星中频信号;再将所述本地载波分别与各半码片对应的卫星中频信号进行混频以剥离载波,进而对剥离载波后的信号进行重采样,得到各半码片对应的重采样数据;
步骤S3d,生成当前卫星对应的B2a信号的伪随机码,利用所述伪随机码对各半码片对应的重采样数据分别进行循环相关,得到各半码片对应的1ms相干积分结果;根据子码符号,分别对各半码片对应的1ms相干积分结果进行5ms的相干累加,得到各半码片对应的5ms相干累加数据;
步骤S3e,根据各半码片对应的5ms相干累加数据,检测是否成功捕获当前卫星发送的B2a信号,并在成功捕获当前卫星发送的B2a信号后,确定当前卫星发送的B2a信号的码相位。
2.根据权利要求1所述的方法,其特征在于,步骤S3a中,所述根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率,包括:
根据当前卫星发送的B1C信号的多普勒频率,利用第一预设公式,计算当前卫星发送的B2a信号的多普勒频率;
其中,所述第一预设公式包括:式中,fdop_B2a表示当前卫星发送的B2A信号的多普勒频率,fB1C表示B1C信号的载波频率,fB2a表示B2a信号的载波频率,fdop_B1C表示当前卫星发送的B1C信号的多普勒频率。
3.根据权利要求1所述的方法,其特征在于,所述B2a信号的码相位搜索范围为以所述B2a信号的码相位边界值为中心的左右各32个半码片构成的码相位范围。
4.根据权利要求1所述的方法,其特征在于,步骤S7包括:
对所述N组相干累加数据进行取模操作,得到对应的N个模值M1,M2,…,Mi,…,MN,确定所述N个模值M1,M2,…,Mi,…,MN中的最大模值Mmax;其中,N表示码相位搜索范围内的半码片个数,i表示码相位搜索范围内的第i个半码片;
判断所述最大模值Mmax是否超过预设门限值:若是,则确定成功捕获当前卫星发送的B2a信号,且所述最大模值Mmax对应的半码片位置即为当前卫星发送的B2a信号的码相位;若否,则确定未成功捕获当前卫星发送的B2a信号。
5.一种B2a信号捕获装置,其特征在于,包括:数字中频信号获取单元、B1C信号捕获单元以及B2a信号捕获单元;
其中,所述数字中频信号获取单元,用于获取北斗天线接收到的北斗卫星信号,并对所述北斗卫星信号依次进行下变频处理和A/D转换,得到数字中频信号;其中,所述北斗卫星信号中包含至少一颗北斗三代卫星发送的B1C信号和B2a信号;
所述B1C信号捕获单元,用于根据所述数字中频信号进行捕获,得到当前卫星发送的B1C信号的多普勒频率和码相位;根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率;根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围;
所述B2a信号捕获单元,用于利用当前卫星发送的B1C信号的多普勒频率和码相位,进行B2a信号捕获;
具体为:根据当前卫星发送的B1C信号的多普勒频率,计算当前卫星发送的B2a信号的多普勒频率;根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围;
所述根据当前卫星发送的B1C信号的码相位,确定所述B2a信号的码相位搜索范围,包括:
根据当前卫星发送的B1C信号的码相位,利用第二预设公式估算当前卫星发送的B2a信号的比特边界值;其中,所述第二预设公式包括:式中,n表示当前卫星发送的B2a信号的比特边界值,L表示主码长度,pB1C表示当前卫星发送的B1C信号的码相位,fs表示采样频率,fCA表示码速率;
根据所述B2a信号的码相位边界值,确定所述B2a信号的码相位搜索范围;
根据所述B2a信号的多普勒频率和中频频率,生成本地载波;
将所述数字中频信号的起始位置依次移动至所述B2a信号的码相位搜索范围内的各半码片处,得到各半码片对应的卫星中频信号;再将所述本地载波分别与各半码片对应的卫星中频信号进行混频以剥离载波,进而对剥离载波后的信号进行重采样,得到各半码片对应的重采样数据;
生成当前卫星对应的B2a信号的伪随机码,利用所述伪随机码对各半码片对应的重采样数据分别进行循环相关,得到各半码片对应的1ms相干积分结果;根据子码符号,分别对各半码片对应的1ms相干积分结果进行5ms的相干累加,得到各半码片对应的5ms相干累加数据;
根据各半码片对应的5ms相干累加数据,检测是否成功捕获当前卫星发送的B2a信号,并在成功捕获当前卫星发送的B2a信号后,确定当前卫星发送的B2a信号的码相位。
6.根据权利要求5所述的装置,其特征在于,所述B2a信号捕获单元包括:本地载波发生器、混频器、重采样模块、伪随机码发生器、相关器、累加器以及捕获检测模块;其中,所述本地载波发生器,用于根据所述B2a信号的多普勒频率和中频频率,生成本地载波;
所述数字中频信号获取单元,还用于将所述数字中频信号的起始位置依次移动至所述B2a信号的码相位搜索范围内的各半码片处,得到各半码片对应的卫星中频信号;
所述混频器,用于将所述本地载波分别与各半码片对应的卫星中频信号进行混频以剥离载波;
所述重采样模块,用于对剥离载波后的信号进行重采样,得到各半码片对应的重采样数据;
所述伪随机码发生器,用于生成当前卫星对应的B2a信号的伪随机码;
所述相关器,用于利用所述伪随机码对各半码片对应的重采样数据分别进行循环相关,得到各半码片对应的1ms相干积分结果;
所述累加器,用于根据子码符号,分别对各半码片对应的1ms相干积分结果进行5ms的相干累加,得到各半码片对应的5ms相干累加数据;
所述捕获检测模块,用于根据各半码片对应的5ms相干累加数据,检测是否成功捕获当前卫星发送的B2a信号,并在成功捕获当前卫星发送的B2a信号后,确定当前卫星发送的B2a信号的码相位。
7.根据权利要求5所述的装置,其特征在于,所述B1C信号捕获单元具体用于:
根据当前卫星发送的B1C信号的多普勒频率,利用第一预设公式,计算当前卫星发送的B2a信号的多普勒频率;
其中,所述第一预设公式包括:式中,fdop_B2a表示当前卫星发送的B2A信号的多普勒频率,fB1C表示B1C信号的载波频率,fB2a表示B2a信号的载波频率,fdop_B1C表示当前卫星发送的B1C信号的多普勒频率。
8.根据权利要求5或6所述的装置,其特征在于,所述B1C信号捕获单元具体还用于:
根据当前卫星发送的B1C信号的码相位,利用第二预设公式估算当前卫星发送的B2a信号的比特边界值;其中,所述第二预设公式包括:式中,n表示当前卫星发送的B2a信号的比特边界值,L表示主码长度,pB1C表示当前卫星发送的B1C信号的码相位,fs表示采样频率,fCA表示码速率;
根据所述B2a信号的码相位边界值,确定所述B2a信号的码相位搜索范围。
9.根据权利要求8所述的装置,其特征在于,所述B2a信号的码相位搜索范围为以所述B2a信号的码相位边界值为中心的左右各32个半码片构成的码相位范围。
10.根据权利要求6所述的装置,其特征在于,所述捕获检测模块具体用于:
对所述N组相干累加数据进行取模操作,得到对应的N个模值M1,M2,…,Mi,…,MN,确定所述N个模值M1,M2,…,Mi,…,MN中的最大模值Mmax;其中,N表示码相位搜索范围内的半码片个数,i表示码相位搜索范围内的第i个半码片;
判断所述最大模值Mmax是否超过预设门限值:若是,则确定成功捕获当前卫星发送的B2a信号,且所述最大模值Mmax对应的半码片位置即为当前卫星发送的B2a信号的码相位;若否,则确定未成功捕获当前卫星发送的B2a信号。
CN201811469341.8A 2018-11-28 2018-11-28 一种B2a信号捕获方法及装置 Active CN109581434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811469341.8A CN109581434B (zh) 2018-11-28 2018-11-28 一种B2a信号捕获方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811469341.8A CN109581434B (zh) 2018-11-28 2018-11-28 一种B2a信号捕获方法及装置

Publications (2)

Publication Number Publication Date
CN109581434A CN109581434A (zh) 2019-04-05
CN109581434B true CN109581434B (zh) 2023-08-18

Family

ID=65925878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811469341.8A Active CN109581434B (zh) 2018-11-28 2018-11-28 一种B2a信号捕获方法及装置

Country Status (1)

Country Link
CN (1) CN109581434B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110068841A (zh) * 2019-05-06 2019-07-30 西安开阳微电子有限公司 一种卫星信号高灵敏度快速捕获方法
CN112578412A (zh) * 2020-11-24 2021-03-30 中国电子科技集团公司第二十研究所 一种兼容B1C信号与B2a信号的捕获方法
CN112558120A (zh) * 2020-12-01 2021-03-26 无锡奇芯科技有限公司 一种提高北斗b3i捕获跟踪灵敏度的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2701471A1 (en) * 2010-04-27 2011-10-27 Ecole De Technologie Superieure Universal acquisition and tracking apparatus for global navigation satellite system (gnss) signals
KR20130011129A (ko) * 2011-07-20 2013-01-30 한국전자통신연구원 Gps 신호 획득 및 추적 장치와 그 방법
CN104536016A (zh) * 2014-11-05 2015-04-22 北京大学 一种gnss新体制信号捕获装置及方法
CN108267755A (zh) * 2017-12-19 2018-07-10 北京遥测技术研究所 一种导航信号的跟踪处理方法
CN108519611A (zh) * 2018-03-01 2018-09-11 上海交通大学 基于gpu的北斗b1c/b1i双频并行多通道协同捕获方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2701471A1 (en) * 2010-04-27 2011-10-27 Ecole De Technologie Superieure Universal acquisition and tracking apparatus for global navigation satellite system (gnss) signals
KR20130011129A (ko) * 2011-07-20 2013-01-30 한국전자통신연구원 Gps 신호 획득 및 추적 장치와 그 방법
CN104536016A (zh) * 2014-11-05 2015-04-22 北京大学 一种gnss新体制信号捕获装置及方法
CN108267755A (zh) * 2017-12-19 2018-07-10 北京遥测技术研究所 一种导航信号的跟踪处理方法
CN108519611A (zh) * 2018-03-01 2018-09-11 上海交通大学 基于gpu的北斗b1c/b1i双频并行多通道协同捕获方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曾超.GNSS接收机快速捕获算法研究.《中国优秀硕士学位论文全文数据库 信息科技辑》.2014,第I136-I237页. *

Also Published As

Publication number Publication date
CN109581434A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
US6407699B1 (en) Method and device for rapidly extracting time and frequency parameters from high dynamic direct sequence spread spectrum radio signals under interference
CN109581434B (zh) 一种B2a信号捕获方法及装置
CN110045397B (zh) 一种基于fpga的l5信号捕获方法及装置
CN104536016B (zh) 一种gnss新体制信号捕获装置及方法
US20010055330A1 (en) Method and apparatus for code phase tracking
US8462616B2 (en) Apparatus and method for estimating a frequency shift and a time shift
CN116482727B (zh) 导航信号跟踪方法、装置及设备、芯片
CN107085222A (zh) 一种 boc 信号捕获方法和卫星导航接收机
CN108897009B (zh) 一种boc导航信号接收机及其码跟踪方法
CN109581433B (zh) 一种l5信号捕获方法、装置及计算机存储介质
CN102854516A (zh) 一种gnss接收机中的载噪比估计方法及其系统
CN109917429B (zh) 一种b1c弱信号的捕获方法、装置及计算机存储介质
US20110050496A1 (en) Energy Domain Based Peak Reconstruction Methods And Apparatuses
CN111665528B (zh) 卫星信号捕捉方法、装置及存储介质
Ye et al. CW interference effects on the performance of GPS receivers
Fortin et al. Implementation strategies for a software‐compensated FFT‐based generic acquisition architecture with minimal FPGA resources
KR20150056279A (ko) 사인 위상 BOC(n,n) 신호 추적 성능 향상을 위한 주변 첨두 제거 방법
JP5321500B2 (ja) 信号捕捉方法
CN110888150B (zh) 一种gnss接收机的抑制多径的方法及系统
CN108627863B (zh) 基于dft的gnss导航电文数据解调方法、装置及gnss接收终端
Ratnam et al. Acquisition of GPS L1 signals using Cooley-tukey FFT algorithm
CN112764063A (zh) 一种实现捕获处理的方法及接收机
CN109633707B (zh) 一种基于预平均处理的变系数匹配滤波的快速捕获方法
CN108226968B (zh) 一种导航信号快速捕获方法
Wei et al. Simulation and analysis of GPS software receiver

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