%Caculate?MI?of?a?and?b?in?the?region?of?the?size（a）

function?[mi]?=?EMI（ab）
%
[MN]?=?size（a）;
%?change?a?b?to?range?1:n
if?min（a（:））<1
????a?=?a-min（a（:））+1;
end
if?min（b（:））<1
????b?=?b-min（b（:））+1;
end
????

%?ncolora?=?max（a（:））-min（a（:））+1;
%?ncolorb?=?max（b（:））-min（b（:））+1;
ncolora?=?max（a（:））;
ncolorb?=?max（b（:））;

hab?=?zeros（ncolorancolorb）;
ha?=?zeros（1ncolora）;
hb?=?zeros（1ncolorb）;

%count?begins
[MN]=size（a）;
for?i=1:M
????for?j=1:N
???????indexx?=??a（ij）;
???????indexy?=??b（ij）;

???????hab（indexxindexy）?=?hab（indexxindexy）+1;
???????ha（indexx）?=?ha（indexx）+1;
???????hb（indexy）?=?hb（indexy）+1;

???end
end

%cal?hahbhab
hsum?=?sum（sum（hab））;
index?=?find（hab~=0）;
p?=?hab/hsum;
Hab?=?sum（sum（-p（index）.*log（p（index））））;

hsum?=?sum（sum（ha））;
index?=?find（ha~=0）;
p?=?ha/hsum;
Ha?=?sum（sum（-p（index）.*log（p（index））））;

hsum?=?sum（sum（hb））;
index?=?find（hb~=0）;
p?=?hb/hsum;
Hb?=?sum（sum（-p（index）.*log（p（index））））;

mi?=?Ha+Hb-Hab;
%?cal?EMI
%?feature?efficiency?coefficient
n?=2/5;
if?Hab?~=?0
????e?=?（mi^n）?/?（?Hab^（1-n））;
else%?not?usual?happen?just?for?test
????e?=?-1;%（mi^n）?/?（?Hab^（1-n）+1e-10）;
end
mi?=?e;
%?%?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????653??2005-10-19?20:10??peizhun\myMI.m

?????文件????????697??2005-10-19?20:09??peizhun\FPMI.m

?????文件???????3338??2005-10-19?20:10??peizhun\GMI.m

?????文件????????228??2005-10-19?20:10??peizhun\gobad.m

?????文件???????1404??2005-10-19?20:10??peizhun\HiMI.m

?????文件???????1034??2005-10-19?20:10??peizhun\MI.m

?????文件???????1218??2005-10-19?20:09??peizhun\EMI.m

?????文件????????439??2005-10-19?20:10??peizhun\oneDimSearch.m

?????文件???????2528??2005-10-19?20:10??peizhun\POWELL.m

?????文件??????35888??2005-10-19?20:10??peizhun\proj.fig

?????文件??????12544??2005-10-19?20:10??peizhun\proj.m

?????文件???????5902??2005-10-19?20:10??peizhun\PSO.m

?????文件????????390??2005-10-19?20:10??peizhun\restore.m

?????文件???????3719??2005-11-27?20:04??peizhun\RMI.m

?????文件?????????67??2005-10-19?20:10??peizhun\test.m

?????文件????????238??2005-10-19?20:10??peizhun\test2.m

?????文件???????7648??2005-10-19?20:10??peizhun\testListBoxlj.fig

?????文件???????7936??2005-11-27?20:02??peizhun\testListBoxlj.m

?????文件????????326??2005-10-19?20:10??peizhun\說明.txt

?????目錄??????????0??2005-11-27?20:08??peizhun

-----------?---------??----------?-----??----

????????????????86197????????????????????20