《基于MATLAB的语音信号处理课程设计(6页).docx》由会员分享，可在线阅读，更多相关《基于MATLAB的语音信号处理课程设计(6页).docx（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-基于MATLAB的语音信号处理课程设计-第 6 页课程设计报告课程设计题目： 基于matlab的语音信号处理 学 号： 学生姓名: 专 业: 通信工程 班 级： 指导教师： 2016年 12 月 28 日目录一 摘要.1二 设计内容 .1三 设计步骤 .2四 实验结果 .4五实验结果分析 .5六. 实验心得 .6七. 参考文献.7八.实验代码 .7一. 摘要MATLAB软件在通信领域中占据着至关重要的作用，通信原理这本书与其更是息息相关，基于MATLAB的语音信号处理，无疑是对通信原理和MATLAB结合使用的能力的综合训练，也是对所学通信知识的有力检验。在本次设计中，采用了windows下的
2、录音机进行录音，并且通过goldwire软件进行了语音文件格式转换，便于文件在MATLAB中顺利运行，然后利用matlab软件对语音信号按照设计要求进行相应的频谱分析及滤波，最后通过回放比较语音信号的差异。二设计内容录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；根据给定的低通滤波器，对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号。三设计步骤 （1）语音信号的采集 利用Windows下的录音机，录制一段自己的话音，时间在1 s内。然后在Matlab软件平台下，利用函数wavread对语音信
3、号进行采样，记住采样频率和采样点数。通过wavread函数的使用，使学生理解采样频率、采样位数等概念。 （2）语音信号的频谱分析 画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性，加深学生对频谱特性的理解。 （3）给定某IIR数字低通滤波器如下： 即滤波器的分子系数b=0.33，0.33，分母系数 a=1,-0.33。 （4）用滤波器对信号进行滤波 用给定的低通滤波器对采集的信号进行滤波，在Matlab中， IIR滤波器利用函数filter对信号进行滤波。 （5）比较滤波前后语音信号的波形及频谱 在一个窗口同时画出滤波前后的波形及频谱。 （6）回放语音信号 在Matl
4、ab中，函数sound可以对声音进行回放。其调用格式：sound(x，fs，bits)；感觉滤波前后的声音变化。四 实验结果 1.语音信号的时域波形 2.语音信号的频谱特性 三滤波器的衰减函数波形四语音信号滤波前后的对比图形五实验结果分析 根据滤波器的衰减函数波形可知，首尾均为一条直线，是由于此滤波器由于系统函数的分子分母系数的原因，使得此低通滤波器的滤波作用不强，几乎无滤波，而刚好与语音信号滤波前后的波形基本没有发生变化的现象前后呼应，做出了合理的解释，在这里也恰恰说明了此滤波器有改进的空间。六实验心得 通过本次课程设计，不仅加深了自己对数字信号处理滤波器设计章节理解，而且对matlab的运
5、用有了新的进步，更增加对通信原理课程的认识，虽然课程设计时间较长，过程也比较吃力，因为这是多门学科的综合运用。不过幸好，从设计过程来看，复习了一遍数字信号处理课程关于数字滤波器的设计的内容，再利用提供的参考文献，加之老师的耐心指导，完成了用Matlab进行语音信号处理课程设计。终归是有所收获。这里涉及matlab的编程，由于自己有一定的matlab软件使用基础，故而使用起来并不是那么费力，但由于滤波器的设计并非完美，所以还可以进一步进行优化设计。相信随着科技的发展，MATLAB在通信领域中必将发挥更大的作用。最后，真诚感谢吴老师的耐心指导，以及设计课程中长时间的寸步不离，非常有责任心。七应收集
6、的资料及主要参考文献1， 数字滤波器设计方法；2， 数字信号处理及MATLAB实现余成波杨如民等编著 清华大学出版社出版3， 数字信号处理教程程佩青清华大学出版社出版 4， 通信原理 第六版 樊昌信 曹丽娜 编著 国防工业出版社八，实验代码clear,clcy,fs,bits=wavread(matlab.mass.wav);sound(y,fs,bits);figure(1); plot(y);title(语音信号的时域波形);Y=fft(y,4096); figure(2); plot(abs(Y);title(语音信号的频谱特性);b=0.33,0.33;a=1,-0.33;figure
7、(3);fk=0:14000/512:14000;wk=2*pi*fk;hk=freqs(b,a,wk);plot(fk/1000,20*log10(abs(hk);xlabel(频率（khz）);ylabel(幅度DB);axis(0,14,-40,5);grid on;title(滤波器的衰减函数波形);x=filter(b,a,y);X=fft(x,4096);figure(4);subplot(2,2,1);plot(y);title(滤波前信号波形);subplot(2,2,2);plot(abs(Y);title(滤波前信号频谱);Subplot(2, 2 ,3);plot(x);title(滤波后信号波形);Subplot(2, 2 ,4);plot(abs(X);title(滤波后信号频谱);sound(x,fs,bits);