语音采样和处理

1、指导教师评定成绩:南昌航空大学科技学院数字信号处理课程设计报告题目语音采集和处理学生姓名： 王定雪专业：自动化班级：108202225学号：25设计时间：2013年1月目录一实践的目的和要求 31. 课程设计基本要求 32. 课程设计实现 3二.设计原理 31、采样定理 32、采样频率 43 、采样位数与采样频率 44 、利用MATLAB对语音信号进行分析和处理 45 、语音信号的滤波示意图 4三.设计步骤 51、语音信号的采集与分析 52、利用余弦函数构造一个高频噪声 53、设计合适的滤波器，滤除高频噪声 6四程序设计和仿真图形 64.1语音信号原程序频谱分析 64.2语音信号加噪与频谱分析

2、 74.3双线性变换法设计低通滤波器对噪音滤波 94.4双线性变换法设计带通滤波器对噪音滤波 11五.心得体会 13参考文献 13课题:语音信号采集与处理一、实践的目的和要求1. 课程设计基本要求 学会MATLAB勺使用，掌握MATLAB勺程序设计方法； 掌握在Windows环境下语音信号采集的方法； 掌握数字信号处理的基本概念、基本理论和基本方法； 掌握MATLABS计FIR和IIR数字滤波器的方法； 学会用MATLAB寸信号进行分析和处理。2. 课程设计实现(1) 语音信号的采集利用windows下的录音机或其他软件，录制一段自己的话音，时间 控制在1秒左右。然后在MATLAB件平台下，禾

3、I用函数 wavread对语音信号进 行采样，记住采样频率和采样点数。通过wavread函数的使用，要求理解采样频 率、采样位数等概念。wavread函数调用格式：y=wavr，读取file 所规定的wav文件，返回采样值放在向量y中。y,fs,nbits=wavread(file)，采样值放在向量 y中，fs表示采样频率(Hz)，nbits表示采样位数。y=wavread(file,N)，读取前N点的采样值放在向量 y中。y=wavread(file,N1,N2) ，读取从N1点到N2点的采样值放在向量 y中。1、设计原理：1、采样定理:在进行模拟与数字信号的转换过程中，当采样大于最高频率的

4、2倍 时，则采样之后的数字信号完整的保留了原始信号中的信息，一般实际 应用中保证采样频率为信号最高频率的 510倍。2、采样频率：采样频率是指计算机每秒钟采样多少个声音样本，是描述声音文件 的音质、音频、衡量声卡、声音文件的质量标准。采样频率越高，即采 样的时间间隔越短，则在单位时间内计算机得到的声音样本数据越多， 对声音波形的表示也越准确。3 、采样位数与采样频率：采样位数即采样值或取样值，用来衡量声音波动变化的参数，是指 声卡在采集和播放声音文件时所使用数字声音信号的二进制位数。采样频率是指录音设备在一秒钟内对声音信号的采样次数，采样频率越高声 音的还原就越真实越自然。采样位数和采样率对于

5、音频接口来说是最为 重要的两个指标。无论采样频率如何，理论上来说采样的位数决定了音 频数据最大的力度范围。采样位数越多则捕捉到的信号越精确。4、利用MATLA对语音信号进行分析和处理：采集语音信号后，利用 MATLA软件平台进行频谱分析；并对所 采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析， 设计合适的滤波器滤除噪声，恢复原信号。5 、语音信号的滤波示意图：二.设计步骤1、语音信号的采集与分析(1) 利用PC机上的声卡和 WINDOWS作系统可以进行数字信号的采集。将话 筒插入计算机的语音输入插口上，启动录音机。按下录音按钮，对话筒说话“ 我 想回家”，说完后停止录音，屏幕左侧将显

6、示所录声音的长度。(2) 以文件名“2”保存入C盘中。可以看到，文件存储器的后缀默认为.wav。 要保存文件时，利用了计算机上的A/D转换器，把模拟的声音信号变成了离散的 量化了的数字信号,放音时,它又通过D/A转换器,把保存的数字数据恢复为原来 的模拟的声音信号。(3) 在Matlab软件平台下可以利用函数 wavread对语音信号进行采样,得到了 声音数据变量x1,x仁wavread('C:22.wav') ;同时把x1的采样频率fs=22050Hz 和数据位Nbits=16Bit放进了 MATALB勺工作空间。图figure 1为原始语音信号 的时域图形。(4) 从图中可

