数字信号处理期末实验-语音信号分析与处理.docVIP

数字信号处理期末实验-语音信号分析与处理 山东建筑大学信电学院课程设计说明书 PAGE 5 语音信号分析与处理 摘 要 用MATLAB对语音信号进行分析与处理，采集语音信号后，在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。 数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换（FT）。离散傅立叶变换（DFT）和数字滤波是数字信号处理的最基本内容。 关键词：MATLAB;语音信号；加入噪声；滤波器；滤波 1. 设计目的与要求 （1）待处理的语音信号是一个在20Hz~20kHz频段的低频信号。 （2）要求MATLAB对语音信号进行分析和处理，采集语音信号后，在MATLAB平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器进行滤除噪声，恢复原信号。 析离散信号和系统的的有力工具。对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。 （5）数字信号滤波器原理和方法： IIR数字滤波器系统函数： 其中H(z)成为N阶IIR数字滤波器系统函数。IIR滤波器设计方法有间接和直接法，间接法是借助于模拟滤波器的设计方法进行的。其步骤是：先设计过度模拟滤波器得到系统函数Ha（s），然后将Ha（s）按某种方法转换成数字滤波器的系统函数H(z)。 利用有限脉冲响应（FIR）滤波器设计滤波器。有限脉冲响应滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。用N表示FIR滤波器单位脉冲响应h（n）的长度，其系统函数H（z）为 H（z）是z-1的N-1次多项式，它在z平面上有N-1个零点，在原点z=0处有一个N-1重极点。因此，H（z）永远稳点。稳定和线性相位是FIR滤波器最突出的优点。 （6）各种不同类型滤波器的性能比较： 巴特沃斯滤波器具有单调下降的幅频特性；切比罗夫滤波器的幅频特性在通带或阻带有等波纹特性，可以提高选择性；贝塞尔滤波器通带内有有较好的线性相位特性；椭圆滤波器的选择性相对前三种是最好的，但通带和阻带内均呈现等波纹幅频特性，相对特性的非线性稍重。 IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR滤波器的低的多。IIR数字滤波器的设计方法是利用模拟滤波器成熟的理论及设计图进行设计的，因而保留了一些典型模拟滤波器的优良的幅度特性。 （7）离散傅立叶变换 其中WN= ，N为DFT变换空间长度。 3.2 信号采集 从网上下载一段wav格式的文件，把文件“000.wav”保存在MATLAB文件夹下的work文件夹中，以.wav格式保存，这是windows操作系统规定的声音文件保存的标准。 [x1,fs]=audioread(000.wav);%把语音信号进行加载入MATLAB仿真软件平台中，采样值放在向量x1中，fs表示采样频率（Hz） x=x1(1:5000,1);%对双声道信号取单声道并取其5000点 X=fft(x,4096);%对信号做4096点FFT变换 调用参数x为被变换的时域序列向量，变换区间长度为4096，当x小于4096时，fft函数自动在x后面补零。函数返回x的4096点DFT变换结果X。当x大于4096时，fft函数计算x前面4096个元素构成的长序列的4096点DFT，忽略x后面的元素。 进行图形分区，首先画出语音信号的时域波形，然后对其进行频谱分析。在MATLAB中利用fft对信号进行快速傅立叶变换，得到信号的频谱特性。 magX=abs(X);%把傅里叶变换后的复数值取模 subplot(2,1,1);%图形分区 plot(x);title(原始信号波形);%绘制波形 f=(0:2047)*fs/2/2048;%单位转换 subplot(2,1,2); plot(f,magX(1:2048));title(原始信号频谱); 其程序如下： [x1,fs]=audioread(000.wav); %读取语音信号 x=x1(1:5000,1);%对双声道信号取单声道并取其5000点 X=fft(x,4096);%对信号做4096点FFT变换 magX=abs(X);%把傅里叶变换后的复数值取模 subplot(2,1,1);%图形分区 plot(x);title(原始信号波形);%绘制波形 f=(0:2047)*fs