基于matlab声音信号的滤波去噪处理毕业论文.doc

《基于matlab声音信号的滤波去噪处理毕业论文.doc》由会员分享，可在线阅读，更多相关《基于matlab声音信号的滤波去噪处理毕业论文.doc（23页珍藏版）》请在三一办公上搜索。

1、基于matlab声音信号的滤波去噪处理毕业论文 基于matlab声音信号的滤波去噪处理摘要滤波器设计在数字信号处理中占有极其重要的地位。FIR数字滤波器和IIR滤波器是滤波器设计的重要组成部分。Matlab功能强大、简单易学、编程效率高,深受广大科技工作者的欢迎。特别是Matlab还具有信号分析工具箱,不需具备很强的编程能力,就可以很方便地进行信号分析、处理和设计。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用MATLAB

2、作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。ABSTRACTThe Design of Analysis and Processing Voice Signal Abstract Speech signal processing is to study the use of d

3、igital signal processing technology and knowledge of the voice signal voice processing of the emerging discipline is the fastest growing areas of information science one of the core technology. Transmission of information through the voice of humanitys most important, most effective, most popular an

4、d most convenient form of exchange of information. Matlab language is a data analysis and processing functions are very powerful computer application software, sound files which can be transformed into discrete data files, then use its powerful ability to process the data matrix operations, such as

5、digital filtering, Fourier transform, when domain and frequency domain analysis, sound playback and a variety of map rendering, and so on. Its signal processing and analysis toolkit for voice signal analysis provides a very rich feature function, use of these functions can be quick and convenient fe

6、atures complete voice signal processing and analysis and visualization of signals, makes computer interaction more convenient . Matlab Signal Processing is one of the important areas of application. The design of voice-processing software for most of the content are numerous, easy to maneuver and so

7、 on, using MATLAB7.0 comprehensive use GUI interface design, various function calls to voice signals such as frequency, amplitude, Fourier transform and filtering, the program interface concise, simple, has some significance in practice. Finally, the speech signal processing further development put

8、forward their own views. 目 录摘要 ABSTRACT 绪论 1.1研究的目的和意义 1.2国内外同行的研究状况 1.3本课题的研究内容和方法 语音信号去噪方法的研究 2.1去噪的原理 2.2去噪的方法 去噪和仿真的研究 3.1语音文件在MATLAB平台上的录入与打开 3.2 原始语音信号频谱分析及仿真 3.3 加噪语音信号频谱分析及仿真 3.4 去噪及仿真 3.5 结合去噪后的频谱图对比两种方式滤波的优缺点 总结 致谢 参考文献 1.绪论1.1研究的目的和意义 语音信号的采集与分析技术是一门涉及面很广的交叉科学，它的应用和发展与语音学、声音测量学、电子测量技术以及数字

9、信号处理等学科紧密联系，语音是人类获取信息的重要来源和利用信息的重要手段。在信号传输过程中，由于实验条件或各种其他主观或客观条件的原因，语音处理系统都不可避免地要受到各种噪声的干扰。噪声不但降低了语音质量和语音的可懂度，而且还将导致系统性能的急剧恶化，严重时使整个系统无法正常工作。 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程

10、度上摆脱了传统非交互式程序设计语言，如C、Fortran的编辑模式，代表了当今国际科学计算软件的先进水平。其强大的数据处理能力可以极大程度上削弱噪声影响还原出真实的语音信号相符度在90%以上。 1.2 国内外同行研究现状 20世纪60年代中期形成的一系列数字信号处理的理论和算法，如数字滤波器、快速傅立叶变换（FFT）等是语音信号数字处理的理论和技术基础。随着信息科学技术的飞速发展，语音信号处理取得了重大的进展：进入70年代之后，提出了用于语音信号的信息压缩和特征提取的线性预测技术（LPC），并已成为语音信号处理最强有力的工具，广泛应用于语音信号的分析、合成及各个应用领域，以及用于输入语音与参考

11、样本之间，时间匹配的动态规划方法；80年代初一种新的基于聚类分析的高效数据压缩技术矢量量化（VQ）应用于语音信号处理中；而用隐马尔可夫模型（HMM）描述语音信号过程的产生是80年代语音信号处理技术的重大发展，目前HMM已构成了现代语音识别研究的重要基石。近年来人工神经网络 ANN 的研究取得了迅速发展，语音信号处理的各项课题是促进其发展的重要动力之一，他的各项成果也体现在语音信号处理的各项技术之中。1.3本课题的研究内容和方法 研究内容 本论文主要介绍的是的语音信号的简单处理。本论文针对以上问题，运用数字信号学基本原理实现语音信号的处理，在matlab7.0环境下综合运用信号提取，幅频变换以及

