CN106556848A - 一种北斗二代卫星b1频点信号的快速捕获方法 - Google Patents

一种北斗二代卫星b1频点信号的快速捕获方法 Download PDF

Info

Publication number
CN106556848A
CN106556848A CN201610931756.7A CN201610931756A CN106556848A CN 106556848 A CN106556848 A CN 106556848A CN 201610931756 A CN201610931756 A CN 201610931756A CN 106556848 A CN106556848 A CN 106556848A
Authority
CN
China
Prior art keywords
signal
frequency
dft
depositor
zero
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
CN201610931756.7A
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.)
Huazhong Institute Of Optoelectronic Technology (china Shipbuilding Industry Corp 717 Institute)
717th Research Institute of CSIC
Original Assignee
Huazhong Institute Of Optoelectronic Technology (china Shipbuilding Industry Corp 717 Institute)
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 Huazhong Institute Of Optoelectronic Technology (china Shipbuilding Industry Corp 717 Institute) filed Critical Huazhong Institute Of Optoelectronic Technology (china Shipbuilding Industry Corp 717 Institute)
Priority to CN201610931756.7A priority Critical patent/CN106556848A/zh
Publication of CN106556848A publication Critical patent/CN106556848A/zh
Pending legal-status Critical Current

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

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

本发明设计一种北斗二代卫星B1频点信号的快速捕获方法,其步骤包含:步骤1、产生所有北斗卫星B1频点的C/A码序列,形成正弦表和余弦表;步骤2、实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号;步骤3、将零频基带数据信号降速率存储、处理;步骤4、捕获北斗二代卫星信号中的有效可见卫星信号。本发明的优点在于,在用较少的资源消耗、较低复杂度的硬件电路实现的基础上,实现对信号的快速捕获。

Description

