|
max_tau = 100; }; s; a5 A0 C3 G3 d2 W- b
Part=128;" R2 J7 d7 s; X, h/ a
data=xlsread('L.xlsx');
" G" F/ i" G) j( ~8 x, x/ ~ loadlibrary('mutual');% a6 k* L# D( z* T
[entropy]=mutual(data);5 r# V! x3 x. E3 n8 v
for i = 1:length(entropy)-1
/ X+ N" ~- w# Q; i3 e u: o/ Z- g8 Y if (entropy(i)<=entropy(i+1))
d$ X ~ n' k& M& C1 @& w: ? tau = i-1; ) f5 n* q! O& C7 m" P U
break;/ Y& I, |' e5 z' T( `( ]
end3 A5 C& K7 J2 N3 V1 O
end
; l" \* \- _6 H. Y3 \" O0 U * E' n9 v+ f# [3 c4 k& Q! o0 @- B1 }
plot(0:length(entropy)-1,entropy)
' l+ `9 `; F5 y6 x# A xlabel('tau');/ G4 o4 T" [/ V& |- D8 g& y1 x# A. q
title('互信息法求tau'); |
|