distance.m

混沌的G-P算法

function [C,D]=Distance(Q,D,m,N,r)  % 由重构矩阵Q计算距离矩阵D
b=1/(N*N);
for k=1:m;

   for i=1:N
      for j=i+1:N % 由于距离矩阵D为对称矩阵，此处将D的下三角及其对角线上的计算过程去掉，减少1倍数据量
        DD(i,j)=max(abs((Q(i,k)-Q(j,k)))); % 重构矩阵第m列的第i行与第j行的距离(|ri-rj|)
        if DD(i,j)>D(i,j)
            D(i,j)=DD(i,j);
        else
        end
      end
   end
%%%%%%%%%%%%%%%%%%%%%%%
for i=1:200
   S=0;       
  for ii=1:N
    for jj=ii+1:N
        S=S+(1&((r(i)-D(ii,jj))>=0)); % 计算Heaviside函数
    end
  end
  c=2*b*S;
  lnc=log(c);
  C(i,k)=lnc;
end
end