基于MATLAB的语音信号分析与处理的课程设计实验报告

1、目 录绪论1正文：一、设计目的及要求2二、设计过程2三、调试分析6四、结果分析与体会11五、参考文献12六、谢辞12绪 论数字信号处理（digital signal processing,dsp）是从20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。它的重要性日益在各个领域的应用中表现出来。简言之，数字信号处理是把信号用数字或符号表示的序列，通过计算机或通用（专用）信号处理设备，用数字的数值计算方法处理（例如滤波、变换、压缩、增强、估计、识别等），以达到提取有用信息便于应用的目的。数字信号处理的应用包括滤波与变换、通信、语音、语言、图像、图形、消费电子、仪器

2、、工业控制与自动化、医疗、军事等。其发展方向也是多方面的，包括数字汇聚、远程会议系统、融合网络、数字图书馆、图像与文本合一的信息检索业务、多媒体通信、个人信息终端等。在本次设计中，就以设计一个数字滤波器系统为目的。滤波器按频率划分有低通、高通、带通、带阻全通等类型。一个数字滤波器可以用一个系统函数来表示，若要得到该系统函数，就要得到以z(-1)升幂排列的传输函数的分子和分母多项式的系数，即向量num、den，而这两个参数由num den=butter(n,wn)、n,wn=buttord(wp,ws,rp,rs)这一函数得出，这是巴特沃思滤波器的设计过程。由于频率响应的周期性，频率变量以数字频

3、率来表示，=2*f/fc，f为模拟域频率，fc为抽样频率，所以数字滤波器设计中必须给出抽样频率。在这次设计中采用了巴特沃思和切比雪夫两种滤波器，这两种滤波器又分别采用了高通、低通、带通、带阻四种频率特性实现，从横向上和纵向上都可以进行对比，表现出各种滤波器的不同特性。 一、设计目的及要求1课程设计目的 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用 matlab 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。2课程设计基本要求 1） 学会 matlab 的使用，掌握 matlab 的程序设计方法。 2） 掌握在 windows 环境下

4、语音信号采集的方法。 3） 掌握数字信号处理的基本概念、基本理论和基本方法。 4） 掌握 matlab 设计 fir 和iir 数字滤波器的方法。 5） 学会用 matlab 对信号进行分析和处理。3课程设计内容录制一段自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，设计数字滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号；最后，用 matlab 设计一信号处理系统界面。二、设计过程1、 打开gui窗口：file new g

5、ui create new gui blank gui(default) ok2、 控件设计：在控件布局设计区放置四个axes控件、六个push button控件、一个button group控件和两个radio button控件、四个static text控件。3、 修改控件属性：选中需要修改属性的控件，双击打开属性查看器，具体设置如下：1) 五个axes控件 tag作用axes1提供坐标画出原始信号波形axes2提供坐标画出原始信号频谱axes3提供坐标画出滤波后信号波形axes4提供坐标画出滤波后信号频谱axes5提供坐标画出滤波器频率响应2) 六个push button控件：strin

6、gtagstringtag低通pushbutton1带阻pushbutton4高通pushbutton2原始信号pushbutton5带通pushbutton3关闭窗口pushbutton63) 一个button group控件和两个radio button控件： stringstyletag请选择滤波器种类button groupuipanel5巴特沃思radio buttonradiobutton1切比雪夫radio buttonradiobutton24) 五个static text控件（这五个static text控件是对坐标轴中图形的说明）：tagstringtext1原始信号波形t

7、ext2原始信号频谱text3滤波后信号波形text4滤波后信号频谱text5滤波器频率响应4、 保存：设置好各个控件的属性，回到gui主窗口保存，给文件命名为hy，同时hy.m文件打开。5、 设置回调函数：在hy.m文件窗口中设置回函数。这里，虽然gui自动生成了回调函数，但是回调函数是空的，需要在hy.m文件中对它进行定义说明。该程序只需要对六个push button控件的回调函数进行定义说明。六个push button控件的回调函数：1）“低通”按键的回调函数function pushbutton1_callback(hobject, eventdata, handles)y,fs,bi

8、ts=wavread('e:hy.wav'); %函数wavread是对语音信号进行 采样，格式是y,fs,nbit=wavread，返回采样值放在向量y中，fs表示采样频率（hz），nbit表示采样位数。fp=1000;fs=2000;rp=0.5;rs=40;fc=40000;%设定通带截止频率（fp）、阻带截止频率（fs）、通带波纹系数（rp）、阻带波纹系数（rs）、抽样频率（fc）。wp=2*fp/fc;ws=2*fs/fc; %将模拟域转化成数字域。if get(handles.radiobutton1,'value') %如果选择radiobutto

9、n1，则制做巴特沃思滤波器。n,wc=buttord(wp,ws,rp,rs); %估算巴特沃思滤波器的阶数n和3db截止频率wc。num,den=butter(n,wc); %求传输函数的分子和分母多项式的系数。else %选择radiobutton2，则制做切比雪夫型滤波器。n,wc=cheb1ord(wp,ws,rp,rs); %估算切比雪夫型滤波器的阶数n和截止频率wc。num,den=cheby1(n,rp,wc); %求传输函数的分子和分母多项式的系数。endx=filter(num,den,y); %根据传输函数的分子和分母多项式的系数得到模拟滤波器，并将y通过该滤波器，得到x。

