[发明专利]一种基于BP神经网络技术的电离层球谐函数改进方法

权利要求书

1.一种基于BP神经网络技术的电离层球谐函数改进方法，其特征在于：包括以下步骤：

S1：获取区域观测点的穿刺点地理经度、穿刺点地理纬度和穿刺点垂直方向的电子含量；

S2：通过电离层球谐函数模型公式建立电离层误差补偿模型，并计算出该穿刺点的垂直方向电子含量以及电离层误差补偿模型的残差值；

S3：利用BP神经网络技术对球谐函数进行模型补偿，建立融合模型；

S4：利用新建立的融合模型对球谐函数模型拟合的电离层垂直总电子含量残差进行预报，进而提高模型补偿精度。

2.根据权利要求1所述的基于BP神经网络技术的电离层球谐函数改进方法，其特征在于：所述步骤S1中，区域观测点的穿刺点地理经度和穿刺点地理纬度分别通过式(1)、式(3)得到：

式(1)中，λ为区域观测点的穿刺点地理经度，λ₀为测站接收机的地理经度，Ψ_pp为地心张角，如式(2)所示，A为卫星方位角，为区域观测点的穿刺点地理纬度；

Ψpp=π2-E-arcsin(RR+Hcos E)---(2)]]>

式(2)中，E为卫星高度角，R为地球半径，H为电离层单层模型的等效高度；

式(3)中，为测站接收机的地理纬度。

3.根据权利要求1所述的基于BP神经网络技术的电离层球谐函数改进方法，其特征在于：所述步骤S1中，区域观测点的穿刺点垂直方向的电子含量由下面方法所得：

假设地球平均半径为R，电离层薄层高度为H，在接收机处卫星的天顶距为z，在穿刺点处的天顶距为z′(天顶距与高度角互为余数)，则有如下关系式：

z′=acsin[RR+Hsinz]---(4)]]>

时钟的基本频率f₀为10.23MHz，GPS上L波段的两种载波L1和L2的频率分别为：

f_L1＝154×f₀＝1575.42MHz(5)

f_L2＝120×f₀＝1227.60MHz (6)

采用伪距观测量时，根据每个历元两个频率上的伪距之差解算出天顶方向自由电子含量VTEC值，其观测方程如下：

Pj-Pi=-40.28fj2VTEC1cos z′+40.28fi2VTEC1cos z′+(BjS-BiS)-(BjR-BiR)---(7)]]>

式(7)中，i,j＝1,2，i≠j；P为伪距观测值，P_j为第j个信号的伪距观测值，P_i为第i个信号的伪距观测值，P_j-P_i为电离层延迟观测量，B^S为伪距观测量的卫星硬件延迟偏差，为第j个信号的伪距观测量的卫星硬件延迟偏差，为第i个信号的伪距观测量的卫星硬件延迟偏差，B^R为伪距观测量的接收机硬件延迟偏差，为第j个信号的伪距观测量的接收机硬件延迟偏差，为第i个信号的伪距观测量的接收机硬件延迟偏差，f_i为载波L1的频率，f_j为载波L2的频率，由此可解出电离层天顶方向总电子含量VTEC的表达式如下：

VTEC=-cos z′40.28fi2fj2fi2-fj2(Pj-Pi-(BjS-BiS)+(BjR-BiR))=-cos z′40.28fi2fj2fi2-fj2(ΔPij-ΔBijS+ΔBijR)---(8)]]>

式(8)中，ΔP_ij是采用伪距观测时的伪距之差，为伪距观测值的卫星相对硬件延迟偏差，为伪距观测值的接收机相对硬件延迟偏差；

采用相位观测量时，其观测方程如下：

式(9)中，λ、N分别为波长和整周模糊度，λ_j为第j个信号的波长，λ_i为第i个信号的波长，N_j为第j个信号的整周模糊度，N_i为第i个信号的整周模糊度，b^S为相位观测值的卫星硬件延迟偏差，为第j个信号的相位观测值的卫星硬件延迟偏差，为第i个信号的相位观测值的卫星硬件延迟偏差，b^R为相位观测值的接收机硬件延迟偏差，为第j个信号的相位观测值的接收机硬件延迟偏差，为第i个信号的相位观测值的接收机硬件延迟偏差；

解方程可得：

式(10)中，L_Δij为相位观测量组合，为载波相位观测值的卫星相对硬件偏差，为载波相位观测值的接收机硬件延迟偏差，Amb_ij为模糊度组合常数。