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

郝天铎 崔琛 龚阳 孙从易

郝天铎, 崔琛, 龚阳, 孙从易. 基于凸优化方法的认知雷达低峰均比波形设计[J]. 雷达学报, 2018, 7(4): 498-506. doi: 10.12000/JR18002
引用本文: 郝天铎, 崔琛, 龚阳, 孙从易. 基于凸优化方法的认知雷达低峰均比波形设计[J]. 雷达学报, 2018, 7(4): 498-506. doi: 10.12000/JR18002
Hao Tianduo, Cui Chen, Gong Yang, Sun Congyi. Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization[J]. Journal of Radars, 2018, 7(4): 498-506. doi: 10.12000/JR18002
Citation: Hao Tianduo, Cui Chen, Gong Yang, Sun Congyi. Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization[J]. Journal of Radars, 2018, 7(4): 498-506. doi: 10.12000/JR18002

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

doi: 10.12000/JR18002
基金项目: 国家部委基金
详细信息
    作者简介:

    郝天铎(1989–),男,内蒙古呼和浩特人,国防科技大学电子对抗学院博士生,主要研究方向为认知雷达信号处理和最优化理论。E-mail: haotianduo0423@126.com

    崔 琛(1962–),男,河北易县人,国防科技大学电子对抗学院教授,博士生导师,主要研究方向为雷达信号处理及雷达对抗技术。E-mail: kycuichen@163.com

    龚 阳(1992–),男,湖北随州人,国防科技大学电子对抗学院博士生,主要研究方向为雷达信号处理。E-mail: 13156527915@163.com

    孙从易(1992–),男,山东青岛人,国防科技大学电子对抗学院硕士生,现就职于96630部队,主要研究方向为MIMO雷达信号处理。E-mail: 1337325128@163.com

    通讯作者:

    郝天铎   haotianduo0423@126.com

Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization

Funds: The National Ministries Foundation
图(7)
计量
  • 文章访问数:  1047
  • HTML全文浏览量:  257
  • PDF下载量:  273
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-02
  • 修回日期:  2018-04-23
  • 刊出日期:  2018-08-28

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

doi: 10.12000/JR18002
    基金项目:  国家部委基金
    作者简介:

    郝天铎(1989–),男,内蒙古呼和浩特人,国防科技大学电子对抗学院博士生,主要研究方向为认知雷达信号处理和最优化理论。E-mail: haotianduo0423@126.com

    崔 琛(1962–),男,河北易县人,国防科技大学电子对抗学院教授,博士生导师,主要研究方向为雷达信号处理及雷达对抗技术。E-mail: kycuichen@163.com

    龚 阳(1992–),男,湖北随州人,国防科技大学电子对抗学院博士生,主要研究方向为雷达信号处理。E-mail: 13156527915@163.com

    孙从易(1992–),男,山东青岛人,国防科技大学电子对抗学院硕士生,现就职于96630部队,主要研究方向为MIMO雷达信号处理。E-mail: 1337325128@163.com

    通讯作者: 郝天铎   haotianduo0423@126.com

摘要: 为了提高雷达发射波形的检测性能,同时使发射机发挥其最大效能,以发射波形的低峰均比(PAR)为约束条件,该文提出了一种信号相关杂波背景下的认知雷达发射波形和接收机滤波器联合优化方法。首先,面向距离扩展目标检测问题,构建关于雷达输出信干噪比(SINR)的优化模型;然后将该模型转化为Rayleigh商形式,给出了接收机权值的解析表达式;在此基础上,通过半正定松弛,将关于发射波形半正定矩阵的非凸问题转化为凸问题,求得发射波形的最优矩阵解;最后,将秩1近似法和最近邻方法相结合,从最优矩阵解中提取出发射波形的最优向量解。该方法在给定PAR取值范围内可使波形的输出SINR达到最大,PAR=2时波形的SINR值与能量约束下优化波形的SINR值相同,并且比PAR=1时所得波形高出约0.5 dB。仿真结果验证了所提方法的有效性。

English Abstract

