H_sq=0; %¼ÆËã(s,q)µÄÁªºÏìØ
for i=1:n
for j=1:n
if Psq(i,j)~=0
H_sq=H_sq-Psq(i,j)*log(Psq(i,j));
end
end
end
I_sq(tau)=H_s+H_q-H_sq; %¼ÆËãtauϵĻ¥ÐÅÏ¢º¯Êý
clear s q; %Çå¿Õ±äÁ¿sºÍq
end
for i=2:length(I_sq)-1
if I_sq(i)<I_sq(i-1)&I_sq(i)<I_sq(i+1)
tau=i;
break;
end
end
plot(1:40,I_sq(1:40),'.-')
xlabel('Lag');
title('»¥ÐÅÏ¢·¨ÇóʱÑÓ');
此互信息法計算程序經過監測人網站團隊進行多次驗證和改進,其計算成果基本正確,但是還有存在有部分問題,歡迎相關學者批評和指正。
本次講座的主要成果為此互信息法計算程序,大家可以仔細研究一下,如果能成果理解本程序,那就說明你對互信息法求延遲時間的相關理論知識的認識較為透徹了,本次講座完成。
以上僅為監測人網站關于混沌時間序列中計算延遲時間的一點簡單,本講座為監測人網站原創,轉載請注明出處。
監測人http://m.mocsil.com 專業監測資源共享網站。