MATLAB课程设计

1、郑州轻工业学院课 程 设 计 任 务 书题目 基于双线性Z变换法的IIR数字带阻滤波器设计 专业、班级 电子信息工程 学号 姓名 主要内容、基本要求、主要参考资料等：主要内容：利用MATLAB软件读取一段含有噪声的 .wav格式的语音信号，然后基于FFT对该信号进行频谱分析；基于含噪语音信号的频谱确定滤波器的参数，利用双线性Z变换法设计IIR数字带阻滤波器，并利用所设计的滤波器对信号进行滤波处理。比较滤波前后语音信号的时域波形及频谱，分析滤波前后的语音变化。基本要求：1、基于含噪语音信号的频谱确定滤波器的参数；2、基于双线性Z变换法分别设计一个巴特沃斯IIR数字带阻滤波器和切比雪夫IIR数字带

2、阻滤波器；3、掌握利用wavread函数读取、播放 .wav格式语音信号的方法；4、对语音信号进行滤波，绘制滤波前后信号的时域波形及频谱；5、回放语音信号，分析滤波前后的语音变化。主要参考资料：1、从玉良.数字信号处理原理及其MATLAB实现M.北京：电子工业出版社.2009.72、胡广书.数字信号处理理论、算法与实现M.北京：清华大学出版社.2003,8完 成 期 限： 2010.6.222010.6.26 指导教师签名： 课程负责人签名： 2013年 6月 22日郑州轻工业学院课程设计说明书题目： 基于双线性Z变换法的IIR数字带阻滤波器设计 姓 名： 院 （系）： 电子信息工程 专业班级

3、： 电信工程 学 号： 指导教师： 蔡超峰 成 绩： 时间： 2015 年 6 月29日至 2015 年 7 月 3 日目 录1引言22课程设计要求42.1课程设计题目42.2课程设计目的42.3设计要求42.4数字滤波器的基本概念52.5数字滤波器设计的基本步骤52.5.1确定指标52.5.2模型逼近52.5.3实现性能分析和计算机仿真53 IIR数滤波器设计的原理与方法字63.1 IIR数字滤波器设计的原理63.2 IIR 数字滤波器设计的基本方法74 IIR带阻滤波器的MATLAB 设计84.1 IIR带阻滤波器的设计流程85 IIR带阻滤波器的仿真结果及波形126 IIR带阻滤波器的仿

4、真结果分析14总结16参考文献181.引言 在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号的处理和分析都是基于滤波器而进行的。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。但是，传统的数字滤波器的设计使用繁琐的公式计算，改变参数后需要重新计算，从而在设计滤波器尤其是高阶滤波器时工作量很大。利用MATLAB信号处理箱(Signal Proc

5、essing Toolbox)可以快速有效地实现数字滤波器的设计与仿真。本文设计一个IIR数字带阻滤波器。IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。在对滤波器实际设计时，整个过程的运算量是很大的。设计完成后对已设计的滤波器的频率响应要进行校核，要得到幅频、相频响应特性，运算量也是很大的。平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候要根据设计要求和滤波效果不断地调整，以达到设计的最优化。在这种情况下，滤波器设计就要进行大量复杂

6、的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成。利用MATLAB强大的计算功能进行计算机辅助设计，可以快速有效地设计数字滤波器，大大地简化了计算量。滤波器的功能是让一定频率范围内的信号通过，而将此频率范围之外的信号加以抑制或使其急剧衰减。当干扰信号与有用信号不在同一频率范围之内，可使用滤波器有效的抑制干扰。 用LC网络组成的无源滤波器在低频范围内有体积重量大，价格昂贵和衰减大等缺点，而用集成运放和RC网络组成的有源滤波器则比较适用于低频，此外，它还具有一定的增益，且因输入与输出之间有良好的隔离而便于级联。由于大多数反映生理信息的光电信号具有频率低、幅度小、易受干扰等特点，因而RC有

