kendallrankcorr.m

加权核函数的静态杂波抑制效果以及残留噪声的高斯性和独立性检验实验结果演示程序。 另外

function [tstar,meanval] = kendallrankcorr(x)
[h,w]=size(x);
tmp=0;
g=0;
% t0=n*(n-1);
% t0=2/t0;
en1=0;
en2=0;

t1=4*h+10;
t1=t1/(9*h*(h-1));
t1=sqrt(t1);

meanval=0;
idx=0;
for coli=1:w-1
 for colj=(coli+1):w
    idx=idx+1;
    for i=1:h-1
      for j=i+1:h
        a1=x(i,coli)-x(j,coli);
        a2=x(i,colj)-x(j,colj);
        tmp=a1*a2;
        if tmp~=0
           en1=en1+1;
           en2=en2+1;
          if tmp>0
            g=g+1;
          else
            g=g-1;
          end
        else
           if a1~=0 en1=en1+1;end
           if a2~=0 en2=en2+1;end               
        end  
      end
   end
   t(idx)=g/sqrt(en1*en2);
   %t(idx)=g*t0;
   meanval=meanval+t(idx);
   tstar(idx)=t(idx)/t1;
   g=0;
   en1=0;
   en2=0;
  end
end

meanval=meanval/idx;

plotdens(tstar);