% filename='wdbc.lz';
filename='ionosphere.lz';
fid=fopen(filename,'rb');
nclass=fread(fid,1,'int');
ndimension=fread(fid,1,'int');
sample=fread(fid,256,'short');
n=sum(sample);
n1=sample(1);
n2=sample(2);
data1=fread(fid,[ndimension,n],'float');
fclose(fid);
meaneachclass(:,1)=mean(data1([1:225]),2);
meaneachclass(:,2)=mean(data1([226:351]),2);
meanallsample=mean(data1,2);
%*****************************求取类间类内离散度矩阵********************************************%
sb=0;
for i=1:2
sb=sb+sample(i)*(meaneachclass(:,i)-meanallsample)*(meaneachclass(:,i)-meanallsample)';
end
sw=0;
for i=1:2
if 1==i
for j=1:225
sw=sw+(meaneachclass(:,i)*data1(:,j))*(meaneachclass(:,i)*data1(:,j))';
end
elseif 2==i
for j=226:351
sw=sw+(meaneachclass(:,i)*data1(:,j))*(meaneachclass(:,i)*data1(:,j))';
end
end
end
cov=inv(sw)*sb;
[x y z]=pcacov(cov*cov/34);
y=x'*data1;
plot(y(1,1:n1),y(2,1:n1),'*',y(1,n1+1:n),y(2,n1+1:n),'*')
LDA.rar （1个子文件） 
LDA.m 1KB
寒泊
- 粉丝: 86
- 资源: 1万+
评论0