sinsignal.m

对低频段的带通信号进行高精度的时延估计

%用sin信号测试MCZT和FICP求时延                   
M=1;                   %循环次数
TimeDelay=zeros(1,M);
for k=1:M
%Ts=1e-7; Tdelay=265e-9;N1=4e+3;N2=400e+3;
%Ts=0.5e-7; Tdelay=265e-9;N1=16e+3;N2=800e+3;
Ts=0.25e-7; Tdelay=265e-9;N1=64e+3;N2=1.6e+6;
t=Ts:Ts:1.024e-4;
f1=1.44e+6;f2=1.44e+6;
x1n=sin(2*pi*f1*t);
%----------------------------------------
t2=t-Tdelay;
  for n=1:6
    if t2(n)<=0     % 信号的延迟
    t2(n)=0;
    end
   end
x2n=sin(2*pi*f2*t2); 
%----------------------------------------------
x1n=awgn(x1n,10,'measured');  %awgn：Add white Gaussian noise to a signal
x2n=awgn(x2n,10,'measured');  %加高斯白噪声8dB
%--------------------------------------------
N=length(x1n);
X1k=MCZT(x1n,N,N1); plot(abs(X1k));
X2k=MCZT(x2n,N,N1);
Rn=FICP( X1k,X2k,N2);  %快速FICP算法
Nn=find(Rn==max(Rn(N:N+300)))-N-1;
TimeDelay(k)=Nn*N1*Ts/N2;
end
plot(TimeDelay,'.');ylim([258e-9,272e-9]); grid on;