%用KPCA進行人臉識別的Matlab代碼?
allsamples=[];%所有訓練圖像?
for?i=1:40?
????for?j=6:10
??????%a=imread（strcat（‘Data\s‘num2str（i）‘\‘num2str（j）‘.pgm‘））;?
??????a=imread（strcat（‘C:\Documents?and?Settings\Administrator\桌面\人臉庫\ORL人臉庫\s‘num2str（i）‘\‘num2str（j）‘.jpg‘））;
??????b=a（1:112*92）;?%?b是行矢量?1×N，其中N＝10304，提取順序是先列后行，即從上到下，從左到右?
??????b=double（b）;?
??????allsamples=[allsamples;?b];??%?allsamples?是一個200?*?10304?矩陣，allsamples?中每一行數據代表一張圖片。
??end?
end?
?
cov_size=size（allsamples1）;?
train_num=cov_size;?
patterns=allsamples;?
for?i=1:cov_size?
??for?j=i:cov_size?
????K（ij）?=（patterns（i:）*patterns（j:）‘+1）^0.7;?
????K（ji）?=?K（ij）;?%核函數是對稱的
??end?
end?
unit?=?ones（cov_size?cov_size）/cov_size;?
K_n?=?K?-?unit*K?-?K*unit?+?unit*K*unit;?%改變核矩陣以實現數據在高維空間去均值
?
[evecsevals]?=?eig（K_n）;?%求去均值后的核矩陣的特征值、特征向量
evals?=?real（diag（evals））;?%取特征值實部
[l1l2]=sort（-1*evals‘）;?
evals=l1*（-1）;?
evecs=evecs（:l2）;?%?取排序后特征值對應的特征向量
for?i=1:cov_size?%因為K為滿秩矩陣
??evecs（:i）?=?evecs（:i）/（sqrt（evals（i）））;?%這樣可以使特征空間中協方差矩陣C的特征向量歸一化
end?
?
?
%??90%的貢獻率?
dsum?=?sum（evals）;?
????dsum_extract?=?0;?
????p?=?0;?
????while（?dsum_extract/dsum?????????p?=?p?+?1;?
????????dsum_extract?=?sum（evals（1:p））;?
????end?
?
?
max_ev=p;?
?
?
unit_train?=?ones（train_numcov_size）/cov_size;?
?
for?i=1:train_num?
??for?j=1:cov_size?
????K_train（ij）?=（patterns（i:）*patterns（j:）‘+1）^0.7;?
??end?
end?
K_train_n?=?K_train?-?unit_train*K?-?K_train*unit?+?unit_train*K*unit;?
features?=?zeros（train_num?max_ev）;?
features?=?K_train_n?*?evecs（:1:max_ev）;?%?訓練樣本在特征空間的投影
?
?
allcoor?=?features;?
?
?
?
accu=0;?
%?測試過程?
for?i=1:40?
????for?j=1:5?%讀入40?x?5?副測試圖像?
????????%a=imread（strcat（‘Data\s‘num2str（i）‘\‘num2str（j）‘.pgm‘））;?
????????a=imread（strcat（‘C:\Documents?and?Settings\Administrator\桌面\人臉庫\ORL人臉庫\s‘num2str（i）‘\‘num2str（j）‘.jpg‘））;
????????b=a（1:10304）;?
????????b=double（b）;?
????????unit_test?=?ones（1cov_size）/cov_size;?
????????K_test?=?zeros（1cov_size）;??
??????for?j=1:cov_size?
????????K_test（j）?=?（b*patterns（j:）‘+1）^0.7;?
??????end?
?
K_test_n?=?K_test?-?unit_test*K?-?K_test*unit?+?unit_test*K*unit;?%?去均值
test_features?=?zeros（1?max_ev）;??
test_features?=?K_test_n?*?evecs（:1:max_ev）;?
?????????
?????????
????????tcoor=test_features;?%計算坐標，是1×p階矩陣?
????????for?k=1:200??
????????????????mdist（k）=norm（tcoor-allcoor（k:））;?%求2范數
????????????end;?
?[distindex2]=sort（mdist）;?
????????class=ceil（?index2（1）/5?）;?
????????if?class==i?
????????????accu=accu+1;?
????????end;?
????end;
end;?
accuracy=accu/200;%輸出識別率
p;

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

?????文件??????10318??1994-04-18?14:17??KPCA人臉識別有效\程序\Data\s1\1.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s1\10.pgm

?????文件??????10318??1994-04-18?14:17??KPCA人臉識別有效\程序\Data\s1\2.pgm

?????文件??????10318??1994-04-18?14:17??KPCA人臉識別有效\程序\Data\s1\3.pgm

?????文件??????10318??1994-04-18?14:17??KPCA人臉識別有效\程序\Data\s1\4.pgm

?????文件??????10318??1994-04-18?14:17??KPCA人臉識別有效\程序\Data\s1\5.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s1\6.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s1\7.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s1\8.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s1\9.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s10\1.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s10\10.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s10\2.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s10\3.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s10\4.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s10\5.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s10\6.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s10\7.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s10\8.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s10\9.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s11\1.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s11\10.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s11\2.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s11\3.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s11\4.pgm

?????文件??????10318??1994-04-18?14:18??KPCA人臉識別有效\程序\Data\s11\5.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s11\6.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s11\7.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s11\8.pgm

?????文件??????10318??1994-04-18?14:07??KPCA人臉識別有效\程序\Data\s11\9.pgm

............此處省略419個文件信息