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('»¥ÐÅÏ¢·¨ÇóʱÑÓ');
此互信息法計(jì)算程序經(jīng)過(guò)監(jiān)測(cè)人網(wǎng)站團(tuán)隊(duì)進(jìn)行多次驗(yàn)證和改進(jìn)，其計(jì)算成果基本正確，但是還有存在有部分問(wèn)題，歡迎相關(guān)學(xué)者批評(píng)和指正。
本次講座的主要成果為此互信息法計(jì)算程序，大家可以仔細(xì)研究一下，如果能成果理解本程序，那就說(shuō)明你對(duì)互信息法求延遲時(shí)間的相關(guān)理論知識(shí)的認(rèn)識(shí)較為透徹了，本次講座完成。
以上僅為監(jiān)測(cè)人網(wǎng)站關(guān)于混沌時(shí)間序列中計(jì)算延遲時(shí)間的一點(diǎn)簡(jiǎn)單，本講座為監(jiān)測(cè)人網(wǎng)站原創(chuàng)，轉(zhuǎn)載請(qǐng)注明出處。
監(jiān)測(cè)人http://m.mocsil.com 專(zhuān)業(yè)監(jiān)測(cè)資源共享網(wǎng)站。
