资源描述：
《实验六数字滤波器设计》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验六数字滤波器设计13通信工程方鹏豪1341737081)基于Butterworth模拟滤波器原型，使用双线性变换法分别设计IIR数字低通滤波器，参数指标：通带截止频率Wp=0.2it，通带波动值Rp=ldB,阻带截止频率Ws=0.3n,阻带最小衰减Rs=20dB。滤波器采样频率Fs=2000Hz。要求：绘出幅频特性、相频特性图并比较结果，绘出滤波器零极点图，写出滤波器的系统函数。代码••%数字滤波器指标Wp=0.2*pi:Ws=0.3*pi:Rp=l;Rs=20;%转换为模拟滤波器指标Fs=2000;T=l/F
2、s;0nigp=(2/T)*tan(Wp/2);%双线性变换法屮频率预畸变0mgs=(2/T)*tan(Ws/2):%模拟原型滤波器计算[N,Omgc]=buttord(Omgp,Omgs,Rp,Rs,’s’)％计算阶数和截止频率[zO,pO,kO]=buttap(N);%归一化原型设计bO=kO*real(poly(zO))%求原型滤波器分子分母系数aO=real(poly(pO));[b,a]=lp21p(b0,aO,Omgc)%变换为模拟滤波器分子分母系数%用双线性变换法计算数字滤波器系数[num,den]=
3、bilinear(b,a,Fs)%数字滤波器设计结果%方法2:直接求模拟滤波器系数[b,a]=butter(N,Omgc,’s’)；[num,den]=bilinear(b,a,Fs);%方法3:直接求数字滤波器系数[N,Wc]=buttord(Wp/pi,Ws/pi,Rp,Rs);[num,den]=butter(N,Wc):%求数字滤波器的频率特性[H,w]=freqz(num,den);db_H=20*logl0((abs(H)+eps)/max(abs(H)));%转换为分贝值subplot(2,2,1);
4、plot(w/pi,abs(II));ylabel(’幅度’)；xlabel(’频率’)；title(’幅度响应’)；axis([0,1,0,1.1]);set(gca,’XTickMode’,’manual’，’XTick’，[0,0.2,0.3,1]);ripple=10"(-Rp/20);rs=l(T(-Rs/20);set(gca,’YTickMode’，’manual’，’YTick’，[0，rs,ripple,1]);gridsubplot(2，2，2);plot(w/pi,angle(II)/pi);
5、ylabel(’相位’)；xlabel(’频率’)；title(’相位响应’)；axis([0,1，-1，1]);set(gca,’XTickMode’，’manual’，’XTick’，[0，0.2,0.3，1]);set(gca,’YTickMode’，’manual’，’YTick’，[-1,0，1]);gridsubplot(2,2,3);plot(w/pi,db_H);ylabelC幅度(dB)’);xlabel(’频率’);title(’幅度响应(dB)’);axis([0,1,-50,5]);set(
6、gca,’XTickMode’，’manual’，’XTick’，[0，0.2，0.3，1]);set(gca,’YTickMode’，’manual’subplot(2,2,4);zplane(num,,’YTick’，[-50,-20,-1,0]);gridden);title(’零极点图’)；0.00120.0002axis([-l.1,1.1,-1.1,运算图形结果:nujn=0.00020.00121.1]);0.00370.00620.00620.0037den=1.0000-3.84716.7860-
7、6.96294.4452-1.75420.3944■0.03882)基于切比雪夫1型模拟滤波器原型设计，使用脉冲响应不变法设计一个11R数字带通滤波器，参数指标：下通带截止频率Wpl=0.3it，上通带截止频率WP2=0.7，通带最大衰减Rp=ldB;下阻带截止频率Wsl=0.2n，上阻带截止频率Ws2=0.8it，阻带最小衰减Rs=20dB。滤波器采样频率Fs=2000Hz。要求：绘出幅频特性、相频特性图并比较结果，绘出滤波器零极点图，写出滤波器的系统函数。代码：%数字滤波器指标Wpl=0.3*pi;Wp2=0.
8、7*pi;Wsl=0.2*pi;Ws2=0.8*pi;Rp=l;Rs=20;%转换为模拟滤波器指标Fs=2000;T=l/Fs;Omgpl=Wpl*Fs;0mgp2=Wp2*Fs;Omgsl=Wsl*Fs;0mgs2=Ws2*Es;0mgp=[0mgpl,0mgp2];0mgs=[0mgsl,0mgs2];bw=0mgp2-0mgpl;w0=sqrt(0mg
