bpsk_rayleigh.m

BPSK over rayleigh channel performance estimation

clear all;
M=[1 2 4 8 16 32 64 128];
Nt=4096;
EbNo=10;
z=10.^(EbNo/10);
for j=1:8
    t=cputime;
    N(j)=Nt/M(j);
    Nor_std_div_t(j)=sqrt((4*z*(1+(M(j)-1)*(0.5-1/pi)))/Nt);
    for i=1:100
        [BER(j,i)]=c10_MCBPSKRayFlatrun(N(j),M(j),z)
    end;
    Avg_BER(j)=sum(BER(j,:))/100;
    var(j)=sum((BER(j,:)-Avg_BER(j)).^2)/100;
    Nor_std_div(j)=sqrt(var(j))/Avg_BER(j);
    e(j)=cputime-t;
    F(j)=var(j)*e(j);
end;
disp(Nor_std_div);
disp(Nor_std_div_t);
disp(e);
disp(F);
semilogy(M,Nor_std_div,'x',M,Nor_std_div_t,M,e,'-.',M,F,'-');