Отчет по лабораторной №4 МСвСПД
.docxМИНИСТЕРСТВО ЦИФРОВОГО РАЗВИТИЯ И МАССОВЫХ КОММУНИКАЦИЙ РОССИЙСКОЙ ФЕДЕРАЦИИ
Федеральное государственное бюджетное образовательное учреждение
высшего образования
«Санкт-Петербургский государственный университет телекоммуникаций
им. проф. М. А. Бонч-Бруевича»
Кафедра сетей связи и передачи данных
Отчет по лабораторной работе №4
«Исследование канала ДСК по методу Монте-Карло в системе Octave»
По дисциплине: «Многофункциональный синтез в системах передачи данных»
Цель работы:
Ознакомиться с общими принципами проведения анализа по методу МонтеКарло на примере проверки правильности модели канала ДСК, реализованной в системе Octave.
Ход выполнения работы:
Вариант 7
Листинг:
Скрипт 5 раз пересылающий через канал ДСК 50000 случайных двоичных цифр.
p0=0.33
N= zeros(8,7);
x= [1 10 50 100 500 1000 5000 10000 50000];
for i1 = 1:8
c_error = 0
for i2 = 1:5
for i3= x(i2):x(i2+1)
b = randint(1);
if (bsc(b,p0) ~= b)
c_error = c_error+1;
endif
endfor
i1
N(i1,i2) = c_error/i3;
endfor
endfor
Расчет стандартного отклонения и математического ожидания.
for i=1:8
N(i,6)= mean(N(i,1:5));
N(i,7) = std(N(i,1:5));
endfor
Построение графиков экспериментальной вероятности ошибки в канале.
subplot (2,1,1)
semilogx(N(:,1:5),"r",N(:,6),"g")
xlabel ("t")
legend("m1","m2","m3","m4","m5","ExValue",4)
Рисунок 1. Графики экспериментальной вероятности ошибки в канале.
Построение графика стандартного отклонения в каждой точке.
subplot(2,1,2)
semilogx(N(:,7),"b")
xlabel ("t")
legend("StDev",4)
Рисунок 2. Графики стандартного отклонения в каждой точке.
Вывод программы в командном окне
>> Lab_3_variant_7
p0 = 0.33000
c_error = 0
i1 = 1
i1 = 1
i1 = 1
i1 = 1
i1 = 1
c_error = 0
i1 = 2
i1 = 2
i1 = 2
i1 = 2
i1 = 2
c_error = 0
i1 = 3
i1 = 3
i1 = 3
i1 = 3
i1 = 3
c_error = 0
i1 = 4
i1 = 4
i1 = 4
i1 = 4
i1 = 4
c_error = 0
i1 = 5
i1 = 5
i1 = 5
i1 = 5
i1 = 5
c_error = 0
i1 = 6
i1 = 6
i1 = 6
i1 = 6
i1 = 6
c_error = 0
i1 = 7
i1 = 7
i1 = 7
i1 = 7
i1 = 7
c_error = 0
i1 = 8
i1 = 8
i1 = 8
i1 = 8
i1 = 8
Вариант 18
Листинг:
Скрипт 5 раз пересылающий через канал ДСК 50000 случайных двоичных цифр.
p0=0.22
N= zeros(8,7);
x= [1 10 50 100 500 1000 5000 10000 50000];
for i1 = 1:8
c_error = 0
for i2 = 1:5
for i3= x(i2):x(i2+1)
b = randint(1);
if (bsc(b,p0) ~= b)
c_error = c_error+1;
endif
endfor
i1
N(i1,i2) = c_error/i3;
endfor
endfor
Расчет стандартного отклонения и математического ожидания.
for i=1:8
N(i,6)= mean(N(i,1:5));
N(i,7) = std(N(i,1:5));
endfor
Построение графиков экспериментальной вероятности ошибки в канале.
subplot (2,1,1)
semilogx(N(:,1:5),"r",N(:,6),"g")
xlabel ("t")
legend("m1","m2","m3","m4","m5","ExValue",4)
Рисунок 1. Графики экспериментальной вероятности ошибки в канале.
Построение графика стандартного отклонения в каждой точке.
subplot(2,1,2)
semilogx(N(:,7),"b")
xlabel ("t")
legend("StDev",4)
Рисунок 2. Графики стандартного отклонения в каждой точке.
Вывод программы в командном окне
>> Lab_3_variant_18
p0 = 0.22000
c_error = 0
i1 = 1
i1 = 1
i1 = 1
i1 = 1
i1 = 1
c_error = 0
i1 = 2
i1 = 2
i1 = 2
i1 = 2
i1 = 2
c_error = 0
i1 = 3
i1 = 3
i1 = 3
i1 = 3
i1 = 3
c_error = 0
i1 = 4
i1 = 4
i1 = 4
i1 = 4
i1 = 4
c_error = 0
i1 = 5
i1 = 5
i1 = 5
i1 = 5
i1 = 5
c_error = 0
i1 = 6
i1 = 6
i1 = 6
i1 = 6
i1 = 6
c_error = 0
i1 = 7
i1 = 7
i1 = 7
i1 = 7
i1 = 7
c_error = 0
i1 = 8
i1 = 8
i1 = 8
i1 = 8
i1 = 8
Вариант 24
Листинг:
Скрипт 5 раз пересылающий через канал ДСК 50000 случайных двоичных цифр.
p0=0.16
N= zeros(8,7);
x= [1 10 50 100 500 1000 5000 10000 50000];
for i1 = 1:8
c_error = 0
for i2 = 1:5
for i3= x(i2):x(i2+1)
b = randint(1);
if (bsc(b,p0) ~= b)
c_error = c_error+1;
endif
endfor
i1
N(i1,i2) = c_error/i3;
endfor
endfor
Расчет стандартного отклонения и математического ожидания.
for i=1:8
N(i,6)= mean(N(i,1:5));
N(i,7) = std(N(i,1:5));
endfor
Построение графиков экспериментальной вероятности ошибки в канале.
subplot (2,1,1)
semilogx(N(:,1:5),"r",N(:,6),"g")
xlabel ("t")
legend("m1","m2","m3","m4","m5","ExValue",4)
Рисунок 1. Графики экспериментальной вероятности ошибки в канале.
Построение графика стандартного отклонения в каждой точке.
subplot(2,1,2)
semilogx(N(:,7),"b")
xlabel ("t")
legend("StDev",4)
Рисунок 2. Графики стандартного отклонения в каждой точке.
Вывод программы в командном окне
>> Lab_3_variant_24
p0 = 0.16000
c_error = 0
i1 = 1
i1 = 1
i1 = 1
i1 = 1
i1 = 1
c_error = 0
i1 = 2
i1 = 2
i1 = 2
i1 = 2
i1 = 2
c_error = 0
i1 = 3
i1 = 3
i1 = 3
i1 = 3
i1 = 3
c_error = 0
i1 = 4
i1 = 4
i1 = 4
i1 = 4
i1 = 4
c_error = 0
i1 = 5
i1 = 5
i1 = 5
i1 = 5
i1 = 5
c_error = 0
i1 = 6
i1 = 6
i1 = 6
i1 = 6
i1 = 6
c_error = 0
i1 = 7
i1 = 7
i1 = 7
i1 = 7
i1 = 7
c_error = 0
i1 = 8
i1 = 8
i1 = 8
i1 = 8
i1 = 8
Вывод:
В ходе проделанной работы был проведен анализ канала ДСК по методу Монте-Карло, в результате которого были получены график экспериментальной вероятности ошибки в канале ДСК и график стандартного отклонения.