clear?all
close?all
%讀入訓練圖像，圖像大小為N=100*100，訓練圖像個數為M
alltrainimg=[];
for?i=1:5
????for?j=1:1:10
??????filename??=?sprintf（‘G:\\faces\\%d\\s%d.bmp‘ij）;
??????a?=?imread（filename）;
??????b=a（1:100*100）;?
??????b=double（b）;
??????alltrainimg?=[alltrainimg;?b];??%?alltrainimg?是一個M?*?N?矩陣，alltrainimg?中每一行數據代表一張圖片
??end
end
%?for?i=1:99
%??????for?j=1:99
%??????????X（ij）=alltrainimg（100i*100+j）;
%??????end
%??end
%??imshow（X[]）;

img_mean=mean（alltrainimg）;?%?求出訓練集對應的平均圖片，為1?*?N?矩陣

M?=?5*10;
X_mean=zeros（M10000）;
for?i=1:M
????X_mean（i:）=alltrainimg（i:）-img_mean;?%?X_mean每一行保存的數據是“每個圖片數據-平均圖片”
end;

?%計算KL變換生成矩陣sigma
?sigma1=X_mean*X_mean‘;???%?M?*?M?階矩陣
?
?%求出sigma的特征值與特征向量，并將其進行排序
??[v?d]=eig（sigma1）;
??d1=diag（d）;
??[d2?index]=sort（d1）;??%以升序排序
??col

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-08-10?16:18??face?recognition\
?????文件????????2313??2011-08-10?16:05??face?recognition\PCAtest.m
?????目錄???????????0??2011-08-03?14:50??face?recognition\Yale_5G\
?????文件?????1221330??2007-07-25?20:30??face?recognition\Yale_5G\allcoor.mat
?????文件?????1221330??2007-07-25?20:30??face?recognition\Yale_5G\base.mat
?????文件??????744959??2011-05-08?15:54??face?recognition\Yale_5G\base1.mat
?????目錄???????????0??2011-08-03?14:50??face?recognition\Yale_5G\face1\
?????文件????????1381??2008-03-18?14:50??face?recognition\Yale_5G\face1\bilinear.m
?????文件?????????306??2008-03-18?15:49??face?recognition\Yale_5G\face1\ch.m
?????文件?????1174815??2008-03-18?15:50??face?recognition\Yale_5G\face1\t.rar
?????文件????????3463??2011-05-08?17:36??face?recognition\Yale_5G\facerec.asv
?????文件????????3463??2011-05-08?17:36??face?recognition\Yale_5G\facerec.m
?????文件????????2003??2007-09-09?21:56??face?recognition\Yale_5G\file.txt
?????文件????????5000??2007-09-09?21:58??face?recognition\Yale_5G\file1.txt
?????文件????????4076??2011-05-08?17:37??face?recognition\Yale_5G\file2.txt
?????文件???????18024??2011-05-08?15:52??face?recognition\Yale_5G\hs_err_pid2268.log
?????文件???????33488??2011-05-08?16:00??face?recognition\Yale_5G\hs_err_pid4016.log
?????文件????????1451??2007-09-09?10:51??face?recognition\Yale_5G\kernel.m
?????文件????????2153??2011-05-08?15:55??face?recognition\Yale_5G\rec.asv
?????文件????????2161??2011-05-08?16:03??face?recognition\Yale_5G\rec.m
?????文件????????8206??2007-09-09?10:51??face?recognition\Yale_5G\svm168.m
?????目錄???????????0??2011-08-03?14:50??face?recognition\Yale_5G\t\
?????文件????????1033??2011-05-08?17:03??face?recognition\Yale_5G\test.m
?????文件??????744959??2011-05-08?15:54??face?recognition\Yale_5G\trainvec.mat
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s1.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s10.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s11.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s2.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s3.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s4.bmp
?????文件????????7478??2008-03-18?15:49??face?recognition\Yale_5G\t\10_s5.bmp
............此處省略339個文件信息