ip_08_07.m

来自「Proakis《contemporarycommunication system」· M 代码 · 共 16 行

% MATLAB script for Illustrative Problem 7, Chapter 8. 
echo on
ep=0.3;
for i=1:2:61
  p(i)=0;
  for j=(i+1)/2:i
    p(i)=p(i)+prod(1:i)/(prod(1:j)*prod(1:(i-j)))*ep^j*(1-ep)^(i-j);
  end
end
pause % Press a key to see the plot
stem((1:2:61),p(1:2:61))
xlabel('n')
ylabel('pe')
title('Error probability as a function of n in simple repetition code')