U=load（‘D:\MATLAB\bin\data.txt‘）;
u=U（:1:7）;%光譜信息
[munu]=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）;?%協(xié)方差
[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（i1: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]‘;
?????[xy]=max（D）;
????class_now（i1）=y;?
end

%uo=U（:8）;
%c=0
%for?i=1:mu
????%if（uo（i1）==class_now（i1））
???????%?c=c+1;
????%end
%end
class=class_now;

[fnamefilepath]=uiputfile（‘*.txt‘‘分類結(jié)果另存為‘）;
fid=fopen（strcat（filepathfname）‘wt‘）;
for?j=1:mu
????fprintf（fid‘%d\n‘class（j1））;???
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1803??2019-06-20?22:59??監(jiān)督分類最大似然法（MATLAB）\最大似然.m
?????目錄???????????0??2020-06-16?10:13??監(jiān)督分類最大似然法（MATLAB）\