|
|
max_tau = 100;
, h. e7 `. ~3 l! C; S Part=128;2 _3 Z2 i2 l- x- n& W
data=xlsread('L.xlsx');+ S9 ?7 E4 k" R h8 y k6 M
loadlibrary('mutual');% G" c* I5 R* R. l8 Y0 H
[entropy]=mutual(data);% s9 L9 F$ S, U* G
for i = 1:length(entropy)-1
% u2 d$ a4 ~7 Z3 n( D if (entropy(i)<=entropy(i+1))
f# q1 N+ J: K' I/ q( x4 j tau = i-1;
' C+ w1 Y+ g7 u& V break;) j9 M: d: q" {. X: @
end
# a! E$ X+ A9 a3 g$ s6 I end
2 C0 `9 y" p5 R8 K( w + ]; D8 [7 ]0 z! w! W2 i" f
plot(0:length(entropy)-1,entropy)7 f1 y" e: e& g$ O
xlabel('tau');, f5 V' ]+ ~7 E
title('互信息法求tau'); |
|
/1 
发表于 2022-3-28 09:14
收藏
淘帖
支持！
反对！