load -ascii svdmat.mat
[u,s,v]=svd(svdmat);

% Reduce reconstruction dimensionality
u_red=u(:,1:2);
s_red=s(1:2,1:2);
v_red=v(:,1:2);

% Reconstruction
recon=u_red*s_red*v_red';

% Similarities
b=s_red*v_red';
temp=sqrt(sum(b.*b));
e = b./[temp;temp];
correla=e'*e;
