clear all;
clc;
%%%%%%%%%%%首先模拟被检测信号产生待采样信号，在本次仿真中，设定被检测信号为BPSK调制，fc=1Mhz
%%%%%%%%%%%基带信号长度为100bit，信号比特率fb=1Mbps，采样频率fs=100Msps
signal_base=randint(1,100,2);
signal_base_bpsk=2*signal_base-1;
T=1/1000000;
t0=0:(T/100):(99*T/100);
for i=1:100,
    if signal_base_bpsk(i)>0
        signal_base_pure(((i-1)*100+1):((i)*100))=sin(2*pi*1000000*t0);
    else
        signal_base_pure(((i-1)*100+1):((i)*100))=-sin(2*pi*1000000*t0);
    end;
end;
signal_base_awgn = awgn(signal_base_pure,0);
%%%%%%%%%将在fs为100Msps下采集得到的数据，进行三倍内插%%%%%%%%%%%%%%

t1=0:(T/100):(100*T-T/100);
s1=signal_base_awgn ;
t2=0:(T/300):(100*T-T/100);
s2=interp1(t1,s1,t2);%%%%3倍线性内插

%%%%%%%%%%分别对内插后数据序列s2，延迟100Ts数据序列s3，延迟250Ts数据序列s4，

s3_temp=[zeros(1,100) s2];
s3=s3_temp(1:2048);
s4_temp=[zeros(1,250) s2];
s4=s4_temp(1:2048);
s6_temp=[zeros(1,120) s2];
s6=s6_temp(1:2048);
s7_temp=[zeros(1,300) s2];
s7=s7_temp(1:2048);
s5=s2(1:2048);
%%%%%%%%%%分别计算s5与s3、s4之间的互相关函数
N=2048;
r35=zeros(1,500);
for m=1:500,
    for n=1:N-m,
   r35(m)=r35(m)+s5(n)*s3(n+m);
    end;
end;
r35=r35/N;
r45=zeros(1,500);
for m=1:500,
    for n=1:N-m,
   r45(m)=r45(m)+s5(n)*s4(n+m);
    end;
end;
r45=r45/N;

r65=zeros(1,500);
for m=1:500,
    for n=1:N-m,
   r65(m)=r65(m)+s5(n)*s6(n+m);
    end;
end;
r65=r65/N;

r75=zeros(1,500);
for m=1:500,
    for n=1:N-m,
   r75(m)=r75(m)+s5(n)*s7(n+m);
    end;
end;
r75=r75/N;