语音信号特征提取程序设计代码Word文件下载.docx

1、n=0:N-1;y=fft（x）;mag=abs（y）;f=（0:length（y）-1）*fs/length（y）;%进行对应的频率转换subplot（212）;plot（f,mag）;原始信号频谱图频率（Hz）y=wavread（%预加重figure（）;原始波形（幅值）xx=double（y）;plot（xx）;预加重波形（幅值）%对指定帧位置进行加窗处理Q=y;N=256;%窗长Hamm=hamming（N）;%加窗（汉明窗）Rect=rectwin（N）;%加窗（矩形窗）frame=60;%需要处理的帧位置M=Q（frame-1）*（N/2）+1）:（frame-1）*（N/2）+N

2、）;Frame=M.*Hamm%加窗后的语音帧frame=M.*Rectplot（Frame）;汉明窗plot（frame）;矩形窗%语音信号分帧plot（y）;原始y=enframe（y,256,80）;帧数分帧%端点检测x1,x2=vad（x）;%调用函数vad实现%浊音， 取13270-13510个点%短时自相关函数temp=y（13271:13510）;Rn1=zeros（1,240）;for nn=1:240 for ii=1:240-nn Rn1（nn）=Rn1（nn）+temp（ii）*temp（nn+ii）; endendfigure（6）;jj=1:240;plot（jj,R

3、n1,b浊音自相关函数短时自相关函数%清音， 取12120-12360个点temp=y（12121:12360）;Rn2=zeros（1,240）; Rn2（nn）=Rn2（nn）+temp（ii）*temp（nn+ii）;plot（jj,Rn2,清音自相关函数分帧功能子程序：function f=enframe（x,win,inc） % F = ENFRAME（X,LEN） splits the vector X（:） up into % frames. Each frame is of length LEN and occupies % one row of the output matr

4、ix. The last few frames of X % will be ignored if its length is not divisible by LEN. % It is an error if X is shorter than LEN. nx=length（x（:）; nwin=length（win）; if （nwin = 1） len = win; else len = nwin; end if （nargin 1） w = win（:） f = f .* w（ones（nf,1）,: End选端点功能子程序：function x1,x2 = vad（x）% x=aud

5、ioread（1.wav%幅度归一化到-1,1x = double（x）;x = x / max（abs（x）;%常数设置FrameLen = 240;FrameInc = 80;amp1 = 10;amp2 = 2;zcr1 = 10;zcr2 = 5;maxsilence = 8; % 6*10ms = 30msminlen = 15; % 15*10ms = 150msstatus = 0;count = 0;silence = 0;%计算过零率tmp1 = enframe（x（1:end-1）, FrameLen, FrameInc）;tmp2 = enframe（x（2:end） ,

6、 FrameLen, FrameInc）;signs = （tmp1.*tmp2）0.02;zcr = sum（signs.*diffs, 2）;%计算短时能量amp = sum（abs（enframe（filter（1 -0.9375, 1, x）, FrameLen, FrameInc）, 2）;%调整能量门限amp1 = min（amp1, max（amp）/4）;amp2 = min（amp2, max（amp）/8）;%开始端点检测x1 = 0;x2 = 0;for n=1:length（zcr） goto = 0; switch status case 0,1 % 0 = 静音,

7、1 = 可能开始 if amp（n） amp1 % 确信进入语音段 x1 = max（n-count-1,1）; status = 2; silence = 0; count = count + 1; elseif amp（n） amp2 | . % 可能处于语音段 zcr（n） zcr2 status = 1; count = count + 1; else % 静音状态 status = 0; count = 0; end case 2, % 2 = 语音段 amp2 | . % 保持在语音段 zcr（n） count = count + 1; else % 语音将结束 silence =

8、 silence+1; if silence maxsilence % 静音还不够长，尚未结束 count = count + 1; elseif count minlen % 语音长度太短，认为是噪声 status = 0; silence = 0; count = 0; else % 语音结束 status = 3; end case 3, break; endend count = count-silence/2;x2 = x1 + count -1;subplot（311）plot（x）axis（1 length（x） -1 1）Speechline（x1*FrameInc x1*FrameInc, -1 1, Color, redline（x2*FrameInc x2*FrameInc, -1 1, subplot（312）plot（amp）;axis（1 length（amp） 0 max（amp）Energyline（x1 x1, min（amp）,max（amp）, line（x2 x2, min（amp）,max（amp）, subplot（313）plot（zcr）;axis（1 length（zcr） 0 max（zcr）ZCRline（x1 x1, min（zcr）,max（zcr）, line（x2 x2, min（zcr）,max（zcr）,