% QPSK Demod Test
clear;
load c:\qpsk\lecroy_1.dir\sc1005.dat;
offset=1;
size=1024*4;
wave=sc1005';
fs=1e9;

% Display wave form

plot(wave(1:100));
grid
title('IF Signal');
%%%%%%pause;

data_fft=10*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
      (wave(offset:offset+size-1))))));
   data_fft=data_fft-max(data_fft);
   nnnn=size;
   ff=-nnnn/2:nnnn/2-1;
   ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
	title(' IF Data Spectrum Before Filtering ');
	axis([400e6, 500e6 -50 0]);
	grid;
	xlabel('FREQUENCY (HZ)');
	ylabel('MAGNITUDE (dB) ');
   %pause;
   
   % generate local carrier
   
   f_lo=(480-50)*1e6;   % carrier
   t=[0:length(wave)-1]/fs;
	carrier_a=sin(2*pi*f_lo*t); 	% carrier sine;

	if_data=wave.*carrier_a;	% down converting

%	% plot spectrum;
data=if_data;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
		 data(offset:offset+size-1)))));
	ff=-nnnn/2:nnnn/2-1;
	ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
	title(' I Data Spectrum After Down Converting');
	grid;
   %%%pause;
   
   axis([0 100e6 -10 30]);
   %%%pause;
   
      
   %************************************
   %		Design A Low-Pass Filter 		*
   %************************************
   
   f_lowpass=[   0,    65/500,   80/500 ,    1.     ];
        m=   [   1.,   1.0 ,       1e-5,       1e-5   ];

n1=80;
b=firls(n1, f_lowpass, m);

a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10(  (abs(amp)+1e-8)  /  (max(abs(amp)) +1e-8)  );
plot(freq,amp);
grid
%%%pause  
axis([0e6 100e6 -50 0]);
title('N1=100 taps FIR low pass filter.');
%%%pause;

if_data=filter(b,a,if_data);

if_50_data=if_data;

   %***************************
   %	Plot Spectrum				*
   %***************************
   data=if_data;
   offset=100;
   data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
      data(offset:offset+size-1)))));
	ff=-nnnn/2:nnnn/2-1;
	ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
	title(' IF Data Spectrum After Low-Pass Filter');
	grid;
   %%%pause;
   axis([0e6 100e6 -50 30]);
	%%%pause;
   
   %***************************
   %		Times 2					*
   %***************************
   
   if2_data=if_data.*if_data;
   
   %***************************
   %	Plot Spectrum				*
   %***************************
   offset=100;
   data=if2_data;
   data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
      		data(offset:offset+size-1)))));
	ff=-nnnn/2:nnnn/2-1;
	ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
   axis([0e6 250e6 -50 30]);
	title(' IF Data Spectrum After X2');
	grid;
   %pause;
   
   %***************************
   %		Times 4					*
   %***************************
   
   if4_data=if2_data.*if2_data;
   
   %***************************
   %	Plot Spectrum				*
   %***************************
   offset=100;
   data=if4_data;
   data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
      		data(offset:offset+size-1)))));
	ff=-nnnn/2:nnnn/2-1;
	ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
   axis([200e6 250e6 -100 -30]);
	title(' IF Data Spectrum After X4');
	grid;
   %pause;
   

   
   
   %*********************************
   %		Design A Low-Pass Filter	*
   %*********************************
      b=fir1(400,[213/500 220/500]);

   a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10(  (abs(amp)+1e-8)  /  (max(abs(amp)) +1e-8)  );
plot(freq,amp);
title(' Band-Pass Filter After X2');
axis([200e6 250e6 -50 0]);
grid
%%pause  

if4_data=filter(b,a,if4_data);