7、以看出在时域环境下，信号呈现出4个不规则的信号峰值。通过freqz(x1)函数绘制原始语音信号的频率响应图 figure 2。(5) 然后对语音信号进行频谱分析，在Matlab中可以利用函数y仁fft(x1,1024);plot(abs(y1(1:512)对信号行快速傅里叶变换,得到信号的频谱图figure 3 ,从图中可以看出对各个频点上的随机信号在频域进行抽样，抽样频率为22050Hz2、利用余弦函数构造一个高频噪声1.在Matlab中人为设计一固定频率5000Hz的噪声干扰信号。噪声信号通常为 随机序列，在本设计中用余弦序列代替。干扰信号构建命令函数为 d=Au*cos(2*pi*500

8、0*t)' x仁wavread('C:22.wav'); x2=x1+d;2给出的干扰信号为一个余弦信号，针对上面的语音信号，采集了其中一段。再 对噪音信号进行频谱变换得到其频谱图 , y2=fft(x2,1024) ；从图中可以看出干 扰信号，在4000Hz和6000Hz频点处有一高峰,其中5000Hz正是本设计所要 利用的。若带噪信号y(m)是纯净语音信号s(m)和平稳加性高斯白噪声n(m)构成的, 即:y(m) =x(m) +n(m);则在频域中表示为 Y( 3 )=X( 3 )+N( 3 );其中 Y( 3 ),X( 3 ),N( 3)分别是 y(m),x(m)

9、,n(m)的傅里叶变换。X( 3 )=Y( 3 )-N( 3); 只要从带噪信号中减去噪声信号即可估计出语音信号 , 再进行傅里叶反变换就可 得到增强的语音。3、设计合适的滤波器，滤除高频噪声由模拟滤波器变换为数字滤波器时 , 采用的是双线性变换法 , 它保留的是从 模拟到数字域的系统函数表示。用双线性变换法设计低通滤波器。在MATLAB中, 可以利用函数 butterworth 设计低通滤波器。数字巴特沃兹低通滤波器函数：functionb,a=digital_lowpass_butter(wp,ws,Rp,As,Fs) N,wn=buttord(wp/pi,ws/pi,Rp,As)b,a=

10、butter(N,wn);时域低通滤波：x=wavread( 2.wav');x 为读取原语音文件的数据，一维数组h=ones(1,220); h为低通滤波器的单位冲激响应y=conv(x,h);时域卷积，低通滤波wavwrite(y,22050, su2.wav );将滤波后的数据保存为.wav文件用MATLAB中的函数freqz画出各滤波器的频率响应。用设计好的带阻滤波 器对含噪语音信号进行滤波,在Matlab中IIR滤波器利用函数filter对信号进行 滤波。在一个窗口同时画出滤波前后的波形及频谱。从图中可以看出,5000Hz看到的高峰消失了 ,语音信号与开始的一样 , 滤波器成

11、功的滤除了干扰信号。利用MATLA对语音信号进行分析和处理，采集语音信号后，利用MATLAB件平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信 号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。四、程序设计及仿真图形4.1、用 Matlab 对原始语音信号进行分析，画出它的时域波形和频谱程序:clear;y,fs,bits=wavread('C:22.wav'); %x 语音数据 fs 采样频率 bits 采样点数 sou nd(y,fs,bits); %话音回放n = len gth (y) ; %求出语音信号的长度Y=fft(y, n);%傅里叶变换s

12、ubplot(2,1,1);plot(y);title('原始信号波形');subplot(2,1,2);plot(abs(Y);title('原始信号频谱')4.2、给原始的语音信号加上一个高频余弦噪声，频率为(22050hz)。对加噪后 的语音进行分析，并画出其 信号时域和频谱图。程序：%加噪声fs=22050;x=wavread('C:22.wav');f=fs*(0:511)/1024;Au=0.03;t=0:1/22050:(le ngth(x)-1)/22050;d=Au*cos(2*pi*5000*t);x2=x+d'd&#

13、39;y2=fft(x2,1024);figure(1);plot(t,x2)title('加噪后的信号');xlabel('time n');ylabel('fuzhi n');figure(2)subplot(2,1,1);plot(f,abs(y2(1:512);%plot(f,abs(y1(1:512);title('原始语音信号频谱');xlabel('Hz');ylabel('fuzhi');subplot(2,1,2);plot(f,abs(y2(1:512);title('加

14、噪后的信号频谱');xlabel('Hz');ylabel('fuzhi');U LINmlimen匣対语音悟号煎谱-5DoJOJOfDOO100M12000151 11 1 11 丿L1X11加輛的信号躺20004000EOOO8000loan12000nu5上Nnj.4.3、设计合适的滤波器，滤除高频噪声，绘出滤波后的信号频域和时域波形: 程序：%用双线性变换法设计巴特沃思低通滤波器双线性变换法低通滤波器fp=1000;fc=1200;As=100;Ap=1;ffs=22050;wc=2*fc/ffs;wp=2*fp/ffs;n,wn=ellipor

15、d(wp,wc,Ap,As);%求阶数 N 和通带截止频率n um,de n=ellip( n,Ap,As,w n);%设计椭圆滤波器freqz( num,de n,256,ffs);%画出频率响应200040DDB00DFrequency (HzJ000010DODInS-仝 apmeBew5LO50200040006000eoro10000Frequency (HR富aJMsa壘亠对噪音进行滤波：x=filter( nu m,de n, y);%滤波X二fft(x,4096);%FFTsubplot(211),plot(x); title('滤波后信号波形');subplo

16、t(212),plot(abs(X); title('滤波后信号频谱');滤波后图形：10探波后信号频语1210364-1II1-IIJ-5001QQ0150020002600300035004QOO454.4、双线性变换法的带通滤波器和对噪音滤波fp=3000,8000;fc二1000,10000;Ap=1;As=100;ffs=22050 wc=2*fc/ffs;wp=2*fp/ffs;n ,w n二ellipord(wp,wc,Ap,As);nu m,de n=ellip( n,Ap,As,w n);freqz (nu m,de n,256,ffs);20004COO60

17、00Frequency (Hz)B000II DODO口-50-100-150-200-250-300O BO20004COO6000Frequency (Hz)BOOQ10000500O -fil进行滤波：x=filter( nu m,de n, y);%滤波滤波后信号波形');X二fft(x,4096);%FFT subplot(211),plot(x); title(' subplot(212),plot(abs(X); title('滤波后信号频谱');滤波后競披启信日获it五、 心得体会：通过这一个星期的课程设计， 我学到了很多的东西 , 不仅巩固了我

18、以前所 学过的知识 , 还让我学到很多在书本上所没有学到过的知识。同时进一步加深了对语音信号的了解和熟练了对 Matlab 的使用 , 让我 对通信原理这门课程有了更加浓厚的兴趣。 因为以前都是基于课本上所学的 理论知识，然而通过这次课程设计之后才能真正理解其意义。在这次课程设计的过程中，我遇到不少的问题 ,比如刚开始 ,要录音的时 候，没有选择正确的频率，导致在运行程序的时候，无法调用声音，经过后 来的改正才可以。还有刚开始由于对滤波器的滤波原理并不是很了解 , 于是 我又翻出学过的数字信号处理课本 , 认真研究起各种滤波器 ,这才使我明白了 大多数滤波器是如何工作 ,不再单单只是懂理论 , 理论与实际相结合是很重要 的 , 只有理论知识是远远不够的 , 只有把所学的理论知识与实践相结合起来 从理论中得出结论。 总的来说 ,通过这次的课程设计我对语音信号有了全面的 认识,对 Matlab 的知识又有了深刻的理解 , 让我感受到只有在充分理解课本 知识的前提下 , 才能更好的应用这个工具。这次