[发明专利]一种有效探测GNSS基准站坐标时间序列不连续点的方法

权利要求书

1.一种有效探测GNSS基准站坐标时间序列不连续点的方法，其特征在于，包括以下步骤：

步骤1，输入GNSS基准站坐标时间序列；

步骤2，将GNSS基准站坐标时间序列中相邻时刻的观测值作差得到GNSS基准站坐标差值序列，将GNSS基准站坐标差值序列使用四分位距方法进行探测得到多个初始不连续点；

步骤3，构造GNSS基准站坐标时间序列的参数估计模型，以步骤1所述的GNSS基准站坐标时间序列作为输入数据，以所述GNSS基准站坐标时间序列的参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、步骤2所述的多个初始不连续点作为未知参数，通过优化求解得到未知参数的最优线性无偏估计值及所述GNSS基准站坐标时间序列的残差序列；

步骤4，根据步骤3所述的残差序列构建检验统计量，对步骤2所述的多个初始不连续点进行统计检验，每次从步骤3所述构造的GNSS基准站坐标时间序列的参数估计模型中剔除不显著的初始不连续点，直到所有的初始不连续点经过统计量检验都显著为止，以所述GNSS基准站坐标时间序列的参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、检验显著的不连续点作为新的未知参数，通过步骤3所述优化求解得到新的未知参数的估计值及残差序列；

步骤5，检查步骤4获取的新的未知参数的最优线性无偏估计值X^new与步骤3获取的未知参数的最优无偏估计值X^initial是否相等，若相等则进入步骤6，否则用X^new替换X^initial重复步骤3～步骤4；

步骤6，输出新的未知参数的最优线性无偏估计值X^new，即为GNSS基准站坐标时间序列中不连续点。