[問題] 有關Matlab中FFT範例的問題
最近論文在做跟FFT有關的東西
在Matlab中FFT的範例如下
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sampling period
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
X = S + 2*randn(size(t));
plot(1000*t(1:50),X(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('t (milliseconds)')
ylabel('X(t)')
Y = fft(X);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1)
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')
查過了幾本Matlab有關FFT的地方,如果沒有這樣做的話振幅都滿離奇的
但為何FFT之後的結果,要乘以2/L之後振幅才會是正確的?
麻煩大家幫忙解答了,感謝
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.240.100.35
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1473743226.A.456.html
推
09/13 13:13, , 1F
09/13 13:13, 1F
→
09/15 01:15, , 2F
09/15 01:15, 2F
→
09/15 01:15, , 3F
09/15 01:15, 3F
→
09/15 01:15, , 4F
09/15 01:15, 4F
推
09/15 09:56, , 5F
09/15 09:56, 5F
推
09/16 16:43, , 6F
09/16 16:43, 6F
→
09/16 16:43, , 7F
09/16 16:43, 7F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章