Tr=100e-6;
Tsub=60e-6;
Bsub=3*10^6;
fs=4*10^6;
fsub=3*10^6;
f0=24*10^6;
M=16;
R=50*10^3;
V=50;
k=Bsub/Tsub;
c=3*10^8;
Nr=fs*Tr;
t1=linspace（-Tr/2Tr/2Nr）;
u=exp（j*pi*k*t1.^2）.*（abs（t1）%?a=[13910135151116148741226];
figure（1）
plot（t1*1e6real（u））;
f=linspace（-fs/2fs/23*Nr+1）;
ufft=（fft（u3*Nr+1））;
figure（2）
plot（f*1e-6abs（fftshift（ufft）））;
%?figure（2）
%?plot（abs（fftshift（ufft）））;
i=0;
x=zeros（3*Nr+1Nr）;
for?freq_del=-Bsub:2*Bsub/（3*Nr）:Bsub;
????i=i+1;
????exf=exp（j*2*pi*freq_del.*t1）;
????exf1=u.*exf;
????ufft1=（fft（exf13*Nr+1））;
????ufft2=ufft1.*conj（ufft）;
????x（:i）=fftshift（abs（ifft（ufft2））.‘）;
end
delay=-Tsub:2*Tsub/（3*Nr）:Tsub;
freq_del=-Bsub:2*Bsub/（3*Nr）:Bsub;
figure（3）
mesh（freq_deldelayx./max（max（x）））;
%?colormap（[.5?.5?