[发明专利]基于相关系数的削弱坐标时间序列中CME影响的方法

权利要求书

1.基于相关系数的削弱坐标时间序列中CME影响的方法，其特征在于，包括步骤：

步骤1，获取GPS测站坐标时间序列观测值，并获取GPS测站的残差坐标时间序列；

步骤2，从GPS网内选定基准站s，计算基准站s与GPS网内普通测站间的公共历元，并确定公共历元个数；所述的普通测站为GPS网内除基准站s外的测站；

步骤3，计算公共历元内基准站s与普通测站p间的相关系数r_sp；

步骤4，根据基准站s与普通测站p间的相关系数获取GPS网内所有普通测站共同作用于基准站s的共模误差其中，S-1为GPS网内参与共模误差计算的普通测站数目；r_sp为公共历元内基准站s与普通测站p间的相关系数；v_p,k和分别为第k个公共历元下普通测站p坐标时间序列的残差和标准差；

步骤5，从基准站s的原始坐标时间序列观测值中对应扣除普通测站共同作用于基准站s的共模误差。

2.如权利要求1所述的基于相关系数的削弱坐标时间序列中CME影响的方法，其特征在于：

所述的GPS测站的残差坐标时间序列采用如下公式获得：

y(t)=a+bt+csin(2πt)+dcos(2πt)+esin(4πt)+fcos(4πt)+Σj=1nggjH(t-Tgj)+vi]]>

其中：

y(t)为时刻t对应的GPS测站坐标观测值；

t表示日坐标解历元，单位：年；

a为GPS测站位置，b为线性速度；

系数c、d用来描述GPS测站的年周期运动，系数e、f用来描述GPS测站的半年周期运动，c、d、e、f为待估计参数，经拟合获得；

为跳变改正项，g_j表示跳变振幅，T_gj表示发生跳变的历元，n_g表示跳变个数，j为跳变编号，H为海维西特阶梯函数（Heaviside step function），在跳变前H值为0，跳变后H值为1；

v_i为时刻t的观测值残差。

3.如权利要求1所述的基于相关系数的削弱坐标时间序列中CME影响的方法，其特征在于：

所述的公共历元内基准站s与普通测站p间的相关系数r_sp为：

rsp=Σk=1N[(mk-m‾)(nk-n‾)]Σk=1N(mk-m‾)2Σk=1N(nk-n‾)2]]>

其中：

N为普通测站p与基准站s间公共历元个数，k表示公共历元编号；

m_k、n_k分别表示基准站s与普通测站p第k个公共历元的残差坐标，可通过残差坐标时间序列获取；

和分别为基准站s与普通测站p在公共历元下残差坐标时间序列均值。