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