12、傅里叶变换、滤波等技术来进行语音信号处理。我所做的工作就是在matlab7.0软件上编写一个处理语音信号的程序，能对语音信号进行采集，并对其进行各种处理，达到简单的语音信号处理的目的。 运行环境 运行环境主要介绍了硬件环境和软件环境。 硬件环境： 处理器：Inter Pentium B950 内存：2G 硬盘空间：460G 显卡：NVIDIA GeForce GT520 操作系统：Window 7（64位旗舰版） 开发环境 开发环境主要介绍了本系统采用的操作系统、开发语言。 操作系统：Windows 7 2 开发环境：Matlab 7.0骤 语音信号的录制。在MATLAB平台上读入语音信号。绘

13、制频谱图并回放原始语音信号。利用MATLAB编程加入一段正弦波噪音，设计滤波器去噪。利用MATLAB编程加入一段随机噪音信号，设计FIR和IIR滤波器去噪，并分别绘制频谱图、回放语音信号。6 通过仿真后的图像以及对语音信号的回放，对比两种去噪方式的优缺点。具体流程图如下所示： 图1.1论文设计流程2.语音信号去噪方法的研究2.1 去噪的原理 采样定理在进行模拟/数字信号的转换过程中，当采样频率fs.大于信号中，最高频率f的2倍时，即：fs. 2f,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的510倍；采样定理又称奈奎斯特定理。 1924年奈奎斯特

14、 Nyquist 就推导出在理想低通信道的最高大码元传输速率的公式: 理想低通信道的最高大码元传输速率 2W*log2 N 其中W是理想低通信道的带宽,N是电平强度 为什么把采样频率设为8kHz?在数字通信中，根据采样定理, 最小采样频率为语音信号最高频率的2倍频带为F的连续信号f t 可用一系列离散的采样值f t1 ,f t1t ，f t12t ，.来表示,只要这些采样点的时间间隔t1/2F，便可根据各采样值完全恢复原来的信号f t 。 这是时域采样定理的一种表述方式。 时域采样定理的另一种表述方式是：当时间信号函数f t 的最高频率分量为fM时,f t 的值可由一系列采样间隔小于或等于1/

15、2fM的采样值来确定,即采样点的重复频率f2fM。图为模拟信号和采样样本的示意图。 时域采样定理是采样误差理论、随机变量采样理论和多变量采样理论的基础对于时间上受限制的连续信号f t （即当t T时,f t 0,这里T T2-T1是信号的持续时间），若其频谱为F（）,则可在频域上用一系列离散的采样值 采样值来表示,只要这些采样点的频率间隔 采样频率，也称为采样速度或者采样率，定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率的倒数是采样周期或者叫作采样时间，它是采样之间的时间间隔。通俗的讲采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，

16、衡量声卡、声音文件的质量标准。采样频率只能用于周期性采样的采样器，对于非周期性采样的采样器没有规则限制。 采样频率的常用的表示符号是 fs。 通俗的讲采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位时间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。采样频率与声音频率之间有一定的关系，只有采样频率高于声音信号最高频率的两倍时，才能把数字信号表示的声音还原成为原来的声音。这就是说采样频率是衡量声卡采集、记录和还原声音文件的质量标准。采样位数和采样率对于音频接口来说是最为重要的两个指标，也

17、是选择音频接口的两个重要标准。无论采样频率如何，理论上来说采样的位数决定了音频数据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dB。采样位数越多则捕捉到的信号越精确。对于采样率来说你可以想象它类似于一个照相机，44.1kHz意味着音频流进入计算机时计算机每秒会对其拍照达441000次。显然采样率越高，计算机摄取的图片越多，对于原始音频的还原也越加精确在我们的日常交流和语音通信系统中，加性宽带噪声严重影响了语音质量和可懂度。从带噪语音中提取原始语音信号的方法很多，在单信道条件下，谱相减算法以其运算量小、原理简单、易于实现并且有不错的增强效果而得到了广泛的应用。谱相减语音增强算法的核心

