|
max_tau = 100;) F) u" Y- u# _4 }' |4 c: A
Part=128;" O- R# O, f5 @ e
data=xlsread('L.xlsx');
$ a; m N+ @2 A/ E" [6 g loadlibrary('mutual');) `% ]; M% U: t9 v. n2 y1 c3 o
[entropy]=mutual(data);0 F# J& B5 G8 l5 p' g% I$ v3 Q4 h
for i = 1:length(entropy)-1 : c/ Q4 { D, N- b9 J ^' o4 T
if (entropy(i)<=entropy(i+1))
3 ^+ g" L+ G; O' H tau = i-1; 5 V/ J9 B4 B% A4 k% {$ ~
break;
$ O5 I3 M5 n- ~4 O) F end* [8 Y& a/ e1 E, T0 N* _( Q- F2 s
end3 g0 x2 m# |% K7 x
# x9 q% Z! k3 G% p0 v0 i
plot(0:length(entropy)-1,entropy)
8 u4 o0 m+ E0 s% e* ]: l! m xlabel('tau');
& P& r0 d0 g6 K+ n title('互信息法求tau'); |
|