[发明专利]一种GNSS时间序列阶跃探测与修复方法

摘要

本发明公开了一种GNSS时间序列阶跃探测与修复方法，针对GNSS时间序列阶跃问题，利用时间序列发生阶跃前后连续时间窗口内的趋势不变特性，通过移动窗口法和时间基函数保系数趋势项逼近方式构建了阶跃探测与修复模型，进而采用稳健最小二乘、广义最小二乘或其他最优估计理论进行阶跃参数估计和方差‑协方差估计，并通过统计假设检验进行阶跃参数显著性检验与阶跃修复。本发明适用于非等时间间隔或数据缺失下多种类型的GNSS时间序列探测与修复，阶跃探测与修复估计理论严密、易于编程，探测结果具有最优统计性质、修复结果具有统计可靠性。

主权项

1.一种GNSS时间序列阶跃探测与修复方法，其特征在于：该方法包括以下步骤：/n(1)获取GNSS基准站的时序文件，读取原始观测值；所述原始观测值指的是时序文件天顶方向的GNSS站坐标时间序列、GNSS组合钟差时间序列或GNSS组合相位观测时间序列；/n(2)根据GNSS时间序列类别和移动窗口法，选取阶跃前后连续时间窗口的移动窗口节点长度，确定时间基函数及其最高阶数；初始化移动窗口序号i＝1；/n(3)基于GNSS时间序列的趋势不变特性，将第i个移动窗口内的GNSS时间序列进行时间基函数保系数趋势项逼近；/n(4)构建基于最优估计理论的GNSS时间序列阶跃探测与修复模型，并进行阶跃参数估计与方差-协方差估计；/n(5)根据参数估计及其方差-协方差估计，构建正态分布统计量，进行GNSS时间序列阶跃参数的显著性检验，若通过显著性检验，则修复GNSS时间序列阶跃项，并令修复后的时间序列替换原时间序列；若未通过显著性检验，则表示GNSS时间序列不存在阶跃项，无需修复；/n(6)将移动窗口序号i增加为i+1，重复上述步骤(3)～步骤(5)，直至完成整个时间序列的阶跃探测和修复。/n