脉冲多普勒雷达

脉冲多普勒雷达是利用运动的物体对发射信号的产生的多普勒效应来测速和测距的，利用回波的延时可以测距，求出多普勒频率可以求出物体的速度

clc;
clear all;
close all;
fs=2e6;%采样率
fc=1e10;%载波
R1=50000;R2=30000;v1=10;v2=3;c=3e8;
Tm=500e-6;%调制周期
T0=5e-6;%脉宽
period_counter=64;%表示经历的周期的数目
t=0:1/fs:period_counter*Tm;
tm=0:1/fs:Tm;
counter_tm=size（tm）;%一个周期采样点的个数
t0=1/fs:1/fs:T0;%一个脉宽采样点
counter_t0=size（t0）;%一个脉宽采样点数
%发射信号
send=0;
for i=1:period_counter
    send=[send ones（1T0*fs） zeros（1（Tm-T0）*fs）];%
end;
figure（1）;
plot（send）;
axis（[0 1e4 -0.5 2.0]）
title（‘发射信号‘）;xlabel（‘采样点数‘）;ylabel（‘幅度‘）;
%回波信号
%有迟延
receive=0;
t_delay1=2*R1/c;
t_delay2=2*R2/c;
delay_array1=zeros（1fix（t_delay1*fs））;%fix:去尾
delay_array2=zeros（1fix（t_delay2*fs））;
receive1=[delay_array1 send];
receive2=[delay_array2 send zeros（1fix（t_delay1*fs）-fix（t_delay2*fs））];
receive=0.3*receive1+receive2;
figure（2）;
plot（receive）;
axis（[0 1e4 -2.0 2.0]）
title（‘固定目标反射的回波‘）;xlabel（‘采样点数‘）;ylabel（‘幅度‘）;
%有多普勒调制
fd1=200;fd2=800;
tt=0:1/fs:period_counter*Tm+t_delay1;
receive_fd1=receive1.*cos（2*pi*fd1*tt）;
receive_fd2=receive2.*cos（2*pi*fd2*tt）;
receive_fd=0.3*receive_fd1+receive_fd2;
figure（3）;
plot（receive_fd）;
axis（[0 5e4 -1.5 1.5]）
title（‘两个移动目标反射的回波‘）;xlabel（‘采样点数‘）;ylabel（‘幅度‘）;

%匹配滤波
cankao=[ones（1T0*fs） zeros（1（Tm-T0）*fs）];
cankao=fliplr（cankao）;
MF=conv（cankaoreceive_fd）;
figure（4）;
plot（MF）;
axis（[0 5e4 -32 32]）
title（‘匹配滤波结果‘）;xlabel（‘采样点数‘）;ylabel（‘幅度‘）;

%数据重排
zz=1000-mod（length（MF）1000）;
MF1=[MF zeros（1zz）];
for m=1:1:（length（MF）/1000）
    for n=1:1:1000
        speed（mn）=MF1（11000*（m-1）+n）;
    end
end
real_speed=real（speed）;
figure（5）;
meshc（real_speed）
title（‘数据重排三维图‘）;xlabel（‘距离门‘）;ylabel（‘速度门‘）;
speed_out=fft（speed）;
figure（6）;
meshc（abs（speed_out））
title（‘FFT处理‘）;xlabel（‘距离门‘）;ylabel（‘速度门‘）;

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

     文件       1882  2010-01-28 19:44  pd_duopule.m

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

                 1882                    1