|
max_tau = 100;
( T+ _& ~6 _5 \: a) z/ L Part=128;8 B" x" S& v; L+ I* g0 h" }
data=xlsread('L.xlsx');3 Q' m, a$ J& P: H, _* S
loadlibrary('mutual');
6 B! s- W% \. o1 J* @6 G, j [entropy]=mutual(data);& |% H0 ]4 _ G7 z" U6 C1 f L
for i = 1:length(entropy)-1
$ }2 I, S# |1 `/ h3 x+ j! I! G if (entropy(i)<=entropy(i+1))
" [6 h. b8 q& W+ W% O) N( L9 H5 K tau = i-1;
D' N' k n+ Q3 _8 C' {1 B" ` break;
! _/ m! R$ g T6 |0 x end
3 L2 j$ S" s1 m- R \+ K- h end
( v/ U" \, c( T2 Z0 N
. \. ~) |3 `# n plot(0:length(entropy)-1,entropy)
; }! R8 p& K/ N xlabel('tau');
( t2 s* i. c! V4 a7 O2 |$ e3 B) C" E title('互信息法求tau'); |
|