ft=['babystar1.wav';'corn1.wav';'babeq1.wav'];%wK鉹t@CׂĂ
fnt=size(ft)(1,1);%t@C̐
d=32;%
c = 1000;%wKf[^Ɏgp鉹t@C̒
nt=1000;%wKf[^̐
Yt=zeros(d,1);%tf[^p̍s
X=zeros(fnt,1);%wKf[^p̍s
for i=1:fnt%wK鉹t@C̐JԂ
  for j=10:nt+10%et@C؂oĕϊ
    Y=log10(abs(fft(wavread(ft(i,:),[j*c+1 (j+1)*c]),d*2)));
    Yt=[Yt,Y([1:d],1)];
    X=[X eye(fnt)(:,i)];
   end
end
Yt = Yt(:,[2:nt*fnt+1]);
X = X(:,[2:nt*fnt+1]);

Rx2 = min_max(Yt);%wKf[^̍őlƍŏl
S   = [d fnt];%wKp[^̐ݒ
NNet = newff( Rx2, S );
NNet.trainParam.epochs = 500;
VV.P = Yt;
VV.T = X;
net = train( NNet, Yt, X, [], [], VV );%wK
res = sim( net, Yt);%wKʂ̌
a=0;
for i=1:fnt
  for j=1+nt*(i-1):nt*i;
    [x,ix] = max(res(:,j));
    if(i==ix)%tf[^Ɠo̓j[̒lł傫ȂĂ邩
      a++;
    end
  end
end
a/(nt*fnt)%̕\
