|
|
max_tau = 100;' R S4 ~8 p! N6 i+ o, I6 \* P
Part=128;
/ t% v9 u8 i. {7 x$ A% t' r data=xlsread('L.xlsx');: B: ~' g% ^& ]. O s
loadlibrary('mutual');
! [; ]& }. p- V$ Y- U$ R9 ~ [entropy]=mutual(data);
@& t" Z) v7 h& h& N2 R X for i = 1:length(entropy)-1
- A, R: R% t7 k0 s if (entropy(i)<=entropy(i+1))
- G0 {, T0 k5 Z+ u6 X' G tau = i-1; ! _/ Q5 [9 s3 q9 a
break;, D) D0 z2 `3 B" ~, b0 H& F+ \* `) D
end. E' t/ g7 {4 j w' ^, u
end
8 U( V! c% }+ G8 ^! E+ B
4 N! K9 \) @7 g2 j plot(0:length(entropy)-1,entropy); W. w( G( ?6 }+ F4 U* Z
xlabel('tau');
/ r* ~# H) u$ E- |% Y title('互信息法求tau'); |
|