CN113820731A - 一种卫星导航接收机无周期长码信号的捕获方法及系统 - Google Patents
一种卫星导航接收机无周期长码信号的捕获方法及系统 Download PDFInfo
- Publication number
- CN113820731A CN113820731A CN202111407144.5A CN202111407144A CN113820731A CN 113820731 A CN113820731 A CN 113820731A CN 202111407144 A CN202111407144 A CN 202111407144A CN 113820731 A CN113820731 A CN 113820731A
- Authority
- CN
- China
- Prior art keywords
- value
- envelope
- correlation
- phase
- accumulation
- 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
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/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)
- Radio Relay Systems (AREA)
Abstract
本发明公开了一种卫星导航接收机无周期长码信号的捕获方法及系统。方法包括以下步骤:S1、快速搜索阶段,使用较短的相干积分时间进行分段相关累加,对分段相关值进行快速傅里叶变换、取包络和后积累处理,对频谱包络取最大得到当前码相位下的多普勒频率初始估计值;S2、虚警验证阶段,对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,同样进行快速傅里叶变换、取包络和后积累处理,使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果。本发明在保证多普勒搜索范围的情况下,能实现较高的检测性能,在不增加相关器资源的情况下实现两级检测。
Description
技术领域
本发明涉及卫星导航技术领域,更具体地说,特别涉及一种卫星导航接收机无周期长码信号的捕获方法及系统。
背景技术
不同于采用短码的民用公开信号,各卫星导航系统的授权信号均使用无周期长码。通常情况下,卫星导航接收机在接收无周期长码时,需要首先接收民用信号,并从民用信号中解调出具体当前的准确时间,但是在某些情况下民用信号可能无法使用,或者为了防止被欺骗,授权信号接收机会对授权信号进行直接捕获。
目前导航信号接收机的捕获模块通常采用频域并行的实现架构,具体如图1所示,其中,Tc表示分段积分时间,M表示相干积累的分段数,N表示FFT点数,L表示后积累次数。为了能够在低虚警概率下实现较高的检测概率,就需要提高检测性能。在输入信号载噪比一定和约束总的积累时长的情况下,检测性能主要取决于分段积分时间Tc。增加分段积分时间可提高检测性能,但同时会减小多普勒频率搜索范围。当接收机动态范围较大时,就需要进行多次捕获,这会严重影响捕获时间。因此,在传统的捕获模块实现架构中,检测性能和多普勒搜索范围是难以兼顾的。
发明内容
本发明的目的在于提供一种卫星导航接收机无周期长码信号的捕获方法及系统,以克服现有技术所存在的缺陷。
为了达到上述目的,本发明采用的技术方案如下:
一种卫星导航接收机无周期长码信号的捕获方法,包括以下步骤:
S1、快速搜索阶段,首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值,并使用多普勒频率的初始估计值对分段相关值进行相位旋转;
S2、虚警验证阶段,首先对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
进一步地,所述步骤S1的搜索阶段具体包括以下步骤:
S10、定义T c1为搜索阶段的分段积分时间,M 1为搜索阶段的相干积累的分段数,将输入的基带复信号和本地码进行相关累加,得到第i次后积累第k个分段的相关值,其中每个分段的相关值的表达式为,表示当前待验证的码相位偏移量;
进一步地,所述步骤S2具体包括:
S23、对积累后的频域包络值取最大得到最大包络值,对应的序号为,将最大包络值与预设的捕获门限比较,若则改变,重复步骤S10~S22,继续下一个相位的搜索,若则判定捕获到信号,输出捕获结果,所述捕获结果包括码相位和多普勒频率,为虚警验证阶段估计得到残留多普勒频率;
进一步地,所述虚警验证阶段估计得到残留多普勒频率的计算公式为:
本发明还提供一种根据上述的卫星导航接收机无周期长码信号的捕获方法的系统,包括:
快速搜索模块,用于首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值,并使用多普勒频率的初始估计值对分段相关值进行相位旋转;
虚警验证模块,用于首先对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
与现有技术相比,本发明的优点在于:本发明在保证多普勒搜索范围的情况下,能够实现较高的检测性能,并且通过复用相关值,在不增加相关器资源的情况下实现两级检测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是现有技术的捕获模块通常采用频域并行的实现架构图。
图2是本发明卫星导航接收机无周期长码信号的捕获方法的流程图。
图3是本发明卫星导航接收机无周期长码信号的捕获系统的原理图。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
针对传统捕获模块在进行无周期长码信号的捕获时,难以兼顾检测性能和多普勒搜索范围的问题,本发明通过基于快速搜索和虚警验证两层处理的实现架构,同时为了降低相关值计算所需的硬件资源,通过缓存相关值在两层处理中实现相关值的复用。
本发明将传统的频域并行捕获过程分为快速搜索和虚警验证两个阶段,即:S1、快速搜索阶段,首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值;S2、虚警验证阶段,首先使用多普勒频率的初始估计值对快速搜索阶段的分段相关值进行相位旋转,并将相位旋转后的多个分段相关值进行累加得到相干积分时间较长的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决。当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
参阅图2所示,其中,T c1和T c2分别表示快速搜索和虚警验证阶段的分段积分时间,M 1和M 2分别表示快速搜索和虚警验证阶段相干积累的分段数,N表示FFT(快速傅里叶变换)点数,L 1和L 2分别表示快速搜索和虚警验证阶段的后积累次数。
以快速搜索阶段的设计参数为例,T c1决定了多普勒搜索范围为,N决定了多普勒频率的估计分辨率,T c1、M 1、N、L 1共同决定了快速搜索阶段的检测性能。这些参数可事先根据捕获模块的设计目标来确定。虚警验证阶段的设计参数T c1、M 1、N、L 1具有类似的作用。
虽然快速搜索和虚警验证两个阶段的设计参数相互间并无约束,但是当T c2=HT c1,M 1=HM 2,L 1= L 2时,两阶段的相关值可以实现最大程度的复用。
本实施例公开了一种卫星导航接收机无周期长码信号的捕获方法,具体包括以下步骤:
步骤S10、定义T c1为搜索阶段的分段积分时间,M 1为搜索阶段的相干积累的分段数,将输入的基带复信号和本地码进行相关累加,得到第i次后积累第k个分段的相关值,其中每个分段的相关值的表达式为,表示当前待验证的码相位偏移量;
所述的步骤S2具体包括:
步骤S23、对积累后的频域包络值取最大得到最大包络值,对应的序号为,将最大值与预设的捕获门限比较,若则改变,重复步骤S10~S22,继续下一个相位的搜索,若则判定捕获到信号,输出捕获结果,所述捕获结果包括码相位和多普勒频率,为虚警验证阶段估计得到残留多普勒频率。
所述虚警验证阶段估计得到残留多普勒频率的计算公式为:
结合图3所示,本发明还提供一种根据上述的卫星导航接收机无周期长码信号的捕获方法的系统,包括:快速搜索模块1,用于首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值,并使用多普勒频率的初始估计值对分段相关值进行相位旋转;虚警验证模块2,用于首先对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
下面通过具体实验对本发明作进一步说明。
下表1分别是输入信号载噪比32dBHz,FFT点数为32点时,采用传统方法和本发明的方法所对应的检测性能。
表1
由上述仿真结果可见,相比传统方法,当总的积累时长均为80ms,虚警概率为0.02%时,本发明的检测性能为93.5%,远优于传统方法的62.0%。传统方法为了实现与本发明相同的检测性能,所需的总积累时长需由80ms增加至140ms,也就是说本发明在保证多普勒搜索范围的情况下,能够实现较高的检测性能。
虽然结合附图描述了本发明的实施方式,但是专利所有者可以在所附权利要求的范围之内做出各种变形或修改,只要不超过本发明的权利要求所描述的保护范围,都应当在本发明的保护范围之内。
Claims (5)
1.一种卫星导航接收机无周期长码信号的捕获方法,其特征在于,包括以下步骤:
S1、快速搜索阶段,首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值,并使用多普勒频率的初始估计值对分段相关值进行相位旋转;
S2、虚警验证阶段,首先对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
2.根据权利要求1所述的卫星导航接收机无周期长码信号的捕获方法,其特征在于,所述步骤S1的搜索阶段具体包括以下步骤:
S10、定义T c1为搜索阶段的分段积分时间,M 1为搜索阶段的相干积累的分段数,将输入的基带复信号和本地码进行相关累加,得到第i次后积累第k个分段的相关值,其中每个分段的相关值的表达式为,表示当前待验证的码相位偏移量;
3.根据权利要求2所述的卫星导航接收机无周期长码信号的捕获方法,其特征在于,所述步骤S2具体包括:
5.根据权利要求1-4任意一项所述的卫星导航接收机无周期长码信号的捕获方法的系统,其特征在于,包括:
快速搜索模块,用于首先使用较短的相干积分时间进行分段相关累加,然后对分段相关值进行快速傅里叶变换、取包络和后积累处理,最后对频谱包络取最大得到当前码相位下的多普勒频率初始估计值,并使用多普勒频率的初始估计值对分段相关值进行相位旋转;
虚警验证模块,用于首先对快速搜索阶段相位旋转后的多个分段相关值进行累加得到更长相干积分时间的分段相关值,然后同样进行快速傅里叶变换、取包络和后积累处理,最后使用频谱包络最大值进行捕获判决,当包络最大值超过捕获门限时,输出当前的码相位和多普勒频率作为信号捕获结果,否则继续下一个码相位的搜索。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111407144.5A CN113820731B (zh) | 2021-11-24 | 2021-11-24 | 一种卫星导航接收机无周期长码信号的捕获方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111407144.5A CN113820731B (zh) | 2021-11-24 | 2021-11-24 | 一种卫星导航接收机无周期长码信号的捕获方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113820731A true CN113820731A (zh) | 2021-12-21 |
CN113820731B CN113820731B (zh) | 2022-06-10 |
Family
ID=78918203
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111407144.5A Active CN113820731B (zh) | 2021-11-24 | 2021-11-24 | 一种卫星导航接收机无周期长码信号的捕获方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113820731B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115113243A (zh) * | 2022-08-29 | 2022-09-27 | 中国人民解放军国防科技大学 | 一种联合频率插值估计的导航信号捕获方法 |
CN115267860A (zh) * | 2022-09-27 | 2022-11-01 | 中国人民解放军国防科技大学 | 一种高动态短突发信号的多相关器组高精度引导方法 |
CN118068384A (zh) * | 2024-04-17 | 2024-05-24 | 湖南跨线桥航天科技有限公司 | 一种基于卫星导航信号测量旋转载体转速的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520423A (zh) * | 2011-12-06 | 2012-06-27 | 成都金本华科技有限公司 | 卫星导航接收机长周期伪随机扩频码的快速捕获电路及其捕获方法 |
US20120250733A1 (en) * | 2011-03-31 | 2012-10-04 | General Dynamics Advanced Information Systems | Method and apparatus for rapid acquisitions of gps signals in space applications |
CN104993844A (zh) * | 2015-06-29 | 2015-10-21 | 和芯星通科技(北京)有限公司 | 一种频域搜索方法及装置 |
CN105049081A (zh) * | 2015-06-05 | 2015-11-11 | 重庆大学 | 一种适应高动态环境的长码扩频信号快速捕获方法 |
CN106253947A (zh) * | 2016-07-29 | 2016-12-21 | 西安电子科技大学 | 基于双层时频分析的低轨卫星直扩系统快速捕获算法 |
CN106772471A (zh) * | 2016-11-16 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 一种基于gpu的长码分段重叠局部相关捕获方法 |
CN112327334A (zh) * | 2020-09-29 | 2021-02-05 | 航天恒星科技有限公司 | 一种低轨卫星辅助的gnss长码信号捕获方法及装置 |
CN112910498A (zh) * | 2021-01-25 | 2021-06-04 | 上海航天电子通讯设备研究所 | Pmf-fft测控信号捕获装置和方法 |
CN112987047A (zh) * | 2021-05-13 | 2021-06-18 | 湖南跨线桥航天科技有限公司 | 基于时分复用技术的卫星导航信号捕获系统及方法 |
-
2021
- 2021-11-24 CN CN202111407144.5A patent/CN113820731B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120250733A1 (en) * | 2011-03-31 | 2012-10-04 | General Dynamics Advanced Information Systems | Method and apparatus for rapid acquisitions of gps signals in space applications |
CN102520423A (zh) * | 2011-12-06 | 2012-06-27 | 成都金本华科技有限公司 | 卫星导航接收机长周期伪随机扩频码的快速捕获电路及其捕获方法 |
CN105049081A (zh) * | 2015-06-05 | 2015-11-11 | 重庆大学 | 一种适应高动态环境的长码扩频信号快速捕获方法 |
CN104993844A (zh) * | 2015-06-29 | 2015-10-21 | 和芯星通科技(北京)有限公司 | 一种频域搜索方法及装置 |
CN106253947A (zh) * | 2016-07-29 | 2016-12-21 | 西安电子科技大学 | 基于双层时频分析的低轨卫星直扩系统快速捕获算法 |
CN106772471A (zh) * | 2016-11-16 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 一种基于gpu的长码分段重叠局部相关捕获方法 |
CN112327334A (zh) * | 2020-09-29 | 2021-02-05 | 航天恒星科技有限公司 | 一种低轨卫星辅助的gnss长码信号捕获方法及装置 |
CN112910498A (zh) * | 2021-01-25 | 2021-06-04 | 上海航天电子通讯设备研究所 | Pmf-fft测控信号捕获装置和方法 |
CN112987047A (zh) * | 2021-05-13 | 2021-06-18 | 湖南跨线桥航天科技有限公司 | 基于时分复用技术的卫星导航信号捕获系统及方法 |
Non-Patent Citations (4)
Title |
---|
XIAOPENG JIAO 等: "High sensitivity GPS acquisition algorithm based on code doppler compensation", 《2012 IEEE 11TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING》, 4 April 2013 (2013-04-04), pages 241 - 245 * |
刘天雄: "《卫星导航系统概论》", 30 November 2018, 中国宇航出版社, pages: 379 - 381 * |
苏悦 等: "GNSS弱信号两级快速傅里叶变换捕获方法", 《中国空间科学技术》 * |
苏悦 等: "GNSS弱信号两级快速傅里叶变换捕获方法", 《中国空间科学技术》, no. 4, 31 August 2014 (2014-08-31) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115113243A (zh) * | 2022-08-29 | 2022-09-27 | 中国人民解放军国防科技大学 | 一种联合频率插值估计的导航信号捕获方法 |
CN115113243B (zh) * | 2022-08-29 | 2022-11-15 | 中国人民解放军国防科技大学 | 一种联合频率插值估计的导航信号捕获方法 |
CN115267860A (zh) * | 2022-09-27 | 2022-11-01 | 中国人民解放军国防科技大学 | 一种高动态短突发信号的多相关器组高精度引导方法 |
CN115267860B (zh) * | 2022-09-27 | 2022-12-27 | 中国人民解放军国防科技大学 | 一种高动态短突发信号的多相关器组高精度引导方法 |
CN118068384A (zh) * | 2024-04-17 | 2024-05-24 | 湖南跨线桥航天科技有限公司 | 一种基于卫星导航信号测量旋转载体转速的方法 |
CN118068384B (zh) * | 2024-04-17 | 2024-07-30 | 湖南跨线桥航天科技有限公司 | 一种基于卫星导航信号测量旋转载体转速的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113820731B (zh) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113820731B (zh) | 一种卫星导航接收机无周期长码信号的捕获方法及系统 | |
CN106646546B (zh) | 一种卫星信号多维快速捕获方法及系统 | |
CN104280750B (zh) | 基于部分匹配滤波fft算法的长码捕获系统及方法 | |
JP5932728B2 (ja) | 無線マイクロフォン信号を検出するための無線センサ及びその方法 | |
CN113050130A (zh) | 卫星通信高动态突发信号的捕获方法及装置 | |
TW200806024A (en) | Detection of presence of television signals embedded in noise using cyclostationary toolbox | |
CN112118196B (zh) | 信号捕获方法及系统 | |
WO2017036325A1 (zh) | 信号同步的方法和信号同步的装置 | |
US20200036557A1 (en) | Transmitter and receiver and corresponding methods | |
CN113972929B (zh) | 一种高动态多普勒下扩频信号的捕获方法 | |
WO2021128782A1 (zh) | 一种nbiot的频域时频同步方法 | |
CN106788584B (zh) | 一种适应于长码字应用的改进pmf-fft捕获方法 | |
CN107907893B (zh) | 基于并行频率与并行码搜索的分段式可配置军码直捕方法 | |
CN116366092A (zh) | 一种多普勒捕获方法、装置和存储介质 | |
CN111624632B (zh) | 一种卫星导航信号捕获方法及装置 | |
CN112738000B (zh) | 一种pss分块互相关检测方法 | |
CN105891853A (zh) | 信号捕获方法及装置 | |
CN103561274A (zh) | 静止摄像头拍摄的运动目标被移除视频时域篡改检测方法 | |
CN111490954B (zh) | 信道脉冲响应的重要时延抽头选择方法及系统 | |
KR100983502B1 (ko) | 직교 주파수 분할 다중 시스템에서의 주파수 오차 검출방법 및 장치 | |
CN114584227B (zh) | 自动化突发信号检测方法 | |
CN108226968B (zh) | 一种导航信号快速捕获方法 | |
CN114944974A (zh) | 一种频偏估计方法、装置及电子设备 | |
US7660371B2 (en) | Normalized auto-correlators | |
CN109633707B (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 |