直升机起降中无线紫外光喷泉码引导方法研究

doi:10.11999/JEIT150004

直升机起降中无线紫外光喷泉码引导方法研究

doi: 10.11999/JEIT150004

基金项目:

国家自然科学基金委员会-中国民航局民航联合研究基金(U1433110)，陕西省科技计划工业公关项目(2014K05-18)，陕西省教育厅产业化培育项目(2013JC09)，西安市科学计划项目(CXY1435(4))和西安市碑林区科技计划项目(GX1302)

计量
- 文章访问数: 1054
- HTML全文浏览量: 65
- PDF下载量: 489
- 被引次数: 14
出版历程
- 收稿日期: 2015-01-05
- 修回日期: 2015-03-31
- 刊出日期: 2015-10-19

Research on Wireless Ultraviolet Guiding Fountain Code Technology of Helicopter Emergency Landing

摘要

摘要: 为了提高直升机在应急无线紫外光通信辅助起降中信标搜寻、定位降落的效率和可靠性，该文提出一种基于喷泉(LT)码的分等级编码通信方案，采用不同等级编码的分级方法设计了不同场景、不同引导阶段的通信策略，分析了高丢包率信道环境中变化删除概率下的误比特率，并与其他传统信道编码进行了比较。实验仿真结果表明，在二进制删除信道下，该分等级编码方法可显著降低LT码的误比特率，提高了编译码性能；该通信策略适用于高丢包率信道环境，能获得更好的通信效果，提高了直升机应急辅助起降的安全性。
- 紫外光通信 /
- 喷泉码 /
- 直升机助降 /
- 引导方法
Abstract: In order to improve the efficiency and the reliability of the helicopter guiding in the beacon researching and locate landing, an encoding communication scheme based on unequal level coding Luby Transform (LT) code is proposed. Different communication strategies are proposed for different scenarios. The different guiding strategies considering the unequal level coding scheme are designed. The bit error rate of the unequal level coding LT code with different communication strategies is analyzed. The coding scheme and communication strategy are simulated by computer. The simulation results show that the bit error rate of the unequal level protection encoding scheme is significantly reduced and the coding performance is improved over binary erasure channel. The strategies can adapt to different channel environments, and the reliability of communication is improved. Thus, helicopter emergency landing safety is improved.
- Ultraviolet communication /
- Fountain code /
- Helicopter landing guiding /
- Guiding method

HTML全文

1. 引言

在外辐射源雷达低空小目标探测中，其背景杂波受传播环境影响而复杂多变，尤其是功率突变的杂波会在距离多普勒谱上形成点状、块状的杂波边缘，极大地影响了低空小目标的恒虚警检测性能。

当前常用的恒虚警算法根据杂波在空间与时序上的分布特性与均匀性可划分为：空域类(也称为滑窗类)和时序类。

空域类恒虚警算法中具有代表性的是单元平均(Cell Averaging)CA-CFAR算法^[1,2]，其采用临近待检测单元的杂波均值来计算检测阈值，在均匀环境下效果较好，但是在多目标及杂波边缘环境下无法消除临近目标与能量突变干扰，使得检测性能恶化。在该算法的基础上，Trunk^[3]提出了最小选择(Smallest Of)SO-CFAR以解决多目标干扰问题，Hansen等人^[4]提出最大选择(Greatest Of)GO-CFAR以解决杂波边缘环境下的能量突变问题，但是这两种算法都在均匀环境下存在一定性能损失，且并不适用于另一种非均匀环境。Smith等人^[5]在以上算法的基础上提出了变指数(Varibility Index)VI-CFAR。该算法先判决当前杂波环境的类型，然后选择与其匹配的其他恒虚警算法^[6,7]，其本质是多种算法的组合。Rohling等人^[8-11]提出有序统计(Order Statistics)OS-CFAR，虽然该算法在多目标、杂波边缘环境下性能得到部分改善，但在均匀环境下性能也有一定损失。

时序类恒虚警算法^[12-16]中，杂波图(Clutter Map)CM-CFAR是该类算法的典型代表，CM-CFAR将杂波空间划分为杂波图单元用以存储杂波数据，每个杂波图单元依靠当前数据和历史数据进行加权迭代更新。由于这类算法在迭代初期采样的数据量少，所以检测性能较差；随着迭代次数增加，检测性能会得到较大的提升。但是这类算法也具有因迭代计算而带来的记忆效应这一不可忽视的缺点。

