CN112649819A - 一种高动态扩频信号捕获装置及捕获方法 - Google Patents

一种高动态扩频信号捕获装置及捕获方法 Download PDF

Info

Publication number
CN112649819A
CN112649819A CN202011328588.5A CN202011328588A CN112649819A CN 112649819 A CN112649819 A CN 112649819A CN 202011328588 A CN202011328588 A CN 202011328588A CN 112649819 A CN112649819 A CN 112649819A
Authority
CN
China
Prior art keywords
data
pseudo code
input
acquisition
module
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
Application number
CN202011328588.5A
Other languages
English (en)
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.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
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 National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN202011328588.5A priority Critical patent/CN112649819A/zh
Publication of CN112649819A publication Critical patent/CN112649819A/zh
Pending legal-status Critical Current

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
    • 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/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明属于扩频信号接收处理的技术领域，具体地说，涉及一种高动态扩频信号捕获装置，包括：FPGA芯片和设置在FPGA芯片上的频率搜索控制器、PRN码发生器、本地振荡器产生模块、第一FFT模块，第二FFT模块、第一复数乘法器、第二复数乘法器、第三复数乘法器、IFFT模块、峰值检测模块、数据处理模块和计数器；所述数据处理模块，用于根据判决结果，获取当前捕获的码相位；并结合下一次捕获得到的码相位，得到捕获时的码相位，完成对高动态扩频信号的捕获。

Description

