f7_4.m

matlab时频分析配套源码第6章的程序

%读入语音文件
fin=fopen('DR4_MLJH0_SX334.ADC','r');
x=fread(fin,'short');
fclose(fin);
%时域波形
figure(1);
plot(x,'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');

%计算短时能量
N=length(x);
%%短时窗长
L=200;
for i=1:N-L
    e(i)=energy(x(i:i+L));
end
figure(2);

%%能量图
plot(e,'LineWidth',2);
xlabel('时间 t');
ylabel('能量 E');
%energy.m
function [e] = energy(s);
e=sum(s.*s)/max(size(s));