Дальше ты имеешь данные (длиной N) в виде: Дата(t) смещение_x(dx) смещение_y(dy) (желательно в виде электронной таблицы).
Ежели данные не в виде электронной таблицы, то составляешь эту самую таблицу из данных.

Даты при этом как вариант в формате дат электронной таблицы, поэтому и годовой период надо прописать в этом же формате: Например:

Код:

A1:="01.01.2023"
A2:="01.01.2024"
A3:=A2-A1

В ячейке A3 будет числовое представление года в единицах электронной таблицы для дат.

После этого рассчитываются периодические значения для каждой даты:

Код:

S=sin(t/T*2*pi)
C=cos(t/T*2*pi)

Получаем два столбца чисел, где t - текущая дата, T - годовой период, pi - 3.14....

Из этих двух столбцов получаем амплитудные значения годовых колебаний:

Код:

Ax=sum(dx*S)/sqrt(N)
Bx=sum(dx*C)/sqrt(N)
Ay=sum(dy*S)/sqrt(N)
By=sum(dy*C)/sqrt(N)

Имея амплитуды годовых колебаний вычленяем их из исходного ряда данных:

Код:

dx = dx0 - Ax*S - Bx*C
dy = dy0 - Ay*S - By*C

Более подробно смотри в спектральном анализе (разложение Фурье).

Ну а с полученным, избавленном от годовой периодики, рядом производишь стандартное построение кубического тренда, либо средствами таблиц, либо непосредственно прописывая формулы.