3

pkg load communications
N=50000;
n=5;
p=(40-4)./100;
SumOfErrors=0;
SumOfProbabilities=0;
M(1:N,1:n)=0;
-------------------------------------------------------------------
P=bsc(M,p);
-------------------------------------------------------------------
for i=1:1:n
for j=1:1:N
SumOfErrors=SumOfErrors+abs(M(j,i)-P(j,i));
w(j,i)=SumOfErrors./j;
if (j==N)
SumOfErrors=0;
endif
end
end
-------------------------------------------------------------------
for i=1:1:n
for j=1:1:N
SumOfProbabilities=SumOfProbabilities+w(j,i);
end
end
mean=SumOfProbabilities./(N.*n);
-------------------------------------------------------------------
x=[10 50 100 500 1000 5000 10000 50000];
y1=[w(10,1), w(50,1), w(100,1), w(500,1), w(1000,1), w(5000,1), w(10000,1), w(50000,1)];
y2=[w(10,2), w(50,2), w(100,2), w(500,2), w(1000,2), w(5000,2), w(10000,2), w(50000,2)];
y3=[w(10,3), w(50,3), w(100,3), w(500,3), w(1000,3), w(5000,3), w(10000,3), w(50000,3)];
y4=[w(10,4), w(50,4), w(100,4), w(500,4), w(1000,4), w(5000,4), w(10000,4), w(50000,4)];
y5=[w(10,5), w(50,5), w(100,5), w(500,5), w(1000,5), w(5000,5), w(10000,5), w(50000,5)];
y=mean;
semilogx(x,y,'o-r',x,y1,x,y2,x,y3,x,y4,x,y5)
set(gca, 'XTick', [0 10 50 100 500 1000 5000 10000 50000]);
xlabel('N');
ylabel('frequency');
title('transmission error rate');
-------------------------------------------------------------------
q=[std(w(:,1)),std(w(:,2)),std(w(:,3)),std(w(:,4)),std(w(:,5)),];
bar(q)
title('standard deviation');