IIR数字滤波器的设计

1、数字信号处理课程设计报告题 目：IIR 数字滤波器的设计学 院：化工过程自动化学院专业班级：学 号：姓 名：指导教师:起止日期：2015年6月22日2015年6月28日1 课程设计的意义与任务要求 11.1 课程设计的意义 11.2 课程设计的任务要求 12 课程设计的理论基础 12.1 数字滤波器简介 22.2 IIR 数字滤波器的设计原理 22.3 IR 数字滤波器的特点 33 MATLA瞅件介绍 33.1 MATLA颤件介绍 33.2 MATLA廊用领域 43.3 MATLAB!关语句 44 课程设计的具体内容 54.1 数字滤波器设计步骤 54.2 脉冲响应不变法和双线性变换法的变换原

2、理和步骤. 错 误 ! 未定义书签。4.2.1 脉冲响应不变法的变换原理和步骤 错 误 ! 未定义书签。4.2.2 双线性变换法的变换原理和步骤 64.3 实验步骤及运行程序 65课程设计的总结与心得 10参考文献 11数字信号处理课程设计1、课程设计的意义与任务要求1.1 课程设计的意义数字滤波器是具有一定传输选择特性的数字信号处理装置，其输入、输出均为 数字信号，实质上是一个由有限精度算法实现的线性时不变离散系统。它的基 本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序 列的频谱或信号波形，让有用频率的信号分量输出。数字滤波器和模拟滤波器 有着相同的滤波概念，根据其频率

3、响应特性可分为低通、高通、带通、带阻等 类型，与模拟滤波器相比，数字滤波器除了具有数字信号处理的固有优点外， 还有滤波精度高、稳定性好、灵活性强等优点。1.2 课程设计的任务要求（1）熟悉用脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理与方法；（2）学会调用MATLABS号处理工具箱中滤波器设计函数设计 IIR数字滤 波器，学会根据滤波需求确定滤波器指标参数。2、课程设计的理论基础利用MATLABJ号处理工具箱中的滤波器设计和分析工具 （FDATool）可以很方便 地设计出符合应用要求的未经量化的IIR数字滤波器。需要将MATLABS计出的 IIR数字滤波器进一步分解和量化，从而获得可

4、用 FPG故现的滤波器系数。IIR数字滤波器的设计方法有两类：间接设计法和直接设计法。间接设计法是借 助模拟滤波器设计方法进行设计的，先根据数字滤波器设计指标设计相应的过 渡模拟滤波器，再将过渡模拟滤波器转换为数字滤波器。直接设计法师在时域 或频域直接设计数字滤波器。由于模拟滤波器设计理论非常成熟，而且有很多性能优良的典型滤波器可供选 择（如，巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器等 ），设计公式和图表 完善，而且许多实际应用需要模拟滤波器的数字仿真，所以间接设计法得到广 泛的应用。而直接设计法要求解联立方程组，必须采用计算机辅助设计。在计 算机普及的今天，各种设计方法都有现成的设计程序（

5、或设计函数）可供调用， 例如利用MATLA防真平台，可以设计不同类型的IIR滤波器。2.1 数字滤波器简介数字滤波器是一种用来过滤时间离散信号的数学系统，通过对抽样数据进行数 学处理来达到频域滤波的目的。可以设计系统的频率响应，让它满足一定的要 求，从而对通过该系统的信号的某些特定的频率成分进行过滤，这就是滤波器 的基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统 是一个离散系统，则滤波器称为数字滤波器。数字滤波器是一个离散时间系统(按预定的算法，将输入离散时间信号(对应数 字频率)转换为所要求的输出离散时间信号的特定功能装置 )。应用数字滤波器 处理模拟信号(对应模拟频率)

6、时，首先须对输入模拟信号进行限带、抽样和模 数转换。数字滤波器输入信号的数字频率 (2九*f/fs,f为模拟信号的频率，fs为采样频率，注意区别于模拟频率)，按照奈奎斯特抽样定理，要使抽样信号的 频谱不产生重叠，应小于折叠频率(ws/2=兀)，其频率响应具有以2冗为间隔的 周期重复特性，且以折叠频率即二九点对称。为得到模拟信号，数字滤波器处 理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、 可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像 信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。2.2 IIR数字滤波器的设计原理IIR数字滤波器的最通