在外辐射源雷达低空小目标探测中，由于目标回波弱、背景杂波复杂多变，空域类算法在该类应用中性能损失严重、虚警概率恶化，传统时序类算法迭代初期性能较差并且还有记忆效应等问题，为此本文提出了动态有序矩阵(Dynamic Order Matrix)DOM-CFAR算法，其将杂波空间划分成有序矩阵，利用雷达初次扫描数据的滑窗中值对有序矩阵中的有序数组进行初始化用以弥补传统时序类算法在迭代初期的性能劣势，并通过雷达扫描数据对有序数组进行动态极值替换用以获取杂波估计中值并弥补记忆效应造成的影响，最后利用杂波估计中值来计算检测阈值。

2. 检测模型

2.1 杂波与目标

在判决雷达回波 $x(t)$ 中是否存在目标时，根据二元假设检测可知

$\left. \begin{aligned} & x(t) = n(t),\qquad\;\;\;{{\rm{H}}_0} \\ & x(t) = n(t) + s(t),{{\rm{H}}_1} \end{aligned} \right\}$

(1)

其中， $s(t)$ 为目标信号， $n(t)$ 为杂波和噪声信号， ${{\rm{H}}_{\rm{0}}}$ 表示只有杂波， ${{\rm{H}}_{\rm{1}}}$ 表示存在目标。

杂波分布模型按照其概率密度函数参数个数可以分为单参和多参分布模型^[17-20]。例如瑞利分布、指数分布等只由1个参数决定其幅度分布规律的为单参分布模型；韦布尔分布、对数正态分布、莱斯分布等由多个参数决定其幅度分布规律的为多参分布模型。本文以杂波数据经平方率检波后服从指数分布为背景，其杂波概率密度函数(Probability Density Function, PDF)为

$f(x) = \frac{1}{\lambda }\exp \left( - \frac{x}{\lambda }\right),0 \le x$

(2)

其中， $\lambda$ 为每个杂波单元总功率。在 ${{\rm{H}}_{\rm{0}}}$ 假设下， $\lambda$ 为杂波和噪声的功率 $\mu$ ；在 ${{\rm{H}}_{\rm{1}}}$ 假设下， $\lambda$ 为杂波、噪声和目标的总功率 $(1 + s)\mu$ , $s$ 表示目标与杂波噪声功率 $\mu$ 的比值。由式(2)可得检测概率(Probability of Detection, PD)与虚警概率(Probability of False Alarm, PFA)为

$\begin{split} {\rm{PD}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_1}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_1}\ \ \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{{(1 + s)\mu }}} \right) \end{split}$

(3)

$\begin{split} {\rm{PFA}} =\,& {\rm{P}} [x > {x_{\rm{t}}}|{{\rm{H}}_{\rm{0}}}] = \int_{{x_{\rm{t}}}}^\infty {f(x){\rm{d}} x} |{{\rm{H}}_{\rm{0}}} \\ =\,& \exp \left( { - \frac{{{x_{\rm{t}}}}}{\mu }} \right) \end{split}$

(4)

在已知参数 $\mu$ 和预设虚警概率 ${\rm{PFA}}$ 时，可由式(4)求得检测阈值T为

$T = - \mu \ln ({\rm{PFA}})$

(5)

2.2 有序矩阵

对于雷达回波的RD谱数据 ${{X}}$ ，每帧数据既有距离信息也有多普勒信息，其表达式为

