2010-11-10 28 views

回答

4

fft可以做到這一點:

Y = fft(signal); 
periodLength = 'second'; %or whatever units your signal was acquired in. 
N = length(Y); 
Y(1) = []; 
power = abs(Y(1:floor(N/2))).^2; 
nyquist = 1/2; 
freq = (1:floor(N/2))/floor(N/2)*nyquist; 
figure; 
plot(freq,power) 
grid on 
xlabel(['cycles/' periodLength]); 
title('Frequency plot'); 

,這給它作爲每個週期的週期數。

period = 1./freq; 
figure; 
plot(period,power) 
grid on 
ylabel('Power') 
xlabel(['Period (' periodLength 's/Cycle)']); 
title('Period Plot'); 
相關問題