一种北斗二代卫星B1频点信号的快速捕获方法
技术领域
本发明涉及北斗卫星导航与定位技术,具体涉及一种北斗二代卫星B1频点信号的快速捕获方法的优化。
背景技术
卫星导航定位是通过接收机天线接收卫星发射的电磁波信号,从接收到的电磁波信号中搜索捕获各可见卫星信号,获得粗略的码相位和载波多普勒,进而将这些捕获参数送入跟踪通道实现对卫星的跟踪,从而获取各种观测量,导航解算模块根据这些观测量实现对接收机的导航定位。因此,实现对卫星信号的成功捕获是卫星导航接收机技术中的最基本的功能,同时也是一项关键技术。
信号捕获的基本过程是使用本地复制产生的C/A码与接收到的卫星信号进行相关运算,根据相关运算的结果判断卫星信号是否捕获成功。根据相关运算实现方式的不同,卫星信号的捕获搜索算法可分为以下三类:线性搜索捕获算法、并行频率搜索捕获算法和并行码相位搜索捕获算法。
线性搜索捕获是指在捕获过程中,对指定卫星信号的码相位和多普勒频移进行扫描式遍历搜索,此方法硬件实现简单,资源消耗少,但是由于需要对所有可能的码相位和多普勒频移进行搜索,导致该方法捕获速度慢,实时性差;并行捕获则对每种相位均设置了相关器,同时做相关运算,将与设定的相关值相同的相位作为捕获成功的相位,和串行捕获相比,并行频率捕获的捕获时间缩短了,需要进行一次傅里叶变换,对硬件资源要求较多,实现起来稍复杂;并行码相位捕获算法时间大大缩短,但是需要进行一次傅里叶变换和一次傅里叶逆变换,电路变得更加复杂,运算量也大大增加,消耗的资源也进一步增多。
发明内容
本发明的目的在于,为了在较少的资源消耗、较低复杂度的硬件电路实现的基础上,实现对信号的快速捕获,本发明提出一种基于并行频率搜索的快速捕获算法。
本发明的技术方案是:
一种北斗二代卫星B1频点信号的快速捕获方法,其步骤包含:
步骤1、产生所有北斗卫星B1频点的C/A码序列,形成正弦表和余弦表:用MATLAB编程产生所有北斗卫星B1频点的C/A码序列,将这些二进制序列存入寄存器一,用MATLAB编程产生正弦波一个周期内的64个离散值,存入寄存器二中,形成正弦表;将正弦波相位滞后90°,取一个周期内的64个离散值,存入寄存器三中,形成余弦表;
步骤2、实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号:通过北斗微带天线接收北斗二代B1频点信号,经过射频前端处理,降频成0.098MHz的近零中频信号,再经过AD正交采样得到零频基带数据,并送入FPGA中,所述零频基带数据由I支路的信号I和Q支路的信号Q组成,信号I存入寄存器四中,该信号I为同相分量信号,信号Q存入寄存器五中,该信号Q为正交信号;在FPGA中进一步进行相位旋转下变频处理,实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号;
步骤3、将零频基带数据信号降速率存储、处理:对零频基带数据信号I和信号Q分别进行累积降速处理,生成降速零频基带数据信号,所述降速零频基带数据信号由I支路的信号I’和Q支路的信号Q’组成;
步骤4、捕获北斗二代卫星信号中的有效可见卫星信号:采用并行频率快速捕获算法来捕获北斗二代卫星信号中的有效可见卫星信号;具体操作时,采用基于循环相关的方式对步骤3所得的信号I’和信号Q’进行并行快速捕获。
进一步的技术方案是:
所述的北斗二代B1频点信号的快速捕获方法,其步骤3中,降速率存储步骤如下:
步骤3.1、利用数字频率合成技术,产生B1频点伪码基频速率时钟使能信号,即时钟频率为2.046MHz的时钟使能信号C1;
步骤3.2、通过使能信号C1对零频基带数据信号I和信号Q分别进行累积降速处理,截取高位,生成2.046MHz数据率的降速零频基带数据信号I’和信号Q’;
步骤3.3、将降速零频基带数据信号I’实时缓存于双端口寄存器六中,降速零频基带数据信号Q’实时缓存于双端口寄存器七中,存储数据点数都为4092点,即2ms时长的降速零频基带数据信号,存储器深度都为4096;
步骤3.4、双端口寄存器六和双端口寄存器七存满之后,重新从起始地址开始存储,覆盖之前存储的信号,这样每个存储的数据有且只包含一个C/A码码片。
所述的北斗二代B1频点信号的快速捕获方法,其步骤4中并行频率快速捕获的步骤如下:
步骤4.1、在双端口寄存器六存满2ms降速零频基带数据信号I’后,从起始位置读取双端口寄存器六中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;同时,在双端口寄存器七存满2ms降速零频基带数据信号Q’后,从起始位置读取双端口寄存器七中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;
步骤4.2、将相关累加结果每隔64点累加出一个点来,并对寄存器六、和寄存器七中的4092点后面补零扩展至4096点,这样共获得I支路64点相关累加值和Q支路64点相关累加值;
步骤4.3、在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器八和双端口寄存器九中,其中,双端口寄存器八存储的I支路的DFT运算结果,双端口寄存器九寄存器存储的是Q支路的DFT运算结果;在获得第二点累加值时,对其进行第二次64点DFT运算,并把结果累加到第一次DFT运算结果中,I支路的结果累加到I支路运算结果中,Q支路的结果累加到Q支路的运算结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果。
步骤4.4、实时计算DFT峰值能量及平均能量,并缓存到寄存器十中;
步骤4.5、重新读取双口寄存器六中的信号I’和双口寄存器七中的信号Q’的数据并与本地伪码序列延迟一个数据点对齐并进行相关累加,获得64点DFT结果,计算DFT峰值能量及平均能量,这样重复2046次,记录64点DFT中峰值能量最大的一次,即可获得伪码初相、DFT峰值能量最大处频率点以及峰值能量和平均能量等信息。
步骤4.6、将峰值能量和平均能量送入捕获判决模块中,比较峰值能量和平均能量,若峰值能量大于或等于16倍的平均能量,则判断捕获到卫星信号,否则,未捕获到卫星信号,捕获通道重新开始新一轮捕获过程。
所述的北斗二代B1频点信号的快速捕获方法,其步骤2中,射频前端将信号频率下变频至0.098MHz,成为近零中频信号。
所述的北斗二代B1频点信号的快速捕获方法,其步骤3中,零频基带数据信号被降速率至2.046MHz,并被按存储器地址顺序实时存储于深度为4096的双口寄存器中,存储时长为2ms,在存储2ms零频基带数据后,回到寄存器的起始地址重新按地址顺序实时存储2ms数据。
所述的北斗二代B1频点信号的快速捕获方法,其步骤4中,在存储2ms零频基带数据信号数据后,在系统时钟的驱动下读取双口寄存器中的数据与本地C/A序列对齐并进行相关累加,每64个相关值累加出一个点,并对双口寄存器中的4092点后面补零扩展至4096点,这样共获得64点相关累加值。
所述的北斗二代B1频点信号的快速捕获方法,其步骤4.3中,在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器中,在获得第二点累加值时,对齐进行第二次64点DFT运算,并把结果累加到第一次结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果。
所述的北斗二代B1频点信号的快速捕获方法,其DFT运算模块包括两个DFT系数存储器,一个DFT系数存储器地址发生器,两个DFT运算结果存储器,一个DFT运算结果存储器地址发生器。
所述的北斗二代B1频点信号的快速捕获方法,所述两个DFT系数存储器都有64个存储单元,每个单元存储一个DFT系数,这两个DFT系数存储器是步骤1中的正弦表和余弦表,DFT系数就是正弦表和余弦表中的离散值。
所述的北斗二代B1频点信号的快速捕获方法,所述步骤4.6中,判决捕获成功的条件是峰值能量数值大于或等于16倍的平均能量数值。
所述的北斗二代B1频点信号的快速捕获方法,其步骤2中所述的载波剥离其过程如下:
步骤2.1 AD正交采样送入FPGA的基带数据信号实时缓存于寄存器四和寄存器五;
步骤2.2 在系统时钟驱动下,同时读取寄存器四及寄存器五中的数据和寄存器二及寄存器三中的离散值,通过相位旋转法进行下变频处理,实现载波剥离。
所述的北斗二代B1频点信号的快速捕获方法,其步骤2中,北斗天线接收北斗卫星的B1频点的射频信号,经硬件下变频处理,信号频率从1561.098MHz下降到了0.098MHz,成为近零中频信号;近零中频信号表达式如下:
式中,d(t)为导航电文,p(t)为为C/A码序列,ωc为包含载波多普勒频率在内的近零中频信号的频率。
所述的北斗二代B1频点信号的快速捕获方法,其步骤2中,近零中频信号的载波剥离实现如下:数控振荡器在相位增量和系统时钟的驱动下产生地址发生器,根据产生的地址读取正弦表、余弦表中,即寄存器二和寄存器三中预先存储的数据,数控振荡器的输出用复数表示为:
相位旋转的过程可以表示成下式:
则相位旋转后的同相分量I、正交分量Q信号分别为:
当选择ωn=-0.098MHz时,则可完全消除载波信息,只留下载波多普勒频率、导航电文和伪码信息。
所述的北斗二代B1频点信号的快速捕获方法,其步骤4.3中,进行DFT运算的实现方法如下:DFT变换式如下:
进一步变换为:
对于固定点数的DFT,变换系数cos(2πnk/N)和-sin(2πnk/N)预先量化后存储于寄存器二和寄存器三中,在做DFT运算时,通过时钟信号产生地址发生器读取相应地址的寄存器中的变换系数来参与DFT运算。
本发明的有益效果在于,本发明充分利用了并行频率搜索时DFT运算所需的DFT系数可以与载波剥离时的数控振荡器的正弦表和余弦表可以共用的特点,在减少硬件资源消耗、较低复杂度的硬件电路实现的基础上,实现了对北斗二代B1频点信号的实时快速捕获。
附图说明
图1是本发明的总体方框图;
图2是在图1技术方案的基础上进一步方案的方框图;
图3是近零中频信号下变频到零频信号流程图;
图4是DFT运算模块。
具体实施方式
结合附图和实施例对本发明作进一步说明如下:
实施例1、是本发明的一个基本实施例。如图1、2所示,一种北斗二代卫星B1频点信号的快速捕获方法,其步骤包含:
步骤1、产生所有北斗卫星B1频点的C/A码序列,形成正弦表和余弦表:用MATLAB编程产生所有北斗卫星B1频点的C/A码序列,将这些二进制序列存入寄存器一,用MATLAB编程产生正弦波一个周期内的64个离散值,存入寄存器二中,形成正弦表;将正弦波相位滞后90°,取一个周期内的64个离散值,存入寄存器三中,形成余弦表;
步骤2、实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号:通过北斗微带天线接收北斗二代B1频点信号,经过射频前端处理,降频成0.098MHz的近零中频信号,再经过AD正交采样得到零频基带数据,并送入FPGA中,所述零频基带数据由I支路的信号I和Q支路的信号Q组成,信号I存入寄存器四中,该信号I为同相分量信号,信号Q存入寄存器五中,该信号Q为正交信号;在FPGA中进一步进行相位旋转下变频处理,实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号;相位旋转下变频处理过程如图3所示;
步骤3、将零频基带数据信号降速率存储、处理:对零频基带数据信号I和信号Q分别进行累积降速处理,生成降速零频基带数据信号,所述降速零频基带数据信号由I支路的信号I’和Q支路的信号Q’组成;
步骤4、捕获北斗二代卫星信号中的有效可见卫星信号:采用并行频率快速捕获算法来捕获北斗二代卫星信号中的有效可见卫星信号;具体操作时,采用基于循环相关的方式对步骤3所得的信号I’和信号Q’进行并行快速捕获。
实施例2、是在实施例1基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,其步骤3中,降速率存储步骤如下:
步骤3.1、利用数字频率合成技术,产生B1频点伪码基频速率时钟使能信号,即时钟频率为2.046MHz的时钟使能信号C1;
步骤3.2、通过使能信号C1对零频基带数据信号I和信号Q分别进行累积降速处理,截取高位,生成2.046MHz数据率的降速零频基带数据信号I’和信号Q’;
步骤3.3、将降速零频基带数据信号I’实时缓存于双端口寄存器六中,降速零频基带数据信号Q’实时缓存于双端口寄存器七中,存储数据点数都为4092点,即2ms时长的降速零频基带数据信号,存储器深度都为4096;
步骤3.4、双端口寄存器六和双端口寄存器七存满之后,重新从起始地址开始存储,覆盖之前存储的信号,这样每个存储的数据有且只包含一个C/A码码片。
实施例3、是在实施例1基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,步骤4中,并行频率快速捕获的步骤如下:
步骤4.1、在双端口寄存器六存满2ms降速零频基带数据信号I’后,从起始位置读取双端口寄存器六中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;同时,在双端口寄存器七存满2ms降速零频基带数据信号Q’后,从起始位置读取双端口寄存器七中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;
步骤4.2、将相关累加结果每隔64点累加出一个点来,并对寄存器六、和寄存器七中的4092点后面补零扩展至4096点,这样共获得I支路64点相关累加值和Q支路64点相关累加值;
步骤4.3、在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器八和双端口寄存器九中,其中,双端口寄存器八存储的I支路的DFT运算结果,双端口寄存器九寄存器存储的是Q支路的DFT运算结果;在获得第二点累加值时,对其进行第二次64点DFT运算,并把结果累加到第一次DFT运算结果中,I支路的结果累加到I支路运算结果中,Q支路的结果累加到Q支路的运算结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果;每个支路的DFT运算过程如图4所示;
步骤4.4、实时计算DFT峰值能量及平均能量,并缓存到寄存器十中;
步骤4.5、重新读取双口寄存器六中的信号I’和双口寄存器七中的信号Q’的数据并与本地伪码序列延迟一个数据点对齐并进行相关累加,获得64点DFT结果,计算DFT峰值能量及平均能量,这样重复2046次,记录64点DFT中峰值能量最大的一次,即可获得伪码初相、DFT峰值能量最大处频率点以及峰值能量和平均能量等信息;
步骤4.6、将峰值能量和平均能量送入捕获判决模块中,比较峰值能量和平均能量,若峰值能量大于或等于16倍的平均能量,则判断捕获到卫星信号,否则,未捕获到卫星信号,捕获通道重新开始新一轮捕获过程。
实施例4、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述步骤2中,射频前端将信号频率下变频至0.098MHz,成为近零中频信号。
实施例5、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述步骤3中,零频基带数据信号被降速率至2.046MHz,并被按存储器地址顺序实时存储于深度为4096的双口寄存器中,存储时长为2ms,在存储2ms零频基带数据后,回到寄存器的起始地址重新按地址顺序实时存储2ms数据。
实施例6、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,步骤4中,在存储2ms零频基带数据信号数据后,在系统时钟的驱动下读取双口寄存器中的数据与本地C/A序列对齐并进行相关累加,每64个相关值累加出一个点,并对双口寄存器中的4092点后面补零扩展至4096点,这样共获得64点相关累加值。
实施例7、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,步骤4.3中,在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器中,在获得第二点累加值时,对齐进行第二次64点DFT运算,并把结果累加到第一次结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果。如图4所示。
实施例8、是在实施例7基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,DFT运算模块包括两个DFT系数存储器,一个DFT系数存储器地址发生器,两个DFT运算结果存储器,一个DFT运算结果存储器地址发生器。所述两个DFT系数存储器都有64个存储单元,每个单元存储一个DFT系数,这两个DFT系数存储器是步骤1中的正弦表和余弦表,DFT系数就是正弦表和余弦表中的离散值。
实施例9、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,其所述步骤4.6中,判决捕获成功的条件是峰值能量数值大于或等于16倍的平均能量数值。
实施例10、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,步骤2中所述的载波剥离其过程如下:
步骤2.1 AD正交采样送入FPGA的基带数据信号实时缓存于寄存器四和寄存器五;
步骤2.2 在系统时钟驱动下,同时读取寄存器四及寄存器五中的数据和寄存器二及寄存器三中的离散值,通过相位旋转法进行下变频处理,实现载波剥离。
实施例11、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。步骤2中,北斗天线接收北斗卫星的B1频点的射频信号,经硬件下变频处理,信号频率从1561.098MHz下降到了0.098MHz,成为近零中频信号;近零中频信号表达式如下:
式中,d(t)为导航电文,p(t)为为C/A码序列,ωc为包含载波多普勒频率在内的近零中频信号的频率。
步骤2中,近零中频信号的载波剥离实现如下:数控振荡器在相位增量和系统时钟的驱动下产生地址发生器,根据产生的地址读取正弦表、余弦表中,即寄存器二和寄存器三中预先存储的数据,数控振荡器的输出用复数表示为:
相位旋转的过程可以表示成下式:
则相位旋转后的同相分量I、正交分量Q信号分别为:
当选择ωn=-0.098MHz时,则可完全消除载波信息,只留下载波多普勒频率、导航电文和伪码信息。
实施例12、是在实施例1基础上进一步的实施例,还可以是在实施例2或3基础上进一步的实施例。所述的北斗二代B1频点信号的快速捕获方法,其步骤4.3中,进行DFT运算的实现方法如下:DFT变换式如下:
进一步变换为:
对于固定点数的DFT,变换系数cos(2πnk/N)和-sin(2πnk/N)预先量化后存储于寄存器二和寄存器三中,在做DFT运算时,通过时钟信号产生地址发生器读取相应地址的寄存器中的变换系数来参与DFT运算。
本发明的权利要求保护范围不限于上述实施例。

