|
max_tau = 100;
h' [! Q3 y7 q0 P, t Part=128;
4 X. I) E0 C6 Z data=xlsread('L.xlsx');4 V, n* e5 m+ e1 U
loadlibrary('mutual');, F$ A8 k' W+ i' h' c$ p
[entropy]=mutual(data);
1 A5 f/ X+ U* o- Z1 g- ~ A+ ] for i = 1:length(entropy)-1 ' g1 ]1 b2 V$ W: |
if (entropy(i)<=entropy(i+1))) h4 B/ |( ~1 q8 e! A
tau = i-1; 2 f8 W. ~" U: b- D- C1 v0 T
break;
! q; v+ H4 r" W* l2 Q end
& s0 y3 V+ j1 `" e8 z& m end
) \( j) c, ?% q z0 S4 \5 W : [: N6 N/ ]2 ]5 f+ `
plot(0:length(entropy)-1,entropy)
$ C5 z9 P- s) g8 j xlabel('tau');
1 j# E+ n' j( _5 T2 N' r title('互信息法求tau'); |
|