基于分数阶傅里叶变换的多辐射源被动直接定位算法

陈芳香 易 伟 周 涛 孔令讲

(电子科技大学信息与通信工程学院 成都 611731)

1 引言

多基站雷达无源定位系统中，雷达本身不发射电磁波，而是借助外辐射源信号或直接截获目标辐射源发射信号进行处理以实现相应的功能；多基站雷达无源定位具有功耗低、隐蔽性好等优势，在军事、民事中得到广泛应用[1–7]。常见的多基站雷达无源定位算法有间接定位算法和直接定位(Direct Position Determination, DPD)算法，也称为两步定位算法和一步定位算法。前者首先从截获的目标观测信号中提取相关参数(如到达角(Angle Of Arrival,AOA)、到达时差(Time Difference Of Arrival,TDOA)等)，然后利用提取到的参数通过数学运算获得目标位置[1–4]。后者则没有参数提取的过程，而是将截获的信号进行简单处理后直接传输至中央处理器并对目标进行定位[5–7]。两步定位算法中，参数提取的过程在每个接收雷达基站独立进行，忽略了所有观测回波信号对应同一个目标这一客观约束，在低信噪比(Signal-to-Noise Ratio, SNR)情况下其定位性能将严重恶化；而DPD算法则能充分利用各接收基站的观测回波信息，具有比两步定位算法更优的定位性能[5]。

针对多基站雷达无源直接定位这一背景，Weiss推导了在发射信号已知和发射信号未知两种情况下的最小二乘估计表达式[5]。该DPD算法只需通过对目标所在2维平面或3维空间区域进行网格搜索，寻找使相应的目标函数最大化的网格点作为辐射源目标的位置估计。对于发射信号已知的情况(例如训练信号或同步信号)，通过直接搜索待估参数使得目标函数最大化即可实现目标的定位；对于更为常见的发射信号未知的情况，根据观测信号构建与目标位置相关的矩阵，并通过一定的数学变换直接搜索该矩阵的最大特征值进而完成参数的估计。一般情况下，前者(DPD-known算法)可以获得最佳的定位性能，而后者(DPD-unknown算法)由于有部分发射信号特征信息的丢失，其定位性能要劣于前者。在很多情况下，辐射源发射信号的模型已知而具体参数未知，例如未知信号参数的线性调频(Linear Frequency Modulation, LFM)信号，可通过估计发射信号参数来进一步提高目标定位精度。

LFM信号具有低截获率、高分辨率、大时宽带宽积等优势，广泛应用于雷达、通信及地震识别等领域[8,9]。因此，对于LFM信号辐射源进行多基站雷达无源直接定位具有重要的理论意义和应用价值。文献[6]为解决单个LFM信号辐射源的定位问题，通过短时傅里叶变换(Short Time Fourier Transform, STFT)推导了联合估计LFM信号参数和目标位置的DPD算法，其定位精度明显优于DPD-unknown算法。但是上述DPD算法无法实现多目标定位，因此在实际场景中对多个LFM信号辐射源的定位问题亟待解决。

为确定监测区域内多个未知信号参数的LFM信号辐射源的位置，本文提出一种基于分数傅里叶变换(FRactional Fourier Transform, FRFT)和DPD算法相结合的多目标定位算法。FRFT作为一种广义的傅里叶变换(Fourier Transform, FT)，在波传播和光学领域是研究、分析和处理信号的一种重要手段[9–12]。因FRFT对LFM信号具有良好能量聚集性，被普遍用来估计LFM信号的相关参数[11,13]。本文首先对观测信号应用FRFT得到其功率谱，并结合目标提取算法进而实现辐射源数量估计；其次利用基本分类方法对估计得到的相关信号参数进行分类，并对感兴趣参数空间进行多次4维网格搜索进而实现多辐射源位置及相应信号参数估计；最后，通过仿真实验验证本文提出算法的定位性能。仿真结果表明该算法能够准确定位监测区域内的多个辐射源目标，且定位精度优于DPD-unknown算法。

2 信号模型与问题导出

3 基于FRFT的多目标直接定位算法

本文提出的基于FRFT的多目标直接定位算法(DPD-FRFT算法)思路：首先，利用FRFT处理接收信号得到其功率谱，并应用目标提取算法提取功率谱局部极值点，进而获得目标数量的估计值；然后，基于FRFT域中相应参数与LFM信号参数的关系式，可以得到相关信号参数的估计值；其次，利用基本分类算法可进一步将多目标问题转化为多个单目标问题；最后，基于4维网格搜索方法便可获得单个目标的位置及相应信号参数估计值。高维似然估计算法是理论最优的，而DPD-FRFT算法作为一种次优的定位算法，通过将高维的多目标定位问题转化为多个低维的单目标定位问题来降低计算量，因此可以认为该算法是在计算复杂度和定位精度之间折中的定位算法。

