LFM调频连续波雷达的MATLAB仿真

%%本振信号 
%------------------------------------------------------------------------

T=1e-4;                                 %调频周期
B=30e6;                                   %调频带宽 30MHz 
K=B/T;                                      %调频斜率 
Fs=2.5*B;Ts=1/Fs;                      %采样频率 
N=T/Ts;
t0=linspace（-T/2T/2N）;
St0=exp（j*pi*K*t0.^2）;
f0=linspace（-Fs/2Fs/2N）;
%plot（t0*1e4real（St0））; 
%xlabel（‘Time in u sec‘）; 
%title（‘Real part of chirp signal‘）; 
%grid on;axis tight;
%figure（2）;
%plot（f0abs（fftshift（fft（St0length（f0）））））

%-------------------------------------------------------------------------

%%目标设定
%-------------------------------------------------------------------------

c=3e8;            %光速
mm=[306090];     %三个目标的距离
tao=2.*mm./c;      %时间延迟
rmax=100;
BSt=zeros（1N）;
for i=1:3
    temp=[];
    k0=ceil（tao（i）/Ts）;
    temp（k0:k0+N-1）=exp（j*pi*K*t0.^2）;
    temp1=temp（1:N）;
    BSt=temp1+BSt;
end
%plot（t0*1e4real（BSt）‘r‘）;

%%加入杂波干扰项
noise=randn（1N）;    %设信噪比0.2
BSt=BSt+noise;
%%混频滤波获取差拍信号
%-------------------------------------------------------------------------
u=real（BSt）.*real（St0）;
%f=linspace（-Fs/2Fs/2N）;
plo