Claims (14)

1.一种北斗二代卫星B1频点信号的快速捕获方法,其特征在于,其步骤包含:
步骤1、产生所有北斗卫星B1频点的C/A码序列,形成正弦表和余弦表:用MATLAB编程产生所有北斗卫星B1频点的C/A码序列,将这些二进制序列存入寄存器一,用MATLAB编程产生正弦波一个周期内的64个离散值,存入寄存器二中,形成正弦表;将正弦波相位滞后90°,取一个周期内的64个离散值,存入寄存器三中,形成余弦表;
步骤2、实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号:通过北斗微带天线接收北斗二代B1频点信号,经过射频前端处理,降频成0.098MHz的近零中频信号,再经过AD正交采样得到零频基带数据,并送入FPGA中,所述零频基带数据由I支路的信号I和Q支路的信号Q组成,信号I存入寄存器四中,该信号I为同相分量信号,信号Q存入寄存器五中,该信号Q为正交信号;在FPGA中进一步进行相位旋转下变频处理,实现对载波的剥离,得到只包含载波多普勒的零频基带数据信号;
步骤3、将零频基带数据信号降速率存储、处理:对零频基带数据信号I和信号Q分别进行累积降速处理,生成降速零频基带数据信号,所述降速零频基带数据信号由I支路的信号I’和Q支路的信号Q’组成;
步骤4、捕获北斗二代卫星信号中的有效可见卫星信号:采用并行频率快速捕获算法来捕获北斗二代卫星信号中的有效可见卫星信号;具体操作时,采用基于循环相关的方式对步骤3所得的信号I’和信号Q’进行并行快速捕获。
2.根据权利要求1所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤3中,降速率存储步骤如下:
步骤3.1、利用数字频率合成技术,产生B1频点伪码基频速率时钟使能信号,即时钟频率为2.046MHz的时钟使能信号C1;
步骤3.2、通过使能信号C1对零频基带数据信号I和信号Q分别进行累积降速处理,截取高位,生成2.046MHz数据率的降速零频基带数据信号I’和信号Q’;
步骤3.3、将降速零频基带数据信号I’实时缓存于双端口寄存器六中,降速零频基带数据信号Q’实时缓存于双端口寄存器七中,存储数据点数都为4092点,即2ms时长的降速零频基带数据信号,存储器深度都为4096;
步骤3.4、双端口寄存器六和双端口寄存器七存满之后,重新从起始地址开始存储,覆盖之前存储的信号,这样每个存储的数据有且只包含一个C/A码码片。
3.根据权利要求1所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤4中并行频率快速捕获的步骤如下:
步骤4.1、在双端口寄存器六存满2ms降速零频基带数据信号I’后,从起始位置读取双端口寄存器六中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;同时,在双端口寄存器七存满2ms降速零频基带数据信号Q’后,从起始位置读取双端口寄存器七中的数据,同时读取寄存器一中的本地C/A码,两者对齐并进行相关累加;
步骤4.2、将相关累加结果每隔64点累加出一个点来,并对寄存器六和寄存器七中的4092点后面补零扩展至4096点,这样共获得I支路64点相关累加值和Q支路64点相关累加值;
步骤4.3、在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器八和双端口寄存器九中,其中,双端口寄存器八存储的I支路的DFT运算结果,双端口寄存器九寄存器存储的是Q支路的DFT运算结果;在获得第二点累加值时,对其进行第二次64点DFT运算,并把结果累加到第一次DFT运算结果中,I支路的结果累加到I支路运算结果中,Q支路的结果累加到Q支路的运算结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果;
步骤4.4、实时计算DFT峰值能量及平均能量,并缓存到寄存器十中;
步骤4.5、重新读取双口寄存器六中的信号I’和双口寄存器七中的信号Q’的数据并与本地伪码序列延迟一个数据点对齐并进行相关累加,获得64点DFT结果,计算DFT峰值能量及平均能量,这样重复2046次,记录64点DFT中峰值能量最大的一次,即可获得伪码初相、DFT峰值能量最大处频率点以及峰值能量和平均能量等信息;
步骤4.6、将峰值能量和平均能量送入捕获判决模块中,比较峰值能量和平均能量,若峰值能量大于或等于16倍的平均能量,则判断捕获到卫星信号,否则,未捕获到卫星信号,捕获通道重新开始新一轮捕获过程。
4.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,所述步骤2中,射频前端将信号频率下变频至0.098MHz,成为近零中频信号。
5.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,所述步骤3中,零频基带数据信号被降速率至2.046MHz,并被按存储器地址顺序实时存储于深度为4096的双口寄存器中,存储时长为2ms,在存储2ms零频基带数据后,回到寄存器的起始地址重新按地址顺序实时存储2ms数据。
6.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤4中,在存储2ms零频基带数据信号数据后,在系统时钟的驱动下读取双口寄存器中的数据与本地C/A序列对齐并进行相关累加,每64个相关值累加出一个点,并对双口寄存器中的4092点后面补零扩展至4096点,这样共获得64点相关累加值。
7.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤4.3中,在系统时钟驱动下每获得一点累加值便对其进行一次64点DFT运算,并把结果存储于双端口寄存器中,在获得第二点累加值时,对齐进行第二次64点DFT运算,并把结果累加到第一次结果中,当做完第64次DFT运算时,即获得了完整的64点DFT运算结果。
8.根据权利要求7所述的北斗二代B1频点信号的快速捕获方法,其特征在于,DFT运算模块包括两个DFT系数存储器,一个DFT系数存储器地址发生器,两个DFT运算结果存储器,一个DFT运算结果存储器地址发生器。
9.根据权利要求8所述的北斗二代B1频点信号的快速捕获方法,其特征在于,所述两个DFT系数存储器都有64个存储单元,每个单元存储一个DFT系数,这两个DFT系数存储器是步骤1中的正弦表和余弦表,DFT系数就是正弦表和余弦表中的离散值。
10.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,所述步骤4.6中,判决捕获成功的条件是峰值能量数值大于或等于16倍的平均能量数值。
11.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤2中所述的载波剥离其过程如下:
步骤2.1 AD正交采样送入FPGA的基带数据信号实时缓存于寄存器四和寄存器五;
步骤2.2在系统时钟驱动下,同时读取寄存器四及寄存器五中的数据和寄存器二及寄存器三中的离散值,通过相位旋转法进行下变频处理,实现载波剥离。
12.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤2中,北斗天线接收北斗卫星的B1频点的射频信号,经硬件下变频处理,信号频率从1561.098MHz下降到了0.098MHz,成为近零中频信号;近零中频信号表达式如下:
式中,d(t)为导航电文,p(t)为为C/A码序列,ωc为包含载波多普勒频率在内的近零中频信号的频率。
13.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤2中,近零中频信号的载波剥离实现如下:数控振荡器在相位增量和系统时钟的驱动下产生地址发生器,根据产生的地址读取正弦表、余弦表中,即寄存器二和寄存器三中预先存储的数据,数控振荡器的输出用复数表示为:
相位旋转的过程可以表示成下式:
则相位旋转后的同相分量I、正交分量Q信号分别为:
I EMBED Equation.3*MERGEFORMAT I.
当选择ωn=-0.098MHz时,则可完全消除载波信息,只留下载波多普勒频率、导航电文和伪码信息。
14.根据权利要求1或2或3所述的北斗二代B1频点信号的快速捕获方法,其特征在于,步骤4.3中,进行DFT运算的实现方法如下:DFT变换式如下:
进一步变换为:
对于固定点数的DFT,变换系数cos(2πnk/N)和-sin(2πnk/N)预先量化后存储于寄存器二和寄存器三中,在做DFT运算时,通过时钟信号产生地址发生器读取相应地址的寄存器中的变换系数来参与DFT运算。
CN201610931756.7A 2016-10-24 2016-10-24 一种北斗二代卫星b1频点信号的快速捕获方法 Pending CN106556848A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610931756.7A CN106556848A (zh) 2016-10-24 2016-10-24 一种北斗二代卫星b1频点信号的快速捕获方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610931756.7A CN106556848A (zh) 2016-10-24 2016-10-24 一种北斗二代卫星b1频点信号的快速捕获方法

