%%%%??讀入圖片???%%%%%%%%%%%%%%%%%%%%%%%%
%?clc;clear;close;
F1=imread（‘SE3_t613.bmp‘）;
G1=imread（‘SE3_t613.bmp‘）;
%%%%%??計算圖像F1和G1的互信息???%%%%%%%%%%%%%%%%%%%%
F1_1=rgb2gray（F1）;
G1_1=rgb2gray（G1）;
Pf_1=imhist（F1_1）./numel（F1_1）;
Pg_1=imhist（G1_1）./numel（G1_1）;
Pfg_1=zeros（256256）;
for?i=1:512
????for?j=1:512
????????Pfg_1（F1_1（ij）+1G1_1（ij）+1）=Pfg_1（F1_1（ij）+1G1_1（ij）+1）+1;
????end
end
num=sum（Pfg_1（:））;
Pfg_1=Pfg_1/num;
W1_1=zeros（1256）;
W2_1=zeros（1256）;
W12_1=zeros（1256）;
for?m=1:256
????if?Pf_1（m）~=0
???????W1_1（m）=Pf_1（m）.*log10（Pf_1（m））;
????end
end
for?m=1:256
????if?Pg_1（m）~=0
???????W2_1（m）=Pg_1（m）.*log10（Pg_1（m））;
????end
end
for?m=1:256
????if?Pfg_1（m）~=0
???????W12_1（m）=Pfg_1（m）.*log10（Pfg_1（m））;
????end
end
N1_1=-sum（W1_1（:））;
N2_1=-sum（W2_1（:））;
N12_1=-sum（W12_1（:））;