功率譜密度和FFT

Question

我知道這可能是一個基本問題，但我在Matlab中應用fft函數時遇到了限制。 我從這個規範的時間序列：

bdf=0.005; 
fHighCut=0.2; 
maxTime=600; 
freq=1/maxTime:df:fHighCut; 
w=2*pi*freq; 
time=linspace(0,600,length(freq)*10); 

我的時間序列lloks像：

它的長度是400秒。我想從它創建PSD但現在看來，我獲得對稱的，以及像：

用以下命令：

timestepFFT=1/time(end); 
freqFFT=(1:length(time))*1/time(end); 
amplitudeEtaSeries=abs(fft(etaSeries)/length(time)); 
powerSpectrumEtaSeries=amplitudeEtaSeries.^2/timestepFFT; 
powerSpectrumEtaSeries(1)=0; 

你能幫幫我嗎？ THanks

Answer 1

請閱讀這篇文章。通過使用fft的功率譜密度計算總是產生對稱的PSD。所以忽略fft輸出的第二部分和細節可以在這裏找到： http://www.mathworks.se/help/signal/ug/psd-estimate-using-fft.html