|
|
max_tau = 100;
; ], k2 |, M, z0 L/ h& D1 d Part=128;) ~/ _6 Y ^! r
data=xlsread('L.xlsx');: e2 }% ^7 l/ A' r8 ?/ O
loadlibrary('mutual');) d# z- s' v. E0 l. z
[entropy]=mutual(data);# Q+ q( m8 t0 C1 Z# u
for i = 1:length(entropy)-1 $ K. `. ?% a7 P0 }& I
if (entropy(i)<=entropy(i+1))
3 Y2 }0 P4 C) @: t: I% i tau = i-1;
. m! m9 T# O3 u: J0 W$ _ break;
' r3 B5 e, C# g8 _ end5 I7 i2 G- d8 s- a
end
( d* a! N& H" J- T3 p; Q7 @% s + a2 K9 a) }# e4 @9 q U$ N C
plot(0:length(entropy)-1,entropy)
, C! _6 F. U* [! W xlabel('tau');
6 J$ ^9 W0 g- V) o) o E1 _; b( ^+ s title('互信息法求tau'); |
|