|
|
max_tau = 100;
) H! G3 ~& m+ m Part=128; I: u- r$ q' y
data=xlsread('L.xlsx');
7 Z* A2 k: O1 u% [ loadlibrary('mutual');+ O( n/ d% |8 r
[entropy]=mutual(data);6 ?+ V/ K% p( z4 y& Z3 w# V4 n2 l; s2 }
for i = 1:length(entropy)-1 + f F) q8 K- N1 U
if (entropy(i)<=entropy(i+1))
$ w" n" d+ N+ @* W tau = i-1; 9 x3 B+ |. B [% M3 H/ W6 Q. I+ t) }
break;6 J5 _7 V2 z+ n( ~6 w" v5 a
end
6 U+ y# O, p7 b" Z3 O+ n end
2 f# _0 \! M C1 ~& W0 ^, w , J5 m% C/ `4 e7 c
plot(0:length(entropy)-1,entropy)
! u4 O' [/ y" N( ^% ?+ f xlabel('tau');
" v; t( L: g- k$ `1 Z1 I title('互信息法求tau'); |
|