返回信息流% data=wavread('20110902_153557.wav');
% datatemp=data(:,1);
用上面两行读入一个.wav文件后，用下面这个50Hz的陷波器滤波，但是波形没有变化，需要修改什么地方呢？参数怎么设置？求懂的指点一二。。。
==================================================
close all；
wp=[40 60]/100；
ws=[49 51]/100； ℅阻带位于通带中心50Hz
rp=3； ℅通带波纹低于3dB
rs=20； ℅阻带衰减低于20dB
fs=200； ℅采样频率为200Hz
[n，wn]=buttord(wp，ws，rp，rs)；
[h]=butter(n，wn，'stop')；
figure(1)，freqz(h，512，fs)；
t=(0：200)/fs；
s=sin(80*pi*t)+sin(100*pi*t)+sin(120*pi*t)； ℅ 混合正弦
信号
sf=filter(h，1，s)； ℅对信号s 进行滤波
figure(2)， subplot(2， 1， 1)， plot(t*200， s)；
ylabel('magnitude')；
xlabel('Hz')；title(' 混合正弦信号')；
figure(2)，subplot(2，1，2)，plot(t*200，sf)；
ylabel('magnitude')；xlabel('Hz')；title(' 滤波后的信号')；
Y=fft(s，512)；
pyy=Y.*conj(Y)/512；
f=1000/512*(0：255)；
figure(3)，subplot(2，1，1)，plot(f/5，pyy(1：256))；
ylabel('magnitude')；
xlabel('Hz')；title(' 滤波前的频谱')；
Y=fft(sf，512)；
pyy=Y.*conj(Y)/512；
f=1000/512*(0：255)；
figure(3)，subplot(2，1，2)，plot(f/5，pyy(1：256))；
ylabel('magnitude')； xlabel('Hz')；title(' 滤波后的频谱')；
附件(540KB) 20110902_153557.wav
这是一条镜像帖。来源：北邮人论坛 / matlab / #8659同步于 2011/12/14
Matlab机器人发帖
【求助】滤工频（50Hz)。急。。。。。。
lzwno1
2011/12/14镜像同步0 回复
订阅后，新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。