CN108226968B - 一种导航信号快速捕获方法 - Google Patents

一种导航信号快速捕获方法 Download PDF

Info

Publication number
CN108226968B
CN108226968B CN201810005376.XA CN201810005376A CN108226968B CN 108226968 B CN108226968 B CN 108226968B CN 201810005376 A CN201810005376 A CN 201810005376A CN 108226968 B CN108226968 B CN 108226968B
Authority
CN
China
Prior art keywords
code
sampling
signal
less
satellite
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
CN201810005376.XA
Other languages
English (en)
Other versions
CN108226968A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810005376.XA priority Critical patent/CN108226968B/zh
Publication of CN108226968A publication Critical patent/CN108226968A/zh
Application granted granted Critical
Publication of CN108226968B publication Critical patent/CN108226968B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, 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

本发明公开了一种导航信号快速捕获方法，其核心是减少二维搜索次数，首先通过组合多颗卫星的扩频码，获得各个卫星大致的码相位和多普勒信息，然后再利用一维搜索进一步确定可见卫星号。本发明实现了混合的一次多颗卫星搜索的方法，简化了传统的对卫星信号逐个进行二维搜索的卫星信号捕获方法，能够进一步减少卫星捕获方法的运算复杂度以及捕获时间。

Description

一种导航信号快速捕获方法
技术领域
本发明属于卫星通信技术领域，具体涉及一种导航信号快速捕获方法的设计。
背景技术
在卫星导航信号处理过程中，接收机在开始信号跟踪之前，首先需要估算出接收信号的载波频率(或者说载波多普勒频移)和码相位这两个参数，然后根据这些信号参数估计值初始化跟踪环路，以帮助接收通道展开对信号的跟踪，而信号捕获的目的正是为了获取所有可见卫星信号的载波频率和码相位的粗略估计值。
为了估算出某个卫星信号的载波频率和码相位这两个参数值，接收机的信号捕获过程一般是通过对该卫星信号的载波频率和码相位这二维进行扫描式搜索来完成的。一旦信号被搜索并得以确认，那么对该信号的捕获过程也就结束。
在对一个卫星信号进行频率和码相位的二维搜索之前，接收机首先要判定这颗卫星有可能可见而值得对其信号进行搜索。这样，卫星信号的捕获实际上是一个关于伪码(卫星号)、频率和码相位的三维搜索过程。
信号捕获的基础是测距码的相关性，所以对于搜索来说，只有当接收机内部复制的载波和测距码信号与接收信号一致时，相关器的输出才会达到最大值。如图1所示，每一个频带对应着一个载波频率搜索值，每一个码相位带对应着一个码相位搜索值，每一个频率带和一个码相位带的交点称为一个基本搜索单元。
卫星信号捕获方法主要有两种，第一种是串行搜索捕获法，串行搜索捕获方法是卫星导航***为解决捕获问题提出的最早的传统捕获方法。该方法是对码相位和多普勒频移分别进行串行搜索。由于码相位搜索的步长通常为半个码片，而一个码元周期较长，所以串行搜索效率低。多普勒频率搜索步长是根据相干积分时间而定，通常为相干积分时间的倒数，其搜索范围由接收机所在的动态状态决定。
第二种是并行搜索捕获法，并行搜索捕获方法是为了解决串行搜索时间长的问题而提出的。通常实现的是码相位采取并行搜索，而多普勒频移采取串行搜索。该方法利用快速傅里叶变换(Fast Fourier Transform，FFT)与快速傅立叶反变换(IFFT)将相关运算从时域转换到频域进行，降低了捕获所用时间，也为实时处理的实现奠定了基础。但是上述两种算法在卫星信号的三维搜索过程中仍然需要较大的计算量。关键的原因是对每颗卫星的捕获都需要进行二维搜索，有多少颗卫星就要进行多少次的二维搜索。二维搜索运算量大，复杂度高，给捕获过程带来负担。
发明内容
本发明的目的是为了解决现有的卫星信号捕获方法需要对卫星信号逐个进行二维搜索，从而造成计算量过大的问题，提出了一种导航信号快速捕获方法。
本发明的技术方案为：一种导航信号快速捕获方法，包括以下步骤：
S1、生成导航***中每颗卫星的C/A码CAk，并以2倍码率对C/A码CAk采样，得到采样C/A码CAsamk；其中k∈{1,2,...,K}，K为导航***中的卫星数量。
S2、将导航***中的K颗卫星分成J组，组间卫星不重复，将每组内的采样C/A码相加得到采样C/A码集合CAsamsetj；其中j∈{1,2,...,J}。
S3、对采样C/A码集合CAsamsetj做FFT运算并对FFT结果取共轭，得到FFT共轭结果
Figure BDA0001538455190000021
S4、接收中频采样信号s。
S5、对中频采样信号s进行下变频得到下变频信号xi；其中i为计数值，初始值为1。
下变频信号xi的计算公式为：
xi＝s·exp[-j2π(fIF+fD-(i-1)Δf)n/fs] (1)
其中fIF为在多普勒频偏为零的情况下中频采样信号的载波频率，fD为对中频采样信号预估的最大多普勒频偏，Δf为对中频采样信号预估的多普勒频移步进，fs为在多普勒频偏为零的情况下中频采样信号的采样率。
S6、以C/A码的2倍码率对下变频信号xi进行打包，即将下变频信号xi中每
Figure BDA0001538455190000022
个采样点累加为1个点，得到序列xpacki，其中fCA为C/A码码率。
S7、对序列xpacki做FFT运算，得到FFT序列XFi。
S8、令计数值i加1。
S9、若i≤I，则返回步骤S5，否则进入步骤S10；其中I为中频采样信号的多普勒频偏估计值个数，计算公式为：
Figure BDA0001538455190000023
S10、设置计数值i为1，设置j为1。
S11、根据步骤S3所得的FFT共轭结果
Figure BDA0001538455190000031
与步骤S7所得的FFT序列XFi计算得到相关运算结果corri j，计算公式为：
Figure BDA0001538455190000032
其中abs[·]为取模运算，ifft(·)为快速傅立叶逆变换。
S12、令计数值i加1。
S13、若i≤I，则返回步骤S11，否则进入步骤S14。
S14、设置计数值i为1，令j加1。
S15、若j≤J，则返回步骤S11，否则进入步骤S16。
S16、设置j为1。
S17、将上标相同的相关运算结果corri j作为一个相关运算结果集合corrj，统计相关运算结果集合corrj中峰均比大于预定门限的所有峰值点rt，并记录下所有的峰值点rt对应的码相位和多普勒频偏组合(pt,dplt)；其中t∈{1,2,...,Tj}，Tj为峰值点总数。
S18、若Tj<1，则进入步骤S19，否则进入步骤S21。
S19、令j加1。
S20、若j≤J，则返回步骤S17，否则捕获阶段结束。
S21、判断采样C/A码集合CAsamsetj是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，返回步骤S19，否则进入步骤S22。
S22、设置t为1。
S23、设有L个集合的存储空间，将采样C/A码集合CAsamsetj中的卫星分成L组，组内的C/A码按步骤S2的方法进行相加，存入上述L个集合的存储空间(若采样C/A码集合CAsamsetj组内卫星数量小于L，则将空出的存储空间复位为零)，得到采样C/A码分组集合
Figure BDA0001538455190000033
其中l∈{1,2,...,L}。
S24、设置l为1。
S25、按(pt,dplt)中的码相位对采样C/A码分组集合
Figure BDA0001538455190000034
进行循环移位，使之与峰值点rt对应的码相位对齐，得到C/A码循环移位集合
Figure BDA0001538455190000035
S26、将C/A码循环移位集合
Figure BDA0001538455190000041
与
Figure BDA0001538455190000042
进行相关运算，得到相关运算结果corrval；其中it为峰值点rt的多普勒频移对应计数值。
S27、比较关系数corrval与峰值点rt的模值大小，若|corrval|<|rt|，则进入步骤S28，否则进入步骤S34。
S28、令l加1。
S29、若l≤L，则返回步骤S25，否则输出错误信息，进入步骤S30。
S30、令t加1。
S31、若t≤Tj，则返回步骤S25，否则输出错误信息。
S32、令j加1。
S33、若j≤J，则返回步骤S17，否则捕获阶段结束。
S34、判断采样C/A码分组集合
Figure BDA0001538455190000043
是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，进入步骤S36，否则进入步骤S35。
S35、将采样C/A码分组集合
Figure BDA0001538455190000044
的卫星继续分成L组并把组内的C/A码按步骤S2的方法进行相加，存入步骤S23所述的L个集合的存储空间(若采样C/A码分组集合
Figure BDA0001538455190000045
组内卫星数量小于L，则将空出的存储空间复位为零)，返回步骤S24。
S36、令t加1。
S37、若t≤Tj，则返回步骤S25，否则进入步骤S38。
S38、令j加1。
S39、若j≤J，则返回步骤S17，否则捕获阶段结束。
本发明的有益效果是：本发明的核心是减少二维搜索次数，首先通过组合多颗卫星的扩频码，获得各个卫星大致的码相位和多普勒信息(此时码相位和多普勒信息还未与卫星号一一对应)，然后再利用一维搜索进一步确定可见卫星号。本发明实现了混合的一次多颗卫星搜索的方法，简化了传统的对卫星信号逐个进行二维搜索的卫星信号捕获方法，能够进一步减少卫星捕获方法的运算复杂度以及捕获时间。
附图说明
图1所示为卫星信号的三维搜索示意图。
图2所示为本发明实施例提供的一种导航信号快速捕获方法流程图。
具体实施方式
现在将参考附图来详细描述本发明的示例性实施方式。应当理解，附图中示出和描述的实施方式仅仅是示例性的，意在阐释本发明的原理和精神，而并非限制本发明的范围。
本发明实施例提供了一种导航信号快速捕获方法，如图2所示，包括以下步骤：
S1、生成导航***中每颗卫星的C/A码CAk，并以2倍码率对C/A码CAk采样，得到采样C/A码CAsamk；其中k∈{1,2,...,K}，K为导航***中的卫星数量。本发明实施例中，设置K＝32。
S2、将导航***中的K颗卫星分成J组，组间卫星不重复，将每组内的采样C/A码相加得到采样C/A码集合CAsamsetj；其中j∈{1,2,...,J}。本发明实施例中，设置J＝8。
S3、对采样C/A码集合CAsamsetj做FFT运算并对FFT结果取共轭，得到FFT共轭结果
Figure BDA0001538455190000051
S4、接收中频采样信号s。
S5、对中频采样信号s进行下变频得到下变频信号xi；其中i为计数值，初始值为1。
下变频信号xi的计算公式为：
xi＝s·exp[-j2π(fIF+fD-(i-1)Δf)n/fs] (1)
其中fIF为在多普勒频偏为零的情况下中频采样信号的载波频率，fD为对中频采样信号预估的最大多普勒频偏，Δf为对中频采样信号预估的多普勒频移步进，fs为在多普勒频偏为零的情况下中频采样信号的采样率。
本发明实施例中，设置fD＝±4kHz，Δf＝500Hz。
S6、以C/A码的2倍码率对下变频信号xi进行打包，即将下变频信号xi中每
Figure BDA0001538455190000052
个采样点累加为1个点，得到序列xpacki，其中fCA为C/A码码率。
本发明实施例中，设置fCA＝1.023MHz。
S7、对序列xpacki做FFT运算，得到FFT序列XFi。
S8、令计数值i加1。
S9、若i≤I，则返回步骤S5，否则进入步骤S10；其中I为中频采样信号的多普勒频偏估计值个数，计算公式为：
Figure BDA0001538455190000061
本发明实施例中，
Figure BDA0001538455190000062
S10、设置计数值i为1，设置j为1。
S11、根据步骤S3所得的FFT共轭结果
Figure BDA0001538455190000063
与步骤S7所得的FFT序列XFi计算得到相关运算结果corri j，计算公式为：
Figure BDA0001538455190000064
其中abs[·]为取模运算，ifft(·)为快速傅立叶逆变换。
S12、令计数值i加1。
S13、若i≤I，则返回步骤S11，否则进入步骤S14。
S14、设置计数值i为1，令j加1。
S15、若j≤J，则返回步骤S11，否则进入步骤S16。
S16、设置j为1。
S17、将上标相同的相关运算结果corri j作为一个相关运算结果集合corrj，统计相关运算结果集合corrj中峰均比大于预定门限的所有峰值点rt，并记录下所有的峰值点rt对应的码相位和多普勒频偏组合(pt,dplt)；其中t∈{1,2,...,Tj}，Tj为峰值点总数。
S18、若Tj<1，则进入步骤S19，否则进入步骤S21。
S19、令j加1。
S20、若j≤J，则返回步骤S17，否则捕获阶段结束。
S21、判断采样C/A码集合CAsamsetj是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，返回步骤S19，否则进入步骤S22。
本发明实施例中，采样C/A码集合CAsamsetj由4颗卫星的C/A码组成，则直接进入步骤S22。
S22、设置t为1。
S23、设有L个集合的存储空间，将采样C/A码集合CAsamsetj中的卫星分成L组，组内的C/A码按步骤S2的方法进行相加，存入上述L个集合的存储空间(若采样C/A码集合CAsamsetj组内卫星数量小于L，则将空出的存储空间复位为零)，得到采样C/A码分组集合
Figure BDA0001538455190000071
其中l∈{1,2,...,L}。
本发明实施例中，设置L＝2.
S24、设置l为1。
S25、按(pt,dplt)中的码相位对采样C/A码分组集合
Figure BDA0001538455190000072
进行循环移位，使之与峰值点rt对应的码相位对齐，得到C/A码循环移位集合
Figure BDA0001538455190000073
S26、将C/A码循环移位集合
Figure BDA0001538455190000074
与
Figure BDA0001538455190000075
进行相关运算，得到相关运算结果corrval；其中it为峰值点rt的多普勒频移对应计数值。
S27、比较关系数corrval与峰值点rt的模值大小，若|corrval|<|rt|，则进入步骤S28，否则进入步骤S34。
S28、令l加1。
S29、若l≤L，则返回步骤S25，否则输出错误信息，进入步骤S30。
S30、令t加1。
S31、若t≤Tj，则返回步骤S25，否则输出错误信息。
S32、令j加1。
S33、若j≤J，则返回步骤S17，否则捕获阶段结束。
S34、判断采样C/A码分组集合
Figure BDA0001538455190000076
是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，进入步骤S36，否则进入步骤S35。
S35、将采样C/A码分组集合
Figure BDA0001538455190000077
的卫星继续分成L＝2组并把组内的C/A码按步骤S2的方法进行相加，存入步骤S23所述的L个集合的存储空间(若采样C/A码分组集合
Figure BDA0001538455190000078
组内卫星数量小于L，则将空出的存储空间复位为零)，返回步骤S24。
S36、令t加1。
S37、若t≤Tj，则返回步骤S25，否则进入步骤S38。
S38、令j加1。
S39、若j≤J，则返回步骤S17，否则捕获阶段结束。
下面将本发明实施例提供的导航信号快速捕获方法与现有的两种卫星信号捕获方法进行算法复杂度对比。
假设某一卫星导航***共有K颗卫星，捕获阶段有I个多普勒频偏估计，有P个码相位估计，比较捕获一颗卫星需要的最大运算量如下：
传统的串行码相位搜索捕获方法的运算量量级是KIP2。
传统的并行码相位搜索捕获方法的运算量量级是KIPlog2P。
采用本发明实施例提供的导航信号快速捕获方法，在二维搜索时分成J组，在卫星PRN确认时分成L组，运算量量级是IJPlog2P。由于J<K(例如GPS中，K＝32,J＝4)，所以本发明的算法复杂度是比传统的两种算法复杂度要低的。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (5)

