巴特沃斯高通数字滤波器1.41.351.31.25Magnitude1.21.151.11.05100.10.20.30.40.50.6Frequency/Hz0.70.80.91
3.3 用脉冲响应不变法设计IIR带通数字滤波器实例。
带通滤波器技术指标为：Wp=[0.25π 0.35π]rad Ap=1dB
T=2s Ws=[ 0.15π 0.40π]rad As=10dB
程序为; T=2;
%采样周期
fs=1T;
%采样频率
Wp=[0.25.*pi 0.35.*pi];
Ws=[0.15.*pi 0.4.*pi]; %设置归一化通带和阻带截止平率
16
Ap= 1; As= 10;
%设置通带最大最小衰减
[N,Wc]=buttord(Wp,Ws,Ap,As,'s'); %调用butter函数确定巴特沃斯滤波器阶数
[B,A]=butter(N,Wc,'bandpass','s'); %调用butter函数设计巴特沃斯滤波器
W=linspace(0,pi,400*pi); %指定一段频率值
[D,C]=impinvar(B,A,fs); %调用脉冲响应不变法¨
Hz=freqz(D,C,W); %·返回频率值| plot(Wpi,abs(Hz)); %绘出巴特沃斯数字滤波器的幅频特性曲线 grid on;
title('巴特沃斯带通滤波器'); xlabel('FrequencyHz'); ylabel('Magnitude'); 得出幅频特性如下图：
17
巴特沃斯带通数字滤波器1.41.21Magnitude0.80.60.40.2000.10.20.30.40.50.6Frequency/Hz0.70.80.91
3.4 用脉冲响应不变法设计IIR带阻数字滤波器实例。
带阻滤波器技术指标为：Wp=[0.15π 0.40π]rad Ap=1dB
T=2s Ws=[ 0.25π 0.35π]rad As=11dB 程序为： T=2;
fs=1T;
Wp=[0.15.*pi 0.40.*pi];
Ws=[0.25.*pi 0.35.*pi]; %设置归一化通带和阻带截止平率 Ap= 1; As= 11;
%设置通带最大最小衰减
18
[N,Wc]=buttord(Wp,Ws,Ap,As,'s'); %调用butter函数确定巴特沃斯滤波器阶数 [B,A]=butter(N,Wc,'stop','s'); %调用butter函数设计巴特沃斯滤波器÷
W=linspace(0,pi,400*pi); %指定一段频率值
[D,C]=impinvar(B,A,fs); %调用脉冲响应不变法·¨
Hz=freqz(D,C,W); %·返回频率响应 plot(Wpi,abs(Hz)); %绘出巴特沃斯数字滤波器的幅频特性曲线 grid on;
title('巴特沃斯带阻滤波器'); xlabel('FrequencyHz'); 得出幅频特性如下图：
19
