[发明专利]一种基于北斗的地基增强系统载波相位平滑伪距方法

权利要求书

1.一种基于北斗的地基增强系统载波相位平滑伪距方法，其特征在于，包括以下步骤：

步骤1，根据北斗卫星导航系统接收机伪距与载波相位观测值计算载波差量，并与计算得到的载波差量门限值进行比较，如果载波差量大于载波差量门限值，说明载波相位发生周跳，则初始化滤波器，否则根据伪距与载波相位观测值计算伪码-载波偏离度值，并与设置的伪码-载波偏离度门限值进行比较；

步骤2，如果伪码-载波偏离度大于门限值，则说明电离层发生异常，此时平滑窗口时间为T1，如果伪码-载波偏离度不大于门限值，说明电离层较为平稳，此时设置平滑窗口时间为T2，利用Hatch滤波器对伪距进行平滑。

2.根据权利要求1所述的方法，其特征在于，步骤1中，根据如下公式计算载波差量temp：

temp＝|λ(φ_k-φ_k-1)-(ρ_k-ρ_k-1)|，

式中，λ为载波波长，φ_k和φ_k-1分别为历元k时刻的载波相位观测值和k-1时刻的载波相位观测值，ρ_k和ρ_k-1分别为k时刻的伪距观测值和k-1时刻的伪距观测值。

3.根据权利要求2所述的方法，其特征在于，步骤1中，根据如下方法计算得到载波差量门限值：

步骤1-1，计算得到每颗卫星每个历元时刻的载波差量；

步骤1-2，按卫星仰角划分载波差量区间，以10°为一个区间进行划分，并分别计算每个区间内载波差量的均值μ_elev与标准差σ_elev；

步骤1-3，利用5阶多项式拟合其它角度均值与标准差，其表达式为：

y＝A₅x⁵+A₄x⁴+A₃x³+A₂x²+A₁x+A₀，

式中，A_i为需要拟合的第i个参数，i取值为0～5，x为输入卫星仰角，其取值分别为5°、15°、25°、35°、45°、55°、65°、75°、85°，y为需要拟合的载波差量均值或者标准差；

步骤1-4，对步骤1-1得到的载波差量进行归一化；

步骤1-5，统计归一化后的载波差量，绘制归一化后载波差量分布直方图，并根据每个间隔点内的载波差量数据个数与载波差量总个数计算其概率密度值，即利用间隔内载波差量数据个数除以载波差量总个数，由此得到其概率密度值；

步骤1-6，计算膨胀系数f；

步骤1-7，根据得到的膨胀系数f以及均值μ_elev与标准差σ_elev获得该组样本数据的载波差量门限值Threshold_elev：

Threshold_elev＝μ_elev±6fσ_elev。

4.根据权利要求3所述的方法，其特征在于，步骤1-4中，采用如下公式对样本参数进行归一化：

式中，x表示载波差量采样数据，表示归一化后的载波差量值。

5.根据权利要求4所述的方法，其特征在于，步骤1-6包括：根据步骤1-5计算得到的概率密度分布情况绘制样本数据的概率密度分布图，并在其中绘制标准正态分布概率密度曲线，然后对标准正态分布的概率密度曲线进行膨胀，直至膨胀所得到的曲线完全包络样本数据的概率密度曲线两侧的样本数据，从而得到膨胀系数f。