pass.m

this code for qpsk with awgn

clear all;
n=20;
samp=32;
data=randint(1,n)
T=10^-3
f=1000;
t=0:1/samp:1-(1/samp);
%map=2*data-1;
%dd=reshape(map,2,[]);
%SQ=dd(1,:)
%SI=dd(2,:)
SQ=data;
SI=data-1;
basis1=(sqrt(2/T))*cos(2*pi*f*t);
basis2=(-1)*(sqrt(2/T))*sin(2*pi*f*t);
S1=basis1*2;
S2=basis2*2;
bpskI=(S1'*SI);
bpskQ=(S2'*SQ);
bpsk=bpskI-bpskQ;
Tx=sign(bpsk)
Tx1=reshape(Tx,1,[])
Txx=(Tx1+1)/2;
%%PENERIMA%%
z=reshape(bpsk,32,[]);
z1=basis1*z;
z2=basis2*z;
Rx1=awgn(z1,80,6.02);
Rx2=awgn(z2,80,6.02);
zE1=sign(Rx1-Rx2);
Rxx=(zE1+1)/2;
%F=reshape(Rxx,1,[]);
BER=biterr(data,Rxx)