function [n,sita,pred]=identi(n0,y,u,t,my,n,m,F,fil,dis,lambda,sita,pred) ty=y(t-1)+my*(u(t)-1); sum=0; uni=ones(1,my); delta=zeros(my,1); delta(y(t))=1; for i=1:my j=i+my*(dis-1)+my*(fil-1)*(m-2); n(j,ty)=lambda*(n(j,ty)+delta(i))+(1-lambda)*n0(j,ty); sum=sum+n(j,ty); end for i=1:my j=i+my*(dis-1)+my*(fil-1)*(m-2); sita(j,ty)=n(j,ty)/sum; end j=y(t)+my*(dis-1)+my*(fil-1)*(m-2); pred(dis+(fil-1)*(m-2))=sita(j,ty);