一维最大熵阈值分割

clear;
I=imread（‘D:\hdtx\000000014.bmp‘）;
count=imhist（I）;
% imhist（I）
[mn]=size（I）;
N=m*n;
       %L=87;     %重新计算t1=21  
       %L=91;    %重新计算t2=80
       L=256;   %重新计算t3=175
p=count/N;
s1=cumsum（p）;
s2=1-s1;
A=zeros（LL）;
B=zeros（LL）;
AA=zeros（LL）;
BB=zeros（LL）;
for i=1:L       %i的取值在计算t1时 i=1:L
                 %i的取值在计算t2时 i=t1:T3    
                 %i的取值在计算t3时 i=t2:256
    for j=1:i
        if （s1（i）>0）
            A（ij）=p（j）/s1（i）;    %计算A灰度级分布
        else
            A（ij）=0;
        end
        for k=i+1:L 
            if （s2（i）>0）
                B（ik）=p（k）/s2（i）;  %计算B灰度级分布
            else
                B（ik）=0;
            end
        end
        if （A（ij）>0）
            AA1（ij）=log10（A（ij））;    %每个分布的熵
        e