[发明专利]适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法

说明书

技术领域

本发明涉及全球导航系统领域，特别是涉及适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法。

背景技术

GNSS电离层建模研究和应用中，主要分为两类电离层模型：一类是经验模型，包括GPS中使用的Klobuchar模型、IRI模型、Bent模型等，经验模型比较复杂且修正效果较差，一般不适用于高精度的电离层延迟建模和GNSS定位中；另一类是基于高精度GNSS双频观测数据构建的拟合电离层模型，根据建模区域的不同，一般常用的电离层模型包括下面三种：多项式模型(POLY)、球谐函数模型(SHF)、三角级数(TSF)模型。然而，现有技术中这三种电离层模型的电离层延迟的拟合精度较低，因此有必要研究出一种拟合精度高的模型。

发明内容

发明目的：本发明的目的是提供一种能够解决现有技术中存在的缺陷的适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法。

技术方案：为达到此目的，本发明采用以下技术方案：

本发明所述的适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法，包括以下步骤：

S1：获取区域观测点的穿刺点地理经纬度信息和相应时刻的电离层垂直总电子含量信息；

S2：通过已知数据转化获得太阳时角差、地磁纬度和地方时的信息；

S3：选定模型参数，构建计算量与所选参数的方程，建立区域电离层垂直总电子含量模型；

S4：通过最小二乘法计算改进的多项式模型参数，完成区域电离层垂直总电子含量模型的构建。

进一步，所述步骤S1中，区域观测点的穿刺点地理经纬度信息通过式(1)和式(2)得到：

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

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

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

进一步，所述步骤S2中，太阳时角差、地磁纬度和地方时分别通过式(4)、(5)和(6)得到：

ΔS＝(λ-λ₀)+(t-t₀)(4)

式(4)中，ΔS为太阳时角差，λ为穿刺点经度，λ₀为测站接收机的地理经度，t为观测时刻点，t₀为建模时间段的中间时刻点；

式(5)中，为地磁纬度，为北磁极的地理纬度，为区域观测点的穿刺点地理纬度，λ_i为区域观测点的穿刺点地理经度，λ_S为北磁极的地理经度；

τ＝UTC+λ/15°(6)

式(6)中，τ为地方时，UTC为协调世界时。

进一步，所述步骤S3中，区域电离层垂直总电子含量模型如式(7)所示：

式(7)中，VTEC为区域电离层垂直总电子含量，a₀₀、a₀₁、a₀₂、a₁₀、a₁₁、a₁₂、a₂₀、a₂₁、a₂₂、b₀₁、b₀₂、b₀₃、b₀₄和b₀₅均为待求系数，ΔS为太阳时角差，为穿刺点距离区域中心纬度的差值，为地磁纬度，τ为地方时。

有益效果：本发明公开了一种适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法，该方法很好的保留了传统多项式拟合模型在空间曲面拟合上的优势，并顾及了电离层延迟信息随时间的非线性振动特性，较好的模拟了区域电离层的时空变化特征，相比传统的多项式模型精度得到了提高。

附图说明

图1为本发明具体实施方式的改进多项式模型与传统的9-DPM模型的建模效果对比图；

图1(a)为前55条检验数据；

图1(b)为后56条检验数据；

图2为本发明具体实施方式的0:00UTC时刻电离层延迟实际空间分布图；

图3为本发明具体实施方式的0:00UTC时刻POLY-SIN模型的改正偏差的空间分布图；