clc;close all;clear all;

%%%%%%%%%%%%%%%% Comparison Between Capacity of SISO and MIMO  %%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%  Yousif Allbadi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

r_M = 2; % Number of Receive Antennas of MIMO

t_M = 2; % Number of Transmit Antennas of MIMO

r_S = 1; % Number of Receive Antennas of SISO

t_S = 1; % Number of Transmit Antennas of SISO

N_M = min(r_M,t_M);

N_S = min(r_S,t_S);

H_M = zeros(r_M,t_M); % H-Matrix of Channel Coefficient of MIMO

H_S = zeros(r_S,t_S); % H-Matrix of Channel Coefficient of SISO

SNR_dB =[0:0.01:35]; % SNR in dB as a Vactor

SNR = 10.^(SNR_dB/10); % Absolute Value of SNR

for k=1:r_M

   for i=1:t_M

       H_M(k,i) = randn(1) + j* randn(1);

   end

end

for k=1:r_S

   for i=1:1

       H_S(k,i) = randn(1) + j* randn(1);

   end

end

[S V D] = svd(H_M*H_M');

[S_S V_S D_S] = svd(H_S*H_S');

Capacity_MIMO = 0;

Capacity_SISO = 0 ;

for k=1:N_M

   lamda(k) = V(k,k);

   Capacity_MIMO = Capacity_MIMO + log2(1+SNR*lamda(k)/t_M); % Capacity of MIMO

end

 for k=1:N_S

   lamda_S(k) = V_S(k,k);

   Capacity_SISO = Capacity_SISO + log2(1+SNR*lamda_S(k/t_S)); % Capacity of SISO

end

%Capacity_SISO = log2(1+SNR); % Shannon Formula

figure(1);

plot(SNR_dB,Capacity_SISO,'-','LineWidth',1.5,'Color','b');

hold on

plot(SNR_dB,Capacity_MIMO,'-','LineWidth',1.5,'Color','r');

legend('SISO','MIMO');

xlabel('SNR (dB)')

ylabel('Capacity(b/s/Hz)')

title('Comparison Between Capacity of SISO and MIMO')

grid;

源码下载地址：

http://page2.dfpan.com/fs/elcc3j5262f192b9166/ 

更多精彩文章请关注微信号：