巴特沃斯数字带阻滤波器

《数字信号处理》课程设计报告

设计课题滤波器设计与实现

专业班级

姓名

学号

报告日期 2012年12月

目录

1. 课题描述 (2)

2. 设计原理 (2)

2.1 滤波器的分类 (3)

2.2 数字滤波器的设计指标 (3)

2.3 巴特沃斯数字带阻模拟滤波器 (3)

2.3.1 巴特沃斯数字带阻滤波器的设计原理 (4)

2.3.2 巴特沃斯数字带阻滤波器的设计步骤 (7)

3. 设计内容 (8)

3.1 用MATLAB编程实现 (10)

3.2 设计结果分析 (10)

4. 总结 (10)

5. 参考文献 (11)

课程设计任务书

1 .课题描述

数字滤波器是由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。使用MATLAB信号处理箱和BW（巴特沃斯）设计低通数字滤波器。

2.设计原理

2.1 滤波器的分类

数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。如果数字滤波器的内部参数不随时间而变化，则称为时不变的，否则为时变的。如果数字滤波器在某一给定时刻的响应与在此时刻以后的激励无关，则称为因果的，否则为非因果的。如果数字滤波器对单一或多个激励信号的响应满足线性条件，则称为线性的，否则为非线性的。应用最广的是线性、时不变数字滤波器。

数字滤波器也可以按所处理信号的维数分为一维、二维或多维数字滤波器。一维数字滤波器处理的信号为单变量函数序列，例如时间函数的抽样值。二维或多维数字滤波器处理的信号为两个或多个变量函数序列。

2.2 数字带阻滤波器的设计指标

阻带频率200~500hz，通带上限频率600hz, 通带下限频率150hz，通带衰减为0.5dB，阻带最大衰减20dB，采样频率2000hz

2.3 巴特沃斯数字带阻模拟滤波器

2.3.1 巴特沃斯数字带阻滤波器的设计原理

巴特沃斯低通滤波器的幅度平方函数|Ha(jΩ)|2用下式表示：

|Ha(jΩ)|=

公式中，N称为滤波器的阶数。在Ω=0时，|Ha（jΩ）||=1;Ω=

Ωc时，|Ha(jΩ)|=1/2，Ωc是3dB截止频率。在Ω=Ωc附近，随Ω加大，幅度迅速下降。幅度特性与Ω和N的关，N越大，通带越平坦，过渡带越宽，过渡带与阻带幅度下降的速度越快，总是频响特性与理想低通滤波器的误差越小。

以S代替jΩ,将幅度平方函数|Ha（jΩ）|2写成s的函数：

Ha（S）Ha(-S)=

复变量S=α+jΩ,此公式表明幅度平方函数有2N个极点，极点Sk 用下公式表示：

Sk=

公式中K=0,1,2……，2N-1。2N个极点等间隔分布在半径为Ωc的圆上，间隔是/N rad。为形成因果稳定的滤波器，2N个极点

中只取S平面左半平面的N个极点构成Ha(S),而右半平面的N个极点构成Ha(-S)。Ha(S)的表达式为

Ha(S)=/ （S-Sk）

设N=3,极点有六个，它们分

S0= S1=-Ωc S2=

S3= S4=Ωc S5=

取S平面的左半平面的极点S0,S1,S2组成系统函数Ha(S),即

Ha(S)=

由于不同的技术指标对应的边界频率和滤波器幅频特性不同，为使设计公式和图表统一将频率归一化。巴特沃斯滤波器采用对3dB截止频率Ωc归一化，归一化后的系统函数为

令=，λ称为归一化频率，称为归一化复变量，这样的巴特沃斯滤波器的归一化低通原型系统函数为

公式中，称为归一化极点，用下公式表示，

k=0,1,…;,N-1

显然， Sk=Ωc*Pk

这样，只要根据技术指标求出阶数N,按照公式求出N个极点，再按得到归一化低通原型系统函数Ga(),如

果给定Ωc，再去归一化，即将，代入中，便可得到期望设计的系统函数Ha(S)。

所使用的MATLAB函数：

1.巴特沃斯滤波器阶数选择函数

（1）[N,Ωc]=buttord(Ωp，Ωs，αp，αs，‘s’)

作用：计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。

说明：Ωp，Ωs，Ωc均为实际模拟角频率。

模拟频率f：每秒经历多少个周期，单位Hz，即1/s，信号的真实频率，可用于模拟信号和数字信号；

模拟角频率Ω：每秒经历多少弧度，单位rad/s，通常只于模拟信号；

数字频率w：每个采样点间隔之间的弧度，单位rad，通常只用于数字信号。

关系：Ω=2pi*f；w = Ω*T=2pi*f/F。(F=1/Ts为采样频率，Ts为采样间隔)

2.完整巴特沃斯滤波器设计函数

（1）格式： [b，a]=butter（N，wc，‘ftype’）

作用：计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。

说明：调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值，一般是调用buttord格式（1）计算N和wc。系数b、a是按照z-1的升幂排列。

（2）格式：[B，A]=butter（N，Ωc，‘ftype’，‘s’）

作用：计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量。

说明：调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB 截止频率(实际角频率)，可调用buttord（2）格式计算N和Ωc。系数B、A按s的正降幂排列。

tfype为滤波器的类型：

◇ftype=high时，高通；Ωc只有1个值。

◇ftype=stop时，带阻；Ωc=[Ωcl,Ωcu]，分别为带阻滤波器的通带3dB下截止频率和上截止频率。

◇ ftype缺省时：若Ωc只有1个值，则默认为低通；若Ωc有2个值，则默认为带通；其通带频率区间Ωcl<Ω <Ωcu。

3. 求离散系统频响特性的函数freqz()

格式：[H,w]=freqz(b,a,N)

说明：b和a分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0～pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512。可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，绘制出系