郝天铎, 崔琛, 龚阳, 孙从易. 基于凸优化方法的认知雷达低峰均比波形设计[J]. 雷达学报, 2018, 7(4): 498-506. doi: 10.12000/JR18002
引用本文: 郝天铎, 崔琛, 龚阳, 孙从易. 基于凸优化方法的认知雷达低峰均比波形设计[J]. 雷达学报, 2018, 7(4): 498-506. doi: 10.12000/JR18002
Hao Tianduo, Cui Chen, Gong Yang, Sun Congyi. Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization[J]. Journal of Radars, 2018, 7(4): 498-506. doi: 10.12000/JR18002
Citation: Hao Tianduo, Cui Chen, Gong Yang, Sun Congyi. Waveform Design for Cognitive Radar Under Low PAR Constraints by Convex Optimization[J]. Journal of Radars, 2018, 7(4): 498-506. doi: 10.12000/JR18002
    • 认知雷达是一种智能化雷达,可以根据对环境和目标信息的认知自适应地调整发射波形,是未来雷达的发展趋势之一[1,2]。认知雷达的发射波形不仅与其所面向的任务有关,而且也受某些约束条件的限制。现阶段,在进行面向检测的波形设计时,大部分方法仅考虑了发射能量的约束,而未对所设计信号的包络加以制约。在工程实际中,为了能够使雷达发射机发挥其最大效能通常要求雷达发射波形具有较低的峰均比(Peak-to-Average power Ratio, PAR)或者恒定包络[3,4]。然而,恒模约束又过于苛刻,往往与雷达的检测或估计性能不能兼顾,与恒模约束相比,PAR是一种更加泛化的约束,恒模波形只是PAR=1时的一种特殊波形,低PAR约束下的波形将具有更高的自由度[57]

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

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

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

      图  1  相关杂波下的信号模型

      Figure 1.  Signal model

      图1中, ${{{s}}} \in {\mathbb{C}^{{N_{{s}}} \times 1}}$ 代表发射波形, ${N_{s}}$ 为发射波形长度; ${{t}} \in {\mathbb{C}^{{N_{{t}}} \times 1}}$ ${{c}} \in {\mathbb{C}^{{N_{c}} \times 1}}$ 分别代表TIR和CIR,假定TIR和CIR长度相同[16],即 ${N_{t}} = {N_{c}}$ ,如果 $c\left( t \right)$ 的时长大于 $t\left( t \right)$ (即 ${N_{c}} > {N_{t}}$ ),本文拟采用的处理方式是通过补零延拓的方法对TIR进行补零处理,使TIR和CIR的采样点个数相同。 ${{n}} \in {\mathbb{C}^{{N_{n}} \times 1}}$ 代表噪声和干扰的总和, $ {{h}}$ 表示接收MMF冲激响应, $ {{x}}$ 代表来自目标和环境的回波, $ {{y}}$ 为滤波器输出后的信号,其长度 ${N_{{y}}}{\rm{ = }}{N_{{h}}}{\rm{ = }}{N_{{s}}} + {N_{t}} - 1$ 。该模型可表示为:

      $$ \begin{align} {{x}} =& {{t}} * {{s}} + {{c}} * {{s}} + {{n}} = {{Ts}} + {{Cs}} + {{n}} \\ = &{{St}} + {{Sc}} + {{n}} = {{{s}}_{{t}}} + {{{s}}_{{c}}} + {{n}} \\ \end{align} $$ (1)
      $$ {{y}} = {{h}} * {{x}} = {{Hx}} = {{Xh}} $$ (2)

      其中,“*”代表卷积运算。 ${{s}_{t}}$ ${{s}_{c}}$ 分别表示回波中的目标信息部分和信号相关杂波部分, $ {{T}}$ , $ {{S}}$ , $ {{C}}$ , $ {{H}}$ $ {{X}}$ 表示对应向量的卷积矩阵,这些矩阵均为Toeplitz阵。以TIR为例,令 ${{t}} = $ $ {\left[ {{{t}}\left( 1 \right)\;\;{{t}}\left( {\rm{2}} \right)\;\; ·\!·\!· \;\;{{t}}\left( {{N_{{t}}}} \right)} \right]^{\rm T}}$ ,则 ${{T}} = \left( {{{{t}}_1}\;\;{{{t}}_{\rm{2}}}\;\; ·\!·\!· \;\;{{{t}}_{{N_{{s}}}}}} \right)$ ,其中, ${{{t}}_i} = {\left( {{{0}}_{i - 1}^{\rm{T}},{{{t}}^{\rm{T}}},{{0}}_{{N_{s}} - i}^{\rm{T}}} \right)^{\rm{T}}}$ , $i = 1,2, ·\!·\!·,{N_{{s}}}$ , ${{{0}}_i}$ 代表i 维零向量, ${\left( \cdot \right)^{\rm T} }$ 代表转置,可得:

      $$ {{T}} = \left[\!\!\! {\begin{array}{*{20}{c}} {{{t}}\left( 1 \right)}&0& ·\!·\!· &0 \\ \vdots &{{{t}}\left( 1 \right)}& \ddots & \vdots \\ {{{t}}({N_{{t}}})}& \vdots & \ddots &0 \\ 0&{{{t}}({N_{{t}}})}& \ddots &{{{t}}\left( 1 \right)} \\ \vdots & \ddots & \ddots & \vdots \\ 0& \cdots &0&{{{t}}({N_{{t}}})} \end{array}}\!\!\!\! \right] \!\!\in {\mathbb{C}^{\left( {{N_{{s}}} + {N_{t}} - 1} \right) \times {N_{{s}}}}} $$

      (3)

      根据式(1)和式(2),并令噪声和杂波的和为:

      $$ {{w}} = {{{s}}_{{c}}} + {{n}} $$ (4)
      $$ {{y}} = {{Hx}} = {{H}}{{{s}}_{{t}}} + {{H}}({{{s}}_{{c}}} + {{n}}) = {{s}}_{{{{t,H}}}} + {{{w}}_{{H}}} $$ (5)

      因此,接收MMF输出信号的SINR可表示为:

      $$ \begin{align} {\rm SINR} & = \frac{{ \mathbb {E}\! \left( {{{s}}_{{{t}},{{H}}}^{\rm{H}}{{{s}}_{{{t}},{{H}}}}} \right)}}{{\mathbb {E}\! \left({{w}}_{{H}}^{\rm{H}}{{{w}}_{{H}}}\right)}} \\ & = \frac{{\mathbb {E}\! \left( {{{{s}}}_{t}^{\rm{H}}{{{{H}}}^{\rm{H}}}{{{H}}}{{{{s}}}_{t}}} \right)}}{\mathbb {E}\! \left( {{{s}}_{{c}}^{\rm{H}}{{{H}}^{\rm{H}}}{{H}}{{{s}}_{{c}}}} \right) + \mathbb {E}\! \left( {{{{n}}^{\rm H}}{{{H}}^{\rm H}}{{{Hn}}} }\right)} \end{align} $$ (6)

      其中, $\mathbb{E}$ 代表求均值。本文主要对确定TIR假设下的优化问题进行分析。

    • 本小节主要采用循环迭代的思想对发射波形 $ {{s}}$ 和接收机滤波器 $ {{h}}$ 进行优化。由于本文主要侧重于最优波形和接收机滤波器的设计,所以在后面的行文中假定TIR和CIR的先验知识可通过文献[17, 18]的方法获得。设 ${{a}} \in {\mathbb{C}^{{N_{a}} \times 1}}$ ,则有 $ {{{a}}^{\rm{H}}}{{a}} = $ $ {\rm{tr}}\left({{a}}{{{a}}^{\rm{H}}}\right)$ , $ {\rm{tr}}$ (·)代表矩阵的迹。因此,根据式(6),在确定TIR假设下目标函数可表示为:

      $$ \begin{align} &\mathop {\max}\limits_{{{s}},{{h}}} \frac{{{{t}}^{\rm{H}}}{{{S}}^{\rm{H}}}{{{H}}^{\rm{H}}}{{{HSt}}}}{\mathbb{E}\!\Bigr( {\rm{tr}}\!\left( {{{{HSc}}{{{c}}^{\rm{H}}}{{{S}}^{\rm{H}}}{{{H}}^{\rm{H}}}} }\right) \Bigr) \!+\! \mathbb{E}\!\Bigr( {{\rm{tr}}\!\left( {{{{Hn}}{{{n}}^{\rm{H}}}{{{H}}^{\rm{H}}}} }\right) }\Bigr)} \\ &\quad = \mathop {\max}\limits_{{{s}},{{h}}} \frac{{{{t}}^{\rm{H}}}{{{S}}^{\rm{H}}}{{{H}}^{\rm{H}}}{{{HSt}}}}{{\rm{tr}}\!\left( {{{{HS}}{{{R}}_{{c}}}{{{S}}^{\rm{H}}}{{{H}}^{\rm{H}}}} }\right) + {\rm{tr}}\!\left({{H}}{{{R}}_{{n}}}{{{H}}^{\rm{H}}}\right)} \end{align} $$ (7)

      其中, ${{{R}}_{n}}$ ${{{R}}_{c}}$ 分别为噪声和CIR的协方差矩阵,可利用其特征值和特征向量进行表示:

      $$ {{{R}}_{{c}}} = \sum\limits_{i = 1}^{{N_{{c}}}} {{\lambda _i}\left( {{{{R}}_{{c}}}} \right){{{u}}_i}\left( {{{{R}}_{{c}}}} \right){{u}}_i^{\rm{H}}\left( {{{{R}}_{{c}}}} \right)} $$ (8)

      其中, ${{{u}}_i}\left( {{{{R}}_{{c}}}} \right)$ 代表矩阵 ${{{R}}_{c}}$ 的第i个特征值 ${\lambda _i}\left( {{{{R}}_{{c}}}} \right)$ 对应的特征向量。

      假定 $ {{h}}$ 已知,从式(7)中提取出 $ {{s}}$ ,根据卷积运算的互易性原则,目标函数可变为:

      $$ \mathop {\max}\limits_{{s}} \frac{{{{{s}}^{\rm{H}}}{{{Ω}} _{{t}}}\left( {{h}} \right)\!{{s}}}}{{{{{s}}^{\rm{H}}}{{{Ω}} _{{c}}}\left( {{h}} \right)\!{{s}}}} $$ (9)

      其中,

      $$ {{{Ω}} _{{t}}}\!\left( {{h}} \right) = {{T}}{{{H}}^{\rm{H}}}{{HT}} $$ (10)
      $$ \begin{align} {{{Ω}} _{{c}}}\left( {{h}} \right) =& \sum\limits_{i = 1}^{{N_{{c}}}} {\lambda _i}\left( {{{{R}}_{{c}}}} \right){{\rm F}^{\rm{H}}}\!\Bigl[ {{{{u}}_i}\left( {{{{R}}_{{c}}}} \right)} \Bigr]\\ & \cdot {{{H}}^{\rm{H}}}{{H}}{\rm F}\!\Bigr[ {{{{u}}_i}\left( {{{{R}}_{{c}}}} \right)} \Bigr]\\ & + \left[ {{\rm tr}\left( {{{H}}{{{R}}_{{n}}}{{{H}}^{\rm{H}}}} \right)\Bigr/{P_{{s}}}} \right]{{{I}}\!_{{N_{{s}}}}} \end{align} $$ (11)

      其中,F代表的是由向量到Toeplitz矩阵的映射, ${P_{s}}$ 代表发射波形的最大能量值, ${{{I}}_{{N_{s}}}}$ 代表 ${N_{{s}}} \times {N_{{s}}}$ 的单位矩阵。

      同理,假定 $ {{s}}$ 已知,提取 $ {{h}}$ ,式(7)可变形为:

      $$ \mathop {\max}\limits_{{h}} \frac{{{{{h}}^{\rm{H}}}{{{{Θ}}}_{{t}}}\left( {{s}} \right){{h}}}}{{{{{h}}^{\rm{H}}}{{{{Θ}}}_{{c}}}\left( {{s}} \right){{h}}}} $$ (12)

      其中,

      $$ {{{{Θ}}}_{{t}}}\left( {{s}} \right) = {{\rm{F}}^{\rm{H}}}\!\left( {{{St}}} \right){\rm{F}}\!\left( {{{St}}} \right) $$ (13)
      $$ \begin{align} {{{{Θ}}}_{{c}}}\left( {{s}} \right) =& \sum\limits_{i = 1}^{{N_{{c}}}} {{\lambda _i}\left( {{{{R}}_{{c}}}} \right){{\rm{F}}^{\rm{H}}}\!\Bigr[ {{{S}}{{{u}}_i}\!\left( {{{{R}}_{{c}}}} \right)} \Bigr]{\rm{F}}\Bigr[ {{{S}}{{{u}}_i}\!\left( {{{{R}}_{{c}}}} \right)} \Bigr]} \\ &+ \sum\limits_{i = 1}^{{N_{{x}}}} {\lambda _i}\left( {{{{R}}_{{n}}}} \right) {{\rm{F}}^{\rm{H}}}\!\Bigr[ {{{{u}}_i}\!\left( {{{{R}}_{{n}}}} \right)} \Bigr]{\rm{F}}\Bigr[ {{{{u}}_i}\!\left( {{{{R}}_{{n}}}} \right)} \Bigr] \end{align} $$ (14)

      $ {{s}}$ 仅受能量约束且 $ {{s}}$ 已知时,式(12)可看做经典的广义Rayleigh商问题[19],可得 $ {{h}}$ 的优化解为:

      $$ {{h}} = \sqrt {{{ P}_{{h}}}} {{{u}}_{{\max}}}\Bigr[ {{{{Θ}}}_{{c}}^{ - 1}\left( {{s}} \right){{{{Θ}}}_{{t}}}\left( {{s}} \right)} \Bigr] $$ (15)

      其中, ${P_{h}}$ 代表对接收机滤波器冲激响应的限制, ${{{{u}}}_{\max}}\left[ {{{Θ}} _{{c}}^{ - 1}\left( {{s}} \right){{{Θ}} _{{t}}}\left( {{s}} \right)} \right]$ 代表矩阵 ${{Θ}} _{{c}}^{ - 1}\left( {{s}} \right){{{Θ}} _{{t}}}\left( {{s}} \right)$ 最大特征值对应的特征向量。同理,当 $ {{h}}$ 已知时,可得式(9)中发射波形的最优解为:

      $$ {{s}} = \sqrt {{P_{{s}}}} {{{u}}_{\max}}\Bigr[ {{Ω} _{{c}}^{ - 1}({{h}}){{Ω} _{{t}}}({{h}})} \Bigr] $$ (16)

      然而,考虑到实际的工程需求,当约束条件变为恒模或者低PAR时,问题就变得复杂起来。

    • 实际雷达系统中,发射波形的总能量是有限的,不失一般性,设发射波形的最大能量不超过 ${N_{{s}}}$ ,则PAR可定义为:

      $$ {\rm{PAR}}({{s}}) \triangleq \frac{{\mathop {\max}\limits_k {{\left| {{{{s}}_k}} \right|}^2}}}{{\displaystyle\frac{1}{{{N_{{s}}}}}{{\sum\limits_{k = 1}^{{N_{{s}}}} {\left| {{{{s}}_k}} \right|^2} }}}} = \mathop {\max}\limits_k {\left| {{{{s}}_k}} \right|^2} $$ (17)

      其中, ${{{s}}_k}$ $ {{s}}$ 的第k个采样。当 $ {{h}}$ 已知时,以式(9)为目标函数的优化问题可表示为:

      $$\left. \begin{align} &\mathop {{\max}}\limits_{{s}} \; \frac{{{{{s}}^{\rm{H}}}{{Ω} _{{t}}}\left( {{h}} \right){{s}}}}{{{{{s}}^{\rm{H}}}{{{{Ω }}}_{{c}}}\left( {{h}} \right){{s}}}} \\ & {\rm{s}.t.}\;\;\;\;{\rm{PAR}}({{s}}) \le \mu ,\;\;\;\mu \in \left[ {1,{N_{{s}}}} \right] \\ & \;\;\;\;\;\;\;\;{\left\| {{s}} \right\|^2} \le {N_{{s}}} \end{align} \!\!\!\!\!\!\!\!\!\!\!\!\!\! \right\} $$ (18)

      式中,当 $\mu = {N_{s}}$ 时,等价于能量约束,PAR就成了冗余约束,式(16)即为最优解;当 $\mu = 1$ 时,等价于恒模约束。由于问题式(18)的目标函数为非凸问题,并且又有PAR的约束,使得该问题成为NP-hard问题。由于式中的 ${{{Ω}} _{{t}}}\left( {{h}} \right)$ ${{{Ω}} _{c}}\left( {h} \right)$ 均为对称矩阵,令 ${{W}} = {{s}}{{{s}}^{\rm{H}}}$ ,此时 $ {{C}}$ 为半正定矩阵,因此可将该问题转化为半正定规划(Semi-Definite Programming, SDP)问题进行求解。令 ${{{{Ω }}}_1} = {{{{Ω}}}_{{t}}}\left( {{h}} \right)$ , ${{{{Ω }}}_2} = {{{{Ω }}}_{{c}}}\left( {{h}} \right)$ ,又由于 ${{{s}}^{\rm{H}}}{Ω} {{s}} = {\rm{tr}}\left({Ω } {{s}}{{{s}}^{\rm{H}}}\right)$ ,则式(18)可转化为:

      $$\left. \!\!\!\!\!\!\! \begin{align} & \mathop {\max}\limits_{{W}} \; \frac{{{\rm {tr}}({{Ω} _1}{{W}})}}{{{\rm {tr}}({{Ω} _2}{{W}})}} \\ & {\rm s.t}.\;\; {\rm {diag}}({{W}}) \le \mu ,\; \\ &\quad\ \ \, {\rm {tr}}({{W}}) \le {N_{{s}}},\;\;\;{{W}} \succeq 0,\;\;{\rm {rank}}({W}) = 1 \\ \end{align} \!\!\!\!\!\!\!\!\!\!\right\} $$ (19)

      其中,diag(·)代表矩阵的对角线元素,rank(·)代表矩阵的秩, ${{W}} \succeq 0$ 表示 $ {{W}}$ 为半正定矩阵。然而,式(19)仍然不是一个凸问题,可以通过SDR[20]技术,先把 ${\rm {rank}}({{W}}) = 1$ 的约束去掉,此时问题变为:

      $$\left. \begin{align} & \mathop {\max}\limits_{{W}} \;\frac{{{\rm {tr}}({{{Ω}} _1}{{W}})}}{{{\rm {tr}}({{{Ω}} _2}{{W}})}} \\ & {\rm s.t.}\;\; {\rm {diag}}({{W}}) \le \mu ,\; \\ & \quad\quad{\rm {tr}}({{W}}) \le {N_{{s}}},\;\;\;{{W}} \succeq 0 \end{align} \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\! \right\} $$ (20)

      接下来,需要对目标函数进行等量代换,使其转化为凸函数。由于 ${\rm {tr}}({{{Ω}} _1}{{W}})$ 的值为一常数,可令 ${\rm {tr}}({{{Ω}} _1}{{W}}){\rm{ = }}q$ , ${{W}'} = {{W}}/q$ ,则有

      $$ {\rm {tr}}({{{Ω}} _1}{{W}'})\,=\, {\rm {tr}}({{{Ω}} _1}{{W}}/q) \,=\, {\rm {tr}}({{{Ω}} _1}{{W}})/q \,=\, 1 $$ (21)

      由于仿射空间并不对除法运算封闭,所以需要对变量q进行处理,令 $p = 1/q$ 。此时,以 ${{W}}'$ p为变量,式(20)可以转化为:

      $$\left. \begin{align} & \mathop {\rm{min}}\limits_{{{W}'},p} \; \rm{tr}\left( {{{{Ω }}}_2}{{{W}'}} \right) \\ & {\rm{s}\rm{.t}\rm{.}}\;\; {\rm{diag}}\left( {{{W}}'} \right) \le \mu p,\;\;\; \\ &\quad\quad {\rm{tr}}\left( {{{W}}'} \right) \le {N_{{s}}}p,\;\;\;{{W}'} \succeq 0 \end{align}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\! \right\} $$ (22)

      式中,由于 ${{W}}'$ 也是半正定矩阵,式(22)也就变成了一个凸问题,该凸问题有成熟的求解方法,实际使用中有相应的工具箱。设 ${{{W}}'\!\!_{\rm opt}}$ ${p_{\rm opt}}$ 为式(22)的最优解,则 ${q_{\rm opt}} = 1/{p_{\rm opt}}$ ,可得式(20)的全局最优解为:

      $$ {{{W}}\!\!_{\rm {opt}}}{\rm{ = }}{{{{W}'}\!\!_{\rm {opt}}}{q_{\rm{opt}}}{\rm{ = }}{{{{{{{W}'}}\!\!_{\rm{opt}}}} / {{p_{\rm{opt}}}}} $$ (23)

      如果 ${{{W}}\!\!_{\rm {opt}}}$ 的秩为1,则可取 ${{{W}}\!\!_{\rm {opt}}}$ 的第1列作为可行解 ${{s}_{\rm {opt}}}$ 。然而,一般情况下,经过SDR后 ${{{W}}\!\!_{\rm {opt}}}$ 的秩并不为1[21]。若 ${\rm {rank}}\left( {{{{W}}\!\!_{\rm opt}}} \right) > 1$ ,根据秩1近似法[20],可以令

      $$ {W}\!\!_{\rm {opt}}\,\!\!\!\!\!^ * \ {\rm{ = }}{\lambda _{\max}}\left( {{{{W}}\!\!_{\rm {opt}}}} \right){{{u}}_{\max}}\left( {{{{W}}\!\!_{\rm {opt}}}} \right){{u}}_{{\max}}^{\rm{H}}\left( {{{{W}}\!\!_{\rm {opt}}}} \right) $$ (24)

      式中, ${{W}}\!\!_{{\rm{opt}\,}}\,\!\!\!\!\!^{\rm{*}}$ 代表 ${{{W}}\!\!_{\rm {opt}}}$ 的最大特征值与特征向量的乘积,在 ${{{W}}\!\!_{\rm {opt}}}$ 的所有特征值与对应特征向量乘积的组合中, ${{W}}\!\!_{\rm {opt}\,}\,\!\!\!\!\!^ * $ ${{{W}}\!\!_{\rm {opt}}}$ 最为近似,此时,就可以选取 ${{{{\tilde {s}}}_{\rm {opt}}}{\rm{ = }}\sqrt {{p_{_{s}}}} {{{u}}_{\max}}({{{W}}\!\!_{\rm {opt}}})$ 作为式(18)的一个候选解。然而,由于 ${{{{\tilde {s}}}_{\rm {opt}}}$ 只是近似解,其元素不一定满足PAR约束,所以需进一步进行优化,将该候选解映射为式(18)的最优可行解,优化问题可表示为:

      $$ \left. \begin{align} &\mathop {\min}\limits_{{s}} \;{\left\| {{{s}} - \sqrt {{p_{_{s}}}} \frac{{{{{u}}_{{\max}}}({{{W}}\!\!_{\rm{opt}}})}}{{{{\left\| {{{{u}}_{{\max}}}({{{W}}\!\!_{\rm{opt}}})} \right\|}_2}}}} \right\|_2} \\ & {\rm{s}}{\rm{.t}\rm{.}}\;\;{\rm PAR}({{s}}) \le \mu ,\;\;\;{\left\| {{s}} \right\|^2} \le {N_{{s}}} \end{align} \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\! \right\} $$ (25)

      此时,式(25)是一个凸问题,该凸问题有成熟的求解方法,实际使用中有相应的CVX工具箱[22]。但直接使用这些成熟的算法时,其复杂度较高。由于式(25)的目标函数和约束条件的形式与文献[15]中的最近邻算法相同,采用该算法可以对问题进行快速求解,每一步迭代的运算复杂度为 ${O}\left( {{N_{s}}\!^2} \right)$ ( ${N_{s}}$ 为序列的码长),远低于凸优化算法 ${ O}\left( {{N_{s}}\!^{4.5}} \right)$ 的运算复杂度。

      综上,基于凸优化的认知雷达低PAR波形设计算法可总结如下:

      步骤1 令 $i = 0$ ,随机产生的复波形 ${{s}^{(i)}}$ 作为初始发射波形,设置门限值 $\tau $ 和最大迭代次数 $\kappa $

      步骤2 令 $i = i + {\rm{1}}$ ,计算 ${{{Θ}} _{c}}\left( {s} \right)$ ${{{Θ}} _{t}}\left( {s} \right)$ ,通过式(15)求得最优接收MMF向量 ${{{h}}^{\left( i \right)}}$

      步骤3 通过求解式(22)得到 ${{W}}\,'\!\!\!_{\rm {opt}} \!\!^{\raisebox{2pt}{\left( i \right)}}$ ,进而得到式(20)中经半正定松弛后的最优矩阵解 ${{W}}\!\!_{\rm {opt}}\,\!\!\!\!^{\raisebox{2pt}{\left( i \right)}}$

      步骤4 如果 ${{W}}\!_{\rm {opt}}\,\!\!\!\!\!^{\raisebox{2pt}{\left( i \right)}}$ 的秩为1,直接提取 ${{W}}\!_{\rm {opt}}\,\!\!\!\!\!^{\raisebox{2pt}{\left( i \right)}}$ 的第1列作为发射波形的可行解 ${{s}}^{\left( i \right)}$ ;否则,采用秩1近似法从 ${{W}}\!_{\rm {opt}}\,\!\!\!\!\!^{\raisebox{2pt}{\left( i \right)}}$ 提取近似解 ${{\tilde {s}}}_{\rm {opt}}^{\left( i \right)}$ ,将其代入式(25)利用最近邻算法得到可行解 ${{s}}^{\left( i \right)}$ ,并计算 ${\rm {SINR}}^{\left( i \right)}$

      步骤5 如果 $\left| {{\rm {SINR}}^{\left( i \right)} - {\rm {SINR}}^{\left( {i - 1} \right)}} \right| < \tau $ 或者迭代的次数大于 $\kappa $ ,迭代停止,否则返回步骤2。

    • 本文算法通过对 $ {{s}}$ $ {{h}}$ 交替迭代进行寻优,对于不同的初始 $ {{s}}$ 或者 $ {{h}}$ ,沿着 $ {{s}}$ $ {{h}}$ 方向的优化结果也是不同的。当初始波形 ${{{s}}_0}$ 已知时,算法会沿着 ${{{s}}_0}$ 的方向求得 ${{h}_1}$ 的最优值,此时的 ${{{h}}_1}$ 仅仅是 ${{{s}}_0}$ 方向上的最优值,而非全局最优值;接下来又会沿着 ${{{h}}_1}$ 的方向得到 ${{{s}}_1}$ 的最优值,如此交替迭代直到满足收敛条件为止。下面给出简要的证明,令 ${\rm {SINR}}^{\left( k \right)}$ 表示第k次迭代后的输出SINR值,则可得

      $$ \begin{align} {\rm {SINR}}^{\left( k \right)} =& {\rm {SINR}}\left( {{{s}_k},{{h}_k}} \right) \le {\rm {SINR}}\left( {{{s}_k},{{h}_{k{\rm{ + }}1}}} \right) \\ \le &{\rm {SINR}}\left( {{{s}_{{k}{\rm{ + }}1}},{{h}_{k{\rm{ + }}1}}} \right) = {\rm {SINR}}^{\left( {k + 1} \right)} \end{align} $$ (26)

      式中,第1个不等式成立的原因是由于 ${{{h}}_{k{\rm{ + }}1}}$ ${{{s}}_k}$ 方向上的最优值,而第2个不等式成立的原因是由于 ${{{s}}_{k{\rm{ + }}1}}$ ${{{h}}_{k{\rm{ + }}1}}$ 方向上的最优值。对于不同的初始发射波形,由于收敛的方向不同,循环迭代所得 $ {{h}}$ & $ {{s}}$ “对”并不是唯一的。在交替迭代过程中,输出SINR是非递减的,所得 $ {{h}}$ & $ {{s}}$ “对”虽然不一定是全局最优点,在大部分情况下,却至少可以保证是沿着初始方向上的局部最优点。

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

    • 假设发射波形和接收滤波器的长度 ${N_{{s}}} = {N_{{h}}} =$ $ 20$ ,约束条件参数 ${P_{{s}}} = {P_{{h}}} = 20$ 。TIR和CIR的长度 ${N_{{t}}} = {N_{{c}}}{\rm{ = }}20$ , TIR的数据 ${{\tilde {t}}}$ 取自文献[23],设 ${{t}} = {\zeta _{{t}}}{{\tilde {t}}}$ , ${\zeta _{t}} = 100$ 。CIR服从均值为 ${{{{0}}}_{{N_{c}}}}$ ,协方差矩阵为 ${{{R}}_{{c}}}$ 的圆对称复高斯分布,即 ${{c}} \sim \mathcal{C}\mathcal{N}\left( {{{{{0}}}_{{N_{{c}}}}},{{{R}}_{{c}}}} \right)$ ,其中, ${{{R}}_{{c}}} \in {\mathbb{C}^{{N_{{c}}} \times {N_{{c}}}}}$ 表示维数为 ${N_{c}} \times {N_{c}}$ 的非奇异Hermitian阵。假设CIR对应的随机过程是宽平稳的,由文献[24]可知,归一化的CIR协方差矩阵可表示为 $ {{{R}}_{({{c}})}} = $ $ {{{U}}\!_{{c}}}{{{Λ}} _{{c}}}{{U}}\!_{{c}}\,\!\!^{\rm{H}}$ ,其中酉矩阵 ${{{U}}\!_c}$ k行第l 列元素为归一化的杂波协方差矩阵,可表示为 ${{{R}}_{({{c}})}} \!=\! {{{U}}\!_{{c}}}{{Λ}_{{c}}}{{U}}\!_{{c}}\,\!\!^{\rm{H}}$ ,则 ${{{R}}_0} = \sigma _{{c}}^2{{{R}}_{({{c}})}}$ ,其中 $\sigma _{c}^2 = 1$ , ${{{Λ}} _{{c}}} \in {\mathbb{C}^{{N_{{c}}} \times {N_{{c}}}}}$ 为对角阵,其对角元素为归一化杂波功率谱密度的采样,其值等于在频率间隔[0,1]上的等间隔采样,代表的是某个距离单元内杂波随时间的变化特性; ${{{U}}\!_{c}}$ 为酉矩阵,第 $k$ 行第 $l$ 列元素为:

      $$ \begin{align} &\frac{1}{{\sqrt {{N_{c}}} }}\exp\left[ {\frac{{ - {\rm j}2{{π}} \left( {k - 1} \right)\left( {{l} - 1} \right)}}{{{N_{c}}}}} \right],\\ & \quad\quad\forall k,l \in \left( {1,{N_{{c}}}} \right) \end{align} $$ (27)

      归一化的TIR和CIR协方差矩阵如图2所示。为便于分析,令噪声为高斯白噪声,其协方差矩阵为 ${{{R}}_{{n}}} = \sigma _{{n}}^2{{{I}}_{{N_{{x}}}}}$ ,其中 $\sigma _{n}^2 = 0.1$

      图  2  确定TIR和随机CIR

      Figure 2.  Determinate TIR and random CIR

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

      图  3  随迭代次数的算法性能分析

      Figure 3.  Algorithm performance analysis with iterations

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

      图  4  不同算法信干噪比随着迭代次数变化的曲线

      Figure 4.  Comparison of the SINR for different algorithms with iterations

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

      图  5  不同算法信干噪比随着CNR变化的曲线

      Figure 5.  Comparison of the SINR for different algorithms with different CNR

    • 图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波形实部和虚部的表示, ${\rm{PAR = }}{N_{s}}$ 代表能量约束下的波形,其中点的分布半径较大,说明波形幅度起伏较大,不能保证所设计波形的恒模特性,不利于实际应用。当PAR=1时,本文方法产生的点位于单位圆上,产生的是恒模波形;而当PAR=2时,点的分布半径与 ${\rm{PAR = }}{N_{s}}$ 相比较小,然而,由前面分析可知,PAR=2时波形的检测性能已经与能量约束下波形较为接近,因此低PAR下产生的波形实用性较强。

      图  6  本文算法不同PAR下的信干噪比性能对比

      Figure 6.  Comparison of the SINR vs different PAR

      图  7  不同PAR波形的实部和虚部

      Figure 7.  Real and imaginary parts of different PAR

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

参考文献 (24)

目录

    /

    返回文章
    返回