(完整word版)《MATLAB课程设计》报告DOC

1、MATLAB 课程设计报告设计题目：基于 MATLABMATLAB 勺语音信号米集与处理学生姓名：学生学号：2010054620100546专业班级：光信息科学与技术答辩时间：20132013 年 1212 月指导教师：冯明库广东技术师范学院电子与信息学院、设计目的及意义1. MATLAB软件功能简介MATLAB勺名称源自Matrix Laboratory,1984年由美国Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB!高性能的数值计算和可视化集成在一起，并提供了大量勺内置函数，从而被广泛勺应用于 科学计算、控制系统和信息处理等领域的分析、仿真和设

2、计工作。MATLAB软件包括五大通用功能，数值计算功能(Nemeric)、符号运算功能(Symbolic) 、数据可视化功能(Graphic) 、数字图形文字统一处理功能 (Notebook) 和建模仿真可视化功能 (Simulink)。其中，符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLA命令窗口。该软件有三大特点，一是功能强大； 二是界面友善、语言自然；三是开放性强。目前，Mathworks公司已推出30多个 应用工具箱。MATLA在线性代数、矩阵分析、数值及优化、数值统计和随机信号 分析、 电路与系统、 系统动力学、 次那好和图像处理、 控制理论分析和系统设计、 过

3、程控制、建模和仿真、通信系统以及财政金融等众多领域的理论研究和工程设 计中得到了广泛应用。MATLABE信号与系统中的应用主要包括符号运算和数值计算仿真分析。由于信号与系统课程的许多内容都是基于公式演算， 而MATLA借助符号数学工具箱提 供的符号运算功能，能基本满足信号与系统课程的需求。例如解微分方程、傅里 叶正反变换、拉普拉斯正反变换和z正反变换等。MATLA在信号与系统中的另一 主要应用是数值计算与仿真分析，主要包括函数波形绘制、函数运算、冲击响应 与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析和零极点图绘制等内容。 数值计算仿真分析可以帮助学生更深入地理解理论知识，

4、 并 为将来使用MATLA进行信号处理领域的各种分析和实际应用打下基础。2.本题目的意义本次课程设计的课题为基于MATLA的语音信号采集与处理，学会运用MATLA的信号处理功能，采集语音信号，并对语音信号进行滤波及变换处 理，观察其时域和频域特性，加深对信号处理理论的理解，并为今后熟练使 用MATLA进行系统的分析仿真和设计奠定基础。此次实习课程主要是为了进一步熟悉对matlab软件的使用，以及学会利用matlab对声音信号这种实际问题进行处理， 将理论应用于实际，加深对它的理解。二、设计任务及指标利用MATLA对语音信号进行分析和处理，采集语音信号后，利用MATLAB软件平台进行频谱分析；并

5、对所采集的语音信号加入干扰噪声，对加入噪声 的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。语音信号的 “ 短时谱”对于非平稳信号,它是非周期的,频谱随时间连续变化,因此由傅里 叶变换得到的频谱无法获知其在各个时刻的频谱特性。如果利用加窗的方法从语 音流中取出其中一个短断,再进行傅里叶变换,就可以得到该语音的短时谱。三、设计过程(1).语音信号的采集及分析基于声卡进行数字信号的采集。 将话筒插入计算机的语音输入插口上,启动录 音机。按下录音按钮，对话筒说话,说完后停止录音。要保存文件时,利用了计算 机上的A/D转换器,把模拟的声音信号变成了离散的量化了的数字信号,放音时,它又通过D/

6、A转换器,把保存的数字数据恢复为原来的模拟的声音信号。在Matlab软件平台下可以利用函数wavread对语音信号进行采样,得到了声音数据 变量x1,同时把x1的采样频率fs=22050Hz和数据位Nbits=16Bit放进了MATALB的工作空间。 图figure 1为原始语音信号的时域图形。 从图中可以看出在时域环 境下,信号呈现出6不规则的信号峰值。通过freqz函数绘制原始语音信号的频率 响应图figure 2。然后对语音信号进行频谱分析,在Matlab中可以利用函数fft对信号行快速傅里叶变换,得到信号的频谱图figure 3，从图中可以看出对各个 频点上的随机信号在频域进行抽样,抽

7、样频率为22050Hz。(2) .给原始信号加上一个高频噪声在Matlab中人为设计一个固定频率5500Hz的噪声干扰信号。噪声信号通常 为随机序列,在本设计中用 正弦序列代替,干扰信号构建命令函数为d=Au*sin(2*pi*5500*t),给出的干扰信号为一个正弦信号， 针对上面的语音信 号,采集了其中一段。再对噪音信号进行频谱变换得到其频谱图,从图中可以看出 干扰信号，在4000Hz和6000Hz频点处有一高峰,其中5500Hz正是本设计所要 利用的。(3) .设计一个滤波器，滤除高频噪声由模拟滤波器变换为数字滤波器时,采用的是双线性变换法,它保留的是从 模拟到数字域的系统函数表示。 在

