实验二. R数字滤波器的设计.docVIP

实验二 IIR数字滤波器的设计 实验目的 学会用MATLAB设计IIR数字滤波器。 实验原理 IIR数字滤波器的设计 （1）脉冲响应不变法变换原理 脉冲响应不变法将模拟滤波器的s平面转换成数字滤波器的z平面，从而将模拟滤波器映射成数字滤波器。IIR滤波器的系统函数为(或z)的有理式，即 ，且一般满足MN （2）双线性变换法变换原理 为克服脉冲响应不变法产生频率响应的混叠失真，可以采用非线性频率压缩方法，使s平面与z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，这就是双线性变换法。 实验内容 问题〔.1〕试设计一个模拟低通滤波器，fp=2400Hz，fs=5000?Hz，Rp=3 dB，Rs=25?dB。分别用巴特沃斯和椭圆滤波器原型，求出其3dB截止频率和滤波器阶数，传递函数，并作出幅频、相频特性曲线。巴特沃斯滤波器s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel(频率 Hz);ylabel(幅度 dB); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel(频率 Hz);ylabel(相角 rad); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 程序（椭圆滤波器 f_p=2400; f_s=5000; R_p=3; R_s=25; % 设计要求指标 [n, fn]=ellipord(f_p,f_s,R_p,R_s,s); % 计算阶数和截止频率 Wn=2*pi*fn; % 转换为角频率 [b,a]=ellip(n,R_p,R_s,Wn,s); % 计算H(s) f=0:100:10000; % 计算频率点和频率范围 s=j*2*pi*f; % s=jw=j*2*pi*f H_s=polyval(b,s)./polyval(a,s); % 计算相应频率点处H(s)的值 figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel(频率 Hz);ylabel(幅度 dB); subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabel(频率 Hz);ylabel(相角 rad); figure(2); freqs(b,a); % 也可用指令freqs直接画出H(s)的频率响应曲线。 问题〔.2〕试设计一个巴特沃斯型数字低通滤波器，设采样率为8000Hz，fp=2100Hz，fs=2500?Hz，Rp=3dB，Rs=25?dB。figure(1); freqz(b,a, 1000, 8000) % 作出H(z)的幅频相频图, freqz(b,a, 计算点数, 采样率) subplot(2,1,1); axis([0 4000 -30 3]) figure(2); % 第二种作图方法 f=0:40:4000; % 计算频率点和频率范围 z=exp(j*2*pi*f./(f_N)); % H_z=polyval(b,z)./polyval(a,z); % 计算相应频率点处H(s)的值 subplot(2,1,1); plot(f, 20*log10(abs(H_z))); % 幅频特性 axis([0 4000 -40 1]); xlabel(频率 Hz);ylabel(幅度 dB); subplot(2,1,2); plot(f, angle(H_z)); % 相频特性 xlabel(频率 Hz);ylabel(相角 rad); 实例〔.3〕试设计一个切比雪夫1型高通数字滤波器，采样率为8000Hz，fp=1000Hz，