发明名称 基于四维UKF的高动态GNSS载波的开环补偿跟踪方法
摘要 基于四维UKF的高动态GNSS载波的开环补偿跟踪方法,它涉及一种全球导航卫星系统载波的跟踪方法。它解决了现有的GNSS载波跟踪方法跟踪精度低、环路跟踪稳定性差、受带宽限制严重而无法适应较大范围的多普勒频移的问题。其方法:捕获数字中频信号的码相位和频率偏移,将频率偏移给UKF滤波器赋初值,将码相位值给本地码跟踪环赋初值并进行本地码跟踪;两路本地载波信号与中频信号相乘后再与本地码跟踪结果进行相关积分,获得两路观测信号后送入UKF滤波器;调整滤波器参数,对高动态的载波信号过程进行跟踪;进行输出补偿,数据解调后进行下一时刻的载波跟踪和数据解调。本发明适用于高动态条件下的信号跟踪过程。
申请公布号 CN101581775A 申请公布日期 2009.11.18
申请号 CN200910072317.5 申请日期 2009.06.18
申请人 哈尔滨工业大学 发明人 孟维晓;陈曦;王文静;韩帅
分类号 G01S5/02(2006.01)I;G01S1/02(2006.01)I 主分类号 G01S5/02(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 代理人 张宏威
主权项 1、基于四维UKF的高动态GNSS载波的开环补偿跟踪方法,其特征是:它由以下步骤完成:步骤一:采用1ms的数据对数字中频信号进行码相位捕获,获得码相位值,采用10ms的数据对所述数字中频信号进行频率捕获,获得频率偏移值,将所述频率偏移值给UKF滤波器赋初值;同时将所得码相位值给本地码跟踪环赋初值并进行本地码跟踪,获得本地码跟踪结果;步骤二:将频率为fc的本地载波信号分别同时经SIN映射和COS映射后获得两路信号,所述两路信号分别与数字中频信号相乘获得信号I和信号Q,所述信号I和信号Q的表达式分别为:<maths id="math0001" num="0001" ><math><![CDATA[ <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mn>2</mn> <mi>P</mi> </msqrt> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>C</mi> <mi>PRN</mi> </msub> <mo>[</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&xi;</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mi>t</mi> <mo>-</mo> <mi>&tau;</mi> <mo>]</mo> <mo>[</mo> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mi>d</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>&Phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>w</mi> </mrow> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>w</mi> <mi>d</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>&Phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>]]></math></maths><maths id="math0002" num="0002" ><math><![CDATA[ <mrow> <mi>Q</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mn>2</mn> <mi>P</mi> </msqrt> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>C</mi> <mi>PRN</mi> </msub> <mo>[</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&xi;</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mi>t</mi> <mo>-</mo> <mi>&tau;</mi> <mo>]</mo> <mo>[</mo> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mi>d</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>&Phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>w</mi> </mrow> <mi>c</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>w</mi> <mi>d</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>&Phi;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>]]></math></maths>所述P为本地载波信号的接收功率,D(t)为周期为20ms的导航电文,CPRN(t)为扩频码,Rb为伪码速率,ξ为高动态引起的码速率偏移,τ为接收扩频码的相位延迟,wc为本地载波频率,wd为本地载波信号的多普勒频移,Φ0为本地载波信号的初始相位,N(t)为高斯白噪声;步骤三:将步骤二获得的信号I和信号Q分别与步骤一中获得的本地码跟踪结果相乘,并对相乘后的结果进行积分清洗,累积积分时间至一个码周期后分别获得信号Ip和信号Qp,所述信号Ip和信号Qp的表达式为:Ip(t)=AD(t)[sin(wdt+Φ0)]+Np(t)Qp(t)=AD(t)[cos(wdt+Φ0)]+Np(t)式中A为信号幅度,Np(t)为经积分清洗后的加性高斯噪声;所述本地码跟踪为同步跟踪;步骤四:将步骤三获得的信号Ip和信号Qp进行二倍相位转换,获得滤波器观测信号IY和QY,所述信号IY和QY的表达式分别为:IY=A2D(t)2sin2(wdt+Φ0)+NYQY=A2D(t)2cos2(wdt+Φ0)+NY步骤五:调整UKF滤波器的参数,对高动态的载波信号过程进行跟踪,所述参数为四个状态变量:载波相位θ、多普勒频移wd、多普勒频移一阶变化率wd1和多普勒频移二阶变化率wd2;步骤六:采用补偿器对UKF滤波器的四个状态变量进行输出补偿,获得补偿后的结果θk,并将补偿后的结果θk代入公式:Fk=Ip,ksin(θk)+Qp.kcos(θk)=D(tk)cos(Δθ)进行数据解调,获得解调结果Dk。
地址 150001黑龙江省哈尔滨市南岗区西大直街92号