>> [M N] = size(X);
R=zeros(M,M);
for i = 1:M
for j = 1:N
R(i,j)=sum(min(X(i,:),X(j,:))')/sum(max(X(i,:),X(j,:))');
end
end
>> R
R =
1.0000 0.4285 0.6066 0.6267 0.3841 0.1602
0.4285 1.0000 0.6807 0.4624 0.7220 0.3011
0.6066 0.6807 1.0000 0.6139 0.6428 0.3579
0.6267 0.4624 0.6139 1.0000 0.5864 0.3593
0.3841 0.7220 0.6428 0.5864 1.0000 0.3788
0.1602 0.3011 0.3579 0.3593 0.3788 1.0000