[发明专利]一种北斗接收机非差观测值噪声评估方法

说明书

技术领域

本发明涉及一种北斗接收机非差观测值噪声评估方法。

背景技术

我国自主研发的北斗区域卫星导航系统已正式建成并开始向亚太地区提供定位、导航、授时等服务，而北斗区域导航系统的各项应用都基于接收机产生的观测值，故观测值的质量将决定上述各项应用的精度。一般使用零基线实验通过载波相位、伪距观测值构成测站间单差或者测站间卫星间双差对北斗接收机观测值噪声进行评估。使用零基线进行北斗接收机观测值噪声分析需要两台接收机共用一个天线进行观测，而同一型号的接收机观测值质量也可能存在一定差异，两台接收机的载波相位、伪距观测值构成单双差的方法，不是直接的评价单一接收机的观测值噪声水平，仅是从侧面反映接收机的观测值噪声。除此之外，零基线单差或双差法完全消除了天线端的噪声，不能准确的反映实际观测值的噪声水平。对于单接收机观测数据噪声的评价，有学者基于单接收机的双频载波相位观测值进行无几何组合，通过滤波进行载波相位的观测值噪声分析，但该组合未完全消除电离层误差的影响，且滤波过程较为复杂。对于接收机伪距观测噪声的评价，一般利用载波相位和伪距组合计算伪距观测值多路径效应，该方法除了会引入了载波相位观测值，增添周跳探测与修复的麻烦，另外部分多路径误差会融入到噪声中。

发明内容

本发明所要解决的技术问题是，针对现有技术不足，提供一种简单易行，且能客观的对北斗单个接收机载波相位和伪距观测值的噪声情况进行评价的方法。

为解决上述技术问题，本发明所采用的技术方案是：一种北斗接收机非差观测值噪声评估方法，包括以下步骤：

1)读取观测值文件的北斗卫星导航系统观测数据，获取B1、B2、B3三个频率的原始载波相位观测值和伪距观测值；

2)通过理论分析，利用步骤1)中的B1、B2、B3三个频率的载波相位观测值获得组合序列；利用步骤1)中B1、B2、B3三个频率的伪距观测值获得组合序列；

3)根据2)中的三频载波相位组合系数对载波相位观测值进行周跳探测与标记；

4)利用所述载波相位观测值组合序列和伪距观测值组合序列分别获得B1、B2、B3各个频率上的载波相位和伪距观测值噪声。

本发明的具体实现步骤如下：

1)读取RINEX观测值文件的北斗卫星导航系统观测数据，获取B1、B2、B3三个频率的载波相位和伪距观测值。根据电离层散射特性，在仅顾及一阶电离层误差的情况下，北斗卫星导航系统的载波相位和伪距观测方程可以表示为如下形式：

其中，下标i代表频率标识，对应取值为1、2、3，分别代表B1、B2、B3；L_i为第i频率上的载波相位观测值，以米为单位；λ_i为第i频率的载波波长；为第i频率的载波相位观测值，单位为周；ρ为卫星至接收机距离，单位为米；c为真空光速；dt为接收机钟差；dT为卫星钟差；d_orb为卫星轨道误差；d_trop为对流层误差；μ_i为第i频率电离层误差系数且μ_i＝f₁²/f_i²；f₁为B1载波频率；f_i为第i频率上的载波频率；d_ion为B1频率电离层延迟误差；N_i为第i频率载波相位模糊度，以周为单位；为第i频率载波相位硬件延迟偏差；为第i频率载波相位观测多路径误差；为第i频率载波相位观测值噪声；P_i为第i频率伪距观测值，单位为米；为第i频率伪距观测多路径误差；为第i频率伪距硬件延迟偏差；为第i频率伪距观测值噪声。由于载波相位硬件延迟偏差与模糊度参数无法分离，因此该偏差项在下面的论述过程中未单独列出。

2)根据理论分析的系数，使用1)中的载波相位观测值进行北斗三频组合。理论上载波相位观测值噪声为波长的1/100，因此不同频率的载波相位观测值噪声与其波长存在一定的比例关系。通过以下关系式获取北斗三频载波相位组合系数：