Progam1

L=10000; % ввод числа испытаний
for snr=1:10; % начало цикла по ОСШ в дБ
sum=0;
for i=1:L; % начало цикла по числу испытаний
x=rand(1,1);
b=bitgenerator(x); % генерирование случайного бита (1 или 0) с помощью х
s=modulator(b); % формирование символа ФМ-2, соответствующего биту b
sigma=10.^(-snr*0.05); % вычисление среднеквадратического отклонения АБГШ для текущего значения ОСШ
n=sigma*randn(1,1); % генерирование значения АБГШ с заданным среднеквадратичным отклонением
y=s+n; % формирование принимаемого сигнала
s_=demodulator(y); % демодуляция
b_=demap(s_); % формирование оценки бита (отображение символа в бит)
e=errordetection(b,b_); % определение наличия ошибки
sum=sum+e; % подсчёт общего числа битовых ошибок
end
BER(snr) = sum/L % вычисление коэффициента битовых ошибок
h=10.^(snr*0.1); % перевод значений ОСШ из дБ в разы
BERtheory(snr)=0.5*erfc(sqrt(0.5*h)); % вычисление значений коэффициента ошибок для теоретической кривой
end % завершение цикла по ОСШ
semilogy(BER(1:10), '-*'); % построение кривой по результатам моделирования
%в логарифмическом масштабе
hold on; % включение режима построения следующей кривой на том же графике
semilogy(BERtheory(1:10), 'm-+'); %построение теоретической кривой в логарифмическом масштабе
hold off; % выключение функции hold on
title('Помехоустойчивость системы радиосвязи с ФМ-2 (Гауссовский канал)');
xlabel('Eb/N0, дБ'); % подпись оси х
ylabel('BER'); % подпись оси у
legend('Моделирование', 'Теоретическая кривая');
grid