监督分类最大似然法MATLAB.zip

最大似然法是经常使用的监督分类方法之一，适用于遥感数据分类，压缩包中的代码可以直接在MATLAB上运行，不提供数据集。

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）; %协方差
[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‘‘分类结果另存为‘）;
fid=fopen（strcat（filepathfname）‘wt‘）;
for j=1:mu
    fprintf（fid‘%d\n‘class（j1））;   
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1803  2019-06-20 22:59  监督分类最大似然法（MATLAB）\最大似然.m
     目录           0  2020-06-16 10:13  监督分类最大似然法（MATLAB）\