лабы / lab2

pkg load communications
function OUT = ZC(IN, P)
OUT = and(IN, rand(size(IN)) > P);
endfunction

p = [1e-2 5e-2 1e-1 2e-1 3e-1 7e-1 9e-1]; % Вероятности ошибки
d1 = randint(1,1e5);
d2 = zeros(1,1e5);
d3 = ones(1,1e5);
% Массивы статистики (5 строк экспериментальных значений и среднее)
statForRandint = zeros(6, length(p));
statForZeros = zeros(6, length(p));
statForOnes = zeros(6, length(p));
% Основной цикл
for i=1:1:length(p)
% 5 экспериментов
for j=1:1:5
[a statForRandint(j, i)] = biterr(d1, ZC(d1, p(i)));
[a statForZeros(j, i)] = biterr(d2, ZC(d2, p(i)));
[a statForOnes(j, i)] = biterr(d3, ZC(d3, p(i)));
endfor
% Вычисление среднего
statForRandint(6,i) = mean(statForRandint(1:5,i));
statForZeros(6,i) = mean(statForZeros(1:5,i));
statForOnes(6,i) = mean(statForOnes(1:5,i));
endfor
% Построение графиков
subplot(2, 3, 1);
h = semilogx(p(:),statForRandint(1:5,:),"b",p,statForRandint(6,:),"r");
grid;
legend ([h(1), h(6)], {"experiments", "mean"}, "location", "northeast");
title("ZCH random value IKTU-98 Boyko, Amosova");
xlabel('Given probability');
ylabel('Experimental probability');
subplot(2, 3, 2);
h = semilogx(p(:),statForZeros(1:5,:),"b",p,statForZeros(6,:),"r");
grid;
legend ([h(1), h(6)], {"experiments", "mean"}, "location", "northeast");
title("ZCH zeros IKTU-98 Boyko, Amosova");
xlabel('Given probability');
ylabel('Experimental probability');
subplot(2, 3, 3);
h = semilogx(p(:),statForOnes(1:5,:),"b",p,statForOnes(6,:),"r");
grid;
legend ([h(1), h(6)], {"experiments", "mean"}, "location", "northeast");
title("ZCH ones IKTU-98 Boyko, Amosova");
xlabel('Given probability');
ylabel('Experimental probability');
subplot(2, 3, 4:6);
h = semilogx(p(:),statForRandint(6,:),"r",p(:),statForZeros(6,:),"b",p(:),statForOnes(6,:),"g");
grid;
legend ({"mean for random values array", "mean for zeros array", "mean for ones array"}, "location", "northeast");
title("ZCH ones IKTU-98 Boyko, Amosova");
xlabel('Given probability');
ylabel('Experimental probability');