基于Matlab的IIR数字滤波器设计(脉冲响应不变法)毕业论文

巴特沃斯高通数字滤波器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