${{X}} = \left[ {\begin{array}{*{20}{c}} {{x_{1,1}}}&{{x_{1,2}}}& ··· &{{x_{1,d}}}& ··· &{{x_{1,D}}}\\ {{x_{2,1}}}&{{x_{2,2}}}& ··· &{{x_{2,d}}}& ··· &{{x_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{r,1}}}&{{x_{r,2}}}& ··· &{{x_{r,d}}}& ··· &{{x_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{x_{R,1}}}&{{x_{R,2}}}& ··· &{{x_{R,d}}}& ··· &{{x_{R,D}}} \end{array}} \right]$

(6)

其中， ${x_{r,d}}$ 表示RD谱 ${{X}}$ 中待检测的第 $(r,d)$ 单元数据，其中 $r$ 为第 $r$ 个距离元、 $d$ 为第 $d$ 个多普勒元。由于本文所提算法利用待检测杂波的历史数据来估计中值 ${y_{{\rm{mid}}}}$ 求解检测阈值T，为了准确地估算估计中值 ${y_{{\rm{mid}}}}$ ，这里利用有序矩阵 ${{M}}$ 来存储相关历史数据，有序矩阵 ${{M}}$ 的表达式为

${{M}} = \left[ {\begin{array}{*{20}{c}} {{{{Y}}_{1,1}}}&{{{{Y}}_{1,2}}}& ··· &{{{{Y}}_{1,d}}}& ··· &{{{{Y}}_{1,D}}}\\ {{{{Y}}_{2,1}}}&{{{{Y}}_{2,2}}}& ··· &{{{{Y}}_{2,d}}}& ··· &{{{{Y}}_{2,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{r,1}}}&{{{{Y}}_{r,2}}}& ··· &{{{{Y}}_{r,d}}}& ··· &{{{{Y}}_{r,D}}}\\ \vdots & \vdots &{\ddots}& \vdots &{\ddots}& \vdots \\ {{{{Y}}_{R,1}}}&{{{{Y}}_{R,2}}}& ··· &{{{{Y}}_{R,d}}}& ··· &{{{{Y}}_{R,D}}} \end{array}} \right]$

(7)

其中， ${{{Y}}_{r,d}}$ 为有序矩阵 ${{M}}$ 中第 $(r,d)$ 单元的有序数组，其长度为 $2K + 1$ ，用于存储 ${x_{r,d}}$ 历史数据，表达式为

${{{Y}}_{r,d}} = \left[ {\begin{array}{*{20}{c}} {{y_1}}\;\;\;{{y_2}}\;\;\; ···\;\;\; {{y_k}}\;\;\; ···\;\;\; {{y_{2K + 1}}} \end{array} } \right]$

(8)

其中， ${y_k}$ 为有序数组 ${{{Y}}_{r,d}}$ 中的第 $k$ 大的数据， $K$ 为本算法的半参考窗值，K值越大，所需的存储空间越多。待检测单元 ${x_{r,d}}$ 的杂波估计中值 ${y_{{\rm{mid}}}}$ 为 ${{{Y}}_{r,d}}$ 中的 ${y_{K + 1}}$ 。

3. DOM-CFAR原理

DOM-CFAR算法实现流程如图1所示，分为初始化、迭代优化与检测3个步骤。

图 1 DOM-CFAR算法实现流程图

下载: 全尺寸图片幻灯片

3.1 初始化

初始化时，以RD谱数据 ${{X}}$ 中待检测单元 ${x_{r,d}}$ 为中心的长度为 $2K + 1$ 的距离多普勒域滑窗进行排序整理，用滑窗中值 $x{'_{r,d}}$ 对有序数组 ${{{Y}}_{r,d}}$ 中的所有数据赋初值，如式(9)所示

${y_k} = x{'_{r,d}},k \in [1,2K + 1]$

(9)

初始化后，有序数组 ${{{Y}}_{r,d}}$ 中的数据 ${y_k}$ 服从 $N({u_0} + \tau ,{\sigma ^2})$ 的高斯分布，其均值与方差为

$\qquad {u_0} + \tau = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{y_k}} }}{{2K + 1}}$

(10)

$\qquad {\sigma ^2} = \frac{{\displaystyle\sum\limits_{k = 1}^{k = 2K + 1} {{{[{u_0} + \tau - {y_k}]}^2}} }}{{2K + 1}}$

(11)

其中， ${u_0}$ 为杂波理论中值， $\tau$ 为 ${y_k}$ 与 ${u_0}$ 的差值， ${\sigma ^2}$ 表示 ${y_k}$ 的方差。

由式(9)可知，初始化后有序数组 ${{{Y}}_{r,d}}$ 中 ${y_k}$ 的数值都为 $x{'_{r,d}}$ ，方差值为0，并且当K值越大，采样滑窗所包含的杂波样本越多， $x{'_{r,d}}$ 越趋近 ${u_0}$ ， $\left| \tau \right|$ 数值越小。

由于本算法初始化参考OS-CFAR的采样策略，所以在迭代初期可获得与其接近的算法性能。

3.2 迭代优化

算法迭代通过将有序数组 ${{{Y}}_{r,d}}$ 中偏离理论中值 ${u_0}$ 的数据替换，保留接近 ${u_0}$ 的数据，从而获取更加准确的估计中值 ${y_{{\rm{mid}}}}$ 。

在数值替换中对极小值 ${y_1}$ 和极大值 ${y_{2K + 1}}$ 进行动态轮询替换，以雷达扫描次数t的奇偶性为标志，当雷达扫描次数t为奇数时，将待检测单元 ${x_{r,d}}$ 赋值给 ${{{Y}}_{r,d}}$ 中的 ${y_1}$ ；当t为偶数时，将 ${x_{r,d}}$ 赋值给 ${y_{2K + 1}}$ 。重新对 ${{{Y}}_{r,d}}$ 进行排序以保持有序性，以被替换极值点为起点，未被替换极值点为终点进行单向冒泡排序。K值越大，排序耗时越多。

由于初始化时 ${{{Y}}_{r,d}}$ 中的数据都相同，所以在算法迭代中需要经过2K+1轮迭代才能将 ${{{Y}}_{r,d}}$ 中的初始数据全部替换，之后 ${{{Y}}_{r,d}}$ 中的 ${y_{{\rm{mid}}}}$ 开始变化，初始数据被全部替换也标志着算法渡过2K+1轮迭代的初期阶段。

3.3 检测阈值计算

通过杂波估计中值 ${y_{{\rm{mid}}}}$ 逆推求解杂波模型参数，即令 $\displaystyle\int_0^{{y_{{\rm{mid}}}}} {f(x){\rm{d}} x} {\rm{ = }}\frac{1}{2}$ 可得杂波参数 $\mu = \dfrac{{{y_{{\rm{mid}}}}}}{{\ln 2}}$ 。

由于预设虚警概率已知，可通过式(4)求解标称因子 $\alpha$ 和检测阈值T为

$\qquad \alpha = - \frac{{\ln ({\rm{PFA}})}}{{\ln 2}}$

(12)

$\qquad T = \alpha {y_{{\rm{mid}}}}{\rm{ = }} - \frac{{{y_{{\rm{mid}}}}\ln ({\rm{PFA}})}}{{\ln 2}}$

(13)

使用本算法进行检测处理时，若杂波数据服从单参分布模型，可直接按原流程处理；若杂波数据服从多参分布模型，则需要先保留其中1个参数为未知状态，将其它参数由参数估计方法(最大似然估计、矩估计等)估算，之后按照原流程估算这个未知参数，计算检测阈值进行检测处理。

3.4 半参考窗值K的取值

半参考窗数值K需要综合考虑内存需求、运算速度、迭代速度、期望检测性能这4个因素进行设置。为便于分析K值对迭代速度和期望检测性能的影响，这里设计如下2个评价指标

${m_{t,K}} = \frac{1}{N}\sum\limits_{i = 1}^N {|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}|}$

