雷达信号处理——恒虚警检测（CFAR）_2dcfar 翻转法

雷达信号处理的流程

雷达信号处理的一般流程：ADC数据——1D-FFT——2D-FFT——CFAR检测——测距、测速、测角。

雷达目标检测

首先要搞清楚什么是检测，检测就是判断有无。雷达在探测的时候，会出现很多峰值，这些峰值有可能是目标、噪声和干扰，那么就需要进行一种判决，来决定它是否属于目标。
雷达目标检测问题是利用信号的统计特性和噪声的统计特性解决在受噪声干扰的观测中，信号有无的判决问题。其数学基础就是统计判决理论，又称假设检验理论。假设检验是进行统计判决的重要工具，信号检测相当于数理统计中的假设检验。

恒虚警检测概念

假设就是检验对象的可能情况或状态。对于雷达检测来说，可以选用两个假设，即目标存在或不存在。用$H_0$和$H_1$分别表示接收机的无信号输入和有信号输入的假设，用$D_0$和$D_1$和分别表示检测器作出无信号和有信号的判决结果。于是接收机的输入与检测器的判决将有四种情况：
(1) $H_0$为真，判为 $D_0$，即接收机无信号输入, 检测器判为无信号, 其概率称为正确不发现概率$P_{an}$；
(2) $H_0$为真, 判为 $D_1$， 即接收机无信号输入, 检测器判为有信号, 其概率称为虚警概率 $P_{fa}$；
(3) $H_1$为真, 判为 $D_0$，即接收机有信号输入, 检测器判为无信号, 其概率称为漏警概率 $P_{la}$；
(4) $H_1$为真, 判为 $D_1$，即接收机有信号输入, 检测器判为有信号, 其概率称为正确检测概率 $P_d$；
我们希望的是正确不发现概率$P_{an}$和真确检测概率 $P_d$越高越好，而虚警概率 $P_{fa}$和漏警概率 $P_{la}$越低越好。
恒虚警检测，又称CFAR，全称是Constant False Alarm Rate Detector，是雷达目标检测的一种常见的手段。CFAR检测算法属于信号检测中的自动检测算法，在雷达信号处理中主要应用的有四种，即CA-CFAR、SO-CFAR、GO-CFAR、WCA-CFAR。

恒虚警检测算法

