[x]=wavread('c101.wav');
fclose(x);
figure(1);
plot(x,'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
N=length(x);
L=200;
for i=1:N-L;
    e(i)=engery(x(i:i+L));
end
figure(2);
plot(e,'LineWidth',2);
xlabel('時間 t');
ylabel('能量 E');
function[e]=energy(s);
e=sum(s.*s)/max(size(s));