实验四IIR数字滤波器设计.docxVIP

实验四 IIR 数字滤波器设计 一、实验目的 1、熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法； 2、熟悉用脉冲响应不变法设计 IIR 数字滤波器的原理与方法。 二、实验内容 1、已知低通滤波器的指标为： 通带边缘频率： 0.4π， Ap=0.5dB； 阻带边缘频率： 0.6π， As=50dB； a、 采用脉冲响应不变法设计巴特沃斯， T=1. 画出幅度响应和数字滤波器 的脉冲响应 h(n)； b、 b、采用脉冲响应不变法设计巴特沃斯， T=1. 画出幅度响应和数字 滤波器的脉冲响应 h(n)。 2、用双线性变换法设计低通滤波器，满足技术指标 wp=0.2π， Ap=0.25dB； ws=0.4π， As=50dB， 并对方波信号进行滤波，画出滤波前后的波形图并进行简要分析 3、设计一个数字高通滤波器 H(z)，它用在结构 xa(t ) A/D H(z) D/A ya(t ) 中，满足下列要求： a、采样速率为 10kHZ； b、阻带边缘频率为 1.5kHZ，衰减为 40dB； c、通带边缘频率为 2kHZ，衰减为 3dB； d、单调的通带和阻带。 4、设计一个带阻滤波器，要求通带上下截止频率为 0.8 π，0.2 π，通带内衰减 不大于 1dB，阻带起始频率为 0.7 π，0.4 π，阻带内衰减不小于 30dB。 设计巴特沃斯带阻滤波器并画出该数字高通滤波器的幅度响应和脉冲响应。 三、实验程序及解释和实验分析及图形 1a、 clear;close all; %清屏 wp=0.4*pi; ws=0.6*pi; Ap=0.5 As=50; T=1; Fs=1/T; OmegaP=wp/T; OmegaS=ws/T; [cs,ds]=afd_butt(OmegaP,OmegaS,Ap,As); [b,a]=impinvar(cs,ds,Fs); [h,w]=freqz(b,a); subplot(2,2,1);plot(w/pi,abs(h)); title( 幅度响应 ); grid; subplot(2,2,2);plot(w/pi,angle(h)); title( 相位响应 ); grid; subplot(2,2,2);plot(w/pi,angle(h)); title( 相位响应 ); grid; subplot(2,2,3); plot(w/pi,20*log(abs(h))); title( 幅度响应 dB);grid; n=[0:1:59]; imp=[1;zeros(59,1)]; y=filter(b,a,imp); subplot(2,2,4);plot(n,y); title( 脉冲响应 ); grid; wp=0.4*pi;%ws=0.6*pi;% wp=0.4*pi; % ws=0.6*pi; % Ap=0.5; % As=50; % T=1; % 1b、 clear;close all; Fs=1/T; 计算通带截止频率原型 计算阻带截止频率原型巴特沃斯原型滤波器计算双线性不变法变换 分析：因为 w=Ω t,w 与Ω呈线性关系， 所以其相位响应图是线性的， 如相位响应 图所示；因为其设计的是低通滤波器，所以会把其高频部分滤掉，留下低频，从 而不是特别陡峭，如脉冲响应图所示。 通带边缘频率设为 0.4*pi 阻带边缘频率设为 0.6*pi 通带最大衰减为 0.5dB阻带最小衰减为 50dB把周期 T 设为 1 OmegaP=(2/T)*tan(wp/2); % OmegaS=(2/T)*tan(ws/2); % [cs,ds]=afd_butt(OmegaP,OmegaS,Ap,As);% [b,a]=bilinear(cs,ds,Fs); % [h,w]=freqz(b,a); 响应 subplot(2,2,1);plot(w/pi,abs(h)); % title( 幅度响应 ); grid; subplot(2,2,2);plot(w/pi,angle(h)); % title( 相位响应 );grid; subplot(2,2,3);plot(w/pi,20*log(abs(h)));% title( 幅度响应 dB); axis([0,1,-80,5]);grid; n=[0:1:59]; imp=[1;zeros(59,1)]; y=filter(b,a,imp); subplot(2,2,4);plot(n,y); title( 脉冲响应 ); % % % % % 计算数字滤波器的 Z 变换频率 画幅频响应曲线图 画出相频响应曲线图 画出取其对数的幅频响应曲线 设置坐标轴，加网格线 设n值 输入矩形波 对 imp 进行滤波 画出 y 的图形 线，而是曲线。由于其设