一种高动态扩频信号捕获装置及捕获方法
技术领域
本发明属于空间飞行器上设备高动态扩频信号捕获的技术领域，具体地说，涉及一种高动态扩频信号捕获装置及捕获方法。
背景技术
由于飞行载体的高速运动，空间飞行器上导航接收机所接收到的导航信号通常会有较大的多普勒频移，而较大的多普勒偏移会给高动态扩频信号捕获带来以下问题：
问题一：大的多普勒偏移增加了扩频信号捕获中多普勒搜索范围，增长了捕获时间。较长的捕获时间会使捕获到的码相位发生码相位漂移，捕获到伪码相位时刻与伪码相位输出时刻并非同一时刻。
其中，码相位漂移是指，由于在捕获过程中，捕获到的时刻和结果输出的时刻是有时间差的，而大的多普勒在搜索时需要消耗更多的时间，而在捕获过程中，捕获码相位是一直在向后搜索，也就是说，码相位是一直在变化的，我们在动态不高的情况下，将结果输出时刻近似于捕获时刻。而在高动态情况下，捕获时间长了，不能近似等同，因此，捕获到的时刻和结果输出时刻的码相位不一样。
问题二：多普勒偏移同时会造成接收信号的伪码速率变化。在导航信号的捕获过程中，使得本地伪码与输入伪码之间速率差距较大，影响捕获结果。
在捕获的过程中，由于大多普勒的影响，输入伪码带有较大的多普勒，导致了输入伪码带有的码元被拉长或者压缩，码元长度发生了变化，而本地伪码的码元宽度是不变的。最佳的捕获方案是输入伪码和本地伪码的速率一致，码元宽度相等。因此，如果输入伪码由于大的多普勒影响导致输入伪码速率和本地伪码速率有很大差距时，难以获得精确的捕获结果。
现有的解决高动态扩频信号捕获的方法主要有PMF-FFT(Partial MatchedFilters-Fast Fourier Transform)基于部分匹配滤波器的快速傅里叶捕获方法)方法，但是，该方法存在消耗资源比较多，捕获时间长的问题。
发明内容
为解决现有技术存在的上述缺陷，本发明提出了一种高动态扩频信号捕获装置，解决高动态扩频信号捕获中伪码相位漂移问题以及本地伪码与输入伪码速率不一致问题；该装置包括：FPGA芯片和设置在FPGA芯片上的频率搜索控制器、PRN码发生器、本地振荡器产生模块、第一FFT模块，第二FFT模块、第一复数乘法器、第二复数乘法器、第三复数乘法器、IFFT模块、峰值检测模块、数据处理模块和计数器；
所述频率搜索控制器，用于实时产生搜索频率值并输出至PRN码发生器；
所述PRN码发生器，用于根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码并输出至第一FFT模块；其中，更改频率值后的本地伪码为伪码序列，其中包括多个伪码；
所述第一FFT模块，用于对更改后频率值后的本地伪码进行FFT变换，将更改后频率值后的本地伪码变换到频域，得到频域数据，并取频域数据的复共轭，输出第一输入数据，并输出至第三复数乘法器；
所述本地振荡器产生模块，用于产生sin数据和cos数据,将sin数据输出至第一复数乘法器；将cos数据输出至第二复数乘法器；
所述第一复数乘法器，用于将第一输入伪码与sin数据进行复乘运算，得到第一输入伪码数据并输出至第二FFT模块；
所述第二复数乘法器，用于将第一输入伪码与cos数据进行复乘运算，得到第二输入伪码数据并输出至第二FFT模块；
所述第二FFT模块，用于将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第二输入数据频域并输出至第三复数乘法器；
所述第三复数乘法器，用于将第一输入数据和第二输入数据频域进行复乘运算，得到第三输入数据并输出至IFFT模块；
所述IFFT模块，用于对第三输入数据进行IFFT变换，得到时域数据，对时域数据中的每个数据进行求绝对值处理，得到多个时域数值，进而得到多个模方并输出至峰值检测模块；
所述峰值检测模块，用于从多个模方中选取其中的最大值，对该最大值进行峰值检测判决,将判决结果输出至数据处理模块；
所述数据处理模块，用于根据判决结果，获取当前捕获的码相位；并结合下一次捕获得到的码相位，得到捕获时的码相位，完成对高动态扩频信号的捕获；
所述计数器，用于统计捕获过程的次数；其中，捕获过程的次数大于或等于2。
作为上述技术方案的改进之一，所述本地振荡器产生模块为载波发生器。
作为上述技术方案的改进之一，所述峰值检测模块的具体判决过程为：
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为当前捕获的码相位P1，并将当前捕获的码相位P1发送至数据处理模块；
如果该最大值小于预先设定的门限阈值，则捕获失败。
作为上述技术方案的改进之一，所述数据处理模块的具体过程为：
根据当前捕获的码相位P1；并结合下一次捕获得到的码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩频信号的捕获。
本发明还提供了一种高动态扩频信号捕获方法，该方法包括：
通过第一次捕获，得到所捕获的第一捕获伪码相位；
通过第二次捕获，得到所捕获的第二捕获伪码相位；
计算两次捕获得出的第一捕获伪码相位和第二捕获伪码相位的码相位差；
根据该码相位差，得到捕获时的码相位，实现高动态扩频信号的捕获。
作为上述技术方案的改进之一，所述通过第一次捕获，得到所捕获的第一捕获伪码相位；其具体过程为：
所述PRN码发生器根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码，并输出至第一FFT模块；
第一FFT模块对更改频率值后的本地伪码进行FFT变换，将更改频率值后的本地伪码变换到频域，得到频域数据，并取频域数据的复共轭，输出第一输入数据；
第二FFT模块将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第一输入数据频域；
所述第三复数乘法器将第一输入数据和第一输入数据频域进行复乘运算，得到第二输入数据；
所述IFFT模块对第二输入数据进行IFFT变换，得到第一时域数据，对第一时域数据中的每个数据进行求绝对值处理，得到多个第一时域数值，进而得到多个模方；
所述峰值检测模块从得到的多个模方中选取其中的最大值，对该最大值进行检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第一次捕获码相位；
如果该最大值小于预先设定的门限阈值，则捕获失败。
作为上述技术方案的改进之一，第二输入数据频域的获取过程具体为：
本地振荡器产生模块产生sin数据和cos数据；
第一复数乘法器将输入伪码与sin数据进行复乘运算，得到第一输入伪码数据；
第二复数乘法器将输入伪码与cos数据进行复乘运算，得到第二输入伪码数据；
第二FFT模块将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第一输入数据频域。
作为上述技术方案的改进之一，所述通过第二次捕获，得到所捕获的第二捕获伪码相位；其具体过程为：
所述PRN码发生器根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码，并输出至第一FFT模块；
第一FFT模块对更改频率值后的本地伪码进行FFT变换，将更改频率值后的本地伪码变换到频域数据，取该频域数据的复共轭，输出第一输入数据；
第二FFT模块将第三输入伪码数据和第四输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第二输入数据频域；
所述第三复数乘法器将第一输入数据和第二输入数据频域进行复乘运算，得到第三输入数据；
所述IFFT模块对第三输入数据进行IFFT变换，得到第二时域数据，对第二时域数据中的每个数据进行求绝对值处理，得到多个第二时域数值，进而得到多个模方；
所述峰值检测模块从得到的多个模方中选取其中的最大值，对该最大值进行峰值检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第二次捕获码相位；
如果该最大值小于预先设定的门限阈值，则捕获失败。
作为上述技术方案的改进之一，所述根据该码相位差，得到捕获时的码相位，实现高动态扩频信号的捕获；其具体过程为：
所述数据处理模块根据判决结果，获取第一次捕获码相位P1；并结合第二次捕获得到的第二次捕获码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩信号的捕获。
本发明与现有技术相比的有益效果是：
1、本发明解决了高动态情况下扩频信号捕获的出现的大的多普勒偏移增加了扩频信号捕获中多普勒搜索范围，增长了捕获时间的问题，以及多普勒偏移同时会造成接收信号的伪码速率变化的问题，以及在导航信号的捕获过程中，使得本地伪码与输入伪码之间速率差距较大，影响捕获结果的问题，在几乎不增加资源的前提下，仅采用二次捕获法，实现了高动态信号的捕获；
2、可以快速的进行高动态扩频信号的捕获，捕获时间短。
附图说明
图1是本发明的一种高动态扩频信号捕获装置的结构示意图；
图2是更改本地伪码的频率值的流程图；
图3是图2的码相位值和多普勒值组成的二维搜索原理图。
具体实施方式
现结合附图对本发明作进一步的描述。
如图1所示，本发明提供了一种高动态扩频信号捕获装置，该装置包括：FPGA芯片和设置在FPGA芯片上的频率搜索控制器、PRN码发生器、本地振荡器产生模块、第一FFT模块，第二FFT模块、第一复数乘法器、第二复数乘法器、第三复数乘法器、IFFT模块、峰值检测模块、数据处理模块和计数器；
所述PRN(pseudo random noise)码发生器，根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码并输出至第一FFT模块；其中，更改频率值后的本地伪码为伪码序列，其中包括多个伪码；
所述频率搜索控制器，用于实时产生搜索频率值，并输出至PRN码发生器；其中，该搜索频率值是在二维搜索过程中所搜索到的频率范围，如图1所示，由频率搜索控制器来产生。
所述第一FFT模块，用于对PRN码发生器输出的更改频率值后的本地伪码x(n)进行FFT变换，将更改频率值后的本地伪码x(n)变换到频域，得到频域数据X(k)，取频域数据X(k)的复共轭，输出第一输入数据X(k)*，并输出至第三复数乘法器；
所述本地振荡器产生模块为载波发生器，用于产生sin数据和cos数据，将sin数据输出至第一复数乘法器；将cos数据输出至第二复数乘法器；
其中，sin数据为不同频率的正弦波；cos数据为不同频率的余弦波；产生sin数据和cos数据是用一个DDS(Direct Digital Synthesizer)发生器；其中，DDS发生器可以通过程序实现，也可以通过IP核实现。
所述第一复数乘法器，用于将第一输入伪码r(n)与载波发生器输出的sin数据进行复乘运算，得到第一输入伪码数据Q，并输出至第二FFT模块；
所述第二复数乘法器，用于将第一输入伪码r(n)与载波发生器输出的cos数据进行复乘运算，得到第二输入伪码数据I，并输出至第二FFT模块；
所述第二FFT模块，用于将第一输入伪码数据Q和第二输入伪码数据I进行整合，得到整合后的输入伪码数据S，并对其进行FFT操作，得到第一输入数据频域R(k)；
所述第三复数乘法器，用于将第一输入数据X(k)*和第一输入数据频域R(k)进行复乘运算，得到第二输入数据L(k)；
所述IFFT模块，用于对第二输入数据L(k)进行IFFT变换，得到第一时域数据l(n)，对第一时域数据为l(n)中的每个数据l(1)，l(2)，l(3)…l(n)进行求绝对值处理，得到多个第一时域数值|l(1)|，|l(2)|，|l(3)|，…|l(n)|，进而得到多个模方|l(1)|2，|l(2)|2，|l(3)|2，…|l(n)|2；
所述峰值检测模块，用于对得到的多个模方|l(n)|2，选取其中的最大值，对该最大值进行峰值检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为当前捕获的码相位P1；
如果该最大值小于预先设定的门限阈值，则捕获失败；
所述数据处理模块，用于根据判决结果，获取当前捕获的码相位P1；并结合下一次捕获得到的码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩频信号的捕获；
所述计数器，用于统计捕获过程的次数；其中，捕获过程的次数大于或等于2。
本发明提供了一种高动态扩频信号捕获方法，该方法包括：
通过第一次捕获，得到所捕获的第一捕获伪码相位P1；
具体地，所述PRN码发生器根据实时产生的搜索频率值,更改PRN码发生器产生的本地伪码的频率值，并输出更改频率值后的本地伪码x(n)；其中，更改频率值后的本地伪码x(n)为伪码序列，其中包括多个伪码；
第一FFT模块对PRN码发生器输出的更改频率值后的本地伪码x(n)进行FFT变换，将更改频率值后的本地伪码x(n)变换到频域，得到频域数据X(k)，取频域数据X(k)的复共轭，输出第一输入数据X(k)*；
第二FFT模块将第一输入伪码数据Q和第二输入伪码数据I进行整合，得到整合后的输入伪码数据S，并对其进行FFT操作，得到第一输入数据频域R(k)；
具体地，本地振荡器产生模块产生sin数据和cos数据；其中，sin数据为不同频率的正弦波；cos数据为不同频率的余弦波；
第一复数乘法器将第一输入伪码r(n)与sin数据进行复乘运算，得到第一输入伪码数据Q；
第二复数乘法器将第一输入伪码r(n)与cos数据进行复乘运算，得到第二输入伪码数据I；
第二FFT模块将第一输入伪码数据Q和第二输入伪码数据I进行整合，得到整合后的输入伪码数据S，并对其进行FFT操作，得到第一输入数据频域R(k)。
所述第三复数乘法器将第一输入数据X(k)*和第一输入数据频域R(k)进行复乘运算，得到第二输入数据L(k)；
所述IFFT模块对第二输入数据L(k)进行IFFT变换，得到第一时域数据l(n)，对第一时域数据为l(n)中的每个数据l(1)，l(2)，l(3)…l(n)进行求绝对值处理，得到多个第一时域数值|l(1)|，|l(2)|，|l(3)|，…|l(n)|，进而得到多个模方|l(1)|2，|l(2)|2，|l(3)|2，…|l(n)|2；
所述峰值检测模块对得到的多个模方|l(n)|2，选取其中的最大值，对该最大值进行峰值检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第一次捕获码相位P1；
如果该最大值小于预先设定的门限阈值，则捕获失败。
如果第一次捕获成功，则进行第二次捕获；通过第二次捕获，得到所捕获的第二捕获伪码相位P2；
具体地，所述PRN码发生器根据实时产生的搜索频率值,更改PRN码发生器产生的本地伪码的频率值，并输出更改频率值后的本地伪码x(n)；其中，更改频率值后的本地伪码x(n)为伪码序列，其中包括多个伪码；
第一FFT模块对PRN码发生器输出的更改频率值后的本地伪码x(n)进行FFT变换，将更改频率值后的本地伪码x(n)变换到频域，得到频域数据X(k)，取频域数据X(k)的复共轭，输出第一输入数据X(k)*；
第二FFT模块将第三输入伪码数据Q’和第四输入伪码数据I’进行整合，得到整合后的输入伪码数据S’，并对其进行FFT操作，得到第二输入数据频域R(k)1；
具体地，本地振荡器产生模块产生sin数据和cos数据；
第一复数乘法器将第二输入伪码r(n)1与sin数据进行复乘运算，得到第一输入伪码数据Q’；
第二复数乘法器将第二输入伪码r(n)1与cos数据进行复乘运算，得到第二输入伪码数据I’；
第二FFT模块将第一输入伪码数据Q’和第二输入伪码数据I’进行整合，得到整合后的输入伪码数据S’，并对其进行FFT操作，得到第二输入数据频域R(k)1。
所述第三复数乘法器，用于将第一输入数据X(k)*和第二输入数据频域R(k)1进行复乘运算，得到第三输入数据L(k)1；
所述IFFT模块对第三输入数据L(k)1进行IFFT变换，得到第二时域l(n)1，对第二时域数据为l(n)1中的每个数据l(1)1，l(2)1，l(3)1…l(n)1进行求绝对值处理，得到多个第二时域数值|l(1)|1，|l(2)|1，|l(3)|1…|l(n)|1，进而得到多个模方|l(1)|1 2，|l(2)|1 2，|l(3)|1 2…|l(n)|1 2；
所述峰值检测模块对得到的多个模方|l(n)1|2，选取其中的最大值，对该最大值进行检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第二次捕获码相位P2；
如果该最大值小于预先设定的门限阈值，则捕获失败。
计算两次捕获得出的第一捕获伪码相位和第二捕获伪码相位的码相位差P2-P1，该码相位差为由于捕获时间和多普勒频偏导致的码相位漂移；
根据该码相位差，得到捕获时的码相位Pinitial＝P1-(P2-P1)，实现高动态扩频信号的捕获。
具体地，所述数据处理模块根据判决结果，获取第一次捕获码相位P1；并结合第二次捕获得到的第二次捕获码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩信号的捕获。
本发明方法采用二次捕获方法，可以估计捕获码相位，解决了高动态条件下伪码相位漂移问题。用根据搜索频点更新本地伪码方法减小了本地伪码和输入伪码速率差，如图2所示，搜索控制器接收到来自多普勒搜索控制的当前频率搜索单元对应的搜索频率值，根据搜索频率值，实时更新伪码频率字，经相位累加器后，查找对应的码表，输出更改频率值后的本地伪码，使得捕获信号频点的本地伪码与输入伪码的伪码速率之差不大于1/2搜索步进。
如图2所示，根据实时产生的搜索频率点对应的搜索频率值更新PRN码发生器产生的本地伪码的频率值，缩小输入伪码的速率和本地伪码的速率之间的差距，获得精确的捕获时的码相位，从而提高捕获结果的精度。搜索控制器接收到来自多普勒控制的当前频率搜索单元对应的多普勒值，PRN码发生器接收到频率搜索控制器发送的更改频率值后的本地伪码，并输出更改频率值后的本地伪码；其中，更改后的本地伪码对应的频率值为当前频率搜索单元对应的多普勒值，根据搜索频点更新伪码频率字，经相位累加器后查找对应码表，输出当前频率搜索单元对应的更改频率值后的本地伪码。
如图3所示，是一个总体二维搜索过程的示意图，该搜索过程就是捕获过程；横坐标值为码相位搜索范围，纵坐标值为多普勒搜索范围，由于输入伪码的码相位和多普勒是未知的，但是可以肯定的是输入的码相位和多普勒在此范围内。扩频信号的捕获过程属于二维搜索，如图3所示，在多普勒这一维度搜索时，设定搜索范围-D到+D，频率步进为d，则搜索的顺序为0-+d、-d-0、+d-+2d、-2d--d、......。在每一个搜索区间，设定本地伪码的频率为此搜索区间的中间值。例如：搜索区间为+3d-+4d，则设定本地伪码频率为+3.5d。这样当搜索到输入伪码所在的多普勒区间时，本地伪码和输入伪码多普勒差距小于1/2d，所以减小了本地伪码和输入伪码的速率差。图3中所示的灰色方格表示捕获过程的初始位置；从这个初始位置开设，对每个多普勒值和对应的码相位开始逐一搜索，从而获取捕获时的码相位。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (9)

