impulse_duopule.m

来自「脉冲多普勒调制 用多普勒雷达」· M 代码 · 共 55 行

clc;
clear all;
close all;
fs=2e8;%采样率
fc=1e9;%载波
f_middle=30e6;
R=8000;v=400;c=3e8;
Tm=80e-6;%调制周期
T0=15e-6;%脉宽
skip_counter=1;%第一次采样和第二次采样之间差
period_counter=64;%表示经历的周期的数目
t=0:1/fs:period_counter*Tm;%
t_tm=0:1/fs:Tm;
counter_tm=size(t_tm);
%R_new=R-v*t;
t1=1/fs:1/fs:T0;
counter_t1=size(t1);
s=cos(2*pi*f_middle*t1);
%t_delay=2*R_new/c;
receive=zeros(period_counter,3*T0*fs);
all=0;
for i=1:period_counter
   R_new=R-v*(i-1)*Tm;
   t_delay=2*R_new/c;
%(i-1)*Tm++(i-1)*Tm
   t_tmp=t_delay+1/fs:1/fs:t_delay+T0;
   delay_array=zeros(1,fix(t_delay*fs));
   t11=1/fs:1/fs:T0;
   signal_delay=cos(2*pi*f_middle*t11-2*pi*fc*t_delay);
    receive(i,1:3*T0*fs)=[zeros(1,T0*fs) signal_delay zeros(1,T0*fs)];

end;


relative_array=zeros(3*T0*fs,period_counter);
for i=1:period_counter
    for j=1:2*T0*fs
relative_array(j,i)=s*receive(i,j+1:j+T0*fs)';
    end;
end;

for jj=1:3*T0*fs
relative_array(jj,:)=fft(relative_array(jj,1:period_counter),period_counter);

end;
figure;
mesh(abs(relative_array));