1.一种导航信号快速捕获方法，其特征在于，包括以下步骤：
S1、生成导航***中每颗卫星的C/A码CAk，并以2倍码率对C/A码CAk采样，得到采样C/A码CAsamk；其中k∈{1,2,...,K}，K为导航***中的卫星数量；
S2、将导航***中的K颗卫星分成J组，组间卫星不重复，将每组内的采样C/A码相加得到采样C/A码集合CAsamsetj；其中j∈{1,2,...,J}；
S3、对采样C/A码集合CAsamsetj做FFT运算并对FFT结果取共轭，得到FFT共轭结果
Figure FDA0002971312500000011
S4、接收中频采样信号s；
S5、对中频采样信号s进行下变频得到下变频信号xi；其中i为计数值，初始值为1；
S6、以C/A码的2倍码率对下变频信号xi进行打包，得到序列xpacki；
S7、对序列xpacki做FFT运算，得到FFT序列XFi；
S8、令计数值i加1；
S9、若i≤I，则返回步骤S5，否则进入步骤S10；其中I为中频采样信号的多普勒频偏估计值个数；
S10、设置计数值i为1，设置j为1；
S11、根据步骤S3所得的FFT共轭结果
Figure FDA0002971312500000012
与步骤S7所得的FFT序列XFi计算得到相关运算结果
Figure FDA0002971312500000013
S12、令计数值i加1；
S13、若i≤I，则返回步骤S11，否则进入步骤S14；
S14、设置计数值i为1，令j加1；
S15、若j≤J，则返回步骤S11，否则进入步骤S16；
S16、设置j为1；
S17、将上标相同的相关运算结果corri j作为一个相关运算结果集合corrj，统计相关运算结果集合corrj中峰均比大于预定门限的所有峰值点rt，并记录下所有的峰值点rt对应的码相位和多普勒频偏组合(pt,dplt)；其中t∈{1,2,...,Tj}，Tj为峰值点总数；
S18、若Tj<1，则进入步骤S19，否则进入步骤S21；
S19、令j加1；
S20、若j≤J，则返回步骤S17，否则捕获阶段结束；
S21、判断采样C/A码集合CAsamsetj是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，返回步骤S19，否则进入步骤S22；
S22、设置t为1；
S23、设有L个集合的存储空间，将采样C/A码集合CAsamsetj中的卫星分成L组，组内的C/A码按步骤S2的方法进行相加，存入上述L个集合的存储空间，得到采样C/A码分组集合
Figure FDA0002971312500000021
其中l∈{1,2,...,L}；若采样C/A码集合CAsamsetj组内卫星数量小于L，则将空出的存储空间复位为零；
S24、设置l为1；
S25、按(pt,dplt)中的码相位对采样C/A码分组集合
Figure FDA0002971312500000022
进行循环移位，使之与峰值点rt对应的码相位对齐，得到C/A码循环移位集合
Figure FDA0002971312500000023
S26、将C/A码循环移位集合
Figure FDA0002971312500000024
与
Figure FDA0002971312500000025
进行相关运算，得到相关运算结果corrval；其中it为峰值点rt的多普勒频移对应计数值；
S27、比较关系数corrval与峰值点rt的模值大小，若|corrval|<|rt|，则进入步骤S28，否则进入步骤S34；
S28、令l加1；
S29、若l≤L，则返回步骤S25，否则输出错误信息，进入步骤S30；
S30、令t加1；
S31、若t≤Tj，则返回步骤S25，否则输出错误信息；
S32、令j加1；
S33、若j≤J，则返回步骤S17，否则捕获阶段结束；
S34、判断采样C/A码分组集合
Figure FDA0002971312500000026
是否只由单颗卫星的C/A码组成，若是则说明捕获到了对应的单颗导航卫星的信号多普勒频偏和码相位，进入步骤S36，否则进入步骤S35；
S35、将采样C/A码分组集合
Figure FDA0002971312500000031
的卫星继续分成L组并把组内的C/A码按步骤S2的方法进行相加，存入步骤S23所述的L个集合的存储空间，返回步骤S24；若采样C/A码分组集合
Figure FDA0002971312500000032
组内卫星数量小于L，则将空出的存储空间复位为零；
S36、令t加1；
S37、若t≤Tj，则返回步骤S25，否则进入步骤S38；
S38、令j加1；
S39、若j≤J，则返回步骤S17，否则捕获阶段结束。
2.根据权利要求1所述的导航信号快速捕获方法，其特征在于，所述步骤S5中下变频信号xi的计算公式为：
xi＝s·exp[-j2π(fIF+fD-(i-1)Δf)n/fs] (1)
其中fIF为在多普勒频偏为零的情况下中频采样信号的载波频率，fD为对中频采样信号预估的最大多普勒频偏，Δf为对中频采样信号预估的多普勒频移步进，fs为在多普勒频偏为零的情况下中频采样信号的采样率。
3.根据权利要求1所述的导航信号快速捕获方法，其特征在于，所述步骤S6中对下变频信号xi进行打包的具体方法为：将下变频信号xi中每
Figure FDA0002971312500000033
个采样点累加为1个点，得到序列xpacki，其中fs为在多普勒频偏为零的情况下中频采样信号的采样率，fCA为C/A码码率。
4.根据权利要求1所述的导航信号快速捕获方法，其特征在于，所述步骤S9中多普勒频偏估计值个数I的计算公式为：
Figure FDA0002971312500000034
其中fD为对中频采样信号预估的最大多普勒频偏，Δf为对中频采样信号预估的多普勒频移步进。
5.根据权利要求1所述的导航信号快速捕获方法，其特征在于，所述步骤S11中相关运算结果corri j的计算公式为：
Figure FDA0002971312500000035
其中abs[·]为取模运算，ifft(·)为快速傅立叶逆变换。
CN201810005376.XA 2018-01-03 2018-01-03 一种导航信号快速捕获方法 Active CN108226968B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810005376.XA CN108226968B (zh) 2018-01-03 2018-01-03 一种导航信号快速捕获方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810005376.XA CN108226968B (zh) 2018-01-03 2018-01-03 一种导航信号快速捕获方法

