clear;
clc;
af=[1.24 1.36 1.38 1.38 1.38 1.4 1.48 1.54 1.56;
    1.72 1.74 1.64 1.82 1.90 1.7 1.82 1.82 2.08];
apf=[1.14 1.18 1.20 1.26 1.28 1.30;
     1.78 1.96 1.86 2.00 2.00 1.96];
 x=[af apf];
 y0=[2*ones(1,9) ones(1,6)];
net=newff(minmax(x),[15,1],{'tansig','purelin'});
 net.trainFcn='trainlm';
 net.trainParam.goal=1e-6;
  net.trainParam.epochs=1000;
  net=train(net,x,y0);
  ytmp=sim(net,x);
  y1=vec2ind(ytmp);
  xt=[1.24 1.28 1.40;1.80 1.84 2.04];
  yttmp=sim(net,xt);
  yt=vec2ind(yttmp)
  figure(1);
  plot(af(1,:),af(2,:),'+',apf(1,:),apf(2,:),'o',xt(1,:),xt(2,:),'*');
  legend('af','apf','檢測樣本')