【雷达原理】雷达杂波抑制方法

Zhi Zhao

已于 2024-02-15 20:42:35 修改

阅读量2k

点赞数 5

分类专栏： # 调频连续波雷达文章标签：雷达原理对消滤波器

于 2023-11-04 22:46:15 首次发布

本文链接：https://blog.csdn.net/weixin_45317919/article/details/133966972

版权

调频连续波雷达专栏收录该内容

6 篇文章 4 订阅

订阅专栏

一、杂波及其特点

1.1 什么是杂波？

杂波是相对目标回波而言的，它是指对雷达正常工作时的目标回波有干扰的回波信号。对于工作在地面的雷达，杂波主要来源于固定建筑物、围栏、树木等静止目标或慢速运动目标产生的回波以及地表反射的回波。

由于单个杂波分量具有随机的相位和幅度，因此合成的总杂波也是随机的，具有类似热噪声的统计特性，并且杂波信号的强度要远大于接收机内热噪声的强度。因此，雷达在强杂波背景下检测目标的能力主要取决于信杂比，而不是信噪比。

杂波抑制的目的就是采用数字信号处理技术滤除或者抑制固定目标杂波和慢速目标杂波，从而保留运动目标回波。

1.2 杂波的频谱特性

雷达工作时可能会碰到的杂波包括地物、海浪、云雨及敌人施放的箔条等。除了孤立的建筑物可认为是固定的点目标外，大多数杂波均属于分布杂波且包含内部运动。当天线不扫描时，固定杂波的功率谱是位于 $nf_{r}$ (n=0，±1，±2，...) 位置上的离散谱线，其中 $f_{r}$ 为脉冲重复频率。当天线扫描时，由于回波脉冲数有限，将引起谱线的展宽。此外，由于天线扫描引起双程天线方向图对回波信号调幅，也会导致功率谱展宽。

一般杂波谱可用高斯函数表示为：

$G(f)=G_{0}e^{-f^2/2\sigma _{s}^2}$

其中， $\sigma _{s}=0.265/n$ ，n为单程天线方向图3dB宽度内的脉冲数。

二、动目标显示(MTI)技术

2.1 对消原理

当固定目标、地杂波等与运动目标处于同一距离单元时，前者的回波通常较强，以至于运动目标的回波被淹没其中，故必须设法对二者进行区分。

因固定目标回波的多普勒频率为零，慢速运动的杂波中所含的多普勒频移也集中在零频附近，它们的回波经相位检波后，输出信号的相位将不随时间变化或者不随时间做缓慢变化，反映在幅度上则为其幅度不随时间变化或者不随时间做缓慢变化。相反运动目标回波经相检输出后，因其相位随时间变化较大，反映在幅度上也是幅度随时间变化较快。

因此若将同一距离单元在相邻重复周期内的脉冲做相减运算，则固定目标回波将被完全对消，慢速杂波也将得到很大程度的衰减，只有运动目标回波得以保留。这样就可以将固定目标、慢速目标与运动目标区别开来。

2.2 数字对消器设计

数字对消器本质上就是一种数字滤波器，根据所设计滤波器的幅频响应特性，可以滤除相应频率的信号分量。常用的MTI滤波器为延迟线对消器，根据对消次数的不同，可分为单延迟线对消器、双延迟线对消器和多延迟线对消器。

1、单延迟线对消器

单延迟线对消器的结构如图1-1所示，滤波器的输出为输入的两个脉冲进行运算得到，该滤波器也称为双脉冲对消器或者一次对消器，其用数学公式可表示为：

$y(t)=x(t)-x(t-T)$

其中，T为脉冲重复周期PRT。

对上式求Z变换，可得： $Y(z)=X(z)-z^{-1}X(z)$

由此可知其系统函数为： $H(z)=1-z^{-1}$

由傅里叶变换与Z变换的关系， $z=e^{jwT}=e^{j2\pi fT}$ ，得到对消器的功率增益为：

$\left | H(w) \right |^2=H(w)H^{*}(w)=(1-e^{-jwT})(1-e^{jwT})=4(sin(\pi fT))^2$

进而得到其幅频响应为：

$\left | H(w) \right |=2\left | sin(\pi fT) \right |$

单延迟线对消器的幅频响应曲线如图1-2所示，取脉冲重复周期 $T=10^{-3}$ ，则脉冲重复频率 $f_{r}=1000$ ，由滤波器的幅频特性曲线可知，在频率 $f=nf_{r}(n=0,\pm 1,\pm 2,...)$ 处有零点，因而该滤波器起到了抑制固定目标回波和慢速运动目标回波的作用。

