[問題]FFT後保留某頻率能量，再作IFFT

[jiahunglin (haha)]

搜尋過IFFT的文章
好像跟我要做的都不太一樣
所以把我的問題po上來，請板上的高手指點
謝謝

我現在有一筆2922*1的資料dis_from_coast
想利用FFT找出變化頻率
參照matlab的demo寫成

relNums=dis_from_coast;
Y = fft(relNums);
Y(1)=[];

n=length(Y);
power = abs(Y(1:floor(n/2))).^2;
nyquist = 1/2;
freq = (1:n/2)/(n/2)*nyquist;

figure(1)
FREQ=freq*365;    %demo是一年一筆資料，我是一天一筆，所以乘以365
plot(FREQ,power)
xlabel('cycles/year')
ylabel('power')
title('Periodogram')

如此可以得到主要變化頻率，也就是power最大的頻率
結果大約是0.99966cycles/year，約一年一次的頻率

如果我現在想要將這個主要頻率的power保留
然後將其他頻率的power歸零
接著再作IFFT得到轉換後的時間序列資料

請問該怎麼作比較好?

感恩

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 210.240.12.220
→ qeagle:除了那一根和對應的負頻保留,其餘設為0,再做ifft就好了
→ qeagle:不過因為LEAKAGE的關係,還原回來的東西可能和你想的不太同