《MATLAB课程设计要点.pdf》由会员分享,可在线阅读,更多相关《MATLAB课程设计要点.pdf(21页珍藏版)》请在三一文库上搜索。
1、郑州轻工业学院 课 程 设 计 任 务 书 题目基于双线性 Z 变换法的 IIR 数字带阻滤波器设计 专业、班级电子信息工程学号姓名 主要内容、基本要求、主要参考资料等: 主要内容 : 利用 MATLAB 软件读取一段含有噪声的.wav 格式的语音信号,然后基于FFT 对 该信号进行频谱分析; 基于含噪语音信号的频谱确定滤波器的参数,利用双线性 Z 变换 法设计 IIR 数字带阻滤波器,并利用所设计的滤波器对信号进行滤波处理。比较滤波前 后语音信号的时域波形及频谱,分析滤波前后的语音变化。 基本要求 : 1、基于含噪语音信号的频谱确定滤波器的参数; 2、 基于双线性 Z 变换法分别设计一个巴特
2、沃斯IIR 数字带阻滤波器和切比雪夫IIR 数字带阻滤波器; 3、掌握利用 wavread函数读取、播放.wav 格式语音信号的方法; 4、对语音信号进行滤波,绘制滤波前后信号的时域波形及频谱; 5、回放语音信号,分析滤波前后的语音变化。 主要参考资料: 1、 从玉良 .数字信号处理原理及其MATLAB 实现M. 北京:电子工业出版社 .2009.7 2、胡广书 .数字信号处理理论、算法与实现M. 北京:清华大学出版社 .2003,8 完成期限:2010.6.222010.6.26 指导 教 师 签名 : 课程负责人签名: 2013年 6 月 22 日 郑州轻工业学院 课程设计说明书 题目:基
3、于双线性 Z 变换法的 IIR 数字带阻滤波器设计 姓名: 院 (系) :电子信息工程 专业班级:电信工程 学号: 指导教师:蔡超峰 成绩: 时间:2015 年 6 月 29 日至 2015 年 7 月 3 日 目 录 1 引言 错误!未定义书签。 2 课程设计要求. 3 2.1 课程设计题目 . 3 2.2 课程设计目的 . 3 2.3 设计要求 . 4 2.4 数字滤波器的基本概念 . 5 2.5 数字滤波器设计的基本步骤 . 5 2.5.1 确定指标 5 2.5.2 模型逼近 6 2.5.3 实现性能分析和计算机仿真 6 3 IIR数滤波器设计的原理与方法字. . 6 3.1 IIR数字
4、滤波器设计的原理 . 6 3.2 IIR 数字滤波器设计的基本方法. 8 4 IIR带阻滤波器的 MATLAB 设计 9 4.1 IIR带阻滤波器的设计流程 . 9 5 IIR带阻滤波器的仿真结果及波形. 10 6 IIR带阻滤波器的仿真结果分析. 14 总结 . 16 参考文献 . 17 1 1. 引言 在现代通信系统中, 由于信号中经常混有各种复杂成分,所以很多信号的处理和分 析都是基于滤波器而进行的。 数字滤波器可以理解为是一个计算程序或算法,将代表输 入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中, 使信号 按预定的形式变化。数字滤波器有多种分类,根据数字滤波器冲
5、激响应的时域特征,可 将数字滤波器分为两种,即无限长冲激响应(IIR) 滤波器和有限长冲激响应(FIR) 滤波 器。但是,传统的数字滤波器的设计使用繁琐的公式计算,改变参数后需要重新计算, 从而在设计滤波器尤其是高阶滤波器时工作量很大。利用MATLAB 信号处理箱 (Signal Processing Toolbox)可以快速有效地实现数字滤波器的设计与仿真。 本文设计一个 IIR 数字带阻滤波器。 IIR 数字滤波器具有无限宽的冲激响应,与模 拟滤波器相匹配,所以 IIR 滤波器的设计可以采取在模拟滤波器设计的基础上进一步变 换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设
6、计法。在对 滤波器实际设计时,整个过程的运算量是很大的。 设计完成后对已设计的滤波器的频率响应要进行校核,要得到幅频、相频响应特性, 运算量也是很大的。平时所要设计的数字滤波器,阶数和类型并不一定是完全给定的, 2 很多时候要根据设计要求和滤波效果不断地调整,以达到设计的最优化。 在这种情况下, 滤波器设计就要进行大量复杂的运算,单纯的靠公式计算和编制简单的程序很难在短时 间内完成。利用MATLAB 强大的计算功能进行计算机辅助设计,可以快速有效地设计数 字滤波器,大大地简化了计算量。 滤波器的功能是让一定频率范围内的信号通过,而将此频率范围之外的信号加以抑 制或使其急剧衰减。 当干扰信号与有
7、用信号不在同一频率范围之内,可使用滤波器有效 的抑制干扰。用 LC 网络组成的无源滤波器在低频范围内有体积重量大,价格昂贵和 衰减大等缺点,而用集成运放和RC 网络组成的有源滤波器则比较适用于低频,此外, 它还具有一定的增益, 且因输入与输出之间有良好的隔离而便于级联。由于大多数反映 生理信息的光电信号具有频率低、幅度小、易受干扰等特点,因而RC 有源滤波器普遍 应用于光电弱信号检测电路中 数字滤波器的实用型式很多,大略可分为有限冲激响应型和无限冲激响应型两类, 可用硬件和软件两种方式实现。 在硬件实现方式中, 它由加法器、乘法器等单元所组成, 这与电阻器、 电感器和电容器所构成的模拟滤波器完
8、全不同。数字信号处理系统很容易 用数字集成电路制成,显示出体积小、稳定性高、可程控等优点。数字滤波器也可以用 软件实现。软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数 3 字滤波计算。 带阻滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰 减到极低水平的滤波器, 与带阻滤波器的概念相对。 一个模拟带阻滤波器的例子是电阻 -电感 -电容电路 (RLC circuit) 。这些滤波器也可以用低通滤波器同高通滤波器组合来产 生. 2. 课程设计要求 2.1 课程设计题目 数字滤波系统的MATLAB 实现与分析 双线性变换法设计数字带阻IIR 滤波器 2.2 课程
9、设计目的 通过本次综合设计, 运用已学的课程知识, 根据题目要求进行软件仿真设计和调试, 对数字信号处理 课程中涉及的滤波器原理和设计方法有一定的感性认识和实践操作 能力,从而加深对本课程知识点的理解,使学生应用知识能力、设计能力、调试能力以 及报告撰写能力等方面有显著提高。学会MATLAB 的使用,掌握 MATLAB 的程序设 计方法;掌握在 Windows 环境下语音信号采集的方法;掌握数字信号处理的基本概念、 4 基本理论和基本方法;掌握MATLAB 设计 FIR 和 IIR 数字滤波器的方法;学会用 MATLAB 对信号进行分析和处理。 2.3 设计要求 本设计的主要内容包括掌握数字信
10、号处理的基础知识和利用MATLAB 软件进行编 程,利用几种常见的数字滤波器对信号进行处理,并对过滤结构进行分析。具体要求包 括以下几个方面: (1)滤波器与数字信号处理的理论知识学习和分析 通过查找、学习相关文献资料, 掌握滤波器与数字信号处理的基本理论知识,掌握 滤波器设计的方法和基础知识,掌握滤波器过滤信号的基本过程和原理,并对相应信号 进行时域和频域分析,理解信号的频谱变化的特点。 (2)利用 MATLAB 对上述过程进行仿真 利用 MATLAB 软件进行相应的设计,实现滤波器对信号的过滤处理。 (3)结果分析 利用滤波器处理相应信号, 对过滤结果进行分析, 对比处理前后的信号频率的衰
11、变, 说明信号过滤的情况好坏,并分析该类型的滤波器的优缺点和适用场合。 5 2.4 数字滤波器的基本概念 滤波器的种类很多,从功能上可分为低通、高通、带通和带阻滤波器,每一种又有 模拟滤波器和数字滤波器两种形式。 如果滤波器的输人和输出都是离散时间信号,则该滤波器的冲击响应也必然是离 散的,这种滤波器称之为数字滤波器。数字滤波器是一种用来过滤时间离散信号的数字 系统,通过对抽样数据进行数学处理来达到频域滤波的目的。 数字滤波器也是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数 字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基木工作原理 是利用离散系统特性对系统输人
12、信号进行加工和变换,改变输入序列的频谱或信号波 形,让有用频率的信号分量通过,抑制无用的信号分量输出。 2.5 数字滤波器设计的基本步骤 数字滤波器的设计要经过如下三个步骤: 2.5.1 确定指标 在设计一个滤波器之前, 必须首先根据工程实际的需要确定滤波器的技术指标。在很多 6 实际应用中,数字滤波器常常被用来实现选频操作。因此,指标的形式一般在频域中给 出幅度和相位响应。 2.5.2 模型逼近 一旦确定了技术指标, 就可利用已学习过的基本原理和关系式,提出一个滤波器模型来 逼近给定的指标体系。这是滤波器设计所要研究的主要问题。 2.5.3 实现性能分析和计算机仿真 以上两步的结果得到的滤波
13、器,通常是以差分方程、系统函数或脉冲响应描述的。 根据这个描述就可以分析其频率特性和相位特性、以验证设计结果是否满足指标要求; 或者利用计算机仿真实现设计的滤波器,再分析滤波结果来判断。 3.IIR数字滤波器设计的原理与方法 3.1 IIR数字滤波器设计的原理 IIR 数字滤波器具有无限持续时间冲激响应, 需要用递归模型来实现, 其系统函数 为: 7 0 1 ( ) ( ) ( ) 1 M k k k N k k k b z Y z H z X z a z (2.1) 设计 IIR 滤波器的任务就是寻求一个因果、物理上可实现的系统函数H (z) ,使 其频率响应 () jw H e 满足所希望
14、得到的频域指标,即符合给定的通带截止频率、阻带截 止、通带衰减和阻带衰减 . 。 利用冲激响应不变法设计数字滤波器时可能会导致频域混叠现象,为了克服这一问 题,需要找到由 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 ,那么 0 () j H e 1。 双线性变换满足以上4 个条件的映射关系,其变换公式为 21 1 s
15、z s T z (2.2 ) 双线性变换的基本思路: 首先将整个平面压缩到一条从/Ts 变换到 2/Ts 的横带里, 然后通过标准的变换关系 sT ze 将横带变换到整个平面上去,这样就得到 8 了平面与平面间的一一对应的单值关系。 3.2 IIR 数字滤波器设计的基本方法 IIR 数字滤波器的设计方法有两类, 一类是借助于模拟滤波器的设计方法设计出模 拟滤波器 , 利用冲激响应不变法或双线性变换法转换成数字滤波器, 再用硬件或软件实 现; 另一类是直接在频域或时域中进行设计, 设计时需要计算机作辅助工具。 随着 MATLAB 软件尤其是MATLAB 的信号处理工作箱的不断完善, 不仅数字滤波
16、器的计算机辅助设计 有了可能 , 而且还可以使设计达到最优化。 IIR 数字滤波器设计的基本步骤如下: (1) 根据任务 , 确定性能指标。在设计一个滤波器之前 , 首先根据工程实际的需要确 定滤波器的技术指标如:边界频率:p ,s ,c ;阻带最小衰减 As 和通带最大 衰减 Rp ; (2) 将数字滤波器的技术指标转换成模拟滤波器指标。利用冲激响应不变法与双线 性变换法进行频率间的转换,主要是边界频率Wp与 Ws 的转换。 (3) 用模拟滤波器设计方法得到模拟滤波器的传输函数Ha ( s) ;可借助巴特沃 斯(Butterworth) 滤波器、切比雪夫 (Chebyshev) 滤波器、椭圆
17、 (Cauer) 滤波器、贝塞 9 尔(Bessel) 滤波器等,这些滤波器都有严格的设计公式、现成曲线和图表供设计人员 使用。 4. IIR带通滤波器的 MATLAB 设计 4.1 IIR带阻滤波器的设计流程 本文设计的 IIR 带阻滤波器是从低通变换过来的, 利用的是双线性变换以及切比雪 夫 II滤波器的原型,其具体的设计流程为上图所示。首先根据题目要求确定带阻滤波 器的技术指标,先要进行频率的预畸变,并且归一化频率,再设计出切比雪夫II 模拟 低通滤波器, 并求出其阶数等相关参数。 其次利用双线性变换法设计数字带通滤波器, 再调用函数进行双线性变换,并求出分子、分母的系数向量。最后通过画
18、图求出其幅频 响应、相频响应、幅度特性曲线与零极点,并画出波形图。最后进行验证,看所设计的 滤波器能否达到要求的指标,若能达到,则说明该滤波器设计符合要求。 3.2 IIR带阻滤波器的设计步骤 (1)根据设计流程,首先确定所要设计的数字带阻滤波器的相关指标: 通带截止频率 wp1 =0.4,wp2 =0.6 ,通带最大衰减 Rp = 2dB; 阻带截止频率 ws1=0.2 和 ws2=0.8,阻带最小衰减Rs = 30dB; 10 取样间隔 T = 0.1s 。 其实现程序如下(程序中pi 代表) : clear;clc;close all ; S,Fs,Bit = wavread(C:Use
19、rsyyyDesktopccSHE.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(xlabel,X(1:150001);title( 原始信号频谱图 ); % 分析音频 信号,找出干扰频率段为8000Hz subplot(2,1,2);plo
20、t(xlabel,X1(1:150001);title( 加噪信号频谱图 ); % 对于带阻滤波器: fp = 5000; fs = 6000; Rp = 3; Rs = 40; % 求数字频率 wp = 2*pi*fp/Fs; ws = 2*pi*fs/Fs; 11 % 预畸变处理 Wp = 2*Fs*tan(wp/2); Ws = 2*Fs*tan(ws/2); % 开始设计模拟低通滤波器 N,Wn = buttord(Wp,Ws,Rp,Rs,s ); B,A = butter(N,Wn,s ); b,a = bilinear(B,A,Fs); h,w = freqz(b,a,1000,w
21、hole ); 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); set(gca, YTickMode , manual ,
22、YTick,-40,-20,-3,0);grid; 12 S2 = filter(b,a,S1); sound(S2,Fs); X2=abs(fft(S2); figure;plot(xlabel,X2(1:150001);title( 滤波以后的频谱 ) 5.IIR带阻滤波器频谱 00.511.522.5 x 10 4 0 500 1000 1500 2000 原 始 信 号 频 谱 图 00.511.522.5 x 10 4 0 2 4 6 x 10 4 加 噪 信 号 频 谱 图 13 00.51 0 0.2 0.4 0.6 0.8 1 幅 频 响 应 00.20.31 -40 -20
23、-3 0 幅 频 响 应 (dB) 00.511.522.5 x 10 4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 x 10 4 滤 波 以 后 的 频 谱 14 6. IIR带阻滤波器的仿真结果分析 根据上述的仿真波形,可以看出: (1)由 ChebyshevI 型、ChebyshevII 型低通原型变换成带通模型的幅频特性可以 用分贝形式表示。对于ChebyshevI 型的幅频响应在通带内为等纹波衰减、阻带内为单 调减小的,且通带(归一化)截止频率在波形上显示为:0.4 与 0.6 ,阻带(归一化) 截止频率在波形上显示为:0.2 与 0.8 ,与设计要
24、求基本一致。对于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 型带阻 滤波器的幅频响应曲线上看,当N 3 时,在 0 (归一化后为0-1 )范围内通带波 动只有两个波谷,阻带则单调减小,这与理论内容一致。对于Cheb
25、yshevII 型带阻滤 波器来说,其幅频响应曲线的通带波动在 0 (归一化后为 0-1 )范围内阻带波动只 有两个波峰,通带则单调减小,这与理论内容一致。因此,对于滤波器阶数这一指标来 15 说,满足要求。 (3)对于分贝化后的幅度特性, ChebyshevI 型带阻滤波器的通带最大衰减大约为 2dB, 阻带衰减最小为 30dB , 而 ChebyshevII 型带阻滤波器的阻带最小衰减大约为30dB , 通带最小衰减大约为2dB,这与带阻滤波器的技术指标相同。因此,在衰减性能上,所 设计的滤波器达到了要求。 (4)对于相频响应来说, ChebyshevI 型带阻滤波器在通带内的相频响应曲线
26、接近 为一条平滑曲线,在阻带内则存在衰减、畸变。ChebyshevII 型带阻滤波器在通带截止 频率以内的相频响应曲线几乎为一条直线,但在通带截止频率处有很大角度的转折,使 相频特性在通带内产生了的畸变,而在阻带内存在很大的畸变。因此,对于双线性变换 法, 其相位特性得不到满足, 必须用其它方法加以校正或采用其它方法重新设计滤波器, 使之满足线性相位。 (5)群延迟是衡量一个滤波器或是整个系统性能指标的一个重要参数。对于利用 这两种低通滤波器原型模型设计的这两种带通滤波器,带阻滤波器的延迟是同等带宽的 低通滤波器延迟的两倍。这个结果是由低通特性向带通特性变换引起的,n 阶低通滤波 器传递函数总
27、是变换为2n 阶带通滤波器传递函数,其滤波器的群延迟特性的中心频率 与理论上的很接近。 16 总结 此次 Matlab 课程设计是在一定的理论基础之上进行的,在先修课程信号与系统 与数字信号处理中,大量有Matlab 设计方面的知识,所以做起来还比较容易,而 且经过了自己的亲身实践,学到了许多实践方面的知识。 首先,在信号滤波系统中, 有时因为模拟滤波器阶数太高,硬件占用空间太大为某 些仪器的实现设置了障碍, 而对于一些窄带情况下的低通滤波器用模拟手段往往很难实 现。在这些情况下,数字滤波器将会是一个很好的解决办法。MATLAB 信号处理工具箱 提供了丰富而简便的设计、实现FIR 和 IIR
28、的方法,使原来繁琐的程序设计简化成函 数的调用,特别是滤波器的表达方式和滤波器之间的相互转换显得十分简便。 其次,IIR 数字滤波器的设计和模拟滤波器的设计有着紧密的关系。通常要先设计 出适当的模拟滤波器, 再通过一定的频带变换把它转换成为所需的数字IIR 滤波器。此 外,任何数字信号处理系统中也还不可避免地用到模拟滤波器,比如 A/D 变换器前的抗 混叠滤波器和 D/A变换器后的平滑滤波器,因此模拟滤波器设计也是很重要的。 最后,在比较设计滤波器的方法上应该明确其技术指标以及某些参数的实际意义。 比如本文用双线性变换法设计数字带阻滤波器时,必须先将频率归一化, 并且进行频率 预畸变,然后设计
29、模拟滤波器,再利用频率变换法将模拟低通变为模拟带通,最后经过 17 双线性变换法将模拟带通变换为数字带通。如果不进行频率预畸变, 那么设计出来的带 通滤波器的幅频特性与相频特性将会产生很严重的畸变(如图5、图 6 所示,为设计时 没有进行频率预畸变的特性) ,使设计的结果不满足给定的要求,在实际中会造成很严 重的危害。 在做本次课程设计的过程中, 我深深地感受到了自己所学到知识的有限,明白了只 学好课本上的知识是不够的, 要通过图书馆和互联网等各种渠道来扩充自己的知识。在 实验过程中我们曾经遇到过问题。但是从中我们学习到了如何对待遇到的困难,进一步 培养了我们一丝不苟的科学态度和不厌其烦的耐心。所有的这些心得会对我以后的学习 和工作有帮助作用,忠心感谢学校给我们提供这次实验机会。 18 参考文献 1、从玉良 . 数字信号处理原理及其MATLAB 实现M. 北京:电子工业出版社 .2009.7 2、胡广书 . 数字信号处理理论、算法与实现M. 北京:清华大学出版社 .2003,8
链接地址:https://www.31doc.com/p-5197524.html