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://jianceren.cn 专业监测资源共享网站。