1.一种高动态扩频信号捕获装置，其特征在于，该装置包括：FPGA芯片和设置在FPGA芯片上的频率搜索控制器、PRN码发生器、本地振荡器产生模块、第一FFT模块，第二FFT模块、第一复数乘法器、第二复数乘法器、第三复数乘法器、IFFT模块、峰值检测模块、数据处理模块和计数器；
所述频率搜索控制器，用于实时产生搜索频率值并输出至PRN码发生器；
所述PRN码发生器，用于根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码并输出至第一FFT模块；其中，更改频率值后的本地伪码为伪码序列，其中包括多个伪码；
所述第一FFT模块，用于对更改后频率值后的本地伪码进行FFT变换，将更改后频率值后的本地伪码变换到频域，得到频域数据，并取频域数据的复共轭，输出第一输入数据，并输出至第三复数乘法器；
所述本地振荡器产生模块，用于产生sin数据和cos数据,将sin数据输出至第一复数乘法器；将cos数据输出至第二复数乘法器；
所述第一复数乘法器，用于将第一输入伪码与sin数据进行复乘运算，得到第一输入伪码数据并输出至第二FFT模块；
所述第二复数乘法器，用于将第一输入伪码与cos数据进行复乘运算，得到第二输入伪码数据并输出至第二FFT模块；
所述第二FFT模块，用于将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第二输入数据频域并输出至第三复数乘法器；
所述第三复数乘法器，用于将第一输入数据和第二输入数据频域进行复乘运算，得到第三输入数据并输出至IFFT模块；
所述IFFT模块，用于对第三输入数据进行IFFT变换，得到时域数据，对时域数据中的每个数据进行求绝对值处理，得到多个时域数值，进而得到多个模方并输出至峰值检测模块；
所述峰值检测模块，用于从多个模方中选取其中的最大值，对该最大值进行峰值检测判决,将判决结果输出至数据处理模块；
所述数据处理模块，用于根据判决结果，获取当前捕获的码相位；并结合下一次捕获得到的码相位，得到捕获时的码相位，完成对高动态扩频信号的捕获；
所述计数器，用于统计捕获过程的次数；其中，捕获过程的次数大于或等于2。
2.根据权利要求1所述的高动态扩频信号捕获装置，其特征在于，所述本地振荡器产生模块为载波发生器。
3.根据权利要求1所述的高动态扩频信号捕获装置，其特征在于，所述峰值检测模块的具体判决过程为：
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为当前捕获的码相位P1，并将当前捕获的码相位P1发送至数据处理模块；
如果该最大值小于预先设定的门限阈值，则捕获失败。
4.根据权利要求3所述的高动态扩频信号捕获装置，其特征在于，所述数据处理模块的具体过程为：
根据当前捕获的码相位P1；并结合下一次捕获得到的码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩频信号的捕获。
5.一种高动态扩频信号捕获方法，其特征在于，该方法基于上述权利要求1-5中任一所述的高动态扩频信号捕获装置实现，该方法包括：
通过第一次捕获，得到所捕获的第一捕获伪码相位；
通过第二次捕获，得到所捕获的第二捕获伪码相位；
计算两次捕获得出的第一捕获伪码相位和第二捕获伪码相位的码相位差；
根据该码相位差，得到捕获时的码相位，实现高动态扩频信号的捕获。
6.根据权利要求5所述的高动态扩频信号捕获方法，其特征在于，所述通过第一次捕获，得到所捕获的第一捕获伪码相位；其具体过程为：
所述PRN码发生器根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码，并输出至第一FFT模块；
第一FFT模块对更改频率值后的本地伪码进行FFT变换，将更改频率值后的本地伪码变换到频域，得到频域数据，并取频域数据的复共轭，输出第一输入数据；
第二FFT模块将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第一输入数据频域；
所述第三复数乘法器将第一输入数据和第一输入数据频域进行复乘运算，得到第二输入数据；
所述IFFT模块对第二输入数据进行IFFT变换，得到第一时域数据，对第一时域数据中的每个数据进行求绝对值处理，得到多个第一时域数值，进而得到多个模方；
所述峰值检测模块从得到的多个模方中选取其中的最大值，对该最大值进行峰值检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第一次捕获码相位；
如果该最大值小于预先设定的门限阈值，则捕获失败。
7.根据权利要求6所述的高动态扩频信号捕获方法，其特征在于，第二输入数据频域的获取过程具体为：
本地振荡器产生模块产生sin数据和cos数据；
第一复数乘法器将输入伪码与sin数据进行复乘运算，得到第一输入伪码数据；
第二复数乘法器将输入伪码与cos数据进行复乘运算，得到第二输入伪码数据；
第二FFT模块将第一输入伪码数据和第二输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第一输入数据频域。
8.根据权利要求5所述的高动态扩频信号捕获方法，其特征在于，所述通过第二次捕获，得到所捕获的第二捕获伪码相位；其具体过程为：
所述PRN码发生器根据搜索频率值更改本地伪码的频率值，并输出更改频率值后的本地伪码，并输出至第一FFT模块；
第一FFT模块对更改频率值后的本地伪码进行FFT变换，将更改频率值后的本地伪码变换到频域数据，取该频域数据的复共轭，输出第一输入数据；
第二FFT模块将第三输入伪码数据和第四输入伪码数据进行整合，得到整合后的输入伪码数据，并对其进行FFT操作，得到第二输入数据频域；
所述第三复数乘法器将第一输入数据和第二输入数据频域进行复乘运算，得到第三输入数据；
所述IFFT模块对第三输入数据进行IFFT变换，得到第二时域数据，对第二时域数据中的每个数据进行求绝对值处理，得到多个第二时域数值，进而得到多个模方；
所述峰值检测模块从得到的多个模方中选取其中的最大值，对该最大值进行峰值检测判决；
如果该最大值大于或等于预先设定的门限阈值，则捕获成功，并将该最大值作为第二次捕获码相位；
如果该最大值小于预先设定的门限阈值，则捕获失败。
9.根据权利要求5所述的高动态扩频信号捕获方法，其特征在于，所述根据该码相位差，得到捕获时的码相位，实现高动态扩频信号的捕获；其具体过程为：
所述数据处理模块根据判决结果，获取第一次捕获码相位P1；并结合第二次捕获得到的第二次捕获码相位P2，计算前后相邻两次捕获得出的码相位差P2-P1，作为捕获高动态扩频信号的捕获时间，得到捕获时的码相位Pinitial＝P1-(P2-P1)，作为捕获时的初始相位值，完成对高动态扩信号的捕获。
CN202011328588.5A 2020-11-24 2020-11-24 一种高动态扩频信号捕获装置及捕获方法 Pending CN112649819A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011328588.5A CN112649819A (zh) 2020-11-24 2020-11-24 一种高动态扩频信号捕获装置及捕获方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011328588.5A CN112649819A (zh) 2020-11-24 2020-11-24 一种高动态扩频信号捕获装置及捕获方法