(14)

$\sigma _{t,K}^2 = \frac{1}{N}\sum\limits_{i = 1}^N {{{(|{y_{t,K,i}} - {{\bar y}_{{\rm{mid}}}}| - {m_{t,K}})}^2}}$

(15)

其中， ${y_{t,K,i}}$ 为杂波估计中值， ${\bar y_{{\rm{mid}}}}$ 为杂波理论中值， ${m_{t,K}}$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的均值， $\sigma _{t,K}^2$ 为 ${y_{t,K,i}}$ 与 ${\bar y_{{\rm{mid}}}}$ 之差绝对值的方差， $t$ 为迭代次数， $K$ 为半参考窗值， $i$ 为蒙特卡罗实验次数， $N$ 为蒙特卡罗实验总次数。均值 ${m_{t,K}}$ 越小，表明 ${y_{t,K,i}}$ 越接近 ${\bar y_{{\rm{mid}}}}$ 、准确性越高；方差 $\sigma _{t,K}^2$ 越小，表明稳定性越高；准确性与稳定性越高计算出的检测阈值更接近理论值，可获得更优的期望检测性能。

利用杂波功率 $\lambda$ 为0 dBm的指数分布杂波数据进行2×10⁴组蒙特卡罗实验，每组实验进行2000次迭代、K取值范围为[1, 100]，对获得的 ${m_{t,K}}$ 和 $\sigma _{t,K}^2$ 取自然对数绘制等高线图谱如图2、图3所示。

图 2 均值等高线

下载: 全尺寸图片幻灯片

图 3 方差等高线

下载: 全尺寸图片幻灯片

由图2、图3可知，两图左上角出现三角形异常区域，这是由算法迭代初期杂波估计中值不变暂稳态导致。K取值在[10,20]范围内时，两图的评价指标处于等高线山谷区域，均值和方差较低(异常区域数据除外)，准确性和稳定性较高，所以应尽量将K取值在[10,20]范围内。

