calculate pseudo inverse.m

来自「协同神经网络的程序」· M 代码 · 共 43 行

clear all
clc

NumOfImage = 34;
width=32;
hight=32;

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

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

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

for k=1:NumOfImage
   %伪逆法
   A(k)=VV(k,:)*Q;
end

for i=1:NumOfImage
   for j=1:NumOfImage 
      a(i,j) = v(i,:)*VV(j,:)';
   end    
end