Publications (1)

Publication Number Publication Date
CN112649819A true CN112649819A (zh) 2021-04-13

Family

ID=75349957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011328588.5A Pending CN112649819A (zh) 2020-11-24 2020-11-24 一种高动态扩频信号捕获装置及捕获方法

Country Status (1)

Country Link
CN (1) CN112649819A (zh)

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5271034A (en) * 1991-08-26 1993-12-14 Avion Systems, Inc. System and method for receiving and decoding global positioning satellite signals
US20020051434A1 (en) * 1997-10-23 2002-05-02 Ozluturk Fatih M. Method for using rapid acquisition spreading codes for spread-spectrum communications
US6407699B1 (en) * 2000-04-14 2002-06-18 Chun Yang Method and device for rapidly extracting time and frequency parameters from high dynamic direct sequence spread spectrum radio signals under interference
JP2003032144A (ja) * 2001-07-17 2003-01-31 Furuno Electric Co Ltd スペクトル拡散信号捕捉装置および方法
US20060013287A1 (en) * 2002-10-15 2006-01-19 Per-Ludvig Normark Spread spectrum signal processing
US20060115022A1 (en) * 2004-03-19 2006-06-01 Ziedan Nesreen I System and method for high dynamic acquisition and tracking of signals from the global positioning system
CN101625404A (zh) * 2008-07-09 2010-01-13 杭州中科微电子有限公司 Gps信号大规模并行快速捕获方法及其模块
CN102565822A (zh) * 2011-10-09 2012-07-11 中国科学院计算技术研究所 对gps l5信号的捕获方法及装置
CN102571137A (zh) * 2012-02-14 2012-07-11 浙江大学 一种全数字直接序列扩频通信系统及其伪码快速捕获方法
CN103529458A (zh) * 2013-10-18 2014-01-22 北京邮电大学 卫星捕获方法及装置
CN103595441A (zh) * 2013-10-30 2014-02-19 中国运载火箭技术研究院 适用于重复使用轨道飞行器的极高动态扩频信号快捕方法
US20140300514A1 (en) * 2012-12-05 2014-10-09 Samsung Electronics Co., Ltd. Method and apparatus for acquiring signal of global navigation satellite system
CN104360357A (zh) * 2014-11-24 2015-02-18 四川九洲电器集团有限责任公司 基于循环方式的北斗卫星信号快速捕获方法及系统
US20150059432A1 (en) * 2013-08-27 2015-03-05 Crystal Instruments Corporation Cross-path phase calibration for high dynamic range data acquisition
US20160116599A1 (en) * 2013-06-05 2016-04-28 Airbus Defence And Space Limited Receiver and method for direct sequence spread spectrum signals
CN105846855A (zh) * 2016-02-28 2016-08-10 中国人民解放军63686部队 一种基于频率引导的大动态扩频信号快速捕获方法
CN106646546A (zh) * 2016-12-20 2017-05-10 南京六九零二科技有限公司 一种卫星信号多维快速捕获方法及系统
CN107135013A (zh) * 2017-05-04 2017-09-05 中国电子科技集团公司第五十四研究所 一种直接序列扩频信号快速捕获方法
CN109100747A (zh) * 2018-07-20 2018-12-28 西安空间无线电技术研究所 一种适用于多路随机接入信号的快速捕获系统及方法
CN110824512A (zh) * 2019-11-26 2020-02-21 中国科学院国家空间科学中心 一种非均匀码片实时延迟多普勒映像数据生成器
CN111694027A (zh) * 2020-06-04 2020-09-22 长沙北斗产业安全技术研究院有限公司 超大动态扩频信号捕获方法与装置
CN111934710A (zh) * 2020-07-06 2020-11-13 南京天际砺剑科技有限公司 一种高动态扩频信号快速捕获算法

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5271034A (en) * 1991-08-26 1993-12-14 Avion Systems, Inc. System and method for receiving and decoding global positioning satellite signals
US20020051434A1 (en) * 1997-10-23 2002-05-02 Ozluturk Fatih M. Method for using rapid acquisition spreading codes for spread-spectrum communications
US6407699B1 (en) * 2000-04-14 2002-06-18 Chun Yang Method and device for rapidly extracting time and frequency parameters from high dynamic direct sequence spread spectrum radio signals under interference
JP2003032144A (ja) * 2001-07-17 2003-01-31 Furuno Electric Co Ltd スペクトル拡散信号捕捉装置および方法
US20060013287A1 (en) * 2002-10-15 2006-01-19 Per-Ludvig Normark Spread spectrum signal processing
US20060115022A1 (en) * 2004-03-19 2006-06-01 Ziedan Nesreen I System and method for high dynamic acquisition and tracking of signals from the global positioning system
CN101625404A (zh) * 2008-07-09 2010-01-13 杭州中科微电子有限公司 Gps信号大规模并行快速捕获方法及其模块
CN102565822A (zh) * 2011-10-09 2012-07-11 中国科学院计算技术研究所 对gps l5信号的捕获方法及装置
CN102571137A (zh) * 2012-02-14 2012-07-11 浙江大学 一种全数字直接序列扩频通信系统及其伪码快速捕获方法
US20140300514A1 (en) * 2012-12-05 2014-10-09 Samsung Electronics Co., Ltd. Method and apparatus for acquiring signal of global navigation satellite system
US20160116599A1 (en) * 2013-06-05 2016-04-28 Airbus Defence And Space Limited Receiver and method for direct sequence spread spectrum signals
US20150059432A1 (en) * 2013-08-27 2015-03-05 Crystal Instruments Corporation Cross-path phase calibration for high dynamic range data acquisition
CN103529458A (zh) * 2013-10-18 2014-01-22 北京邮电大学 卫星捕获方法及装置
CN103595441A (zh) * 2013-10-30 2014-02-19 中国运载火箭技术研究院 适用于重复使用轨道飞行器的极高动态扩频信号快捕方法
CN104360357A (zh) * 2014-11-24 2015-02-18 四川九洲电器集团有限责任公司 基于循环方式的北斗卫星信号快速捕获方法及系统
CN105846855A (zh) * 2016-02-28 2016-08-10 中国人民解放军63686部队 一种基于频率引导的大动态扩频信号快速捕获方法
CN106646546A (zh) * 2016-12-20 2017-05-10 南京六九零二科技有限公司 一种卫星信号多维快速捕获方法及系统
CN107135013A (zh) * 2017-05-04 2017-09-05 中国电子科技集团公司第五十四研究所 一种直接序列扩频信号快速捕获方法
CN109100747A (zh) * 2018-07-20 2018-12-28 西安空间无线电技术研究所 一种适用于多路随机接入信号的快速捕获系统及方法
CN110824512A (zh) * 2019-11-26 2020-02-21 中国科学院国家空间科学中心 一种非均匀码片实时延迟多普勒映像数据生成器
CN111694027A (zh) * 2020-06-04 2020-09-22 长沙北斗产业安全技术研究院有限公司 超大动态扩频信号捕获方法与装置
CN111934710A (zh) * 2020-07-06 2020-11-13 南京天际砺剑科技有限公司 一种高动态扩频信号快速捕获算法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
方科: "高动态低信噪比下扩频信号捕获算法研究", 电视技术, vol. 37, no. 13, 31 December 2013 (2013-12-31), pages 2 *
王鹏;陈国瑛;: "高动态导航接收机直扩信号捕获方法研究", 现代防御技术, no. 02, 15 April 2015 (2015-04-15), pages 2 - 3 *

