数字信号处理 IIR数字滤波器的设计 课程设计 吉林建筑大学

吉林建筑大学电气与计算机学院数字信号处理课程设计报告设计题目：IIR数字滤波器的设计专业班级：学生姓名：学号：指导教师：高晓红吕卅设计时间：2016.12.12－2016.12.16目录第1章概述 (2)1.1设计的作用、目的 (2)1.2设计任务及要求 (2)1.3设计内容 (2)第2章 IIR数字滤波器的设计 (3)2.1设计原理 (3)1. 数字低通滤波器的设计原理 (3)2. 变换方法的原理 (5)2.2设计步骤 (6)第3章数字低通滤波器MATLAB编程及幅频特性曲线 (7)3.1 MATLAB语言编程 (8)3.2 幅频特性曲线 (9)第4章总结 (11)参考文献 (13)第1章 概述1.1设计的作用、目的课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。

本次课程设计一方面通过MATLAB 仿真设计内容，使我们加深对理论知识的理解，同时增强其逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充。

1.2设计任务及要求用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，并观察所设计数字滤波器的幅频特性曲线。

通过课程设计各环节的实践，应使学生达到如下要求：（1）掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器以及窗函数法设计FIR 数字滤波器的原理、具体方法及计算机编程。

（2）观察双线性变换法、脉冲响应不变法及窗函数法设计的滤波器的频域特性，了解各种方法的特点。

（3） 用MATLAB 画出三种方法设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查结果是否满足要求。

1.3设计内容设计题目：IIR 数字滤波器的设计设计内容：已知通带截止频率kHz f p 2.0=,通带最大衰减dB P 1=α, 阻带截止频率kHz f s 3.0=，阻带最小衰减dB s 25=α，T=1ms ，按照以上技术要求，用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，并观察所设计数字滤波器的幅频特性曲线。

第2章 IIR数字滤波器的设计2.1设计原理1.数字低通滤波器的设计原理数字滤波器可以分为两大类：一类是经典滤波器，即一般的滤波器特点是输入信号中有用的频率成分和希望滤去的频率成分各占不同的频率带，通过一个合适的选频滤波器达到绿宝锝目的，这种滤波器主要有维纳滤波器尔曼滤波器，自适应滤波器等。

数字滤波器是一个离散时间系统。

应用数字滤波器处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。

数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍。

数字滤波器的频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率（即二分之一抽样频率点）呈镜像对称。

为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。

数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。

数字滤波器在语声信号处理、图像信号处理、医学生物信号处理以及其他应用领域（如通信、雷达、声纳、仪器仪表和地震勘探等）都得到了广泛的应用。

数字滤波器有低通、高通、带通、带阻和全通等类型。

它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。

如果数字滤波器的内部参数不随时间而变化，则称为时不变的，否则为时变的。

如果数字滤波器在某一给定时刻的响应与在此时刻以后的激励无关，则称为因果的，否则为非因果的。

如果数字滤波器对单一或多个激励信号的响应满足线性条件，则称为线性的，否则为非线性的。

应用最广的是线性、时不变数字滤波器。

数字滤波器可以按所处理信号的维数分为一维、二维或多维数字滤波器。

一维数字滤波器处理的信号为单变量函数序列，例如时间函数的抽样值。

二维或多维数字滤波器处理的信号为两个或多个变量函数序列。

例如，二维图像离散信号是平面坐标上的抽样值。

与模拟滤波器相同，数字滤波器从功能上看可分为低通、高通、带通和带阻几类。

因为它们的单位脉冲响应是非因果且无限长的，所以实际上理想滤波器是不可能实现的。

与模拟滤波器不同的是数字滤波器的传输函数都是以2π为周期的，滤波器的低通频带处于2π的整数倍处，而高频频带处于π的奇数倍附近。

理想滤波器就是一个让输入信号中的某些有用频谱分量无任何变化的通过，同时又能完全抑制另外那些不需要的成分的具有某种选择性的器件、网络或以计算机硬件支持的计算机程序。

根据对不同信号的处理可将滤波器分为模拟滤波器和数字滤波器两种。

模拟滤波器和数字滤波器的概念相同，只是信号的形式和实现滤波方法不同。

数字滤波器是指输入输出都是数字信号的滤波器。

滤波器的滤波原理局势根据信号与噪声占据不同的频带，将噪声的频率放在滤波器的阻带中，而由于阻带的响应为零，这样就滤去了噪声。

一个理想滤波器将。

是一个无法实现的非因果系统，我们只能用一个稳定的因果系统函数去逼近工程需要的确定的性能要求。

数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。

IIR滤波器的特征是，具有无限持续时间冲激响应。

这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。

FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。

数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。

随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。

（1）确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。

在很多实际应用中，数字滤波器常常被用来实现选频操作。

因此，指标的形式一般在频域中给出幅度和相位响应。

幅度指标主要以两种方式给出。

第一种是绝对指标。

它提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。

第二种指标是相对指标。

它以分贝值的形式给出要求。

在工程实际中，这种指标最受欢迎。

对于相位响应指标形式，通常希望系统在通频带中具有线性相位。

运用线性相位响应指标进行滤波器设计具有如下优点：①只包含实数算法，不涉及复数运算；②不存在延迟失真，只有固定数量的延迟；③长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。

（2）逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。

通常采用理想的数字滤波器模型。

之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。

（3）性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。

根据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。

（4）滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。

