基于凸优化方法的认知雷达低峰均比波形设计

郝天铎 崔 琛 龚 阳 孙从易

①(国防科技大学电子对抗学院 合肥 230037)

②(96630部队 北京 102206)

1 引言

认知雷达是一种智能化雷达，可以根据对环境和目标信息的认知自适应地调整发射波形，是未来雷达的发展趋势之一[1,2]。认知雷达的发射波形不仅与其所面向的任务有关，而且也受某些约束条件的限制。现阶段，在进行面向检测的波形设计时，大部分方法仅考虑了发射能量的约束，而未对所设计信号的包络加以制约。在工程实际中，为了能够使雷达发射机发挥其最大效能通常要求雷达发射波形具有较低的峰均比(Peak-to-Average power Ratio,PAR)或者恒定包络[3,4]。然而，恒模约束又过于苛刻，往往与雷达的检测或估计性能不能兼顾，与恒模约束相比，PAR是一种更加泛化的约束，恒模波形只是PAR=1时的一种特殊波形，低PAR约束下的波形将具有更高的自由度[5–7]。

在信号相关杂波背景下进行检测波形设计时，通常选择最大化信干噪比(Signal-to-Interferenceplus-Noise Ratio, SINR)作为优化准则[8,9]。然而，基于该准则进行发射波形和接收机滤波器联合设计时，现有大部分方法主要将匹配滤波器(Matched Filter, MF)作为最佳接收机，而在相关杂波背景下，采用失配滤波器(MisMatched Filter, MMF)可以获得更优的性能。文献[10]对杂波背景下发射波形和MMF的联合设计问题进行了研究，该方法可以保证SINR逐渐递增并收敛于一个定值。但其只对发射能量进行了约束，所求波形的包络不受限制。为了增加恒模或者低PAR约束，文献[11–13]采用凸优化方法，通过半正定松弛(Semi-Definite Relaxation, SDR)技术，对恒模和低PAR波形的优化问题进行了研究。然而，上述文献中均未考虑信号相关杂波对波形设计的影响。文献[14]弥补了这方面的不足，在信号相关杂波背景下，将恒模波形的求解引入到认知雷达中，通过SDR高斯随机化的方法，得到了针对扩展目标检测的最优恒模波形，该波形通过较小的SINR损失即可换来发射机效能的提高，但仍然存在以下不足：(1)只进行了恒模约束下的波形设计，与恒模约束相比，PAR约束下的波形将具有更高的自由度和更好的实用性；(2)采用SDR得到最优波形矩阵解后，通过高斯随机化的方法从该矩阵解中提取最优波形向量解时，由于高斯随机化方法具有一定的随机性，需要取较多的随机高斯变量才能以较高概率得到最优可行解，因此算法运算量偏大。

综上，本文将波形的约束条件松弛为低PAR，同时考虑信号相关杂波的影响，并且兼顾算法的运算速率，提出了一种面向扩展目标检测的低PAR发射波形和接收机MMF的联合设计方法。该方法可同时兼顾波形的检测性能和放大器效能，在给定的PAR范围内可使波形的输出SINR达到最大。此外，本文采用秩1近似法结合最近邻方法[15]，取代高斯随机化方法，进一步降低了算法运算量(详见第4节分析)。仿真结果验证了所提方法理论分析的正确性和算法的有效性。

2 信号模型

本文考虑相关杂波背景下针对扩展目标检测的低PAR波形设计。目标的散射特性用目标冲激响应(Target Impulse Response, TIR)表示，相关杂波用杂波冲激响应(Clutter Impulse Response,CIR)表示。为方便起见，本文对单输入单输出(Single Input Single Output, SISO)雷达进行研究，主要在时域对离散的时间信号进行分析，信号模型如图1所示。

3 波形设计方法

3.1 确定TIR假设下的问题描述

3.2 低PAR波形设计

4 算法性能分析

4.1 算法收敛性分析

4.2 算法复杂度分析

本文主要与文献[14]算法的运算复杂度进行对比。文献[14]采用SDR高斯随机算法，由于高斯分布的随机性，产生的随机变量越多，得到最优解的概率也就越大，该文献平均需要K=10000个随机变量才能得到较好的可行解，其算法复杂度为：。本文的算法复杂度为：。可以看出，本文算法复杂度的第2项要远小于文献算法的第2项。然而，值得说明的是，虽然本文算法在第2步中提高了运算速率，但主要运算量集中在第1步，因此总体而言在算法复杂度上与文献[14]相比只是略有提高，而本文算法的优势在于波形包络的自由度得到了提升。