8、滤波器的选取时,由于设计方法的侧重点不同,作出比较是困难的。如果FIR滤波器情况下，最优的设计将是椭圆滤波器。用双线 性变换法设计低通滤波器。的本设计是用双线性变换法设计BW带阻型滤波器。在MATLA中,可以利用函数butterworth,设计FIR滤波器,利用MATLAB中的函数freqz画出各滤波器的频率响应。用设计好的带阻滤波器对含噪语音信号进行滤 波,在Matlab中FIR滤波器利用函数fftfilt对信号进行滤波,FIR滤波器利用函 数filter对信号进行滤波。 在一个窗口同时画出滤波前后的波形及频谱。 从图中 可以看出,5500Hz看到的高峰消失了,语音信号与开始的一样,滤波器成

9、功的 滤除了干扰信号。利用MATLAB寸语音信号进行分析和处理，采集语音信号后， 利用MATLAB软件平台进行频谱分析；并对所采集的语音信号加入干扰噪声， 寸加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。四、结论及分析1用MATLA对原始语音信号进行分析，画出它的时域波形和频谱程序：件.wav);fs=22050;语音信号采样频率为22050 x1=wavread(C:Documentsandsound(x1,22050);%播放语音信号y1=fft(x1,1024);%对信号做1024点FFT变换f=fs*(0:511)/1024;%将0到511，步长为1的序列的值与除以

10、1024的值，赋值给ffigure(1)%创建图形窗1plot(x1)%做原始语音信号的时域图形title(原始语音信号);xlabel(time n);%x轴的名字是“time n”ylabel(fuzhi n);%y轴的名字是“fuzhi n”%读取语音信号的数据，赋给变量x1，fs相乘并figure(2)freqz(xl)%title(频率响应图)figure)subplot(2,1,1);%plot(abs(y1(1:512)%title(原始语音信号FFT频谱)subplot(2,1,2);plot(f,abs(y1(1:512); %abs线title(原始语音信号频谱)xlabe

11、l(Hz);ylabel(幅值);绘制原始语音信号的频率响应图创建两行一列绘图区间的第1个绘图区间做原始语音信号的FFT频谱图是绝对值，plot是直角坐标下线性刻度曲时域波形和频谱:图1原始语音信号Figure 2SB图2语音信号频率响应图图3原始语音信号FFT与信号频谱2.给原始的语音信号加上一个高频余弦噪声，频率为5500hz。对加噪后的语音进行分析，并画出其信号时域和频谱图。程序：fs=22050;xlabel(Hz);%x轴的名字是“time n”ylabel(fuzhi);% y轴的名字是“fuzhi”subplot(2,1,2);%创建两行一列绘图区间的第2个绘图区间plot(f,

12、abs(y2(1:512);%做加噪后的语音信号的频谱图（abs是绝对值，plot是直角坐标下线性刻度曲线）title（加噪后的信号频谱）;xlabel(Hz);%x轴的名字是“time n”ylabel(fuzhi);% y轴的名字是“fuzhi”%读取语音信号的数据，赋给变量x1f=fs*(0:511)/1024;%将0到511，步长为1的序列的值与fs相乘并除以1024的值，赋值给ft=0:1/fs:(length(x1)-1)/fs; %将0到x1的长度减1后的值除以fs的值，且步长为1/fs的值，的序列的值，赋予tAu=0.03;%噪声幅值d=Au*sin(2*pi*5500*t);

13、 %所加的噪声是正弦信号x2=x1+d;%将正弦信号噪声加在语音信号上sound(x2,22050);%播放语音信号y1=fft(x1,1024);%对信号y1做1024点FFT变换y2=fft(x2,1024);%对信号y2做1024点FFT变换figure(1);%创建图形窗1plot(t,x2);%做加噪后的信号时域图形 （绘出t对x2的线性图）title（加噪后的信号);xlabel(time n);%x轴的名字是“time n”ylabel(fuzhi n);%y轴的名字是“fuzhi n”figure(2)%创建图形窗2subplot(2,1,1);%创建两行一列绘图区间的第1个绘

14、图区间plot(f,abs(y1(1:512);%做原始语音信号的频谱图桌 面 声 音 文件.wav);title（原始语音信号频谱）;x1=wavread(C:Documents and时域波形和频谱:图4加噪后的信号时域图图5原始信号与加噪后信号频谱对比与原始信号对比，区别：先原始信号没加噪音之前0到2000有幅值，在4000到6000之间没有幅值，但是在加了噪音之后4000到6000之间出现最大幅值12,超出正常值。如图5原始信号与加噪后信号频谱对比图所示3.设计合适的滤波器，滤除高频噪声，绘出滤波后的信号频域和时域波形： 程序：%用双线性变换法设计巴特沃思低通滤波器wp=0.25*pi

15、;%通带截止频率ws=0.3*pi;%阻带截止频率Rp=1;%通带最大衰减（db）Rs=15;%阻带最大衰减（db）Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2);%将模拟指标转换成数字指标N,Wn=buttord(wp1,ws1,Rp,Rs,s); %选择滤波器的最小阶数（估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率WCZ,P,K=buttap(N);%创建Butterworth低通滤波器原型Bap,Aap=zp2tf(Z,P,K);%将零极点增益转换为普遍分子，分母b,a=lp2

16、lp(Bap,Aap,Wn);%将普遍的分子和分母转换为以Wn为截止频率bz,az=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换（模拟转换为数字）H,W=freqz(bz,az);%求频率响应figure(9)plot(W*Fs/(2*pi),abs(H)%绘制Butterworth低通滤波器频率响应曲线grid %添加图格xlabel(频率Hz)ylabel（频率响应幅度）title(Butterworth)f1=filter(bz,az,z);%利用上面已做好的滤波器（已bz和az为特征的滤波器）对x2信号进行滤波，赋值给f1figure(10)sub

17、plot(2,1,1) %间plot(t,z)%title(滤波前的时域波形);subplot(2,1,2) %间plot(t,f1);title(滤波后的时域波形);sound(f1,22050); % p=length(f1);F0=fft(f1,p); f=0:fs/p:fs*(s-1)/p;figure(11)y2=fft(z,1024);%subplot(2,1,1);%间plot(f,abs(y4);title(滤波前的频谱) xlabel(Hz);%xylabel(fuzhi);%ysubplot(2,1,2) F1=plot(f,abs(F0);title(滤波后的频谱) xl

18、abel(Hz);ylabel(fuzhi);创建两行一列绘图区间的第1绘图区做加噪后的信号时域图形创建两行一列绘图区间的第2绘图区播放语音信号对信号y2做1024点FFT变换创建两行一列绘图区间的第1绘图区轴的名字是“Hz” 轴的名字是“fuzhi图7滤波前后时域波形图时域波形和频谱:7 Figure 9-r xEile Edit丑“Ins er t L&ols Itsktop jjtindoii Help O B| fe|O|DS|s0Butterworth1JIII b I 1f1 1 I I 11111riid1F4i1111, _ _ _ _ _ x_ _k41R11hdi1

19、1btl1F11btl11111i11_ _ _ _ _ _X. . . =二二二111111i.1 1.7i i i i i i i丨丨丨- - - - - f -111i111厂门1. - - _-T- -111i 1 i:1|1丨U丨11111110.2200D1D00Q120006D00 BOOD频率/Hzo.ao.s图6低通滤波器频率响应曲线*/ FigiiEE 1Dh xEile Edit野“Ins er t L&ols Itsktop jStindon HlpDO a fe gQ紳倉盪Id目B Q0.40.20-D2Q.2Q.3Q.50.S0.7as ag1图8滤波前后频

20、谱图与原噪音信号对比，区别：通过对比分析可知，滤波后的输出波形和原始语音加噪声信号的图形发生了一些变化，在 3000-6000HZ3000-6000HZ 之间的信号消失，出现 0-10000-1000 和7000-80007000-8000 之间的信号。滤波后的输出波形明显在对应时间幅度比 原语音加噪声信号的要小，而且滤波的效果也与滤波器的选择有关， 可以看出滤波器的性能差异以及参数的调节会对滤波器产生一定的影响，并且通过回放可以发现滤波前后的声音有变化低通滤波后，已很接近原来的声音，人耳几乎辨别不出。从频谱图中我们还可以看出声 音的能量信号主要集中在低频部分，说明高频语音信号被滤出，滤波 器

21、达到既定要求。五、设计体会正所谓“纸上得来终觉浅，觉知此事要躬行。”学习任何知识，仅从理论上去 求知，而不去实践、探索是不够的。因此在学期末来临之际，我们迎来了MATLAB课程设计。通过为期一周半的MATLA課程设计，我对MATLA这个仿真软件有了更进一步的认识和了解。在这一周半时间里，我通过自己摸索，查阅资料，并且在老师 的指导下完成了：语音信号的采集及分析；给原始信号加上一个高频噪声；设计 一个滤波器，滤除高频噪声；并最终将课程设计报告总结完毕。在整个设计过程中我懂得了许多东西，也培养了独立思考和设计的能力，树 立了对知识应用的信心， 相信会对今后的学习工作和生活有非常大的帮助， 并且 提

22、高了自己的动手实践操作能力，使自己充分体会到了在设计过程中的成功喜 悦。虽然这个设计做的不怎么好，但是在设计过程中所学到的东西是这次课程设 计的最大收获和财富，使我终身受益。在没有做课程设计以前，觉得课程设计只是对知识的单纯总结，但是通过这 次课程设计发现自己的看法有点太片面，课程设计不仅是对前面所学知识的一种 检验，也是对自己能力的一种提高，通过这次课程设计使自己明白了原来的那点 知识是非常欠缺的，要学习的东西还很多，通过这次课程设计，明白学习是一个 长期积累的过程，在以后的工作和生活中都应该不断的学习，努力提高自己的知 识和综合素质。希望以后像这样的课程设计在多一点。参考文献：1姚东.MA

23、TLAB及在电子信息课程中的应用M.第二版.北京：北京电子工业出 版社,20002张智星.MATLAE程序设计与应用北京：清华大学出版社，20023陈怀琛.MATLAB及其在理工课程中的应用指南.西安：西安电子科技大学出版社，1999.104肖伟、刘忠. MATLAB程序设计与应用M.北京：清华大学出版社2005x1=wavread(C:Documents and件.wav）;%读取语音信号的数据，赋给变量x1，播放语音信号对信号做1024点FFT变换将0到511，步长为1的序列的值与fs相乘并 除以1024的值，赋值给f创建图形窗1做原始语音信号的时域图形轴的名字是“time n” 轴的名字

24、是“fuzhi n”绘制原始语音信号的频率响应图创建两行一列绘图区间的第1个绘图区间做原始语音信号的FFT频谱图title（原始语音信号FFT频谱）subplot(2,1,2);plot(f,abs(y1(1:512); %abs title（原始语音信号频谱）xlabel(Hz); ylabel(幅值);fs=22050;x1=wavread(C:Documents and件.wav）;读取语音信号的数据，赋给变量x1程序清单：fs=22050; %语音信号采样频率为22050sound(x1,22050);%y1=fft(x1,1024);%f=fs*(0:511)/1024;%figur

25、e(1)%plot(x1)%title（原始语音信号);xlabel(time n);%xylabel(fuzhi n);%yfigure(2)freqz(x1)%title（频率响应图）figure(3)subplot(2,1,1);%plot(abs(y1(1:512)%是绝对值，plot是直角坐标下线性刻度曲f=fs*(0:511)/1024;%将0到511，步长为1的序列的值与fs相乘 并除以1024的值，赋值给ft=0:1/fs:(length(x1)-1)/fs; %将0到x1的长度减1后的值除以fs的值，且步长为1/fs的值，的序列的值，赋予tAu=0.03;%噪声幅值d=Au*

26、sin(2*pi*5500*t); %所加的噪声是正弦信号x2=x1+d;%将正弦信号噪声加在语音信号上sound(x2,22050);%播放语音信号y1=fft(x1,1024);%对信号y1做1024点FFT变换y2=fft(x2,1024);%对信号y2做1024点FFT变换figure(1);%创建图形窗1plot(t,x2);%做加噪后的信号时域图形（绘出t对x2的线性图）title（加噪后的信号);xlabel(time n);%x轴的名字是“time n”ylabel(fuzhi n);%y轴的名字是“fuzhi n”figure(2)%创建图形窗2subplot(2,1,1);

27、%创建两行一列绘图区间的第1个绘图区间plot(f,abs(y1(1:512);%做原始语音信号的频谱图title（原始语音信号频谱）;xlabel(Hz);%x轴的名字是“time n”ylabel(fuzhi);% y轴的名字是“fuzhi”subplot(2,1,2);%创建两行一列绘图区间的第2个绘图区间plot(f,abs(y2(1:512);%做加噪后的语音信号的频谱图（abs是绝对值，plot是直角坐标下线性刻度曲线）title（加噪后的信号频谱）;xlabel(Hz);%x轴的名字是“time n”ylabel(fuzhi);% y轴的名字是“fuzhi”%用双线性变换法设计巴特沃思低通滤波器wp=0.25*pi;%通带截止频率ws=0.3*pi;%阻带截止频率title（滤波前的时域波形）;Rp=1;Rs=15;Fs=22050;%通带最大衰减（db）%阻带最大衰减（db）Ts=1/Fs;wp1=2/Ts*tan(wp/2);ws