18、是噪声检测和谱减规则。在分析了语音增强算法理论的基础上，本文首先研究了语音激活检测算法。对基于短时能量和短时过零率双门限法语音激活检测的噪声估计算法做了研究及仿真，同时还研究了一种基于最小子带能量的噪声估计方法。然后，通过分析经典谱减法的原理及其一般改进形式，研究了一种基于噪声残差的谱相减改进算法和一种可以不以噪声是零均值的高斯分布为前提的谱减法改进算法。最后通过大量的仿真实验，验证了所研究的几种改进算法都能有效地提高增强效果。在进入FIR滤波器前，首先要将信号通过A/D器件进行模数转换，使之成为8bit的数字信号，一般可用速度较高的逐次逼进式A/D转换器，不论采用乘累加方法还是分布式算法设计

19、FIR滤波器，滤波器输出的数据都是一串序列，要使它能直观地反应出来，还需经过数模转换，因此由FPGA构成的FIR滤波器的输出须外接D/A模块。FPGA有着规整的内部逻辑阵列和丰富的连线资源，特别适合于数字信号处理任务，相对于串行运算为主导的通用DSP芯片来说，其并行性和可扩展性更好，利用FPGA乘累加的快速算法，可以设计出高速的FIR数字滤波器。有限长单位冲激响应（FIR）滤波器有以下特点： 位冲激响应h n 在有限个n值处不为零 0处收敛，极点全部在z 0处结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中（例如频率抽样结构）也包含有反馈的递归部分。 设FIR滤波器的单位冲激响应h

20、n 为一个N点序列，0 n N 1，则滤波器的系统函数为 H z h n *z-n （2-3）就是说，它有（N1）阶极点在z 0处，有（N1）个零点位于有限z平面的任何位置。FIR滤波器基本结构FIR滤波器有以下几种基本结构： 横截型式的系统的差分方程表达式为 y n h m x n-m 2-4）很明显，这就是线性移不变系统的卷积和公式，也是x n 的延时链的横向结构，称为横截型结构或卷积型结构，也可称为直接型结构。将转置定理用于可得到的转置直接型结构。 FIR滤波器的横截型结构 级联型其中N/2表示取N/2的整数部分。若N为偶数，则N1为奇数，故系数B2K中有一个为零，这是因为，这时有奇数个

21、根，其中复数根成共轭对必为偶数，必然有奇数个实根。画出N为奇数时，FIR滤波器的级联结构，其中每一个二阶因子用图4-11的横型结构。 这种结构的每一节控制一对零点，因而再需要控制传输零点时，可以采用它。但是这种结构所需要的系数B2k（I 0，1，2，k， 1，2，N/2）比卷积型的系数h n 要多，因而所需的乘法次数也比卷积型的要多。 nfinite Impulse Response 数字滤波器，又名“无限脉冲响应数字滤波器”，或“递归滤波器”。递归滤波器，也就是IIR数字滤波器，顾名思义，具有反馈，一般认为具有无限的脉冲响应。IIR数字滤波器的设计利用分析工具 FDATool 可以很方便地设

22、计出符合应用要求的未经量化的IIR数字滤波器。需要将MATLAB设计出的IIR数字滤波器进一步分解和量化，从而获得可用FPGA实现的滤波器系数。 IIR数字滤波器的设计步骤由于采用了级联结构，因此如何将滤波器的每一个极点和零点相组合，从而使得数字滤波器输出所含的噪声最小是个十分关键的问题。为了产生最优的量化后的IIR数字滤波器，采用如下步骤进行设计。 首先计算整体传递函数的零极点； 选取具有最大幅度的极点以及距离它最近的零点，使用它们组成一个二阶基本节的传递函数； 对于剩下的极点和零点采用与相类似的步骤，直至形成所有的二阶基本节。 通过上面三步法进行的设计可以保证IIR数字滤波器中N位乘法器产

23、生的量化舍入误差最小。 获得最优IIR数字滤波器系数为了设计出可用FPGA实现的数字滤波器，需要对上一步分解获得的二阶基本节的滤波器系数进行量化，即用一个固定的字长加以表示。量化过程中由于存在不同程度的量化误差，由此会导致滤波器的频率响应出现偏差，严重时会使IIR滤波器的极点移到单位圆之外，系统因而失去稳定性。为了获得最优的滤波器系数，采用以下步骤进行量化。 计算每个系数的B.查找出每个系数绝对值中的最大值； 计算比此绝对值大的最小整数； 取反获得负整数； 计算需要表示此整数的最小位数； 计算用于表示系数值分数部分的余下位数。 除了系数存在量化误差，数字滤波器运算过程中有限字长效应也会造成误差

