实验报告3滤波器.docVIP

大连理工大学实验报告 学院（系）： 信息与通信工程 专业： 通信工程 班级： 电通1101 姓 名： 殷青 学号： 201181227 实验时间： 2014.5.21 实验室： 创C221 指导教师签字： 成绩： 实验三 IIR FIR数字滤波器设计 实验目的和要求 （见预习报告） 主要仪器设备 MATLAB 实验步骤与操作方法 （见预习报告） 四、实验结果与分析 IIR数字滤波器设计 用buttord和butter函数，直接，通带内衰减不大于1dB，阻带起始频率为，阻带内衰减不小于15dB，观察其频谱响应的特点。 MATLAB源码： %butttordbutter design 巴特沃兹高通滤波器 Wp=0.3;Ws=0.2;Rp=1;Rs=15; [N,Wn]=buttord(Wp,Ws,Rp,Rs); [Bz,Az]=butter(N,Wn,high); figure; freqz(Bz,Az,512); % Hf=abs(H); % Hx=angle(H); % clf % figure(1) % plot(w,Hf) % title(幅频特性曲线) % figure(2) % plot(w,Hx) % title(离散系统相频特性曲线) 实验结果：   给定带通滤波器的技术指标：通带上下截止频率为，通带内衰减不大于3dB，阻带上下起始频率为，阻带内衰减不小于18dB。用buttord和butter函数，巴特沃斯 Wp1=0.3*pi; Wp2=0.4*pi; Ws1=0.2*pi; Ws2=0.5*pi; Ap=3; As=18; %Fs=1000; Wp=[Wp1 Wp2];Ws=[Ws1 Ws2]; [n,Wn]=buttord(Wp/pi,Ws/pi,Ap,As); [b,a]=butter(n,Wn,bandpass);%直接设计滤波器 freqz(b,a,512,1000); （3）切比雪夫 %切比雪芙滤波器 Wp1=0.3*pi; Wp2=0.4*pi; Ws1=0.2*pi; Ws2=0.5*pi; Ap=3; As=18; %Fs=1000; %Chebyshew iLPF设计 figure; Wp=[Wp1 Wp2];Ws=[Ws1 Ws2]; [n,Wn]=cheb1ord(Wp/pi,Ws/pi,Ap,As); [b,a]=cheby1(n,Ap,Wn,bandpass);%直接设计滤波器 freqz(b,a,512,1000); %%Chebyshew LPF设计 figure; [n1,Wn1]=cheb2ord(Wp/pi,Ws/pi,Ap,As); [b1,a1]=cheby2(n1,Ap,Wn1,bandpass);%直接设计滤波器 freqz(b1,a1,512,10000);  实验结果： 巴特沃斯 切比雪夫I型 切比雪夫II型 分析：对比巴特沃斯和切比雪夫I型，可以发现巴特沃斯滤波器在通带内幅度响应平坦，而切比雪夫I型在通带内有波纹，切比雪夫II型在阻带有波纹。 3.用双线性变换法模拟滤波器原型，抽样频率为 。 MATLAB源码： %双线性法的模拟滤波器原型设计巴特沃斯低通滤波器 fp=100;fst=300;Ap=3;As=20; Fs=1000;Ts=1/Fs; wp=2*pi*fp/Fs;ws=2*pi*fst/Fs; Bp=2*tan(wp/2)/Ts; %模拟滤波器指标 Bs=2*tan(ws/2)/Ts; [N,Wn]=buttord(Bp,Bs,Ap,As,s); [z,p,k]=buttap(N); %设计模拟低通原型滤波器 [b,a]=zp2tf(z,p,k); %把滤波器零极点模型转化为传递函数模型 [B,A]=lp2lp(b,a,Wn); [Bz,Az]=bilinear(B,A,Fs); figure; [H,w]=freqz(Bz,Az,512); % Hf=abs(H); % Hx=angle(H); % clf % figure(1) % plot(w,Hf) % title(幅频特性曲线) % figure(2) % plot(w,Hx) % title(离散系统相频特性曲线) 实验结果：  4.用双线性变换法模拟滤波器原型和设计法（butter）两种方法，设计Fs=2000Hz，试设计一个为此系统使用的带通数字滤波器。要求： （1）通带范围为300～400Hz，在