GPS测站坐标时间序列周期性探测方法及系统

摘要

本发明提供了一种GPS测站坐标时间序列周期性探测方法及系统，包括步骤：步骤1，获取GPS测站坐标时间序列观测值；步骤2，剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差；步骤3，采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型，并设计谐函数模型矩阵；步骤4，采用最小二乘法解算谐函数模型矩阵，探测GPS测站坐标时间序列周期性。本发明可能处理含有时变周期性信号的GPS测站坐标时间序列，并保证探测准确性。

主权项

GPS测站坐标时间序列周期性探测方法，其特征是，包括步骤：步骤1，获取GPS测站坐标时间序列观测值；步骤2，剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差；步骤3，采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型，并设计谐函数模型矩阵其中，y表示GPS测站坐标时间序列观测值；$A=\left[\begin{array}{cc}1& t_1\\ 1& t_2\\ ...& ...\\ 1& t_m\end{array}\right],$$A_k=\left[\begin{array}{cc}\cos {\omega }_k{t}_1& \sin {\omega }_k{t}_1\\ \cos {\omega }_k{t}_2& \sin {\omega }_k{t}_2\\ .& .\\ .& .\\ .& .\\ \cos {\omega }_k{t}_m& \sin {\omega }_k{t}_m\end{array}\right],$t_i为观测历元，i＝1,2,...,m，m为观测历元数量，ω_k表示可变周期性信号的频率；$x_k=\left[\begin{array}{c}{a}_k\\ b_k\end{array}\right],$a_k、b_k分别为正弦函数分量、余弦函数分量；$x=\left[\begin{array}{c}y\left(0\right)\\ a_0\end{array}\right],$y(0)为GPS测站参考历元坐标，a₀为GPS测站的线性速度；k表示谐函数编号，q为谐函数数量；步骤4，采用最小二乘法解算谐函数模型矩阵，获得可变周期性信号的频率估计值，并采用假设检验方法验证频率估计值；步骤5，基于假设检验验证后的频率估计值，采用最小二乘方法解算谐函数模型矩阵获得正弦函数分量a_k、余弦函数分量b_k和线性速度a₀。