数字信号处理课程设计题目巴特沃斯高通数字滤波器老师陈忠泽老师学院电气工程学院班级电子信息工程081班学号20084470110姓名何依阳二0一一年五月目录：一、IIR数字高通滤波器的设计1、数字滤波器的概述2、数字滤波器的设计步骤3、设计方法4、IIR巴特沃斯数字高通滤波器的实例计算二、软件仿真工具及实现环境简介1、计算机辅助设计方法2、MATLAB直接设计IIR巴特沃斯数字高通滤波器三、滤波器结构对数字滤波器性能指标的影响分析1、IIR系统的基本网络结构（1）直接型（2）级联型（3）并联型四、有限字长运算在网络结构中对数字滤波器的影响1、运算量化效应对数字滤波器的影响2、参数的字长对数字滤波器性能指标的影响2.1、系数量化对数字滤波器的影响五、运用MATLAB的辅助工具FDATOOL画出系统函数图像六、设计心得IIR数字高通滤波器的设计一、IIR数字高通滤波器的设计1、数字滤波器的概述所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。
2、数字滤波器的设计步骤设计一个IIR数字滤波器主要包括下面5个步骤：（1）确定滤波器要求的规范指标。
（2）选择合适的滤波器系数的计算（如图一流程图所示）。
（3）用一个适当的结构来表示滤波器（实现结构）。
（4）有限字长效应对滤波器性能的影响分析。
（5）用软件或硬件来实现滤波器。
确定数字巴特沃斯高通滤波器指标推导出归一化模拟巴特沃斯低通滤波器指标计算出归一化模拟巴特沃斯低通滤波器去归一化推导出模拟巴特沃斯高通滤波器双线性变换推导出数字巴特沃斯高通滤波器图一流程图本次设计的IIR数字滤波器系数的计算是根据已知的模拟滤波器的特性转换到等价的数字滤波器。
本次设计用双线性变换法得到数字滤波器。
而且，双线性变换法得到的数字滤波器保留了模拟滤波器的幅度响应特性。
3、设计方法频率变换法设计思想：1、从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S 域变换到Z 域，而得到所需类型的数字滤波器。
2、先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z 域内经数字频率变换为所需类型的数字滤波器。
4、IIR 巴特沃斯数字高通滤波器的实例计算（1）设计要求设计一个巴特沃斯数字高通滤波器,要求通带截止频率0.32p ωπ=rad ,通带衰减不大于1dB,阻带截止频率0.16p ωπ=rad，阻带衰减不小于15dB（2）确定数字高通的技术指标：通带截止频率0.32p ωπ=rad，通带最小衰减1p a =dB 阻带截止频率0.16p ωπ=rad，阻带最大衰减15s a =dB（3）将巴特沃斯高通数字滤波器的技术指标转换成巴特沃斯高通模拟滤波器的设计指标：令T=2s，预畸变校正得到的模拟边界频率：1tan 0.5494/,12ph p p rad s dB ωαΩ===；1tan0.2566/,152sh s s rad s dB ωαΩ===。
（4）模拟低通滤波器的技术指标计算如下：对通带边界频率归一化，由于本设计的为1dB 截止频率，所以1,1p c p dBλλα===把1p λ=和-sh Ω带入p p hss hλλ-Ω=Ω求得归一化巴特沃斯低通滤波器的阻带截止频率为：2.141p p hs shλλ-Ω==Ω，15s α=dB（5）设计归一化巴特沃斯模拟滤波器G（p）。
10.8757sp k ===：2.141s s ppλλλ==：l g 3.1349l g s p s pk N λ==：所以取N=4，根据巴特沃斯归一化低通滤波器参数表(见附录)可得归一化模拟低通原型系统函数G(p)为：4321() 2.613 3.4142 2.61311G p p p p p =++++221(0.76541)( 1.84781)p p p p =++++（6）利用频率变换公式p p p hsλΩ=将G(p)转换成模拟高通()HP H s ：p ()(p |p p ha sH s G λΩ==）4432234S +2.613S+3.4142S +2.6131S +S ph ph ph ph =ΩΩΩΩ把0.5494p hΩ=代入此式可得：42340.09110.4333 1.0305 1.4356ss s s s =++++用双线性变换法将模拟高通()a H s 转换成数字高通()H z ：令zz T S 11112--+-=,11123412341210.36471.45872.18811.45870.3647()()|12.05781.85450.78950.1331a z s zz z z z H z H s z z z z -----------=+-+-+==-+-+二、软件仿真工具及实现环境简介1、计算机辅助设计方法在优秀科技应用软件MATLAB 的信号处理工具箱中提供了一整套模拟，数字滤波器的设计命令和运算函数，方便准确，简单容行使得设计人员除了可按上述传统设计步骤快速的进行较复杂高阶选频滤波器的计算、分析外，还可通过原型变换直接进行各种典型数字滤波器设计，即应用MATLAB 设计工具从模拟原型直接变换成满足原定频域指标要求有数字滤波器。
2、MATLAB 直接设计IIR 巴特沃斯数字高通滤波器MATLAB 编程如下：fs=5000;wp=800*2/fs;ws=400*2/fs;rp=1;rs=15;Nn=128;[N,wn]=buttord(wp,ws,rp,rs)[b,a]=butter(N,wn,'high')freqz(b,a,Nn,fs)MATLAB 运算结果如下：N=4wn=0.2388b=0.3647-1.4587 2.1881-1.45870.3647 a= 1.0000-2.0578 1.8545-0.78950.1331三、滤波器结构对数字滤波器性能指标的影响分析1、IIR系统的基本网络结构IIR系统的基本网络结构有三种，即直接型、级联型和并联型。
（1）直接型N阶差分方程如下：∑∑==-+-=MiNiiiinyainx bny01) ()()(对应的系统函数为∑∑=-=--=Ni ii Mi i i z a z b z H 11)(设M=N=2，按照差分方程可以直接画出网络结构如下图（a）所示。
图中第一部分系统函数用)(1z H 表示，第二部分用)(2z H 表示，那么)()()(21z H z H z H =，当然也可以写成)()()(12z H z H z H =,按照该式，相当于将下图(a)中两部分流图交换位置，如下图（b）所示。
该图中节点变量w1=w2，因此前后两部分的延时支路可以合并，形成如下图(c)所示的网络结构流图，我们将下图(c)所示的这类流图称为IIR 直接型网络结构。
M=N=2时的系统函数为2211221101)(------++=z a z a z b z b b z H 对照下图(c)的各支路的增益系数与)(z H 分母分子多项式的系数可见，可以直接按照)(z H画出直接型结构流图。
IIR 网络直接型结构由bz 和az 写出数字滤波器系统函数：123412340.36471.45872.18811.45870.3647()12.05781.85450.78950.1331z z z z H z z z z z---------+-+=-+-+有H（z）写出差分方程如下：() 2.0578(1)1.8545(2)0.7895(3)0.1331(4)y n y n y n y n y n =---+---0.3647()1.4587(1)2.1881(2)1.4587(3)0.3647(4)x n x n x n x n x n +--+---+-直接型网络结构如下图：z x(n)y(n)0.3642.0578-1.85450.7895-0.1331-0.14582.18811.45870.3647z z z 直接型网络结构图在后面的分析中我们将发现，直接型系统对滤波器的性能控制作用不明显，极点对系数的变化不灵敏，易出现不稳定或较大误差，而且运算的累积误差较大。
因此，在设计时一般不选用直接型。
（2）级联型级联型结构是将系统传递函数H(z)写成具有实系数的二阶节的乘积。
将分子和分母多项式分解为各自的根，然后将一对复数共轭根（或者任意两个实数根）组合成二阶多项式。
在直接型表示的系统函数H(z)中，分子、分母均为多项式，且多项式的系数一般为实数。
现将分子、分母多项式分别进行因式分解，得到：∏∏=-=---=N r r Mr r z d z C Az H 1111)1()1()(上式中,A 是常数;Cr 和dr 分别表示H(z)的零点和极点。
由于多项式的系数是实数，Cr 和dr 是实数或者是共轭成对的复数，将共轭成对的零点（极点）放在一起，形成一个二阶多项式，其系数仍为实数；再将分子、分母均为实系数的二阶多项式放在一起，形成一个二阶网络)(z H j 。
)(z H j 如下式：2j 21j 12j 211j 0j j 1)(------++=z a z a z z z H βββ上式中)(z H j 表示一个一阶或二阶的数字网络的子系统函数，每个)(z H j 的网络结构均采用前面介绍的直接型网络结构，如下图所示，H(z)则由k 个子系统级联构成。
一阶和二阶直接型网络结构级联型的MATLAB 的表示与实现：将数字滤波器系统函数H（z）：123412340.36471.45872.18811.45870.3647()12.05781.85450.78950.1331z z z z H z z z z z---------+-+=-+-+用MATLAB 转换成级联型的程序如下：bz=[0.3647-1.4587 2.1881-1.45870.3647];az=[1.0000-2.0578 1.8545-0.78950.1331];[S,G]=tf2sos(bz,az)freqz(bz,az)说明：[S,G]=tf2sos(bz,az):实现直接型到级联型的变换。
B 和A 分别为直接型系统函数的分子和分母多项式系数向量。
返回L 级二阶级联型结构的系数矩阵S 和增益常数G。
MATLAB 运算结果如下：S = 1.0000-2.1807 1.1988 1.0000-0.89760.22721.0000-1.81900.8341 1.0000-1.16020.5859G =0.3647由S 和G 写出数字滤波器级联型系统函数：1212121212.1807 1.198811.8190.8341()0.3647()()10.89760.227211.16020.5859z z z z H z z z z z---------+-+=-+-+级联型网络结构如下图：-0.2961x(n)y(n)-0.8976-0.2272-2.18071.1988z 1z -z 1z --1.0003-0.123780.06189z 1z -z 1z -级联型网络结构图在级联型结构中，每一级分子的系数确定一对零点，分母的系数确定一对极点，因为子网络的零极点也即整体网络的零极点，所以整个系统的零极点都可以准确的由每一级的系数来调整和控制，这样便于调整滤波器的频率响应性能，其灵敏度特性优于直接型和正准型结构。