Publications (1)

Publication Number Publication Date
CN106556848A true CN106556848A (zh) 2017-04-05

Family

ID=58443160

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610931756.7A Pending CN106556848A (zh) 2016-10-24 2016-10-24 一种北斗二代卫星b1频点信号的快速捕获方法

Country Status (1)

Country Link
CN (1) CN106556848A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110007322A (zh) * 2019-04-17 2019-07-12 南京航空航天大学 基于相干降采样的北斗b1i信号捕获方法
CN110320539A (zh) * 2018-03-30 2019-10-11 华为技术有限公司 一种应用于卫星定位系统的比特同步方法以及相关装置
CN112596082A (zh) * 2020-11-16 2021-04-02 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) 北斗b3i信号的载波频率牵引方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2554412A1 (en) * 2004-01-28 2005-08-11 Qualcomm Incorporated Rapid acquisition methods and apparatus for gps signals background
CN102901973A (zh) * 2012-09-25 2013-01-30 郑州威科姆科技股份有限公司 基于北斗卫星的信号实时并行快速捕获方法
CN104297773A (zh) * 2014-02-27 2015-01-21 北京航天时代光电科技有限公司 一种高精度北斗三频sins深组合导航系统
CN104483684A (zh) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 一种快速捕获北斗d1卫星导航系统弱信号的方法
CN105866803A (zh) * 2016-03-23 2016-08-17 沈阳航空航天大学 基于fpga北斗二代卫星导航接收机基带信号快捕算法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2554412A1 (en) * 2004-01-28 2005-08-11 Qualcomm Incorporated Rapid acquisition methods and apparatus for gps signals background
CN102901973A (zh) * 2012-09-25 2013-01-30 郑州威科姆科技股份有限公司 基于北斗卫星的信号实时并行快速捕获方法
CN104297773A (zh) * 2014-02-27 2015-01-21 北京航天时代光电科技有限公司 一种高精度北斗三频sins深组合导航系统
CN104483684A (zh) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 一种快速捕获北斗d1卫星导航系统弱信号的方法
CN105866803A (zh) * 2016-03-23 2016-08-17 沈阳航空航天大学 基于fpga北斗二代卫星导航接收机基带信号快捕算法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320539A (zh) * 2018-03-30 2019-10-11 华为技术有限公司 一种应用于卫星定位系统的比特同步方法以及相关装置
CN110320539B (zh) * 2018-03-30 2024-01-30 华为技术有限公司 一种应用于卫星定位系统的比特同步方法以及相关装置
CN110007322A (zh) * 2019-04-17 2019-07-12 南京航空航天大学 基于相干降采样的北斗b1i信号捕获方法
CN110007322B (zh) * 2019-04-17 2022-06-10 南京航空航天大学 基于相干降采样的北斗b1i信号捕获方法
CN112596082A (zh) * 2020-11-16 2021-04-02 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) 北斗b3i信号的载波频率牵引方法及装置
CN112596082B (zh) * 2020-11-16 2023-12-26 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) 北斗b3i信号的载波频率牵引方法及装置

