使用matlab fft 做震動訊號分析問題!!!!

matlab fft

louis80392 2018-12-15 00:52:44 ‧ 2121 瀏覽

大家好!! 請問當我有一筆1024個數據的震動訊號時，我應該要如何使用fft轉換成頻譜?
我有試過別人提供的代碼

t = 0:0.001:0.6;
x = sin(2*pi*50*t)+sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(1000*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
Y = fft(y,512);
Pyy = Y.* conj(Y) / 512;
f = 1000*(0:256)/512;
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')

可是他們都是自己創立一個函數再做fft，我照著他們的做沒辦法跑出圖形，請問有人可以幫我解答一下，我用加速計測量到的訊號，應該經過哪些步驟來轉換成頻譜??