3.1 分数傅里叶变换简介

3.2 目标数量估计

3.3 信号参数分类

其中，S是Ts,f的量纲归一化因子；n0是LFM信号起始时间的序号，在接收信号中表现为由发射时间及传输时延决定的发射信号的起始样本点编号；f1表示在得到u0,p0后即可得到的初始频率的部分估计值，f2表示与时延相关的初始频率的部分估计值。由式(24)可知，初始频率和调频斜率的估计是基于量纲归一化方法的。

3.4 目标位置估计

4 仿真实验及结果分析

本节通过Monte Carlo仿真实验对比分析本文提出的DPD-FRFT算法、文献[5]中提出的DPD-unknown算法与DPD-known算法的性能。显然，前文推导的高维最大似然估计器在理论上是最优的，相比DPD-FRFT算法、DPD-unknown算法与DPD-known算法有最小的估计误差。然而由于计算的限制，实施该高维最大似然估计器花费的时间过长、代价过大，因此用DPD-known算法替代该算法。可以认为，DPD-known算法由于利用了发射信号完整的无误差的信息，因此是这几种定位算法误差的下界。

如图3所示，在一个2维笛卡尔坐标平面上，3个辐射源的位置坐标为 (-2，0)km , ( 2，0)km,(0，2)km。5个接收基站放置在以远点为圆心、半径为6 km 的圆弧上。SNR变化区间为-15～ 1 5 dB，间隔 3 dB。

辐射源均发射LFM信号，其信号参数如表1所示。

表1 信号参数设置Tab.1 Signal parameter setting

4.1 目标数量仿真

观测信号与式(5)具有相同形式，该信号由3个经过一定时延的LFM信号和噪声叠加得到。对该信号进行FRFT，得到的3维功率谱如图4所示。

由图4可知，当SNR=5 dB时，观测信号的功率谱中明显有3个局部极值点(红色椭圆标出)；但当SNR=–10 dB时，虽然也可分辨出功率谱的3个极值点，但由于受噪声影响并不明显。这表明应用FRFT可有效实现对LFM信号辐射源数量的估计，但受信噪比影响严重。其具体关系如图5所示。

由图5可知，在SNR大于0 dB左右时目标数量能够得到较好的估计。

4.2 定位性能仿真

本文用均方根误差(Root Mean Square Error,RMSE)衡量定位精度，其定义如下：

其中，||·||表示斐波那契范数；N是Monte Carlo实验次数，这里N=200;表示在第i次实验中得到的某个辐射源位置的估计值。明确一点，利用图像膨胀算法提取到的观测信号功率谱的极值点只对应潜在的目标，该过程是基于恒虚警率(Constant False Alarm Rate, CFAR)的，而目标数量的估计值是基于这些可能的目标，并且定位误差处于一定范围内。

仅考虑单目标情况下DPD-unknown算法、DPD-known算法和本文提出的ST-DPD-FRFT算法(针对单个目标的DPD-FRFT算法)的性能比较。需要注意的是，该种对比方式使得DPD-unknown算法有一定的优势；但若ST-DPD-FRFT算法的定位性能优于DPD-unknown算法，那么该算法同多目标DPD算法将具有更大的优势。

由图6可知，在高信噪比下(约为 1 0 dB)，4种定位算法都取得了较好的定位精度；随着信噪比降低至 5 dB，算法的定位性均出现不同程度的恶化；且由于受目标数量的影响，DPD-FRFT算法对单个目标的定位精度要劣于ST-DPD-FRFT算法。在所有算法中DPD-known算法的定位精度最优，而DPD-FRFT算法只有在高信噪比(0 dB以上)情况下其性能才接近DPD-known算法。这主要是由于在低信噪比情况下，DPD-FRFT算法对LFM信号参数的估计误差增大，从而影响其定位精度。为了验证这一点，定义信号参数估计的准确率(Correct Rate of Estimation, CRE)：

其中，Nc表示在N次Monte Carlo实验中信号参数的估计误差小于8%的次数。

图7表示目标发射信号的初始频率和调频斜率的估计准确率随信噪比的变化曲线。由图7可知，在信噪比低于0dB时信号参数的估计准确率开始降低。进一步对比图5和图7可知，目标个数的估计准确性也取决于信号参数的估计准确率。这一结果表明目标发射信号参数的估计准确性直接影响了辐射源目标的定位性能。

5 结论

针对多个未知LFM信号辐射源，本文提出的DPD-FRFT算法在控制计算量的同时，可有效联合估计出目标位置和信号参数。该算法有益探索了多目标、多基站无源直接定位问题，获得了较好的精度定位。