![](https://csdnimg.cn/release/download_crawler_static/87956635/bg8.jpg)
图2 cheby1滤波器设计流程图
程序为:
[y,fs,bit]=wavread('xm.wav');
fp=4300;
fr=4000;
ap=1;
as=40;
[n,fn]=cheb1ord(fp/(fs/2),fr/(fs/2),ap,as,'z');
[b,a]=cheby1(n,ap,fn,'high');
freqz(b,a);
[h,w]=freqz(b,a);
subplot(221);plot(w*fs/(2*pi),abs(h));title('cheby1 高通滤波');
y1=filter(b,a,y);
subplot(222);plot(y);title('原始语音信号');
subplot(223);plot(y1);title('IIR 滤波后语音信号');
X=fft(y1,4096);
subplot(224);plot(abs(X));title('滤波后信号频谱');
显示为: