ksw.m

来自「遗传算法就的开始株型紧凑可执行：早晨从看见着在科学考察偶支持」· M 代码 · 共 35 行

function y=ksw(t,mingrayvalue,maxgrayvalue,hist1,HT)
   
    %计算最佳直方图熵（KSW熵）
    
    Pt=0;
    for i=mingrayvalue:t
        Pt=Pt+hist1(i+1);
    end
    
    Ht=0;
    for i=mingrayvalue:t
        if hist1(i+1)==0
           temp=0;
        else
           temp=hist1(i+1)*log(1/hist1(i+1));
        end
        Ht=Ht+temp;
    end
    
    if Pt==0 | Pt==1           %   or(Pt==0,Pt==1)
        temp1=0;
    else 
        temp1=log(Pt*(1-Pt))+Ht/Pt+(HT-Ht)/(1-Pt);
    end
    
    if temp1 < 0
        H=0;
    else
        H=temp1;
    end
        
    
    y=H;