pca 人脸识别+orl人脸数据库

修改过的pca人脸识别程序+orl 人脸库，点击facrec就能看见结果，路径是相对路径，

% eigface.m
function [] = eigface（）
% calc xmeansigma and its eigen decomposition
allsamples=[];%所有训练图像
for i=1:40
for j=1:5
a=imread（strcat（‘.\ORL人脸数据库\s‘num2str（i）‘\‘num2str（j）‘.bmp‘））;
% imshow（a）;
b=a（1:112*92）; % b 是行矢量 1×N，其中N＝10304，提取顺序是先列后行，即从上到下，从左到右
b=double（b）;
allsamples=[allsamples; b]; % allsamples 是一个M * N 矩阵，allsamples 中每一行数据代表一张图片，其中M＝200
end
end
samplemean=mean（allsamples）; % 平均图片，1 × N
for i=1:200 xmean（i:）=allsamples（i:）-samplemean; % xmean 是一个M × N 矩阵，xmean每一行保存的数据是“每个图片数据-平均图片”
end;
% 获取特征值及特征向量
sigma=xmean*xmean‘; % M * M 阶矩阵
[v d]=eig（sigma）;
d1=diag（d）;
% 按特征值大小以降序排列
dsort = flipud（d1）;
vsort = fliplr（v）;
%以下选择90%的能量
dsum = sum（dsort）;
dsum_extract = 0;
p = 0;
while（ dsum_extract/dsum < 0.9）
p = p + 1;
dsum_extract = sum（dsort（1:p））;
end
p = 199;
% （训练阶段）计算特征脸形成的坐标系
%while （i<=p && dsort（i）>0）
% base（:i） = dsort（i）^（-1/2） * xmean‘ * vsort（:i）; % base 是N×p 阶矩阵，除以dsort（i）^（1/2）是对人脸图像的标准化，详见《基于PCA 的人脸识别算法研究》p31
% i = i + 1; % xmean‘ * vsort（:i）是小矩阵的特征向量向大矩阵特征向量转换的过程
%end
base = xmean‘ * vsort（:1:p） * diag（dsort（1:p）.^（-1/2））;
% 生成特征脸
figure 
for （k=1:p）
temp = reshape（base（:k） 11292）;
subplot（1pk）
imshow（temp）
newpath = [‘.\test\‘ int2str（k） ‘.bmp‘];
imwrite（mat2gray（temp） newpath）;
end
avg = reshape（samplemean 11292）;
imwrite（mat2gray（avg） ‘.\test\average.bmp‘）;
% 将模型保存
save（‘.\model.mat‘ ‘base‘ ‘samplemean‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1638  2014-04-14 23:07  eigface.m
     文件        2784  2014-04-14 22:52  FaceRec.m
     目录           0  2014-04-14 23:09  ORL人脸数据库\
     目录           0  2014-04-11 12:42  ORL人脸数据库\s1\
     目录           0  2014-04-11 12:42  ORL人脸数据库\s10\
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\1.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s10\1.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\10.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s10\10.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\2.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s10\2.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\3.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s10\3.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\4.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s10\4.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\5.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s10\5.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\6.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s10\6.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\7.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s10\7.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\8.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s10\8.pgm
     文件       11384  2010-03-25 19:23  ORL人脸数据库\s10\9.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s10\9.pgm
     目录           0  2014-04-11 12:42  ORL人脸数据库\s11\
     文件       11384  2010-03-25 19:25  ORL人脸数据库\s11\1.bmp
     文件       10318  1994-04-18 14:18  ORL人脸数据库\s11\1.pgm
     文件       11384  2010-03-25 19:25  ORL人脸数据库\s11\10.bmp
     文件       10318  1994-04-18 14:07  ORL人脸数据库\s11\10.pgm
     文件       11384  2010-03-25 19:25  ORL人脸数据库\s11\2.bmp
............此处省略1004个文件信息