X = imread('face01.png'); X = double(X) ./ 255.0; [U,S,V] = svd(X); imshow(U(:,1) * S(1,1) * V(:,1)' ); imshow(U(:,1:5)* S(1:5, 1:5) * V(:,1:5)'); imshow(U(:,1:10)* S(1:10, 1:10) * V(:,1:10)');