plot faces and eigenface

main.m

close all;
clear;
clc;
X=[];
for i=1:10,
    x=imread([int2str(i) '.pgm']);
    [w,h]=size(x);
    x=reshape(x,w*h,1);
    x=double(x);
    X=[X x];
end

for i=1:10,
    subplot(2,5,i);
    I=X(:,i);
    I=uint8(I);
    u=reshape(I,w,h);
    imshow(u);
end
hold on;
figure;
[U,S]=pca(X);
for i=1:10,
    subplot(2,5,i);
    I=U(:,i);
    minI=min(I);
    maxI=max(I);
    I=(I-minI)/(maxI-minI);
    I=I*255;
    I=uint8(I);
    u=reshape(I,w,h);
    imshow(u);
    colormap(jet(256));
    title(sprintf('%f',S(i,i)));
end