1、实验一 语音信号的短时分析一、实验目的1.在理论学习的基础上，进一步地理解和掌握语音信号短时分析的意义，短时时域分析的基本方法。2.进一步理解和掌握语音信号短时平均能量函数及短时平均过零数的计算方法和重要意义。二、实验原理及方法一定时宽的语音信号，其能量的大小随时间有明显的变化。其中清音段（以清音为主要成份的语音段），其能量比浊音段小得多。短时过零数也可用于语音信号分析中，发浊音时，其语音能量约集中于 3kHz 以下，而发清音时，多数能量出现在较高频率上，可认为浊音时具有较低的平均过零数，而清音时具有较高的平均过零数，因而，对一短时语音段计算其短时平均能量及短时平均过零数，就可以较好地区分其中
2、的清音段和浊音段，从而可判别句中清、浊音转变时刻，声母韵母的分界以及无声与有声的分界。这在语音识别中有重要意义。三、实验仪器微型计算机，Matlab 软件环境四、实验步骤1.上机前用 Matlab 语言完成程序编写工作。2.程序应具有加窗（分帧）、计算 、 以及绘制曲线等功能。3.上机实验时先调试程序，通过后进行信号处理。4.对录入的语音数据进行处理，并显示运行结果。5.依据曲线对该语音段进行所需要的分析，并作出结论。6.改变窗的宽度（帧长），重复上面的分析内容。五、操作步骤所用语音信号文件名为“shop.wav“,拷贝到 MATLAB 工作目录。Matlab 编程实验步骤：1新建 M 文件，
3、扩展名为“.m”，编写程序；2选择 File/Save 命令，将文件保存在 F 盘中；3在 Command Window 窗中输入文件名，运行程序；Matlab 部分函数语法格式：读 wav 文件： x=wavread(filename) 数组 a 及 b 中对应元素相乘： a.*b创建图形窗口命令： figure绘图函数： plot(x) 坐标轴： axis(xmin xmax ymin ymax)坐标轴注解： xlabel() ylabel()图例注解： legend( )一阶高通滤波器： y=filter(1-0.09375,1,x)voicebox 工具箱介绍：分帧函数：f=enfra
4、me(x,len,inc) x 为输入语音信号，len 指定了帧长，inc 指定帧移，函数返回为 nlen 的一个矩阵，每一行都是一帧数据。voicebox 工具箱安装方法： 把工具箱拷贝入.MatLabtoolbox 目录下。菜单 file-set path-“Add folder“ 导入 voicebox； 菜单 file-Preference-General-“Updata Toolbox Path Cache六、实验结果七、思考题1语音信号短时平均能量及短时平均过零数分析的主要用途是什么？可以从语音中区别出浊音来，因为浊音时短时平均能量的值要比清音时短时平均能量的值大很多；可以用来区别
5、声母和韵母的分界、无声和有声的分界等最为一种超音段信息，用于语音识别中。浊音时能量集中于较低频率段内，具有较低的过零率，而清音时能量集中于较高频率段内，具有较高的过零率2窗的宽度（帧长）的改变，对 的特性产生怎样的影响？窗口的长度，无论什么形状的窗口，窗口序列的长度 N 将起决定性的作用。N 选得太大，滤波器的通带变窄，波形的振幅变化细节就看不出来，反之，如果 N 太小，则滤波器的通带变宽，信号得不到足够的平均八、程序x,z=wavread(shop.wav);figure;sound(x,z);subplot(4,1,1);plot(x);axis(1 length(x) -1 1);yla
6、bel(Speech);enhance=filter(1-0.9375,1,x);%高频预加重FrameLen=240;FrameInc=80;yframe=enframe(x,FrameLen,FrameInc);amp1=sum(abs(yframe),2);subplot(4,1,2);plot(amp1);axis(1 length(amp1) 0 max(amp1);ylabel(Amp);legend(amp1=x );amp2=sum(abs(yframe.*yframe),2);subplot(4,1,3);plot(amp2);axis(1 length(amp2) 0 ma
7、x(amp2);ylabel(Energy);legend(amp1=x*x);tmp1=enframe(x(1:end-1),FrameLen,FrameInc);tmp2=enframe(x(2:end),FrameLen,FrameInc);signs=(tmp1.*tmp2)0.02;zcr=sum(signs.*diffs,2);subplot(4,1,4);plot(zcr);axis(1 length(zcr) 0 max(zcr);ylabel(ZCR);legend(zcr);建立分帧函数文件function f=enframe(x,win,inc)nx=length(x);nwin=length(win);len = win;nf = fix(nx-len+inc)/inc);f=zeros(nf,len);indf= inc*(0:(nf-1).;inds = (1:len);f(:) = x(indf(:,ones(1,len)+inds(ones(nf,1),:);
