|
max_tau = 100;
% o9 ]& @! o* v6 z* W1 B/ ^3 q, C Part=128;
2 I- r! D4 i- A; N0 ` data=xlsread('L.xlsx');
: V% M3 n. I) W: a& T- I# O loadlibrary('mutual');
! G6 t, m$ L# T {4 H [entropy]=mutual(data);1 g- |1 ]9 J- Q. N- l2 N- V" {" H
for i = 1:length(entropy)-1
a, p" \+ @8 A0 c; Z/ | if (entropy(i)<=entropy(i+1))
2 @$ ?/ W& e) _: O+ A& T/ e tau = i-1;
$ N% M- ^+ ~+ R" ?* ?$ o* G: i break;4 V0 ]4 W- S0 ~6 |
end; n4 Q7 h4 Y3 {6 `
end
7 |( C' W/ h8 H( E% K 2 O8 W; L9 |# a, @
plot(0:length(entropy)-1,entropy)! c( L; y2 ]5 l9 F% |' O: m& e
xlabel('tau');! B/ t4 W4 q7 G* |6 {4 y8 [" K# \
title('互信息法求tau'); |
|