%?input（‘\nPulse?radar?compression?processing:?\n?‘）;
clear;
close?all;


????T=10e-6;????????????????????????????????????????%pulse?duration?10us
????B=30e6;?????????????????????????????????????????%chirp?frequency?modulation?bandwidth?30MHz
????Rmin=8500;Rmax=11500;??????????????????????????%range?bin
????R=[90001000010200];????????%position?of?ideal?point?targets
????%R=[12000?];
????RCS=[1?1?1??];??????????????????????????????%radar?cross?section
????%RCS=[1?];
%Parameter
C=3e8;??????????????????????????????????????????????%propagation?speed
K=B/T;??????????????????????????????????????????????%chirp?slope
Rwid=Rmax-Rmin;?????????????????????????????????????%receive?window?in?meter
Twid=2*Rwid/C;??????????????????????????????????????%receive?window?in?second
Fs=5*B;Ts=1/Fs;?????????????????????????????????????%sampling?frequency?and?sampling?spacing
Nwid=ceil（Twid/Ts）;????????????????????????????????%receive?window?in?number

%Gnerate?the?echo??????
t=linspace（2*Rmin/C2*Rmax/CNwid）;?????????????????%receive?window
??????????????????????????????????????????????????%open?window?when?t=2*Rmin/C
??????????????????????????????????????????????????%close?window?when?t=2*Rmax/C????????????????????????????
M=length（R）;???????????????????????????????????????%number?of?targets????????????????????????????????????????
td=ones（M1）*t-2*R‘/C*ones（1Nwid）;

SNR=[11010010000];

for?i=1:1:4
Srt1=RCS*（exp（1i*pi*K*td.^2）.*（abs（td）n=sqrt（0.5*SNR（i））*（randn（size（Srt1））+1i*randn（size（Srt1）））;
Srt=Srt1+n;

%Digtal?processing?of?pulse?compression?radar?using?FFT?and?IFFT
Nchirp=ceil（T/Ts）;??????????????????????????????????%pulse?duration?in?number
Nfft=2^nextpow2（Nwid+Nwid-1