clc;
clear;
cls1_data=[2.93?6.634;?2.53?7.79;?3.57?5.65;3.16?5.47];%第一個類的訓(xùn)練集
cls2_data=[2.58?4.44;?2.16?6.22;?3.27?3.52];%第二個類的訓(xùn)練集
%求期望
E_cls1=mean（cls1_data）;%第一類數(shù)據(jù)的期望矩陣
E_cls2=mean（cls2_data）;%第二類數(shù)據(jù)的期望矩陣
E_all=mean（[cls1_data;cls2_data]）;%所有訓(xùn)練集的期望矩陣
%%%%%%%%%%%%%%%%%%%%分類前畫圖%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%計算類間離散度矩陣：
x1=E_all-E_cls1;
x2=E_all-E_cls2;
Sb=4*x1‘*x1/7+3*x2‘*x2/7;%%%%%%%%%為什么不用式（3）？？？
%計算類內(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: