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-17 21:51
fs本身存储的就是采样频率吧,为什么要定义新的采样频率Fs?
shawon
2楼-- · 2019-07-17 16:55
是不是截止频率选的不合适啊?
凌晨3点睡
3楼-- · 2019-07-17 21:34
FS如果是8192乘以2呢
陶醉了的猫
4楼-- · 2019-07-18 03:45
 精彩回答 2  元偷偷看……
往事只能回首
5楼-- · 2019-07-18 04:01
Fs/2为折叠频率,IIR数字滤波器的频谱分析范围为[0:Fs/2]
陶醉了的猫
6楼-- · 2019-07-18 08:01
善良的小包子 发表于 2012-6-14 15:29
fs本身存储的就是采样频率吧,为什么要定义新的采样频率Fs?

哪儿?第二行的?

一周热门 更多>