数字滤波器实验.docVIP

实验二 IIR数字滤波器的设计实验报告 实验目的? 掌握双线性变换法及冲激响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 2）观察双线性变换及冲激响应不变法设计的滤波器的频域特性，了解双线性变换法及冲激响应不变法的特点。 3）熟悉Butterworth滤波器、Chebyshev滤波器的频率特性。 二、实验内容 1、fp=0.2KHz,Ap=1dB, fs=0.3KHz,As=25dB,T=1ms；分别用冲激响应不变法及双线性变换法设计一Butterworth数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 首先使用matlab tool箱的FDA进行低通IIR滤波器设计 使用FDA进行仿真如图： Butterworth低通数字滤波器： 使用脉冲响应不变法设计Butterworth低通数字滤波器： fp=200; fs=300; Fs=1000; Rp=1; Rs=25; T=1/Fs; %设计指标 W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率 [N,Wn]=buttord(W1p,W1s,Rp,Rs,s); %确定butterworth的最小介数N和频率参数Wn [z,p,k]=buttap(N); %设计模拟低通原型的零极点增益参数 [bp,ap]=zp2tf(z,p,k); %将零极点增益转换成分子分母参数 [bs,as]=lp2lp(bp,ap,Wn*pi*Fs);%将低通原型转换为模拟低通 [bz,az]=impinvar(bs,as,Fs); %用脉冲响应不变法进行模数变换 sys=tf(bz,az,T); %给出传输函数H(Z) [H,W]=freqz(bz,az,512,Fs); %生成频率响应参数 subplot(2,1,1); plot(W,20*log10(abs(H))); %绘制幅频响应 grid on; %加坐标网格 xlabel(频率/Hz); ylabel(振幅/dB); subplot(2,1,2); plot(W,abs(H)); grid on; xlabel(频率/Hz); ylabel(振幅/H);  Fp=200hz的衰减0.7732db；Fp=300hz的衰减25.0043；IIR滤波器的节数N=9 脉冲响应不变法得到的Butterworth低通数字滤波器满足要求； 使用双线性变换法设计Butterworth低通数字滤波器：  fs=1000; wp=2*pi*200/fs; ws=2*pi*300/fs; Rp=1; Rs=25; Ts=1/fs; Wp=2/Ts*tan(wp/2);Ws=2/Ts*tan(ws/2); %按频率转换公式进行转换 [N,Wn]=buttord(Wp,Ws,Rp,Rs,s); %计算模拟滤波器的最小阶数 [z,p,k]=buttap(N);%设计模拟原型滤波器 [Bap,Aap]=zp2tf(z,p,k); %零点极点增益形式转换为传递函数形式 [b,a]=lp2lp(Bap,Aap,Wn); %低通转换为低通滤波器的频率转化 [bz,az]=bilinear(b,a,fs); %运用双线性变换法得到数字滤波器传递函数 [H,f]=freqz(bz,az,500,fs) subplot(2,1,1); plot(f,20*log10(abs(H))); title(频率响应); grid on; xlabel(频率/Hz); ylabel(振幅/dB); subplot(2,1,2); plot(f,abs(H)); grid on; xlabel(频率/Hz); ylabel(振幅/H); -20*log10(abs(H(201))) -20*log10(abs(H(301))) N  Fp=200hz的衰减0.5977db；Fp=300hz的衰减25.0000；IIR滤波器的节数N=6 双线性变换法得到的Butterworth低通数字滤波器满足要求； 脉冲响应不变法的优点：能够模拟滤波器保持线性相位关系， 缺点：BW型模拟低通滤波器不是带限的，S-Z的变换不是一对一的对应关系脉冲响应不变法设计的滤波器存