Publications (2)

Publication Number Publication Date
CN108226968A CN108226968A (zh) 2018-06-29
CN108226968B true CN108226968B (zh) 2021-05-11

Family

ID=62642716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810005376.XA Active CN108226968B (zh) 2018-01-03 2018-01-03 一种导航信号快速捕获方法

Country Status (1)

Country Link
CN (1) CN108226968B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110068841A (zh) * 2019-05-06 2019-07-30 西安开阳微电子有限公司 一种卫星信号高灵敏度快速捕获方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101464508A (zh) * 2008-12-19 2009-06-24 苏州莱迪斯特电子有限公司 一种gps的c/a码信号的捕获方法
CN101515820A (zh) * 2008-02-22 2009-08-26 上海思谋微电子技术有限公司 定位接收机单通道多卫星并行搜索方法
CN101819278A (zh) * 2010-03-25 2010-09-01 北京航空航天大学 一种gps接收机的信号捕获方法
JP2012256991A (ja) * 2011-06-08 2012-12-27 Nec Aerospace Syst Ltd シンボル同期捕捉システム及びその方法
CN103078660A (zh) * 2013-01-06 2013-05-01 中国电子科技集团公司第十研究所 减少大动态下扩频码捕获时间的方法
CN103197328A (zh) * 2013-03-14 2013-07-10 中国人民解放军国防科学技术大学 一种基于2ms数据存储的高动态卫星接收机快速捕获方法
JP5483750B2 (ja) * 2009-11-30 2014-05-07 古野電気株式会社 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末
CN104101883A (zh) * 2013-04-11 2014-10-15 中国科学院微电子研究所 一种全球卫星导航信号的快速捕获方法及装置
CN105717522A (zh) * 2016-02-23 2016-06-29 电子科技大学 “北斗二代”b1频段弱信号捕获方法
CN106814375A (zh) * 2017-01-24 2017-06-09 中国电子科技集团公司第五十四研究所 一种抗转发式欺骗的捕获方法和接收机
CN107153206A (zh) * 2016-10-13 2017-09-12 南京航空航天大学 一种基于fft的北斗卫星信号快速捕获方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101515820A (zh) * 2008-02-22 2009-08-26 上海思谋微电子技术有限公司 定位接收机单通道多卫星并行搜索方法
CN101464508A (zh) * 2008-12-19 2009-06-24 苏州莱迪斯特电子有限公司 一种gps的c/a码信号的捕获方法
JP5483750B2 (ja) * 2009-11-30 2014-05-07 古野電気株式会社 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末
CN101819278A (zh) * 2010-03-25 2010-09-01 北京航空航天大学 一种gps接收机的信号捕获方法
JP2012256991A (ja) * 2011-06-08 2012-12-27 Nec Aerospace Syst Ltd シンボル同期捕捉システム及びその方法
CN103078660A (zh) * 2013-01-06 2013-05-01 中国电子科技集团公司第十研究所 减少大动态下扩频码捕获时间的方法
CN103197328A (zh) * 2013-03-14 2013-07-10 中国人民解放军国防科学技术大学 一种基于2ms数据存储的高动态卫星接收机快速捕获方法
CN104101883A (zh) * 2013-04-11 2014-10-15 中国科学院微电子研究所 一种全球卫星导航信号的快速捕获方法及装置
CN105717522A (zh) * 2016-02-23 2016-06-29 电子科技大学 “北斗二代”b1频段弱信号捕获方法
CN107153206A (zh) * 2016-10-13 2017-09-12 南京航空航天大学 一种基于fft的北斗卫星信号快速捕获方法
CN106814375A (zh) * 2017-01-24 2017-06-09 中国电子科技集团公司第五十四研究所 一种抗转发式欺骗的捕获方法和接收机

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A fast algorithm for long code direct acquisition;Haoxiang Liang等;《 2013 International Workshop on Microwave and Millimeter Wave Circuits and System Technology》;20140515;全文 *
A Fast Time Acquisition Method for PSK-2DSS Systems;LIN Jing-ran等;《2011 International Conference on Computational Problem-Solving (ICCP) 》;20111201;全文 *
GPS信号C/A码快速精密捕获技术研究;李新山等;《计算机应用研究》;20140430;第31卷(第4期);全文 *
Two dimension direct GPS signal acquisition based on unsymmetrical FFT;Hu Lu;《 2009 4th IEEE Conference on Industrial Electronics and Applications》;20090630;全文 *
一种基于FFT的高动态GPS信号快速捕获方法;刘晓莉等;《***仿真学报》;20070531;第19卷(第10期);全文 *

