function?[evectorstrain_kpcatest_kpca]?=KPCA（traintest）
rbf_var=10000;
threshold?=?90;
patterns=zscore（train）;?%訓練數據標準化
test_patterns=zscore（test）;?%測試數據標準化
train_num=size（patterns1）;?%train_num是訓練樣本的個數
test_num=size（test_patterns1）;%test_num是測試樣本的個數
cov_size?=?train_num;?%cov_size是訓練樣本的個數
for?i=1:cov_size
for?j=i:cov_size
K（ij）?=?exp（-norm（patterns（i:）-patterns（j:））^2/rbf_var）;?%核函數?rbf_var?？？
K（ji）?=?K（ij）;
end
end
unit?=?ones（cov_size?cov_size）/cov_size;%cov_size是樣本的個數
K_n?=?K?-?unit*K?-?K*unit?+?unit*K*unit;%?中心化核矩陣
[evectors_1evaltures_1]?=?eig（K_n/cov_size）;
[xindex]=sort（real（diag（evaltures_1）））;?%sort每行按從小到大排序，x為排序后結果，index為索引
evals=flipud（x）?;%?flipud函數實現矩陣的上下翻轉
index=flipud（index）;
evectors=evectors_1（:index）;
%?for?i=1:cov_size
%?evecs（:i）?=?evectors（:i）/（sqrt（evectors（:i）））;
%?end
train_eigval?=?100*cumsum（evals）./sum（evals）;
index?=?find（train_eigval?>threshold）;

train_kpca?=?zeros（train_num?index（1））;?%train_num是訓練樣本的個數
tra

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1651??2016-12-22?19:10??kpca+elm\KPCA.m
?????文件?????????351??2016-12-22?16:15??kpca+elm\LEM_test.m
?????文件?????????434??2016-12-22?16:05??kpca+elm\LEM_train.m
?????文件?????????706??2017-03-24?15:49??kpca+elm\Main.m
?????文件??????376589??2016-12-16?17:53??kpca+elm\ORL.mat
?????文件?????????227??2017-03-24?15:49??kpca+elm\tiqu.m
?????目錄???????????0??2017-03-24?15:48??kpca+elm\