[問題] 請問關於消除60Hz雜訊?已回收
小弟目前程式為:
clear
Data = load('101.txt');
data1 = Data(:,2);
x=data1(1:2000)/100;
x=x-mean(x);
Fs=360;
for k=1:length(x)
noise(k,1)=sin(2*pi*60*k/Fs);
end
xn=x+noise;
Y=fft(xn,1024);
f=Fs*(0:511)/1024;
a=abs(Y);
plot(f,a(1:512));
xlabel('Frequency (Hz)')
==========以上可以show出60Hz的頻譜圖==========
接下來想要把60Hz給消掉,show出60Hz消失的頻譜圖
Y(172)=0;
Y(844)=0;
.
.
.
想請問各位前輩我資料為零的點數有算錯嗎?
後面程式應該還要加上甚麼才能show出60Hz消失
的頻譜圖呢?
附上101號資料:http://www.megaupload.com/?d=Z1DBBOQA
m(_ _)m...感恩!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.117.57.152
推
05/29 12:18, , 1F
05/29 12:18, 1F
→
05/29 12:19, , 2F
05/29 12:19, 2F
→
05/29 12:20, , 3F
05/29 12:20, 3F