4. 仿真分析

将DOM-CFAR与最优检测器(opt), CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行4×10³组蒙特卡罗仿真实验，其中背景杂波功率 $\lambda$ 为20 dBm, DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，预设虚警概率PFA为10^–4。

4.1 均匀环境下的性能分析

对均匀环境下单目标分别进行10次、200次、1000次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，其检测结果如图4、图5、图6所示。

图 4 均匀杂波下迭代10次

下载: 全尺寸图片幻灯片

图 5 均匀杂波下迭代200次

下载: 全尺寸图片幻灯片

图 6 均匀杂波下迭代1000次

下载: 全尺寸图片幻灯片

由图4可知，DOM迭代初期检测性能接近OS，优于SO，较CA, GO略差。

由图5可知，经过200次迭代DOM性能已优于参照对比算法。

由图6可知，经过1000次迭代DOM性能趋近于最优检测器性能。

4.2 多目标环境下的性能分析

在目标参考前窗设置相同信杂比的第2个目标形成拥有多目标的背景杂波环境进行10次、200次迭代优化，其中信杂比取值范围为[0 dB, 30 dB]，检测结果如图7、图8所示。

图 7 多目标下迭代10次

下载: 全尺寸图片幻灯片

图 8 多目标下迭代200次

下载: 全尺寸图片幻灯片

由图7可知，DOM迭代初期拥有与OS相近的多目标检测优势。

由图8可知，DOM检测性能不受邻近目标的影响，且经过200次迭代后，其检测性能已经优于参照对比算法。

4.3 杂波边缘环境下的性能分析

设置前500个距离单元杂波功率为20 dBm，后500个距离单元杂波功率为30 dBm，形成杂波边缘环境，分别进行10次、200次迭代检测，其结果如图9、图10。

图 9 杂波边缘下迭代10次

下载: 全尺寸图片幻灯片

图 10 杂波边缘下迭代200次

下载: 全尺寸图片幻灯片

由图9可知，DOM迭代初期的虚警概率在杂波边缘的低功率区相较于CA, GO和CM更接近预设值，在杂波边缘的高功率区相较于CA, SO更接近预设值。

由图10可知，DOM经过200次迭代，虚警概率在杂波边缘的低功率区更接近预设值，在杂波边缘的高功率区相较于CA, SO和OS更接近预设值。

4.4 记忆效应分析

DOM与CM都是利用杂波历史数据计算检测阈值的时序类恒虚警算法，算法处理时每帧数据(无论是否有目标能量信息)都参与迭代并影响后续检测，从而导致算法的记忆效应。

由于DOM的参考数据为估计中值 ${x_{{\rm{DOM}}}}(t)$ , CM的参考数据为估计均值 ${x_{{\rm{CM}}}}(t)$ ，为了统一评价指标进行如下处理

