[問題] 關於FFT繪圖的問題已回收
我有兩筆資料,一筆是時間,一筆是對應時間的data,每一筆有1146項
我想把data轉FFT去觀察頻率,我已經寫好程式
但我不太確定把時間刻度轉成頻率上的刻度
我這樣寫不知道正不正確,轉出來的圖總覺得怪怪的
資料內容:
data=[1.0124 1.0368 1.0057 1.0368 1.0124 ..........]
time=[0.0660 0.4960 0.9260 1.2680 1.6260 ..........]
------------------------------------------------------------------
程式:
data=load('data.txt');
time=load('time.txt');
N=length(time);
dt=(time(length(time))-time(1))/N;
df=1/(time(length(time))-time(1));
f=0:df:(N-1)*df;
yp=fft(data);
yp=fftshift(yp);
subplot(2,1,1);
plot(time,data, '-');
title('Sinusoidal signals');grid on
xlabel('Time (seconds)'); ylabel('Amplitude');
axis tight
subplot(2,1,2);
plot(f, abs(yp), '.-b'); grid on
xlabel('Frequency)');
ylabel('Magnitude (Linear)');
轉出來的圖
http://yfrog.com/18fftj
請高手指點一下
謝謝大家
matlab新手
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 60.249.132.173
→
06/12 21:53, , 1F
06/12 21:53, 1F
→
06/13 10:44, , 2F
06/13 10:44, 2F
→
06/13 10:48, , 3F
06/13 10:48, 3F
討論串 (同標題文章)
以下文章回應了本文 (最舊先):
完整討論串 (本文為第 1 之 5 篇):