% Công suất phát theo dB
PdB = 0:1:10;
P = 10.^(PdB/10);
% Độ lợi kênh truyền trung bình
Omega = [1 2 3]; %
N =
length(Omega);
% Lý thuyết
g_ = Omega(:).*P;
gMRC_a =
sum(g_,1);
% Mô phỏng
N = 10^4; % So luong phep thu
h1 =
sqrt(Omega(1)/2).*(randn(length(PdB),N) + 1i*randn(length(PdB),N));
h2 =
sqrt(Omega(2)/2).*(randn(length(PdB),N) + 1i*randn(length(PdB),N));
h3 =
sqrt(Omega(3)/2).*(randn(length(PdB),N) + 1i*randn(length(PdB),N));
% MRC
gMRC_s =
P(:).*sum(abs(h1).^2 + abs(h2).^2 + abs(h3).^2,2)/N;
% Vẽ
figure(1);
set(gcf,'color','White');
semilogy(PdB,gMRC_a,'-',PdB,gMRC_s,'o');
legend('Lý thuyết','Mô phỏng');
xlabel('Tỷ lệ tín hiệu trên nhiễu theo dB');
ylabel('Tỷ lệ tín hiệu trên nhiễu');