clc;
clear?all;
close?all;
warning?off;
fc=100e6;????????????????????
Tr=11e-6;
fd=10e6;?
td=2e-6;?
Tp=0.6e-6;
fs=400e6;
t=-Tp/2:1/fs:Tp/2;
N=length（t）;
B=40e6;
lfm=exp（j*2*pi*（fc*t+0.5*B/Tp*t.^2））
figure（1）
plot（-Tp/2+Tp/N:Tp/N:Tp/2lfm）;
xlabel（‘脈沖寬度（s）‘）;?
ylabel（‘幅度‘）;
title（‘信號波形‘）;
axis（[-Tp/2Tp/2-22]）;
grid?on;
fft_lfm=fftshift（fft（lfmN））;
figure（2）
subplot（211）plot（-fs/2+fs/