实验五IIR数字滤波器设计与滤波附思考题程序.docVIP

试验五 IIR数字滤波器设计和滤波 1.试验目标 （1）加深对信号采样了解， （2）掌握滤波器设计方法； （3）复习低通滤波器设计。 2．试验原理 现在，设计IIR数字滤波器通用方法是先设计对应低通滤波器，然后再经过双线性变换法和频率变换得到所需要数字滤波器。模拟滤波器从功效上分有低通、高通、带通及带阻四种，从类型上分有巴特沃兹（Butterworth）滤波器、切比雪夫（Chebyshev）I型滤波器、切比雪夫II型滤波器、椭圆（Elliptic）滤波器和贝塞尔（Bessel）滤波器等。 经典模拟低通滤波器指标以下：分别为通带频率和阻带频率，分别为通带和阻带容限（峰波纹值）。在通带内要求,有时指标由通带最大衰减和阻带最小衰减给出，定义以下： 和 第二种常见指标是用参数和A表示通带和阻带要求，图所表示： 二者之间关系为：和，依据这多个参数可导出另外两个参数d，k，分别称为判别因子和选择性因子。 BUTTERWORTH低通滤波器：幅度平方函数定义为，N为滤波器阶数，为截止频率。当初，有，为3DB带宽。 BUTTERWORTH低通滤波器系统函数有以下形式： 由模拟滤波器设计IIR数字滤波器，必需建立好s平面和z平面映射关系。使模拟系统函数变换成数字滤波器系统函数，通常采取冲激对应不变法和双线性变换法。冲激对应不变法存在频谱混叠现象，双线性变换法消除了这一线象，在IIR数字滤波器设计中得到了更广泛应用。 s平面和Z平面映射关系为,将和待入数字频率和等效模拟频率之间映射关系：，因为二者不是线性关系，所以称为预畸变。 3.试验内容及其步骤 试验步骤： （1）给定数字滤波器幅度对应参数。 （2）用预畸变公式将数字滤波器参数变换为对应等效模拟滤波器参数。 （3）采取模拟滤波器设计方法设计等效模拟滤波器 （4）采取双线性变换公式把等效模拟滤波器映射为所期望数字滤波器。 其中第三步中模拟滤波器设计步骤为： 首先，依据滤波器指标求选择因子k和判别因子d 其次，确定满足技术所需滤波器阶数N, 再次，设3db截止频率 最终由表查出归一化巴特沃斯滤波器系数。 设计举例： 例1 设计一个模拟巴特沃特低通滤波器，它在30rad/s处含有1dB或愈加好波动，在50rad/s处含有最少30dB衰减。求出级联形式系统函数，画出滤波器幅度响应、对数幅度响应、相位响应和脉冲响应图。 MATLAB参考程序： Wp=30;Ws=50;Rp=1;As=30; %技术指标 Ripple=10^(-Rp/20); Attn=10^(-As/20); [b,a]=afd_butt(Wp,Ws,Rp,As) %巴特沃兹低通滤波器子程序 [db,mag,pha,w]=freqs_m(b,a,50); %计算幅频响应 [ha,x,t]=impulse(b,a); %计算模拟滤波器单位脉冲响应 figure(1);clf; subplot(2,2,1);plot(w,mag);title(Magnitude Response); xlabel(Analog frequency in rad/s); ylabel(H); axis([0,50,0,1.1]);grid; subplot(2,2,2);plot(w,db);title(Magnitude in dB); xlabel(Analog frequency in rad/s); ylabel(decibels); axis([0,50,-40,5]) grid subplot(2,2,3);plot(w,pha/pi);title(Phase Response); xlabel(Analog frequency in rad/s); ylabel(radians); axis([0,50,-1.1,1.1]) grid subplot(2,2,4);plot(t,ha);title(Impulse Response); xlabel(time in seconds); ylabel(ha(t)); axis([0,max(t)+0.05,min(ha),max(ha)+0.025]); grid %巴特沃兹模拟滤波器设计子程序 function[b,a]=afd_butt(Wp,Ws,Rp,As); if Wp=0 error(Passband edge must be larger than 0) end if Ws=Wp error(Stopband edge must be larger than Passed edge) end if (Rp=0)|(As0) error(PB ri