信息科学与工程学院 数字信号处理实验
[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 -
相关推荐: