%用KPCA進行人臉識別的Matlab代碼?

clear?
clc?
allsamples=[];%所有訓練圖像?
for?i=1:40?
????for?j=6:10
??????a=imread（strcat（‘Data\s‘num2str（i）‘\‘num2str（j）‘.pgm‘））;?
??????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?
??evecs（:i）?=?evecs（:i）/（sqrt（evals（i）））;?
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‘））;?
????????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:））;?
????????????end;?
?[distindex2]=sort（mdist）;?
????????class=ceil（?index2（1）/5?）;?
????????if?class==i?
????????????accu=accu+1;?
????????end;?
????end;
end;?
accuracy=accu/200?%輸出識別率
p

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

?????文件???????2345??2009-06-18?10:33??程序\kpca.asv

?????文件???????2352??2009-06-23?11:10??程序\kpca.m

?????文件??????86528??2009-06-25?09:53??程序\程序說明.doc

?????文件??????10318??1994-04-18?14:18??程序\Data\s9\1.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s9\10.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s9\2.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s9\3.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s9\4.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s9\5.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s9\6.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s9\7.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s9\8.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s9\9.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s8\1.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s8\10.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s8\2.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s8\3.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s8\4.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s8\5.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s8\6.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s8\7.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s8\8.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s8\9.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s7\1.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s7\10.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s7\2.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s7\3.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s7\4.pgm

?????文件??????10318??1994-04-18?14:18??程序\Data\s7\5.pgm

?????文件??????10318??1994-04-18?14:07??程序\Data\s7\6.pgm

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