首先，指定一个距离多普勒单元，称为检测单元(CUT, cell under test)，为了判断检测单元处是否存在目标，需要设定一个阈值$T$，当检测单元处的检测值$Y$大于$T$时，判断为有目标，否则判断为没有目标：
d ( Y ) = { 有目标， Y ≥ T 无目标， Y < T . d(Y)= \left\{

$$\begin{array}{ll} 有目标，Y\geq T\\ 无目标，Y< T. \end{array}$$ \right. d(Y)={有目标，Y≥T无目标，Y<T.​
具体的做法是先在测试单元的周围放置保护单元，以消除信号在测试单元周围的多余能量，然后再保护单元的周围选择参考单元：

图1 恒虚警检测的单元分布

上图展示的是同时在距离维和多普勒维进行恒虚警检测的原理，如果只是在距离维或多普勒维做CFAR，检测单元的左右两侧均有设定好数量的保护单元和参考单元。
恒虚警检测选择均值类自适应门限 CFAR 检测技术，其中单元平均方式(CA-CFAR)、两侧单元平均 选大方式(GO-CFAR)、两侧单元平均选小方式 (SO-CFAR) 三种方法在主瓣多目标的背景下，检测性能都有所下降，容易出现漏警和漏检小信号源。单限切除恒虚警概率方法
(E-CFAR)，在主瓣多目标的背景下，具有较好的检 测性能。E-CFAR 方法是在进行单元平均之前设置切 除上限，切除信号的强分量，从而避免对小信号源 漏检，防止对信源数目的估计出现偏差。
如果是CA-CFAR，则取左右两边的参考单元的均值再取平均值，然后与$\alpha$相乘（$\alpha$计算公式如下方公式），再将结果与阈值相比较。如果大于阈值则代表有目标，否则认为无目标。

图2 一维恒虚警检测的阈值计算方法

从上图可以看到，由于待检单元左侧还有几个单元，所以做完恒虚警检测之后会使得序列变短，这个和深度学习中的卷积类似。要想使得不变短的话可以采取补零法或者翻转法来实现。
从图中也可以看到，根据阈值的选择不同，有不同的检测器，其中CA-CFAR是取的参考单元的均值，SO-CFAR取的是参考单元的最小值，而GO-CFAR是取的参考单元的最大值。
其中最重要的是门限因子$\alpha$的计算方法：
$$\alpha =N(P_{fa}^{-\frac{1}{N}}-1)$$

恒虚警检测matlab编程

不同的CFAR检测器适用于不同的场景，CA-CFAR算法适用于均匀环境中的单目标或彼此相距较远的多目标环境，在邻近多目标环境会发生目标的遮蔽效应，在杂波边缘会发生杂波边缘效应。
下面是不进行补零操作的CA-CFAR的程序：

function [position, threshold, start_cell, stop_cell] = cacfar(signal, Pfa, ref_num, guard_num)
% ======>INPUT:
% signal: Data of signal(include signal and noise).[DATATYPE: row vector]
% Pfa: Probability of false alarm.[DATATYPE: scalar]
% ref_num: Number of reference cell.[DATATYPE: scalar]
% guard_num: Number of guard cell.[DATATYPE: scalar]
% ======>OUTPUT:
% position: positions of target.[DATATYPE: row vector]
% threshold: CFAR threshold of input signal.[DATATYPE: row vector]
position = [];
left_num = guard_num + ref_num;
start_cell = left_num + 1;
stop_cell = length(signal) - left_num;
N = 2*ref_num;
alpha = N * (Pfa ^ (-1/N) - 1);
threshold = zeros(1, stop_cell - start_cell + 1);
for ii = start_cell : stop_cell
    tmp_data = [signal(ii-left_num : ii-guard_num-1), ...
        signal(ii+guard_num+1 : ii+left_num)];
    mean(tmp_data);
    tmp = mean(tmp_data) * alpha;
    % threshold(ii - left_num)
    tmp;
    threshold(ii - left_num) = tmp;
    if tmp < signal(ii)
        position = [position, ii];
    end
end
return
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

下面是进行补零操作的CA-CFAR的程序：

function [position, threshold] = MYCACFAR(signal, Pfa, ref_num, guard_num)
% 该函数用于进行CA-FAR恒虚警检测，在信号首位两端进行补零，使得输出的阈值序列和信号长度一致
% ======>INPUT:
% signal: Data of signal(include signal and noise).[DATATYPE: row vector]
% Pfa: Probability of false alarm.[DATATYPE: scalar]
% ref_num: Number of reference cell.[DATATYPE: scalar]
% guard_num: Number of guard cell.[DATATYPE: scalar]
% ======>OUTPUT:
% position: positions of target.[DATATYPE: row vector]
% threshold: CFAR threshold of input signal.[DATATYPE: row vector]
position = []; %存储高于检测门限的信号区域的位置
left_num = guard_num + ref_num; %检测单元两边的单元数
sigLen = length(signal); %信号长度
tempSig = zeros(1,sigLen+2*left_num); %补零后的信号向量
tempSig((left_num + 1) :(sigLen + left_num)) = signal; %对中间部分进行信号的赋值，首位两端各有left_num个零值
N = 2*ref_num; %参考单元总数
alpha = N * (Pfa ^ (-1/N) - 1); %门限因子，横虚警概率越小，参考单元数越多，因子越大，检测阈值也越大
threshold = zeros(1, sigLen); %用于储存检测阈值
for ii = 1 : 1 :sigLen
    centNum = ii+left_num; %对应的检测单元下标
    tmp_data = [tempSig((centNum-left_num):(centNum-guard_num-1)), ...
        tempSig((centNum+guard_num+1):(centNum+left_num))]; %取检测单元的值
    tmp = mean(tmp_data) * alpha; %与门限因子相乘得到检测阈值
    threshold(ii) = tmp; %进行存储
    if tmp < signal(ii) %阈值小于信号
        position = [position, ii];  %存储被检测到的下标
    end
end
return
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

下面是进行翻转补充首位的操作的CA-CFAR的程序：

function [position, threshold] = MYCACFAR2(signal, Pfa, ref_num, guard_num)
% 该函数用于进行CA-FAR恒虚警检测，在信号首位两端进行翻转，使得输出的阈值序列和信号长度一致
% ======>INPUT:
% signal: Data of signal(include signal and noise).[DATATYPE: row vector]
% Pfa: Probability of false alarm.[DATATYPE: scalar]
% ref_num: Number of reference cell.[DATATYPE: scalar]
% guard_num: Number of guard cell.[DATATYPE: scalar]
% ======>OUTPUT:
% position: positions of target.[DATATYPE: row vector]
% threshold: CFAR threshold of input signal.[DATATYPE: row vector]
position = []; %存储高于检测门限的信号区域的位置
left_num = guard_num + ref_num; %检测单元两边的单元数
sigLen = length(signal); %信号长度
tempSig = zeros(1,sigLen+2*left_num); %补零后的信号向量
tempSig(1:left_num) = signal(1:left_num);
tempSig(sigLen + left_num+1:end) = signal(sigLen-left_num+1:end);
tempSig((left_num + 1) :(sigLen + left_num)) = signal; %对中间部分进行信号的赋值，首位两端各有left_num个零值
N = 2*ref_num; %参考单元总数
alpha = N * (Pfa ^ (-1/N) - 1); %门限因子，横虚警概率越小，参考单元数越多，因子越大，检测阈值也越大
threshold = zeros(1, sigLen); %用于储存检测阈值
for ii = 1 : 1 :sigLen
    centNum = ii+left_num; %对应的检测单元下标
    tmp_data = [tempSig((centNum-left_num):(centNum-guard_num-1)), ...
        tempSig((centNum+guard_num+1):(centNum+left_num))]; %取检测单元的值
    tmp = mean(tmp_data) * alpha; %与门限因子相乘得到检测阈值
    threshold(ii) = tmp; %进行存储
    if tmp < signal(ii) %阈值小于信号
        position = [position, ii];  %存储被检测到的下标
    end
end
return
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31

下面是仿真结果：

图2 单目标脉冲压缩结果及CA-CFAR阈值图

对目标峰值部分局部放大后

图3 局部放大

附录

二维CFAR 仿真代码地址：https://github.com/tooth2/2D-CFAR