pulses.m

%SANJEET KUMAR
%M.Tech(Telecommunication) 
%NIT Calicut

% for the given g T (t)= (1-cos2?t/T)
echo on
T=1;
M=4;
N=100;
dT=T/(N-1);
t=0:dT:T;
Es=T/2
fc=6/T;
t=0:dT:T;
u0=(.5-.5*cos(2*pi*t/T)).*sqrt(Es).*cos(2*pi*fc*t);
u1=(.5-.5*cos(2*pi*t/T)).*sqrt(Es).*cos(2*pi*fc*t+2*pi/M);
u2=(.5-.5*cos(2*pi*t/T)).*sqrt(Es).*cos(2*pi*fc*t+4*pi/M);
u3=(.5-.5*cos(2*pi*t/T)).*sqrt(Es).*cos(2*pi*fc*t+6*pi/M);
figure(1);
subplot(6,1,1);
plot(t,u0);
title('QPSK Constant Envelope Signal Waveforms');
xlabel('time,t(s)-->');
ylabel('u0(t)-->');
subplot(6,1,2);
plot(t,u1);
xlabel('time,t(s)-->');
ylabel('u1(t)-->');
subplot(6,1,3);
plot(t,u2);
xlabel('time,t(s)-->');
ylabel('u2(t)-->');
subplot(6,1,4);
plot(t,u3);
xlabel('time,t(s)-->');
ylabel('u3(t)-->');
x=(abs(fft(u0))).*(abs(fft(u0)));
u4=sqrt(2*Es/T)*cos(2*pi*fc*t);
y=(abs(fft(u4))).*(abs(fft(u4)));
subplot(6,1,5);
plot(x);