lr4_3

clear all
close all
clc
N = 10^6; % число символов для оценки ПУ
Eb_N0_dB = [0:20]; % диапазон значений ОСШ
nTx = 2; % число передающих антенн
ip = rand(1,N)>0.5; % формирование битового потока
s = 2*ip-1; % отображение битов в символы ФМ-2
for ii = 1:length(Eb_N0_dB) % цикл по ОСШ
% формирование коэффициентов замираний и АБГШ
h = 1/sqrt(2)*[randn(nTx,N) + j*randn(nTx,N)];
n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)];
% Разнесенная передача с ОС
sr = (1/sqrt(nTx))*kron(ones(nTx,1),s);
hEff = h.*exp(-j*angle(h));
% добавление АБГШ и замираний
y1 = sum(h.*sr,1) + 10^(-Eb_N0_dB(ii)/20)*n;
y2 = sum(hEff.*sr,1) + 10^(-Eb_N0_dB(ii)/20)*n;
% Комбинирование сигналов
y1Hat = y1./sum(h,1); % (без ОС)
y2Hat = y2./sum(hEff,1); % (с ОС)
% приемник с жесткими решениями
ip1Hat = real(y1Hat)>0; % (без ОС)
ip2Hat = real(y2Hat)>0; % (с ОС)
% подсчет числа ошибок
nErr1(ii) = size(find([ip- ip1Hat]),2); % (без ОС)
nErr2(ii) = size(find([ip- ip2Hat]),2); % (с ОС)
end
simBer1 = nErr1/N; % (без ОС)
simBer2 = nErr2/N; % (с ОС)
EbN0Lin = 10.^(Eb_N0_dB/10); theoryBer = 0.5.*(1-sqrt(EbN0Lin./(EbN0Lin+1)));
p = 1/2 - 1/2*(1+1./EbN0Lin).^(-1/2); theoryBer_nRx2 = p.^2.*(1+2*(1-p));
close all; figure;
semilogy(Eb_N0_dB,theoryBer,'k','LineWidth',2); hold on;
semilogy(Eb_N0_dB,theoryBer_nRx2,'mx-','LineWidth',2);
semilogy(Eb_N0_dB,simBer1,'gs--','LineWidth',2);
semilogy(Eb_N0_dB,simBer2,'cs-.','LineWidth',2);
axis([min(Eb_N0_dB) max(Eb_N0_dB) 10^-5 0.5]); grid on;
legend('формула (M_{Tx}=1,N_{Rx}=1)','формула (M_{Tx}=1,N_{Rx}=2, MRC)',...
'ИМ (M_{Tx}=2,N_{Rx}=1) (без ОС)','ИМ (M_{Tx}=2,N_{Rx}=1) (с ОС)');
xlabel('E_{b}/N_{0}, дБ'); ylabel('BER (вероятность ошибки)');
title('Помехоустойчивость разнесенной передачи с ОС для MISO 2*1');