Re: [問題] FFT與IFFT圖形疑問?已回收
ts=0.001; %取樣週期 0.001秒
fs=1/ts; %取樣頻率
t=0:ts:0.6;
x=sin(2*pi*50*t);
xn=x+2*randn(size(t)); %疊加隨機噪音
Y=fft(xn);
y= real(ifft(Y));
f=-fs/2:fs/(length(Y)-1):fs/2;
subplot(311);plot(f,fftshift(Y));
subplot(312);plot(xn);
subplot(313);plot(y);
※ 引述《ice1229 (ice1229)》之銘言：
: 首先這是小弟現在的程式
: t=0:0.001:0.6; %採樣週期為0.001秒
: x=sin(2*pi*50*t);
: xn=x+2*randn(size(t)); %疊加隨機噪音
: Y=fft(xn,512);
: Y(26)=0;
: Y(488)=0;
: y= ifft(Y,512);
: f=1024*(0:256)/512; %設頻率為軸座標,1024是採樣頻率
: subplot(311);plot(f,Y(1:257));
: subplot(312);plot(xn);
: subplot(313);plot(y);
: 目標是首先加入50Hz後fft過後==>完成第一圖形
: Second,把50Hz消除後==>第二圖形
: 最後IFFT inverse回到原始訊號!!==>第三圖形
: 這個是我目前run出來的圖形 http://tinyurl.com/2awfv59
: 感覺上在第一個圖方面好像沒問題 恩恩~50Hz有出來!
: 到了第二個圖,那樣的話算是有消成功嗎?
: 第三個圖,咦,圖形應該是長這樣嗎= =?
: 有大大可以指點迷津嗎? 或是幫我修個程式,感覺圖形也沒有很漂亮= =!
: 感恩(跪....
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 120.126.194.135
※ 編輯: ABCDEFGH1 來自: 120.126.194.135 (05/24 18:43)
※ 編輯: ABCDEFGH1 來自: 120.126.194.135 (05/24 21:37)
※ 編輯: ABCDEFGH1 來自: 120.126.194.135 (05/24 23:01)
討論串 (同標題文章)
本文引述了以下文章的的內容：
完整討論串 (本文為第 2 之 3 篇)：