need help mimo Vs. BER MATLAB code

mkhan said:

Hi;

****************************************************************
what i need is the mimo simulation that shows the increasing in capacity when using siso, simo, miso and mimo...
****************************************************************


Here is attached my personal matlab files (remember capmimo is capmimo.m file, put .m bcz EDA does not allow the .m extension) for the capacity comparison and a tutorial report i wrtoe in the beginning, i did this when i started working on MIMO. I hope this will give u some help. If you need some more help, please let me know.

I will recommend you to go through the book "Introduction to space time wireless communications" By Arogyasawami Paulraj et al, I have gone thorugh and have many MIMO related books with me, but the top three, I recommend are as follows.

No. 1 I already mentioned Paulraj's book.

No. 2, Wireless Communications by Andrea Goldsmith, but be careful the e-copy available on EDA board is evaluation copy and some concepts are not explained as they are explained in original hard copy that was printed by cambridge press in 2005.

No. 3, Fundamentals of Wireless Communications by David Tse and P. Viswanath, this is an excellent reference as well.

Please press helped me !!! this will help me.

Arif Khan

Click to expand...

mkhan said:

Hi,

I hope this time u will find the code. No sorry , what shud i do ?

Added after 2 minutes:

No no, i dont know i cant attach my m file, so i m just copying it here. what u can do is just copy this in an empty new m file in ur matlab editor. Ok , i hope u got what i m saying.

here is the code

******************************************************
clear all
clc
%Shannon capacity
snr=0;
for i = 1:10
snr = snr +2;
c=(log(1+10^(snr/10)))/log(2);
x(i)=snr;
y(i)=c;
end
figure
plot(x,y,'-','LineWidth',1.5)
hold on

% capacity of MIMO Link with NR=2, NT=2
NR=2;
rand('state',456321)
snr=0;
for i=1:100;
snr=snr+2;
for j=1:10000;
c(j)=(NR*log(1+(10^(snr/10))*abs(normrnd(0,1)))/log(2));
end
yy(i)=mean(c);
xx(i)=snr;
end
plot(xx,yy,':','LineWidth',1.5)

% capacity of MIMO Link with NR=3, NT=3
NR=3;
rand('state',456321)
snr=0;
for i=1:100;
snr=snr+2;
for j=1:10000;
c(j)=(NR*log(1+(10^(snr/10))*abs(normrnd(0,1)))/log(2));
end
yy(i)=mean(c);
xx(i)=snr;
end
plot(xx,yy,'-.','LineWidth',1.5)

% capacity of MIMO Link with NR=4, NT=4
NR=4;
rand('state',456321)
snr=0;
for i=1:100;
snr=snr+2;
for j=1:10000;
c(j)=(NR*log(1+(10^(snr/10))*abs(normrnd(0,1)))/log(2));
end
yy(i)=mean(c);
xx(i)=snr;
end
plot(xx,yy,'--','LineWidth',1.5)
xlabel('SNR(dB)')
ylabel('Capacity (bit/s/Hz)')
grid on
legend('Shannon Capacity','MIMO, NT=NR=2','MIMO, NT=NR=3','MIMO, NT=NR=4',2)
title('MIMO Capacity')
print -deps -tiff -r300 capmimo
>>
********************************************************************

Dont forget to press helped me !!!!

Arif Khan

Click to expand...