多普勒脉冲雷达回波仿真

多普勒脉冲雷达回波仿真。产生回波，对回波进行距离压缩，进行两脉冲对消，观察运动、盲速与静止目标的对消情况。

clear all;clc;close all;
fc=3e9;                 %载波频率
PRF=2000;       
Br=5e6;                 %带宽
fs=10*Br;               %采样频率
Tp=5e-6;                %脉宽
Kr=Br/Tp;               %频率变化率
c=3e8;                  %光速
lamda=c/fc;             %波长
Tr=1/PRF;               %脉冲重复周期
N_mc=1.5/60*PRF;        %脉冲个数
t=0:1/fs:15*Tp+Tp;      %采样时间
N_r=length（t）;          %采样点数
N_target=5;             %目标个数
Rmax=c/2*15*Tp;                             %目标最大距离
R_t=Rmax*abs（rand（1N_target））;             %目标的距离
RCS_t=10*（exp（i*2*pi*rand（1N_target）））;    %目标RCS，幅度为10，相位在（02pi）之间随机分布
Vmax=lamda*PRF/2;                           %目标最大速度
v=Vmax*（（1+rand（1N_target））/2）;            %目标速度
%% 生成目标矩阵
sr=zeros（N_mcN_r）;
for i=1:N_mc
    ta=（i-1）*Tr;
    sri=0;
    for k=1:N_target
        tao=2*（R_t（k）-v（k）.*（ta+t））/c;
        srj=RCS_t（k）.*rectpuls（t-tao-Tp/2Tp）.*exp（-1j*2*pi*fc*tao+1j*pi*Kr.*（t-tao-Tp/2）.^2）;
        sri=sri+srj;
    end
    sr（i:）=sri;
end
%% 距离压缩前的回波
tm=（1:N_mc）/PRF;
R=c*t/2;
figure（1）;
image（Rtm255*abs（sr）/max（max（abs（sr））））  
figure（2）;
plot（t*c/2abs（sr（1:）））         
%% 距离压缩
st=rectpuls（t-Tp/2Tp）.*exp（1i*pi*Kr*（t-Tp/2）.^2）;
stf=conj（fft（st））;
for i=1:N_mc
    sr（i:）=ifft（fft（sr（i:））.*stf）;            
end
figure（3）;
image（Rtm255*abs（sr）/max（max（abs（sr））））                
figure;
plot（t*c/2abs（sr（1:）））                       
             
sr=fft（sr[]1）;
V=linspace（0PRF50）*lamda/2;
figure;image（RV255*abs（sr）/max（max（abs（sr））））  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1624  2017-05-22 11:46  echo.m

----------- ---------  ---------- -----  ----

                 1624                    1