|
max_tau = 100;' B; ^5 h# t) Q, }/ w4 H- {: D1 {
Part=128;; C7 G. a5 U! w) h3 x+ m! r
data=xlsread('L.xlsx');- z9 T( T; J; g: T
loadlibrary('mutual');7 N! M2 g; z' ]4 U4 t9 O
[entropy]=mutual(data);
* n# ^! V8 G9 m! }$ D for i = 1:length(entropy)-1 8 V: ?' }6 s. A/ X- C% G, }3 {
if (entropy(i)<=entropy(i+1))& ]& F. n% B; i- N4 E
tau = i-1;
' C1 s# B4 o1 `! v } break;
# n4 z ], _ I# \ end
4 a7 w+ D0 ]7 _8 w* i: b end4 T Q, e* M/ y; z
& j& P, d6 M/ W/ V plot(0:length(entropy)-1,entropy)
4 ?$ Y) O8 t _" v. l xlabel('tau');: a: p. I1 J7 b$ @0 r) l9 g m
title('互信息法求tau'); |
|