Similar Documents

Publication Publication Date Title
CN106646546B (zh) 一种卫星信号多维快速捕获方法及系统
CN110501729B (zh) 一种基于fpga分步码相位细化的gnss信号的捕获方法
CN112118196B (zh) 信号捕获方法及系统
CN105158779A (zh) 一种改进的pmf_fft的pn码捕获方法
CN102914782A (zh) 一种适用于gps弱信号的快速捕获方法
JP2005031073A (ja) Ｇｐｓ相関ピーク信号の探索方法及びこれのためのシステム。
CN113972929B (zh) 一种高动态多普勒下扩频信号的捕获方法
CN106526631B (zh) 北斗b1频点卫星信号载波频率高精度估计方法
CN104111465A (zh) 连续波干扰信号的估计器、估计方法、消除器及消除方法
CN105807264A (zh) 雷达脉冲重复频率检测与初始脉冲到达时间的估计方法
CN109239743B (zh) 一种卫星信号捕获方法及装置
CN102707297A (zh) 一种compass信号的快速捕获方法
CN102798871B (zh) 基于伪码重构的伪码捕获方法及装置
CN101788677A (zh) 全球卫星导航系统信号捕获扫频系统及其捕获方法
CN102655419A (zh) 一种捕获扩频信号的判决函数计算方法
CN117214928A (zh) 一种伪码跟踪方法、装置、终端设备及存储介质
CN112649819A (zh) 一种高动态扩频信号捕获装置及捕获方法
CN106772470A (zh) 卫星导航多频接收机多频点时延检测和校正方法、装置
CN109633707B (zh) 一种基于预平均处理的变系数匹配滤波的快速捕获方法
CN110958035A (zh) 一种基于三角包络拟合的频率二次捕获装置及方法
CN110007322A (zh) 基于相干降采样的北斗b1i信号捕获方法
CN103869343B (zh) 捕获gps中l2c信号的方法及装置
CN112764063A (zh) 一种实现捕获处理的方法及接收机
CN113671547B (zh) 一种改进的高动态捕获方法、装置、设备及存储介质
CN110061760B (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