设计一个数字带通滤波器

西安邮电大学

数字信号处理课内上机报告

专业班级： 通工1002班 学生姓名： 王敏 学号(班内序号): 03101069（31号）

2013 年 6 月 20 日

——————————————————————————

装

订

线————————————————————————————————

报告份数：

实验内容：

设计一个数字带通滤波器，通带范围为0.3πrad到0.4πrad，通带最大衰减为αp = 3dB，阻带最小衰减为αs = 18dB;0.2π一下和0.5π以上范围为阻带；1）采用巴特沃斯性，分别利用

A）脉冲响应不变法；

B）双线性变换法；

设计之。

脉冲响应不变法设计：

Matlab 代码：

Wph = 0.4*pi;Wpl = 0.3*pi;Wsh = 0.5*pi;Wsl = 0.2*pi;Ap = 3;As = 18;%给定滤波器参数

Wp = [Wpl Wph];Ws = [Wsl Wsh];%通带范围，阻带范围

[N,Wn] = buttord(Wp,Ws,Ap,As,'s');%计算巴特沃斯模拟带通滤波器的阶数

[Bt,At] = butter(N,Wn,'s');%设计巴特沃斯模拟滤波器

[Bz,Az] = bilinear(Bt,At,1);%双线性变换为数字带通滤波器

%[Bz,Az] = impinvar(Bt,At,1);%脉冲响应不变法变换为数字带通滤波器

[H,W] = freqz(Bz,Az);%求滤波器的频率响应

figure(1);

subplot(2,1,1);

plot(W/pi,abs(H));%绘制出数字带通滤波器的幅频幅度特性

title('数字带通滤波器幅频响应');

xlabel('归一化频率');ylabel('幅值');grid on;

subplot(2,1,2);

plot(W/pi,angle(H));%绘制出数字带通滤波器的幅频相位特性

title('数字带通滤波器幅频响应');

xlabel('归一化频率');ylabel('相位');grid on;

figure(2);

plot(W/pi,20*log(abs(H)));%绘制出数字带通滤波器的幅频衰减的性title('数字带通滤波器幅频特性（dB）');

xlabel('归一化频率');ylabel('幅频特性（dB）');grid on;

脉冲响应不变法：

双线性变换法：

指导教师评语：

实 验 成 绩: 指导(辅导)教师 :

—————————————————————————

装

订

线————————————————————————————————