Вот собственно эти две формулы используются, где как я понял N - число выборок за период, и к - номер грамоники (у меня N=80 к=1). Тогда Cos(2*PI*k/N) = Cos(2*PI/80) = 0,999999. Меня вполне устроит такая точность, если я приму это за 1. Тогда vi = x[i] + 2*vi_1 - vi_2; Xk = vi_1*vi_1 - 2*vi_1*vi_2 + vi_2*vi_2; или если брать только положит.корни то Xk = vi_1 - vi_2;
И чтобы получить амплитуду надо еще умножить на 2/N.