资源描述:
《特征提取》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、“这下全完啦”高兴语调短时能量“这下全完啦”悲伤语调短时能量“这下全完啦”高兴语调短时平均能量“这下全完啦”悲伤语调短时平均能量“这下全完啦”高兴语调短时过零率“这下全完啦”悲伤语调短时过零率“这下全完啦”高兴语调清音短时自相关函数悲伤语调清音短时自相关函数“这下全完啦”高兴语调浊音短时自相关函数悲伤语调浊音短时自相关函数>>closeall;clearall;>>[filename,filepath]=uigetfile('.wav','Openwavfile');y=wavread([filepath,filename]);>>[x,fs,bit]=wavrea
2、d('AH.wav');%读语音数据,数据放入了X,fs代表采样率,bit代表采样位数>>framelength=240;%设置帧长,%30msunder8khzframenumber=fix(length(x)/(framelength));%总的数据帧数totaltime=length(x)/fs;%该段语音总的时间长度sp1=x(fix(1.1*fs):(fix(1.1*fs)+framelength+framelength-1));%1.1*fs=8.8秒处为浊音始点,取两帧sp2=x(0.8*fs:(0.8*fs+framelength+frameleng
3、th-1));%0.8*fs=6.4秒处为清音始点,取两帧d=0;%初始化,浊音波峰判断R1=0;%初始化,浊音一个基音周期内帧数%计算每帧的能量fori=1:framenumber;E(i)=0;%短时能量初始化Z(i)=0;%短时过零率初始化M(i)=0;%短时平均幅度初始化K(i)=0;%短时零能比初始化forj=(framelength*(i-1)+1):framelength*i;E(i)=E(i)+x(j)*x(j);%第i帧短时能量endforj=(framelength*(i-1)+2):framelength*i;Z(i)=Z(i)+abs(sig
4、n(x(j)-0.2)-sign(x(j-1)));%第i帧短时过零率endforj=(framelength*(i-1)+1):framelength*i;M(i)=M(i)+abs(x(j));endM(i)=M(i)/framelength;%第i帧短时平均幅度endfori=1:framelength;R(i)=0;forj=1:framelength;R(i)=R(i)+sp1(j)*sp1(i+j);%浊音的短时自相关函数endif(R(i)>=110.11)R(i)=110.11;%修正最高波峰最大值为110.11d=d+1;%当波峰为110.11时,
5、自加一endif(d==1)%以第一个波峰110.11为起始点,下一个波峰110.11为终点R1=R1+1;%计算一个基音周期内帧数endendfori=1:framelength;r(i)=0;forj=1:framelength;r(i)=r(i)+sp2(j)*sp2(i+j);%清音的自相关函数endend%思考为什么取sign(x(j)-0.2)%把数据点转化为时间长度fori=1:framenumber;anyfrequency(i)=totaltime*i/framenumber;endfori=1:length(x);anytime(i)=i*tot
6、altime/length(x);endT=R1/fs;%浊音段语音的基音周期,单位为秒(s)figure(1)subplot(2,1,1);plot(anytime,x);title('语音图');%语音图subplot(2,1,2);plot(anyfrequency,E);title('短时能量');%短时能量figure(2)subplot(2,1,1);plot(anytime,x);title('语音图');%语音图subplot(2,1,2);plot(anyfrequency,M);title('短时平均幅度');%短时平均幅度figure(3)su
7、bplot(2,1,1);plot(anytime,x);title('语音图');%语音图subplot(2,1,2);plot(anyfrequency,Z);title('短时过零率');%短时过零率figure(4)subplot(3,1,1);plot(anytime,x);title('语音图');%语音图subplot(3,1,2);plot(anyfrequency,Z);title('短时过零率');%短时过零率subplot(3,1,3);plot(anyfrequency,E);title('短时能量');%短时能量figure(5)sub