《信号处理实验4网上资料.docVIP

/blog/static/2758818820090433142767/ 3．设计数字滤波器，画出其频率响应曲线 各滤波器的性能指标： (1)低通滤波器性能指标　fp＝1000Hz, fs＝1200Hz, As＝40dB, Ap＝1dB。 (2)高通滤波器性能指标　fs＝4000Hz, fp＝4300Hz, As＝40dB, Ap＝1dB (3)带通滤波器性能指标　fp1＝1200 Hz, fp2＝3 000 Hz, fs1＝1000 Hz, fs2＝3200 Hz, As＝40dB, Ap＝1dB。 要求： (1)频率变换法设计IIR滤波器： 可以利用函数butter、cheby1、cheby2和ellip等设计。 (2)(选做)设计IIR滤波器：用双线性变换法设计上面要求的3种滤波器。bilinear (3)(选做)设计FIR滤波器：用窗函数法设计上面要求的3种滤波器。可以利用函数fir1设计FIR滤波器。 (4)函数freqz画出各滤波器的频率响应。 低通滤波器： 1(频率变换法设计IIR滤波器)butter函数 [y,fs,bit]=wavread(1); fp=1000; fr=1200; ap=1; as=40; [n,fn]= buttord(fp/(fs/2),fr/(fs/2),ap,as,z); [b,a]=butter(n,fn); freqz(b,a); [h,w]=freqz(b,a); subplot(221);plot(w*fs/(2*pi),abs(h));title(butter低通滤波); y1=filter(b,a,y); y2=fftfilt(b,y); subplot(222);plot(y);title(原始语音信号); subplot(223);plot(y1);title(IIR滤波后语音信号); subplot(224);plot(y2);title(FIR滤波后语音信号); sound(y1,fs,bit); 2(频率变换法设计IIR滤波器)cheby1函数 [y,fs,bit]=wavread(1); fp=1000; fr=1200; ap=1; as=40; [n,fn]=cheb1ord(fp/(fs/2),fr/(fs/2),ap,as,z); [b,a]=cheby1(n,ap,fn); freqz(b,a); [h,w]=freqz(b,a); subplot(221);plot(w*fs/(2*pi),abs(h));title(cheby1低通滤波); y1=filter(b,a,y); y2=fftfilt(b,y); subplot(222);plot(y);title(原始语音信号); subplot(223);plot(y1);title(IIR滤波后语音信号); subplot(224);plot(y2);title(FIR滤波后语音信号); sound(y1,fs,bit); 3(频率变换法设计IIR滤波器)cheby2函数 [y,fs,bit]=wavread(1); fp=1000; fr=1200; ap=1; as=40; [n,fn]=cheb2ord(fp/(fs/2),fr/(fs/2),ap,as,z); [b,a]=cheby2(n,ap,fn); freqz(b,a); [h,w]=freqz(b,a); subplot(221);plot(w*fs/(2*pi),abs(h));title(cheby2低通滤波); y1=filter(b,a,y); y2=fftfilt(b,y); subplot(222);plot(y);title(原始语音信号); subplot(223);plot(y1);title(IIR滤波后语音信号); subplot(224);plot(y2);title(FIR滤波后语音信号); sound(y1,fs,bit); 得出来的信号波形如下： 4(频率变换法设计IIR滤波器)ellip函数 [y,fs,bit]=wavread(1); fp=1000; fr=1200; ap=1; as=40; [n,fn]=ellipord(fp/(fs/2),fr/(fs/2),ap,as,z); [b,a]=ellip(n,ap,as,fn); freqz(b,a); [h,w]=freqz(b,a); subplot(221);plot(w*fs/(2*pi),abs(h));title(ellip低通滤波); y1=filter(b,a,y); y2=fftfilt(b,y); subplot(222);plot(y);title(原始语音信号); subplot(223);plot(y1);tit