文档库

最新最全的文档下载
当前位置:文档库 > 数字信号处理实验六 IIR数字滤波器的设计实验报告

数字信号处理实验六 IIR数字滤波器的设计实验报告

1、基于Butterworth 型模拟滤波器原型使用冲激不变转换方法设计数字滤波器,要求具有下面的参数指标:

通带截止频率:π2.0=p w 通带波动值:dB R p 1=

阻带截止频率:π3.0=s w 阻带波动值:dB A s 15=

clear all ;

Rp=1;

As=15;

wp=0.2*pi;

ws=0.3*pi;

T=1;

wap=wp/T;

was=ws/T; %得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式

[N,wn]=buttord(wap,was,Rp,As,'s'); %计算模拟滤波器的最小阶数

[bs,as]=butter(N,wn,'s'); %设计出所需的模拟低通滤波器 [bz,az]=impinvar(bs,as,1/T);

%冲激不变法AF 到DF Rip=10^(-Rp/20);

Atn=10^(-As/20);

[H,w]=freqz(bz,az,512);

mag=abs(H);

db=20*log10(mag/max(mag));

figure;

subplot(2,1,1);

plot(w/pi,mag);

title('幅频特性');

xlabel('w(/pi)');

ylabel('|H(jw)|');

axis([0,1,0,1.1]);

set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 0.5 1]); set(gca,'YTickMode','manual','YTick',[0 Atn Rip 1]);

grid on ;

subplot(2,1,2);

plot(w/pi,db);

title('幅频特性(db) ');

xlabel('w(/pi)');

ylabel('dB');

axis([0,1,-40,5]);

set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 0.5 1]); set(gca,'YTickMode','manual','YTick',[-40 -As -Rp 0]);

grid on ;