一种适用于全球区域的gnss电离层延迟三维建模方法
【技术领域】
[0001] 本发明属于卫星导航、无线电波传播以及空间电离层领域，具体涉及一种适用于 全球区域的GNSS电离层延迟三维建模方法。
【背景技术】
[0002] 经过 20 多年的发展，GNSS(Global Navigation Satellite System)不仅使用在 军事领域，也广泛地应用在民用领域。它以精确、定时、不受地域限制等优点，成为在陆地交 通、航海和航空中通用的导航工具，具有全能性、全球性、全天候、连续性和实时性的导航、 定位和定时功能，能为用户提供精密的三维坐标、速度和时间。GNSS从一问世起就不是一个 单一星座系统，目前包含了美国的GPS、俄罗斯的GL0NASS、中国的BDS和欧盟的Galileo等 系统。
[0003] 电离层会对穿过其中的GNSS卫星信号产生附加延迟。用户终端通过传输时间进 行测距，电离层产生的时间延迟在终端被转化为距离误差。电离层延迟所引起的误差占总 测距误差的绝大部分，是终端定位解算中的主要误差源。精确可靠的电离层时延修正已成 为进一步提升GNSS系统服务性能亟待解决的主要技术问题之一。
[0004] 常用的全球区域电离层延迟二维建模方法忽略了电离层电子密度在垂直方向的 变化，假设电子分布在某一指定高度的薄层上，并且假设某交叉点处的电离层电子浓度总 含量（TEC)是各向同性的，使用投影函数将视线方向的TEC(STEC，Slant TEC)与垂直方向 的TEC(VTEC，Vertical TEC)进行转换。上述模型非常有利于描述电离层电子密度的整层 变化，对分析穿越整个电离层的GNSS信号受到电离层的影响是非常有效的，大大简化了数 据处理的过程，在电离层活动较为平静的中炜度地区，上述假设是基本成立的；但是，对于 电离层活动剧烈的赤道地区或者是"赤道异常"双峰结构的边缘地区，交叉点南北两侧的电 离层TEC变化梯度在低高度角时具有较大的差异，若采用简单的投影函数描述视线与垂直 方向电离层TEC之间的关系，将会带来较大的误差。
【发明内容】
[0005] 有鉴于此，本发明提供了一种适用于全球区域的GNSS电离层延迟三维建模方法， 利用经验正交函数对电离层电子密度在垂直方向的分布进行表征，有效利用了经验正交函 数优良的数学特性，提高了全球区域GNSS电离层延迟三维建模的整体精度与可靠性。
[0006] -种适用于全球区域的电离层延迟的三维建模方法，包括如下步骤：
[0007] 步骤1、获取表征电离层电子垂直分布的经验正交函数，具体过程为：
[0008] S1、使用国际参考电离层模型IRI2012对与观测数据同一天的电离层电子剖面进 行采样，得到一天中不同时刻电离层在不同高度上的电子密度采样数据，构成电离层剖面 时间序列；
[0009] S2、求取电离层剖面时间序列的协方差矩阵；
[0010] S3、对协方差矩阵进行特征值分解得到各个特征值对应的经验正交函数；
[0011] 步骤2、构建电离层三维模型，具体为：
[0012] 用球谐函数表征电离层电子密度的水平方向分布，挑选步骤1得到的部分经验正 交函数，用挑选出的各个经验正交函数的线性组合表示电离层电子密度的垂直方向分布； 对得到的电离层电子密度进行积分即得到电子含量TEC的表达式，其中，积分上、下限分别 为地面GNSS接收机的位置和空间GNSS卫星位置；
[0013] 电离层电子密度的积分形式即为电离层延迟三维模型；其中的经验正交函数的系 数以及球谐函数的系数组成待求解的电离层三维模型的系数；
[0014] 步骤3、解算电离层延迟三维模型系数，具体为：
[0015] Sl、获取X个具有不同接收机位置和/或卫星位置的电子含量TEC的观测量；
[0016] S2、将各个电子含量TEC观测量与具有相应的地面GNSS接收机位置和空间GNSS 卫星位置分别作为积分上、下限的电子含量TEC表达式联立，得到具有X个方程的非线性方 程组；
[0017] S3、采用非线性估计方法对所述非线性方程组进行解算，得到电离层三维模型的 系数后，建立适用于全球区域的GNSS电离层延迟三维模型。
[0018] 较佳的，所述步骤1中获得电离层剖面时间序列的方法为：将一天的时间划分为 多个时间段，在每个时间段内，获取电离层剖面不同高度上的电子密度数据，则按顺序排列 的所有时间段的电子密度数据即为电离层剖面时间序列；其中不同高度上的电子密度数据 对应的经、炜度相同或者不同。
[0019] 较佳的，所述步骤2中挑选经验正交函数的方法为：将特征值按从大到小顺序排 列，将前k个特征值对应的经验正交函数进行线性组合用以表示电离层电子密度垂直方向 的分布；其中，k的取值由精度要求和计算速度要求决定。
[0020] 较佳的，所述非线性估计方法为：第一次迭代时，对电离层延迟三维模型的系数赋 初值，分别代入到非线性方程组的各个方程中，得到电子含量TEC的计算值后，与对应的观 测量求差值，得到当前的电子含量误差后再对各电离层延迟三维模型系数分别求偏导，得 到各电离层延迟三维模型系数的误差值；将该误差值与本次迭代的电离层延迟三维模型系 数求和值，作为下一次迭代的电离层延迟三维模型系数值，再次得到电子含量TEC的计算 值后，继续与对应的观测量求差值，以此类推，进行多次迭代，直到电离层延迟三维模型系 数的误差值满足要求后，停止迭代，当前的电离层延迟三维模型系数即为最终的电离层三 维模型方程系数。
[0021] 本发明具有如下有益效果：
[0022] (1)本发明的一种适用于全球区域的GNSS电离层延迟三维建模方法，首先采用国 际参考电离层模型IRI2012获取用于表征电离层电子垂直分布的经验正交函数，基于球谐 函数和经验正交函数使用预先已经精确估计的电离层TEC观测量建立电离层模型，构建适 用于全球区域电离层延迟的精确三维建模方法。
[0023] 相对于现有方法，该方法使用经验正交函数表征电离层电子密度的垂直分布，不 仅克服了电离层延迟二维模型的局限，而且可以有效的监测电离层大尺度结构（如赤道异 常、中炜槽及不同的电离层扰动等）的变化信息，对于监测和研究电离层不同尺度不均匀 性、电离层环境、电离层时空变化、电离层异常变化及电离层全球范围的监测均有重要的科 学意义和应用价值。
【具体实施方式】
[0024] 下面举实施例对本发明进行详细描述。
[0025] 本发明提供了一种适用于全球区域的电离层延迟精确三维建模方法，包含表征电 离层电子垂直分布的经验正交函数获取、电离层三维模型观测方程构建、电离层延迟三维 模型系数解算三大步骤，具体步骤如下：
[0026] 步骤1、表征电离层电子垂直分布的经验正交函数获取
[0027] 本方法使用国际参考电离层模型IRI2012获取电离层延迟三维建模所需的经验 正交函数，对电子密度剖面在不同地区、不同时间以及不同高度上进行采样，获取一个时间 序列的协方差矩阵，采用特征值分解的方法获取此协方差矩阵的特征值和特征向量，具体 包括如下步骤：
[0028] (1)获取与观测数据同时期的电离层剖面时间序列
[0029] 使用国际参考电离层模型IRI2012对与观测数据同一天的电子剖面进行采样，即 分析一天的数据，时间分辨率为T小时，剖面水平经炜度分辨率为L° X L°，高度分辨率 为H米，共得到N个电子密度样本采样数据，其中经炜度分辨率和高度分辨率根据用户需求 进行设定。令各样本采样点的电子密度为Xi(z_j)，其中，i = 1，2,…N ;j = 1，2,…M。则N 个电子密度样本采样数据可表示成如下的矩阵形式：
[0031] 这样就得到与观测数据同时期的电离层剖面时间序列，其中每一行表示一天中获 得的相同高度上的电子密度数据，每一列表示同一时间段中获得的不同高度上的电子密度 数据，其中，不同高度上的数据有可能取自相同的经炜度。
[0032] (2