24、，因此对滤波器中乘法器、加法器及寄存器的数据宽度要也进行合理的设计，以防止产生极限环现象和溢出振荡。 与FIR数字滤波器的设计不同，IIR滤波器设计时的阶数不是由设计者指定，而是根据设计者输入的各个滤波器参数（截止频率、通带滤纹、阻带衰减等），由软件设计出满足这些参数的最低滤波器阶数。在MATLAB下设计不同类型IIR滤波器均有与之对应的函数用于阶数的选择。 ”。利用MATLAB中的wavread命令来读入（采集）语音信号，将它赋值给某一向量。Wavread函数，可得出信号的采样频率为22500，并且声音是单声道的。利用Sound函数可以清晰的听到“毕业设计声音文件”的语音。采集数据并画出波形

25、图。在MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。下面介绍这些函数。函数FFT用于序列快速傅立叶变换。函数的一种调用格式y fft x 。其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT。且和x相同长度。若x为一矩阵，则y是对矩阵的每一列向量进行FFT。如果x长度是2的幂次方，函数fft执行高速基2FFT算法；否则fft执行一种混合基的离散傅立叶变换算法，计算速度较慢。函数FFT的另一种调用格式为y fft x,N 式中，x，y意义同前，N为正整数。函数执行N点的FFT。若x为向量且长度小于N，则函数将x补零至长度N。若向量

26、x的长度大于N，则函数截短x使之长度为N。若x 为矩阵，按相同方法对x进行处理。经函数fft求得的序列y一般是复序列，通常要求其幅值和相位。MATLAB提供求复数的幅值和相位函数：abs，angle，这些函数一般和 FFT同时使用。函数abs x 用于计算复向量x的幅值，函数angle x 用于计算复向量的相角，介于和之间，以弧度表示。用MATLAB工具箱函数fft进行频谱分析时需注意：（1）函数fft返回值y的数据结构对称性一般而言，对N点的x n 序列的FFT是N点的复数序列，其点n N/2+1对应Nyquist频率，作频谱分析时仅取序列X k 的前一半，即前N/2点即可。X k 的后一半

27、序列和前一半序列时对称的。频率计算 若N点序列x n n 0,1,N-1 是在采样频率 下获得的。它的FFT也是N点序列，即X k k 0,1,2,N-1 ，则第k点所对应实际频率值为f k*f /N.绘出了语音信号的波形频谱图。x,fs,bits wavread wangqingtian.wav ;sound x,fs,bits ;X fft x,4096 ;magX abs X ;angX angle X ;subplot 221 ;plot x ;title 原始信号波形 ;subplot 222 ;plot X ; title 原始语音信号采样后的频谱图 subplot 223 ;pl

28、ot magX ;title 原始信号幅值 ; subplot 224 ;plot angX ;title 原始信号相位 ;程序运行可以听到声音，得到的图形为：atlab函数randn：产生正态分布的随机数或矩阵的函数产生均值为0，方差 2 1，标准差 1的正态分布的随机数或矩阵的函数。用法：Y randn n 返回一个n*n的随机项的矩阵。如果n不是个数量，将返回错误信息。Y randn m,n 或 Y randn m n 返回一个m*n的随机项矩阵。Y randn size A 返回一个和A有同样维数大小的随机数组。randn返回一个每次都变化的数量y,fs,bits wavread ch

29、ushi.wav ;sound y,fs n length y y_p fft y,n ;f fs* 0:n/2-1 /n;figure 1 subplot 2,1,1 ;plot y ;title 原始语音信号采样后的时域波形 ;xlabel 时间轴 ylabel 幅值A subplot 2,1,2 ;plot f,abs y_p 1:n/2 ;title 原始语音信号采样后的频谱图 ;xlabel 频率Hz ;ylabel 频率幅值 ;L length y noise 0.1*randn L,2 ;y_z y+noise;sound y_z,fs n length y ;y_zp fft

