U=load('D:\MATLAB\bin\data.txt');
u=U(:,1:7);%光谱信息
[mu,nu]=size(u);
A=load('C:\Users\86159\Desktop\MY\1d.txt'); %墙
B=load('C:\Users\86159\Desktop\MY\2d.txt'); %陶瓷
C=load('C:\Users\86159\Desktop\MY\\3d.txt'); %纸板
D=load('C:\Users\86159\Desktop\MY\4d.txt'); %泡沫塑料
E=load('C:\Users\86159\Desktop\MY\5d.txt'); %仙人掌
F=load('C:\Users\86159\Desktop\MY\6d.txt'); %植物
G=load('C:\Users\86159\Desktop\MY\7d.txt'); %枯萎植物
MA=mean(A);MB=mean(B);MC=mean(C);MD=mean(D);ME=mean(E);MF=mean(F);MG=mean(G);%均值
EA=cov(A); EB=cov(B); EC=cov(C); ED=cov(D);EE=cov(E);EF=cov(F);EG=cov(G); %协方差
[ma,~]=size(A);[mb,~]=size(B);[mc,~]=size(C);[md,~]=size(D);[me,~]=size(E);[mf,~]=size(F);[mg,~]=size(G);
PA=ma/(ma+mb+mc+md+me+mf+mg);
PB=mb/(ma+mb+mc+md+me+mf+mg);
PC=mc/(ma+mb+mc+md+me+mf+mg);
PD=md/(ma+mb+mc+md+me+mf+mg);
PE=me/(ma+mb+mc+md+me+mf+mg);
PF=mf/(ma+mb+mc+md+me+mf+mg);
PG=mg/(ma+mb+mc+md+me+mf+mg);
%判别
for i=1:mu
%选择第i个向量U
%计算Da Db Dc Dd
X=U(i,1:7);
Da=-0.5*((X-MA)/(EA)*(X-MA)')-0.5*log(det(EA))+log(PA);
Db=-0.5*((X-MB)/(EB)*(X-MB)')-0.5*log(det(EB))+log(PB);
Dc=-0.5*((X-MC)/(EC)*(X-MC)')-0.5*log(det(EC))+log(PC);
Dd=-0.5*((X-MD)/(ED)*(X-MD)')-0.5*log(det(ED))+log(PD);
De=-0.5*((X-ME)/(EE)*(X-ME)')-0.5*log(det(EE))+log(PE);
Df=-0.5*((X-MF)/(EF)*(X-MF)')-0.5*log(det(EF))+log(PF);
Dg=-0.5*((X-MG)/(EG)*(X-MG)')-0.5*log(det(EG))+log(PG);
D=[Da Db Dc Dd De Df Dg]';
[x,y]=max(D);
class_now(i,1)=y;
end
%uo=U(:,8);
%c=0
%for i=1:mu
%if(uo(i,1)==class_now(i,1))
% c=c+1;
%end
%end
class=class_now;
[fname,filepath]=uiputfile('*.txt','分类结果另存为');
fid=fopen(strcat(filepath,fname),'wt');
for j=1:mu
fprintf(fid,'%d\n',class(j,1));
end
