1_短时能量幅度_audition_

%%短时能量 [x,Fs]=audioread('E:\Matlab\bin\homework\voice.wav'); % 读入数据文件 wlen=20; inc=10; % 给出帧长和帧移 win=rectwin(wlen); % 给出海宁窗 N=length(x); % 信号长度 X=enframe(x,win,inc)'; % 分帧 fn=size(X,2); % 求出帧数 time=(0:N-1)/Fs; % 计算出信号的时间刻度 for i=1 : fn u=X(:,i); % 取出一帧 u2=u.*u; % 求出能量 En(i)=sum(u2); % 对一帧累加求和 end subplot(2,1,1); plot(time,x); % 画出时间波形 title('语音波形'); ylabel('幅值'); xlabel(['时间/s' 10 '(a)']); frameTime=(((1:fn)-1)*inc+wlen/2)/Fs; % 求出每帧对应的时间 subplot(2,1,2); plot(frameTime,En) % 画出短时能量图 title('短时能量'); ylabel('幅值'); xlabel(['时间/s' 10 '(b)']);