clear all
EbN0=12;
M=2;

n=log2(M)
nb_symbole=n*1e5;
hmod=pskmod(M,pi/4);
hmod.SymbolOrder='gray';
hdemod=pskdemod(M,pi/4);
hdemod.SymbolOrder='gray';

input = randi([0 M-1], nb_symbole,1);
xm = modulate(hmod, input);
SNR=EnB0+10+log10(n)
xmnoise=awgn(xm,SNR,'measured');
scatterplot(xm);grid
scatterplot(xmnoise);grid
ydemod = demodulate(hdemod, xmnoise);
[ydemod(1:10) input(1:10)]
%nb_symbole*n==nb_bits !
TEB=biterr(ydemod,input)/(nb_symbole*n)
BERtheory = berawgn(EbN0,'psk',M,'nondiff')
BERerfc_BPSK = 1/2.*erfc(sqrt(10^(EbN0/10)))