CN110007322A - 基于相干降采样的北斗b1i信号捕获方法 - Google Patents
基于相干降采样的北斗b1i信号捕获方法 Download PDFInfo
- Publication number
- CN110007322A CN110007322A CN201910306689.3A CN201910306689A CN110007322A CN 110007322 A CN110007322 A CN 110007322A CN 201910306689 A CN201910306689 A CN 201910306689A CN 110007322 A CN110007322 A CN 110007322A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- sampled
- concerned
- code
- 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.)
- Granted
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/246—Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
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
本发明提供一种基于相干降采样的北斗B1I信号捕获方法,通过相干降采样的采样方式减少了参与FFT的点数,使采样频率不再局限在固定的频率点上,减少了捕获算法的计算复杂度,结构简单,在不影响信噪比的同时有效提高了捕获效率,能够实现北斗B1I信号的快速捕获。
Description
技术领域
本发明属于卫星导航技术领域,涉及一种卫星信号的捕获方法,具体涉及一种基于相干降采样的北斗B1I信号捕获方法,适用于北斗B1I信号的捕获。
背景技术
北斗卫星导航系统是中国自行研制的全球卫星导航系统,致力于为全球用户提供稳定可靠的卫星导航服务。信号捕获是北斗系统软件接收机中的关键技术,用于搜索所有可见卫星并获得相应卫星的伪码相位和多普勒频移的粗略值。
传统的捕获方法有串行频率捕获法、并行码相位捕获法和并行频率捕获法三种。三种算法的核心思想都是利用北斗伪随机码的强自相关性实现对信号的捕获。串行频率捕获法利用时域相关器对信号的载波频率和伪码相位进行扫描式的二维搜索,算法不具有并行性,因为对于每一个伪码相位和载波频率的组合,其积分时间是固定的,且伪码相位和载波频率的组合有很多种,所以捕获时间很长;并行频率捕获法通过FFT(Fast FourierTransformation,快速傅氏变换)运算一次就完成对一个码单元内所有频率的搜索,所以只需要对伪码相位进行遍历,算法的并行性不够高,需要循环遍历的次数多,捕获效率低;并行码相位法是基于FFT和信号卷积的原理,将对码相位和载波频率的搜索转换为仅对频率的搜索,只对载波频率进行遍历。该算法是三种传统算法中并行性最强、捕获效率最高的方法,但是其结构复杂,采样频率的选择有限且参与FFT的点数相对较多,影响了卫星信号的捕获速度。
传统的并行频率捕获算法在实现了码相位剥离后,仍以原始的采样频率对信号进行采样。但是被剥离了码相位后的卫星信号带宽减小到10kHz左右,仍保持较高的采样频率会使参与FFT运算的数据量非常大,不利于信号的快速捕获。因此本文提出了一种新的捕获算法,在输入中频信号剥离伪码和载波后通过相干降采样的方式对其进行采样,使参与FFT的点数减少,不影响信号信噪比的同时提高了捕获效率。
发明内容
针对上述问题,本发明提出一种基于相干降采样的北斗B1I信号捕获方法,将其应用于软件接收机中,能够在节约资源的同时有效提高捕获速度。
本发明是采用如下的技术方案实现的:一种基于相干降采样的北斗信号捕获算法,包括以下步骤:
步骤一:软件接收机读取经过射频前端电路滤波、放大处理后的数字中频信号;
步骤二:设定捕获卫星的PRN码(pseudo random noise code,伪随机噪声码),将该卫星作为捕获对象;
步骤三:生成本地载波信号,将数字中频信号与本地载波信号混频;
步骤四:根据设定的PRN码生成相应的本地伪码信号,将步骤三得到的混频信号与本地伪码信号在某相位处对齐并进行相关累加;
步骤五:对累加结果相干降采样;
步骤六:对相干降采样结果进行离散傅里叶变换,记录峰值能量;
步骤七:将本地伪码信号移1位后重复执行步骤四到步骤六。由于北斗B1I信号的伪随机码的码速率为2.046MHz,码长为2046个码片,所以如此循环2046次,可以遍历完所有的伪码相位并得到2046个峰值能量。遍历得到的所有峰值,将峰值结果中的最大值作为检测值;
步骤八:将检测值和设定的门限值进行比较,若高于门限值则捕获成功,否则执行步骤九;
步骤九:改变本地载波信号频率,重复执行步骤三到步骤八,直到捕获到卫星信号或搜索完所有可能的频率点。
步骤十:若在所有的频率点上都捕获不到信号,改变设定的PRN码,重复执行步骤二到步骤九,直到捕获完所有可见卫星。
优选的,步骤一中,数字中频信号表达式为:
其中,n表示采样后得到的离散信号的数据点中的第n个采样数据;j是PRN码;Aj是信号振幅;Cj是C/A码;Dj是导航电文比特;fIF是载波中心频率;是多普勒频移;是初始载波相位;w(n)是高斯白噪声。
优选的,步骤三中,将数字中频信号sj(n)与本地载波信号SIF混频,过程表示如下:
其中,fIF表示载波中心频率;本地载波信号SIF=exp[-i(2πfIF)n];i表示虚数单位。
优选的,步骤四中,将步骤三得到的混频信号和本地伪码信号C(n-m)在某相位处对齐并进行相关累加,过程表示如下:
m代表本地伪码的码相位,m=0,1,2...2045。
优选的,步骤五中,对累加结果相干降采样结果表示如下:
k表示相干降采样结果的第k个数据点,M表示相干降采样累加点数。
优选的,步骤六中,离散傅里叶变换的步骤包括:
首先对无限长序列xj(k)进行加窗处理,得到信号该过程表示为:
其中,Ts表示采样周期,L表示FFT点数;
根据时域相乘等于频域卷积原理,将信号转为频域信号X(f),再对频域信号X(f)进行离散加窗处理得到信号X′(kk),过程表示为:
对信号X′(kk)取模,得到取模结果为:
将得到的信号X(kk)的峰值作为检测值;
其中,f表示将时域信号变到频域信号之后的频率,kk表示对频域信号离散加窗之后得到的信号的数据点中的第kk个数据,Δf表示采样每一小格的频率值,分别表示I支路和Q支路的分量。
有益效果:本发明通过相干降采样的采样方式减少了参与DFT的点数,使采样频率不再局限在固定的频率点上,减少了捕获算法的计算复杂度,结构简单,在不影响信噪比的同时有效提高了捕获效率,能够实现北斗B1I信号的快速捕获。
附图说明
图1是基于相干降采样的北斗B1I信号捕获算法的结构示意图。
图2是本发明中相干降采样模块的示意图。
图3是基于相干降采样的北斗B1I信号捕获算法的流程框图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
基于相干降采样的并行频率捕获算法的原理图如图1所示。具体实施时,接收天线可以将北斗卫星发射的电磁波信号转换为电压信号或电流信号,随后将信号送入射频前端进行处理。在经过射频前端模块中前置滤波器和前置放大器的滤波放大后,信号变为频率较低、载噪比较高的中频信号,再经过A/D转换器的采样,得到的数字中频信号表达式为:
其中,n表示采样后得到的离散信号的数据点中的第n个采样数据;j是PRN码;Aj是信号振幅;Cj是C/A码;Dj是导航电文比特;fIF是载波中心频率;是多普勒频移;是初始载波相位;w(n)是高斯白噪声。对该信号的处理过程为:
第一步:将数字中频信号与本地载波信号相乘,实现对信号的下变频,过程表示如下:
其中,fIF表示载波中心频率,本地载波信号SIF=exp[-i(2πfIF)n];i表示虚数单位。
第二步:将得到的零中频数字信号和本地伪码信号C(n-m)相乘,完成伪码的剥离,过程表示如下:
在信号捕获的过程中,可以通过改变m值来调整本地生成的伪码相位。当m=0时,本地伪码和输入信号的相位对齐,即可实现伪码剥离。m代表本地伪码的码相位,因为北斗B1信号码长为2046个码片,m=0,1,2...2045。C表示本地伪码序列。
第三步:对已经实现载波剥离和伪码剥离后的信号进行相干降采样,将采样频率减小为原始采样频率的1/M。首先将采集到的数据点数补零至M的整数倍,然后对数据进行相干降采样,得到的降采样结果表示如下:
k表示相干降采样结果的第k个数据点,M表示相干降采样累加点数。
相干降采样的运算过程如图2所示。以处理采集4092个数据点为例,首先将数据补零至4096个数据点,然后以M=64对原始信号进行相干降采样,即可得到64点相干降采样结果。
第四步:对降采样结果进行离散傅里叶变换得到峰值,并将结果取模作为检测值;
首先对无限长序列xj(k)进行加窗处理,得到信号该过程表示为:
其中,Ts表示采样周期,L表示FFT点数。
根据时域相乘等于频域卷积原理,将信号转为频域信号X(f),再对频域信号X(f)进行离散加窗处理得到信号X′(kk),过程表示为:
对信号X′(kk)取模,得到取模结果为:
将得到的信号X(kk)的峰值作为检测值。
其中,f表示把时域信号变到频域信号之后的频率,kk表示对频域信号离散加窗之后得到的信号的数据点中的第kk个数据,Δf表示采样每一小格的频率值,分别表示I支路和Q支路的分量。
第五步:将检测值与之前设定的门限值进行比较。若结果大于设定的门限值,说明捕获成功。若没有通过门限判决,那么改变本地伪码相位重新进行检测。当检测完全部2046个码片后,若还没有通过门限判决的结果,那么改变本地载波频率重新进行检测。若搜索完全部频率点仍没有大于门限值的结果,那么改变PRN码,进行下一颗星的搜索。
基于相干降采样的并行频率捕获算法的流程框图如图3所示,包括如下步骤:
步骤一:软件接收机读取经过射频前端电路的滤波、放大处理后的数字中频信号;
步骤二:设定捕获卫星的PRN码,将该卫星作为捕获对象。
步骤三:生成本地载波信号,将数字中频信号与本地载波信号混频,通过数字下变频模块将数字中频信号变为零中频数字信号。
步骤四:根据设定的PRN码生成相应的本地伪码信号,将零中频数字信号累积降速至本地伪码信号的频率,再与本地伪码信号在某相位处对齐并进行相关累加。
步骤五:对累加结果相干降采样。
步骤六:对相干降采样结果进行离散傅里叶变换,记录峰值能量。
步骤七:将本地伪码信号移1位后重复执行步骤四到步骤六。由于北斗B1I信号的伪随机码的码速率为2.046MHz,码长为2046个码片,所以如此循环2046次,可以遍历完所有的伪码相位并得到2046个峰值能量。遍历得到的所有峰值,将峰值结果中的最大值作为检测值。图中的count是循环次数。
步骤八:将检测值和设定的门限值进行比较,若高于门限值则捕获成功,否则执行步骤九。
步骤九:改变本地载波信号频率,重复执行步骤三到步骤八,直到捕获到卫星信号或搜索完所有可能的频率点。
步骤十:若在所有的频率点上都捕获不到信号,改变设定的PRN码,重复执行步骤二到步骤九,直到捕获完所有可见卫星。
除上述实施例外,本发明还可以有其他实施方式。凡本技术领域的普通技术人员在不脱离本发明原理的前提下,通过改进和润饰所形成的等同替换或等效变换,均落在本发明要求的保护范围。
Claims (6)
1.基于相干降采样的北斗B1I信号捕获方法,其特征在于,包括如下步骤:
步骤一:软件接收机读取经过射频前端电路滤波、放大处理后的数字中频信号;
步骤二:设定捕获卫星的PRN码,将该卫星作为捕获对象;
步骤三:生成本地载波信号,将数字中频信号与本地载波信号混频;
步骤四:根据设定的PRN码生成相应的本地伪码信号,将步骤三得到的混频信号与本地伪码信号在某相位处对齐并进行相关累加;
步骤五:对累加结果相干降采样;
步骤六:对相干降采样结果进行离散傅里叶变换,记录峰值能量;
步骤七:将本地伪码信号移1位后重复执行步骤四到步骤六,如此循环2046次,遍历完所有的伪码相位并得到2046个峰值能量,将峰值结果中的最大值作为检测值;
步骤八:将检测值和设定的门限值进行比较,若高于门限值则捕获成功,否则执行步骤九;
步骤九:改变本地载波信号频率,重复执行步骤三到步骤八,直到捕获到卫星信号或搜索完所有可能的频率点。
步骤十:若在所有的频率点上都捕获不到信号,改变设定的PRN码,重复执行步骤二到步骤九,直到捕获完所有可见卫星。
2.根据权利要求1所述的基于相干降采样的北斗B1I信号捕获方法,其特征在于,步骤一中,数字中频信号表达式为:
其中,n表示采样后得到的离散信号的数据点中的第n个采样数据;j是PRN码;Aj是信号振幅;Cj是C/A码;Dj是导航电文比特;fIF是载波中心频率;是多普勒频移;是初始载波相位;w(n)是高斯白噪声。
3.根据权利要求2所述的基于相干降采样的北斗B1I信号捕获方法,其特征在于,步骤三中,将数字中频信号sj(n)与本地载波信号SIF混频,过程表示如下:
其中,fIF表示载波中心频率;本地载波信号SIF=exp[-i(2πfIF)n];i表示虚数单位。
4.根据权利要求3所述的基于相干降采样的北斗B1I信号捕获方法,其特征在于,步骤四中,将步骤三得到的混频信号和本地伪码信号C(n-m)在某相位处对齐并进行相关累加,过程表示如下:
m代表本地伪码的码相位,m=0,1,2...2045。
5.根据权利要求4所述的基于相干降采样的北斗B1I信号捕获方法,其特征在于,步骤五中,对累加结果相干降采样结果表示如下:
k表示相干降采样结果的第k个数据点,M表示相干降采样累加点数。
6.根据权利要求5所述的基于相干降采样的北斗B1I信号捕获方法,其特征在于,步骤六中,离散傅里叶变换的步骤包括:
首先对无限长序列xj(k)进行加窗处理,得到信号该过程表示为:
其中,Ts表示采样周期,L表示FFT点数;
根据时域相乘等于频域卷积原理,将信号转为频域信号X(f),再对频域信号X(f)进行离散加窗处理得到信号X′(kk),过程表示为:
对信号X′(kk)取模,得到取模结果为:
将得到的信号X(kk)的峰值作为检测值;
其中,f表示将时域信号变到频域信号之后的频率,kk表示对频域信号离散加窗之后得到的信号的数据点中的第kk个数据,Δf表示采样每一小格的频率值,分别表示I支路和Q支路的分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306689.3A CN110007322B (zh) | 2019-04-17 | 2019-04-17 | 基于相干降采样的北斗b1i信号捕获方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306689.3A CN110007322B (zh) | 2019-04-17 | 2019-04-17 | 基于相干降采样的北斗b1i信号捕获方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110007322A true CN110007322A (zh) | 2019-07-12 |
CN110007322B CN110007322B (zh) | 2022-06-10 |
Family
ID=67172445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910306689.3A Active CN110007322B (zh) | 2019-04-17 | 2019-04-17 | 基于相干降采样的北斗b1i信号捕获方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110007322B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115061165A (zh) * | 2022-08-18 | 2022-09-16 | 中国人民解放军国防科技大学 | 基于低损耗降采样策略的导航扩频信号捕获方法及系统 |
CN115657093A (zh) * | 2022-12-29 | 2023-01-31 | 成都奇芯微电子有限公司 | 基于捕获数据存储的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102759739A (zh) * | 2011-04-26 | 2012-10-31 | 中国科学院微电子研究所 | 快速联合捕获装置及其方法 |
CN104360357A (zh) * | 2014-11-24 | 2015-02-18 | 四川九洲电器集团有限责任公司 | 基于循环方式的北斗卫星信号快速捕获方法及系统 |
CN105301608A (zh) * | 2015-10-30 | 2016-02-03 | 西安烽火电子科技有限责任公司 | 高动态接收机及其对北斗卫星信号b1频点的捕获方法 |
CN106556848A (zh) * | 2016-10-24 | 2017-04-05 | 华中光电技术研究所(中国船舶重工集团公司第七七研究所) | 一种北斗二代卫星b1频点信号的快速捕获方法 |
CN106814375A (zh) * | 2017-01-24 | 2017-06-09 | 中国电子科技集团公司第五十四研究所 | 一种抗转发式欺骗的捕获方法和接收机 |
CN107153206A (zh) * | 2016-10-13 | 2017-09-12 | 南京航空航天大学 | 一种基于fft的北斗卫星信号快速捕获方法 |
-
2019
- 2019-04-17 CN CN201910306689.3A patent/CN110007322B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102759739A (zh) * | 2011-04-26 | 2012-10-31 | 中国科学院微电子研究所 | 快速联合捕获装置及其方法 |
CN104360357A (zh) * | 2014-11-24 | 2015-02-18 | 四川九洲电器集团有限责任公司 | 基于循环方式的北斗卫星信号快速捕获方法及系统 |
CN105301608A (zh) * | 2015-10-30 | 2016-02-03 | 西安烽火电子科技有限责任公司 | 高动态接收机及其对北斗卫星信号b1频点的捕获方法 |
CN107153206A (zh) * | 2016-10-13 | 2017-09-12 | 南京航空航天大学 | 一种基于fft的北斗卫星信号快速捕获方法 |
CN106556848A (zh) * | 2016-10-24 | 2017-04-05 | 华中光电技术研究所(中国船舶重工集团公司第七七研究所) | 一种北斗二代卫星b1频点信号的快速捕获方法 |
CN106814375A (zh) * | 2017-01-24 | 2017-06-09 | 中国电子科技集团公司第五十四研究所 | 一种抗转发式欺骗的捕获方法和接收机 |
Non-Patent Citations (1)
Title |
---|
曾庆喜等: "基于分段折叠算法的北斗的CB21码快速捕获", 《中国惯性技术学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115061165A (zh) * | 2022-08-18 | 2022-09-16 | 中国人民解放军国防科技大学 | 基于低损耗降采样策略的导航扩频信号捕获方法及系统 |
CN115657093A (zh) * | 2022-12-29 | 2023-01-31 | 成都奇芯微电子有限公司 | 基于捕获数据存储的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110007322B (zh) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100747B (zh) | 一种适用于多路随机接入信号的快速捕获系统及方法 | |
CN105005057B (zh) | 一种北斗导航系统d1导航电文的捕获方法 | |
CN103441777B (zh) | 一种星载扩频接收机及其提高载波捕获精度的方法 | |
CN109061693B (zh) | 一种适用于p码捕获的改进pmf-fft方法 | |
CN110071738A (zh) | 基于多伪码支路接收机的扩频信号解扩跟踪方法 | |
CN105743612B (zh) | 实时盲解调大频移短时突发信号的方法 | |
CN110007322A (zh) | 基于相干降采样的北斗b1i信号捕获方法 | |
CN105553506B (zh) | 一种长码扩频信号的快速捕获方法及装置 | |
CN108169772A (zh) | 一种加窗fft的卫星信号捕获方法 | |
CN103941269A (zh) | 用于卫星导航系统的pn码捕获方法 | |
CN105974448A (zh) | 捕获卫星信号的方法和装置 | |
CN103197328B (zh) | 一种基于2ms数据存储的高动态卫星接收机快速捕获方法 | |
CN109782310A (zh) | 一种BOC(n,n)调制信号快速捕获方法 | |
CN107247276B (zh) | 基于重叠多块补零算法的弱信号捕获方法 | |
CN105242286A (zh) | 一种基于小波域滤波码捕获方法 | |
CN117270002B (zh) | 一种新体制北斗信号无模糊度快速捕获方法、系统及设备 | |
CN107907893B (zh) | 基于并行频率与并行码搜索的分段式可配置军码直捕方法 | |
CN103760578A (zh) | 一种gnss卫星导航信号的无模糊跟踪方法 | |
CN109581434A (zh) | 一种B2a信号捕获方法及装置 | |
CN102394669B (zh) | 一种用于高增益扩频通信系统的快速伪码捕获方法 | |
CN107621643A (zh) | 一种适用于导航信号质量评估的相关域参数精确解算方法 | |
CN109633707B (zh) | 一种基于预平均处理的变系数匹配滤波的快速捕获方法 | |
CN104297767B (zh) | 一种自适应的多次重采样导航卫星捕获方法及其实现系统 | |
CN103869343B (zh) | 捕获gps中l2c信号的方法及装置 | |
CN105846855A (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 |