clear all; close all;
M=[12 1 6 2 4 3];
N=[1 12 2 6 3 4];
SNR_1=linspace(0,30,31);
SNR = 10.^(SNR_1/10);
C= zeros(6,31);
for i=1:6
    if N(i)>M(i)
        C(i,:) = M(i)*log2(1+(N(i)*SNR)/M(i));
    elseif M(i)>=N(i)
        C(i,:) = N(i)*log2(1+SNR);
    end
end
figure (1);
plot(SNR_1,C(1,:),'ko-');
hold on;
plot(SNR_1,C(2,:),'rs-');
hold on;
plot(SNR_1,C(3,:),'gd-');
hold on;
plot(SNR_1,C(4,:),'bh-');
hold on;
plot(SNR_1,C(5,:),'c*-');
hold on;
plot(SNR_1,C(6,:),'mx-');
hold on;
grid on
title('Independent Rayleigh fading channel, CDI/CSI mode')
xlabel('SNR in dB')
ylabel('MIMO capacity (bits/transmission)')
axis([0 30 0 32]);
legend('(M=12, N=1)','(M=1, N=12)','(M=6, N=2)','(M=2, N=6)','(M=4, N=3)','(M=3, N=4)');