专利名称：基于四维ukf的高动态gnss载波的开环补偿跟踪方法
技术领域：
本发明涉及一种全球导航卫星系统载波的跟踪方法。
技术背景全球定位导航卫星系统在陆地、海洋、航空和航天等领域都有着广泛的应 用。GNSS (全球导航卫星系统)接收机的信号层处理部分主要由捕获和跟踪 模块组成。GNSS接收机在完成导航信号的捕获后，便转入跟踪流程，在捕获 所得频偏的基础上进一步准确跟踪信号的相位、频率变化以实现导航电文的正 确解调处理。捕获只在接收机开机阶段工作，而信号跟踪要贯穿整个接收过程。 后者性能的好坏直接决定着接收机能否实现准确定位。GNSS信号的跟踪就是要准确估计出信号的频率、相位，利用本地NCO 实现载波剥离，获得导航电文。传统的跟踪方法大多采用硬件跟踪技术，能够 适应中低动态环境。采用最为广泛的是Costas环路结构，其中的鉴相器和环路 滤波器是影响跟踪能力的关键。通常采用的算法为二阶PLL环路滤波法，其原 理简单，硬件实现也很容易。但其所能跟踪的频率变化只为十几赫兹，如果频 率变化较大则需要FLL跟踪环路作辅助。虽然FLL可在一定范围内提高环路 的跟踪能力，但是这种能力的提高需要增大环'路带宽，由此会引进更多的噪声， 导致PLL跟踪精度下降。因此FLL辅助PLL会受到环路带宽的限制。由于近 年来DSP和嵌入式处理器的运行速度大大提高，使得在软件接收机上方便灵 活的修改和设计新的算法成为可能。本发明即是针对软件接收机设计的一种载 波信号跟踪方法。导航定位技术在高动态领域的应用越来越多，例如对高速运动的飞机或 其它机动目标进行定位跟踪。但是与静态、中低动态环境相比，高动态环境使 GNSS信号产生了较大的多普勒频移，并且伴随着较高的频率变化率，传统的 锁相锁频环由于受到环路带宽的限制难以实现高动态信号的跟踪。为了能够使 接收机在高动态环境下保持正常工作，必须改进已有算法或采用新的算法。现代GNSS系统中的GPS导航定位系统应用最为广泛，Galileo导航系统和我国BD导航系统采用的信号也均采用扩频调制，在信号跟踪部分所采用的 方法流程一致，仅仅是参数不同。目前对于高动态条件下的信号跟踪技术已经有了一些进展，主要在于以下 几个方面。 一种是提高PLL的阶数，采用三阶环路滤波器，构成四阶锁相环 路。此方法参数设计较为复杂，其辅助端FLL的频偏跟踪范围为i0.5iOfe，如 果Av^在调整过程中超过0.5kHz，环路就会失锁，受到这一限制，锁相锁频环 无法实现高动态信号的跟踪。另一种方法是寻求其他算法代替FLL的作用， 比较有代表性的是采用扩展卡尔曼滤波算法(EKF)去辅助PLL。 EKF虽然能 够跟踪高动态信号，估计出载波信号的频偏，但是其受初值影响较大，如果选 择不好会造成滤波器发散。而且EKF的高动态信号跟踪能力受到调整时间的 限制，调整时间间隔越长跟踪效果越差。由于扩频码周期限定了调整时间间隔 为lms，其跟踪能力表现得很不稳定，会严重影响数据解调。另外还有人尝试 将UKF使用到闭环当中，取代PLL的作用。由于UKF跟踪的频率范围没有 限制，所以不需要其它辅助。但是研究中发现由于闭环结构的限制使得UKF 只能采用三维模型表示相位误差信号，跟踪精度低。
发明内容
' 本发明是为了解决现有的GNSS载波跟踪方法跟踪精度低和环路跟踪稳 定性差、受带宽限制严重、无法适应较大范围的多普勒频移的问题，提出基于 四维UKF的高动态GNSS载波的开环补偿跟踪方法。四维UKF的高动态GNSS载波的开环补偿跟踪方法，它由以下步骤完成 步骤一采用lms的数据对数字中频信号进行码相位捕获，获得码相位 值，釆用10ms的数据对所述数字中频信号进行频率捕获，获得频率偏移值， 将所述频率偏移值给UKF滤波器赋初值；同时将所得码相位值给本地码跟踪 环赋初值并进行本地码跟踪，获得本地码跟踪结果；步骤二将频率为/。的本地载波信号分别同时经SIN映射和COS映射后获得两路信号，所述两路信号分别与数字中频信号相乘获得信号/和信号2，所述信号/和信号2的表达式分别为 '/(0 = ^D(/)C觸[(1 +脾-r][si勿+ 00) + sin(2wc, ++ 00)]+聰 ,=^Z)(OC』+柳-r][cos(w/ + 。o) + cos(2Wcf ++ 0>0)] + ，) 所述尸为本地载波信号的接收功率，D(0为周期为20ms的导航电文，C層(/) 为扩频码，^为伪码速率，《为高动态引起的码速率偏移，"T为接收扩频码的 相位延迟，Wc为本地载波频率，、&为本地载波信号的多普勒频移，0。为本地载波信号的初始相位，iV(/)为高斯白噪声；步骤三将步骤二获得的信号/和信号。分别与步骤一中获得的本地码跟踪结果相乘，并对相乘后的结果进行积分清洗，累积积分时间至一个码周期后分别获得信号/p和信号^ ，所述信号/p和信号^的表达式为=相(,)[si勿+ (Do)] + ^(,)& (,)=卿)[co—/ + <D0)] + (/)式中A为信号幅度， (f)为经积分清洗后的加性高斯噪声； 所述本地码跟踪为同步跟踪；步骤四将步骤三获得的信号/p和信号&进行二倍相位转换，获得滤波 器观测信号^和& ，所述信号八和込的表达式分别为/y =^D")2 sin2(>y + <D0) + A^= J2Z)0)2 COS 2(W/ + 0。) + TVy步骤五调整UKF滤波器的参数，对高动态的载波信号过程进行跟踪， 所述参数为四个状态变量载波相位e、多普勒频移^、多普勒频移一阶变化 率^和多普勒频移二阶变化率w〗； ，步骤六采用补偿器对UKF滤波器的四个状态变量进行输出补偿，获得 补偿后的结果《，并将补偿后的结果《代入公式A = /p,fc si,) + gp fc co赋)=) COS,进行数据解调，获得解调结果A。步骤五中调整UKF滤波器参数的方法是UKF滤波器的状态变量为载 波相位P、多普勒频移v^、多普勒频移一阶变化率^和二阶变化率^，调整 时间间隔&为lms，当前时刻为/fc， UKF滤波器各状态变量与第/k-l时刻的关 系为<formula>formula see original document page 7</formula>式中W为系统噪声;观测向量<formula>formula see original document page 7</formula>式:与状态向量的关系为<formula>formula see original document page 7</formula>
式中R为测量噪声；步骤六所述补偿器对UKF滤波器的四个状态变量的输出补偿是根据公
<formula>formula see original document page 7</formula>
进行的。有益效果:本发明能够实现高动态载波信号的跟踪;本发明采用四维UKF 实现高动态环境下的载波跟踪，不仅跟踪信号频率，同时跟踪载波相位；本发 明不会受到类似FLL和PLL所用带宽的限制，能够适应较大范围的多普勒频 移。而且跟踪精度和环路跟踪稳定性能也明显优于采用闭环形式的EKF辅助 PLL的方法，环路跟踪稳定性好。与采用闭环结构相比，本发明使用的UKF 使用四维模型将结构上的开环引入到滤波器内部，通过UKF自身的调整代替 外部反馈，提高了跟踪精度。
图l是本发明的方法的信号处理流程示意图；图2是本发明的UKF相位 跟踪能力效果图，其中横坐标为时间(ms)，纵坐标为相位(rad);图3是本 发明的UKF频率跟踪能力效果图，其中横坐标为时间(ms)，纵坐标为频率 (他)；图4是本发明的UKF频率一阶变化率跟踪能力效果图，其中横坐标为时间(ms)，纵坐标为频率一阶变化率(论")；图5是本发明UKF频率二阶 变化率跟踪能力效果图，其中横坐标为时间(ms)，纵坐标为频率二阶变化率(历"2);图6是本发明的相位跟踪误差曲线图，其中横坐标为时间(ms)， 纵坐标为相位误差(rad);图7是频率跟踪误差曲线图，其中横坐标为时间(ms)， 纵坐标为频率误差(他)；图8是频率一阶突化率跟踪误差曲线图，其中横坐 标为时间(ms)，纵坐标为频率一阶变化率误差(他")；图9是频率二阶变化 率跟踪误差曲线图，其中横坐标为时间(ms)，纵坐标为二阶变化率误差(Hz"2);图IO为EKF辅助PLL、 UKF闭环跟踪与本发明的效果对比图，其 中横坐标为加加速度(g/s)，纵坐标为均方根误差(论)，曲线21为EKF辅 助PLL的相位跟踪曲线，曲线22为UKF闭环跟踪曲线，曲线23为本发明的 跟踪曲线；图11为EKF辅助PLL、 UKF闭环跟踪与UKF开环跟踪动态频率 跟踪效果对比图，其中横坐标为加加速度(g/s),纵坐标为均方根误差(/fe)， 图中曲线31为EKF辅助PLL频率跟踪曲线，曲线32为UKF闭环频率跟踪曲 线，曲线33为本发明的频率跟踪曲线。 '具体实施方式
具体实施方式
一结合图1说明本具体实施方式
，基于四维UKF的高动 态GNSS载波的开环补偿跟踪方法，它由以下步骤完成步骤一采用lms的数据对数字中频信号进行码相位捕获，获得码相位 值，采用10ms的数据对所述数字中频信号进行频率捕获，获得频率偏移值， 将所述频率偏移值给UKF滤波器赋初值；同时将所得码相位值给本地码跟踪 环赋初值并进行本地码跟踪，获得本地码跟踪结果；本步骤采用先粗捕再精捕的方法提高捕获精度，使UKF滤波器的初始值 与初始频偏更加接近，加快滤波器的收敛速度。捕获所得码偏给本地码跟踪环 进行码跟踪，复现本地伪码。步骤二将频率为/。的本地载波信号分别同时经SIN映射和COS映射后 获得两路信号，所述两路信号分别与数字中频信号相乘获得信号/和信号2， 所述信号/和信号g的表达式分别为/(,) = V^ZP(f )C顺[(l +柳-"[si勿+ 4>0) + sin(2w/ + >V + Oo)] + , 2") = #Z)(0C』1 +柳-r][cos(w/ + O0)+cos(2Wcf + w乂 + 00)] +竭所述尸为本地载波信号的接收功率，D(0为周期为20ms的导航电文，C，(O 为扩频码，A为伪码速率，《为高动态引起的码速率偏移，r为接收扩频码 的相位延迟，We为本地载波频率，V^为本地载波信号的多普勒频移，。。为本地载波信号的初始相位，A^)为高斯白噪声；步骤三将步骤二获得的信号/和信号^分别与步骤一中获得的本地码跟踪结果相乘，并对相乘后的结果进行积分清洗，累积积分时间至一个码周期后分别获得信号/p和信号& ，所述信号/P和信号&的表达式为;(f)=麟)[si勿+ 。。)〗+ #p (/)& W =鄉)[cos(iv/ + O0)〗+ (,)式中A为信号幅度， W为经积分清洗后的加性高斯噪声； 所述本地码环为同步跟踪； 本步骤通过解扩提高信噪比；步骤四将步骤三获得的信号/p和信号^进行二倍相位转换，获得滤波 器观测信号和&，所述信号/y和&的表达式分别为A = fiXO2 sin 2(wJ + d>0) + W 込=j cos 2(wJ +①o) + ^本步骤通过步骤三累积积分时间至一个码周期(即1ms)后，起到滤除二 倍频信号的作用；并进行UKF滤波器前端观测信号的预处理；本步骤的目的 是因为步骤三获得的/,和A两路信号不能直接送入滤波器，因为这两路信号 还属于BPSK调制信号，直接送入UKF滤波器会由于存在180°的相位翻转 使得UKF滤波器无法正确识别载波相位，导致环 路失锁，因此在信号进入UKF 滤波器之间做如下处理 込=《—《=^D(02[cos20v + O0)-sin2(>V + a)。)] +ATy 使/>(02=1，这样就可以消除BPSK调制的影响，使载波跟踪环路正常工作。 步骤五调整UKF滤波器的参数，对高动态的载波信号过程进行跟踪，所述参数为四个状态变量载波相位0、多普勒频移^、多普勒频移一阶变化率^和多普勒频移二阶变化率<;步骤六采用补偿器对UKF滤波器的四个状态变量进行输出补偿，获得补偿后的结果《，并将补偿后的结果《代入公式= sin(《)+仏* C0S(A ) = ) C0S(Ae) 进行数据解调，获得解调结果A后，进行下一时刻的载波跟踪和数据解调；步骤六所述UKF滤波器的四个状态变量分别为信号的载波相位^、多普 勒频移W、多普勒频移一阶变化率乂和二阶变化率^。当A^足够小时，cos(A0大于O且接近l。由于导航电文采用BPSK调制，Z)("为土l，那么设置判决门限为O,就可以进行数据的正确解调。步骤五中调整UKF滤波器参数的方法是UKF滤波器的状态变量为载波相位0、多普勒频移W、多普勒频移一阶变化率"和二阶变化率w〗，调整时间间隔&为lms，当前时刻为h UKF滤波器各状态变量与第;fc-l时刻的关系为 '<formula>formula see original document page 10</formula>式中W为系统噪声；观测向量<formula>formula see original document page 10</formula>及与状态向量的关系为sin 2(《+ 力/2 + i*2/6 + w么<ir3/24) cos 2(《+ + A2/6 +《fc &3/24)_式中R为测量噪声；将这两个方程应用到UKF滤波器中，适当调整滤波器参数，进行预测和 更新，就可以对高动态的载波信号过程进行跟踪。UKF滤波器不仅可以给出 解调所需的当前载波相位，还可以给出当前载波的多普勒频移以及其一阶变化 率和二阶变化率。但由于UKF的滤波估计具有一定的滞后性，直接利用估计出来的相位进行数据解调误差较大，因此需要输出信号补偿。步骤六所述补偿器对UKF滤波器的四个状态变量的输出补偿是根据公式<formula>formula see original document page 11</formula>本发明所述四维就是利用信号的相位、，频率、频率一阶变化率和二阶变 化率四个参量来表示载波相位变化信息。 一般的中低动态信号仅采用前面三个参量，原因是EKF在保证理论没有误差的条件下只能到三维，而其他采用UKF 四维的却没有解决数据解调问题，而且根据它们建立的模型不准确，跟踪精度 也没有本发明采用补偿方法的精度高。所以四维既充分表现了载波信号的动态 性，也在UKF的合理设计之内，还能够充分利用UKJF的非线性跟踪能力。本发明采用的环路不再使用PLL，在很大程度上减少了参数设计所带来 的困难。同时用二倍相位信号作为输入，克服了一般环路对导航数据处理的困 难。由于UKF滤波器的估计值是时间段内的均值，数据解调则需要某一时刻 的瞬时值，相比而言滤波器输出值具有一定的滞后性，所以在输出端创新性的 采用一定的补偿方式更准确输出载波相位和频率。下面是采用具体实施方式
一的方法与其它现有方法作对比，说明本发明的 效果，以GPS Ll信号为例根据JPL实验室规定的GPS高动态信号环境的特 点，机动目标的运动速度至少达到1000m/s，加加速度达到100g/s，持续时间 0.5s，加速度斜升50g。由于运动速度要大于1000m/s,那么载体运动所引起的 载波多普勒频移为5.25KHz (GPSL1信号的发射频率为1.5752GHz)，超过了 FLL所能跟踪的范围，造成环路失锁。而其他几种方法与本发明(Openloop) 的动态性能对比如图2所示。图2中的载体运动的加加速度采用lOOg/s的阶 跃形式，持续时间为0.5s。初始运动速度为^00m/s。图2、 3、 4、 5中的曲线 表示的是载体运动引起的多普勒频移和与之基本重合的跟踪结果，足见其跟踪效果好。图6、 7、 8、 9中的曲线表示的是跟踪误差曲线，可见相位误差在零 附近有很小抖动。频率误差由于跟踪的滞后性在阶跃信号出现和消失的地方各 有一个尖峰，这个尖峰高度在2Hz左右，而稳态时的抖动很小，在lHz以下。 综上所述，显示的结果表明本发明所采用的^法可以很好的解决高动态环境下 的载波信号跟踪问题，跟踪精度足以实现正确解调。图10和图11表示的是EKF辅助PLL，以及UKF在闭环跟踪和开环跟踪 当中，相位误差与频率误差的对比。从图中明显看到，当加加速度在5g/s以 上时EKF就无法实现跟踪。原因就是EKF对非线性函数的处理采用泰勒展开 式，舍去了高阶项。对比采用的是精度较高的二阶EKF，具有到二阶变化率 的识别能力，但其对高动态机动目标的加加速度所产生的载&相位三阶变化率 无法正确识别，这种被忽略了的高阶项引起的大的误差导致了 EKF滤波器发 散，无法进行跟踪，而三阶EKF由于求导计算难度太大，计算复杂，目前未 被采用。此外，由于载波信号始终处于一种高速变化当中，环路滤波器PLL 无法达到一个稳态收敛的过程，这将导致环路滤波器的作用不明显，甚至会由 于设计不当引入其他干扰。具体描述本发明的方法步骤，比较性能参数 & (/) = V^"(/)C尸柳[(1 + , - r〗cos[(w + wrf )/+0>0〗+尊)(1) 其中P为信号的接收功率，Z)W为导航数据周期为20ms.， ^为伪码速率， 《为高动态引起的伪码多普勒频移，r为接收扩频码的相位延迟，w,为多普勒 频移，外)为高斯白噪声。GPSL1信号采用的PRN码是码周期长度为1023码 片的C/A码，码频率为1.023MHz。由于伪码信号的频率较低，高动态引起的 多普勒变化《很小，采用传统的码延迟锁相环(DLL)可以进行码跟踪。信号 正常捕获之后进行载波剥离和解扩，再经过低通滤波器得到如下信号形式 W) = ,0cos[(Awrf)"O>0] + iV(0 (2) 其中A为信号幅度，A^是由捕获所得载波频率进行载波剥离后的频率差 值，该值随着时间会快速变化，伴随着较大的一阶和二阶变化率。跟踪环路设 法消除Av^从而实现数据解调。此外，信号的可靠跟踪需要解扩， 一方面去掉 伪码调制所引起的相位翻转，另一方面利用伪码的自相关值提高信噪比。GPSLl信号采用的扩频码周期为lms，所以载波跟踪环路需要每lms调整一次。从UKF的闭环跟踪和开环跟踪对比可见，整个环路结构都属于四阶环路。 一个是由反馈单元构成的闭环形式，UKF采用三维模型；另一个是用算法实 现内部反馈的开环形式，UKF也因此采用四维模型。图中对比可知本发明 提出的UKF开环补偿跟踪方法的跟踪精度明k高于采用UKF闭环形式的跟踪 环路。本方法跟踪得到的相位误差和频率误差随着动态性的提高而略微增大， 而闭环形式却始终保持同一较高数值，这对于中低动态环境来讲是一种损失。 同时这里也证明了本发明采用的方法具有一定的自适应能力。
权利要求
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>ξ</mi> <mo>)</mo></mrow><msub> <mi>R</mi> <mi>b</mi></msub><mi>t</mi><mo>-</mo><mi>τ</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>Φ</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>Φ</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>ξ</mi> <mo>)</mo></mrow><msub> <mi>R</mi> <mi>b</mi></msub><mi>t</mi><mo>-</mo><mi>τ</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>Φ</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>Φ</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。
2、根据权利要求1所述的基于四维UKF的高动态GNSS载波的开环补 偿跟踪方法，其特征在于，步骤五中调整U^F滤波器参数的方法是UKF滤 波器的状态变量为载波相位^、多普勒频移 、多普勒频移一阶变化率^和二 阶变化率^，调整时间间隔^为lms，当前时刻为/fc， UKF滤波器中上述状态 变量满足关系式<formula>formula see original document page 3</formula>式中W为系统噪声;、观测向:<formula>formula see original document page 3</formula>与上述状态向量的关系为:.<formula>formula see original document page 3</formula>及式是R为测量噪声。
3、根据权利要求1所述的基于四维UKF的高动态GNSS载波的开环补 偿跟踪方法，其特征在于，步骤六所述补偿器对UKF滤波器的四个状态变量 的输出补偿是根据公式 '<formula>formula see original document page 3</formula>进行的。
全文摘要
基于四维UKF的高动态GNSS载波的开环补偿跟踪方法，它涉及一种全球导航卫星系统载波的跟踪方法。它解决了现有的GNSS载波跟踪方法跟踪精度低、环路跟踪稳定性差、受带宽限制严重而无法适应较大范围的多普勒频移的问题。其方法捕获数字中频信号的码相位和频率偏移，将频率偏移给UKF滤波器赋初值，将码相位值给本地码跟踪环赋初值并进行本地码跟踪；两路本地载波信号与中频信号相乘后再与本地码跟踪结果进行相关积分，获得两路观测信号后送入UKF滤波器；调整滤波器参数，对高动态的载波信号过程进行跟踪；进行输出补偿，数据解调后进行下一时刻的载波跟踪和数据解调。本发明适用于高动态条件下的信号跟踪过程。
文档编号G01S5/02GK101581775SQ20091007231
公开日2009年11月18日 申请日期2009年6月18日 优先权日2009年6月18日
发明者孟维晓, 王文静, 曦 陈, 帅 韩 申请人:哈尔滨工业大学