数字信号处理第二次实验报告..docVIP

数字信号处理实验报告第二次实验 ——IIR数字滤波器的设计 姓名:陈桐 学号 实验日期:2006年11月14日 实验目的: (1) 掌握双线形变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线形变换法及脉冲响应不变法设计低通，高通和带通IIR数字滤波器的计算机编程。 (2) 观察双线形变换法及脉冲响应不变法设计的滤波器的频域特性，了解双线形变换法及脉冲响应不变法的特点。 (3) 熟悉巴特沃思滤波器，切比雪夫滤波器和椭圆滤波器的频率特性。 实验原理: 实验中有关变量的定义：fc通带边界频率，fr阻带边界频率，tao通带波动，at最小阻带衰减，fs采样频率，t采样周期。 设计一个数字滤波器一般包括以下两步： a.按照任务要求，确定滤波器性能指标 b.用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求 数字滤波器的实现： 对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。 用双线形变换法设计IIR数字滤波器的过程： a.将设计性能指标中的关键频率点进行“预畸” b.利用“预畸”得到的频率点设计一个模拟滤波器。 c.双线形变换，确定系统函数 实验内容: fc=0.3kHz,δ=0.8dB,fr=0.2kHz,At=20dB,T=1ms;设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。 程序: clc; clear; figure(2); %预畸变模拟边界频率 wc=2*1000*tan(2*pi*300/(2*1000)); wr=2*1000*tan(2*pi*200/(2*1000)); %阻带单调，通带波纹，方法为切比雪夫 [N,wn]=cheb1ord(wc,wr,0.8,20,s); %高通 [B,A]=cheby1(N,0.8,wn,high,s); [num,den]=bilinear(B,A,1000); %频谱 [h,w]=freqz(num,den); f=w/pi*500; %画图 plot(f,20*log10(abs(h))); axis([0,500,-80,10]); grid on; xlabel(频率/Hz); ylabel(幅度/dB); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 图形为: 由图易得，该滤波器完全满足要求。δ趋近于0，f=200Hz时，幅度约为-30dB，满足At=20dB的要求。 fc=0.2kHz, δ=1dB,fr=0.3kHz,At=25dB,T=1ms;分别用脉冲响应不变法及双线形变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 程序为: clc; clear; fs=1000; wp=2*pi*200;ws=2*pi*300 [n,wn]=buttord(wp,ws,1,25,s); % 设计Butterworth低通模拟滤波器 % 通带衰减小于1分贝、阻带（0.6pi）衰减大于25分贝 % n: 阶数 % wn:低通模拟滤波器的3分贝通带截止频率 [B1,A1] = BUTTER(n,wn,s); % 设计Butterworth低通模拟滤波器 % B,A：模拟滤波器传递函数分子、分母多项式系数 [b1,a1]=impinvar(B1,A1,fs); % 用脉冲响应不变法将模拟低通滤波器转换为数字滤波器b,a：分子、分母多项式系数 % 双曲变换的预畸变 [B2,A2]=butter(n,2000*tan(pi/5),s); % 用双曲变换法将模拟低通滤波器转化为数字滤波器b,a:分子,分母多项式系数 [b2,a2]=bilinear(B2,A2,fs); [h1,w1]=freqz(b1,a1,100,whole); % h: 在单位圆上采样100点的频率响应 % w：频率矢量 c1=20*log10(abs(h1)+0.0001); [h2,w2]=freqz(b2,a2,100,whole); % h: 在单位圆上采样100点的频率响应 % w：频率矢量 c2=20*log10(abs(h2)+0.0001); plot(w1(1:50), c1(1:50),-.); grid on; hold on; plot(w2(1:50), c2(1:50),-); grid on; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%