124
FFErrVectorCurrent=FFErrVector(1:i);
[I1,J1]=min_matrixel_search(ErrMatrixCurrent);
ErrCurrentMin=ErrMatrixCurrent(I1,J1);
ErrTreshold=(1+0.05)*ErrCurrentMin;
anniter=1:i;
AnniterMin=min(anniter(min(ErrMatrixCurrent,[],2)<=ErrTreshold));
I0=AnniterMin;
ErrMatrixClip=ErrMatrixCurrent(AnniterMin,:);
[AbsoluteErrMin,J0]=min(ErrMatrixClip);
best_NEURON=n_NEURON(I0);
best_Q=n_Q(J0);
if gui
cla(H1)
yGlMin=min(min(min(ErrMatrixCurrent)),min(FFErrVectorCurrent));
yGlMax=max(max(max(ErrMatrixCurrent)),max(FFErrVectorCurrent));
dQ=yGlMax-yGlMin;
if dQ<1e-8
set(H1,'ylim',[yGlMin-1e-8 yGlMax+1e-8])
else
set(H1,'ylim',[yGlMin-0.1*dQ yGlMax+0.1*dQ])
end
if i<2
set(H1,'xlim',[min(n_NEURON(1:i))-1 max(n_NEURON(1:i))+1])
else
dn=max(n_NEURON(1:i))-min(n_NEURON(1:i));
set(H1,'xlim',[min(n_NEURON(1:i))-0.1*dn
max(n_NEURON(1:i))+0.1*dn])
end
xIntv=get(H1,'xlim');