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是訓練樣本的個數
train_kpca=[K_n?*?evectors（:1:index（1））];
unit_test?=?ones（test_numcov_size）/cov_size;%cov_size是訓練樣本的個數
K_test?=?zeros（test_numcov_size）;?%test_num是測試樣本的個數，?cov_size是訓練樣本的個數
for?i=1:test_num?%test_num是測試樣本的個數
for?j=1:cov_size%cov_size是訓練樣本的個數
K_test（ij）?=?exp（-norm（test_patterns（i:）-patterns（j:））^2/rbf_var）;
end
end
K_test_n?=?K_test?-?unit_test*K?-?K_test*unit?+?unit_test*K*unit;
test_kpca?=?zeros（test_num?index（1））;%test_num是測試樣本的個數
test_kpca?=?[K_test_n?*?evectors（:1:index（1））];
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1651??2016-12-22?19:10??kpca?lda?mds降維\KPCA.m
?????文件????????1454??2016-12-28?15:45??kpca?lda?mds降維\LDA.m
?????文件?????????800??2016-12-28?16:38??kpca?lda?mds降維\MDS.m
?????文件?????????898??2017-03-24?16:00??kpca?lda?mds降維\Main.m
?????文件??????376589??2016-12-16?17:53??kpca?lda?mds降維\ORL.mat
?????文件?????????833??2016-12-28?15:49??kpca?lda?mds降維\PCA.m
?????文件?????????230??2017-03-24?16:00??kpca?lda?mds降維\tiqu.m
?????目錄???????????0??2017-03-24?16:01??kpca?lda?mds降維\