5 实验仿真

5.1 算法的有效性验证

令PAR=1，初始波形为随机的相位编码波形，比较算法中波形最优矩阵解和波形最优向量解的检测性能，如图3(a)所示。可以看出，随着迭代次数的增加，二者均递增收敛，且后者的检测性能要略差于前者，这是因为从矩阵解转化为向量解时，还需要再进行一次PAR的约束，缩小了可行解的求解区域。图3(b)验证了不同初始波形对算法的影响。其中，优化复高斯包络初始波形和优化随机相位编码波形是通过文献[10]中的算法1，将其进行一步优化而得到的初始波形。可以看出，优化初始波形的收敛速度和初始SINR都要更优；此外，线性调频(Linear Frequency Modulation, LFM)信号的性能较差。这表明算法对初始值是敏感的，为了提高收敛速率可选择性能较好的波形作为初始发射波形。

图4将本文算法与现有算法进行了对比。设门限τ=10-6，最大迭代次数κ=1000，分别与文献[14]中的恒模波形、文献[10]中能量约束下的最优波形、恒模波形(恒模波形3)、文献[4]中频域设计的最优波形以及LFM信号进行比较。可以看出，在恒模波形中，文献[14]和本文波形性能相当，恒模波形3的性能最差，这是因为其直接提取能量约束下最优波形的相位而成为恒模波形，并非最优的恒模波形；此外，本文算法性能要优于文献[4]，这是因为文献[4]频域波形不含有相位信息，变换到时域时会带来检测性能的下降；文献[10]算法是能量约束下的最优波形，本文产生的恒模波形的SINR与其相差约0.5 dB，这是因为能量约束下有着更高的包络自由度，此外，本文算法要明显好于LFM信号，而且本文方法可设计PAR在1～50之间的任意波形，这说明了本文方法对波形包络设计的灵活性。

图5在不同杂噪比(Clutter-to-Noise Ration, CNR)下对各种算法进行比较。图5(b)是图5(a)的放大。由于在PAR=1时，本文方法与文献[14]性能较为接近，故在该仿真中不与文献[14]作比较。可以看出，随着CNR的增大，本文方法产生的恒模波形检测性能优于文献[4]和恒模波形3，验证了本文算法的有效性。

5.2 算法的性能与PAR之间的关系

图6对不同PAR(1.0, 1.3, 1.6, 2.0, 2.6, 3.0)下优化波形的输出SINR进行了比较，图中每条实线对应着不同的PAR，随着箭头的指向PAR的值由1.0增大到3.0。可以看出，随着PAR的增大，式(18)的可行解区域也变大，因此输出SINR不仅收敛速度更快，而且最终收敛值也会相应变大。当PAR＞2时，得到的SINR最终收敛于同一值，接近能量约束下SINR的上限。图7为不同PAR波形实部和虚部的表示，代表能量约束下的波形，其中点的分布半径较大，说明波形幅度起伏较大，不能保证所设计波形的恒模特性，不利于实际应用。当PAR=1时，本文方法产生的点位于单位圆上，产生的是恒模波形；而当PAR=2时，点的分布半径与相比较小，然而，由前面分析可知，PAR=2时波形的检测性能已经与能量约束下波形较为接近，因此低PAR下产生的波形实用性较强。

6 结论

针对传统雷达波形设计方法中发射机发射功率不被充分利用的问题，以PAR为约束条件，从时域出发，在信号相关杂波背景下，提出了一种低PAR发射波形和接收机联合设计方法。相比现有方法，本文方法可以实现任意PAR下的联合设计。理论分析和仿真实验表明，通过将非凸问题转化为凸问题，可以有效提高相关杂波背景下距离扩展目标的检测性能；随着PAR的增大，波形性能逐渐接近能量约束下的曲线；相比于高斯随机法，通过将秩1近似法和最近邻算法相结合，可以在波形性能保持相当的情况下使算法复杂度进一步降低。值得注意的是，本文分析了确定TIR和随机CIR假设下的优化问题，对其他形式的TIR和CIR的组合，算法同样适用。