2、双延迟线对消器

双延迟线对消器的结构如图2-1所示，它由两个单延迟线对消器级联而成。双延迟线对消器经常被称为三脉冲对消器或者二次对消器，其数学表达式如下：

$y(t)=x(t)-2*x(t-T)+x(t-2T)$

对上式求Z变换，可得： $Y(z)=X(z)-2z^{-1}X(z)+z^{-2}X(z)$

由此可知其系统函数为： $H(z)=1-2z^{-1}+z^{-2}=(1-z^{-1})^2$ ，从系统函数的表达式可以看出其等效于两个一次对消器级联。

由此得到二次对消器的功率增益为：

$\left | H(w) \right |^2=\left | H_{1}(w) \right |^2\cdot \left | H_{1}(w) \right |^2=16(sin(\pi fT))^4$

式中， $\left | H_{1}(w) \right |^2$ 为单延迟线对消器的功率增益。

进而得到其幅频响应为：

$\left | H(w) \right |=4\left | sin(\pi fT) \right |^2$

双延迟线对消器的幅频响应曲线如图2-2所示，可以看出，双延迟线对消器比单延迟线对消器具有更好的响应(更深的凹口和更平坦的通带响应)。单延迟线对消器的频率响应较差，原因在于其阻带没有宽的凹口。而双延迟线对消器无论在阻带还是通带上都比单延迟线对消器有更好的频率响应，因此得到了更广泛的应用。

3、多延迟线对消器

多延迟线对消器的结构如图3-1所示，它是由多个单延迟线对消器级联而成，N延迟线对消器的数学表达式为：

$y(t)=\sum_{n=0}^{N}w_{n}x(t-nT)$

式中，N为对消器的次数，对消器的系数 $w_{n}$ 为二项式系数，用下式计算：

$w_{n}=(-1)^nC_{N}^{n}=(-1)^n\frac{N!}{(N-n)!n!},n=0,1,...,N$

N次对消器的系统函数可表示为：

$H(z)=(1-z^{-1})^N=\sum_{n=0}^{N}w_{n}z^{-n}$

其幅频响应和相频响应分别为：

$\varphi (w)=N(\frac{\pi }{2}-wT/2)=N(\frac{\pi }{2}-\pi fT)$

可见，相位响应φ(w)与w是线性关系，所以对消器是一种线性相位滤波器，回波信号通过它后，相位关系不产生非线性变化。

三、MATLAB仿真

3.1 对消效果验证

由第二节的内容可知，MTI技术实质上就是设计一个滤波器，通过对多组脉冲回波的同一个距离单元加权求和即可实现。对于静止目标在该距离单元频谱固定，加权求和后将固定的频谱相互抵消；对于运动目标在该距离单元频谱存在多普勒频移，加权求和后目标频谱不会完全抵消。

通过MATLAB仿真多组具有静止目标和运动目标的回波信号，对比不同次数的对消器对杂波的抑制效果。
案例及对应的代码会在博文【雷达测速原理及实现方法】写完后更新。

2024年2月15日更新：根据最新建立的目标回波模型，添加对消效果的验证

首先产生256组线性调频连续波信号，设置3个目标，其距离分别为500米，1000米，2000米，速度分别为0m/s，8m/s，-5m/s，目标回波信号的幅度均为1，利用上述目标参数产生模拟回波信号，并添加信噪比为10dB的高斯白噪声作为信号的底噪，对该信号先在距离维上进行FFT，然后进行对消处理，效果如图4-1所示。

由图4-1可知，一维FFT结果有3个尖峰，分别代表了距离为500米、1000米、2000米的3个目标，经过一次对消器器处理后，在500米处的静止目标被抑制掉。

对比2脉冲对消器和3脉冲对消器对静止目标的抑制效果，如图4-2所示。

2脉冲对消器和3脉冲对消器均能抑制静止目标的回波信号，对模拟目标回波信号中静止目标的抑制效果极好，可以对消干净。

3.2 代码

clc;
clear;
close all;

%% 单延迟线对消器
PRT = 1e-3;
PRF = 1/PRT;
f = (-2000:1:2000);
H1 = 2*abs(sin(pi*f*PRT));
H1 = H1/max(H1);
H1_d = mag2db(H1);
figure;
subplot(2,1,1);plot(f/PRF,H1);
xlabel('归一化频率f/fr');ylabel('幅度');title('延迟线对消器的幅频响应');
subplot(2,1,2);plot(f/PRF,H1_d);
xlabel('归一化频率f/fr');ylabel('幅度(dB)');ylim([-50 0]);

