《利用matlab处理语音信号及滤波》由会员分享，可在线阅读，更多相关《利用matlab处理语音信号及滤波（17页珍藏版）》请在金锄头文库上搜索。
1、word数字信号处理课程设计摘要：语音信号处理是研究用数字信号处理技术和语音学知识对语音信号进展处理的新兴学科，是目前开展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。Matlab语言是一种数据分析和处理功能十分强大的计算机应用软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！本设计要求自己录制一段自己的语音后，在MATLAB软件中采集语音信号、回放语音信号并画出语音信号的时域波形和频谱图。再在Matlab中分别设计不同形式的FIR数字滤波器。之后对采
2、集的语音信号经过不同的滤波器低通、高通、带通后，观察不同的波形，并进展时域和频谱的分析。比照处理前后的时域图和频谱图，分析各种滤波器对于语音信号的影响。最后分别收听进展滤波后的语音信号效果，做到了解在怎么样的情况下该用怎么样的滤波器。目录4445563.1.1原语音信号时域、频域图667891015161617171.设计容：1首先录制好一段自己的语音。2用Matlab分别设计好3种类型的滤波器指标自己确定：低通型、高通型、带通型。3用Matlab将语音信号进展采样，并分别将其通过所设计的3种滤波器。4用Matlab自带的语音返回函数收听滤波后的语音信号，分析并比拟其与原语音信号的差异。2.设
3、计原理：语音信号时一种非平稳的时变信号，它带着各种信息。在语音编码、语音合成、语音识别和语音增强等语音处理中无一例外需要提取语音中包含的各种信息。语音信号分析的目的就在于方便有效的提取并表示语音信号所携带的信息。语音信号处理可以分为时域和变换域等处理方法，其中时域分析是最简单的方法，直接对语音信号的时域波形进展分析，崎岖的特征参数主要有语音的短时能量，短时平均过零率，短时自相关函数等。信号提取：通过图形用户界面上的菜单功能按键采集电脑上的一段音频信号，完成音频信号的频率，幅度等信息的提取，并得到该语音信号的波形图。信号调整：在设计的用户图形界面下对输入的音频信号进展各种变化，如变化幅度、改变频
4、率等操作，以实现对语音信号的调整。信号的傅里叶表示在信号的分析和处理中起着重要的作用。因为对于线性系统来说，可以很方便地确定其对正弦或复指数和的响应，所以傅里叶分析方法能完善地解决许多信号分析和处理问题。另外，傅里叶表示使信号的某些特性变得更明显，因此，它能更深入地说明信号的各项红物理现象。由于语音信号时随着时间变化的，通常认为，语音是一个受准周期脉冲或随机噪声源激励的线性系统的输出。输出频谱是声道系统频率响应与激励源频谱的乘积。身份到系统的频率响应与激励源都是随时间变化的，因此一般标准的傅里叶表示虽然适用于周期与平稳随机信号的表示，但不能直接用于语音信号。由于语音信号可以认为在短时间，近似不
5、变，因而可以采用短时分析法。1.信号变换：在用户图形界面西啊对采集的语音信号进展Fourier等变换，并画出变换前后的频谱图和倒谱图。2.信号滤波：滤除语音信号中的噪音局部，可以采用抵用滤波、高通滤波、带通滤波，并比拟各种滤波后的效果。3设计过程：3.1实验程序源代码原语音信号存放在e:下：3.1.1原语音信号的时域、频域图x1,fs,bits=wavread(e:txwz.wav);%sound(x1,fs,bits);figure(1);plot(x1);%做原始语音信号的时域图形title(原始语音信号);xlabel(时间 t);ylabel(音量 n);figure(2);y1=ff
6、t(x1);%做length(x1)点的FFTy1=fftshift(y1); %平移，是频率中心为0derta_fs = fs/length(x1); %设置频谱的间隔，分辨率plot(-fs/2:derta_fs: fs/2-derta_fs,abs(y1);%画出原始语音信号的频谱图title(原始语音信号的频谱);grid on;fc1=4000;N1=2*pi*0.9/(0.1*pi)wc1=2*pi*fc1/fs;if rem(N1,2)=0 N1=N1+1;endWindow= boxcar (N1+1); %长度为N1的矩形窗Windowb1=fir1(N1,wc1/pi,Wi
7、ndow);figure(3);freqz(b1,1,512);title(低通滤波器的频率响应);x1_low = filter(b1,1,x1);%对信号进展低通滤波figure(4);plot(x1_low);title(信号经过低通滤波器(时域);figure(5);plot(-fs/2:derta_fs:fs/2-derta_fs,abs(fftshift(fft(x1_low);title(信号经过低通滤波器频域);fc2=4000;N2=2*pi*3.1/(0.1*pi)wc2=2*pi*fc1/fs;N2=N2+mod(N2,2);Window=hanning(N2+1);b2
8、=fir1(N2,wc2/pi,high,Window);figure(6);freqz(b2,1,512);%数字滤波器频率响应title(高通滤波器的频率响应);x1_high = filter(b2,1,x1);%对信号进展高通滤波figure(7);plot(x1_high);title(信号经过高通滤波器(时域);figure(8);plot(-fs/2:derta_fs:fs/2-derta_fs,abs(fftshift(fft(x1_high);title(信号经过高通滤波器频域)f1=4000;f2=8000; %带通滤波器的通带围w1=2*pi*f1/fs;w3=w1+0.
9、1*pi;w2=2*pi*f2/fs;w4=w2-0.1*pi;w=(w1+w3)/2,(w2+w4)/2;B=0.1*pi;N3=ceil(2*0.9*pi/B);N3=N3+mod(N3,2);Window=boxcar(N3+1);b3=fir1(N3,w/pi,stop,Window);%带通滤波器figure(9);freqz(b3,1,512);%数字滤波器频率响应title(带通滤波器的频率响应);x1_daitong = filter(b3,1,x1);%对信号进展带通滤波figure(10);plot(x1_daitong);title(信号经过带通滤波器(时域);figur
10、e(11);plot(-fs/2:derta_fs:fs/2-derta_fs,abs(fftshift(fft(x1_daitong);title(信号经过带通滤波器频域);3.1.5语音信号的回放分别执行sound(x1,fs,bits);%原始信号sound(x1_low,fs,bits);%经过低通滤sound(x1_high,fs,bits); %经过高通滤波sound(x1_daitong,fs,bits);%经过带通滤波原始语音信号的时域图形：原始语音信号频谱：低通滤波器的频率响应：信号经过低通滤波后的时域波形：信号经过低通滤波后的频域波形高通滤波器的频率响应：信号经过高通滤波后
11、的时域波形：信号经过高通滤波后的频域波形：带通滤波器的频率响应：信号经过带通滤波后的时域波形：信号经过带通滤波后的频域波形：1在高通与带通滤波器的设计时老是报错，但同样的用法在低通滤波器中就可以实现b2=fir1(N2,wc2/pi,high,Window);? Error using = fir1The window length must be the same as the filter length.其要求在fir函数中所选用的窗长要和滤波器长度一致。但在参考书上指出，滤波器阶数必须为窗长加1。经上网查询后，原来高通、带阻滤波器的阶数应该控制为奇数，因为如果阶数为偶数，如此在点必有一零
12、点，这对于高通带阻来说是不允许的，故取阶数为奇数，而你FIR1滤波器阶数为M+1阶，所以你的M必须为偶数，所以可以将程序中去窗长算法由原程序的: N2=2*pi*0.9/(0.1*pi);if rem(N2,2)=0N2=N2+1;End和: N3=2*pi*0.9/(0.1*pi);if rem(N3,2)=0N3=N3+1;End改为了：N2=N2+mod(N2,2);和: B=0.1*pi;B为过渡带宽N3=ceil(2*0.9*pi/B);N3=N3+mod(N3,2);2.在设计高通滤波器时先是使用的矩形窗，用矩形窗验证出来的结果中低频语音分量依旧很强，不能将其全部抑制在0，之后换窗
13、，选着了最小衰减可以到达53dB的海明窗，再次试验，非常成功！经过回放三个不同类型滤波器输出的语音信号，并与原语音信号比照得到了如下结论。语音高频成分音质非常锋利，齿音中，声音有些暗淡。语音低频成分音质沉稳，空间感觉强，语音浑厚。语音中频成分音质有力度，有通透感。4体会与收获以往都是通过课本来感性的认知语音信号，通过本次的课程设计，让我对语音信号有了一个较为实际的认识。于此同时，让我再次把数字信号处理与数字滤波器的设计方法重新进展了复习和学习。而对于Matlab，也再次让我感受到了其功能的强大。最为重要的是，本次课程设计让我重新审视了学习的过程：只去做实验是不行的，首先还是要思考，遇到了问题查书籍，百度搜索也只是一种手段，更加重要的是想，再理解，只有这样才能真正的做好实验。 /
