%Název: Stanovení poètu uživatelù v CDMA systému
%Autor: Jan Paseèný
%Datum: 28.5.2008
%Program pro MATLAB k bakaláøské práci FEKT VUT Brno

clear all;
N=60;                       %poèet uživatelù
alfa=0.375;                 %voice activity factor
delta=32;
%W=1.25*10^6;               %
%R=8*10^3;                  %èást pro komplexnìjší zadání vstupních velièin
%ebn0=5;                    %
%delta=(W/R)/ebn0;          %


hold off;
axis([20 60 -4 -1]);
hold on;
grid on;
xlabel('Poèet uživatelù na sektor'), ylabel('Pr(BER>0.001');
set(gca,'YTick',-4:1:-1);
set(gca,'YTickLabel',{'0.0001','0.001','0.01','0.1'});

for delta=32 %:15:47    %pokud chceme zobrazit souèasnì více závislostí
%for delta              %èást pro komplexnìjší zadání vstupních velièin 
  for NS=20:1:N
    s=0;
    for k=0:1:NS-1, 
        s=s+nchoosek(NS,k)*(alfa^k)*((1-alfa)^(NS-1-k))*0.5*erfc((delta-k-0.274*NS)/sqrt(0.156*NS));
    end;
    P(NS)=s;
  end;  

  plot(20:1:60,log10(P(20:1:N)));
end;


