说明：  for K=1:10000 T=randn(M,M)+j*randn(M,M); T=0.707*T; I=eye(M); eigen=eig(T*T'); %extract eigenvalues gamma=zeros(M,1); r=M;%set rank = number of antennas (full rank) p=1;%initial count for i=1:r mu=getmu(r,SNR,T,p,M);%determine mu value (see eqn. (1.30)) gamma(i)=mu-(M/SNR)*(1/eigen(i)); %calculate gamma (see eqn (1.31)) if gamma(i)<0 gamma(i)=0;% if gamma < 0, set it to zero i.e. discard it p=p+1;% increment count mu=0;%clear register else mu=0; % if gamma >0 store it and clear register end end %solve eqn (1.26)using the determinant form a=I+(SNR/M)*diag(gamma).*diag(eigen); a=det(a); y(K)=log2(a); end