10、sound(x,fc,bits); %将滤波后的信号x进行回放。x=fft(x); %将x信号进行快速傅里叶变换。axes(handles.axes3);plot(x); %在axes3坐标轴上画出x信号的波形图。axes(handles.axes4);plot(abs(x); %在axes4坐标轴上画出x信号的频谱图。h,f=freqz(num,den,256,fc); %求滤波器的频率响应。axes(handles.axes5);plot(f,abs(h),'k'); %在axes5坐标轴上以黑线画出滤波器的频率响应。2）“高通”按键的回调函数function pushbu

11、tton2_callback(hobject, eventdata, handles)y,fs,bits=wavread('e:hy.wav');fp=2000;fs=1000;rp=0.5;rs=40;fc=40000;wp=2*fp/fc;ws=2*fs/fc;if get(handles.radiobutton1,'value')n,wc=buttord(wp,ws,rp,rs);num,den=butter(n,wc,'high'); %返回n阶高通滤波器。else n,wc=cheb1ord(wp,ws,rp,rs);num,den=c

12、heby1(n,rp,wc,'high');endx=filter(num,den,y);sound(x,fc,bits);x=fft(x);axes(handles.axes3);plot(x);axes(handles.axes4);plot(abs(x);h,f=freqz(num,den,256,fc);axes(handles.axes5);plot(f,abs(h),'k');3）“带通”按键的回调函数function pushbutton3_callback(hobject, eventdata, handles)y,fs,bits=wavread

13、('e:hy.wav');fp=3000,8000;fs=1000,10000;rp=0.5;rs=40;fc=40000;wp=2*fp/fc;ws=2*fs/fc;if get(handles.radiobutton1,'value')n,wc=buttord(wp,ws,rp,rs);num,den=butter(n,wc); %wc为双元素向量，返回2n阶带通滤波器。else n,wc=cheb1ord(wp,ws,rp,rs);num,den=cheby1(n,rp,wc);endx=filter(num,den,y);sound(x,fc,bits)

14、;x=fft(x);axes(handles.axes3);plot(x);axes(handles.axes4);plot(abs(x);h,f=freqz(num,den,256,fc);axes(handles.axes5);plot(f,abs(h),'k');4）“带阻”按键的回调函数function pushbutton4_callback(hobject, eventdata, handles)y,fs,bits=wavread('e:hy.wav');fp=1000,10000;fs=3000,8000;rp=0.5;rs=40;fc=40000

15、;wp=2*fp/fc;ws=2*fs/fc;if get(handles.radiobutton1,'value')n,wc=buttord(wp,ws,rp,rs);num,den=butter(n,wc,'stop'); %返回2n阶带阻滤波器。else n,wc=cheb1ord(wp,ws,rp,rs);num,den=cheby1(n,rp,wc,'stop');endx=filter(num,den,y);sound(x,fc,bits);x=fft(x);axes(handles.axes3);plot(x);axes(handl

16、es.axes4);plot(abs(x);h,f=freqz(num,den,256,fc);axes(handles.axes5);plot(f,abs(h),'k');5）“原始信号”按键的回调函数function pushbutton5_callback(hobject, eventdata, handles)y,fs,bits=wavread('e:hy.wav');sound(y,fs,bits);y=fft(y);axes(handles.axes1);plot(y);axes(handles.axes2);plot(abs(y);6）“关闭窗口”

17、按键的回调函数function pushbutton6_callback(hobject, eventdata, handles)close6、 保存修改后的hy.m文件，单击gui主窗口工具栏中的“激活运行”按钮，在gui界面中按下“原始信号”按键就可以看到原始信号的波形和频谱，并播放原始信号；选择滤波器种类（巴特沃思或切比雪夫），然后按下“高通”、“低通”、“带通”、“带阻”其中一个按键，就可以看到原信号经过滤波器后的信号的波形与频谱，并播放该信号。三、调试分析图1 控件布局图2 单击“激活运行”后界面图3 巴特沃思高通滤波器图4 切比雪夫高通滤波器图5 巴特沃思低通滤波器图6 切比雪夫低

18、通滤波器 图7 巴特沃思带通滤波器 图8 切比雪夫带通滤波器 图9 巴特沃思带阻滤波器图10 切比雪夫带阻滤波器图3图10为巴特沃思和切比雪夫型滤波器在不同频率特性下产生的信号波形与频谱，同时有原信号作对比。可以看出，巴特沃思滤波器的频率特性曲线无论在通带还是在阻带内，都是频率的单调函数，切比雪夫型滤波器的幅频特性在通带内是等波纹的，阻带内是单调的。因为在回调函数中，高通和低通、带通和带阻滤波器的截止频率、波纹系数、抽样频率设定的值是相同的，所以在下图中纵向对比可以看出，信号通过高通和低通滤波器之后频谱相对于原信号是互补的，如果对于高通滤波器在某一频率段上有幅值，那么对于低通滤波器该频率幅值为

19、零。同理于带通和带阻滤波器。横向上对比可以看出不同滤波器选频特性也不同，特别是在图7图10中选频边界处的频谱处理可以看出，切比雪夫滤波器比巴特沃思滤波器的选频特性好。四、结果分析与体会滤波器的作用就是滤除掉阻带内有频率分量，保留通带内的频率分量。在图3图10的波形图和频谱图中可以看出，经过滤波器处理后，不仅滤除掉了原信号的部分频率分量，原信号的幅度也会有所衰减，而低通和带阻滤波器产生的信号相对于原信号变化不大，是因为设定的选频范围是原信号的频率段中较集中的部分，所以与此同时，高通和带通滤波器产生的信号相对于原信号的幅度就会大大衰减，其原因也是在于高通和带通滤波器的选频范围是原信号频率范围中较少的一部分。在信号回放时，低通和带阻滤波器回放的音频信号与原信号相比，音色有一点低沉；而高通和带通滤波器回放的音频信号，几乎已经听不到什么声音，这也说明了这两种滤波器产生的信号频率只占原信号频率范围中较少的部分。所以根据各个滤波器的选频范围及滤波后信号的音效，可大致判断原信号的频率多集中在1000hz