一种GNSS载波相位平滑伪距处理方法

摘要

本发明属于全球卫星导航系统(GNSS)技术领域，公开了一种GNSS载波相位平滑伪距处理方法。本发明通过对平滑常数的解算，通过当前历元的载波相位、下一历元的粗码伪距和载波相位、元载波相位平滑伪距初始值获取的下一历元外推伪距，计算出下一历元的载波相位平滑伪距值。重复上述过程，实时获取GNSS载波相位平滑伪距值。采用本发明提高了载波相位平滑伪距精度，进而大大提高了GNSS定位精度。

主权项

1.一种GNSS载波相位平滑伪距处理方法，其特征在于：包括如下步骤：步骤一、GNSS接收机实时输出粗码伪距和载波相位；步骤二、由GNSS接收机输出的当前历元的粗码伪距，获得载波相位平滑伪距初始值，所述当前历元为接收机输出的第一个粗码伪距；步骤三、根据步骤一获得的当前历元粗码伪距和下一历元载波相位，与步骤二获得的载波相位平滑伪距初始值，获取下一历元的外推伪距P_ex(k)；步骤四、获取下一历元的载波相位平滑伪距值P_sm(k)P_sm(k)＝w(k)×P(k)+[1-w(k)]×P_ex(k)其中，w(k)为前历元k的伪距权重，如果平滑常数M值大于历元k，则伪距权重取w(k)＝1／k，如果平滑常数M值不大于历元k，则伪距权重取w(k)＝1／M；$M=\mathrm{INT}\left[\frac{1}{2}\sqrt[3]{\frac{{\sigma }_{\mathrm{DLL}}^2}{2T_m^2{I}_d^2}}\right]$其中，M为平滑常数，为GNSS粗码伪距的设备测量误差的方差，T_m为数据采样周期，I_d为电离层延迟的发散变化率；P(k)为步骤一获得的当前历元的粗码伪距，P_ex(k)为骤三获取的当前历元外推伪距；步骤五、根据步骤一获得的下一历元粗码伪距和下下一历元载波相位，与步骤四获得的载波相位平滑伪距值，获得下下一历元的外推伪距；步骤六、重复步骤四与步骤五，实时获取GNSS载波相位平滑伪距值。