for i = 1:M
for j = N:-1:1
norm = 0;
for r = 1:R
for t = 1:T
X = x(i,r,t);
for n = 1:N
if (n ~= j)
X = X - coupling(i,n) *
alpha(n,r) *
s (n,pad+t-tau(n,r));
end
end
Y = alpha(j,r)*s(j,pad+t-tau(j,r));
c(i,j) = c(i,j) + X*Y;
norm = norm + Y^2;
end
end
c(i,j) = c(i,j)/norm;
end
end