|
回复 #2 UCLAN 的帖子
clear;
[fid,message]=fopen();
data=fread(fid,'float32');
fclose(fid);
da=data
fs=6000
Rp = 1; Rs = 10;
Wp = [1500 2500]/(fs/2); Ws = [1400 2600]/(fs/2);
[n,Wn] = buttord(Wp,Ws,Rp,Rs)
[b,a] = butter(n,Wn);
[h,w]=freqz(b,a)
figure(1)
subplot(211);
plot(fs*w/2/pi,abs(h));
xlabel('频率(HZ)');
ylabel('幅值');
grid on
sf=filter(b,a,da)
subplot(212);
plot(sf)
axis('tight')
xlabel('时间 (seconds)');
ylabel('幅值');
这是一个Butterworth带通滤波器,希望对你能有帮助,具体参数还是要进一步按照不同的要求修改
[ 本帖最后由 zhangnan3509 于 2007-12-15 22:31 编辑 ] |
评分
-
1
查看全部评分
-
|