td=ones(M,1)*t-2*R'/C*ones(1,Nwid);
Srt1=RCS*(exp(1i*pi*K*td.^2).*(abs(td) Nchirp=ceil(T/Ts); Srw=fft(Srt,Nfft); Srw1=fft(Srt1,Nfft); t0=linspace(-T/2,T/2,Nchirp); St=exp(1i*pi*K*t0.^2); Sw=fft(St,Nfft); Sot=fftshift(ifft(Srw.*conj(Sw))); Sot1=fftshift(ifft(Srw1.*conj(Sw))); N0=Nfft/2-Nchirp/2; Z=abs(Sot(N0:N0+Nwid-1)); figure subplot(211) plot(t*1e6,real(Srt)); axis tight; xlabel('us');ylabel('幅度') title(['线性调频信号压缩前']); subplot(212) plot(t*C/2,Z) xlabel('Range in meters');ylabel('幅度 ') title(['线性调频信号压缩后']); Nfft=2^nextpow2(Nwid+Nwid-1);
