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

实验二IIR数字滤波器的设计 一、 实验目的 学会用MATLAB设计IIR数字滤波器。 二、 实验原理 1IR数字滤波器的设计 （1）脉冲响应不变法变换原理 脉冲响应不变法将模拟滤波器的s平面转换成数字滤波器的z平面，从而将 容E ,且一般满足M5N模拟滤波器映射成数字滤波器。IIR滤波器的系统函数为z（或z）的有理式，即 容E ,且一般满足M5N H （乙）= k = l （2）双线性变换法变换原理 为克服脉冲响应不变法产生频率响应的混叠失真，可以釆用非线性频率压缩 方法，使s平面与z平面建立了一一对应的单值关系，消除了多值变换性，也就 消除了频谱混叠现象，这就是双线性变换法。 三、实验内容 问题（2.1）试设计一个模拟低通滤波器，fp=2400Hz, fs=5000 Hz, Rp=3 dB, Rs=25 dB。 分别川巴特沃斯和椭恻滤波器原型，求出其3dB截止频率和滤波器阶数，传递函数，并作 出幅频、相频特性曲线。 程序：（巴特沃斯滤波器） Lp=2400; Ls=5000; R_p=3; R_s=25; % 设计要求指标 %转换为角频率%计算H（s）%计算频率点和频率范围% s=jw=j*2*pi*f%计算相应频率点处 %转换为角频率 %计算H（s） %计算频率点和频率范围 % s=jw=j*2*pi*f %计算相应频率点处H⑸的值 fb,a]=butter(n, Wn, s); f=0:100:1(X)00; s=j*2*pi*f; H_s=polyval(b,s)./polyval(a,s); figure(l); subplot(2,1,1); plot(f, 20*log 10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); %相频特性xlabclC频率 Hz);ylabclC幅度 dB1); subplot(2,l,2); plot(f, angle(H_s)); %相频特性 xlabelC频率 Hzr);ylabelC相角 rad); figure(2); freqs(b,a); %也J用指令freqs直接1丽出H(s)的频率响应曲线。 ■200 .103Frequency (rad/$)000004— 4— (?g?p)力左 ■200 . 103 Frequency (rad/$) 00000 4— 4— (?g?p)力左 程序(椭圆滤波器) Lp=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⑸的值 figure(l); subplot(2,1,1); plot(f, 20*log 10(abs(H_s))); % 幅频特性 axis([0 10000 -40 1]); xlabel(频率 Hzjylabel(幅度 dB‘)； subplot(2,1,2); plot(f, angle(H_s)); % 相频特性 xlabe频率 Hz);ylabelC相角 rad); figure(2); freqs(b,a); %也可丿LI指令freqs直接画出H(s)的频率响应ill]线。 问题(2.2)试设计一个巴特沃斯型数字低通滤波髀，设采样率为8000Hz, fp=2100Hz, fs=2500 Hz, Rp=3dB, Rs=25 dBo 程序： f_N=8000; % 釆样率 f_p=2100; f_s=2500; R_p=3; R_s=25; % 设计要求指标 Ws=Ls/(fLN/2); Wp=Lp/(f_N/2); % 计算归一化频率 [n, Wn]=buttord(Wp,Ws,R_p,R_s); % 计算阶数和截止频率 figure(l); freqz(b,a, 1000, 8000) subplot(2丄 1); axis([0 4000 -30 3]) figure(2); f=0:40:4000; %作出H（z）的輛频相频图，freqz（b,a，计算点数，采样率）%第二种作图方法%计算频率点和频率范围z=e xp(j * %作出H（z）的輛频相频图，freqz（b,a，计算点数，采样率） %第二种作图方法 %计算频率点和频率范围 H_z=polyval(b,z)./polyval(a,z); %计算相应频率点处H⑸的值