|
max_tau = 100;% \! x) K& n! y6 O Q9 {- N% J
Part=128;
+ ], V" [6 Q( P, C: u9 A+ c6 u1 \ data=xlsread('L.xlsx');
. b& O2 ?; v) @6 ~1 Y" c# m loadlibrary('mutual');
7 e; q& X( d3 z" @/ m& G/ K [entropy]=mutual(data);4 R1 B8 l n- F8 g7 J0 T
for i = 1:length(entropy)-1
" x0 N2 ]$ }" i) H/ D! w0 O/ L if (entropy(i)<=entropy(i+1))5 \. \% h, X" C9 E q) c
tau = i-1;
! A; z9 o* Y1 U ? break; c9 S, e8 {. Y2 O) j
end7 d. E$ b L" o8 }! t8 _
end8 U% d! h) O9 Q6 m( ] `
$ N0 w5 {9 c, D0 N
plot(0:length(entropy)-1,entropy)
% n' I2 f' E& q- j4 k# u5 @* u xlabel('tau');
& N! U, w/ z1 G9 M" U0 r& E6 C title('互信息法求tau'); |
|