clear all;
close all;
clc;
M=1000;
Nt=[2,4,6,8];
Nr=[2,2,2,2];
s={'b-o','b-*','b-square','b-hexagram'};
SNR=[0:2:20]; %dB
figure(1);
title('接收天线数不变,改变发射天线数系统容量仿真分析')
xlabel('SNR [dB]');
ylabel('Capacity');
grid on
hold on;
for n = 1:length(Nt)
t = Nt(n);
r = Nr(n);
for m = 1:M
H = raylrnd(1,r,t);
for l = 1:length(SNR)
rou = 10^(SNR(l)/10);
% Find the capacity while CSIT unknown
CU(m,l) = log2(real(det(eye(r)+rou*H*H'/t)));
end
end
C_unknown(:,n) = mean(CU)';
for l = 1:length(SNR)
[cdf_u,co_u] = hist(CU(:,l),100); %绘制直方图,用来显示数据集的分布情况