%讀入語音文件
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));