В статье предложен метод обнаружения выбросов в погрешностях выработки времени спутниковыми часами (далее – погрешности часов), который основан на взвешенном варианте описания данных с использованием опорных векторов (SVDD – support vector data description) в сочетании с модифицированным медианным абсолютным отклонением (MAD – median absolute deviation), что позволяет устранить ограничения, присущие общепринятому методу MAD. В качестве весов для SVDD применяется локальная достижимая плотность, которая в полной мере учитывает локальные характеристики всех данных. Путем построения минимальной гиперсферы в многомерном пространстве признаков одномерные данные преобразуются в расстояния от векторов до центра гиперсферы, что увеличивает несоответствие между обычными данными и выбросами. Модифицированный метод MAD предполагает проверку расстояния от вектора текущей эпохи до центра гиперсферы на наличие грубых погрешностей, что позволяет выяснить, является ли таковым текущее смещение часов. Для имитационных экспериментов были задействованы точные данные о погрешностях часов спутников BDS-3 с различными интервалами выборки, предоставленные Немецким центром исследования Земли (GFZ). Сравнение результатов применения методов MAD и комбинированного SVDD-MAD (WS-MAD) показывает, что последний способен обнаруживать даже небольшие погрешности в гладко распределенных частотных данных часов и значительные – в данных смещения часов при наличии тренда. Кроме того, обнаружение выбросов и прогнозирование погрешности часов спутников с различными орбитами и типами часов позволяют заключить, что метод WS-MAD повышает их точность в большей степени для спутников МЕО и IGSO, чем GEO. Для тех же МЕО-спутников указанный эффект для водородных часов проявляется сильнее, чем для часов на рубидии.
This paper proposes a clock bias gross error detection method, which combines weight- ed support vector data description (SVDD) with modified median absolute deviation (MAD) to address the limitations of the traditional MAD method. The method uses the local reachable density as the weighting factor of SVDD, which fully considers the local characteristics of each data. By constructing the minimum hypersphere in the high-dimen- sional feature space, the one-dimensional data are transformed into the distance from the vector point to the center of the hypersphere in the high-dimensional space. This trans- formation increases the discrepancy between normal data and gross errors. The modified MAD method probes the distance from the current epoch vector point to the center of the hypersphere for gross errors, and thus determines whether the current clock bias is a gross error. The precision clock bias data of BDS-3 with different sampling intervals provided by GFZ were used for simulation experiments. By comparing the results of the MAD method and the WS-MAD method, it is found that the WS-MAD method can detect the small gross errors in the smooth clock bias frequency data and more gross errors in the clock bias data with trend term floating. The fitting and prediction analyses on satellites with different orbits and clock types show that the WS-MAD method improves the fitting and prediction accuracy of MEO and IGSO satellites better than that of GEO satellites. For the same MEO satellites, the enhancement effect of hydrogen clocks is better than that of rubidium clocks.
The authors declare that there are no conflicts of interest present.