IIR低通滤波的结果分析

2019-07-17 14:46发布

ti, simsun, sans-serif">明明的低通为啥滤波后是空白?高通的也是空白。  求解释。详细点~clear;close all
[z,fs,bits]=wavread('e:1234.wav')
y=z(1:8192)
Y=fft(y)
fp=1000;fc=1200;As=100;Ap=1;Fs=8000
wc=2*fc/Fs;wp=2*fp/Fs
[N,ws]=buttord(wc,wp,Ap,As)
[b,a]=butter(N,ws)
figure(1)
freqz(b,a,512,Fs)
x=filter(b,a,z)
X=fft(x,8192)
figure(2)
subplot(2,2,1);plot(abs(Y));axis([0,1000,0,1.0])
title('滤波前信号频谱')
subplot(2,2,2);plot(abs(X));axis([0,1000,0,1.0])
title('滤波后信号频谱')
subplot(2,2,3);plot(z)
title('滤波前信号波形')
subplot(2,2,4);plot(x)
title('滤波后信号波形')
b_large_AOwX_74d80000b8201262.jpg
b_large_QtcB_2af2000012a81261.jpg
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。
该问题目前已经被作者或者管理员关闭, 无法添加新回复
9条回答
陶醉了的猫
1楼-- · 2019-07-18 11:49
往事只能回首 发表于 2012-6-14 15:30
Fs/2为折叠频率,IIR数字滤波器的频谱分析范围为[0:Fs/2]

高通低通都这样?为什么换成高通的时候就出错误了呢。。。。低通时候Fs/2没问题。高通就报错了
往事只能回首
2楼-- · 2019-07-18 12:42
 精彩回答 2  元偷偷看……
善良的小包子
3楼-- · 2019-07-18 14:28
第二行读文件的时候就已经给出了fs作为采样频率,你在第5行又定义了新的采样频率Fs,不过这个问题不足以导致你图中的情况,我猜测是你的滤波器要求过高了,改变下滤波器各参数试一试

一周热门 更多>