capacityuncorrelatedaveragehome2.m

MIMO v&#7899 i Zero bu&#7897 c equalizer

%% This file will plot the average Capacity of MIMO Channels.
%% The plots will be Capacity vs SNR, Capacity VS MRx Capacity vs. NTx, 
%% Capacity comparisons with othe channels
close all;
clear all;
type=['o' '+' 's' 'd'];
colour=['k' 'r' 'b' 'k'];
disp('1- Capacity vs SNR input parameters');
SNRArraydB=input('Input the SNR Range in dBs >> ');
SNRArray=10.^(SNRArraydB/10);
it=input('Input the number of antenna arrangments up to four >> ');
%it=3;
NumberOfRelizations=input('Input the number of Relizations >> ');
%AR1=zeros(1,length(SNRArray));
%AR2=zeros(1,length(SNRArray));
%AR3=zeros(1,length(SNRArray));
figure;
subplot(2,2,1);hold;
for ii=1:it
    AR(ii,:)=zeros(1,length(SNRArray));
    MRx=input('Input the number of receive antennas for the &ii Arrangment >> ');
    NTx=input('Input the number of transmit antennas 1st Arrangment >> ');
    clear H;
    for jj=1:NumberOfRelizations
    H=sqrt(0.5)*randn(MRx,NTx)+i*sqrt(0.5)*randn(MRx,NTx);
    for snr=1:length(SNRArray)
    Cmimo1st(snr)=log2(det(eye(MRx)+SNRArray(snr)/NTx*H*H'));
end
AR(ii,:)=AR(ii,:)+Cmimo1st;
end
AR(ii,:)=AR(ii,:)/NumberOfRelizations;
plot(SNRArraydB,AR(ii,:),colour(ii),SNRArraydB,AR(ii,:),type(ii));
title(sprintf('Capacity vs. SNR for different {\it N^{Tx}X M^{Rx}}');
xlabel('SNR dB');
ylabel('Capacity bps/Hz');
legend(sprintf('%2.0fX%2.0f',[NTx MRx]));
end