Also Published As

Publication number Publication date
CN108226968A (zh) 2018-06-29

Similar Documents

Publication Publication Date Title
CN106646546B (zh) 一种卫星信号多维快速捕获方法及***
CN101132191B (zh) 一种gnss接收机基带信号处理方法
RU2510134C2 (ru) Способ и система для обнаружения сигнала с расширенным спектром
US7526015B2 (en) Parallel correlator implementation using hybrid correlation in spread-spectrum communication
CN102162852A (zh) 大规模频率偏移下微弱gnss信号捕获方法及其***
TW201512689A (zh) 用於處理組合式導航信號的方法、裝置及製品
CN100501441C (zh) 定位全球定位***相关峰值信号的方法和***
CN113009522B (zh) 多普勒频率残差校正的长时间相干积分捕获算法模块
CN109655847B (zh) 一种适于动态信号的快速捕获方法
CN102798871B (zh) 基于伪码重构的伪码捕获方法及装置
CN104215979A (zh) 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法
CN107247276B (zh) 基于重叠多块补零算法的弱信号捕获方法
CN113820731B (zh) 一种卫星导航接收机无周期长码信号的捕获方法及***
CN106526631B (zh) 北斗b1频点卫星信号载波频率高精度估计方法
CN103760578B (zh) 一种gnss卫星导航信号的无模糊跟踪方法
CN101072045A (zh) 基于频率折叠的短周期扩频码时频并行搜索方法
CN109581434B (zh) 一种B2a信号捕获方法及装置
CN108226968B (zh) 一种导航信号快速捕获方法
CN110018500B (zh) 一种基于圆周移位的北斗卫星信号捕获方法
US8149937B2 (en) Signal receiver and method for estimating residual doppler frequencies of signals thereof
KR102155083B1 (ko) 이진 오프셋 반송파 신호를 위한 비 모호 상관 함수 생성 방법 및 장치
CN109921824B (zh) 一种高阶boc调制信号无模糊捕获方法
CN110958035A (zh) 一种基于三角包络拟合的频率二次捕获装置及方法
CN115113243A (zh) 一种联合频率插值估计的导航信号捕获方法
CN104993844B (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