Similar Documents

Publication Publication Date Title
CN102282773B (zh) 用于捕捉扩频信号的方法和系统
US8301678B2 (en) Method and device for fast correlation calculation
CN105301608B (zh) 高动态接收机及其对北斗卫星信号b1频点的捕获方法
CN104360357A (zh) 基于循环方式的北斗卫星信号快速捕获方法及系统
CN105717522B (zh) “北斗二代”b1频段弱信号捕获方法
CN110045397B (zh) 一种基于fpga的l5信号捕获方法及装置
CN106556848A (zh) 一种北斗二代卫星b1频点信号的快速捕获方法
CN106788470B (zh) 基于扩频技术的高数据速率突发通信快速捕获方法及装置
CN101952735B (zh) 多相位码发生器和gnss接收机
CN102928854A (zh) 一种基于匹配滤波器的gps捕获单元设计方法
CN100581068C (zh) 一种多星座卫星信号的混合串并行相关捕获方法
CN102353968A (zh) 一种基于fpga的gps信号捕获方法及其系统
CN108169772A (zh) 一种加窗fft的卫星信号捕获方法
CN103698784A (zh) 一种p码捕获方法及装置
CN104849732B (zh) 一种二进制偏移载波射频导航信号跟踪方法
CN106772475A (zh) 一种基于fpga捕获微弱gps导航信号的电路及方法
CN106526631B (zh) 北斗b1频点卫星信号载波频率高精度估计方法
CN103941269A (zh) 用于卫星导航系统的pn码捕获方法
CN103439718A (zh) 一种高阶boc调制信号的无模糊跟踪单元
CN101320084A (zh) 卫星定位信号的实时处理方法
CN106093981A (zh) 基于优化并行码相位搜索的gps捕获电路
JP2009519448A (ja) 信号処理装置及び信号処理方法
CN101788677A (zh) 全球卫星导航系统信号捕获扫频系统及其捕获方法
CN201152895Y (zh) 一种gps信号相关器
CN104833990B (zh) 基于fpga的北斗卫星导航系统基带信号预处理方法和装置

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170405

WD01 Invention patent application deemed withdrawn after publication