b=[1,1,0,0,0,1,0,1]
fs=40000;
f0=4000;
br=2000;
M=fs/br;
xm=kron(b,ones(1,M));
N=max(size(b))*M;
n=0:N-1;
E0=0;
Sask=1*cos(2*pi*n*fc/fs).*xm;
plot(xm);
pause;
plot(Sask);
%srt(2Ei/T)=1
fi=pi
Sbpsk=cos(2*pi*n*fc/fs+fi);
pause
plot(Sbpsk)