利用Matlab实现Butterworthchebyshev滤波器的模拟.docVIP

数字信号处理实验 科目 数字信号处理基础 实验项目 学生 学号 指导老师 实验日期 一：实验目的： 掌握Matlab的操作； 学会利用Matlab实现Butterworth chebyshev滤波器的模拟； 二：实验原理： Butterworth chebyshev模拟滤波器的理论知识和设计原理； 冲激响应不变法 双线性变换法的理论知识 三：实验内容： 设计满足下列指标的数字低通滤波器（6个图） Wp=0.2π Rp=12 dB; Ws=0.5π Rs=20 dB; fs=1 KHz; 利用Butterworth chebyshevI型设计出模拟低通，采用冲激响应不变法、双线性变换法转换成。 直接设计出Butterworth chebyshevI型LPDF 一、 1、双线性Butterworth Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T; Lp=2*Fs*tan(Wp/2); Ls=2*Fs*tan(Ws/2); [N,Ln]= buttord(Lp,Ls,Rp,Rs,s); [B,A]=butter(N,Ln,s); [b,a]=bilinear(B,A,Fs); [h,w]=freqz(b,a,256); h1=20*log10(abs(h)); plot(w/pi,h1); grid; xlabel(Digital Frequency in pi units); ylabel(Gain in DB); axis([0 1 -50 10]); 2、双线性Chebyshev Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T; Lp=2*Fs*tan(Wp/2); Ls=2*Fs*tan(Ws/2); [N,Wn]= cheb1ord(Lp,Ls,Rp,Rs,s) [B,A]=cheby1(N,Rp,Wn,s); [b,a]=bilinear(B,A,Fs); [h,w]=freqz(b,a,256); h1=20*log10(abs(h)); plot(w/pi,h1); grid; xlabel(Digital Frequency in pi units); ylabel(Gain in DB); axis([0 1 -50 10]); 3、冲激响应Butterworth Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T; Lp=Wp*Fs; Ls=Ws*Fs; [N,Ln]= buttord(Lp,Ls,Rp,Rs,s); [B,A]=butter(N,Ln,s); [b,a]=impinvar(B,A,Fs); [h,w]=freqz(b,a,256); h1=20*log10(abs(h)); plot(w/pi,h1); grid; xlabel(Digital Frequency in pi units); ylabel(Gain in DB); axis([0 1 -50 10]); 4、冲激响应Chebyshev Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T; Lp=Wp*Fs; Ls=Ws*Fs; [N,Ln]= cheb1ord(Lp,Ls,Rp,Rs,s); [B,A]=cheby1(N,Rp,Wn,s); [b,a]=impinvar(B,A,Fs); [h,w]=freqz(b,a,256); h1=20*log10(abs(h)); plot(w/pi,h1); grid; xlabel(Digital Frequency in pi units); ylabel(Gain in DB); axis([0 1 -50 10]); 二、 1、直接设计Butterworth Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T; Lp=Wp/pi; Ls=Ws/pi; [N,Ln]= buttord(Lp,Ls,Rp,Rs); [B,A]=butter(N,Ln); [h,w]=freqz(B,A,256); h=20*log10(abs(h)); plot(w/pi,h); grid; xlabel(Digital Frequency in pi units); ylabel(Gain in DB); axis([0 1 -50 10]); 2、直接设计Chebyshev Wp=0.2*pi; Ws=0.5*pi; Rp=1; Rs=20; T=0.001; Fs=1/T;