1.一种北斗接收机非差观测值噪声评估方法,其特征在于,包括以下步骤:
1)读取观测值文件的北斗卫星导航系统观测数据,获取B1、B2、B3三个频率的载波相位观测值和伪距观测值;
2)通过理论分析,利用步骤1)中的载波相位观测值求取B1、B2、B3三个频率的载波相位观测值组合序列;利用步骤1)中的伪距观测值,求取B1、B2、B3三个频率的伪距观测值组合序列;
3)通过2)中理论分析获得的三频载波相位系数对步骤1)中的载波相位观测值进行周跳探测与标记;
4)利用所述载波相位观测值组合序列和伪距观测值组合序列分别获得B1、B2、B3各个频率上的载波相位和伪距观测值噪声。
2.根据权利要求1所述的北斗接收机非差观测值噪声评估方法,其特征在于,所述步骤1)中,B1、B2、B3三个频率的载波相位观测值 和伪距观测值P的观测方程表示为如下形式:;
;
其中,下标i代表频率标识,i=1、2、3,分别对应B1、B2、B3频率;Li为第i频率上的载波相位观测值,以米为单位;λi为第i频率的载波波长; 为原始观测值文件中第i频率的载波相位观测值,单位为周;ρ为卫星至北斗接收机距离,单位为米;c为真空光速;dt为北斗接收机钟差;dT为卫星钟差;dorb为卫星轨道误差;dtrop为对流层误差;μi为第i频率的电离层误差系数且μi=f12/fi2;f1为B1载波频率;fi为第i频率上的载波频率;dion为B1频率的电离层延迟误差;Ni为第i频率载波相位模糊度,以周为单位; 为第i频率载波相位硬件延迟偏差; 为第i频率载波相位观测多路径误差; 为第i频率载波相位观测值噪声;Pi为原始观测值文件中第i频率伪距观测值,单位为米; 为第i频率伪距硬件延迟偏差; 为第i频率伪距观测多路径误差; 为第i频率伪距观测值噪声。由于载波相位硬件延迟偏差与模糊度参数无法分离,因此该偏差项在下面的论述过程中未单独列出。
3.根据权利要求2所述的北斗接收机非差观测值噪声评估方法,其特征在于,B1、B2、B3三个频率的载波相位观测值组合序列为:其中,Nc为B1、B2、B3三个频率组合观测值模糊度,以米为单位; 为B1、B2、B3三个频率载波相位组合观测值多路径误差; 为B1、B2、B3三个频率载波相位组合观测值噪声;B1、B2、B3三个频率伪距观测值组合序列为:其中, 为组合伪距观测硬件延迟偏差; 为三频伪距组合观测值多路径误差;
为B1、B2、B3伪距组合观测值噪声。
4.根据权利要求3所述的北斗接收机非差观测值噪声评估方法,其特征在于,通过以下关系式进行载波相位观测值周跳探测:△Nc(k)=Lc(k)-Lc(k-1)
;
其中,△Nc(k)为周跳检测量;k表示观测历元,即时刻;Lc(k)、Lc(k-1)为第k历元与k-1两相邻历元三频载波相位组合观测值; 为周跳判断阈值; 为B1频率上的载波相位观测值精度,设为2毫米;若第k历元△Nc大于 则认为该历元出现周跳,并对该历元进行标记。
5.根据权利要求4所述的北斗接收机非差观测值噪声评估方法,其特征在于,通过以下关系式计算B1、B2、B3三个频率的载波相位观测值噪声和伪距观测值噪声:
其中, 为B1频率载波相位观测值噪声;n为载波相位观测量或伪距观测量统计个数; 为剔除模糊度值的载波相位观测值残差; 分别为B2、B3频率上的载波相位观测值噪声;σB1/P、σB2/P、σB3/P分别为B1、B2、B3频率上的伪距观测值噪声;
为剔除硬件延迟偏差值的伪距观测值残差。
6.根据权利要求5所述的北斗接收机非差观测值噪声评估方法,其特征在于, 的计算过程为:设北斗接收机的载波相位观测值统计时长为t小时,北斗接收机的采样间隔为s秒,将所述统计时长分成每m分钟一段,则对每 个Lc计算一个平均值,然后用所述 个Lc减去所述平均值,得到m分钟的 个 依此类推,一共得到个 序列,记为 的计算与上面的过程类似:设北斗接收机的伪距
观测值统计时长为t小时,北斗接收机的采样间隔为s秒,将所述统计时长分成每m分钟一段,则对每 个Pc计算一个平均值,然后用所述 个Pc减去所述平均值,得到m分钟的 个 依此类推,一共得到 个 序列,记为 此处的即为上面公式中的载波相位观测量或伪距观测量统计个数n。
7.根据权利要求6所述的北斗接收机非差观测值噪声评估方法,其特征在于,m=3~
15。