求sena.IMG的一阶熵，二阶熵和差分熵的matlab程序

求sena.IMG的一阶熵，二阶熵和差分熵的matlab程序

fig=fopen（‘SENA.img‘）;
A=fread（fig）;

[MN]=size（A）;
temp=zeros（1256）;
for m=1:M;
    for n=1:N;
    
        if A（mn）==0;
   i=0;
  else
   i=A（mn）;
  end
        temp（i）=temp（i）+1;
    end
end
temp=temp./（M*N）;
firstorder=0;

for i=1:length（temp）
    if temp（i）==0;
        firstorder=firstorder;
    else
        firstorder=firstorder-temp（i）*log2（temp（i））;
    end
end
firstorder
   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        417  2010-12-01 10:56  一阶熵 二阶熵 差分熵matlab程序\a.m

     文件        610  2010-12-01 10:52  一阶熵 二阶熵 差分熵matlab程序\b.m

     文件        447  2010-12-01 10:56  一阶熵 二阶熵 差分熵matlab程序\c.m

     文件      65536  2010-11-27 14:21  一阶熵 二阶熵 差分熵matlab程序\SENA.IMG

     目录          0  2010-12-01 15:39  一阶熵 二阶熵 差分熵matlab程序

----------- ---------  ---------- -----  ----

                67010                    5