[发明专利]一种顾及周期信号的GNSS测站速度稳健估测方法

权利要求书

1.一种顾及周期信号的GNSS测站速度稳健估测方法，其特征在于：包括以下步骤，

步骤1，获取单个测站同一分量m个多源GNSS时间序列，m表示多源GNSS时间序列的来源个数；

步骤2，获取各单一来源GNSS时间序列的周期信号，得到的信号周期集合记作其中，n₁，…，n_m为各GNSS时间序列中的周期信号个数；p表示信号周期，其上标表示获得该信号周期的GNSS时间序列编号，下标为周期编号；

步骤3，筛选信号周期，使所有周期均满足以下条件，得到速度样本窗口列表1，...，速度样本窗口列表m，

其中，表示信号周期，i＝1，…，n_j′，j＝1，…，m，i和j为循环变量，分别用于表示任一周期和任一序列来源；T_j表示各序列跨度，单位为年；b为预先设置的崩溃点；

设n_j′表示筛选后剩余的信号周期数量，筛选后速度样本窗口列表记为

步骤4，按下式计算各单一来源序列速度样本，得到速度样本列表1，...，速度样本列表m，

其中，j＝1，…，m，分别为第j个序列的第k个和第l个历元，分别对应与历元下的坐标；与历元的时间间隔与对应的速度样本窗口最接近；

将获得的速度样本列表1，...，速度样本列表m分别相应记作o₁，…，o_m为各序列获得的速度样本个数；

步骤5，计算各单一来源序列速度样本的绝对中位差MAD^j，并估计标准差σ^j，

σ^j＝1.4826·MAD^j

其中，j＝1，…，m，median{ }表示计算集合的中位数；

步骤6，精化速度样本，使多源序列获取的所有速度样本均满足以下条件，得到精化的速度样本窗口列1，...，精化的速度样本窗口列m，

将精化后速度样本记为o₁′，…，o_m′为各序列精化后速度样本个数；

步骤7，计算各单一来源序列精化后速度样本的中位数及其速度样本方差VR^j，

VR^j＝(1.4826·MAD^j′)²

其中，j＝1，…，m，MAD^j′为精化后速度样本的绝对中位差；

步骤8，计算各单一来源序列速度中位数的权值w^j，

其中，j＝1，…，m，VR^h表示任一个速度样本方差值；

步骤9，计算单测站单分量最终速度估值及不确定性如下，

所得结果用于测站周边地理环境监测。

2.根据权利要求1所述顾及周期信号的GNSS测站速度稳健估测方法，其特征在于：预设b＝0.15。

3.根据权利要求1所述顾及周期信号的GNSS测站速度稳健估测方法，其特征在于：当n₁′，…，n_m′中任一值小于5，假设n_a′＜5，a∈{1，…，m}，且对应的周期信号集合内均不满足时，在筛选后速度样本窗口列表中加入一个取值为365的元素。

4.根据权利要求1或2或3所述顾及周期信号的GNSS测站速度稳健估测方法，其特征在于：所述测站周边地理环境监测，为测站周边地震灾害监测、滑坡灾害监测或水文环境监测。

5.一种顾及周期信号的GNSS测站速度稳健估测系统，其特征在于：用于实现如权利要求1-4任一项所述的一种顾及周期信号的GNSS测站速度稳健估测方法。