语音信号的短时分析matlab仿真,包括短时过零率,短时能量,短时自相关等【包含程序操作录像,代码中文注释】

clear all; clc; close all; [x,Fs]=audioread('C3_2_y.wav'); % 读入数据文件 wlen=200; inc=100; % 给出帧长和帧移 win=hanning(wlen); % 给出海宁窗 N=length(x); % 信号长度 time=(0:N-1)/Fs; % 计算出信号的时间刻度 En=STEn(x,win,inc); %短时能量 Mn=STMn(x,win,inc); %短时平均幅度 Zcr=STZcr(x,win,inc); %短时过零率 %此处和上述3个参数不同，返回的不是向量而是矩阵，因为一帧信号得到的不是一个数值 X=enframe(x,win,inc)'; % 分帧 xn=X(:); Ac=STAc(X); %计算短时自相关 Ac=Ac(:); Amdf=STAmdf(X); %计算短时幅度差 Amdf=Amdf(:); fn=length(En); % 求出帧数 figure(1) subplot 311; plot(time,x,'b'); axis tight% 画出时间波形 title('(a)语音波形'); ylabel('幅值'); xlabel(['时间/s' 10 ]); frameTime=FrameTimeC(fn,wlen,inc,Fs); % 求出每帧对应的时间 subplot 312; plot(frameTime,Mn,'b') % 画出短时幅度图 title('(b)短时幅度'); ylabel('幅值'); xlabel(['时间/s' 10 ]); subplot 313; plot(frameTime,En,'b') % 画出短时能量图 title('(c)短时能量'); ylabel('幅值'); xlabel(['时间/s' 10 '(b)']); figure(2) subplot 211; plot(time,x,'b'); axis tight% 画出时间波形 title('(a)语音波形'); ylabel('幅值'); xlabel(['时间/s' 10 ]); subplot 212; plot(frameTime,Zcr,'b') % 画出短时过零率图 title('(b)短时过零率'); ylabel('幅值'); xlabel(['时间/s' 10 ]); figure(3) subplot 211; plot(xn,'b'); % 画出时间波形 title('(a)语音波形'); ylabel('幅值'); xlabel(['点数' 10 ]); subplot 212; plot(Ac,'b') % 画出短时自相关图 title('(b)短时自相关'); ylabel('幅值'); xlabel(['点数' 10 ]); figure(4) subplot 211; plot(xn,'b'); % 画出时间波形 title('(a)语音波形'); ylabel('幅值'); xlabel(['点数' 10 ]); subplot 212; plot(Amdf,'b') % 画出短时幅度差 title('(b)短时幅度差'); ylabel('幅值'); xlabel(['点数' 10 ]);