随机信号处理基础 - matlab仿真

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);