Спектральный анализ временного ряда Matlab

[vintaj~]

Задача: выявить сезонную составляющую(СС) временного ряда. СС представляет собой гармоническую функцию, с помощью МБПФ строю спектр сигнала, но получается довольно странная картина. Сезонная составляющая отсутствует, хотя должна обязательно быть. Помогите найти ошибку:

]B=[605281;588766;588854;619040;616270;617959;619601;622267;602953;599164;621322;623076;618577;621612;623093;602032;596682;623036;626237;627384;624973;626504;604649;595122;624765;628500;629024;629145;627528;611012;604645;627193;619767;618282;616177;618971;603108;595998;617560;620233;620839;623284;622506;612085;605421;631246;628480;625417;627146;630114;615449;612750;632746;641805;646257;653437;651664;630002;623895;641005;636995;631881;0;0;0;0];
FftL=124;
Fd=64;
signal=B;
X=abs(fft(signal,FftL));%преобразование Фурье
Pxx=X.*conj(X)./FftL;% Нормировка спектра по амплитуде
F=0:Fd/FftL:Fd/2-1/FftL;% Массив частот вычисляемого спектра Фурье
figure% Создаем новое окно
subplot(2,1,1);% Выбор области окна для построения
plot(F,X(1:length(F)));% Построение спектра Фурье сигнала

]

[ura_111]

Для начала, я бы советовал подсчитать в другом месте...

Посмотри такое тебе надо?

https://www.youtube.com/watch?v=JUY52vCuj5g

[ura_111]

И расскажи, что такое "сезонную составляющую" и как ты видешь что её нет?

[vintaj~]

что то подобное. Данных много поэтому, хотелось бы получить решение непосредственно методом преобразования Фурье.

[ura_111]

А там что вывело?

Что "сезонную составляющую" есть или нет? Или ты не пробывал?

[vintaj~]

Какая то гармоническая функция или сумма гармоник. Зашумленная разумеется. Ну спектральный анализ говорит, что эта функция sin(0*x*2pi) чего быть не должно

[ura_111]

Когда что то не получается - надо исключить что в самих данных нет ошибки.

[ura_111]

Покажи (пусть даже схематично), что должно получиться...

[vintaj~]

Судя по графику какая то периодичность присутствует. Хотя... возможно, что слишком мало данных

[ura_111]

А нули в конце нужно учитывать?