30、y_z,n ;f fs* 0:n/2-1 /n;figure 2 subplot 2,1,1 ;plot y_z ;title 加噪语音信号时域波形 ;xlabel 时间轴 ylabel 幅值A subplot 2,1,2 ;plot f,abs y_zp 1:n/2 ;title 加噪语音信号频谱图 ;xlabel 频率Hz ;ylabel 频率幅值 ;初始信号的时域波形、频谱图与加噪后语音信号的时域波形、频谱图分别如图3.5、3.6所示。图3.5 图3.6通过对两张图片的对比，很明显可以看加噪后的语音信号时域波形比原始语音信号浑浊了许多，在时间轴上可以明显看出00.5S的幅值增大了；通过对

31、原始语音信号的频谱图与加噪后的语音信号频谱图的对比，也可以看出在频率5000Hz以后的频率幅值发生了明显的增加。再通过对原始语音信号的回放效果与加噪后的语音信号回放的效果的对比，人耳可以明显辨别出两种语音信号不一样了，加噪后的语音信号在听觉上比原始语音信号要浑浊很多，而且还有吱吱嘎嘎的混杂音。3.4 去噪及仿真（1）FIR滤波器法去噪通过对上一节中加噪语音信号和原始语音信号频谱图对比可以知道，噪音大部分是Hz大于5000的部分，故设计低通滤波器进行滤波处理。接下来我们要用设计的FIR低通滤波器对上一节中加噪语音信号进行滤波处理。用自己设计的FIR数字低通滤波器对加噪的语音信号进行滤波时，在Ma

32、tlab中，FIR滤波器利用函数fftfilt对信号进行滤波。函数fftfilt用的是重叠相加法实现线性卷积的计算。调用格式为：y fftfilter h,x,M 。其中，h是系统单位冲击响应向量；x是输入序列向量；y是系统的输出序列向量；M是有用户选择的输入序列的分段长度，缺省时，默认的输入向量的重长度M 512。用设计好的FIR数字低通滤波器对加噪语音信号的滤波程序：y,fs,bits wavread chushi.wav ;sound y,fs n length y y_p fft y,n ;f fs* 0:n/2-1 /n;figure 1 subplot 2,1,1 ;plot y

33、;title 原始语音信号采样后的时域波形 ;xlabel 时间轴 ylabel 幅值A subplot 2,1,2 ;plot f,abs y_p 1:n/2 ;title 原始语音信号采样后的频谱图 ;xlabel 频率Hz ;ylabel 频率幅值 ;L length y noise 0.1*randn L,2 ;y_z y+noise;sound y_z,fs n length y ;y_zp fft y_z,n ;f fs* 0:n/2-1 /n;figure 2 subplot 2,1,1 ;plot y_z ;title 加噪语音信号时域波形 ;xlabel 时间轴 ylabel

34、 幅值A subplot 2,1,2 ;plot f,abs y_zp 1:n/2 ;title 加噪语音信号频谱图 ;xlabel 频率Hz ;ylabel 频率幅值 ;Ft 5000;Fp 1000;Fs 1200;wp 2*Fp/Ft;ws 2*Fs/Ft;rp 1;rs 50;p 1-10. -rp/20 ;s 10. -rs/20 ;fpts wp ws;mag 1 0;dev p s;n21,wn21,beta,ftype kaiserord fpts,mag,dev ;b21 fir1 n21,wn21,Kaiser n21+1,beta ;h,w freqz b21,1 ;pl

35、ot w/pi,abs h ;title FIR低通滤波器 ;x fftfilt b21,y_z ;X fft x,n ;figure 4 ;subplot 2,2,1 ;plot f,abs y_zp 1:n/2 ;title 滤波前信号的频谱 ;subplot 2,2,2 ;plot f,abs X 1:n/2 ;title 滤波后信号的频谱 ;subplot 2,2,3 ;plot y_z ;title 滤波前信号的时域波形 subplot 2,2,4 ;plot x ;title 滤波后信号的时域波形 sound x,fs,bits 而后得到图像如下：图3.7分析：从以上四图可以很明显

36、和直观的看出原始语音信号和加噪语音信号时域波形和频谱图的区别。加噪后的语音信号的时域波形比原始语音信号要模糊得多，频谱图则是在频率5000Hz以后出现了明显的变化。再通过滤波前的信号波形和频谱图的对比，可以明显看出滤波后的波形开始变得清晰了，有点接近原始信号的波形图了。滤波后信号的频谱图也在5000Hz以后开始逐渐接近原始语音信号的频谱图。再从对语音信号的回放，人耳可以明显辨别出加噪后的语音信号比较浑浊，还有很明显嘎吱嘎吱的杂音在里面。滤波后，语音信号较加噪后的信号有了明显的改善，基本可以听清楚了，而且杂音也没有那么强烈，但是声音依然没有原始语音信号那么清晰脆耳。（2）IIR滤波器法去噪同样，

