[发明专利]一种大地测量非等距时序噪声的自适应差分估计方法

权利要求书

1.一种大地测量非等距时序噪声的自适应差分估计方法，其特征在于：包括以下步骤：

(1)读取不同空间大地测量技术观测值文件或IGS站坐标文件中的观测值数据，并生成观测值时间序列；

(2)选取时序移动窗口节点长度作为差分算子长度，确定观测值时序逼近的时间基函数；

(3)基于所选时间基函数构建差分算子估计方程，利用最小二乘方法和最小值搜索准则确定自适应差分算子，并通过差分算子矩阵计算差分序列；

设移动窗口节点长度为K的第i个窗口的差分算子d_i＝[d_i,1…d_i,K]^T，构建差分算子估计方程：

令差分序列满足零均值期望特性，同时考虑非零基准向量e_i＝[e_i1…e_iK]，则可建立差分算子估计方程：

其中：

F_p(t_i+k)为时间基函数；

(4)通过改变差分算子长度获得不同长度下的差分算子和相应的差分序列，利用峰度和偏态检验确定最佳的自适应差分算子及其差分序列；

(5)利用误差传播律进行观测值时序的噪声估计。

2.根据权利要求1所述的大地测量非等距时序噪声的自适应差分估计方法，其特征在于：步骤(2)设移动窗口节点长度为K的第i个窗口内观测值序列的移动时间窗口为t∈[t_i+1,t_i+K]，选取时间基函数F_p(t)；

在移动时间窗口t∈[t_i+1,t_i+K]内，观测值时间序列y(t)可采用时间基函数F_p(t)表示为逼近函数形式：

式中，y(t)表示窗口内的观测值时序，F_j(t)表示j阶时间基函数，p表示时间基函数的最高阶数，β_j表示F_j(t)的确定性系数，ε_t表示观测噪声。

3.根据权利要求2所述的大地测量非等距时序噪声的自适应差分估计方法，其特征在于：所述时间基函数F_p(t)的确定包括以下步骤：

①对移动时间窗口t∈[t_i+1,t_i+K]标准化为[-1,1]，即

其中，窗口内节点序号k∈[1,K]，τ_i+k为t_i+k的归一化时间；

②在以下两种时间基函数中任意选择一种作为时间基函数F_p(t_i+k)，p满足p≤K-2：

A)时间基函数F_p(t_i+k)取p阶一般多项式T_p(τ_i+k)，其具体形式为：

B)时间基函数F_p(t_i+k)取p阶Bernoulli多项式B_p(τ_i+k)，其具体形式为：

其中，组合数b_q表示q阶Bernoulli数，其递推公式为：