matlab 线性调频信号匹配滤波

% 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