function?[Train_FeatureTest_Feature]=GetFeature_KFDA（Train_SampleTest_SampleSample）;
%使用核Fisher方法對樣本進行特征提取

sn_ori=size（Sample2）;
Train_Num=5;
[Dim_orisncn]=size（Sample）;
Test_Num=sn-Train_Num;
Train_count=Train_Num*cn;
Test_count=Test_Num*cn;

[Train_Feature2Test_Feature2]=GetProjection（SampleTrain_NumTrain_SampleTest_Sample）;
%訓練樣本，測試樣本分別二維變成三維
[dimn]=size（Train_Feature2）;
T=ones（dimTrain_Numcn）;
for?i=1:cn
????T1=[];
????for?j=1:Train_Num
????????T1=[T1?Train_Feature2（:（i-1）*Train_Num+j）];
????end
????T（::i）=T1;
end
Train_Feature3=T;

T=ones（dimTest_Numcn）;
for?i=1:cn
????T1=[];
????for?j=1:Test_Num
????????T1=[T1?Test_Feature2（:（i-1）*Test_Num+j）];
????end
????T（::i）=T1;
end
Test_Feature3=T;

%在訓練樣本的KPCA特征矩陣上求得散布矩陣Sb、Sw、St
[SbSwSt]=GetScatter（Train_Feature3）;

[VS]=eig（Sw）;
s=diag（S）;
[sindex]=sort（s‘descend‘）;
Rank_Sw=rank（Sw）;
%取Sw的零空間
k=1;
for?i=（Rank_Sw+1）:size（Sw1）
????P1（:k）=V（:index（i））;
????k=k+1;
end
Sb1=P1‘*Sb*P1;
Rank_SbI=rank（Sb1）;
[V1S1]=eig（Sb1）;
s1=diag（S1）;
[s1index1]=sort（s1‘descend‘）;
for?i=1:Rank_SbI
????Z1（:i）=V1（:index1（i））;
end
W1=P1*Z1;
%取P1的正交補空間
for?i=1:Rank_Sw
????P2（:i）=V（:index（i））;
end
Sb2=P2‘*Sb*P2;
St2=P2‘*St*P2;
[V2S2]=eig（Sb2St2）;
s2=diag（S2）;
[s2index2]=sort（s2‘descend‘）;
k=1;
for?j=1:size（s21）
????if?s2（j）>1e-6
????Z2（:k）=V2（:index2（j））;
????k=k+1;
????end
end
W2=P2*Z2;
W=[W1W2];

ProjectSet=W;

%降維，由于?rank（（Sw^-1）*Sb）?<=?rank（Sb）?<=?K-1，所以非零的特征根的個數最多為K-1。所以可以到?K-1?維。
ProjectSet=ProjectSet（:1:sn-1）;

%特征提取
Train_Feature=ProjectSet‘*Train_Feature2;
Test_Feature=ProjectSet‘*Test_Feature2;
%訓練樣本，測試樣本分別二維變成三維
[dimn]=size（Train_Feature）;
T=ones（dimTrain_Numcn）;
for?i=1:cn
????T1=[];
????for?j=1:Train_Num
????????T1=[T1?Train_Feature（:（i-1）*Train_Num+j）];
????end
????T（::i）=T1;
end
Train_Feature=T;

T=ones（dimTest_Numcn）;
for?i=1:cn
????T1=[];
????for?j=1:Test_Num
????????T1=[T1?Test_Feature（:（i-1）*Test_Num+j）];
????end
????T（::i）=T1;
end
Test_Feature=T;

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

?????文件???????2187??2008-06-15?17:45??KFDA\GetFeature_KFDA.m

?????文件????????704??2008-06-15?16:03??KFDA\GetProjection.m

?????文件????????800??2008-06-13?21:38??KFDA\GetScatter.m

?????文件???????2620??2008-06-13?16:49??KFDA\kpca_calc.m

?????文件????????738??2008-06-12?21:34??KFDA\kpca_map.m

?????文件???????1024??2006-03-27?22:52??KFDA\kpca_matrix.m

?????文件????????567??2008-06-15?17:51??KFDA\main.m

?????文件???????1201??2008-06-14?11:27??KFDA\ReadImage.m

?????文件????????758??2008-06-15?11:00??KFDA\Recognition.m

?????目錄??????????0??2008-06-15?17:56??KFDA

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

????????????????10599????????????????????10