e4_1.rar_MIMO容量和SNR_SNR和MIMO容量_mimo_mimo 平均 容量_平均 容量

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;