K-L变换算法matlab实现

clear;
load  MALE.txt
load  FEMALE.txt
data1=[MALE];
data2=[FEMALE];
data=[MALE;FEMALE];
      k1=mean（data1）;          %第一类均值
      k2=mean（data2）;    %第二类均值
      k=mean（data）;

R1=cov（data1（:1）data1（:2））;           %协方差矩阵
R2=cov（data2（:1）data2（:2））; 
R=cov（data（:1）data（:2））;
SW=0.5*（R1+R2）;
[VD] = eig（SW）; 
%[V1D1] = eig（R1）; 
%[V2D2] = eig（R2）; 
sb=0.5*（（k1-k）‘*（k1-k）+（k2-k）‘*（k2-k））
J1（:1）=V（:1）‘*sb*V（:1）/D（11）;
J2（:2）=V（:2）‘*sb*V（:2）/D（22）;
i=1;
if J1（:1）    i=2;
end
                    %namida=D（ii）
w=V（i:）