funevalueci.m

来自「相关积分计算matlab程序」· M 代码 · 共 23 行

function CI=FunEvalueCI(N,m,R,PS)

% PS为所构造的N×m维相空间矩阵
% Xi与Xj的距离

h1=0;              % Heaviside函数值累计求和
for i=1:N
    for j=i+1:N
        D1 =norm(PS(i,:) - PS(j,:),2);
%         D = 0;
%         for k=1:m
%             D=D+(x(i,k)-x(j,k))^2;
%         end
%         D1=D^0.5;
        if (R>D1)
            h=1;
        else
            h=0;
        end
        h1=h1+h;
    end
end
CI=2*h1/(N*(N-1));