X1=[x1;x2;x3]; X2=X1'; k1=mean(x1);%第一类均值 k2=mean(x2);%第二类均值 k3=mean(x3);%第三类均值 k=mean(X1);%总均值向量 R1=cov(x1(:,:,:));%x1类协方差矩阵 R2=cov(x2(:,:,:));%x2类协方差矩阵 R3=cov(x3(:,:,:));%x3类协方差矩阵 %总的类内离散度矩阵 SW=(200/600)*R1+(190/600)*R2+(210/600)*R3; [V,D]=eig(SW);%求取SW的本征值和本征向量 %类间离散度矩阵 sb=(200/600)*(