%%============================================
%%         協(xié)同識(shí)別系統(tǒng)
%%============================================
clc
%訓(xùn)練
Par=[ones(1,4) 1 1];


Time0=clock;

l=64;w=64;

S3=example8Im2Vec('example81sample1.bmp');
S4=example8im2vec('example81sample2.bmp');
S0=example8im2vec('example81sample3.bmp');


v3=S3-mean(S3);v4=S4-mean(S4);v0=S0-mean(S0);
V=[v3/norm(v3) v4/norm(v4) v0/norm(v0)];

P=pinv(V);



%測(cè)試
Test=example8im2vec('example81sample4.bmp');

q=Test;
q=q-mean(q);
q=q/norm(q);

Order0=P*q;
V;

[N,M]=size(V);

for i=1:M
    p(i,1)=Order0(i);
end

B=1;
C=1;

steps=30;

for n=1:steps                                            
       
       D=0;
       for j=1:M
           D=D+p(j,n)^2;
       end
       
       D=(B+C)*D;
       gamma=1/D;

       for k=1:M
           p(k,n+1)=p(k,n)+gamma*(Par(k)-D+B*(p(k,n)^2))*p(k,n);
       end
end


end

L=l;W=w;