Mô phỏng xác suất dừng của hệ thống maxial ratio combining

% 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');