7、源滤波器普遍应用于光电弱信号检测电路中数字滤波器的实用型式很多，大略可分为有限冲激响应型和无限冲激响应型两类，可用硬件和软件两种方式实现。在硬件实现方式中，它由加法器、乘法器等单元所组成，这与电阻器、电感器和电容器所构成的模拟滤波器完全不同。数字信号处理系统很容易用数字集成电路制成，显示出体积小、稳定性高、可程控等优点。数字滤波器也可以用软件实现。软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数字滤波计算。带阻滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带阻滤波器的例子是电阻-电感-电容电路(RLC

8、circuit)。这些滤波器也可以用低通滤波器同高通滤波器组合来产生.2.课程设计要求2.1课程设计题目 数字滤波系统的MATLAB实现与分析双线性变换法设计数字带阻IIR滤波器2.2课程设计目的通过本次综合设计，运用已学的课程知识，根据题目要求进行软件仿真设计和调试，对数字信号处理课程中涉及的滤波器原理和设计方法有一定的感性认识和实践操作能力，从而加深对本课程知识点的理解，使学生应用知识能力、设计能力、调试能力以及报告撰写能力等方面有显著提高。学会MATLAB的使用，掌握MATLAB的程序设计方法；掌握在Windows环境下语音信号采集的方法；掌握数字信号处理的基本概念、基本理论和基本方法；

9、掌握MATLAB设计FIR和IIR数字滤波器的方法；学会用MATLAB对信号进行分析和处理。2.3设计要求本设计的主要内容包括掌握数字信号处理的基础知识和利用MATLAB软件进行编程，利用几种常见的数字滤波器对信号进行处理，并对过滤结构进行分析。具体要求包括以下几个方面：（1）滤波器与数字信号处理的理论知识学习和分析通过查找、学习相关文献资料，掌握滤波器与数字信号处理的基本理论知识，掌握滤波器设计的方法和基础知识，掌握滤波器过滤信号的基本过程和原理，并对相应信号进行时域和频域分析，理解信号的频谱变化的特点。（2）利用MATLAB对上述过程进行仿真利用MATLAB软件进行相应的设计，实现滤波器对

10、信号的过滤处理。（3）结果分析利用滤波器处理相应信号，对过滤结果进行分析，对比处理前后的信号频率的衰变，说明信号过滤的情况好坏，并分析该类型的滤波器的优缺点和适用场合。2.4数字滤波器的基本概念滤波器的种类很多，从功能上可分为低通、高通、带通和带阻滤波器，每一种又有模拟滤波器和数字滤波器两种形式。如果滤波器的输人和输出 都是离散时间信号，则该滤波器的冲击响应也必然是离散的，这种滤波器称之为数字滤波器。数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器也是具有一定传输选择特性的数字信号处理装置，其输入、输出均为数字信号，实质上是一个由有限精

11、度算法实现的线性时不变离散系统。它的基木工作原理是利用离散系统特性对系统输人信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输出。2.5数字滤波器设计的基本步骤数字滤波器的设计要经过如下三个步骤:2.5.1确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。2.5.2模型逼近一旦确定了技术指标，就可利用已学习过的基本原理和关系式，提出一个滤波器模型来逼近给定的指标体系。这是滤波器设计所要研究的主要问题。2.5.3实现性能分析

12、和计算机仿真以上两步的结果得到的滤波器，通常是以差分方程、系统函数或脉冲响应描述的。根据这个描述就可以分析其频率特性和相位特性、以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。3.IIR数字滤波器设计的原理与方法3.1 IIR数字滤波器设计的原理IIR 数字滤波器具有无限持续时间冲激响应,需要用递归模型来实现,其系统函数为: (2.1)设计IIR 滤波器的任务就是寻求一个因果、物理上可实现的系统函数H (z) ，使其频率响应满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止、通带衰减和阻带衰减.。利用冲激响应不变法设计数字滤波器时可能会导致频