%% 双延迟线对消器
H2 = 4*abs(sin(pi*f*PRT)).^2;
H2 = H2/max(H2);
H2_d = mag2db(H2);
figure;plot(f/PRF,H1_d);hold on;plot(f/PRF,H2_d);ylim([-50 0]);
xlabel('归一化频率f/fr');ylabel('幅度(dB)');title('延迟线对消器的幅频响应');
legend('单延迟线对消器','双延迟线对消器');

clc;
clear;
close all;

%% LFM信号参数
B = 25e6;               % 带宽
PRT = 800e-6;           % 脉冲重复周期
PRF = 1/PRT;            % 脉冲重复频率

T_chirp = 600e-6;       % 脉冲宽度
u = B/T_chirp;          % 调频斜率
fs = 3e6;               % 采样率
NumADC = round(T_chirp*fs);         % 单个脉冲的采样点数
NumChirp = 256;                     % 脉冲数
c = physconst('LightSpeed');        % 光速
f0 = 10e9;                          % 载频
Lambda = c/f0;                      % 波长

R_max = fs*c/(2*u);                 % 最大测量距离
V_max = Lambda*PRF/4;               % 最大测量速度

%% 目标参数
tarNum = 3;
tar_R0 = [500,1000,2000];           % 目标距离
tar_V0 = [0,8,-5];                  % 目标速度
Rcs = [0.1,1,5];
SNR = 10;

%% 模拟信号
Phase_t = @(f0,u,t) 2*pi*(f0*t+1/2*u*t.^2);        % LFM信号的相位表达式

Tt = (0:1/fs:(T_chirp-1/fs))';
Signal_Tx = exp(1j*Phase_t(f0,u,Tt));           % 发射信号
Signal_Rx = zeros(NumADC,NumChirp);             % 多脉冲的目标回波信号
for jj = 1:NumChirp
    Signal_temp = 0;
    Rt = Tt+(jj-1)*PRT;
    for kk = 1:tarNum
        tar_R = tar_R0(kk)+tar_V0(kk)*Rt;
        tao = 2*tar_R/c;                % 目标回波的时延
        Ar = 1;                         % 目标回波幅度
        Signal_temp = Signal_temp + Ar.*exp(-1j*Phase_t(f0,u,(Tt-tao)));        % 接收信号
    end
    Signal_Rx(:,jj) = awgn(Signal_temp,SNR,'measured');
end

% 混频
Signal_Mix = Signal_Tx.*Signal_Rx;
rawData = Signal_Mix;
clear Signal_Tx;
clear Signal_Rx;
clear Signal_Mix;

%% 距离维FFT
Nfft1 = 2^ceil(log2(NumADC));               % FFT点数
R_point = (fs/Nfft1)*c/(2*u);               % 距离点精度
delta_R = c/(2*B);                          % 距离分辨率

win1 = hamming(NumADC);       % 加窗
fft_Data = zeros(Nfft1,NumChirp);
for ii = 1:NumChirp
    fft_Data(:,ii) = fft(rawData(:,ii).*win1,Nfft1);
end

figure(101);
mesh(mag2db(abs(fft_Data)));xlabel('脉冲数');ylabel('距离维');title('1维FFT结果');
view([90 0]);

cancel_Data = zeros(Nfft1,NumChirp);
for ii = 2:NumChirp
    cancel_Data(:,ii) = fft_Data(:,ii)-fft_Data(:,ii-1);
end

figure(102);
mesh(mag2db(abs(cancel_Data)));xlabel('脉冲数');ylabel('距离维');title('2脉冲对消结果');
view([90 0]);

figure;
subplot(2,1,1);plot(mag2db(abs(fft_Data(:,10))));xlabel('距离维');ylabel('幅度');title('1维FFT结果');
subplot(2,1,2);plot(mag2db(abs(cancel_Data(:,10))));xlabel('距离维');ylabel('幅度');title('2脉冲对消结果');

%% 对消效果验证
% 2脉冲对消
H1 = [1,-1];
Cancel_Data1 = Cancellation(fft_Data,H1);

% 3脉冲对消
H2 = [1,-2,1];
Cancel_Data2 = Cancellation(fft_Data,H2);

figure;
subplot(2,1,1);plot(mag2db(abs(Cancel_Data1(:,10))));xlabel('距离维');ylabel('幅度');title('2脉冲对消结果');
subplot(2,1,2);plot(mag2db(abs(Cancel_Data2(:,10))));xlabel('距离维');ylabel('幅度');title('3脉冲对消结果');

对消函数Cancellation可以去下面的链接下载：

基于MTI滤波器的LFMCW雷达杂波抑制实现方法