function??[Out]=facerec（Trainset）

Slect=Trainset;
n=1;
M=45;
N=length（Slect）;
allsamples=[];%所有訓練圖像
for?i=1:1:2
????for?j=1:N
????????filename??=?sprintf（‘D:\\zzh\\work\\zzh\\baseline_V8\\matlab\\Yale_5G\\face1\\%d_s%d.bmp‘iSlect（j））;
??????a=imread（filename）;?????
??????a=a‘;??
??????b=a（1:80*80）;??%?b是行矢量?1×N，其中N＝10000，提取順序是先列后行，即從上到下，從左到右
??????b=double（b）;
??????allsamples=[allsamples;?b]??%?allsamples?是一個M?*?N?矩陣，allsamples?中每一行數據代表一張圖片，其中M＝200
????end
end

samplemean=mean（allsamples）;?%?平均圖片

xmean=[];
for?i=1:1:M?
????xmean（i:）=allsamples（i:）-samplemean;?%?xmean是一個M?×?N矩陣，xmean每一行保存的數據是“每個圖片數據-平均圖片”
end;
xmean

sigma=xmean*xmean‘;???%?M?*?M?階矩陣
sigma;
[v?d]=eig（sigma）;
d1=diag（d）;
[d2?index]=sort（d1）;?%以升序排序
cols=size（v2）;%?特征向量矩陣的列數
vsort=[];

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3406??2008-03-18?15:19??PCA-SVM\facerec.m

?????文件???????5000??2007-09-09?21:58??PCA-SVM\file1.txt

?????文件???????2003??2007-09-09?21:56??PCA-SVM\file.txt

?????文件???????1451??2007-09-09?10:51??PCA-SVM\kernel.m

?????文件???????1618??2007-09-09?11:00??PCA-SVM\rec.asv

?????文件???????1646??2007-09-09?11:01??PCA-SVM\rec.m

?????文件???????1032??2007-09-09?21:53??PCA-SVM\test.m

?????文件???????1991??2008-01-23?08:56??PCA-SVM\WavePCA.m

?????文件???????9883??2008-01-23?10:52??PCA-SVM\facerec.asv

?????文件???????8206??2007-09-09?10:51??PCA-SVM\svm168.m

?????目錄??????????0??2009-05-16?03:08??PCA-SVM

-----------?---------??----------?-----??----

????????????????36236????????????????????11