十四、IIRDF双线性变换法设计和滤波.docxVIP

实验题目IIRDF的双线性变换法设计与滤波　实验时间　实验台号　实验目的双线性变换法的原理；滤波的实现；　实验仪器PC台式机，软件MATLAB实验预习实验预习：1、IIR滤波器结构；2、滤波器指标的含义；实验原理：为了克服冲激响应不变法多值映射这一缺点，双线性变换法首先将整个s平面映射到平面中一个带宽为的横带上，然后通过变换关系把平面映射到z平面上。1、系统函数：2、模拟角频和数字角频的关系为： 3、设计方法：方法一：①按一定规律将数字滤波器指标转换为模拟滤波器指标；②根据转换后的指标求出模拟滤波器的阶数和截止频率；③利用阶数设计出模拟滤波器函数；④利用冲激响应不变法函数bilinear（）把模拟滤波器变换成数字滤波器；方法二：①根据数字滤波器指标直接调用设计函数求数字滤波器阶数；②根据阶数，调用直接设计数字滤波器函数获得系统函数分子、分母系数。（注，系统函数直接射时采用的均是双线性变换法）实验步骤应用函数bilinear();buttord();butter()；题目14-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滤波器); fc=1000;ap=1;as=25;fp=100;fs=300;wp=2*pi*fp/fc;ws=2*pi*fs/fc;anp=2*fc*tan(wp/2);ans=2*fc*tan(ws/2);[N,anc]=buttord(anp,ans,ap,as,s);[b,a]=butter(N,anc,s);[B2,A2]=bilinear(b,a,fc)[H2,w]=freqz(B2,A2,whole);subplot(2,1,2);plot(w*fc/2/pi,20*log10(abs(H2)));grid on;axis([0,1000,-100,0]);ylabel(H2幅值dB);title(双线性变换法设计的数字低通IIR滤波器); %结果14-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; %结果14-2实验数据结果及分析14-114-2思考题教师评语　　成绩:　日期: