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:）