clc;
clear;
cls1_data=[2.93?6.634;?2.53?7.79;?3.57?5.65;3.16?5.47];%第一個(gè)類(lèi)的訓(xùn)練集
cls2_data=[2.58?4.44;?2.16?6.22;?3.27?3.52];%第二個(gè)類(lèi)的訓(xùn)練集
%求期望
E_cls1=mean（cls1_data）;%第一類(lèi)數(shù)據(jù)的期望矩陣
E_cls2=mean（cls2_data）;%第二類(lèi)數(shù)據(jù)的期望矩陣
E_all=mean（[cls1_data;cls2_data]）;%所有訓(xùn)練集的期望矩陣
%%%%%%%%%%%%%%%%%%%%分類(lèi)前畫(huà)圖%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?i=1:4
?????plot（cls1_data（i1）cls1_data（i2）‘.r‘）;
?????hold?on;
end;
plot（E_cls1（1）E_cls1（2）‘^r‘）;
hold?on;
for?i=1:3
?????plot（cls2_data（i1）cls2_data（i2）‘*b‘）;
?????hold?on;
end;
plot（E_cls2（1）E_cls2（2）‘^b‘）;
plot（E_all（1）E_all（2）‘vc‘）;
hold?on;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%計(jì)算類(lèi)間離散度矩陣：
x1=E_all-E_cls1;
x2=E_all-E_cls2;
Sb=4*x1‘*x1/7+3*x2‘*x2/7;%%%%%%%%%為什么不用式（3）？？？
%計(jì)算類(lèi)內(nèi)離散度矩陣
y1=0;
for?i=1:4
????y1=y1+（cls1_data（i:）-E_cls1）‘*（cls1_data（i:）-E_cls1）;
end;
y2=0;
for?i=1:3
????y2=y2+（cls2_data（i:）-E_cls2）‘*（cls2_data（i:）-E_cls2）;
end;
Sw=4*y1/7+3*y2/7;%%%%%%%%%為什么不用式（3）？？？
%求最大特征值和特征向量
[VL]=eig（inv（Sw）*Sb）;

[ab]=max（max（L））;
newspace=V（:b）;%最大特征值所對(duì)應(yīng)的特征向量
new_cls1_data=cls1_data*newspace;%訓(xùn)練后的數(shù)據(jù)集
new_cls2_data=cls2_data*newspace;%訓(xùn)練后的數(shù)據(jù)集
%%%%%%%%%%%%%%%%%%畫(huà)圖代碼%%%%%%%%%%%%%%%%%

k=newspace（2）/newspace（1）;
plot（[06][06*k]‘-c‘）;%畫(huà)出最大特征值對(duì)應(yīng)的特征向量，即樣本所組成的線性空間所投影的子空間
axis（[0?7?-2?9]）;
hold?on;
%畫(huà)出樣本投影到子空間點(diǎn)
for?i=1:4
????temp=cls1_data（i:）;
????newx=（temp（1）+k*temp（2））/（k*k+1）;
????newy=k*newx;
????plot（newxnewy‘*r‘）;
end;

for?i=1:3
????temp=cls2_data（i:）;
????newx=（temp（1）+k*temp（2））/（k*k+1）;
????newy=k*newx;
????plot（newxnewy‘ob‘）;
end;
%預(yù)測(cè)
prediction=[4.81?3.46];
result=prediction*newspace;
temp=new_cls1_data-[result;result;result;result];
temp2=new_cls2_data-[result;result;result];
if（min（abs（temp））>min（abs（temp2）））
????output=‘該樣本屬于不合格產(chǎn)品‘
else
????output=‘該樣本屬于合格產(chǎn)品‘
end;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????2044??2010-08-07?17:15??ＭＡＴＬＡＢ實(shí)現(xiàn)的ＬＤＡ算法及結(jié)果截圖\LDA.m

?????文件???????5335??2010-08-07?17:17??ＭＡＴＬＡＢ實(shí)現(xiàn)的ＬＤＡ算法及結(jié)果截圖\LDA_result.fig

?????目錄??????????0??2010-08-07?17:21??ＭＡＴＬＡＢ實(shí)現(xiàn)的ＬＤＡ算法及結(jié)果截圖

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

?????????????????7379????????????????????3