OFDM雷达多目标运动参数的近似最大似然估计

黄 瑞 杜小勇 胡卫东

(国防科技大学电子科学学院 长沙 410073)

1 引言

由于脉冲雷达发射的大功率信号在城市环境下容易干扰电台等其他通讯设备，可以借助通信信号来探测目标。近年来出现的多载波调制(Multi-Carrier Modulation, MCM)技术[1]引起了人们的关注，由此发展而来的正交频分复用(Orthogonal Frequency Division Multiplex, OFDM)[2]技术采用多路正交子载波进行信号调制，具有频率分集和波形分集的潜力，在通信、雷达一体化的发展背景下有着重要的研究价值。

相位编码OFDM信号具有多普勒高分辨力，用于雷达动目标检测时可以更精确地估计目标速度。不同的相位编码信号具有不同的自相关特性。霍夫曼编码可以降低信号自相关函数的整体旁瓣水平，但其包络峰均比(Peak-to-Mean Envelope Power Ratio, PMEPR)明显增大[3]。研究表明，barker码序列具有较好的综合性能[4]，本文将以13位barker码作为相位编码序列。

最近10年，国内外在OFDM信号特性与波形设计上取得了较好的研究成果[5–7]，因此OFDM雷达的信号处理问题受到了研究者的关注。张卫等人[8]利用Keystone变换在信号子载波域、快时间域和慢时间域进行联合解耦合处理，解决了目标的跨距离-多普勒单元走动问题，进而可估计出匀速运动下的多目标参数信息，但计算量较大；Lellouch等人[9]利用回波与发射信号的载频相位信息得到了点目标距离及径向速度估计，但要求目标在一个距离门内运动，即不发生越距离单元走动现象。除此以外，还需进一步研究OFDM雷达回波处理中面临的一些特殊问题，如速度补偿和多普勒解模糊。

在信号处理领域，最大似然估计是一种渐进有效估计量，但是对于多测量的非线性模型而言其计算量较大，不利于实际应用。为了提高计算效率，本文借鉴文献[10]中MIMO雷达信号处理的思路，结合OFDM信号多载波正交结构的特点，对信号进行通道分离，形成多通道信号。通过相关处理得到不同子载波上的距离像；利用Keystone变换进行速度补偿并解多普勒模糊，对同一载波的相同距离单元进行脉冲多普勒处理，得到每个子载波对应的多普勒频谱；进一步在子载波域作相参积累，得到距离-多普勒2维谱。通过谱峰搜索和CLEAN技术[11]的运用，从中提取出峰值位置对应的时延和多普勒参数。将其作为初值，结合观测数据的似然函数，利用牛顿迭代法获得更精确的参数估计，本文称之为近似最大似然估计。近似最大似然估计量可构成复合假设检验的重要环节，提升检测器的目标检测性能。论文组织结构如下：第2节给出了相位编码OFDM信号的回波模型；第3节给出了目标距离和速度参数的最大似然估计模型；为了提高目标运动参数估计的运算效率，第4节提出一种基于通道分离的近似最大似然估计算法；第5节利用仿真实验验证了算法的性能；第6节是总结。

2 相位编码OFDM信号模型

3 最大似然估计

4 基于通道分离的近似最大似然估计

直接利用式(11)求最大似然估计时计算量较大，在此考虑提高运算效率的参数估计方法。式(11)中相位项的子载波k分别与慢时间和快时间相耦合。针对不同的子载波，多脉冲联合处理时目标将在距离单元和多普勒单元上走动。鉴于此，文献[8]利用Keystone变换在信号子载波域、快时间域和慢时间域进行联合解耦合处理，解决了目标的跨距离多普勒单元走动问题，计算量约为 3K3Nlog2(KN)。事实上，式(11)中的相位项体现了不同子载波的回波信号具有不同的频率偏移，因此可以将多载波正交结构的OFDM雷达信号进行分离，通过多通道接收的方式增大距离分辨单元，避免目标的跨距离单元走动和3维Keystone变换。当时，对式(5)进行通道分离，即用参考信号与回波混频并经过低通滤波处理，得到各子载波通道上的信号

利用CLEAN技术从原信号中减去重构的子信号，然后对剩余信号继续上述的操作，直到CFAR检测不出峰值为止。算法处理流程如图1所示。

5 仿真实验

仿真参数设置如表1所示。根据表中数据并结合前文的分析可知，目标产生速度模糊的阈值是vuna=15 m/s ，出现跨距离单元走动的速度临界值是vt≈4.68m/s ，速度分辨率 Δvr=0.06m/s，且满足( 2vTr)/c≪tc。 由于|fd+fdk|≪Δf/2，多普勒频偏导致的子载频间串扰可忽略不计。在高斯白噪声背景下，2个目标距离为R1=10.0 km,R2=1.8 km,R3=17.0 km，径向速度v1=20 m/s,v2=24 m/s,v3=34 m/s，对应的目标散射强度分别为A01=5,A02=2,A03=1。虚警率Pfa=10-3，蒙特卡洛仿真100次。

信号经过通道分离和相关处理后，在子载波-多普勒平面的投影如图2所示。由图2可以看出，目标的多普勒频移随子载波的变化而变化，因此不能直接进行子载波域的相参积累。

经过Keystone变换和CLEAN处理后，相应子载波-多普勒平面的投影如图3所示。可见，多普勒频移与子载波之间的耦合得到校正。

对Keystone变换的数据在子载波域进行相参积累，获得信号的距离-多普勒2维谱如图4所示。其中图4(a)是聚焦于第1个目标的结果，图4(b)是剔除前两个目标后聚焦于第3个目标的结果。结果表明，本文所述的补偿方法能在快时间域、慢时间域以及子载波域将目标的回波能量积累起来，有利于目标检测和后续的参数估计。

表1 仿真参数Tab.1 Simulation parameters

下面考察近似最大似然估计量的性能。当v=20 m/s, S NRi在0～20 dB之间变化时，分别利用文献[8]所述的联合Keystone变换法、通道分离与Keystone变换组合法以及本文所提的近似最大似然估计法进行速度估计，估计量的均方根误差(RMSE)曲线如图5(a)所示。图5(b)给出了各种算法对应的计算量随脉冲数的变化情况。

由图5(a)可知，文献[8](点划线)与基于通道分离和Keystone相结合的方法(实线)的速度估计的RMSE曲线变化趋势类似。在相同的估计精度下，与上述两种方法相比，本文基于通道分离的近似最大似然估计法(点线)的输入信噪比S NRi改善约4 dB。同时，本文算法得到的RMSE非常接近Cramer-Rao下限。由图5(b)可知，在相同脉冲数下，本文算法的计算量较文献[8]的算法大幅降低。由于使用了牛顿迭代法，与基于通道分离和Keystone相结合的方法相比，计算量有所增加，但仍然远小于文献[8]中算法的计算量。因此，本文提出的基于通道分离的近似最大似然估计算法在估计精度和计算复杂度上具有综合优势。

6 总结

本文将OFDM通信信号应用在雷达动目标探测中，在通信、雷达一体化的发展背景下有着重要的应用前景。本文参考多载频MIMO雷达通道分离得到目标高分辨距离信息的方法，将OFDM信号的多载波正交结构与脉冲多普勒处理相结合，并借助Keystone变换解决了多普勒偏移问题。为了得到更好的估计精度，利用牛顿迭代法对似然函数进行优化，得到了基于通道分离的近似最大似然估计方法。仿真结果验证了算法的综合性能。今后还可针对机动目标相干化处理以及参数估计问题展开研究，扩展OFDM雷达的应用范围。