CN101866008A - 多卫星导航下的反射信号接收机的相关器 - Google Patents

多卫星导航下的反射信号接收机的相关器 Download PDF

Info

Publication number
CN101866008A
CN101866008A CN201010193964A CN201010193964A CN101866008A CN 101866008 A CN101866008 A CN 101866008A CN 201010193964 A CN201010193964 A CN 201010193964A CN 201010193964 A CN201010193964 A CN 201010193964A CN 101866008 A CN101866008 A CN 101866008A
Authority
CN
China
Prior art keywords
industrial computer
data processing
sign indicating
indicating number
pattern
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
Application number
CN201010193964A
Other languages
English (en)
Other versions
CN101866008B (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201010193964A priority Critical patent/CN101866008B/zh
Publication of CN101866008A publication Critical patent/CN101866008A/zh
Application granted granted Critical
Publication of CN101866008B publication Critical patent/CN101866008B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种多卫星导航下的反射信号接收机的相关器，该反射信号接收机的相关器包括有模式配置单元(1)、反射信号处理单元(2)和直射信号处理单元(3)；所述反射信号处理单元(2)由载波解调器(201)、延迟CA码发生器(202)、延迟CA码相关运算器(203)和1ms相干累加器(204)构成。本发明的相关器是在FPGA芯片，采用VerilogHDL语言编程来实现多导航卫星下反射信号的相关处理和多种模式配置。配置有本发明相关器的反射信号接收机能够为后续物理状态的反演提供更加合理的、灵活的数据源，并为星载装置的研发提供技术基础。

Description

多卫星导航下的反射信号接收机的相关器
技术领域
本发明涉及一种适用于接收机中的相关器设计，更特别地说，是指一种能同时接收处理多种导航卫星的反射信号，并通过多种模式配置来实现不同功能的反射信号的接收与处理。
背景技术
目前已建成的全球卫星导航系统包括美国的GPS(Global PositioningSystem)、俄罗斯的GLONASS(Global Orbiting Navigation Satellite System)，以及中国的北斗一号(BD-1)区域卫星导航系统，随着GPS、GLONASS现代化计划的实施，将提供更多的信号源；正在建设的系统包括中国北斗二号(BD-2)卫星导航系统、欧盟Galileo系统、印度的卫星导航定位系统、日本的准天顶卫星系统(QZSS，Quasi-Zenith Satellite System)等。
不同的全球卫星导航系统在卫星数量、运行轨道和工作参数方面虽然存在差异，但在用于微波遥感的信号源结构、L波段频率、信号功率等方面相当接近。因此，GNSS-R遥感探测机理适于多系统信号源的兼容使用。
在2008年11月第29卷第6期中“遥测遥控”公开的《GNSS-R延迟映射接收机相关器设计》，如图1所示，讲解了一种延迟映射接收机的总体结构；同时也公开了串行工作模式的相关器和并行工作模式的相关器。
现有延迟映射接收机只能处理GPS单系统，并且模式配置简单，对于不同场合的应用以及操作的灵活性不能实现多种模式的配置。
发明内容
为了解决目前接收与处理导航卫星反射信号装置的单一性，模式配置简单，对于不同场合的应用难以灵活性操作的缺陷，本发明提出了一种能够兼容多卫星导航下对接收的反射信号进行模式配置的相关器。该反射信号接收机的相关器采用在大容量FPGA芯片上通过VerilogHDL编程设计进行模式配置得以实现，为后续物理状态的反演提供更加合理和灵活的数据支持。
本发明的一种多卫星导航下的反射信号接收机的相关器，所述反射信号接收机主要由左旋天线、右旋天线、射频前端、FPGA相关器、DSP处理器、数据处理工控机、任务监控工控机组成；数据处理工控机向FPGA相关器发送模式选择控制指令，同时用于存储FPGA相关器输出的相关累加值
Figure BSA00000137769300021
m表示每一个反射通道的相关累加值的个数，n表示反射通道号；DSP处理器向FPGA相关器发出卫星码相位信息F＝{numPn，SVn，SVPn}，numPn表示卫星反射信号CA码相对于直射信号的码相位延迟数，SVn表示卫星编号，SVPn表示卫星直射信号跟踪后的CA码码相位，n表示反射通道号；其特征在于：该反射信号接收机的相关器包括有模式配置单元(1)、反射信号处理单元(2)和直射信号处理单元(3)；所述反射信号处理单元(2)由载波解调器(201)、延迟CA码发生器(202)、延迟CA码相关运算器(203)和1ms相干累加器(204)构成；
模式配置单元1中存储的模式有：
模式一：当模式配置单元1收到由数据处理工控机下发的0x11H指令后，一方面向数据处理工控机发送应答指令0x12H，另一方面向直射信号处理单元3发出启动指令，记为J1，然后直射信号处理单元3向数据处理工控机输出GPS原始数据；
模式二：当模式配置单元1收到由数据处理工控机下发的0x21H指令后，一方面向数据处理工控机发送应答指令0x22H，另一方面向直射信号处理单元3发出启动指令，记为J2，然后直射信号处理单元3向数据处理工控机输出Galileo原始数据；
模式三：当模式配置单元1收到由数据处理工控机下发的0x31H指令后，一方面向数据处理工控机发送应答指令0x32H，另一方面向直射信号处理单元3发出启动指令，记为J3，然后直射信号处理单元3向数据处理工控机输出GLONASS原始数据；
模式四：当模式配置单元1收到由数据处理工控机下发的0x41H指令后，一方面向数据处理工控机发送应答指令0x42H，另一方面向直射信号处理单元3发出启动指令，记为J4，然后直射信号处理单元3向数据处理工控机输出BD原始数据；
模式五：当模式配置单元1收到由数据处理工控机下发的0x5EH指令后，一方面向数据处理工控机发送应答指令0x5FH，另一方面向直射信号处理单元3发出启动指令，记为J5，然后直射信号处理单元3向数据处理工控机输出GPS、Galileo、GLONASS和BD混合的原始数据；
模式六：当模式配置单元1收到由数据处理工控机下发的0x81H指令后，一方面向数据处理工控机发送应答指令0x82H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I1；
模式七：当模式配置单元1收到由数据处理工控机下发的0x91H指令后，一方面向数据处理工控机发送应答指令0x92H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I2；
模式八：当模式配置单元1收到由数据处理工控机下发的0xA1H指令后，一方面向数据处理工控机发送应答指令0xA2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I3；
模式九：当模式配置单元1收到由数据处理工控机下发的0xB1H指令后，一方面向数据处理工控机发送应答指令0xB2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I4；
模式十：当模式配置单元1收到由数据处理工控机下发的0xC1H指令后，一方面向数据处理工控机发送应答指令0xC2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I5；
模式十一：当模式配置单元1收到由数据处理工控机下发的0xD1H指令后，一方面向数据处理工控机发送应答指令0xD2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I6；
模式十二：当模式配置单元1收到由数据处理工控机下发的0xE1H指令后，一方面向数据处理工控机发送应答指令0xE2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I7；
模式十三：当模式配置单元1收到由数据处理工控机下发的0xF1H指令后，一方面向数据处理工控机发送应答指令0xF2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I8；
直射信号处理单元(3)对接收的F3＝{D51-1，D52-1，D53-1}信息进行处理后，分别得到GPS、Galileo、GLONASS和BD的导航信息的运算值D3和载波信号cn(t)；载波解调器(201)对直射信号处理单元(3)输出的载波信号cn(t)根据载波解调关系式dn
Figure BSA00000137769300041
对反射信号F1＝{D51-2，D52-2，D53-2}进行解调，得到载波剥离后的反射信号解调值dn(t)；式中un(t)表示量化后的数字中频信号，cn(t)表示载波信号，t表示序列值，如1，2，3，……，∞，n表示反射通道号；载波信号cn(t)的电平值由{+1，+2，-1，-2}四种组成；
延迟CA码产生器(202)是在延迟映射接收机定位后，第一方面根据DSP处理器赋予的直射处理信息F＝{numPn，SVn，SVPn}通过相位选择关系
Figure BSA00000137769300042
且Mn＝SVPn-numPn，获得延迟CA码的码相位Pn；N表示CA码最大相位值1022；
延迟CA码产生器(202)第二方面依据Pn通过相位运算得到产生CA码的两组十位移位寄存器的初始值reg1和reg2，然后根据reg1和reg2的值并在SVn的调节下，通过改变抽头的方式，产生在Pn相位点后卫星号为SVn的延迟CA码，该延迟CA码用集合表达为
Figure BSA00000137769300043
所述延迟码CA是一个周期为1023的CA码；
延迟CA码产生器(202)中的时钟clk_2为2倍的延迟CA码码速率，时钟clk_20为20倍的延迟CA码码速率；用clk_2的时钟速率读取延迟CA码，产生CA码码间隔为1/2的延迟CA码用于测风模式下相关累计值的输出；用clk_20的时钟速率读取延迟CA码，产生CA码码间隔为1/20的延迟CA码用于测高模式下相关累计值的输出；通过移位寄存器的移位操作，每个反射通道分别产生m组延迟CA码，每组CA码的依次延迟的间隔由累加模式选择确定；
延迟CA码相关运算器(203)根据延迟CA码电平为{+1，-1}的特点，在反射信号相关处理时应用相关逻辑式
Figure BSA00000137769300044
完成反射信号的相关运算；式中anm(t)表示延迟CA码量化值；dn(t)表示载波剥离后的反射信号；znm(t)表示反射信号相关处理后的第n反射通道的第m点的相关值；
1ms相干累加器(204)根据反射信号相关累加式
Figure BSA00000137769300045
得到反射通道m点1ms相关累加值集合
Figure BSA00000137769300046
q11表示第1反射通道中第1个相关累加点，q1m表示第1反射通道中第m个相关累加点，qn1表示第n个反射通道中第1个相关累加点，qnm表示第n个反射通道中第m个相关累加点，同时理可得矩阵其它集合中的元素的物理意义；M为1ms的相关值点的个数。
本发明是一种接收和处理多导航卫星反射信号的相关器的优点在于：
(1)本发明能够同时接收处理四种导航卫星的反射信号，为后续的物理反演提供更加灵活详尽的反演数据。
(2)本发明的模式配置单元能够选择输出原始采集数据和反射信号相关累加值的类型，有效的减少数据处理工控机的存储空间，并为后续的物理反演工作的灵活操作提供实时性保障。
(3)本发明的延迟CA码发生器可在设定的任意相位处连续的产生导航卫星的扩频码，在反射通道相关中省去了大量的移位寄存器，在硬件实现中能过节省大量的逻辑资源。
(4)本发明的载波解调器和延迟CA码相关运算器，可通过逻辑关系式实现普通乘法器功能，在硬件实现中节省大量的逻辑资源，缩短了系统延迟时间。
(5)本发明的1ms的相关累累加器，输出1ms相关点的累加值为后续的物理反演带来更大的灵活性。
附图说明
图1是实时GNSS反射信号接收机的总体结构。
图2是本发明的多导航卫星系统的反射信号接收处理相关器结构图。
图2(A)是本发明反射信号处理单元的结构框图。
具体实施方式
下面将结合附图对本发明做进一步的详细说明。
本发明主要是在现有GPS反射信号处理系统的基础上，通过大容量FPGA芯片，采用VerilogHDL语言编程来实现多导航卫星下反射信号的相关处理和多种模式配置。本发明多卫星导航下的反射信号接收机的相关器包括有载波解调器、延迟CA码发生器和延迟CA相关运算器，采用硬件实现时与普通的相关器相比节省了大量的逻辑资源、缩短了系统的延迟时间。内嵌有本发明相关器的接收机能够为后续物理状态的反演提供更加合理的、灵活的数据源，并为星载装置的研发提供技术基础。FPGA芯片型号为Xilinx公司的XC5VLX220和XC5VFX100T。
如图1所示的实时GNSS反射信号接收机的总体结构，该接收机主要由左旋天线、右旋天线、射频前端、FPGA相关器、数字信号(DSP)处理器、通用PC(数据处理工控机、任务监控工控机)等组成。指向天顶的右旋天线接收直射卫星信号，指向天底的左旋天线接收海面散射信号；即FPGA相关器包含多个并行通道，直射通道与右旋天线相连，直射信号经相关处理后在DSP中通过通常的闭合码和载波环路获得伪距和多普勒观测值，并得到导航解；反射通道与左旋天线相连，根据选定卫星的参数(码相位和载波频率)，DSP配置对应反射通道并且根据直射卫星多普勒频移和码相位信息调整反射通道的配置，从而通过开环方式获得不同时延和多普勒频移的相关功率曲线。
为了适用于多卫星导航下、能够满足延迟映射接收机对多种天线接收的信号进行处理，如图2所示，本发明提出一种通过模式配置来实现对不同天线接收信号的相关值的获取，因此对于相关值获取所采用的FPGA相关器(如图2A所示)包括有模式配置单元1、反射信号处理单元2和直射信号处理单元3；所述反射信号处理单元2由载波解调器201、延迟CA码发生器202、延迟CA码相关运算器203和1ms相干累加器204构成。
在本发明中，借助现有接收机中的DSP处理器向FPGA相关器的延迟CA码发生器202输入的信息有F＝{numPn，SVn，SVPn}，numPn表示卫星反射信号CA码相对于直射信号的码相位延迟数，SVn表示卫星编号，SVPn表示卫星直射信号跟踪后的CA码码相位，n表示反射通道号。
在本发明中，借助现有接收机中的数据处理工控机向FPGA相关器发送模式选择控制指令，同时用于存储FPGA相关器输出的相关累加值
Figure BSA00000137769300061
m表示每一个反射通道的相关累加值的个数，n表示反射通道号。
如图2所示，多种天线(也称天线组)包括第一右旋天线11、第二右旋天线12、第三右旋天线13、第一左旋天线21、第二左旋天线22和第三左旋天线23；第一右旋天线11与第一左旋天线21构成第一天线对；第二右旋天线12与第二左旋天线22构成第二天线对；第三右旋天线13与第三左旋天线23构成第三天线对。在本发明中，第一天线对用于接收GPS、Galileo的直射和反射信号，且工作频率为1570～1580MHz，放大器增益为26dB。第二天线对用于接收GLONASS的直射和反射信号，且工作频率为1560～1620MHz，放大器增益为26dB。第三天线对用于接收BD(北斗)的直射和反射信号，且工作频率为1600～1620MHz，放大器增益为26dB。
频率变换单元5第一方面对第一右旋天线11、第二右旋天线12和第三右旋天线13接收的卫星直射信号进行多级下变频，分别输出模拟中频信号f11、f12和f13；f11表示GPS和Galileo直射信号，f12表示GLONASS直射信号，f13表示BD直射信号；
频率变换单元5第二方面对第一左旋天线21、第二左旋天线22和第三左旋天线23接收的卫星反射信号进行多级下变频，分别输出模拟中频信号f21、f22和f23；f21表示GPS和Galileo反射信号，f22表示GLONASS反射信号，f23表示BD反射信号。
第一模数转换单元51对接收到的f11和f21进行2bit量化后输出数字量化直射信号D51-1和数字量化反射信号D51-2。
第二模数转换单元52对接收到的f12和f22进行2bit量化后输出数字量化直射信号D52-1和数字量化反射信号D52-2。
第三模数转换单元53对接收到的f13和f23进行2bit量化后输出数字量化直射信号D53-1和数字量化反射信号D53-2。
在本发明中，模式配置单元1中存储的模式有：
模式一：当模式配置单元1收到由数据处理工控机下发的0x11H指令后，一方面向数据处理工控机发送应答指令0x12H，另一方面向直射信号处理单元3发出启动指令，记为J1，然后直射信号处理单元3向数据处理工控机输出GPS原始数据。
模式二：当模式配置单元1收到由数据处理工控机下发的0x21H指令后，一方面向数据处理工控机发送应答指令0x22H，另一方面向直射信号处理单元3发出启动指令，记为J2，然后直射信号处理单元3向数据处理工控机输出Galileo原始数据。
模式三：当模式配置单元1收到由数据处理工控机下发的0x31H指令后，一方面向数据处理工控机发送应答指令0x32H，另一方面向直射信号处理单元3发出启动指令，记为J3，然后直射信号处理单元3向数据处理工控机输出GLONASS原始数据。
模式四：当模式配置单元1收到由数据处理工控机下发的0x41H指令后，一方面向数据处理工控机发送应答指令0x42H，另一方面向直射信号处理单元3发出启动指令，记为J4，然后直射信号处理单元3向数据处理工控机输出BD原始数据。
模式五：当模式配置单元1收到由数据处理工控机下发的0x5EH指令后，一方面向数据处理工控机发送应答指令0x5FH，另一方面向直射信号处理单元3发出启动指令，记为J5，然后直射信号处理单元3向数据处理工控机输出GPS、Galileo、GLONASS和BD混合的原始数据。
在本发明中，模式配置单元1向直射信号处理单元3下发的启动指令写成集合表达为D1-3＝{J1，J2，J3，J4，J5}。当模式配置单元1收到除模式一至模式五之外的指令，则延迟映射接收机继续执行上一条有效指令，不改变原有模式。
模式六：当模式配置单元1收到由数据处理工控机下发的0x81H指令后，一方面向数据处理工控机发送应答指令0x82H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I1。
模式七：当模式配置单元1收到由数据处理工控机下发的0x91H指令后，一方面向数据处理工控机发送应答指令0x92H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I2。
模式八：当模式配置单元1收到由数据处理工控机下发的0xA1H指令后，一方面向数据处理工控机发送应答指令0xA2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I3。
模式九：当模式配置单元1收到由数据处理工控机下发的0xB1H指令后，一方面向数据处理工控机发送应答指令0xB2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I4。
模式十：当模式配置单元1收到由数据处理工控机下发的0xC1H指令后，一方面向数据处理工控机发送应答指令0xC2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I5。
模式十一：当模式配置单元1收到由数据处理工控机下发的0xD1H指令后，一方面向数据处理工控机发送应答指令0xD2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I6。
模式十二：当模式配置单元1收到由数据处理工控机下发的0xE1H指令后，一方面向数据处理工控机发送应答指令0xE2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I7。
模式十三：当模式配置单元1收到由数据处理工控机下发的0xF1H指令后，一方面向数据处理工控机发送应答指令0xF2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I8。
在本发明中，模式配置单元1向延迟CA码发生器202下发的相关通道累加值配置指令写成集合表达为D1-2＝{I1，I2，I3，I4，I5，I6，I7，I8}。当模式配置单元1收到除模式六至模式十三之外的指令，则延迟映射接收机继续执行上一条有效指令，不改变原有模式。
在本发明中，模式配置单元1用于实现原始采集数据(如图1所示)的类型选择和相关累加值模式的选择。所述相关累加值是指对反射信号的原始数据进行相关器的处理后得到。但是，在延迟接收机的测风模式下，延迟CA码发生器202中用于计算相关累加值的延迟CA码码间隔为1\2；在延迟接收机的测高模式下，延迟CA码发生器202中用于计算相关累加值的延迟CA码间隔为1\20。
采用本发明设计的FPGA相关器应用在延迟映射接收机上、进行的直射信号和反射信号处理关系为：
频率变换单元5第一方面对第一右旋天线11、第二右旋天线12和第三右旋天线13接收的卫星直射信号进行多级下变频，分别输出模拟中频信号f11、f12和f13；f11表示GPS和Galileo直射信号，f12表示GLONASS直射信号，f13表示BD直射信号；
频率变换单元5第二方面对第一左旋天线21、第二左旋天线22和第三左旋天线23接收的卫星反射信号进行多级下变频，分别输出模拟中频信号f21、f22和f23；f21表示GPS和Galileo反射信号，f22表示GLONASS反射信号，f23表示BD反射信号。
第一模数转换单元51对接收到的f11和f21进行2bit量化后输出数字量化直射信号D51-1和数字量化反射信号D51-2。
第二模数转换单元52对接收到的f12和f22进行2bit量化后输出数字量化直射信号D52-1和数字量化反射信号D52-2。
第三模数转换单元53对接收到的f13和f23进行2bit量化后输出数字量化直射信号D53-1和数字量化反射信号D53-2。
直射信号处理单元3对接收的F3＝{D51-1，D52-1，D53-1}信息采用2006年7月第27卷第4期在《遥测遥控》公开的“GPS接收机基带信号处理模块的FPGA实现”文献进行处理，分别得到GPS、Galileo、GLONASS和B D的导航信息的运算值D3和载波信号cn(t)。
载波解调器201对直射信号处理单元3输出的载波信号cn(t)根据载波解调关系式
Figure BSA00000137769300091
对反射信号F1＝{D51-2，D52-2，D53-2}进行解调，得到载波剥离后的反射信号解调值dn(t)。式中un(t)表示量化后的数字中频信号，cn(t)表示载波信号，t表示序列值，如1，2，3，……，∞，n表示反射通道号。载波信号cn(t)的电平值由{+1，+2，-1，-2}四种组成。
延迟CA码产生器202是在延迟映射接收机定位后，第一方面根据DSP处理器赋予的直射处理信息F＝{numPn，SVn，SVPn}通过相位选择关系
Figure BSA00000137769300101
且Mn＝SVPn-numPn，获得延迟CA码的码相位Pn；N表示CA码最大相位值1022(起始相位值从0开始)。
延迟CA码产生器202第二方面依据Pn通过相位运算得到产生CA码的两组十位移位寄存器的初始值reg1和reg2，然后根据reg1和reg2的值并在SVn的调节下，通过改变抽头的方式，产生在Pn相位点后卫星号为SVn的延迟CA码，该延迟CA码用集合表达为
Figure BSA00000137769300102
所述延迟码CA是一个周期为1023的CA码。
延迟CA码产生器202中的时钟clk_2为2倍的延迟CA码码速率，时钟clk_20为20倍的延迟CA码码速率。用clk_2的时钟速率读取延迟CA码，产生CA码码间隔为1/2的延迟CA码用于测风模式下相关累计值的输出；用clk_20的时钟速率读取延迟CA码，产生CA码码间隔为1/20的延迟CA码用于测高模式下相关累计值的输出。通过移位寄存器的移位操作，每个反射通道分别产生m组延迟CA码，每组CA码的依次延迟的间隔由累加模式选择确定。
延迟CA码相关运算器203根据延迟CA码电平为{+1，-1}的特点，在反射信号相关处理时应用相关逻辑式完成反射信号的相关运算。式中anm(t)表示延迟CA码量化值；dn(t)表示载波剥离后的反射信号；znm(t)表示反射信号相关处理后的第n反射通道的第m点的相关值。1ms相干累加器204根据反射信号相关累加式
Figure BSA00000137769300104
得到反射通道m点1ms相关累加值集合
Figure BSA00000137769300105
q11表示第1反射通道中第1个相关累加点，q1m表示第1反射通道中第m个相关累加点，qn1表示第n个反射通道中第1个相关累加点，qnm表示第n个反射通道中第m个相关累加点，同时理可得矩阵其它集合中的元素的物理意义。M为1ms的相关值点的个数。
本发明的一种多卫星导航下的反射信号接收机的相关器，该相关器是在FPGA芯片，采用VerilogHDL语言编程来实现多导航卫星下反射信号的相关处理和多种模式配置。配置有本发明相关器的反射信号接收机能够为后续物理状态的反演提供更加合理的、灵活的数据源，并为星载装置的研发提供技术基础。

Claims (6)

1.一种多卫星导航下的反射信号接收机的相关器，所述反射信号接收机主要由左旋天线、右旋天线、射频前端、FPGA相关器、DSP处理器、数据处理工控机、任务监控工控机组成；数据处理工控机向FPGA相关器发送模式选择控制指令，同时用于存储FPGA相关器输出的相关累加值m表示每一个反射通道的相关累加值的个数，n表示反射通道号；DSP处理器向FPGA相关器发出卫星码相位信息F＝{numPn，SVn，SVPn}，numPn表示卫星反射信号CA码相对于直射信号的码相位延迟数，SVn表示卫星编号，SVPn表示卫星直射信号跟踪后的CA码码相位，n表示反射通道号；其特征在于：该反射信号接收机的相关器包括有模式配置单元(1)、反射信号处理单元(2)和直射信号处理单元(3)；所述反射信号处理单元(2)由载波解调器(201)、延迟CA码发生器(202)、延迟CA码相关运算器(203)和1ms相干累加器(204)构成；
模式配置单元1中存储的模式有：
模式一：当模式配置单元1收到由数据处理工控机下发的0x11H指令后，一方面向数据处理工控机发送应答指令0x12H，另一方面向直射信号处理单元3发出启动指令，记为J1，然后直射信号处理单元3向数据处理工控机输出GPS原始数据；
模式二：当模式配置单元1收到由数据处理工控机下发的0x21H指令后，一方面向数据处理工控机发送应答指令0x22H，另一方面向直射信号处理单元3发出启动指令，记为J2，然后直射信号处理单元3向数据处理工控机输出Galileo原始数据；
模式三：当模式配置单元1收到由数据处理工控机下发的0x31H指令后，一方面向数据处理工控机发送应答指令0x32H，另一方面向直射信号处理单元3发出启动指令，记为J3，然后直射信号处理单元3向数据处理工控机输出GLONASS原始数据；
模式四：当模式配置单元1收到由数据处理工控机下发的0x41H指令后，一方面向数据处理工控机发送应答指令0x42H，另一方面向直射信号处理单元3发出启动指令，记为J4，然后直射信号处理单元3向数据处理工控机输出BD原始数据；
模式五：当模式配置单元1收到由数据处理工控机下发的0x5EH指令后，一方面向数据处理工控机发送应答指令0x5FH，另一方面向直射信号处理单元3发出启动指令，记为J5，然后直射信号处理单元3向数据处理工控机输出GPS、Galileo、GLONASS和BD混合的原始数据；
模式六：当模式配置单元1收到由数据处理工控机下发的0x81H指令后，一方面向数据处理工控机发送应答指令0x82H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I1；
模式七：当模式配置单元1收到由数据处理工控机下发的0x91H指令后，一方面向数据处理工控机发送应答指令0x92H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I2；
模式八：当模式配置单元1收到由数据处理工控机下发的0xA1H指令后，一方面向数据处理工控机发送应答指令0xA2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I3；
模式九：当模式配置单元1收到由数据处理工控机下发的0xB1H指令后，一方面向数据处理工控机发送应答指令0xB2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I4；
模式十：当模式配置单元1收到由数据处理工控机下发的0xC1H指令后，一方面向数据处理工控机发送应答指令0xC2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I5；
模式十一：当模式配置单元1收到由数据处理工控机下发的0xD1H指令后，一方面向数据处理工控机发送应答指令0xD2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I6；
模式十二：当模式配置单元1收到由数据处理工控机下发的0xE1H指令后，一方面向数据处理工控机发送应答指令0xE2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I7；
模式十三：当模式配置单元1收到由数据处理工控机下发的0xF1H指令后，一方面向数据处理工控机发送应答指令0xF2H，另一方面向延迟CA码发生器202发送相关通道累加值配置指令，记为I8；
直射信号处理单元(3)对接收的F3＝{D51-1，D52-1，D53-1}信息进行处理后，分别得到GPS、Galileo、GLONASS和B D的导航信息的运算值D3和载波信号cn(t)；
载波解调器(201)对直射信号处理单元(3)输出的载波信号cn(t)根据载波解调关系式对反射信号F1＝{D51-2，D52-2，D53-2}进行解调，得到载波剥离后的反射信号解调值dn(t)；式中un(t)表示量化后的数字中频信号，cn(t)表示载波信号，t表示序列值，如1，2，3，……，∞，n表示反射通道号；载波信号cn(t)的电平值由{+1，+2，-1，-2}四种组成；
延迟CA码产生器(202)是在延迟映射接收机定位后，第一方面根据DSP处理器赋予的直射处理信息F＝{numPn，SVn，SVPn}通过相位选择关系
Figure FSA00000137769200031
且Mn＝SVPn-numPn，获得延迟CA码的码相位Pn；N表示CA码最大相位值1022；
延迟CA码产生器(202)第二方面依据Pn通过相位运算得到产生CA码的两组十位移位寄存器的初始值reg1和reg2，然后根据reg1和reg2的值并在SVn的调节下，通过改变抽头的方式，产生在Pn相位点后卫星号为SVn的延迟CA码，该延迟CA码用集合表达为
Figure FSA00000137769200032
所述延迟码CA是一个周期为1023的CA码；
延迟CA码产生器(202)中的时钟clk_2为2倍的延迟CA码码速率，时钟clk_20为20倍的延迟CA码码速率；用clk_2的时钟速率读取延迟CA码，产生CA码码间隔为1/2的延迟CA码用于测风模式下相关累计值的输出；用clk_20的时钟速率读取延迟CA码，产生CA码码间隔为1/20的延迟CA码用于测高模式下相关累计值的输出；通过移位寄存器的移位操作，每个反射通道分别产生m组延迟CA码，每组CA码的依次延迟的间隔由累加模式选择确定；
延迟CA码相关运算器(203)根据延迟CA码电平为{+1，-1}的特点，在反射信号相关处理时应用相关逻辑式
Figure FSA00000137769200033
完成反射信号的相关运算；式中anm(t)表示延迟CA码量化值；dn(t)表示载波剥离后的反射信号；znm(t)表示反射信号相关处理后的第n反射通道的第m点的相关值；
1ms相干累加器(204)根据反射信号相关累加式
Figure FSA00000137769200034
得到反射通道m点1ms相关累加值集合
Figure FSA00000137769200035
q11表示第1反射通道中第1个相关累加点，q1m表示第1反射通道中第m个相关累加点，qn1表示第n个反射通道中第1个相关累加点，qnm表示第n个反射通道中第m个相关累加点，同时理可得矩阵其它集合中的元素的物理意义；M为1ms的相关值点的个数。
2.根据权利要求1所述的多卫星导航下的反射信号接收机的相关器，其特征在于：通过大容量FPGA芯片，采用VerilogHDL语言编程来实现多导航卫星下反射信号的相关处理和多种模式配置。
3.根据权利要求2所述的多卫星导航下的反射信号接收机的相关器，其特征在于：FPGA芯片型号为Xilinx公司的XC5VLX220和XC5VFX100T。
4.根据权利要求1所述的多卫星导航下的反射信号接收机的相关器，其特征在于：所述右旋天线包括有第一右旋天线(11)、第二右旋天线(12)和第三右旋天线(13)；所述左旋天线包括有第一左旋天线(21)、第二左旋天线(22)和第三左旋天线(23)；第一右旋天线(11)与第一左旋天线(21)构成第一天线对；第二右旋天线(12)与第二左旋天线(22)构成第二天线对；第三右旋天线(13)与第三左旋天线(23)构成第三天线对；第一天线对用于接收GPS、Galileo的直射和反射信号，且工作频率为1570～1580MHz，放大器增益为26dB；第二天线对用于接收GLONASS的直射和反射信号，且工作频率为1560～1620MHz，放大器增益为26dB；第三天线对用于接收北斗的直射和反射信号，且工作频率为1600～1620MHz，放大器增益为26dB。
5.根据权利要求1或4所述的多卫星导航下的反射信号接收机的相关器，其特征在于：频率变换单元(5)第一方面对第一右旋天线(11)、第二右旋天线(12)和第三右旋天线(13)接收的卫星直射信号进行多级下变频，分别输出模拟中频信号f11、f12和f13；f11表示GPS和Galileo直射信号，f12表示GLONASS直射信号，f13表示BD直射信号；
频率变换单元(5)第二方面对第一左旋天线(21)、第二左旋天线(22)和第三左旋天线(23)接收的卫星反射信号进行多级下变频，分别输出模拟中频信号f21、f22和f23；f21表示GPS和Galileo反射信号，f22表示GLONASS反射信号，f23表示BD反射信号。
6.根据权利要求1所述的多卫星导航下的反射信号接收机的相关器，其特征在于：模式配置单元(1)用于实现原始采集数据的类型选择和相关累加值模式的选择；所述相关累加值是指对反射信号的原始数据进行相关器的处理后得到；但是，在延迟接收机的测风模式下，延迟CA码发生器(202)中用于计算相关累加值的延迟CA码码间隔为1\2；在延迟接收机的测高模式下，延迟CA码发生器(202)中用于计算相关累加值的延迟CA码间隔为1\20。
CN201010193964A 2010-06-08 2010-06-08 多卫星导航下的反射信号接收机的相关器 Expired - Fee Related CN101866008B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010193964A CN101866008B (zh) 2010-06-08 2010-06-08 多卫星导航下的反射信号接收机的相关器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010193964A CN101866008B (zh) 2010-06-08 2010-06-08 多卫星导航下的反射信号接收机的相关器

Publications (2)

Publication Number Publication Date
CN101866008A true CN101866008A (zh) 2010-10-20
CN101866008B CN101866008B (zh) 2012-10-10

Family

ID=42957809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010193964A Expired - Fee Related CN101866008B (zh) 2010-06-08 2010-06-08 多卫星导航下的反射信号接收机的相关器

Country Status (1)

Country Link
CN (1) CN101866008B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221700A (zh) * 2011-03-30 2011-10-19 无锡北斗卫导科技有限公司 一种利用反射信号增强直射信号的方法
CN106443734A (zh) * 2016-09-18 2017-02-22 广州知春里网络科技有限公司 一种在arm平台实现厘米级精度gps定位的系统
CN107015063A (zh) * 2017-03-21 2017-08-04 中国科学院国家天文台 宽带多信道数字相关接收机及接收方法
CN111045038A (zh) * 2019-12-27 2020-04-21 中国科学院国家空间科学中心 一种gnss-r接收机反射通道数量的优化方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003012473A1 (en) * 2001-08-01 2003-02-13 Roke Manor Research Limited Passive moving object detection system and method using signals transmitted by a mobile telephone station
CN101344589A (zh) * 2008-08-29 2009-01-14 北京航空航天大学 基于gnss反射信号的空间飞行器探测装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003012473A1 (en) * 2001-08-01 2003-02-13 Roke Manor Research Limited Passive moving object detection system and method using signals transmitted by a mobile telephone station
CN101344589A (zh) * 2008-08-29 2009-01-14 北京航空航天大学 基于gnss反射信号的空间飞行器探测装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《北京航空航天大学学报》 20090930 杨东凯等 软件定义的GNSS反射信号接收机设计 1048-1051 1-6 第35卷, 第9期 *
《航空学报》 20090131 杨东凯等 基于导航卫星反射信号的空间飞行器探测定位方法 143-147 1-6 第30卷, 第1期 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221700A (zh) * 2011-03-30 2011-10-19 无锡北斗卫导科技有限公司 一种利用反射信号增强直射信号的方法
CN106443734A (zh) * 2016-09-18 2017-02-22 广州知春里网络科技有限公司 一种在arm平台实现厘米级精度gps定位的系统
CN107015063A (zh) * 2017-03-21 2017-08-04 中国科学院国家天文台 宽带多信道数字相关接收机及接收方法
CN111045038A (zh) * 2019-12-27 2020-04-21 中国科学院国家空间科学中心 一种gnss-r接收机反射通道数量的优化方法

Also Published As

Publication number Publication date
CN101866008B (zh) 2012-10-10

Similar Documents

Publication Publication Date Title
Jumaah et al. Technique for order performance by similarity to ideal solution for solving complex situations in multi-criteria optimization of the tracking channels of GPS baseband telecommunication receivers
CN101246212B (zh) 一种可接收处理导航卫星反射信号的相关装置
CN100594390C (zh) 全球定位系统讯号的抑制方法及其全球定位系统接收器
CN1543714B (zh) 码相位搜索扩频信号的方法
CN102109604B (zh) Gps/galileo导航基带处理芯片及导航接收机
CN101802634B (zh) Gnss信号处理器
CN100360953C (zh) 带有校正多径信号效应装置的射频信号接收机及启动该接收机的方法
CN101866008B (zh) 多卫星导航下的反射信号接收机的相关器
CN1737603A (zh) 确定射频信号接收机位置的方法及实现该方法的接收机
CN1109413C (zh) 在gps接收机中消除多普勒复制谐波的方法
CN105319569A (zh) 一种基于多核dsp的惯性/卫星深组合信息处理硬件平台
Hurskainen et al. Multicore software-defined radio architecture for GNSS receiver signal processing
WO2007068661A1 (en) Signal processor and signal processing method
CN201716419U (zh) 机载多普勒/延迟映射接收机
CN117269997A (zh) 基带信号处理电路、接收机及基带信号处理方法
Ucar et al. A chip scale atomic clock driven receiver for multi-constellation GNSS
Guruprasad FPGA-based GNSS receiver design for reflectometry applications
CN100438360C (zh) 处理扩频信号的方法、设备和接收扩频信号的接收机
CN105068095B (zh) 一种提高卫星伪距精度的跟踪系统与方法
CN101571584B (zh) 一种可配置的反射信号功率处理系统
CN202102114U (zh) Gps和glonass多模卫星导航(gnss)接收机的快速搜索引擎
Raasakka et al. GNSS baseband processing in a multi-core platform
CN103914625A (zh) 基于cordic算法的无线信号测向方法
Miralles et al. FFT Algorithms Implementation for Real Time GNSS Receivers in Embedded Processors
CN102096081B (zh) 一种载波与伪随机码剥离电路

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121010

Termination date: 20150608

EXPY Termination of patent right or utility model