clear all;
close all;
clc;
%AD采樣率為20倍符號速率；
fs=1228800*20;
%讀取數據，16量化，
datain=read_binary ('PILOTSIG.dat',50000);
% datain=read_binary ('SKYFSIG.dat',50000);
datain=datain';
%中頻采樣后剩余載波頻率為 fc=3.7280e6Hz
fc=3.7280e6;%-100*4;
%fc=3.72755e6;%鎖相環跟蹤后的真實的載波頻率;
length_wave=length(datain);
t=[0:length_wave-1]/fs;
phi0=0;%初始相位
lo_cos=cos(2*pi*fc*t+phi0);
lo_sin=sin(2*pi*fc*t+phi0);
hi_data=datain.*lo_cos';	% down converting
hq_data=datain.*lo_sin';	% down converting

%低通濾波器，
b=fir1(30,0.7e6/fs);
ii_data=filter(b,1,hi_data);
qq_data=filter(b,1,hq_data);

%采樣頻率將為4倍符號速率
fs=fs/5;
i_data=ii_data(100:5:end)';
q_data=qq_data(100:5:end)';
clear datain lo_cos lo_sin t ii_data qq_data hi_data hq_data;

offset_n=3; % 1~4之間
n_ratio=4;

N1=length(i_data)/n_ratio-10;
for kk=1:N1
    baseband_data(kk)=i_data(offset_n+n_ratio*kk)+1i*q_data(offset_n+n_ratio*kk);
end;
figure,
plot(baseband_data,'*');


I_data=real(baseband_data(500:end));
Q_data=imag(baseband_data(500:end));

%讀取短PN序列，該序列已仿真產生并驗證
load PNDATA_align.mat M_PN_I M_PN_Q;


%任意選64個符號數據與短碼序列做滑動相關，求出峰值位置
symbol_shift_pos=700;
c_i_data=I_data(symbol_shift_pos:symbol_shift_pos+63);
c_q_data=Q_data(symbol_shift_pos:symbol_shift_pos+63);

[PeakValue1, PeakPos1 ,corrvalue1] = Pilot_PN_Search(c_i_data,c_q_data,M_PN_I,M_PN_Q);
figure,
plot(corrvalue1);