|
max_tau = 100;' k A8 t$ M1 s- b& k
Part=128;
8 A- y& N9 q0 q8 Z5 @ data=xlsread('L.xlsx');6 T' ^3 S( @1 ?/ M# X
loadlibrary('mutual');, v$ ?5 @, ?" E5 A/ y4 r( H
[entropy]=mutual(data); l& f3 g# R/ r- N! x: Y* a4 v
for i = 1:length(entropy)-1 * F1 U6 g+ h' L% U1 r
if (entropy(i)<=entropy(i+1))
0 a" z1 C/ e& w0 ~& W; ]& Y* a tau = i-1; + X1 \ p7 L4 B2 S: d
break;9 ~; {! q. T' F6 R7 _* A- y
end* E$ j& K* [6 t, b" ^7 F7 M
end
6 l2 ~) n! a- y6 D+ b 0 }8 s& B, [8 [# i `
plot(0:length(entropy)-1,entropy), l- h& U1 `: R M% B3 B- `% F- y
xlabel('tau');- E9 ^0 N! T) S0 T. p$ {
title('互信息法求tau'); |
|