$\quad {s_{{\rm{DOM}}}}(t) = \frac{{|{x_{{\rm{DOM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{DOM}}}}}}$

(16)

$\quad {s_{{\rm{CM}}}}(t) = \frac{{|{x_{{\rm{CM}}}}(t) - {P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}|}}{{{P_{{\rm{dB}}}} \cdot {{\bar x}_{{\rm{CM}}}}}}$

(17)

其中， ${P_{{\rm{dB}}}}$ 为杂波功率值， ${\bar x_{{\rm{DOM}}}}$ 为理论中值， ${\bar x_{{\rm{CM}}}}$ 为理论均值。评价指标 $s$ 的数值越小，参考值越接近理论值，估算的检测阈值更准确。

先在均匀杂波无目标的环境下进行2000次迭代使DOM和CM的性能稳定，从第2001次迭代开始，连续引入50次单目标回波，仿真目标能量参与迭代对性能的影响，其中目标能量信杂比分别为4 dB, 8 dB, 16 dB，结果如图11所示。

图 11 目标干扰对性能的影响

下载: 全尺寸图片幻灯片

由图11可知，当目标干扰出现并参与迭代，对DOM的影响有一定延时，当目标干扰消失时，DOM会以更快的速度重新优化，故其拥有更优的抗目标干扰能力。

4.5 算法复杂度与耗时

由于算法中每个待检测单元都需要一组长度为2K+1的有序数组用于存储其历史数据，因此空间复杂度(使用大O表示法)为O(MK)，其中M由待检测单元总个数决定。在迭代优化中动态极值替换的逻辑计算量很小几乎可以忽略，对运算速度起决定影响的是单向冒泡排序，因此其运算复杂度为O(K)。

以大O表示法为评价指标对各个算法的空间、运算复杂度以及在硬件为i7-7700HQ CPU条件下处理10帧实测RD谱数据的串行、并行计算耗时进行分析，结果如表1所示，其中N由滑窗大小决定、K由半参考窗值决定、M由待检测单元总个数决定。

表 1 算法复杂度及耗时比较

算法名称	空间复杂度	运算复杂度	串行耗时(s)	并行耗时(s)
CA	1	MN	51.337	2.776
SO	1	MN	73.588	3.260
GO	1	MN	75.176	3.289
OS	1	MNlgN	449.651	7.632
CM	M	M	0.087	0.073
DOM	MK	MK	88.490	5.560

下载: 导出CSV

| 显示表格

由表1及前文中各算法在均匀、多目标和杂波边缘环境下的性能对比可知，本算法与现有大部分算法相比，虽然存储空间需求较大、处理速度优势不突出，但在性能上得到了较大提高。同时根据已进行的雷达低空小目标探测实验可知，其算法复杂度对于目前的GPU硬件条件来说是可行的。

5. 实测数据验证

使用武汉大学电波传播实验室的外辐射源雷达系统^[21]在有风力发电机组环境下的实测数据进行验证分析，探测条件如图12所示。

图 12 雷达探测环境

下载: 全尺寸图片幻灯片

实测数据经多径杂波和多普勒扩展杂波抑制后的RD谱中，风机叶片转动的多普勒扩展杂波残余仍然较强，形成杂波边缘环境，如图13所示。

图 13 实测RD谱

下载: 全尺寸图片幻灯片

在实测RD谱上添加100个信杂比为15 dB的随机目标则形成拥有多目标的杂波边缘环境。将DOM-CFAR与CA-CFAR, SO-CFAR, GO-CFAR, OS-CFAR和CM-CFAR进行实测验证，其中DOM-CFAR的半参考窗值K为15, CA-CFAR, SO-CFAR, GO-CFAR和OS-CFAR的保护单元为4、滑窗长度为30, CM-CFAR的遗忘因子为1/128，迭代总次数为400，预设虚警概率为10^–6，结果如表2所示。

由表2数据可知，在迭代初期，DOM相较于参照对比算法拥有较优的检测性能，而虚警性能相对较差。但是DOM通过迭代优化后，检测概率会稳步提升，虚警概率会迅速下降，从而获得更优的检测性能与虚警性能。

表 2 实测的检测概率与虚警概率

CFAR	CA	SO	GO	OS	CM迭代次数				DOM迭代次数
CFAR	CA	SO	GO	OS	10	100	200	400	10	100	200	400
检测概率(%)	72.4	51.3	78.8	86.8	1.3	62.1	76.4	93.0	88.0	98.3	98.1	98.2
虚警概率(10^–5)	2.68	9.58	1.93	22.06	0.13	2.93	4.25	4.45	20.93	5.08	2.85	1.28

下载: 导出CSV

| 显示表格

6. 结束语

本文提出了一种基于动态有序矩阵的外辐射源恒虚警检测算法，利用历史杂波数据来动态替换有序矩阵中的极值，使得杂波估计中值向理论中值收敛，通过估计中值求解检测阈值。仿真与实测结果表明，该算法在迭代初期拥有较优的检测性能，通过多次迭代可以快速达到检测概率与虚警概率的最优性能，在均匀杂波、多目标和杂波边缘的复杂环境中检测性能稳定，并且相较于同属时序类算法的CM-CFAR，拥有更优的抗目标干扰能力。

参考文献(48)

陶白江, 孙可, 曾丁. 建立我国的直升机空中120救援系统势在必行[J].中华损伤与修复杂志(电子版), 2014, 9(3): 247-249.

Tao Bai-jiang, Sun Ke, and Zeng Ding. Establishment of a helicopter air 120 rescue system is imperative[J]. Chinese Journal of Injury Repair and Wound Healing, 2014, 9(3): 247-249.

张璐晶. 孟祥凯: 加快低空空域管理改革[J]. 中国经济周刊, 2014, 14(10): 48-49.

Zhang Lu-jing and Meng Xiang-kai. Speed up low-altitude airspace management reform[J]. China Economic Weekly, 2014, 14(10): 48-49.

宗合. 真高1000米以下空域实行分类管理有序放开[N]. 中国航空报, 2014-11-25001.

Zong He. 1000 meters below the airspace of category management and orderly release[N]. China Aviation News, 2014-11-25001.

于祥明. 低空空域改革倒计时系列政策或年底前揭盅[N]. 上海证券报, 2014-11-22003.

Yu Xiang-ming. Low-altitude airspace reform countdown[N]. Shanghai Securities News, 2014-11-22003.

曾天翔.飞机事故及其原因统计分析[J]. 航空标准化与质量, 1998, 26(6): 37-43.

Zeng Tian-xiang. Statistical analysis of aircraft accidents and their causes[J]. Aeronautic Standardization Quality, 1998, 26(6): 37-43.

邹元振, 孙文胜. 民用直升机飞行事故分析及对策[J]. 直升机技术, 2007, 6(4): 58-64.

Zou Yuan-zhen and Sun Wen-sheng. Civil helicopter accident analysis and countermeasures[J]. Helicopter Technique, 2007, 6(4): 58-64.

刘孙波, 施平安, 陈聚和. 舰载直升机起降安全性分析及其评估[J]. 舰船电子工程, 2014, 34(2): 106-109.

Liu Sun-bo, Shi Ping-an, and Chen Ju-he. Analysis and evaluation of ship-based helicopter take-off and landing safety[J]. Ship Electronic Engineering, 2014, 34(2): 106-109.

唐大全, 毕波, 王旭尚, 等. 自主着陆/着舰技术综述[J]. 中国惯性技术学报, 2010, 18(5): 550-555.

Tang Da-quan, Bi Bo, Wang Xu-shang, et al.. Summary on technology of automatic landing/carrier landing[J]. Journal of Chinese Inertial Technology, 2010, 18(5): 550-555.

Yang Song-pu and Wang Yang-zhu. Research on visual navigation technology of unmanned aerial vehicle landing[C]. 2013 IEEE International Conference on Information and Automation (ICIA), Yinchuan, August, 2013: 463-467.

Saripalli S, Montgomery J F, and Sukhatme G. Visually guided landing of an unmanned aerial vehicle[J]. IEEE Transactions on Robotics and Automation, 2003, 19(3): 371-380.

Scherer S, Chamberlain L, and Singh S. Autonomous landing at unprepared sites by a full-scale helicopter[J]. Robotics and Autonomous Systems, 2012, 60(12): 1545-1562.

Kaloshin G A. Scientific fundamentals of laser aids for low visibility approach landing of aircraft[C]. The 4th Pacific Rim Conference on Lasers and Electro-Optics, Chiba Japan, 2001: 15-19.

赵太飞, 吴鹏飞, 宋鹏. 无线紫外光直升机辅助起降通信技术研究[J]. 激光杂志, 2014, 35(10): 9-13.

Zhao Tai-fei, Wu Peng-fei, and Song Peng. Research on the wireless ultraviolet communication technology in helicopter launching and landing guidance[J]. Laser Journal, 2014, 35(10): 9-13.

Gou Liang, Zhang Geng-xin, Bian Dong-ming, et al.. Relay scheme based on distributed Luby transform codes for InterPlaNetary internet[J]. China Communications, 2013, 10(10): 1-11.

慕建君, 焦晓鹏, 曹训志. 数字喷泉码及其应用的研究进展与展望[J]. 电子学报, 2009, 37(7): 1571-1577.

Mu Jian-jun, Jiao Xiao-peng, and Cao Xun-zhi. A survey of digital fountain code an its application[J]. Acta Electronica Sinica, 2009, 37(7): 1571-1577.

Byers J W, Luby M, and itzenmacher M. A digital fountain approach to reliable distribution of bulk data[C]. Proceedings of SIGCOMM 98 Conference on Applications, Columbia Canada, 1998: 56-67.

Luby M. LT codes[C]. 43rd Annual IEEE Symposium on Foundations of Computer Science in Proc (FOCS), Vancouver Canada, 2002: 271-282.

Shokrollahi A. Raptor codes[J]. IEEE Transactions on Information, 2006, 52(6): 2551-2567.

Zhang Wu-ping, Jing Xiao-rong, Zhang Zu-fan, et al.. Image transmission with UEP-LT over MIMO channels[C]. Wireless and Optical Communication Conference (WOCC), Chongqing, 2013: 105-108.

Rajkumarsingh B and Basant S. LT codes with block duplication and ring type constellation for unequal error protection[C]. 2014 IEEE Region 10 Symposium, Kuala Lumpur Malaysia, 2014: 336-341.

Barbot N, Torkestani S S, Sahuguede S, et al.. LT codes performance over indoor mobile wireless optical channel[C]. IET International Symposium on Communication Systems, Networks and Digital Signal Processing (CSNDSP), Poznan Poland, 2012: 18-20.

Mladenov T, Nooshabadi S, and Kiseon Kim. Strategies for the design of raptor decoding in broadcast/multicast delivery systems[J]. IEEE Transactions on Consumer Electronics, 2010, 56(2): 423-428.

Apavatjrut A, Goursaud C, Jaffres-Runser, et al.. Toward increasing packet diversity for relaying LT fountain codes in wireless sensor networks[J]. IEEE Communications Letters, 2011, 15(1): 52-54.

Hussain I, Land I, Chan T H, et al.. A new design framework for LT codes over noisy channels[C]. 2014 IEEE International Symposium on Information Theory (ISIT), Honolulu USA, 2014: 2162-2166.

范喜全, 匡镜明. 一种复杂环境下的战术通信信道仿真方法[J]. 系统仿真学报, 2008, 20(9): 2502-2504.

Fan Xi-quan and Kuang Jing-ming. Simulation method of tactical communication channel under complex environment[J]. Journal of System Simulation, 2008, 20(9): 2502-2504.

梁大为, 鲍振武. 无线光通信性能影响因素分析[J]. 电子测量技术, 2006, 29(2): 21-22.

Liang Da-wei and Bao Zhen-wu. Analyses of influential factors of wireless optical communication performance[J]. Journal of Electronic Measurement and Instrument, 2006, 29(2): 21-22.

李华, 宋光辉, 吴鹏超. 国外森林防火通信技术现状与展望[J]. 中国林业, 2008, 10(19): 50-51.

Li Hua, Song Gang-hui, and Wu Peng-chao. Forest fire situation and prospects of foreign communications technology[J]. Forest Research, 2008, 10(19): 50-51.

柯熙政. 紫外光自组织网络理论[M]. 北京: 科学出版社, 2011: 3-43.

Ke Xi-zheng. UV Adhoc Communication Network Theory[M]. Beijing: Science Press, 2011: 3-43.

Jing Yue, Lin Zi-huai, Vucetic B, et al.. The design of degree distribution for distributed fountain codes in wireless sensor networks[C]. 2014 IEEE International Conference on Communications (ICC), Sydney NSW, 2014: 5796-5801.

Zhu Hong-peng, Zhang Geng-xin, and Li Guang-xia. A novel degree distribution algorithm of LT codes[C]. IEEE International Conference on Communication Technology, Hangzhou, 2008: 221-224.

Li Lu-ying, Li Zong-yan, and Wang Wen-bo. Adaptive iteration for fountain decoding[J]. The Journal of China Universities of Posts and Telecommunications, 2010, 17(2): 22-25.

Chen Gang, Feras Abou-Galala, and Xu Zheng-yuan. Experimental evaluation of LED-based solar blind NLOS communication links[J]. Optical Society of America, 2008, 16(19): 59-68.

Zuo Yong, Xiao Hou-fei, Wu Jian, et al.. Closed-form path loss model of non-line-of-sight ultraviolet single-scatter propagation [J]. Optics Letters, 2013, 38(12): 2116-2118.

El-Shimy M A and Hranilovic S. Spatial-diversity imaging receivers for non-line-of-sight solar-blind UV communications[J]. Journal of Lightwave Technology, 2015, 33(11): 2246-2255.

施引文献

期刊类型引用(5)

1.	刘言，刘宁波，黄勇，王中训. 利用相位特征筛选参考单元的改进CFAR方法. 烟台大学学报(自然科学与工程版). 2023(03): 371-378 . 百度学术
2.	杜改丽，封治华. 一种改进现有机载雷达探测性能的多普勒处理新技术. 现代雷达. 2022(02): 75-83 . 百度学术
3.	邹俊杰，程丰，万显荣. 外源雷达空时联合恒虚警检测分析与实验. 雷达科学与技术. 2022(04): 415-420+428 . 百度学术
4.	周子铂，王彬彬，张朝伟，刘建卫，徐颖鑫，王志会. 基于迭代对消的外辐射源雷达目标检测方法. 雷达科学与技术. 2022(05): 555-564 . 百度学术
5.	杨威，崔恒荣. 77GHz近程雷达CFAR算法研究. 通信技术. 2021(09): 2125-2131 . 百度学术