13、域混叠现象，为了克服这一问题，需要找到由s平面到z平面的另外的映射关系，这种关系应保证：1) s平面的整个j轴仅映射为z平面单位圆上的一周；2) 若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的；3) 这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)；4)如果G(j0)=1，那么1。双线性变换满足以上4个条件的映射关系，其变换公式为 （2.2）双线性变换的基本思路：首先将整个平面压缩到一条从/Ts变换到2/Ts 的横带里，然后通过标准的变换关系将横带变换到整个平面上去，这样就得到了平面与平面间的一一对应的单值关系。3.2 IIR 数字滤波器设计的基本方法I

14、IR 数字滤波器的设计方法有两类,一类是借助于模拟滤波器的设计方法设计出模拟滤波器,利用冲激响应不变法或双线性变换法转换成数字滤波器,再用硬件或软件实现；另一类是直接在频域或时域中进行设计,设计时需要计算机作辅助工具。随着MATLAB 软件尤其是MATLAB 的信号处理工作箱的不断完善,不仅数字滤波器的计算机辅助设计有了可能,而且还可以使设计达到最优化。IIR 数字滤波器设计的基本步骤如下:(1) 根据任务,确定性能指标。在设计一个滤波器之前,首先根据工程实际的需要确定滤波器的技术指标如：边界频率：p ，s ，c ；阻带最小衰减As 和通带最大衰减Rp；(2) 将数字滤波器的技术指标转换成模拟

15、滤波器指标。利用冲激响应不变法与双线性变换法进行频率间的转换，主要是边界频率Wp与Ws 的转换。(3) 用模拟滤波器设计方法得到模拟滤波器的传输函数Ha ( s) ；可借助巴特沃斯(Butterworth) 滤波器、切比雪夫(Chebyshev) 滤波器、椭圆(Cauer) 滤波器、贝塞尔(Bessel) 滤波器等，这些滤波器都有严格的设计公式、现成曲线和图表供设计人员使用。4. IIR带通滤波器的MATLAB 设计4.1 IIR带阻滤波器的设计流程本文设计的IIR带阻滤波器是从低通变换过来的，利用的是双线性变换以及切比雪夫II滤波器的原型，其具体的设计流程为上图所示。首先根据题目要求确定带阻

16、滤波器的技术指标，先要进行频率的预畸变，并且归一化频率，再设计出切比雪夫II模拟低通滤波器，并求出其阶数等相关参数。其次利用双线性变换法设计数字带通滤波器，再调用函数进行双线性变换，并求出分子、分母的系数向量。最后通过画图求出其幅频响应、相频响应、幅度特性曲线与零极点，并画出波形图。最后进行验证，看所设计的滤波器能否达到要求的指标，若能达到，则说明该滤波器设计符合要求。3.2 IIR带阻滤波器的设计步骤（1）根据设计流程，首先确定所要设计的数字带阻滤波器的相关指标：通带截止频率wp1 =0.4，wp2 =0.6，通带最大衰减Rp = 2dB；阻带截止频率ws1=0.2 和ws2=0.8，阻带最

