ISAR单特显点法相位补偿

%%demo of LFM pulse radar %================================================================== function LFM_radar(T,B,Rmin,Rmax,R,RCS) if nargin==0 T=10e-6; %pulse duration 10us B=30e6; %chirp frequency modulation bandwidth 30MHz Rmin=10000;Rmax=15000; %range bin R=[10500,11000,12000,12008,13000,13002]; %position of ideal point targets RCS=[1 1 1 1 1 1]; %radar cross section end %================================================================== %%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/C,2*Rmax/C,Nwid); %receive window %open window when t=2*Rmin/C %close window when t=2*Rmax/C M=length(R); %number of targets td=ones(M,1)*t-2*R'/C*ones(1,Nwid); Srt=RCS*(exp(j*pi*K*td.^2).*(abs(td) %%demo of LFM pulse radar %================================================================== function LFM_radar(T,B,Rmin,Rmax,R,RCS) if nargin==0 T=10e-6; %pulse duration 10us B=30e6; %chirp frequency modulation bandwidth 30MHz Rmin=10000;Rmax=15000; %range bin R=[10500,11000,12000,12008,13000,13002]; %position of ideal point targets RCS=[1 1 1 1 1 1]; %radar cross section end %================================================================== %%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/C,2*Rmax/C,Nwid); %receive window %open window when t=2*Rmin/C %close window when t=2*Rmax/C M=length(R); %number of targets td=ones(M,1)*t-2*R'/C*ones(1,Nwid); Srt=RCS*(exp(j*pi*K*td.^2).*(abs(td)