clear?all;

f0=10e9;
B=300e6;
fs=f0*10;
Tt=1e-7;
tt=linspace（0TtTt*fs）;
k=B/Tt;
u=1*（exp（1j*pi*（2*（f0-B/2）*tt+k*tt.^2）））;%postive?chirp

M=length（u）;

F=8;
K=100;
df=F/K/M;
T=1;
N=100;
uamp=abs（u）;
phas=angle（u）;
t=0:M-1;
tscale1=[0?0:M-1?M-1];
dphas=[NaN?mod（diff（phas）2*pi）]/2/pi;

%?plot?the?signal?parameters
figure（1）?clf?hold?off?
subplot（311）
plot（tscale1[0?abs（uamp）?0]‘linewidth‘1.5）
ylabel（‘?Amplitude?‘）
axis（[-inf?inf?0?1.2*max（abs（uamp））]）

subplot（312）
plot（t?phas‘linewidth‘1.5）
axis（[-inf?inf?-inf?inf]）
ylabel（‘?Phase?[rad]?‘）

subplot（313）
plot（tdphas*ceil（max（t））‘linewidth‘1.5）
axis（[-inf?inf?-inf?inf]）
xlabel（‘?\itt?/?t_b?‘）
ylabel（‘?\itf?*?Mt_b?‘）

dtau=ceil（T*M）/N;
tau=round（[0:1:N]*dtau）;
f=[0:1:K]*df;?
f=[-fliplr（f）?f];
mat1=spdiags（u‘0M+ceil（T*M）M）;
u_padded=[zeros（1ceil（T*M））uzeros（1ceil（T*M））];
cidx=[1:M+ceil（T*M）];
ridx=round（tau）‘;
index?=?cidx（ones（N+11）:）?+?ridx（:ones（1M+ceil（T*M）））;
%?mat2?=?sparse（u_padded（index））;?
mat2?=?（u_padded（index））;?
uu_pos=mat2*mat1;
e=exp（-j*2*pi*f‘*t）;
a_pos=