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
???

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????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