matlab实现LDA

clc;
clear;
cls1_data=[2.93 6.634; 2.53 7.79; 3.57 5.65;3.16 5.47];%第一个类的训练集
cls2_data=[2.58 4.44; 2.16 6.22; 3.27 3.52];%第二个类的训练集
%求期望
E_cls1=mean（cls1_data）;%第一类数据的期望矩阵
E_cls2=mean（cls2_data）;%第二类数据的期望矩阵
E_all=mean（[cls1_data;cls2_data]）;%所有训练集的期望矩阵
%%%%%%%%%%%%%%%%%%%%分类前画图%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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）？？？
%计算类内离散度矩阵
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: