语音信号的处理与滤波

广西科技大学课程设计用纸1广西科技大学数字信号处理课程设计题目：语音信号的处理与滤波(难度系数：0.8)系别：计算机工程学院专业：通信工程班级：学号：姓名：指导教师：日期：广西科技大学课程设计用纸2目录摘要：...................................................2一．设计内容.............................................3二．设计目的.............................................4三．基本步骤.............................................3四、相关原理知识.........................................5五、实现过程.............................................91．录制声音...........................................92．分别取8000个和16000个数据进行频谱分析，得到幅度和相位谱，比较二者异同并分析原因.........................113．滤波器的设计......................................134．对声音进行滤波....................................145．把处理后的所有数据存储为声音文件，与原始声音进行比较.....................................................15六、心得体会............................................16七、参考文献............................................17摘要：本次设计主要是利用MATLAB工具对自己的语音进行录制，然后设计IIR滤波器，并通过所设计的滤波器进行语音信号滤波以及分析等，广西科技大学课程设计用纸3对以前所学的知识进行一个系统的归纳和进一步理解，并且在实际应用操作当中体会我们《数字信号处理》这门课程的相关知识点。一、设计内容1.设计题目一:语音信号的处理与滤波（难度系数：0.8）1．熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数。2．在MATLAB环境中，使用声音相关函数录制2秒左右自己的声音，抽样率是8000Hz/s。（考虑如何解决一个实际问题：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象，如何排除对这些无效点的采样？）3．画出语音信号的时域波形;然后对语音信号进行FFT,得到信号的频谱特性,然后分别取8000个和16000个数据进行频谱分析，得到幅度和相位谱，比较二者异同并分析原因。4．针对电话信道（最高3500Hz），设计一个FIR或IIR滤波器进行滤波，把抽样率转变为7000Hz/s，并进行频谱分析，得到幅度和相位谱。5.用自己设计的IIR滤波器分别对采集的信号进行滤波比较滤波前后语音信号的波形及频谱..画出滤波前后的波形和频谱.加深对滤波器的了解.6．把处理后的所有数据储存为声音文件，与原始声音进行比较。7.保存滤波后的语音二、设计目的综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。广西科技大学课程设计用纸4三、基本步骤1．语音信号的采集熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数，在MATLAB环境中，使用声音相关函数录制3秒左右自己的声音。（考虑如何解决一个实际问题：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象，如何排除对这些无效点的采样？）2．语音信号的频谱分析要求首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性，从而加深对频谱特性的理解。3．设计数字滤波器并画出其频率响应给出各滤波器的性能指标：(1)低通滤波器性能指标：fb＝1000Hz，fc＝1200Hz，As＝100dB，Ap＝1dB。(2)高通滤波器性能指标：fc＝4800Hz，fb＝5000HzAs＝100dB，Ap＝1dB。(3)带通滤波器性能指标：fp1＝1200Hz，fp2＝3000Hz，fsc1＝1000Hz，fsc2＝3200Hz，As＝100dB，Ap＝1dB。要求用窗函数法和双线性变换法设计上面要求的3种滤波器。在Matlab中，使用函数butter，cheby1等函数设计IIR滤波器；利用Matlab中的函数freqz画出各滤波器的频率响应。4．用滤波器对信号进行滤波要求用自己设计的各种滤波器分别对采集的信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。5．比较滤波前后语音信号的波形及频谱要求在一个窗口同时画出滤波前后的波形及频谱。6．回放语音信号在Matlab中，使用相关声音函数对声音进行回放。辨别滤波前广西科技大学课程设计用纸5后的声音有何变化。四、相关原理知识1.录音原理Matlab中有两个函数可以用于实现录音功能，一个是waverecord函数调用windows音频设备进行录音，主要用于实时语音采集。第二个是使用Matlab中的wavread函数也可以进行语音采集。不过wavread函数只能读取计算机中已经有的“.wav”文件，也就是说必须事先使用windows的录音机等设备把有关的语音资料录制并存储为“.wav”文件。显然wavread函数不能用于实时语音信号采集。另外语音采集的频率和精度不能改变。回放时也有两个函数可以选用。Sound函数和wavplay，只不过后者的功能没有前者强大。不过在本设计中主要是对已有的信号进行处理，所以选用wavread录放。同时，由于信号的时域描述只能反映信号的幅值随时间的变化情况，除只有一个频率分量的简谐波外一般很难明确揭示信号的频率组成和各频率分量的大小，因此，我们要把它转换成频域来分析，它能够提供比时域信号波形更直观，丰富的信息。2、滤波器的设计原理及设计方法FIR滤波器的设计问题在于寻求一系统函数)(zH，使其频率响应)(jeH逼近滤波器要求的理想频率响应)(jdeH，其对应的单位脉冲响应)(nhd。1、用窗函数设计FIR滤波器的基本方法设计思想：从时域从发，设计)(nh逼近理想)(nhd。设理想滤波器)(jdeH的单位脉冲响应为)(nhd。以低通线性相位FIR数字滤波器为例。广西科技大学课程设计用纸6deeHnhenheHjnjddjnndjd)(21)()()()(nhd一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器h(n)，最直接的方法是截断)()()(nwnhnhd，即截取为有限长因果序列，并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即2/)1()()()(Nanwnhnhd用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，这个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。2、典型的窗函数（1）矩形窗(RectangleWindow))()(nRnwN其频率响应和幅度响应分别为：21)2/sin()2/sin()(NjjeNeW，)2/sin()2/sin()(NWR（2）三角形窗(BartlettWindow)121,122210,12)(NnNNnNnNnnw其频率响应为：212])2/sin()4/sin([2)(NjjeNNeW广西科技大学课程设计用纸7（3）汉宁(Hanning)窗，又称升余弦窗)()]12cos(1[21)(nRNnnwN其频率响应和幅度响应分别为：)]12()12([25.0)(5.0)()()]}12()12([25.0)(5.0{)()21(NWN（4）汉明(Hamming)窗，又称改进的升余弦窗)()]12cos(46.054.0[)(nRNnnwN其幅度响应为：)]12()12([23.0)(54.0)(NWN（5）布莱克曼(Blankman)窗，又称二阶升余弦窗)()]14cos(08.0)12cos(5.042.0[)(nRNnNnnwN其幅度响应为：)]14()14([04.0)]12()12([25.0)(42.0)(NWNWNWN（6）凯泽(Kaiser)窗10,)())]1/(21[1()(020NnINnInw其中：β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，β越大，过渡带越宽，阻带越小衰减也越大。I0(·)是第一类修正零阶贝塞尔函数。若阻带最小衰减表示为ssA10log20，β的确定可采用下述经验公式：广西科技大学课程设计用纸850)7.8(1102.05021)21(07886.0)21(5842.02104.0ssssssAAAAAA若滤波器通带和阻带波纹相等即δp=δs时，滤波器节数可通过下式确定：136.1495.7FANs式中：22psF3．利用窗函数设计FIR滤波器的具体步骤如下：（1）按允许的过渡带宽度△ω及阻带衰减AS，选择合适的窗函数，并估计节数N：其中A由窗函数的类型决定。（2）由给定的滤波器的幅频响应参数求出理想的单位脉冲响应)(nhd。（3）确定延时值（4）计算滤波器的单位取样响应)(nh，)()()(nwnhnhd。（5）验算技术指标是否满足要求。3.IIR数字滤波器IIR数字滤波器的系统是Z的有理数,可表示为系统的设计就是要确定系数,以使滤波器达到要求.设计方法为,用模拟滤波器理论来设计数字滤波器.在IIR数字滤波器的设计中较多的采用这种方法.用计算机助设计,优化技术设计.针对电话信道我们可以设计一个IIR低通滤波器来滤我们录制的语音,用fftfilt函数可以把我们录制的语音放到设计好的滤波器进行滤波.然后angle函数求出滤波后语音的相位谱,用abs求出幅度谱.然后再用wavwrite对滤波后的语音保存为声音文件.广西科技大学课程设计用纸9五、实现过程1．录制自己的声音文件（wav格式），内容为“数字信号课程设计，邓国望”程序如下：fs=8000;%抽样频率channel=1;%1信道t=3;%设定录音时间为3秒钟fprintf('按任意键后开始%d秒录音：',t);pause;%文字提示%暂停命令fprintf('录音中...');%提示当前正在录音当中x=wavrecord(t*fs,fs,channel,'double');%录制自己的声音信号fprintf('录音结束\n');%提示录音已经结束wavwrite(x,fs,'F:\新建文件夹\dgw.wav');%存储录音信号fprintf('按任意键后回放:');pause;wavplay(x,fs);%播放录制的声音下面读取原始声音，并且画出它的时域图形，程序如下：y=wavread('F:\新建文件夹\dgw.wav',[120000]);%读取声音plot(y);%画出声音时域图形tit