雷达发射线性调频信号的模糊函数图以及等高线图，速度模糊，时延模糊图

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