clear;
% AWGN CCMC capacity 
snr_dB = -10:1:25;
snr = 10.^(snr_dB/10);
C = log2(1+snr);
plot(snr_dB,C,'-d','color','k','MarkerSize',10);
grid on
xlabel('SNR (dB)','FontSize', 16);
ylabel('Bandwidth Eficiency (bit/s/Hz)','FontSize', 16);
legend('AWGN: 1x1','Location','northwest');
axis([snr_dB(1) snr_dB(end) 0 10])
ax = gca; ax.FontSize = 16; grid on;
 
hold on

% SISO Rayleigh CCMC capacity 
N=5000;
h=sqrt(1/2)*(randn(N,1)+1i*randn(N,1));
factor = abs(h).^2;

% update the following line and put your codes here
C2 = ones(1,length(snr));
%%%%%

plot(snr_dB,C2,'-+','DisplayName','Rayleigh: 1x1','color','b','MarkerSize',10);

hold off;
