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/
更多精彩文章请关注微信号: