基于语音信号去噪处理的FIR滤波器设计要点.pdf
《基于语音信号去噪处理的FIR滤波器设计要点.pdf》由会员分享,可在线阅读,更多相关《基于语音信号去噪处理的FIR滤波器设计要点.pdf(24页珍藏版)》请在三一文库上搜索。
1、摘要 本次课程设计的题目是基于语音信号去噪处理的FIR 滤波器设计,主要是 利用 MATLAB 设计一个 FIR 数字带通滤波器,对一段噪声环境下的语音信号进行 滤波。在设计过程中, 首先要录制一段语音信号, 并对录制好的信号进行时域和 频域分析; 然后对原始的语音信号进行加噪处理,对加噪后的信号进行分析, 绘 出时域和频域的分析图; 最后利用设计出的FIR 带通滤波器,针对语音信号的性 质选取一种适合的窗函数设计滤波器进行滤波,最后对仿真结果进行分析。 关键词 : FIR 带通滤波器;语音信号;MATLAB 仿真;加噪;滤波 目录 前言 1 第一章 基本原理 . 2 1.1 语音处理中的采样
2、原理 2 1.2 数字滤波器的设计 2 1.3 窗函数法 4 第二章 FIR 带通滤波器设计. 6 2.1 FIR 滤波器简介 6 2.2 FIR 带通滤波器设计要求 7 2.3 设计方法 7 2.4 设计步骤 9 第三章 FIR 带通滤波器的软件仿真 10 3.1 程序流程图 . 10 3.2 仿真结果及分析 . 11 总结. 15 参考文献 . 16 附录. 17 致谢. 22 一、 1 前言 在信号处理过程中,所处理的信号往往混有噪音,从接收到的信号中 消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪 音信号的不同特性,提取有用信号的过程称为滤波。实现滤波功能的系统 被称
3、为滤波器。在近代电信设备和各类控制系统中,数字滤波器应用极为 广泛。 语音处理是最早应用数字滤波器的领域之一,也是最早推动数字信号 处理理论发展的领域之一。该领域主要包括5 个方面的内容:第一,语音 信号分析。即对语音信号的波形特性,统计特性,模型参数等进行分析计 算。第二,语音合成。即利用专用硬件或在通用计算机上运行软件来产生 语音。第三,语音识别。即利用专用硬件或计算机识别人的讲话,或者识 别说话的人。第四,语音增强。即从噪音或者干扰中提取被掩盖的语音信 号。第五,语音编码。主要用于语音数据的压缩,目前已经建立了一系列 语音编码国际标准,大量用于通信和音频处理。 MATLAB 是一种面向科
4、学和工程计算的语言,它集数值分析、矩阵运算、信 号处理和图形显示于一体, 具有编程效率高、 调试手段丰富、 扩充能力强等特点。 MATLAB 的信号处理工具箱具有强大的函数功能,它不仅可以用来设计数字滤波 器,还可以使设计达到最优化,是数字滤波器设计的强有力工具。 2 第一章基本原理 1.1 语音处理中的采样原理 因为录制的语音信号是模拟信号, 要想使用数字滤波器对叠加了噪声的信号 进行滤波,则在设计数字滤波器之前首先要进行模数转换,将模拟信号转换为数 字信号。 在进行模数转换的过程中,当最高采样频率fs 大于信号中最高频率f 的 2 倍时,即: fsmax2fmax, 采样之后的数字信号可以
5、完整地保留原始信号中的信 息,一般实际应用中保证采样频率为信号最高频率的510倍;采样定理又称奈 奎斯特定理。 频带为 F 的连续信号 f(t)可用一系列离散的采样值f(t1),f(t1t) , f(t1 2t) ,. 来表示, 只要这些采样点的时间间隔t 1/2F,便可根据各采 样值完全恢复原来的信号f(t)。这是时域采样定理的一种表述方式。时域采样 定理的另一种表述方式是:当时间信号函数f(t)的最高频率分量为fm 时,f(t) 的值可由一系列采样间隔小于或等于1/2fm 的采样值来确定 , 即采样点的重复频 率 f 2fm。 采样频率,也称为采样速度或者采样率, 定义了每秒从连续信号中提
6、取并组 成离散信号的采样个数,它用赫兹(Hz)来表示。采样频率的倒数是采样周期或 者叫采样时间, 它是采样之间的时间间隔。 通俗的讲采样频率是指计算机每秒钟 采集多少个声音样本,是描述声音文件的音质、音调,衡量声卡、声音文件的质 量标准。 采样位数和采样率对于音频接口来说是最为重要的两个指标,也是选择音频 接口的两个重要标准。 无论采样频率如何, 理论上来说采样的位数决定了音频数 据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dB。采样位数 越多则捕捉到的信号越精确 1 。 1.2 数字滤波器的设计 数字滤波器 (Digital Filter,简称为 DF)是指用来对输入信号进行滤
7、 波的硬件和软件。所谓数字滤波器,是指输入、输出均为数字信号,通过一定运 3 算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数 字滤波器和模拟滤波器相比, 因为信号的形式和实现滤波的方式不同,数字滤波 器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等 优点。一般用两种方法来实现数字滤波器: 一是采用通用计算机,把滤波器所要 完成的运算编成程序通过计算机来执行,也就是采用计算机软件来实现; 二是采 用实际专用的数字处理硬件。 数字滤波器根据其冲击响应函数的时域特性,可分为两种,即无限冲击响应 IIR 滤波器和有限冲击响应FIR 滤波器。 IIR 滤波器
8、的特性是,具有无限持续时 间冲击响应。 这种滤波器一般需要用递归模型来实现,因而有时也称之为递归滤 波器。FIR 滤波器的冲击响应只能延续一定的时间,在工程实际中可以采用递归 的方法来实现,也可以采用非递归的方式来实现。 数字滤波器的设计方法有多种, 如双线性变换法,窗函数法,插值逼近法和切比雪夫逼近法等等。随着MATLAB 软件尤其是 MATLAB 信号处理工作的不断完善,不仅数字滤波器的计算机辅助设 计有了可能,而且还可以使设计达到最优化 2 。 数字滤波器的设计基本步骤如下: (1)确定指标 在设计一个滤波器之前, 必须首先根据工程实际的需要确定滤波器的技术指 标。在很多实际应用中,数字
9、滤波器常常被用来实现选频操作。因此,指标的形 式一般在频域中给出幅度和相位响应。幅度指标主要用以下两种方式给出。第一 种是绝对指标。它提供对幅度响应函数的要求,一般应用于FIR 滤波器的设计。 第二种指标是相对指标。 它以分贝值的形式给出要求。 运用线性相位响应指标进 行滤波器设计具有如下优势: 只包含实数算法, 不涉及复数运算。 不存在延 时失真,只有固定数量的延迟。长度为N的滤波器,计算量为N/2 数量级。 (2)逼近 确定了技术指标后, 就可以建立一个目标的数字滤波器模型。通常采用理想 的数字滤波器模型。 之后,利用数字滤波器的设计方法,设计出一个实际滤波器 模型来逼近给定的目标。 (3
10、)性能分析和计算机仿真 4 上两步的结果是得到以差分或者系统函数或者冲击响应描述的滤波器。根据 这个描述就可以分析器频率特性和相位特性,以验证设计结果是否满足指标要 求:或者利用计算机仿真实现设计的滤波器,再分析滤波器结果来判断 3 。 1.3 窗函数法 数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间 序列转化为代表输出信号的数字时间序列,并在转化过程中, 使信号按预定的形 式变化。数字滤波器有多种分类, 根据数字滤波器冲激响应的时域特征,可将数 字滤波器分为两种,即无限长冲激响应(IIR )滤波器和有限长冲激响应(FIR) 滤波器。FIR 数字滤波器的单位脉冲响应是有限长
11、序列。它的设计问题实质上是 确定能满足所要求的转移序列或脉冲响应的常数问题,设计方法主要有窗函数 法、频率采样法和等波纹最佳逼近法等 2 。 用窗函数设计滤波器首先要对滤波器提出性能指标。一般是给定一个理想的 频率响应,使所设计的 FIR 滤波器的频率响应去逼近所要求的理想的滤波器的响 应。窗函数法设计的任务在于寻找一个可实现(有限长单位脉冲响应) 的传递函 数。用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列 4 。主 要设计步骤为: (1)通过傅里叶逆变换获得理想滤波器的单位脉冲响应)(n hd 。 从时域出发, 截取有限长的一段冲击响应作为H(z) 的系数,冲击响应长度 N
12、就是系统函数H(z) 的阶数。只要N 足够长,截取的方法合理,总能满足频域的 要求。一般这种时域设计、 频域检验的方法要反复几个回合才能成功。要设计一 个线性相位的FIR 数字滤波器,首先要求理想频率响应 )( jw d eH 。 )( jw d eH 是 w 的周期函数,周期为 2 ,可以展开成傅氏级数: eheH nj d j d n)()( (1-1) 其中 )(nhd 是与理想频响对应的理想单位抽样响应序列。但不能用来作为设 计 FIR DF用的 h(n) ,因为 )(nhd 一般都是无限长、非因果的,物理上无法实现。 为了设计出频响类似于理想频响的滤波器,可以考虑用h(n) 来近似
13、)(nhd 5 。 窗函数的基本思想: 先选取一个理想滤波器 (它的单位抽样响应是非因 5 果、无限长的),再截取(或加窗)它的单位抽样响应得到线性相位因果FIR 滤 波器,这种方法的重点是选择一个合适的窗函数和理想滤波器 6 。 (2)由性能指标确定窗函数W(n)和窗口长度 N 。 设 x(n) 是一个长序列, w(n) 是长度为 N的窗函数,用 w(n)截断 x(n) ,得到 N点序列 xn(n) ,即 )()()(nnxn xn (1-2) 在频域上则有 jjj de 2 1 eWeXX N (1-3) (3) 求得实际滤波器的单位脉冲响应h(n) , h(n) 即为所设计 FIR 滤波
14、器系 数向量。 )()()(nnnh hd (1-4) 由此可见,窗函数w(n)不仅仅会影响原信号x(n) 在时域上的波形,而且也 会影响到频域内的形状 4 。 6 第二章 FIR 带通滤波器设计 2.1 FIR 滤波器简介 数字滤波器包括 FIR(有限单位脉冲响应)滤波器与IIR (无限单位脉冲响 应)滤波器两种。在现代信号处理技术中,例如数据传输、雷达接收以及一些要 求较高的电子系统, 都越来越多地要求信道具有线性的相位特性。在这方面,FIR 滤波器具有独到的优点, 它可以在幅度特性随意设计的同时,保证精确、 严格的 线性相位特性 5 。 FIR 滤波器的单位脉冲响应h (n) 是有限长的
15、(0nN-1) , 其 z 变换为 1 z 的(N-1)阶多项式: z n N on nh zX zY zH 1 )( )( )( )( (2-1 ) 可得 FIR 滤波器的系统差分方程为: ) 1()1() 1()1()0()0()(NnxNbnxbxbny )()()()( 1 0 nxnbmnxmb N m (2-2) 因此, FIR 滤波器又称为卷积滤波器。FIR 滤波器的频率响应表达式为: 1 0 )()( N n njj enheH (2-3) 信号通过FIR 滤波器不失真条件是在通带内具有恒定的幅频特性和线性相位特 性。理论上可以证明:当FIR 滤波器的系数满足下列中心对称条件:
16、 )1()(nNhnh(2-4) 或者)1()(nNhnh(2-5) 时,滤波器设计在逼近平直幅频特性的同时,还能获得严格的线性相位特性。线 性相位 FIR 滤波器的相位滞后和群延迟在整个频带上是相等且不变的。对于一个 N 阶的线性相位 FIR 滤波器,群延迟为常数, 即滤波后的信号简单地延迟常数个 时间步长。这一特性使通带频率内信号通过滤波器后仍保持原有波形形状而无相 7 位失真 7 。 FIR 滤波器的设计任务是选择有限长度的h(n) ,使传输函数)( j eH满足技 术要求。 FIR 滤波器的设计方法有多种,如窗函数法、频率采样法及其它各种优 化设计方法,本次设计使用窗函数法设计FIR
17、带通滤波器。 2.2 FIR 带通滤波器设计要求 利用 MATLAB 仿真软件系统结合窗函数法设计一个数字带通FIR 滤波器。 设计要求: (1)用所设计的滤波器对受噪声影响的信号进行滤波,画出滤波后语音信号的 时域波形图和频谱图; (2)对滤波前后的信号进行对比,分析信号的变化;回放语音信号,并与原始 语音信号对比。 技术指标: 低端阻带截止频率 fc11000 Hz 低端通带截止频率 fb11200 Hz 高端通带截止频率 fb23000 Hz 高端阻带截止频率 fc23200 Hz 通带衰减系数 ap1dB 阻带衰减系数 as100 dB 2.3 设计方法 设计 FIR 数字滤波器的方法
18、通常有三种:窗函数法, 频率抽样法, 等纹波逼 近法。本次课程设计讨论的是第一种窗函数法。这种方法也叫傅里叶级数法。 一般是先给定所要求的理想滤波器频率响应)( jw d eH,导出)(nhd,我们知 道理想滤波器的冲击响应)(nhd是无限长的非因果序列, 而我们要设计的是h(n) 是有限长的 FIR 滤波器,所以要用有限长序列h(n) 来逼近无限长序列)(nhd, 设: dweeHn jwnjw d h )( 2 1 )((2-6) 8 常用的方法是使用有限长的窗函数w(n) 来截取)(nhd即: )()()(nnwnh hd (2-7) 根据在时域是相乘关系,在频域则是卷积关系: deWe
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 语音 信号 处理 FIR 滤波器 设计 要点
链接地址:https://www.31doc.com/p-5207651.html