数字滤波器的设计实验VIP专享

实验二 IIR 数字滤波器的设计 实验内容及步骤: 数字滤波器的性能指标：通带临界频率fp、阻带临界频率fr；通带内的最大衰减 Ap；阻带内的最小衰减 Ar；采样周期 T； (1)、fp=0.3KHz,Ap=0.8dB, fr=0.2KHz,Ar=20dB,T=1ms；设计一 Chebyshev 高通滤波器；观察其通带损耗和阻带衰减是否满足要求。 程序如下： fp=300; fr=200; Ap=0.8; Ar=20; T=0.001;fs=1/T; wp=2*pi*fp*T; wr=2*pi*fr*T; Wp=2/T*tan(wp/2); Wr=2/T*tan(wr/2); [N,Wn]=cheb1ord(Wp,Wr,Ap,Ar,'s'); [B,A] = cheby1(N,Ap,Wn,'high','s'); [num,den]=bilinear(B,A,1/T); [h,w]=freqz(num,den); plot(w*fs/(2*pi),20*log10(abs(h))); %衰减及频率都用归一化的1为单位显示 axis([0,500,-30,0]); title('Chebyshev高通滤波器'); xlabel('频率'); ylabel('衰减'); grid on; 根据下图知道通带损耗与阻带衰减满足要求 (2)、fp=0.2KHz,Ap=1dB, fr=0.3KHz,Ar=25dB,T=1ms；分别用脉冲响应不变法及双线性变换法设计一Butterworth 数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 程序如下： fp=200; fr=300; Ap=1;Ar=25; T=0.001;fs=1/T; wp=2*pi*fp*T; wr=2*pi*fr*T; Wp=2/T*tan(wp/2); Wr=2/T*tan(wr/2); [N,Wn]=buttord(Wp,Wr,Ap,Ar,'s'); [B,A] = butter(N,Wn,'s'); [num1,den1]=impinvar(B,A,1/T); %脉冲响应不变法得出设计的传递函数 [num2,den2]=bilinear(B,A,1/T); %双线性变换法得出设计的传递函数 [h1,w]=freqz(num1,den1); plot(w*fs/(2*pi),20*log10(abs(h2)),w*fs/(2*pi),20*log10(abs(h1)),'r.');grid on; %衰减及频率都用归一化的1为单位显示 axis([0,500,-30,0]); title('Butterworth低通滤波器（红线—脉冲响应不变法 蓝线—双线性变换法）'); xlabel('ƵÂÊ'); ylabel('Ë¥¼õ'); grid on; 优缺点： 采用脉冲响应不变法 优点： 1.h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好 2 线性相位模拟滤波器转变为线性相位数字滤波器 缺点： 1.对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大 2 不能用来设计高通和带阻滤波器。只适用于限带的低通、带通滤波器 采用双线性变换法 优点：1 避免了频率响应的混迭失真现象 2 在特定数字滤波器和特定模拟...