%This?program?is?writen?by?zhangming?on?2009-5-10
%this?is?my?homeworke?of?machine?leaning
%my?e-mail:?zm_fred@163.com
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all
clc
close?all
%選擇訓練樣本路徑
TrainDatabasePath?=?uigetdir（‘D:\ORL人臉數據庫‘?‘Select?training?database?path‘?）;
%選擇測試樣本路徑
TestDatabasePath?=?uigetdir（‘D:\ORL人臉數據庫‘?‘Select?test?database?path‘）;
%生成測試數據集
%T?=?CreateDatabase（TrainDatabasePath）;
TrainFiles?=?dir（TrainDatabasePath）;
class_Number?=?0;
s=size（TrainFiles1）;
%計算訓練樣本個數
for?i?=?1:size（TrainFiles1）
????if?not（strcmp（TrainFiles（i）.name‘.‘）|strcmp（TrainFiles（i）.name‘..‘）|strcmp（TrainFiles（i）.name‘README‘））
????????class_Number?=?class_Number?+?1;?%?Number?of?all?images?in?the?training?database
????end
end????????????????????????????????????????%class_Number=40??????
%讀取訓練集第一個圖像，得到其尺寸，初始化圖像矩陣
str?=?strcat（TrainDatabasePath‘\s1\1.bmp‘）;
img?=?imread（str）;
[irowicol]?=?size（img）;
A=zeros（irowicol2）;
%合計有Train_Number人的目錄，每個目錄存放某人的10幅人臉圖像取前n副圖像做訓練樣本
n=5;
for?i?=?1?:?class_Number
????topdir?=?int2str（i）;
????topdir?=?strcat（TestDatabasePath‘s‘topdir）;
????for?j=?1?:?n
?????%for?j=?1?:2?:9
????????str?=?int2str（j）;
????????str?=?strcat（topdir‘‘str‘.pgm‘）;
????????img?=?imread（str）;
????????A（::n*（i-1）+j）?=?img;?
????end
end
%得到訓練樣本的個數Train_Number
Train_Number?=?size（A3）;
%求所有訓練樣本向量的均值?M
M=mean（A3）;

%求圖像散布矩陣Gt
Gt=zeros（icolicol）;
for?i?=?1?:?Train_Number
????temp?=?A（::i）-M;
????Gt?=??Gt?+?temp‘*temp;
end
Gt=Gt/Train_Number;
d?=?10;
[V?D]?=?eigs（Gtd）;
%ＰＣＡ空間投影圖像
%V?=?orth（V）;
ProjectedImages?=?zeros（irowd）;
for?i?=?1?:?Train_Number
????ProjectedImages（::i）?=?A（::i）*V;????
end
%已識別測試集的數目Recnum
Recnum=0
for?i?=?1?:?class_Number
????topdir?=?int2str（i）;
????topdir?=?strcat（TestDatabasePath‘s‘topdir）;
????for?j=?n+1:10
?????%for?j=?2?:2?:10
????????str?=?int2str（j）;
????????TestImage?=?strcat（topdir‘‘str‘.pgm‘）;?%測試樣本路徑
?????????%測試樣本原始數據
????????im?=?imread（TestImage）;
????????%測試樣本的投影
????????ProjectedTestImage?=?double（im）*V;
???????
????????gc=[];
?