巴特沃斯数字带通滤波器

1、巴特沃斯数字带通滤波器 数字信号处理课程设计报告设计课题滤波器设计与实现专业班级姓名学号报告日期2012年12月目录1. 课题描述2. 设计原理2.1 滤波器的分类2.2 数字滤波器的设计指标2.3 巴特沃斯数字带阻模拟滤波器2.3.1巴特沃斯数字带通滤波器的设计原理2.3.2巴特沃斯数字带通滤波器的设计步骤3. 设计内容3.1 用MATLAB编程实现3.2 设计结果分析4总结5.参考文献课程设计任务书题目滤波器设计与实现学生姓名设计內容与要求一、设计内容：设计巴特沃斯数字带通滤波器，通带频率200500hz，阻带上限频率600h乙阻带卜限频率150hz，通带衰减取大0.5dB，阻带最小衰减4

2、0dB，采样频率2000hz，画出幅频、相频响应曲 线，并设计信号验证滤波器设计的正确性。二、设计要求1设计报告一律按照规定的格式，使用A4纸，格式、封面统一给出 模版。2报告内容(1)设计题目及要求(2 )设计原理(包括滤波器工作原理、涉及到的MATLAB函数的 说明)(3 )设计内容(设计思路，设计流程、仿真结果)(4)设计总结(收获和体会)(5 )参考文献(6 )程序清单祀【卜时间2012年12月3日至2011年12月11日指导教师签名2011 年 12 月 2日系(教研室)主任 签名年月日学生签名年月日1 课题描述数字滤波器是由数字乘法器、加法器和延时单元组成的一种算法或 装置。 数字

3、滤波器的功能是对输入离散信号的数字代码进行运算处 理，以达到改变信 号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波器 已可用计算机软件实现，也可用大规模集成数字硬件实时实现。使用 MATLAB信号处理箱和BW （巴特沃斯）设 计低通数字滤波器。2.设计原理2.1滤波器的分类 数字滤波器有低通、高通、带通、带阻和全通等类型。它 可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。如果数 字滤波器的内部参数不随时间而变化，则称为时不变的，否则为时变的。如 果数字滤波器在某一给定时刻的响应与在此时刻以后的激励无尖，则称为因果 的，否则为非因果的。如果数字滤波器对单一或多个激

4、励信号的响应满足线 性条件，则称为线性的，否则为非线性的。应用最广的是线性、时不变数字 滤波器。数字滤波器也可以按所处理信号的维数分为一维、二维或多维数字滤波 器。一维数字滤波器处理的信号为单变量函数序列，例如时间函数的抽样 值。二维或多维数字滤波器处理的信号为两个或多个变量 函数序列。2.2 数字带通滤波器的设计指标阻带频率150600hz 通带上限频率500,通带下限频率200hz，通带衰 减为0.5dB，阻带最大衰减40dB 采样频率2000hz2.21巴特沃斯原理实际的滤波电路往往难以达到理想的要求，如要同时在幅频和相频响应两方 面都满足要求就更为困难。因此，只有根据不同的实际需要，寻

5、求最佳的近似 理想特性。例如，可以主要着眼于幅频响应，而不考虑相频响应；也可以从满 足相频响应出发，而把幅频响应居于次要位置。介绍一种最简单也是最常用的 滤波电路一一巴特沃斯滤波电路(又叫最平幅度滤波电路)。这种 滤波电路对幅 频响应的要求是：在小于截止频率CCO的范围内，具有最平幅度的响应，而在 CCOCO 后，幅频响应迅速下降。2.3.2巴特沃斯带通数字滤波器的设计步骤(1) 确定滤波器的指标即：通带上截止频率。通带下截止频率。阻带上截止频率，阻带下截止频率。以及通带内最大衰减和阻带最小衰减。(2) 求出模拟带通滤波器指标(3) 模拟归一化低通滤波器技术指标(4) 设计模拟低通滤波器(5)

6、 将归一化模拟低通妆化为模拟带通(6) 利用双线性变换法将Ha(s)转化为数字带通滤波器H(Z).(7) 作图显示滤波器的幅频特性和相位特性。H,W = freqz(B 乙 Az); subplot(2,2,1);plot(W*ft/(2*pi),abs(H);gridon;xlabel(1 频率);ylabel(，幅值 *);3设计内容3.1用MATLAB编程实现ft=2000;fpl=150; fph=600;wp1= fpl *2*pi;%临界频率采用模拟角频率表 示wph= fph*2*pi; wp= wp1 ,wph;%临界频率采用模拟角频率表 示wpb=wp/ ft;%求数字频率r

7、p=0.5; rs=40; fsl=200;fsh=500; ws1= fsl *2*pi;%临界频率米用模拟角频率表不wsh= fsh *2*pi; ws= ws1, wsh;%临界频率采用模拟角频率表示wsb=ws/ ft; 率%求数字频OmegaP=2* ft *tan(wpb/2);OmegaS=2* ft*tan(wsb/2);%选择滤波器的最小阶数N,Wn=buttord(OmegaP,OmegaS, rp, rs,变后获得 0 ); %此处是代入经预畸 的归一化模拟频率参数Bt,At=butter(N,Wn, );% 设计一个 N 阶的巴特沃 思模拟滤波器Bz,Az=biline

8、ar(Bt,At, ft);%双线性变换为数字滤波器%求解数字滤波器的频率响应title(1数字滤波器幅频响应)；subplot(2,2,2); plot(W*ft/(2*pi),angle(H);grid on; xlabel( * 频率/Hz );ylabel( 相位 *); t=0:50;y=sin(50 t+1/3 pi)+sin(300 t+pi); subplot(2,2,3);plot(t,y);grid on;y1 =filter(Bz,Az,y); subplot(2,2,4); plot(t,y1);grid on;仙：血 Tcdi 上$刪HrpX礙1（，谶鑼辩做斤设计疇盧沃斯带阻滤波器时，由于通带为一定的 频率段。因此被滤掉的部分为平缓先线条。4总结本次课程实验中，让我慢慢了解了程序的编译，对于程序的编译也渐渐地掌握了其中的方法与规律。通过这次试验让我收获颇大，在 实验设计中，要先确定需要的函数，然后根据所给条