passband.asv

this code for qpsk with awgn

clear all;clf
samp=32;
data=randint(1,20,2,1)
T=10^-3;
f=10;
t=0:1/samp:1-(1/samp);
SQ=data;
SI=1-data;
basis1=(sqrt(2/T))*cos(2*pi*f*t);
basis2=(-1)*(sqrt(2/T))*sin(2*pi*f*t);
S1=basis1;
S2=basis2;
qpskI=(S1'*SI);
qpskQ=(S2'*SQ);
qpsk=qpskI+qpskQ;    %%%%%keluaran pemancar
subplot()

qpsknoise=awgn(qpsk,-50,6.02)    %%%%penambahan noise AWGN (diumpamakan di kanal)

z1=(basis1*qpsknoise)
z2=(basis2*qpsknoise)
z=sign(z1-z2)
z=(1-z)/2

%menghitung BER
[N,BER]=biterr(z,data)

%Tx=sign(bpsk)
%Tx1=reshape(Tx,1,[]);
%Txx=(Tx1+1)/2;
%plot(0:length(bpsk)-1,bpskI); 
%grid on;
%sket=SQ-SI;
%scatterplot(sket);
%%Penerima%%
%z=reshape(bpsk,32,[]);
%z1=basis1*z;
%z2=basis2*z;
%zE1=sign([z1;z2]*(1/128));

%%% Noise %%%%%%

%snr= 10; %%%dapat divariasi nilai snr 0,3,6,10 %%%%
%Rx1=awgn(bpsk,snr,'measured')
%Rx=sign(Rx1)
%Rx=reshape(Rx,1,[]);
%Rxx=(Rx+1)/2;
%BER=biterr(Txx,Rxx)
%subplot (2,1,1)
%plot(t,S1)
%axis([0 1 -3 3])
%subplot (2,1,2)
%plot(t,S2)
%axis([0 1 -3 3 ])