37、也设计一个IIR低通滤波器对加噪语音信号进行内部处理，程序如下：Ft 8000;Fp 1000;Fs 1200;wp 2*pi*Fp/Ft;ws 2*pi*Fs/Ft;fp 2*Ft*tan wp/2 ;fs 2*Fs*tan wp/2 ;n11,wn11 buttord wp,ws,1,50,s ; %求低通滤波器的阶数和截止频率b11,a11 butter n11,wn11,s ; %求S域的频率响应的参数num11,den11 bilinear b11,a11,0.5 ; %双线性变换实现S域到Z域的变换h,w freqz num11,den11 ; %根据参数求出频率响应plot w*

38、8000*0.5/pi,abs h ;legend 用butter设计 ;grid;y,fs,nbits wavread BYSJ.wav ;n length y ; %求出语音信号的长度noise 0.01*randn n,2 ; %随机函数产生噪声s y+noise; %语音信号加入噪声S fft s ; %傅里叶变换z11 filter num11,den11,s ;sound z11 ;m11 fft z11 ; %求滤波后的信号subplot 2,2,1 ;plot abs S ,g ;title 滤波前信号的频谱 ;grid;subplot 2,2,2 ;plot abs m11

39、,r ;title 滤波后信号的频谱 ;grid;subplot 2,2,3 ;plot s ;title 滤波前信号的波形 ;grid;subplot 2,2,4 ;plot z11 ; title 滤波后的信号波形 ; 得到图像如下：图3.8 通过程序运行，我们可以听出，滤波效果很不好，失真现象严重，并且出现了严重的寄生震荡。3.5 结合去噪后的频谱图对比两种方式滤波的优缺点IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断

40、累积，有时会产生微弱的寄生振荡。 （1）IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。FIR滤波器则要求较低。（2）IIR滤波器运算误差大，有可能出现极限环振荡，FIR相比之下运算误差较小，不会出现极限环振荡。（3）IIR幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上； （4）与FIR滤波器的设计不同，IIR滤波器设计时的阶数不是由设计者指定，而是根据设计者输入的各个滤波器参数（截止频率、通带滤纹、阻带衰减等），由软件设计出满足这些参数的最低滤波器阶数。在MATLAB下设计不同类型IIR滤波器均有与之对应的函数用于阶数的选择。 （5）IIR单

41、位响应为无限脉冲序列FIR单位响应为有限的 （6）FIR幅频特性精度较之于iir低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变。这是很好的性质。 （7）IIR滤波器有噪声反馈，而且噪声较大，FIR滤波器噪声较小。FIR幅频特性精度较之于iir低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变。这是很好的性质。4.总结本文对语音信号处理系统的设计作了详细的介绍，采用一系列图像分析和处理技术，实现了语音信号的基本处理的功能，经过测试运行，本设计圆满的完成了对语音信号的读取与打开较好的完成了对语音信号的频谱分析，通过fft变换，得出了语音信号的频谱

42、图；在滤波这一块，课题主要是从滤波器入手来设计滤波器，基本实现了滤波与课题的要求十分相符在此论文撰写过程中，要特别感谢我的导师的指导与督促，同时感谢的谅解与包容。没有的帮助也就没有今天的这篇论文。求学历程是艰苦的，但又是快乐的。感谢我的老师，谢谢他在这四年中为我们全班所做的一切，他不求回报，无私奉献的精神很让我感动，再次向表示由衷的感谢。在这四年的学期中结识的各位生活和学习上的挚友让我得到了人生最大的一笔财富。在此，也对他们表示衷心感谢。本文参考了大量的文献资料，在此，向各学术界的前辈们致敬！高等教育出版社5皇甫堪,陈建文,楼生强.现代数字信号处理M 北京电子工业出版社2003. 11蔡旭辉，刘卫国，蔡立燕.MATLAB基础与应用教程 M .人民邮电出版社.2009.12尚涛，谢龙汗，杜如虚.MATLAB工程计算及分析 M .清华大学出版社.2011.宿州学院毕业论文设计语音信号采集信号加噪语音信号录入语音信号变换语音信号滤波效果显示对比- 1 - 11 -