17、小衰减Rs = 30dB；取样间隔T = 0.1s。其实现程序如下（程序中pi代表）：clear;clc;close all;S,Fs,Bit = wavread(C:UsersyyyDesktopccSHE.wav); %读取音频信号S1,Fs,Bit = wavread(C:UsersyyyDesktopccBandstop.wav); %读取音频信号(含有噪声)sound(S,Fs);sound(S1,Fs);X = abs(fft(S); X1 = abs(fft(S1); xlabel = 0:Fs/(length(S)-1):Fs/2;subplot(2,1,1);plot(xla

18、bel,X(1:150001);title(原始信号频谱图); %分析音频信号，找出干扰频率段为8000Hzsubplot(2,1,2);plot(xlabel,X1(1:150001);title(加噪信号频谱图); % 对于带阻滤波器：fp = 5000;fs = 6000;Rp = 3;Rs = 40; % 求数字频率wp = 2*pi*fp/Fs;ws = 2*pi*fs/Fs;% 预畸变处理Wp = 2*Fs*tan(wp/2);Ws = 2*Fs*tan(ws/2);% 开始设计模拟低通滤波器N,Wn = buttord(Wp,Ws,Rp,Rs,s);B,A = butter(N,

19、Wn,s);b,a = bilinear(B,A,Fs); h,w = freqz(b,a,1000,whole);h = h(1:501);w = w(1:501);mag = abs(h);db = 10*log10(mag+eps)/max(mag);figure;subplot(1,2,1);plot(w/pi,mag);title(幅频响应);axis(0,1,0,1.1);subplot(1,2,2);plot(w/pi,db);title(幅频响应(dB);axis(0,1,-40,1);set(gca,XTickMode,manual,XTick,0,0.2,0.3,1);se

20、t(gca,YTickMode,manual,YTick,-40,-20,-3,0);grid; S2 = filter(b,a,S1);sound(S2,Fs);X2=abs(fft(S2);figure;plot(xlabel,X2(1:150001);title(滤波以后的频谱)5.IIR带阻滤波器频谱6. IIR带阻滤波器的仿真结果分析根据上述的仿真波形，可以看出：（1）由ChebyshevI 型、ChebyshevII型低通原型变换成带通模型的幅频特性可以用分贝形式表示。对于ChebyshevI 型的幅频响应在通带内为等纹波衰减、阻带内为单调减小的，且通带（归一化）截止频率在波形上显

21、示为：0.4与0.6，阻带（归一化）截止频率在波形上显示为：0.2与0.8，与设计要求基本一致。对于ChebyshevII型的幅频响应来说，在通带内单调减小、阻带内为等纹波衰减的，而通带、阻带（归一化）截止频率为：0.4与0.6，0.2与0.8，均满足设计要求。（2）此次设计的IIR带阻滤波器的阶数可由函数N,OmegaC=cheb1ord(Wp,Ws,Rp,Rs,s)或N,OmegaC=cheb2ord(Wp,Ws,Rp,Rs,s)求得，N=3。则从ChebyshevI 型带阻滤波器的幅频响应曲线上看，当N3时，在（归一化后为0-1）范围内通带波动只有两个波谷，阻带则单调减小，这与理论内容一

22、致。对于ChebyshevII 型带阻滤波器来说，其幅频响应曲线的通带波动在（归一化后为0-1）范围内阻带波动只有两个波峰，通带则单调减小，这与理论内容一致。因此，对于滤波器阶数这一指标来说，满足要求。（3）对于分贝化后的幅度特性，ChebyshevI 型带阻滤波器的通带最大衰减大约为2dB，阻带衰减最小为30dB，而ChebyshevII 型带阻滤波器的阻带最小衰减大约为30dB，通带最小衰减大约为2dB，这与带阻滤波器的技术指标相同。因此，在衰减性能上，所设计的滤波器达到了要求。（4）对于相频响应来说，ChebyshevI 型带阻滤波器在通带内的相频响应曲线接近为一条平滑曲线，在阻带内则存

23、在衰减、畸变。ChebyshevII 型带阻滤波器在通带截止频率以内的相频响应曲线几乎为一条直线，但在通带截止频率处有很大角度的转折，使相频特性在通带内产生了的畸变，而在阻带内存在很大的畸变。因此，对于双线性变换法，其相位特性得不到满足，必须用其它方法加以校正或采用其它方法重新设计滤波器，使之满足线性相位。（5）群延迟是衡量一个滤波器或是整个系统性能指标的一个重要参数。对于利用这两种低通滤波器原型模型设计的这两种带通滤波器，带阻滤波器的延迟是同等带宽的低通滤波器延迟的两倍。这个结果是由低通特性向带通特性变换引起的，n阶低通滤波器传递函数总是变换为2n阶带通滤波器传递函数，其滤波器的群延迟特性的中心频率与理论上的很接近。总结此次Matlab课程设计是在一定的理论基础之上进行的，在先修课程信号与系统与数字信号处理中，大量有Matlab设计方面的知识，所以做起来还比较容易，而且经过了自己的亲身实践，学到了许多实践方面的知识。首先，在信号滤波系统中，有时因为模拟滤波器阶数太高，硬件占用空间太大为某些仪器的实现设置了障碍，而对于