匹配滤波处理不同距离维回波信号

• 可以参考之前写的两篇博客
  • 线性调频信号
  • 原理说明与仿真

(1) 原理说明

• 假设目标与雷达的相对距离为 $R$，雷达发射信号 $s(t)$ ,传播速度为光速 $c$ ，则经过时间 $\frac{R}{c}$ 后电磁波到达目标，照射到目标上的电磁波可写成：$s(t-\frac{R}{c})$。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为 $\sigma \cdot s(t-\frac{R}{c})$，其中$\sigma$ 为目标的雷达散射截面。再经过时间$\frac{R}{c}$后，被雷达接收天线接收的信号为$\sigma \cdot s(t-2\frac{R}{c})$。
   

• 如果将雷达天线和目标看作一个系统，便得到如图的等效，而且这是一个 LTI（线性时不变）系统。

  

  等效 LTI 系统的冲激响应可写成：
  $$h(t)=\sum _{i=1}^M{\sigma }_i\delta (t-{\tau }_i)$$

  其中 $M$ 表示目标个数，${\sigma }_i$ 表示目标的散射特性，${\tau }_i$ 是光速在雷达与目标之间往返一次的时间
  $${\tau }_i=\frac{2R_i}{c}$$
  $R_i$ 为第 $i$ 个目标与雷达的相对距离。雷达发射信号 $s(t)$ 经过该LTI系统，得输出信号（即雷达的回波信号）$sr(t)$
  $$sr(t)=s(t)\ast h(t)=\sum _{i=1}^M{\sigma }_{i}s(t-{\tau }_i)$$

(2) 代码仿真

• matlab 代码

  %% 产生线性调频波
  TimeWidth = 10e-6;      %脉冲持续时间10us
  T = 12e-6;              %选取T在12us
  BandWidth = 20e6;       %线性调频信号的频带宽度20MHz
  Fs = 50e6;Ts = 1/Fs;    %采样频率，注意需要满足奈奎斯特频率
  sample_dot_num = round(TimeWidth * Fs);%表示采样点的个数
  
  f0 = 0;%初始频率
  f1 = f0 + BandWidth;%终止频率
  t=0:1/Fs:TimeWidth;%根据结束时间生成时间序列
  signal = chirp_signal(t,f0,f1);
  plot(t*1e6,signal,'LineWidth',2);
  %% 进行传播并接收回波
  c = 3e5;  %光速单位Km/s
  R1 = 80;  %第一个物体80Km
  R2 = 120; %第二个物体120Km
  R3 = 180; %第三个物体180Km
  %计算时间延迟
  tau1 = 2*R1/c;
  tau2 = 2*R2/c;
  tau3 = 2*R3/c;
  
  %产生回波信号
  Nmax = round(1.1*tau3*Fs);
  signal_total = zeros(1, Nmax);
  signal_total(fix(tau1*Fs):fix(tau1*Fs)+length(signal)-1) = signal;
  signal_total(fix(tau2*Fs):fix(tau2*Fs)+length(signal)-1) = signal;
  signal_total(fix(tau3*Fs):fix(tau3*Fs)+length(signal)-1) = signal;
  snr = -20; %添加信噪比为-20dB
  signal_twn = awgn(signal_total,snr, 'measured');
  plot((0:Nmax-1)*Ts*1e3,signal_twn);
  
  mf = fliplr(signal);                    %输入序列的复共轭（匹配滤波器的脉冲响应）
  sr_noise_mf = abs(conv(mf,signal_twn)); %加入噪声后的卷积反馈
  
  ty = (0:length(signal_twn)-length(signal)-1)*Ts*1e3; %注意时间这里要减去因为卷积增加的一部分
  plot(ty,sr_noise_mf(1:length(ty)))
  plot(ty*c/2/1000,sr_noise_mf(1:length(ty)))
  

• 仿真结果：