自相关和周期自相关函数matlab

clc;clear all;close all;
%% ‘Barker 13 elements‘
% r=10;% signal is sampled r times in each tb （r/tb sampling rate）
% u_amp=ones（113）;% signal amplitude vector
% u_phase=[0 0 0 0 0 1 1 0 0 1 0 1 0];% signal phase vector
% u = [ones（1r*5）-（ones（1r*2））ones（1r*2）-ones（1r）ones（1r）-ones（1r）ones（1r）];%13 bit
% df=.0308;
% K=50;% number of points in the Doppler axis for ambiguity plot is K
%% ‘FrankP1P2 M elements‘
% r=10;% signal is sampled r times in each tb （r/tb sampling rate）
% M=8;
% u_amp=ones（1M*M）;% signal amplitude vector
% for i=1:M
%     for j=1:M
% %         phi（ji）=2*pi/M*（i-1）*（j-1）;% ‘Frank M elements‘
%         phi（ji）= （-pi/M/2）*（2*i-1-M）*（2*j-1-M）;% ‘P2 M elements‘
% %         phi（ji）= （-pi/M）*（M-（2*j-1））*（（j-1）*M+（i-1））;% ‘P1 M elements‘
%     end
% end
% nn=0;
% for ii=1:M
%     for jj=1:M
%         nn=nn+1;
%         phi2（nn）=phi（iijj）;
%     end
% end
% j=sqrt（-1）;
% u_phase=phi2;
% u_basic=u_amp.*exp（j*u_phase）;
% m_basic=length（u_basic）;
% ud=diag（u_basic）;
% ao=ones（rm_basic）;
% m=m_basic*r;
% u=reshape（ao*ud1m）;    % u_basic with each eleement repeated r times
% df=.0312;% maximal displayed Doppler is F/Mtb
K=50;% number of points in the Doppler axis for ambiguity plot is K
%% ‘P3P4 M elements‘
% r=10;% signal is sampled r times in each tb （r/tb sampling rate）
% M=64;
% u_amp=ones（1M）;% signal amplitude vector
% for i=1:M
% %     phi（i）= （pi/M）*（i-1）*（i-1）;% ‘P3 M elements‘
%     phi（i）= （pi/M）*（i-1）*（i-1）-pi*（i-1）;%‘P4 M elements‘
% end
% j=sqrt（-1）;
% u_phase=phi;
% u_basic=u_amp.*exp（j*u_phase）;
% m_basic=length（u_basic）;
% ud=diag（u_basic）;
% ao=ones（rm_basic）;
% m=m_basic*r;
% u=reshape（ao*ud1m）;    % u_basic with each eleement repeated r times
% df=.0312;% maximal displayed Doppler is F/Mtb
% K=50;% number of points in the Doppler axis for ambiguity plot is K
%% ‘T1T2T3T4‘
r=1;
m=2;
k=4;
T=16e-3;
ttt=1;
deltaf=1e3;
for tt = linspace（0T120）;
    jj = floor（k*tt/T）;
%     phase（ttt）= mod（（（2*pi/m）*floor（（（k*tt - jj*T）*（jj*m/T）））） 2*pi）; % T1
%     phase（tt