function CorrMatrix_sqrt=getCorrMatrix(rho,AntennaConfig)

CorrMatrix=rho*ones(str2num(AntennaConfig),str2num(AntennaConfig));
for loop=1:str2num(AntennaConfig)
    CorrMatrix(loop,loop)=1;
end

[v,d]=eig(CorrMatrix);  % !Note: although svd can get the same result, eig is right in definition.
CorrMatrix_sqrt=v*sqrt(d)*v';