湖南大学数字信号处理实验三 - 图文

信息科学与工程学院 数字信号处理实验

[N,Wc]=cheb1ord(wp,ws,Ap,As,'s'); [b,a]=cheby1(N,Ap,Wc,'s');

[B,A]=impinvar(b,a,fs); %脉冲响应不变法 freqz(B,A,'whole'); title('Chebyshev1'); figure %Chebyshev2

[N,Wc]=cheb2ord(wp,ws,Ap,As,'s'); [b,a]=cheby2(N,Ap,Wc,'s');

[B,A]=impinvar(b,a,fs); %脉冲响应不变法 freqz(B,A,'whole'); title('Chebyshev2'); figure %Elliptic

[N,Wc]=ellipord(wp,ws,Ap,As,'s'); [b,a]=ellip(N,Ap,As,Wc,'s');

[B,A]=impinvar(b,a,fs); %脉冲响应不变法 freqz(B,A,'whole'); title('Elliptic');

－ 29 －

信息科学与工程学院 数字信号处理实验

－ 30 －

信息科学与工程学院 数字信号处理实验

－ 31 －

信息科学与工程学院 数字信号处理实验

双线性变换法

%'Butter low-pass filter' %数字低通滤波器技术指标 fs=10000; T=1/fs;

vp=2*pi*1000/fs; vs=2*pi*1400/fs; Ap=2; As=20;

%模拟低通滤波器技术指标 wp=(2/T)*tan(vp/2); ws=(2/T)*tan(vs/2);

[N,Wc]=buttord(wp,ws,Ap,As,'s'); [b,a]=butter(N,Wc,'s'); w=linspace(100,1500,500)*2*pi; [B,A]=bilinear(b,a,fs);%双线性变换法 freqz(B,A,'whole');

title('Butter low-pass filter'); figure %Chebyshev1

－ 32 －