capacityuncorrelatedaveragehome.m

来自「MIMO v&#7899 i Zero bu&#7897 c equalizer」· M 代码 · 共 55 行

%% 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;
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 >> ');
it=3;
NumberOfRelizations=input('Input the number of Relizations >> ');
AR1=zeros(1,length(SNRArray));
AR2=zeros(1,length(SNRArray));
AR3=zeros(1,length(SNRArray));
%for ii=1:it
    MRx=4;%input('Input the number of receive antennas 1st Arrangment >> ');
    NTx=4;%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
AR1=AR1+Cmimo1st;
end
    MRx=4;%input('Input the number of receive antennas 2nd Arrangment >> ');
    NTx=2;%input('Input the number of transmit antennas 2nd 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)
    Cmimo2nd(snr)=log2(det(eye(MRx)+SNRArray(snr)/NTx*H*H'));
end
AR2=AR2+Cmimo2nd;
end
    MRx=2;%input('Input the number of receive antennas 3rd Arrangment >> ');
    NTx=4;%input('Input the number of transmit antennas 3rd 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)
    Cmimo3rd(snr)=log2(det(eye(MRx)+SNRArray(snr)/NTx*H*H'));
end
AR3=AR3+Cmimo3rd;
end
AR1=AR1/NumberOfRelizations;
AR2=AR2/NumberOfRelizations;
AR3=AR3/NumberOfRelizations;
%end
figure;
subplot(2,2,1);plot(SNRArraydB,AR1,'k',SNRArraydB,AR2,'b',SNRArraydB,AR3,'r');