БПФ на данных временных рядов - PullRequest
Новая
       31

БПФ на данных временных рядов

0 голосов
/ 13 апреля 2019
  1. Какой должна быть частота выборки для набора данных, который собирается каждый месяц и непрерывно в течение 265 лет?
  2. Как определить частоту повторения, увидев БПФ?

Я взял период выборки за 30 * 24 * 60 * 60 секунд, потому что данные записываются только раз в месяц!Когда я строю БПФ, я получаю доминирующие частоты как 3,28x10 ^ -8 Гц.Теперь, как мне это интерпретировать?Как сказать, как часто появляется этот сигнал?

Набор данных находится здесь . 

datax=VarName1;
datay=VarName2;
datay_mean=datay-mean(datay);
class(datax)
dataxx=datetime(datax);
class(dataxx)
figure;
plot(dataxx,datay)
datetick('x', 'keepticks','keeplimits')
title('actual plot of temp in years')

Y = fft(datay_mean);
L = length(datay_mean); 

T = 30*24*60*60 
Fs = 1/(T);

f_adj1=(0:L-1)/L;
figure;
plot(f_adj1,abs(Y/L))                        
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('normalized (freq)')
ylabel('|P1(f)|')

f_adj2=(0:L-1)*Fs/L;
figure;
plot(f_adj2,abs(Y/L))                         
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('(freq)')
ylabel('|P1(f)|')
...