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