7、用的方法是借助于模拟滤波器的设计方法。模拟滤波器 设计已经有了相当成熟的技术和方法，有完整的设计公式，还有比较完整的图 表可以查询，因此设计数字滤波器可以充分利用这些丰富的资源来进行。对于IIR数字滤波器的设计具体步骤如下：(1)按照一定的规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的 技术指标。(2)根据转换后的技术指标设计模拟低通滤波器G(s)(G(s)是低通滤波器的传递函数)。再按照一定的规则将G(s)转换成H(z)(H(z)是数字滤波器的传递函数)。若 设计的数字滤波器是低通的，上述的过程可以结束，若设计的是高通、带通或 者是带阻滤波器，那么还需要下面的步骤：将高通、带通或带阻

8、数字滤波器的技术指标转换为低通模拟滤波器的技术指标， 然后设计出低通G(s),再将G(s)转换为H(z)。2.3 IIR数字滤波器的特点(1) IIR 数字滤波器的系统函数可以写成封闭函数的形式。(2) IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算 结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准 型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理， 使误差不断累积，有时会产生微弱的寄生振荡。(3) IIR 数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、 契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工

9、作量比较 小，对计算工具的要求不高。在设计一个 IIR数字滤波器时，我们根据指标先 写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数 字滤波器的公式。(4) IIR 数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准 网络。3、关于 MATLAB3.1 MATLAB软件介绍MATLAB矩阵实验室)是MATrix LABoratory 的缩写，是一款由美国TheMathWorks公司出品的商业数学软件。MATLAB1一种用于算法开发、数据可视 化、数据分析以及数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB®