一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。

以低通滤波器特性为例，频率响应有通带、过渡带及阻带三个范围。

在通带内：1|)(|Ap -1≤≤ωj e H c ||ωω≤在阻带内：st j A e ≤|)(H |ω c st ωωω≤≤||其中Wc 为通带截止频率，Ap 为通带误差，Ast 为阻带误差。

与模拟滤波器类似，数字滤波器按频率特性划分为低通、高通、带通、带阻、全通等类型，由于数字滤波器的频率响应是周期性的，周期为2Pi 。

2.变换方法的原理利用模拟滤波器成熟的理论及其设计方法来设计IIR 数字低通滤波器是常用的方法。

2.2.1设计过程按照数字滤波器技术指标要求设计一个过度模拟低通滤波器，在按照一定的转换关系将转换成数字低通滤波器的系统函数。

由此可见，设计的关键问题就是找到这种转换关系，将s 平面上的转换成z 平面上的。

为了保证转换后的稳定且满足技术指标要求，对转换关系提出两点要求：（1）因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。

（2）数字滤波器的频率响应模仿模拟滤波器的频响特性，s 平面的虚轴映射为z 平面的单位圆，响应的频率之间呈线性关系。

一个无限脉冲响应系统的传递函数具有可选择极点和可选择零点（而有线脉冲响应系统只具有可选择零点以及位于原点的固定极点）。

模拟系统同样具有极点和零点，甚至经常只有极点（多项式滤波器或全几点滤波器）。

由此得出结论是，通过把PN 图从s 域映像变换到z 域中，可以用一个无限脉冲（IIR ）系统对模拟滤波器进行数字仿真。

目前，IIR 数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。

模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便。

将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再由冲击响应不变法或双线性变换将H(s)转换为所需的H(z)。

s - z 映射的方法有：冲激响应不变法、阶跃响应不变法、双线性变换法等。

实际中有三种广泛应用的滤波器原型，既巴特沃斯低通滤波器、切比雪夫低通滤波器、椭圆低通滤波器。

2.2.2巴特沃斯低通滤波器这种滤波器的特征是其通带和阻带都有平坦的幅度响应。

N 阶低通滤波器的平方幅度响应的表达式为： N ca j j j H A 222)(11|)(|)(ΩΩ+=Ω=Ω 通带： 使信号通过的频带阻带：抑制噪声通过的频带过渡带：通带到阻带间过渡的频率范围Ωc ：截止频率2.2.3脉冲响应不变法设计IIR 数字滤波器利用模拟滤波器成熟的理论及其设计方法来设计IRR 数字低通滤波器是通常的方法。

从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应模仿模拟滤波器的冲击响应，即h(n)是ha(t)的采样值。

设T 为采样周期，变换过程为：)()()()(z H n h t h s H a a −→−−→−−→−利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可以从不同的角度出发。

脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t)，即将ha(t)进行等间隔采样，使h(n)正好等于ha(t)的采样值，满足h(n)=ha(nT)式中,T 是采样周期。

如果令Ha(s)是h 其中T 为采样间隔。

如果以Ha(s)及H(z)分别表示ha(t)的拉氏变换及h(n)的Z 变换，则a(t)的拉普拉斯变换，H(z)为h(n)的Z 变换，利用采样序列的Z 变换与模拟信号的拉普拉斯变换的关系。

2.2设计步骤求模拟滤波器的系统函数Ha （S ）(1)确定阶数N ： 1101101.01.0--=p s sp K αα fpfs∏∏=22sp λsp Ksp N λlg lg =（2）查表确定极点，写出Ga （P ）（3）将Ga （P ）进行归一化，求出Ha （S ）（4）公式：①脉冲响应不法 ∑=-=Ni iis s A z H 1)()()(1t u e A t h ts Ni i a i ∑==)()(1nT u e A n h nT sNi i i ∑==∑=--=Ni t siz e A i 111H(z)②双线性变换法)T 21tan(T 21Ω=ΩT T1111T 2j Ω-Ω-+-=Ωj j e e1111T 2s --+-=z zS-+=T 2ST 2z第3章数字低通滤波器MATLAB编程及幅频特性曲线3.1 MATLAB语言编程1.脉冲响应不变法Fs=1000;Wp=400*pi;Ws=600*pi;As=25;Ap=1;[N,Wc]=buttord(Wp,Ws,Ap,As,'s');[B,A]=butter(N,Wc,'s');W=(0:pi:1000);[D,C]=impinvar(B,A,Fs);W=(0:0.001:pi)Hz=freqz(D,C,W);plot(W,abs(Hz)/abs(Hz(1)));grid;title('脉冲响应不变法');xlabel('频率');ylabel('幅值\dB');2.双线性变换法wp=400*pi;ws=600*pi;Ap=1;As=25;T=0.001;fs=1/T;[N,Wn]=buttord(wp,ws,Ap,As,'s');[b,a]=butter(N,Wn,'s');[numd,dend]=bilinear(b,a,fs);w=linspace(0,pi,500);h=freqz(numd,dend,w);plot(w/pi,20*log10(abs(h))); axis([0 1 -50 0]);grid;title('双线性变换法');xlabel('归一化频率');ylabel('幅度/dB');3.2 幅频特性曲线第4章总结通过这次数字信号处理的设计，让我感受到“温故而知新”是的重要性，我们总是以为知识学过了，就不用再花费时间去复习，其实，这是错误的。