十三、IIRDF冲激响应不变法设计和滤波.docx

实验题目IIRDF的冲激响应不变法设计与滤波　实验时间　实验台号　实验目的冲激响应不变法的原理；滤波的实现；　实验仪器PC台式机，软件MATLAB实验预习实验预习：1、IIR滤波器结构；2、滤波器指标的含义；　实验原理：1、系统函数：2、冲激响应不变法原理：将模拟滤波器的单位冲激响应加以等时间间隔抽样使得。将进行拉普拉斯逆变换得到 ，再将进行抽样得到，将进行Z变换得到。其中： 3、模拟角频和数字角频的关系为： 4、设计方法：方法一：①按一定规律将数字滤波器指标转换为模拟滤波器指标；②根据转换后的指标求出模拟滤波器的阶数和截止频率；③利用阶数设计出模拟滤波器函数；④利用冲激响应不变法函数impinvar（）把模拟滤波器变换成数字滤波器；方法二：①根据数字滤波器指标直接调用设计函数求数字滤波器阶数；②根据阶数，调用直接设计数字滤波器函数获得系统函数分子、分母系数。（注，系统函数直接射时采用的均是双线性变换法）实验步骤应用函数impinvar();buttord();butter()题目13-1：分别用冲激响应不变法和直接法设计一个巴特沃斯低通数字滤波器。该滤波器通带截止频率为，阻带截止频率，通带最大衰减，阻带最小衰减。抽样周期T=1ms？fc=1000;ap=1;as=25;fp=100;fs=300;wp=2*pi*fp/fc;ws=2*pi*fs/fc;Wanp=wp*fc;Wans=ws*fc;[N,Wanc]=buttord(Wanp,Wans,ap,as,s);[b,a]=butter(N,Wanc,s);[B1,A1]=impinvar(b,a,fc)[H1,w]=freqz(B1,A1,whole);plot(w*fc/2/pi,20*log10(abs(H1)));grid on;axis([0,1000,-40,0]);ylabel(H1幅值dB);title(冲激响应不变法设计的数字低通IIR滤波器); %结果13-1fc=1000;ap=1;as=25;fp=100;fs=300;wp=2*fp/fc;ws=2*fs/fc;[N,Wn]=buttord(wp,ws,ap,as);[B,A]=butter(N,Wn);[H,W]=freqz(B,A,whole);figure(2);subplot(2,1,1);plot(W*fc/2/pi,20*log10(abs(H)));grid on;axis([0,1000,-100,0]);ylabel(H幅值dB);title(直接设计的数字低通IIR滤波器);subplot(2,1,2);plot(W*fc/2/pi,angle(H));xlabel(频率(Hz));ylabel(相位（弧度）);grid on; %结果13-2实验数据结果及分析13-113-2思考题教师评语　　成绩:　日期: