|
|
max_tau = 100;
& V: v8 P3 w) U; V. \9 N Part=128;
" O+ r! f1 y9 l; N2 a7 U: z7 M data=xlsread('L.xlsx');
~- g+ E$ N& C& N: h) ~ loadlibrary('mutual');
; S: |) Q; G; T5 C( x7 c [entropy]=mutual(data);& j: U; N3 r2 Q4 f, m# M
for i = 1:length(entropy)-1
$ J I/ @5 I: @: E if (entropy(i)<=entropy(i+1))8 _/ U' w8 r, u/ k
tau = i-1; + }; `1 F; x% q7 b0 b9 c' B
break;
" H7 \7 m4 n' h end
" \" J: u* ^, i4 O" W0 { end" ]/ D1 M. X7 }/ Y4 e5 u
* b7 `; j% M) g1 Z8 n, M% @# {8 o
plot(0:length(entropy)-1,entropy)* w, p0 z9 T7 n. w& a
xlabel('tau');
) A6 M6 N5 ?3 ^ p z title('互信息法求tau'); |
|
/1 
发表于 2022-3-28 09:14
收藏
淘帖
支持！
反对！