设计一个数字带通滤波器.doc

实验总成绩：
——————————————————————————装 订 线—————————————————————————————————

报告份数：

西安邮电大学

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

专业班级： 通工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;
脉冲响应不变法：

双线性变换法：

——————————————————————————装 订 线—————————————————————————————————
指导教师评语：

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

- 5 -