语音信号的滤波处理八：.doc

1 设计目的及内容 1.1 设计目的 (1)掌握数字信号处理的基本概念，基本理论和基本方法。 (2)熟悉离散信号和系统的时域特性。 (3)掌握序列快速傅里叶变换方法。 (4)学会MATLAB的使用，掌握MATLAB的程序设计方法。 (5)掌握利用MATLAB对语音信号进行频谱分析。 (6)掌握滤波器的网络结构。 (7)掌握MATLAB设计IIR、FIR， A为幅度因子，a为衰减因子，为模拟角频率，其中n为学号(例如，王墨同学n=23) ② 要求输入采样频率fs（根据程序处理需要指定范围）后，在时域演示信号波形、采样脉冲及采样后信号；在频域演示不同采样频率下对应信号的频谱。 1.2.2 设计题部分 (1)选择Windows系统的C:\WINDOWS\Media\Windows XP 关机.wav，并对该信号进行采样； (2)语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图； (3)产生噪声信号加到语音信号中，得到被污染的语音信号，并回放语音信号； (4)污染信号的频谱分析，画出被污染的语音信号时域波形和频谱； (5)根据有关的频谱特性，采用间接法设计IIR数字滤波器，并画出相应滤波器的幅频图（设计3个IIR滤波器） a. 模拟滤波器类型：切比雪夫滤波器（低通、带通、高通） b．总体要求：MATLAB原程序+仿真波形+技术指标 (6)根据有关的频谱特性，采用直接法设计FIR数字滤波器，并画出相应滤波器的幅频，相频图（设计3个FIR滤波器） a. 滤波器类型：Blackman窗（低通、带通、高通） b．总体要求：MATLAB原程序+仿真波形+技术指标+窗函数 (7)用自己设计的这些滤波器分别对被不同噪声污染的信号进行滤波； (8)分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化； (9)回放语音信号。 2 课程设计基本原理 2.1设计思路 (1)语音信号声音选自XP系统盘目录下C:\WINDOWS\Media\Windows XP 关机.wav (2)根据不同的滤波器类型选择不同（频率）的噪声信号加到语音信号中，得到被污染的语音信号。 (3)分别设计切比雪夫（低通、带通、高通）以及Blackman窗滤波器（低通、带通、高通）对被污染的语音信号滤波，滤掉相应的噪音信号，得到符合要求的语音信号。 如图2.1所示图为切比雪夫滤波器或Blackman窗滤波器设计流程： 图2.1 总体方案设计流程图 2.2 设计原理 设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统，并使系统函数H(z)具有指定的频率特性。数字滤波器从实现的网络结构或者从单位冲激响应分类，可以分成无限长单位冲激响应（IIR）数字滤波器和有限长单位冲激响应（FIR）数字滤波器。本次数字信号处理课程设计选用切比雪夫数字滤波器（IIR）以及Blackman窗数字滤波器（FIR）。 2.2.1 IIR数字滤波器的结构与设计 IIR滤波器系统函数的极点可以在单位圆内的任何位置，实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，又由于IIR数字滤波器能够保留一些模拟滤波器的优良特性，因此应用很广。设计IIR数字滤波器的方法主要有基于冲激响应不变法的IIR数字滤波器设计，基于双线性Z变换法的IIR数字滤波器设计，数字高通及带通IIR滤波器设计，基于MATLAB函数直接设计IIR数字滤波器。 2.2.1.1 基于双线性Z变换法的IIR数字低通滤波器设计 利用冲激响应不变法设计数字滤波器时会产生频域混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服这一问题，需要找到由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那么，双线性Z变换满足以上4个条件的映射关系，其变换公式为： 双线性Z变换的基本思路是，首先将整个s平面压缩到一条从–π/Ts到π/Ts的带宽为2π/Ts的横带里，然后通过标准的变换关系将横带变换成整个z平面上去，这样就得到s平面与z平面间的一一对应的单值关系。 在MATLAB中双线性Z变换可以通过bilinear函数实现，其调用格式为[Bz,Az]=bilinear(B,A,Fs)。 其中B,A为模拟滤波器传递函数G(s)的分子分母多项式的系数向量，而Bz,Az为数字滤波器的传递函数H(z)的分子分母多项式的系数向量。 2.2.1.2 数字高通及带通IIR滤波器设计 除了低通