|
|
max_tau = 100;/ x0 v6 {9 M1 N k0 H2 ^
Part=128;2 i1 h7 r$ U. F& g# R& G
data=xlsread('L.xlsx');
, V! t: k( \( P# M- L ?) T loadlibrary('mutual');
+ I1 B$ F1 s: o/ ^. | [entropy]=mutual(data);
9 H# { [% r' w$ y for i = 1:length(entropy)-1 ( K. ~$ @/ q" s
if (entropy(i)<=entropy(i+1))$ W! T) N; |! _7 |7 J4 E, X
tau = i-1;
8 c5 ~* n g6 T2 O" P. I) | break;! @7 A+ f$ O- v; d! o9 `
end
' @( D1 d! ~* j9 O8 Z G end
* v+ s+ g4 f% n7 D7 |+ q
5 G* c2 ~8 V1 H" n, P* G plot(0:length(entropy)-1,entropy)
* s4 s7 ]3 |) \" O7 [. t xlabel('tau');
8 }8 Q! J& N5 }! f title('互信息法求tau'); |
|