10、;可以用来创建用户界面及与调用 其它语言(包括C, C+林口 FORTRAN编写的程序。1MATLABF仅在数值运算上继续保持着相对其他同类软件的绝对优势,而且还开 发了自己的符号运算功能。用户只要学会了 MATLA就可以方便解决诸如矩阵变 换及运算、多项式运算、微积分运算、线性与非线性方程求解、常微分方程求 解、偏微分方程求解、插值与拟合、统计及优化等问题。而且MATLAB发展到今天，其可扩充性和可开发性起着不可估量的作用。MATLA冰身就像一个解释系统，对其中的函数程序的执行以一种解释执行的方 式进行，具最大好处是 MATLA瑞全成了一个开放的系统，用户可以方便的与FORTRANC等语言进

11、行连接，以充分利用各种资源。用户只需将已有的EXE文件转换成ME双件，就可以方便地调用有关程序和子程序。3.2 MATLAB应用领域MATLA且一种应用于科学计算领域的高级语言，它的主要功能包括数值计算和 符号计算功能、绘图功能、编程语言及应用工具箱。尽管MATLASE要用于数值 运算，但利用为数众多的附加工具箱（Toolbox）它也适合不同领域的应用，例 如控制系统设计与分析、图像处理、信号处理与通讯、金融建模和分析等。另 外还有一个配套软件包 Simulink ,提供了一个可视化开发环境，常用于系统模 拟、动态/嵌入式系统开发等方面。MATLABM以进行矩阵运算、绘制函数和数 据、实现算法

12、、创建用户界面、连接其他编程语言的程序等，主要应用于工程 计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分 析等领域。2MATLAB1一个包含大量计算算法的集合。 其拥有600多个工程中要用到的数学 运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都 是科研和工程计算中的最新研究成果，而且经过了各种优化和容错处理。在通 常情况下，可以用它来代替底层编程语言，如C和C+。在计算要求相同的情况下，使用MATLAB勺编程工作量会大大减少。MATLAB勺这些函数集包括从最 简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数 所能解决的问题其大致包

13、括矩阵运算和线性方程组的求解、微分方程及偏微分 方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问 题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数 组操作以及建模动态仿真等。33.3 MATLAB相关语句convert转换量化滤波器的结构copyobj*为量化器、量化滤波器或量化 FFT创建一个独立的复制品disp*显示量化器、量化滤波器或量化 FFTeps*返回量化器、量化滤波器或量化 FFT的量化阶数filter*将一量化滤波器应用于数据，可访问状态和滤波信息 freqz*计算量化滤波器的频响特性get*返回量化器、量化滤波器或量化FFT的属性imp

14、z*计算量化滤波器的单位冲激响应isallpass测试量化滤波器是否为全通特性isfir测试是否为FIR滤波器islinphase测试是否为线性相位noperations*量化器、量化滤波器或量化 FFT的量化操作次数noverflows*返回最后一次FFT或IFFT运算的溢出次数qfft 量化 FFTqreport*显示应用量化器、量化滤波器或量化FFT的结果quantizer*构建一个量化器range*返回量化器的数值范围reset*将一个或多个量化器、量化滤波器或量化FFT复位setbits*将量化器、量化滤波器或量化 FFT的数据设置为Format形式 4、设计的具体内容4.1 数字滤

15、波器设计步骤用MATLA进行模拟原型的数字滤波器的设计，一般步骤如下：(1)按一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术 指标。(2)根据转换后的技术指标使用滤波器阶数选择函数， 确定最小阶数N和固有 频率 Wn根据选用的模拟低通滤波器的类型可分为：buttord , cheblord ,cheb2ord , ellipord 等函数；(3)运用最小阶数 N产生模拟滤波器原型，模拟低通滤波器的创建函数有： buttap , cheb1ap, cheb2ap, ellipap , besselap 等；(4)运用固有频率 Wn把模拟滤波器原型转换成模拟低通、高通、带通、带阻滤

16、波器,可分别用函数lp21P , lp2hp , lp2bp , lp2bs ;(5)运用脉冲响应不变法和双线性不变法把模拟滤波器转换成数字滤波器，分别用函数impinva和bilinear 来实现。低通Chebyshevl型数字滤波器的设计： 设计中需要限定其通带上限临界频率 Wp阻带临界滤波频率 Ws在通带内的最 大衰减Rp,阻带内的最小衰减Rso 44.2 脉冲响应不变法和双线性变换法的变换原理与步骤4.2.1 脉冲响应不变法的变换原理与步骤 设计巴特沃斯低通滤波器，要求通带边界频率fp=2.1KHZ,通带最大衰减Rp=0.5dB;阻带边界频率fs=8KHZ,阻带最小衰减Rs=30 dB

17、,采样频率为Fs= 20 KHZ。从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应模仿模拟滤波器的冲 击响应，即h(n)是ha(t)的采样值。tnT h(n) ZT H(z)设T为采样周期，变换过程：Ha(s) ilt ha(t)如果模拟滤波器的系统函数只有单阶极点，且分母的阶数高于分子阶数，用脉冲响应不变法求数字滤波器的系统函数有简便方法：将Ha(S)展成部分分式的NHa(s)并联形式，再利用下述变换公式直接写出H(z)AkH(z) N TAks Skk 1 1 eSkTz 14.2.2 双线性变换法的变换原理和步骤5(1)保证S平面压缩到S1平面的宽为2nt/T的横带内(2)保证低频部分

18、基本对应Ctan(-T) Ctan(-)彳 tan(2) T tan(-y)根据要求，确定数字滤波器指标。如是模拟频率临界点，则要先转变成数字频率，以便预畸变处理。将数字指标转换成与 Ha(s)对应的模拟性能指标。设计模拟滤波器的系统函数 Ha(s)。将映射关系代入Ha(s)中得数字滤波器系1统函数 H(z) ° H(z) Ha(吼 2L Ha,1TT 1 z1T 1 z由于数字滤波器传输函数只与频域的相对值有关，故在设计时可先将滤波器设计指标进行归一化处理。设采样频率为Fs,归一化频率的计算公式是：归_化频率实际模拟角频率(rad/s) 实际数字频率实际模频率(Hz)“FsFs/2

19、4.3 实验步骤及运行程序一、利用典型法设计数字滤波器的步骤：1、将设计指标归一化处理。如果采用双线性变换法，还需进行预防变。2、根据归一化频率，确定最小阶数N和频率参数 Wn可供选用的阶数选择函数有：buttord , cheblord , cheb2ord, ellipord 等。3、运用最小阶数N设计模拟低通滤波器原型。模拟低通滤波器的创建函 数有:buttap, cheblap, cheb2ap, ellipap 和 besselap ,这些函数输出的是零 极点增益形式，还要用zp2tf函数转换成分子分母多项式形式。如果想根据最 小阶数直接设计模拟低通滤波器原型，可用 butter,

20、chebyl, cheby2, ellip, bessel等函数，只是注意要将函数中的 Wn设为1。4、根据第2步的频率参数 Wn将模拟低通滤波器原型转换成模拟低通、 高通、带通、带阻滤波器，可用函数分别是：lp21p, lp2hp, lp2bp, lp2bs 。5、运用脉冲响应不变法或双线性变换法把模拟滤波器转换成数字滤波器， 调用的函数是impinvar和bilinear 。脉冲响应不变法适用于采样频率大于4倍截止频率的锐截止低通带通滤波器，而双线性变换法适合于相位特性要求不 高的各型滤波器。6、根据输出的分子分母系数，用tf函数生成H (z)的表达式，再用freqz 函数验证设计结果。用

21、冲激响应不变法设计的巴特沃斯数字低通滤波器的M程序如下：fp=2100;fs=8000;Fs=20000;Rp=0.5;Rs=30;T=1/Fs;%设计指标W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率N,Wn=buttord(W1p,W1s,Rp,Rs, 's');%!定butterworth的最小介数解口频率参数Wn z,p,k=buttap(N); %设计模拟低通原型的零极点增益参数bp,ap=zp2tf(z,p,k); %将零极点增益转换成分子分母参数bs,as=lp21P(bp,ap,Wn*pi*Fs);%将低通原型转换为模拟低通11bz,az=imp

22、invar(bs,as,Fs);% sys=tf(bz,az,T); % H,W=freqz(bz,az,512,Fs); % subplot(2,1,1);plot(W,20*log10(abs(H); %grid on;%xlabel('频率/Hz');ylabel('振幅/dB');subplot(2,1,2);plot(W,abs(H);grid on;xlabel('频率/Hz');ylabel('振幅/dB');运行后的波形如下用脉冲响应不变法进行模数转换给出传递函数H(z)生成频率响应参数绘制幅频响应加坐标网格010

23、00 2000 3000 4000 5000 6M0 7000 8000 9000 10000频率，出-401505 一01000 2000 3000 4000 5000 6000 7000 SOOO 9000 10000频率上E运行结果:N=4bz= 0.0000 0.0999 0.1914 0.0252az= 1.0000 -1.4336 1.0984 -0.4115 0.0627fp=2.1KHZ,Rs=30 dB,采样可以看出：通过编程，结果十分直观，题目中要求的通带边界频 通带最大衰减Rp=0.5dB;阻带边界频率fs=8KHZ,阻带最小衰减 频率为Fs= 20 KHZ,等设计指标都

24、达到，滤波器设计符合要求。 二、用双线性变换法设计椭圆数字低通滤波器的 M程序如下： fs=20000;wp=2*pi*2100/fs;ws=2*pi*8000/fs;Rp=0.5;Rs=30;Ts=1/fs;Wp=2/Ts*tan(wp/2);Ws=2/Ts*tan(ws/2); N,Wn=ellipord(Wp,Ws,Rp,Rs, 's'); z,p,k=ellipap(N,Rp,Rs);Bap,Aap=zp2tf(z,p,k);b,a=lp21P(Bap,Aap,Wn);bz,az=bilinear(b,a,fs);H,f=freqz(bz,az,512,fs);subp

25、lot(2,1,1);plot(f,20*log10(abs(H);grid on;xlabel('频率 /Hz');ylabel('振幅/dB');subplot(2,1,2);plot(f,abs(H);grid on;xlabel('频率/Hz');ylabel('振幅/dB');oJI-Jn i mfvJSwo o 20 no OJLI 3O no 6z o/H ram领 o on 4DRno9C00010U0 2000 3口 W 4000 5000 印加 7000 2000 9000 10000频率/H工运行结果N=2b