plot(if4_data(1000:3100));
grid;
title('IF Waveform X4');
%%pause;


	%***************************
   %		sign  					*
   %***************************
   
   if4_data=sign(if4_data);
   
   %***************************
   %	Plot Spectrum				*
   %***************************
   offset=100;
   data=if4_data;
   data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
      		data(offset:offset+size-1)))));
	ff=-nnnn/2:nnnn/2-1;
	ff=ff/nnnn*2*fs/2;
   plot(ff,data_fft);
   %axis([200e6 250e6 -100 -30]);
	title(' IF Data Spectrum After X4');
	grid;
   %pause;
   

   
   
if4_data=filter(b,a,if4_data);

plot(if4_data);
grid;
title('IF Waveform X4 After Sign 1');
%pause;



for ii=1:10
   step=300;
plot(if4_data(100+(ii-1)*step:ii*step+100));
grid;
title('IF Waveform after 1 sign ');
%%pause;
end; 



%***********************
% count cycle
%***********************

carrier_ref=if4_data(1000:5000-100);
length_ref=length(carrier_ref);
cycle_count=0;
for i=1:length_ref-1;
   if carrier_ref(i)*carrier_ref(i+1) <=0
      cycle_count=cycle_count+1;
   end
end

carrier_cycle=(length_ref-1)/fs/cycle_count*2;
carrier_frequency=1/carrier_cycle;

%*******************
%	Generate LO	   *
%*******************
length_if=length(if4_data);
t=[0:length_if-1]/fs;
sin_lo=sin(2*pi*(carrier_frequency+1e3)*t);

for ii=1:20
   step=200;
plot(if4_data(500+(ii-1)*step:ii*step+500));
hold on;
aaa=5;
plot(sin_lo(500+aaa+(ii-1)*step:ii*step+500+aaa),'r');
hold off;
grid;
title('Carrier and Signal Waveform ');
%pause;
end; 

true_carrier_frequency=carrier_frequency/4;

   %*******************
	%	Generate LO	   *
   %*******************
   for i_phase=34-10:34+10;
   
   length_wave=length(if_50_data);
   t=[0:length_wave-1]/fs;
   phase_offset=pi*(i_phase*5)/180;
   df=-0e3;
   lo_cos=cos(2*pi*(df+true_carrier_frequency)*t+phase_offset);
   lo_sin=sin(2*pi*(df+true_carrier_frequency)*t+phase_offset);
   
   %************************
   %	Check Carrier Phase
   %***********************
   
   %for i=1:10
   %  step=400;
	%plot(if_50_data(1+(i-1)*step:i*step)*8);
	%hold on;
	%aaa=-1;
   %plot(lo_sin(1+aaa+(i-1)*step:i*step+aaa),'r');
   %plot(lo_cos(1+aaa+(i-1)*step:i*step+aaa),'g');
	%hold off;
	%grid;
	%title('***** Check Carrier Phase ***** ');
	%%pause
	%end

   
   
   %****************
   %	Down Convert LO
   %***************
   
   i_data=if_50_data.*lo_cos;
   q_data=if_50_data.*lo_sin;
   
   
   %*********************************
   %		Design A Low-Pass Filter	*
   %*********************************
      b=fir1(100,25/500);

   a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10(  (abs(amp)+1e-8)  /  (max(abs(amp)) +1e-8)  );
plot(freq,amp);
title(' Low Pass Filter');
%axis([200e6 250e6 -50 0]);
grid
%%%pause  

i_data=filter(b,a,i_data);
q_data=filter(b,a,q_data);


for i=1:4
   step=500;
plot(i_data(100+(i-1)*step:100+i*step));
hold on;
plot(q_data(100+(i-1)*step:100+i*step),'r');
hold off;
grid;
title(' Demod I Data');
axis([1 500 -.1 .1]);
pause;
end



sample_i_data=reshape(i_data(1+500:4000+500),20,200);
sample_q_data=reshape(q_data(1+500:4000+500),20,200);

for ii=1:20
   
   i_out=sample_i_data(ii,:);
   q_out=sample_q_data(ii,:);
   
   iq=i_out+1i*q_out;
   
   plot(iq,'*');
   grid;
   %pause
end



end  %phase