一种基于双参数混合正则化的电离层掩星反演方法
【技术领域】
[0001] 本发明涉及卫星大地测量及大气探测领域，具体是涉及一种基于双参数混合正则 化的电离层掩星反演方法。
【背景技术】
[0002] 电离层是人类生存的地球近地空间环境的一部分，与人类活动息息相关。充分了 解电离层，是人类对自身生存环境认识和利用的重要基础。它的活动和特性直接影响到无 线电通信、卫星精密导航定位以及无线电技术的航天测控及其他相关领域。因此，监测与研 宄电离层中的各种现象，从而揭示现象背后的物理机制，具有重要的科学、经济、军事意义。
[0003] 电离层探测有着悠久的历史，全球定位系统（Global Positioning System, GPS) 作为当前探测电离层的有效手段已得到广泛应用，通过其获取总电离层含量（Total Electron Content, TEC)和电子密度Ne两个重要物理参量来揭示时空变化规律和特征。最 近发展起来的LEO (Low Earth Orbit) GPS掩星电离层反演技术是薪新的、有效的电离层探 测技术，且仍在不断发展和完善之中。
[0004] GPS掩星电离层反演技术是GPS无线电掩星技术的重要组成部分之一，迄今为止， 国内外许多学者为此进行了大量的研宄。对于它，基于总电子含量的Abel积分变换反演方 法是最为常用的，该方法本质上是一种线性迭代算法，反演过程是从电子密度廓线顶层递 推至底层。在该过程中，如上一层的反演结果存有一定的误差，则此误差在递推过程中得到 放大，从而影响反演的精度。为获取高精度的电离层电子密度，寻求新的方法迫在眉睫。
【发明内容】
[0005] 针对上述存在的问题，本发明的目的是提供一种基于双参数混合正则化方法反演 电离层电子密度，以减小测量误差的影响，从而提高电子密度的反演精度。
[0006] 为达到上述目的，本发明采用如下的技术方案：
[0007] -种基于双参数混合正则化的电离层掩星反演方法，具体步骤如下：
[0008] (1)数据预处理。进行坐标转换以及GPS周跳的探测与修复；
[0009] (2)利用GPS精密星历内插出所需时刻的卫星坐标；
[0010] (3)利用GPS双频观测值计算传播路径上的TEC ;
[0012] 式中，月和為分别为两个频率上相位平滑码伪距观测值，Bs和以分别为码观测中 卫星和接收机的硬件延迟偏差。
[0013] (4)双参数混合正则化方法反演电离层电子密度
[0014] 总电子含量（TEC)为电离层电子密度沿信号传播路径的线积分，表示为：
[0015] TECray= / rayne(l)dl (2)
[0016] 如图2所示的电离层掩星观测几何示意图，AC为斜向TEC路径，D为此路径TEC对 应的掩星切点，位于掩星切点与GPS卫星一侧的DB为非掩星侧，⑶一侧为掩星侧，图2中 B、C点距球心距离相等，低轨卫星（LE0)在B和C点都观测GPS信号，获得TECa^PTEC AC。 若假设电子密度为局部对称分布，则AB段的TEC在B和C点观测是时相同，则有
[0018] 其中，ne为电子密度，r为积分路径上的点到地心的距离，r(l为掩星切点到地心的 距离，即碰撞参数，为LE0卫星到地心的距离。
[0019] 对（3)式进行离散化处理，得
[0021] 其中，A为线性转换矩阵，h为第j层的电子密度。在离散化形式下，即使j = i 时，A也是一个病态矩阵，其谱不集中。为获得问题的稳定近似解，需对（4)式进行正则化 过程，即求解下述的极小化问题
[0023] 其中，a和0分别为大于0的正则化参数；第一个正则化项||A?-可看作是 引入了对参数n的约束，％是参数的先验值；L为离散化后与正则算子相对应的正则矩阵， 在此取为二阶正则矩阵
[0024]
[0025] 第二个正则化项取为||?|(。
[0026] 本发明反演中采用正则化思想，用正则解+辦rk/TEC+afW)作为精 确解的近似。正则化参数的选取采用广义交叉校检准则进行计算，然后再进行优化组合，以 此得到最佳的正则化参数。
[0027] 本发明技术有益效果：
[0028] 本发明包括：（1)对卫星观测数据进行预处理。（2)利用双频GPS载波相位观测值 计算传播路径上的TEC。（3)通过双参数混合正则化方法反演电离层电子密度。为确定最 佳正则化参数，本发明采用UPRE法、L曲线法和GCV法的优化组合来确定。提高了计算电 离层电子密度的精度，可大大减小测量误差的影响，从而提高电子密度的反演精度。
[0029] 在空间环境探测、通信导航等技术领域具有重要应用价值。
【附图说明】
[0030] 图1是本发明实施例1的电离层掩星反演流程图；
[0031] 图2是本发明实施例1的电离层掩星观测几何图；
[0032] 图3是图1中的双参数混合正则化反演流程图；
[0033] 图4是第一例利用不同方法反演电离层掩星电子密度剖面与测高仪站的比较个 例图；
[0034] 图5是第二例利用不同方法反演电离层掩星电子密度剖面与测高仪站的比较个 例图。
【具体实施方式】
[0035] 下面结合附图和具体实施例对本发明作进一步的描述，然而本发明的范围并不限 于下述实施例。
[0036] 为了验证本发明基于双参数混合正则化的电离层掩星反演方法的可行性以及相 对于Tikhonov正则化和Abel积分变换反演方法的优越性，选取了 2011年8月6日COSMIC 电离层掩星数据进行反演。并将反演剖面与测高仪JicamarCa(76° W，11.9° S)站探测结 果进行比较。
[0037] (1)获取2011年8月6日COSMIC电离层掩星观测数据以及测高仪Jicamarca站 数据；
[0038] (2)数据预处理。GPS周跳的探测与修复；
[0039] (3)利用GPS双频观测值计算出随时间变化传播路径上的TEC
[0041] 式中，戶t和為分别为两个频率上相位平滑码伪距观测值，B s和B 8分别为码观测中 卫星和接收机的硬件延迟偏差。
[0042] (4)计算出掩星碰撞参数及其相应的炜度、经度和高度；
[0043] (5)将坐标转换至地固参考系下，并计算在该坐标系下碰撞参数及其相应的炜度、 经度和高度；
[0044] (6)根据LE0、GPS卫星位置信息以及掩星事件表，将随时间变化的TEC转换为随碰 撞参数变化的TEC。如图2所示的电离层掩星观测几何示意图，AC为斜向TEC路径，D为此 路径TEC对应的掩星切点，位于掩星切点与GPS卫星一侧的DB为非掩星侧，⑶一侧为掩星 侦牝图2中B、C点距球心距离相等，低轨卫星（LE0)在B和C点都观测GPS信号，获得TECab 和TECa。。若假设电子密度为局部对称分布，则AB段的TEC在B和C点观测是时相同，则有
[0046] 其中，ne为电子密度，r为积分路径上的点到地心的距离，r ^为掩星切点到地心的 距离，即碰撞参数，为LE0卫星到地心的距离。
[0047] (7)利用双参数混合正则化方法反演电离层电子密度廓线。对（2)式进行离散化 处理，得
[0049] 其中，A为线性转换矩阵，rij为第j层的电子密度。在离散化形式下，即使j=i 时，A也是一个病态矩阵，其谱不集中。为获得问题的稳定近似解，需对（3)式进行正则化 过程，即求解下述的极小化问题
[0051] 其中，a和0分别为大于〇的正则化参数；第一个正则化项|I(?-可看作是 引入了对参数n的约束，％是参数的先验值；L为离散化后与正则算子相对应的正则矩阵， 在此取为二阶正则矩阵
[0052]
[0053] 第二个正则化项取为。
[0054] 本发明反演中采用正则化思想；
[0055] 用正则解
作为精确解的近似。正则化参 数的选取采用广义交叉校检准则进行计算，然后再进行优化组合，以此得到最佳的正则化 参数。
[0056] (8)图4和图5给出了不同时间利用不同方法反演结果对比个例图，本发明反演的 电子密度与Tikhonov正则化和Abel积分变换反演结果的对比情况，其结果证明了本发明 的可行性和优越性。
[0057] (9)通过以上计算方案可得如下结论，本发明基于双参数混合正则化方法可大大 减小测量误差的影响，从而提高电子密度的反演精度。
[0058] 本发明涉及的所有参数符号的具体含义为：
[0059] 月、戶:分别为两个频率&和f2上相位平滑码伪距观测值；
[0060]Bs、BK分别为码观测中卫星和接收机的硬件延迟偏差；
[0061] ne为电离层电子密度；
[0062] r为积分路径上的点到地心的距离；
[0063] &为掩星切点到地心的距离，即碰撞参数；
[0064] r为LE0卫星到地心的距离；
[0065]A为线性转换矩阵；
[0066] 11』为第j层的电子密度。
[0067] 以上所述仅为本发明的较佳实施例而已，凡在本发明的精神和原则之内所作的任 何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种基于双参数混合正则化的电离层掩星反演方法，其特征在于，包括w下步骤： (1) 对卫星观测数据进行预处理，将其进行坐标转换W及GI^S周跳的探测与修复，为后 处理提局精度提供保障； (2) 利用GI^S精密星历内插出所需时刻的卫星坐标； (3) 利用GPS双频观测值计算传播路径上的TEC;(1) 式中，聲和4分别为两个频率上相位平滑码伪距观测值，BS和B"分别为码观测中卫星 和接收机的硬件延迟偏差； (4) 双参数混合正则化方法反演电离层电子密度； 总电子含量TEC为电离层电子密度沿信号传播路径的线积分，表示为： TECrw=/rayDeWdl似 通过电离层掩星观测得到几何示意图，选定B、C双参数，AC为斜向TEC路径，D为此路 径TEC对应的掩星切点，位于掩星切点与GPS卫星一侧的DB为非掩星侧，CD-侧为掩星侧， B、C点距球屯、距离相等，低轨卫星（LEO)在B和C点都观测GPS信号，获得TECw和TECAC; 若假设电子密度为局部对称分布，则AB段的TEC在B和C点观测是时相同，则有其中，n。为电子密度，r为积分路径上的点到地屯、的距离，r。为掩星切点到地屯、的距离， 即碰撞参数，ru。为LEO卫星到地屯、的距离； 对（3)式进行离散化处理，得(4) 其中，A为线性转换矩阵，rij.为第j层的电子密度；在离散化形式下，即使jSi时，A也是一个病态矩阵，其谱不集中；为获得问题的稳定近似解，需对（4)式进行正则化过程， 即求解下述的极小化问题其中，a和P分别为大于0的正则化参数；第一个正则化项|Z(n-所)|可看作是引入 了对参数n的约束，n。是参数的先验值；L为离散化后与正则算子相对应的正则矩阵，在此 取为二阶正则矩阵第二个正则化项取为HI:; 上述反演中采用正则化思想； 用正则解6 =(心^ +af￡ +片。、方TEC+aZ7Zn。）作为精确解的近似；正则化参数的 选取采用广义交叉校检准则进行计算，然后再进行优化组合，W此得到最佳的正则化参数。
【专利摘要】本发明公开了一种基于双参数混合正则化的电离层掩星反演方法，其特征在于，包括：(1)对卫星观测数据进行预处理，获取观测值；(2)利用GPS精密星历内插出所需时刻的卫星坐标；(3)利用GPS双频观测值计算传播路径上的TEC；(4)双参数混合正则化方法反演电离层电子密度；本发明提高了计算电离层电子密度的精度，可大大减小测量误差的影响，从而提高电子密度的反演精度。
【IPC分类】G06F19/00
【公开号】CN104933316
【申请号】CN201510382191
【发明人】汤俊, 李长春
【申请人】华东交通大学
【公开日】2015年9月23日
【申请日】2015年7月2日