发明名称 多卫星导航下的反射信号接收机的相关器
摘要 本发明公开了一种多卫星导航下的反射信号接收机的相关器,该反射信号接收机的相关器包括有模式配置单元(1)、反射信号处理单元(2)和直射信号处理单元(3);所述反射信号处理单元(2)由载波解调器(201)、延迟CA码发生器(202)、延迟CA码相关运算器(203)和1ms相干累加器(204)构成。本发明的相关器是在FPGA芯片,采用VerilogHDL语言编程来实现多导航卫星下反射信号的相关处理和多种模式配置。配置有本发明相关器的反射信号接收机能够为后续物理状态的反演提供更加合理的、灵活的数据源,并为星载装置的研发提供技术基础。
申请公布号 CN101866008B 申请公布日期 2012.10.10
申请号 CN201010193964.4 申请日期 2010.06.08
申请人 北京航空航天大学 发明人 杨东凯;吴红甲;张凤元;李明里
分类号 G01S19/13(2010.01)I;G01S19/35(2010.01)I 主分类号 G01S19/13(2010.01)I
代理机构 北京永创新实专利事务所 11121 代理人 李有浩
主权项 1.一种多卫星导航下的反射信号接收机的相关器,所述反射信号接收机主要由左旋天线、右旋天线、射频前端、FPGA相关器、DSP处理器、数据处理工控机、任务监控工控机组成;数据处理工控机向FPGA相关器发送模式选择控制指令,同时用于存储FPGA相关器输出的相关累加值<maths num="0001"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>q</mi><mn>11</mn></msub></mtd><mtd><msub><mi>q</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>21</mn></msub></mtd><mtd><msub><mi>q</mi><mn>22</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>q</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>q</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mi>nm</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>m表示每一个反射通道的相关累加值的个数,n表示反射通道号;DSP处理器向FPGA相关器发出卫星码相位信息F={numP<sub>n</sub>,SV<sub>n</sub>,SVP<sub>n</sub>},numP<sub>n</sub>表示卫星反射信号CA码相对于直射信号的码相位延迟数,SV<sub>n</sub>表示卫星编号,SVP<sub>n</sub>表示卫星直射信号跟踪后的CA码码相位,n表示反射通道号;其特征在于:该反射信号接收机的相关器包括有模式配置单元(1)、反射信号处理单元(2)和直射信号处理单元(3);所述反射信号处理单元(2)由载波解调器(201)、延迟CA码发生器(202)、延迟CA码相关运算器(203)和1ms相干累加器(204)构成;模式配置单元(1)中存储的模式有:模式一:当模式配置单元(1)收到由数据处理工控机下发的0x11H指令后,一方面向数据处理工控机发送应答指令0x12H,另一方面向直射信号处理单元(3)发出启动指令,记为J<sub>1</sub>,然后直射信号处理单元(3)向数据处理工控机输出GPS原始数据;模式二:当模式配置单元(1)收到由数据处理工控机下发的0x21H指令后,一方面向数据处理工控机发送应答指令0x22H,另一方面向直射信号处理单元(3)发出启动指令,记为J<sub>2</sub>,然后直射信号处理单元(3)向数据处理工控机输出Galileo原始数据;模式三:当模式配置单元(1)收到由数据处理工控机下发的0x31H指令后,一方面向数据处理工控机发送应答指令0x32H,另一方面向直射信号处理单元(3)发出启动指令,记为J<sub>3</sub>,然后直射信号处理单元(3)向数据处理工控机输出GLONASS原始数据;模式四:当模式配置单元(1)收到由数据处理工控机下发的0x41H指令后,一方面向数据处理工控机发送应答指令0x42H,另一方面向直射信号处理单元(3)发出启动指令,记为J<sub>4</sub>,然后直射信号处理单元(3)向数据处理工控机输出BD原始数据;模式五:当模式配置单元(1)收到由数据处理工控机下发的0x5EH指令后,一方面向数据处理工控机发送应答指令0x5FH,另一方面向直射信号处理单元(3)发出启动指令,记为J<sub>5</sub>,然后直射信号处理单元(3)向数据处理工控机输出GPS、Galileo、GLONASS和BD混合的原始数据;模式六:当模式配置单元(1)收到由数据处理工控机下发的0x81H指令后,一方面向数据处理工控机发送应答指令0x82H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>1</sub>;模式七:当模式配置单元(1)收到由数据处理工控机下发的0x91H指令后,一方面向数据处理工控机发送应答指令0x92H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>2</sub>;模式八:当模式配置单元(1)收到由数据处理工控机下发的0xA1H指令后,一方面向数据处理工控机发送应答指令0xA2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>3</sub>;模式九:当模式配置单元(1)收到由数据处理工控机下发的0xB1H指令后,一方面向数据处理工控机发送应答指令0xB2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>4</sub>;模式十:当模式配置单元(1)收到由数据处理工控机下发的0xC1H指令后,一方面向数据处理工控机发送应答指令0xC2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>5</sub>;模式十一:当模式配置单元(1)收到由数据处理工控机下发的0xD1H指令后,一方面向数据处理工控机发送应答指令0xD2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>6</sub>;模式十二:当模式配置单元(1)收到由数据处理工控机下发的0xE1H指令后,一方面向数据处理工控机发送应答指令0xE2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>7</sub>;模式十三:当模式配置单元(1)收到由数据处理工控机下发的0xF1H指令后,一方面向数据处理工控机发送应答指令0xF2H,另一方面向延迟CA码发生器(202)发送相关通道累加值配置指令,记为I<sub>8</sub>;直射信号处理单元(3)对接收的F<sub>3</sub>={D<sub>51-1</sub>,D<sub>52-1</sub>,D<sub>53-1</sub>}信息进行处理后,分别得到GPS、Galileo、GLONASS和B D的导航信息的运算值D<sub>3</sub>和载波信号c<sub>n</sub>(t);D<sub>51-1</sub>表示第一模数转换单元(51)对接收到的f<sub>11</sub>进行2bit量化后输出的数字量化直射信号,D<sub>52-1</sub>表示第二模数转换单元(52)对接收到的f<sub>12</sub>进行2bit量化后输出的数字量化直射信号,D<sub>53-1</sub>表示第三模数转换单元(53)对接收到的f<sub>13</sub>进行2bit量化后输出的数字量化直射信号;f<sub>11</sub>表示GPS和Galileo直射信号,f<sub>12</sub>表示GLONASS直射信号,f<sub>13</sub>表示BD直射信号;载波解调器(201)对直射信号处理单元(3)输出的载波信号c<sub>n</sub>(t)根据载波解调关系式<maths num="0002"><![CDATA[<math><mrow><msub><mi>d</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>u</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>c</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>2</mn><msub><mi>u</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>c</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mn>2</mn></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>u</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>c</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>2</mn><msub><mi>u</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>c</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mn>2</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>对反射信号F<sub>1</sub>={D<sub>51-2</sub>,D<sub>52-2</sub>,D<sub>53-2</sub>}进行解调,得到载波剥离后的反射信号解调值d<sub>n</sub>(t);式中u<sub>n</sub>(t)表示量化后的数字中频信号,c<sub>n</sub>(t)表示载波信号,t表示序列值,n表示反射通道号;载波信号c<sub>n</sub>(t)的电平值由{+1,+2,-1,-2}四种组成;D<sub>51-2</sub>表示第一模数转换单元(51)对接收到的f<sub>21</sub>进行2bit量化后输出的数字量化反射信号D<sub>51-2</sub>,D<sub>52-2</sub>表示第二模数转换单元(52)对接收到的f<sub>22</sub>进行2bit量化后输出的数字量化反射信号D<sub>52-2</sub>,D<sub>53-2</sub>表示第三模数转换单元(53)对接收到的f<sub>23</sub>进行2bit量化后输出的数字量化反射信号;f<sub>21</sub>表示GPS和Galileo反射信号,f<sub>22</sub>表示GLONASS反射信号,f<sub>23</sub>表示BD反射信号;延迟CA码发生器(202)是在反射信号接收机定位后,第一方面根据DSP处理器赋予的直射处理信息F={numP<sub>n</sub>,SV<sub>n</sub>,SVP<sub>n</sub>}通过相位选择关系<maths num="0003"><![CDATA[<math><mrow><msub><mi>P</mi><mi>n</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>M</mi><mi>n</mi></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>M</mi><mi>n</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>N</mi><mo>-</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>M</mi><mi>n</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><msub><mi>M</mi><mi>n</mi></msub><mo>&lt;</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>且M<sub>n</sub>=SVP<sub>n</sub>-numP<sub>n</sub>,获得延迟CA码的码相位P<sub>n</sub>;N表示CA码最大相位值1022;延迟CA码发生器(202)第二方面依据延迟CA码的码相位P<sub>n</sub>通过相位运算得到产生CA码的两组十位移位寄存器的初始值reg1和reg2,然后根据reg1和reg2的值并在SV<sub>n</sub>的调节下,通过改变抽头的方式,产生在所述延迟CA码的码相位P<sub>n</sub>的相位点之后的卫星号为SV<sub>n</sub>的延迟CA码,该延迟CA码用集合表达为<img file="FSB00000799234100032.GIF" wi="643" he="80" />所述延迟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}的特点,在反射信号相关处理时应用相关逻辑式<maths num="0004"><![CDATA[<math><mrow><msub><mi>z</mi><mi>nm</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mo>-</mo><msub><mi>d</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>a</mi><mi>nm</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>d</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>a</mi><mi>nm</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>完成反射信号的相关运算;式中a<sub>nm</sub>(t)表示延迟CA码量化值;d<sub>n</sub>(t)表示载波剥离后的反射信号解调值;z<sub>nm</sub>(t)表示反射信号相关处理后的第n反射通道的第m点的相关值;1ms相干累加器(204)根据反射信号相关累加式<img file="FSB00000799234100034.GIF" wi="290" he="122" />得到反射通道m点1ms相关累加值集合<maths num="0005"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>q</mi><mn>11</mn></msub></mtd><mtd><msub><mi>q</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>21</mn></msub></mtd><mtd><msub><mi>q</mi><mn>22</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>q</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>q</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>q</mi><mi>nm</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>q<sub>11</sub>表示第1反射通道中第1个相关累加点,q<sub>1m</sub>表示第1反射通道中第m个相关累加点,q<sub>n1</sub>表示第n个反射通道中第1个相关累加点,q<sub>nm</sub>表示第n个反射通道中第m个相关累加点,同理可得矩阵其它集合中的元素的物理意义;M为1ms的相关值点的个数。
地址 100191 北京市海淀区学院路37号