balancequicknetwork.m

协同神经网络的程序

clear all
clc

NumOfImage = 34;
width=32;
hight=32;
maxGray = 1;

for i = 1: NumOfImage
%   i
   R=imread(num2str(i),'bmp');
   B=(R(:))';
   v(i,:)=double(B);
   vmean(i)=mean(v(i,:));
   v(i,:)=v(i,:)-mean(v(i,:));
   v(i,:)=v(i,:)/norm(v(i,:));
end
%size(v)
%pause

%we get the value of pinv v
VV = (pinv(v))';
%size(VV)

%Q is the row parament
R=imread('q','bmp');
Q=double(R(:));
Qmean= mean(Q);
Q=Q-mean(Q);
Q=Q/norm(Q);

%A is the order parameter
for k=1:NumOfImage
   A(k)=VV(k,:)*Q;
end

max = 1;
flops(0);
tic;
for k=1:NumOfImage
   A(k)=VV(k,:)*Q;
   if(A(k)>A(max)) max = k;end
end
flops
toc

% q is the recognized pattern;Output the value of q and draw it;      
q=v(max,:);  
q = reshape(q,hight,width);
q = mat2gray(q);
imshow (q);