《数字语音处理》实验指导书2015.docVIP

《数字语音处理》 实验指导书 主编 王都生 大连民族大学信息与通信工程学院 2015年10月 目 录 实验一 汉明窗加窗的时域波形以及频谱的显示...........................................3? 实验二?WAV格式文件的录制及其短时能量、短时平均幅度分析.............5? 实验三 语音端点检测......................................................................................7? 实验四 语音信号的倒谱的编程实现.............................................................11? 实验五 基于短时自相关函数的基音周期估计.............................................15 实验六 基于短时自相关函数的基音周期估计续.........................................19? 实验七 声道冲激响应及激励信号的分离实现.............................................23?  实验一 汉明窗加窗的时域波形以及频谱的显示 一、实验环境 1．Windows XP，Windows 98，Windows2000 都行。 2．Matlab 二、实验目的 1． 掌握加窗时域波形显示程序的编写。 2．了解清音和浊音的时域波形特征。 3．了解清音和浊音的频谱特征。 三、实验原理 1．窗函数 a 汉明窗函数 b 矩形窗函数 2. 给输入语音信号加窗 3.?MATLAB?程序设计? 分为两步，对于浊音的时域波形显示及其频谱，对于清音的时域波形及其频谱的计算和显示。? 其中，语音信号必须放在该程序所在的文件夹下。? ％汉明窗的短时傅里叶变换 ?[y,Fs,bits]=wavread(1.wav)?? %读出信号，采样率和采样位数。 ?y=y(:,1)?? %这里假设声音是双声道，只取单声道作分析。 ?y=y(4000:4400)? %分别这里选择浊音信号和清音信号分别进行讨论。 y=y.*hamming(length(y))?? figure?? subplot(2,1,1)?? plot(y)?? sigLength=length(y)?? Y=zeros(Fs,1)?? Y?=?20*log10(abs(fft(y,sigLength))) subplot(2,1,2) plot(Y)? % 一点对应 1HZ xlabel(Frequency(Hz))? ylabel(20lg|Xn(ejw)|)? 四、 实验步骤 1．汉明窗加窗的时域波形及其短时频谱显示的源代码。 ①录音 wavrecord (N,Fs,CH) ②读语音 wavread(“sound.wav”) ③截取清浊音 ④加窗 y=y1.*hamming(length(y1)) 2．利用编写的源代码分别显示浊音和清音的时域波形及频谱。 求频谱 20lg|xn(e jw )| 五、实验报告要求 总结语音信号的频谱特征，并指出加矩形窗和汉明窗的浊音信号、清音信号的频谱特征。 实验二 WAV格式文件的录制及其短时能量、短时平均幅度分析 一、实验环境 1、Windows XP，Windows 98，Windows 2000 都行。 2、Matlab 二、实验目的 1、掌握利用MATLAB 录制 WAV格式的语音文件编程。 2、掌握清音和浊音的短时能量和短时平均幅度计算的编程。 3、了解清音和浊音在不同的窗函数下的短时能量和短时平均幅度的特点。 4、了解在不同窗长的情况下，语音信号的短时能量及平均幅度的变化的规律。 三、实验原理 1、短时能量定义为： 其中，N为窗长。可见，短时能量为一帧样点值的加权平方和。 若令：， 则 短时能量可以理解为：语音信号各个样点值平方，然后样点通过一个冲激响应为 h(n)的滤波器，输出为由短时能量构成的时间序列。因此，冲激响应 h(n)的选择或者说窗函数的选择直接影响着短时能量的计算。 短时平均幅度函数 Mn用来衡量语音幅度的变化