M=1000;
n_bins=round(M/10);
Nt=[2:2:8];
SNR=[10];
figure(1);
title('发送端未知CSI情况下MIMO系统容量CCDF','FontSize',14)
xlabel('容量（bit/s/Hz）','FontSize',12);
ylabel('Pr（容量>=横坐标）','FontSize',12);
hold on;
text(1.8,0.8,'1×1','FontSize',12);
text(3,0.9,'2×2','FontSize',12);
text(6,0.8,'2×4','FontSize',12);
text(8,0.7,'2×6','FontSize',12);
text(10.4,0.2,'2×8','FontSize',12);
text(13,0.1,'4×4','FontSize',12);
text(19,0.1,'6×6','FontSize',12);
text(24,0.1,'8×8','FontSize',12);
%画1发1收时SISO系统容量
for n=1:1
for m=1:M
H=(randn(1,1)+1i*randn(1,1))/sqrt(2);
rho=10^(SNR/10);
CU(m,n)=log2(real(det(eye(1)+rho*H*H')));
end
[cdf,c]=hist(CU(:,n),n_bins);
plot(c,1-(cumsum(cdf))/M)
end
%画收发天线个数为2×2，2×4，2×8时的容量
for n=2:length(Nt)
N=Nt(n);
for m=1:M
H=(randn(2,N)+1i*randn(2,N))/sqrt(2);
rho=10^(SNR/10);
CU(m,n)=log2(real(det(eye(2)+rho*H*H'/2)));
end
[cdf,c]=hist(CU(:,n),n_bins);
plot(c,1-(cumsum(cdf))/M)
end
%
for n=1:length(Nt)
N=Nt(n);
for m=1:M
H=(randn(N,N)+1i*randn(N,N))/sqrt(2);
rho=10^(SNR/10);
CU(m,n)=log2(real(det(eye(N)+rho*H*H'/N)));
end
[cdf,c]=hist(CU(:,n),n_bins);
plot(c,1-(cumsum(cdf))/M)
end
grid on;
