ISSN 
下载:
-
现代电子战中，各类干扰技术快速发展，从箔条类的无源干扰到噪声压制、假目标欺骗的有源干扰[1-2]，使得现代雷达的生存环境越来越复杂，严重降低了雷达的生存性能。为了保证雷达在干扰环境下正常工作，雷达抗干扰技术得到了广泛关注与研究[3-5]。
波形设计作为一种典型的雷达抗干扰手段，主要分为抗距离欺骗干扰和抗速度欺骗干扰两方面。针对距离欺骗干扰，文献[6]利用波形分集技术，在不同脉冲重复时间(pulse repetition time, PRT)发射相互正交信号，并利用匹配滤波抑制距离欺骗干扰，但是该方法采用的正交波形会导致较高的距离旁瓣。文献[7]从低截获抗干扰角度出发，通过对线性调频−频率编码复合波形设计滤波器抑制距离旁瓣，进一步提高了线性调频−频率编码复合波形的低截获抗干扰能力。文献[8]针对切片重构干扰和频谱弥散干扰导致的距离维假目标，通过最小化发射信号与干扰信号的互相关积分电平(integrated level, IL)和发射信号的自相关积分旁瓣电平(integrated sidelobe level, ISL)的加权和，在恒模约束下，设计发射波形，有效抑制了这两种干扰类型。文献[9]针对间歇采样干扰，通过建立最小化目标信号与接收滤波器的ISL和干扰信号与接收滤波器的IL的加权和的优化准则，在PAR约束下，采用交替方向乘子法(alternation direction method of multipliers, ADMM)联合设计发射波形与接收滤波器实现间歇采样干扰抑制，但是该算法无法保证收敛性。另一方面，针对速度欺骗干扰，文献[10]通过设计脉冲初相有效地对抗速度欺骗干扰。该方法考虑恒模约束，建立最小化干扰阻带内能量的优化问题，利用共轭梯度法求解最优脉间初相，极大地改善了脉冲多普勒雷达在速度欺骗干扰环境下的目标探测性能。在此研究基础上，文献[11-12]分别利用修正牛顿法和WeCAN算法求解文献[10]中的优化问题。但是，文献[11]的工作仅涉及脉间初相，未考虑脉间幅度设计，没有充分利用幅度这一自由度；同时，相比离散量化相位，连续相位设计在实际工程应用中较为困难。
对此，本文增加脉冲幅度这一自由度，同时考虑离散相位，研究脉间波形幅相联合设计抗速度欺骗干扰方法。通过将最小化阻带内干扰能量与目标旁瓣能量的加权和作为优化准则，在离散量化相位和PAR约束下，构建幅度−相位联合设计抗速度欺骗干扰的优化问题。然后，根据交替方向惩罚法(alternating direction penalty method, ADPM)框架，提出IADPM算法求解该优化问题，具体地，通过引入辅助变量将原优化问题分解为一系列子问题，每个子问题利用不同算法求出近似解，进而不断循环迭代，求出最优幅度和相位，最终实现脉间幅度相位联合设计抗速度欺骗干扰。
-
假设脉冲多普勒雷达在一个CPI内发射
$N$ 个不同幅度和初相的脉冲信号，第$n$ 个脉冲发射信号为：$${s_n}\left( t \right) = {x_n}{{\rm{e}}^{{\rm{j}}{\theta _n}}}u\left( {t - n{T_r}} \right)\quad n = 0,1, \cdots ,N - 1$$ (1) 式中，
${x_n}$ 和${\theta _n}$ 分别为第$n$ 个脉冲信号的幅度和初相；${T_r}$ 为脉冲重复时间；$u\left( t \right)$ 定义为：$$u\left( t \right) = \left\{ \begin{aligned} &1\quad 0 < t \leqslant {T_p} \\ & 0\quad {\text{其他}} \end{aligned} \right.$$ (2) 式中，
${T_p}$ 为脉冲宽度。假设在某一距离单元内，有
$Q$ 个不同多普勒频率的目标，因此，第$n$ 个脉冲接收到的目标回波为：$${T_n}\left( t \right) = {x_n}{{\rm{e}}^{{\rm{j}}{\theta _n}}}\sum\limits_{q = 1}^Q {{\alpha _{Tq}}\exp \left( {{\rm{j}}2{{\text π}}n{f_{Tq}}} \right)u\left( {t - n{T_r} - {\tau _0}} \right)} $$ (3) 式中，
${\tau _0}$ 为目标双程时延；${\alpha _{Tq}}$ 为第$q$ 个目标的回波幅度；${f_{Tq}} = {{2{v_{Tq}}{T_r}} / \lambda }$ 为第$q$ 个目标的归一化多普勒频率；${v_{Tq}}$ 为相应的目标速度；$\lambda $ 为雷达工作波长。数字射频存储器(digital radio frequency memory, DRFM)干扰机截获到雷达信号后需要对信号进行分选识别，调制虚假多普勒信息形成干扰信号后转发出去，一般假设DRFM干扰机发射的脉冲信号滞后真实雷达信号
$i$ 个PRT[10]，如图1所示。
图 1 干扰脉冲示意图
假设在时延
${\tau _0}$ 处，DRFM干扰机调制了$P$ 个不同多普勒频率的虚假目标，则第$n$ 个脉冲接收到的干扰回波为：$$\begin{split} & \qquad\qquad\qquad {J_n}\left( t \right) = {x_{n - i}}{{\rm{e}}^{{\rm{j}}{\theta _{n - i}}}} \times \\ & \sum\limits_{p = 1}^P {{\alpha _{Jp}}\exp \left( {{\rm{j}}2{{\text π}}\left( {n - i} \right){f_{Jp}}} \right)u\left( {t - \left( {n - i} \right){T_r} - {\tau _0}} \right)} \end{split} $$ (4) 式中，
${\alpha _{Jp}}$ 和${f_{Jp}}$ 分别为第$p$ 个虚假目标的幅度和归一化多普勒频率。因此，第
$n$ 个脉冲的接收回波为：$$\begin{split} & \qquad\;\; {y_n}\left( t \right) = {T_n}\left( t \right) + {J_n}\left( t \right) + {\omega _n}\left( t \right) =\\ & {x_n}{{\rm{e}}^{{\rm{j}}{\theta _n}}}\sum\limits_{q = 1}^Q {{\alpha _{Tq}}\exp \left( {{\rm{j}}2{{\text π}}n{f_{Tq}}} \right)u\left( {t - n{T_r} - {\tau _0}} \right)} + \\ & \qquad {x_{n - i}}{{\rm{e}}^{{\rm{j}}{\theta _{n - i}}}}\sum\limits_{p = 1}^P {\alpha _{Jp}}\exp \left( {{\rm{j}}2{{\text π}}\left( {n - i} \right){f_{Jp}}} \right) \times \\ &\qquad\qquad u\left( {t - \left( {n - i} \right){T_r} - {\tau _0}} \right) + {\omega _n}\left( t \right) \end{split} $$ (5) 式中，
${\omega _n}\left( t \right)$ 表示加性高斯白噪声信号。经过匹配滤波处理后，在时延
${\tau _0}$ 处，第$n$ 个脉冲的输出为：$$\begin{split} & \qquad\quad\;\; y\left( n \right) = {x_n}\sum\limits_{q = 1}^Q {{\alpha _{Tq}}\exp \left( {{\rm{j}}2{{\text π}}n{f_{Tq}}} \right)}+ \\ & {x_{n - i}}{{\rm{e}}^{{\rm{j}}\left( {{\theta _{n - i}} - {\theta _n}} \right)}}\sum\limits_{p = 1}^P {{\alpha _{Jp}}\exp \left( {{\rm{j}}2{{\text π}}\left( {n - i} \right){f_{Jp}}} \right)} + \upsilon \left( n \right) \end{split} $$ (6) 式中，
$\upsilon \left( n \right)$ 为${\omega _n}\left( t \right)$ 的滤波输出值，满足均值为0，方差为${\sigma ^2}$ 的高斯分布。 -
利用文献[10]中的多通道处理技术，获取目标和干扰信号的多普勒先验信息。假设干扰信号与目标信号的多普勒频率均在阻带范围内，设置阻带范围为：
$$\Omega = \mathop { \mathop \bigcup \limits^{{N_s}} }\limits_{k = 1} \left[ {{f_{k1}},{f_{k2}}} \right]$$ (7) 式中，
$ \bigcup\limits $ 为取并集操作；${N_s}$ 为阻带个数；${f_{k1}}$ 和${f_{k2}}$ 分别为第$k$ 个阻带的下限和上限频率。根据文献[13]，忽略干扰回波幅度
${\alpha _{Jp}}$ ，可以推导出$P$ 个干扰信号在阻带内的能量为：$${C_J} = {{{s}}^{\rm{H}}}{{{R}}_J}{{s}}$$ (8) 式中，
$ {(\cdot)}^{{\rm{H}}}$ 表示共轭转置操作；${{s}}$ 表示为：$$\begin{split} &\qquad\quad {{s}} = {\left[ {{s_0},{s_1}, \cdots ,{s_{N - 1}}} \right]^{\rm{T}}} =\\ & \left[ {{x_{ - i}}{{\rm{e}}^{{\rm{j}}{\varphi _0}}},{x_{1 - i}}{{\rm{e}}^{{\rm{j}}{\varphi _1}}},} \right. \cdots ,{\left. {{x_{N - 1 - i}}{{\rm{e}}^{{\rm{j}}{\varphi _{N - 1}}}}} \right]^{\rm{T}}} \end{split} $$ (9) 式中，
$ {(\cdot)}^{{\rm{T}}}$ 为转置操作。矩阵${{{R}}_J}$ 的第$\left( {m,n} \right)$ 个元素表示为：$$\begin{split} & {{{R}}_J}\left( {m,n} \right) = \sum\limits_{k = 1}^{{N_s}} {\sum\limits_{p = 1}^P {{{\rm{e}}^{ - {\rm{j2{\text π} }}\left( {m - n} \right)f_{Jp}^{}}}\int_{{f_{k1}}}^{{f_{k2}}} {{{\rm{e}}^{{\rm{j2{\text π} }}f\left( {m - n} \right)}}{\rm{d}}f} } } = \\ & \left\{ \begin{array}{l} \displaystyle\sum\limits_{k = 1}^{{N_s}} {\displaystyle\sum\limits_{p = 1}^P {\frac{{{{\rm{e}}^{{\rm{j2{\text π} }}\left( {{f_{k2}} - f_{Jp}^{}} \right)\left( {m - n} \right)}} - {{\rm{e}}^{{\rm{j2{\text π} }}\left( {{f_{k1}} - f_{Jp}^{}} \right)\left( {m - n} \right)}}}}{{{\rm{j}}2{{\text π}}\left( {m - n} \right)}}} } \quad m \ne n \\ P\displaystyle\sum\limits_{k = 1}^{{N_s}} ({{f_{k2}} - {f_{k1}}}) \qquad\qquad\qquad\qquad m = n \end{array} \right. \end{split} $$ (10) 相似地，
$Q$ 个目标信号在阻带内的旁瓣能量可以推导为：$${C_T} = {{{x}}^{\rm{T}}}\left( {{{{R}}_{T1}} + {{{R}}_{T2}}} \right){{x}} = {{{x}}^{\rm{T}}}{{{R}}_T}{{x}}$$ (11) 式中，
${{x}} = {\left[ {{x_0},{x_1}, \cdots ,{x_{N - 1}}} \right]^{\rm{T}}}$ ；${{{R}}_{T1}}$ 和${{{R}}_{T2}}$ 中的第$\left( {m,n} \right)$ 个元素分别表示为：$$\begin{split} & {{{R}}_{T1}}\left( {m,n} \right) = \sum\limits_{k = 1}^{{N_s}} {\sum\limits_{q = 1}^Q {{{\rm{e}}^{ - {\rm{j2{\text π} }}\left( {m - n} \right)f_{Tq}^{}}}\int_{{f_{k1}}}^{{f_{Tq}} - \Delta f} {{{\rm{e}}^{{\rm{j2{\text π} }}f\left( {m - n} \right)}}{\rm{d}}f} } } = \\ & \left\{ \begin{array}{l} \displaystyle\sum\limits_{k = 1}^{{N_s}} {\displaystyle\sum\limits_{q = 1}^Q {\frac{{{{\rm{e}}^{{\rm{j2{\text π} }}\left( { - \Delta f} \right)\left( {m - n} \right)}} - {{\rm{e}}^{{\rm{j2{\text π} }}\left( {{f_{k1}} - f_{Tq}^{}} \right)\left( {m - n} \right)}}}}{{{\rm{j}}2{\rm{{\text π} }}\left( {m - n} \right)}}} } \quad m \ne n \\ \displaystyle\sum\limits_{k = 1}^{{N_s}} {\displaystyle\sum\limits_{q = 1}^Q ({{f_{Tq}} - \Delta f - {f_{k1}}} }) \qquad\qquad\qquad\qquad m = n \end{array} \right. \\[-33pt] \end{split} $$ (12) $$\begin{split} & {{{R}}_{T2}}\left( {m,n} \right) = \sum\limits_{k = 1}^{{N_s}} {\sum\limits_{q = 1}^Q {{{\rm{e}}^{ - {\rm{j2{\text π} }}\left( {m - n} \right)f_{Tq}^{}}}\int_{{f_{Tq}} + \Delta f}^{{f_{k2}}} {{{\rm{e}}^{{\rm{j2{\text π} }}f\left( {m - n} \right)}}{\rm{d}}f} } } = \\ & \left\{ \begin{array}{l} \displaystyle\sum\limits_{k = 1}^{{N_s}} {\displaystyle\sum\limits_{q = 1}^Q {\frac{{{{\rm{e}}^{{\rm{j2{\text π} }}\left( {{f_{k2}} - f_{Tq}^{}} \right)\left( {m - n} \right)}} - {{\rm{e}}^{{\rm{j2{\text π} }}\Delta f\left( {m - n} \right)}}}}{{{\rm{j2{\text π} }}\left( {m - n} \right)}}} } \quad m \ne n \\ \displaystyle\sum\limits_{k = 1}^{{N_s}} {\displaystyle\sum\limits_{q = 1}^Q ({{f_{k2}} - {f_{Tq}} - \Delta f} } )\qquad\qquad\qquad m = n \end{array} \right. \\[-33pt] \end{split} $$ (13) 式中，
$\Delta f = {1 / N}$ 为归一化多普勒频率间隔。令：
$$\begin{split} & \qquad\qquad\qquad {{\varphi }} = {\left[ {{\varphi _0},{\varphi _1}, \cdots ,{\varphi _{N - 1}}} \right]^{\rm{T}}} \\ & {{a}} = {\left[ {{a_0},{a_1}, \cdots ,{a_{N - 1}}} \right]^{\rm{T}}} = {\left[ {{x_{ - i}},{x_{1 - i}}, \cdots ,{x_{N - 1 - i}}} \right]^{\rm{T}}} \end{split} $$ (14) 因此有：
$${{{C}}_J} = {{{s}}^{\rm{H}}}{{{R}}_J}{{s}} = {\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)^{\rm{H}}}{{{R}}_J}\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)$$ (15) $${{{C}}_T} = {{{x}}^{\rm{T}}}{{{R}}_T}{{x}} = {\left( {{{J}}\left| {{s}} \right|} \right)^{\rm{T}}}{{{R}}_T}{{J}}\left| {{s}} \right| = {\left( {{{Ja}}} \right)^{\rm{T}}}{{{R}}_T}{{Ja}}$$ (16) 式中，
$ |\cdot|$ 为取模操作；$ \odot $ 为哈达玛积；${{J}}$ 为循环移位矩阵，定义为：$${{J}} = \left[ {\begin{array}{*{20}{c}} {{{{{\text{0}}}}_{\left( {N - i} \right) \times i}}}&{{I}}_{\left( {N - i} \right) \times {\left( {N - i} \right)}} \\ {{{{I}}_{i\times i}}}&{{{{{\text{0}}}}_{i \times \left( {N - i} \right)}}} \end{array}} \right]$$ (17) 为了抑制欺骗干扰，需要使干扰信号在阻带范围内能量尽可能小，同时也要保证目标信号在阻带范围内的旁瓣能量尽可能小。同时，为了避免波形的非线性失真，考虑PAR约束，可以构造以下优化问题：
$${{\cal{P}}_0}\left\{ \begin{array}{l} \mathop {\min }\limits_{{s}} f\left( {{s}} \right) = {\mu _1}{{{s}}^{\rm{H}}}{{{R}}_J}{{s}} + {\mu _2}{\left( {{{J}}\left| {{s}} \right|} \right)^{\rm{H}}}{{{R}}_T}\left( {{{J}}\left| {{s}} \right|} \right) \\ {\rm{s.t.}}{\rm{PAR}}\left( {{s}} \right) \leqslant {\gamma ^2} \\ \quad \arg \,{s_n} \in \left\{ {0,\dfrac{1}{M}2{{\text π}}, \cdots ,\dfrac{{M - 1}}{M}2{{\text π}}} \right\}\quad \forall n \end{array} \right.$$ (18) 式中，
$\,{\mu _1}$ 和$\,{\mu _2}$ 为权重因子，满足$\,{\mu _1} + {\mu _2} = 1$ ，用来折中干扰能量和目标旁瓣能量；$ \mathrm{arg}(\cdot)$ 表示取相位操作；${s_n}$ 为矢量${{s}}$ 中的第$n$ 个元素；$M$ 为离散相位个数；$\gamma \in \left[ {1,\sqrt N } \right]$ 控制最大PAR值；其中${\rm{PAR}}\left( {{s}} \right)$ 定义为[14]：$${\rm{PAR}}\left( {{s}} \right) = \frac{{\mathop {\max }\limits_n {{\left| {{s_n}} \right|}^2}}}{{\dfrac{1}{N}\displaystyle\sum\limits_{n = 0}^{N - 1} {{{\left| {{s_n}} \right|}^2}} }}$$ (19) -
本节提出了IADPM算法求解问题
${{\cal{P}}_0}$ 。该算法的核心思想是通过引入辅助变量将原问题转换为一系列子问题，通过对每个子问题求取非精确解并不断交替循环迭代，最终求解包含幅度与相位的序列${{s}}$ 。 -
将式(15)～式(16)带入问题
${{\cal{P}}_0}$ 中，问题${{\cal{P}}_0}$ 可进一步写为：$${{\cal{P}}_1}\left\{ \begin{split} & \mathop {\min }\limits_{{{a}},{{\varphi }}} {\mu _1}{\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)^{\rm{H}}}{{{R}}_J}\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right) + {\mu _2}{\left( {{{Ja}}} \right)^{\rm{H}}}{{{R}}_T}\left( {{{Ja}}} \right) \\ & {\rm{s.t.}}\quad 0 \leqslant {a_n} \leqslant \gamma \quad \forall n \\ & \quad {{a}} \in {\mathbb{R}^N} \\ & \quad {\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)^{\rm{H}}}\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right) = N \\ & \quad {\varphi _n} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\} \end{split} \right.$$ (20) 式中，
${\mathbb{R}^N}$ 为$N$ 维实向量集合。基于文献[15]中ADPM框架，本小节提出了IADPM算法求解问题
${{\cal{P}}_1}$ 。首先，引入辅助变量${{z}}$ ，问题${{\cal{P}}_1}$ 可重写为：$${{\cal{P}}_2}\left\{ \begin{aligned} & \mathop {\min }\limits_{{{a}},{{\varphi }}} {\mu _1}{\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)^{\rm{H}}}{{{R}}_J}\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)\; + {\mu _2}{\left( {{{Ja}}} \right)^{\rm{H}}}{{{R}}_T}\left( {{{Ja}}} \right) \\ & {\rm{ s.t.}}\quad{{z}} = {{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}} \\ & \quad \,\left| {{z_n}} \right| \leqslant \gamma ,\forall n \\ & \quad {{{z}}^{\rm{H}}}{{z}} = N \\ & \quad {\varphi _n} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\} \end{aligned} \right.$$ (21) 进而，构造增广拉格朗日函数：
$$\begin{split} & L\left( {{{a}},{{\varphi }},{{z}}} \right) = {\mu _1}{\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)^{\rm{H}}}{{{R}}_J}\left( {{{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)+ \\ & \qquad\qquad {\mu _2}{\left( {{{Ja}}} \right)^{\rm{H}}}{{{R}}_T}\left( {{{Ja}}} \right) +\\ & \Re \left\{ {{{{\lambda }}^H}\left( {{{z}} - {{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right)} \right\} + \frac{\rho }{2}{\left\| {{{z}} - {{a}} \odot {{\rm{e}}^{{\rm{j}}{{\varphi }}}}} \right\|^2} \end{split} $$ (22) 式中，
${{\lambda }}$ 为拉格朗日乘子；$\rho $ 为惩罚因子；$ \Re \{\cdot\}$ 为取实部操作；$ \Vert \cdot\Vert $ 表示二范数。根据ADPM框架，通过交替更新
${{z}},{{a}},{{\varphi }}$ 来最小化$L\left( {{{a}},{{\varphi }},{{z}}} \right)$ ，进而求得原问题的解。假设${{z}},{{a}},{{\varphi }}$ 在第$t$ 次迭代的解分别为${{{z}}^{\left( t \right)}},{{{a}}^{\left( t \right)}},{{{\varphi }}^{\left( t \right)}}$ ，则提出的IADPM算法的迭代规则如下：$${{{z}}^{\left( t \right)}} = \arg \mathop {\min}\limits_{{{z}} \in {\Omega _{{z}}}} L\left( {{{z}},{{{a}}^{\left( {t - 1} \right)}},{{{\varphi }}^{\left( {t - 1} \right)}},{{\lambda }}_{}^{\left( {t - 1} \right)},{\rho ^{\left( {t - 1} \right)}}} \right)$$ (23) $${\varOmega _{{z}}} = \left\{ \begin{array}{l} \left| {{z_n}} \right| \leqslant \gamma ,{{{z}}^H}{{z}} = N \\ \arg {z_n} \in \left\{ {0,\dfrac{1}{M}2{{\text π}}, \cdots ,\dfrac{{M - 1}}{M}2{{\text π}}} \right\},\forall n \end{array} \right\}$$ (24) $${{{a}}^{\left( t \right)}} = \arg \mathop {\min }\limits_{{\Omega _{{a}}}} {\mkern 1mu} L\left( {{{{z}}^{\left( t \right)}},{{a}},{{{\varphi }}^{\left( {t - 1} \right)}},{{\lambda }}_{}^{\left( {t - 1} \right)},{\rho ^{\left( {t - 1} \right)}}} \right)$$ (25) $${\varOmega _{{a}}} = \left\{ {{{a}} \in {\mathbb{R}^N},0 \leqslant {a_n} \leqslant \gamma ,\forall n} \right\}$$ (26) $${{{\varphi }}^{\left( t \right)}} = \arg \mathop {\min }\limits_{{{\varphi }} \in {\Omega _{{\varphi }}}} {\mkern 1mu} L\left( {{{{z}}^{\left( t \right)}},{{{a}}^{\left( t \right)}},{{\varphi }},{{\lambda }}_{}^{\left( {t - 1} \right)},{\rho ^{\left( {t - 1} \right)}}} \right)$$ (27) $${\varOmega _{{\varphi }}} = \left\{ {{\varphi _n} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\},\forall n} \right\}$$ (28) $$\rho _{}^{\left( t \right)} = \left\{ \begin{aligned} & \rho _{}^{\left( {t - 1} \right)},{\rm{if}}\left\| {{{{z}}^{\left( t \right)}} - \left( {{{{a}}^{\left( t \right)}} \odot {{\rm{e}}^{{\rm{j}}{{{\varphi }}^{\left( t \right)}}}}} \right)} \right\| \leqslant \\ & \quad {{\rm{\delta }}_{{\rm{1c}}}}\left\| {{{{z}}^{\left( {t - 1} \right)}} - \left( {{{{a}}^{\left( {t - 1} \right)}} \odot {{\rm{e}}^{{\rm{j}}{{{\varphi }}^{\left( {t - 1} \right)}}}}} \right)} \right\| \\ & \rho _{}^{\left( {t - 1} \right)}{{\rm{\delta }}_{{\rm{2c}}}}\quad {\rm{else}} \end{aligned} \right.$$ (29) $${{\lambda }}_{}^{\left( t \right)} = {{\lambda }}_{}^{\left( {t - 1} \right)} + \rho _{}^{\left( t \right)}\left( {{{{z}}^{\left( t \right)}} - \left( {{{{a}}^{\left( t \right)}} \odot {{\rm{e}}^{{\rm{j}}{{{\varphi }}^{\left( t \right)}}}}} \right)} \right)$$ (30) 式中，
${{\rm{\delta }}_{{\rm{1c}}}}$ 和${{\rm{\delta }}_{{\rm{2c}}}}$ 为接近1的两个常数，且满足$0 < {{\rm{\delta }}_{{\rm{1c}}}} < 1$ 和${{\rm{\delta }}_{{\rm{2c}}}} > 1$ 。首先，固定
${{{a}}^{\left( {t - 1} \right)}},\;{{{\varphi }}^{\left( {t - 1} \right)}},\;{{\lambda }}_{}^{\left( {t - 1} \right)},\;{\rho ^{\left( {t - 1} \right)}}$ ，更新${{{z}}^{\left( t \right)}}$ 。忽略与${{z}}$ 的无关项，式(23)可以重写为：$$\begin{split} & \mathop {\min }\limits_{{z}} \frac{\rho }{2}{{{z}}^{\rm{H}}}{{z}} + \Re \left( {{{{z}}^{\rm{H}}}{{\lambda }} - \rho {{{z}}^{\rm{H}}}{{s}}} \right) \\ & \quad\quad\quad\;{\rm{s.t.}}\quad\left| {{z_n}} \right| \leqslant \gamma ,\forall n \\ & \quad\quad\quad\quad{\left\| {{z}} \right\|^2} = N \\ & \arg {z_n} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\} \end{split} $$ (31) 式(31)可以用经典的Power Method[16]方法求解，首先定义：
$${{b}} = {{s}} - {{{\lambda }} / \rho }$$ (32) ${z_n}$ 的相位满足$\arg {z_n} = \arg b_n^{}$ ，由于本文考虑离散量化相位，则选取最接近最优值的离散量化相位即可。接下来计算$\left| {{z_n}} \right|$ ，假设矢量${{b}}$ 中的元素满足$\left| {b_1^{}} \right| \geqslant \cdots \geqslant \left| {b_N^{}} \right|$ ，其非零元素个数为$m$ ，当$m{\gamma ^2} \leqslant N$ 时，式(31)的解为：$$\left| {{z_n}} \right| = \left\{ \begin{split} & \gamma \quad n = 1,2, \cdots ,m \\ & \sqrt {\frac{{N - m{\gamma ^2}}}{{N - m}}} \quad n = m + 1,m+2, \cdots ,N \end{split} \right.$$ (33) 当
$m{\gamma ^2} > N$ 时，式(31)的解为：$$\left| {{z_n}} \right| = \min \left\{ {\beta \left| {b_n^{}} \right|,\gamma } \right\}\quad n = 1,2, \cdots ,N$$ (34) 式中，
$$\beta \in \left\{ {\beta \left| \begin{array}{l} \displaystyle\sum\limits_n {\min \left\{ {{\beta ^2}{{\left| {b_n^{}} \right|}^2},{\gamma ^2}} \right\}} - N = 0 \\ \beta \in \left[ {0,\dfrac{\gamma }{{\min \left\{ {\left| {b_n^{}} \right|\left| {,\left| {b_n^{}} \right| \ne 0} \right.} \right\}}}} \right] \end{array} \right.} \right\}$$ (35) 由于函数
$f\left( \beta \right) = \displaystyle\sum\limits_n {\min \left\{ {{\beta ^2}{{\left| {{b_n}} \right|}^2},{\gamma ^2}} \right\}} - N$ 在区间上单调递增，且$f\left( 0 \right) \leqslant 0$ ，因此可用二分法求得$\,\beta $ 。接着，固定
${{{z}}^{\left( t \right)}},{{{\varphi }}^{\left( {t - 1} \right)}},{{\lambda }}_{}^{\left( {t - 1} \right)},{\rho ^{\left( {t - 1} \right)}}$ ，更新${{{a}}^{\left( t \right)}}$ 。忽略${{a}}$ 的无关项，式(25)可以写为：$$\begin{split} &\quad \; \mathop {\min }\limits_{{a}} {\mu _1}{\left( {{{Da}}} \right)^{\rm{H}}}{{{R}}_J}{{Da}} + {\mu _2}{{{a}}^{\rm{H}}}{{{J}}^{\rm{H}}}{{{R}}_T}{{Ja}} - \\ & \Re \left\{ {{{{\lambda }}^{\rm{H}}}{{Da}}} \right\} - \rho \Re \left\{ {{{{z}}^{\rm{H}}}{{Da}}} \right\} + \frac{\rho }{2}{{{a}}^{\rm{H}}}{{{D}}^{\rm{H}}}{{Da}} \\ & \quad \;\quad \;\quad \; {\rm{s.t.}}\quad{{a}} \in {{\mathbb{R}}_ {a_n}^N} \geqslant 0,\forall n \end{split} $$ (36) 式中，矩阵
${{D}}$ 定义为：$${{D}} = {\rm{diag}}\left( {{{\rm{e}}^{{\rm{j}}{\varphi _0}}},{{\rm{e}}^{{\rm{j}}{\varphi _1}}}, \cdots ,{{\rm{e}}^{{\rm{j}}{\varphi _{N - 1}}}}} \right)$$ (37) 式中，
${\rm{diag}}\left( \cdots \right)$ 表示将数据形成对角矩阵。利用坐标下降(coordinate descent, CD)[17]算法求解
${{a}}$ ，固定${a_0}, \cdots ,{a_{p - 1}},{a_{p + 1}}, \cdots ,{a_{N - 1}}$ ，式(36)中目标函数可以表示为关于${a_p}$ 的表达式为：$$\begin{split} & {\mu _1}{\left( {{{Da}}} \right)^{\rm{H}}}{{{R}}_J}{{Da}} + {\mu _2}{{{a}}^{\rm{H}}}{{{J}}^{\rm{H}}}{{{R}}_T}{{Ja}} - \Re \left\{ {{{{\lambda }}^{\rm{H}}}{{Da}}} \right\}- \\ & \quad \rho \Re \left\{ {{{{z}}^{\rm{H}}}{{Da}}} \right\} + \frac{\rho }{2}{{{a}}^{\rm{H}}}{{{D}}^{\rm{H}}}{{Da}} = ua_p^2 + v{a_p} + l \end{split} $$ (38) 式中，
$$\begin{split} &\;\; v = {\mu _1}\left( {{{\rm{e}}^{ - {\rm{j}}{\varphi _p}}}\sum\limits_{n \ne p} {{a_n}{{\rm{e}}^{{\rm{j}}{\varphi _n}}}{r_{pn}}} + {{\rm{e}}^{{\rm{j}}{\varphi _p}}}\sum\limits_{n \ne p} {{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}{r_{np}}} } \right) +\\ & {\mu _2}\left( {\sum\limits_{n \ne p} {{a_n}{{r'}_{pn}}} + \sum\limits_{n \ne p} {{a_n}{{r'}_{np}}} } \right) - \Re \left( {{{\bar \lambda }_p}{{\rm{e}}^{{\rm{j}}{\varphi _p}}}} \right) - \rho \Re \left( {{{\bar z}_p}{{\rm{e}}^{{\rm{j}}{\varphi _p}}}} \right) \end{split} $$ (39) $$u = {\mu _1}{r_{pp}} + {\mu _2}{r'_{pp}} + 0.5\rho $$ (40) $$\begin{array}{l} l = {\mu _1}\displaystyle\sum\limits_{m \ne p} {\displaystyle\sum\limits_{n \ne p} {{a_m}{{\rm{e}}^{{\rm{j}}{\varphi _m}}}{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}{r_{nm}}} } + {\mu _2}\displaystyle\sum\limits_{m \ne p} {{a_m}\displaystyle\sum\limits_{n \ne p} {{a_n}{{r'}_{nm}}} } - \\ \Re \left( {\displaystyle\sum\limits_{n \ne p} {{{\bar \lambda }_n}{a_n}{{\rm{e}}^{{\rm{j}}{\varphi _n}}}} } \right) - \rho \Re \left( {\displaystyle\sum\limits_{n \ne p} {{{\bar z}_n}{a_n}{{\rm{e}}^{{\rm{j}}{\varphi _n}}}} } \right) + \displaystyle\sum\limits_{n \ne p} {a_n^2}\\[-22pt] \end{array} $$ (41) 式中，
${r_{nm}}$ 和${r'_{nm}}$ 分别为矩阵${{{R}}_J}$ 和${{{R}}_T}$ 中的第$\left( {n,m} \right)$ 个元素。$ \overline{(\cdot)}$ 表示取共轭操作，由于$u > 0$ ，可得${a_p}$ 最优解为：$$a_p^{} = \left\{ \begin{aligned} & - \dfrac{v}{{2u}}\quad \; - \dfrac{v}{{2u}} > 0 \\ & 0\quad\quad\quad {\text{其他}} \end{aligned} \right.$$ (42) 最后，固定
${{{z}}^{\left( t \right)}},{{{a}}^{\left( t \right)}},{{\lambda }}_{}^{\left( {t - 1} \right)},{\rho ^{\left( {t - 1} \right)}}$ ，更新${{{\varphi }}^{\left( t \right)}}$ 。忽略${{\varphi }}$ 的无关项，关于${{\varphi }}$ 的优化问题可写为：$$\begin{split} & \mathop {\min }\limits_{{\varphi }} {\mu _1}{\left( {{{Da}}} \right)^{\rm{H}}}{{{R}}_J}{{Da}} - \Re \left\{ {{{{\lambda }}^{\rm{H}}}{{Da}}} \right\} - \rho \Re \left\{ {{{{z}}^{\rm{H}}}{{Da}}} \right\} \\ & \qquad {\rm{s.t.}}\quad {\varphi _n} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\}\quad \forall n \end{split} $$ (43) 同样采取CD算法求解
${{\varphi }}$ 。固定${\varphi _0}, \cdots ,{\varphi _{p - 1}},$ ${\varphi _{p + 1}}, $ $ \cdots ,{\varphi _{N - 1}}$ ，式(38)关于${\varphi _p}$ 的表达式可以重写为：$$ {\mu _1}{\left( {{{Da}}} \right)^{\rm{H}}}{{RDa}} \!-\! \Re \left\{ {{{{\lambda }}^{\rm{H}}}{{Da}}} \right\} \!-\! \rho \Re \left\{ {{{{z}}^{\rm{H}}}{{Da}}} \right\} = c + \Re \left( {{a_p}{{\rm{e}}^{{\rm{j}}{\varphi _p}}}g} \right) $$ (44) 式中，
$$\begin{split} & c = {\mu _1}\left( {\sum\limits_{m \ne p} {\sum\limits_{n \ne p} {{a_m}{{\rm{e}}^{{\rm{j}}{\varphi _m}}}{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}{r_{nm}}} } + x_p^2{r_{pp}}} \right) -\\ & \qquad \frac{1}{2}\sum\limits_{n \ne p} {{{\bar \lambda }_n}{x_n}{{\rm{e}}^{{\rm{j}}{\varphi _n}}}} - \frac{1}{2}\sum\limits_{n \ne p} {{\lambda _n}{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}} - \\ & \qquad \frac{\rho }{2}\sum\limits_{n \ne p} {{{\bar z}_n}{a_n}{{\rm{e}}^{{\rm{j}}{\varphi _n}}}} - \frac{\rho }{2}\sum\limits_{n \ne p} {{z_n}{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}} \end{split} $$ (45) $$g = 2{\mu _1}\sum\limits_{n \ne p} {{a_n}{{\rm{e}}^{ - {\rm{j}}{\varphi _n}}}{r_{np}}} - {\bar \lambda _p} - \rho {\bar z_p}$$ (46) 因此，式(44)关于
${\varphi _p}$ 的优化问题可以重写为：$$\begin{split} & \qquad\quad \mathop {\min }\limits_{{\varphi _p}} \Re \left( {{a_p}{{\rm{e}}^{{\rm{j}}{\varphi _p}}}g} \right) \\ & {\rm{s.t.}}\quad{\varphi _p} \in \left\{ {0,\frac{1}{M}2{{\text π}}, \cdots ,\frac{{M - 1}}{M}2{{\text π}}} \right\} \end{split} $$ (47) 可以求得：
$${\varphi _p} = f\left( {{{\text π}} - \arg g} \right)$$ (48) 式中，
$f\left( \psi \right)$ 表示最接近$\psi $ 的离散量化相位值。${{z}},{{a}},{{\varphi }},{{\lambda }},\rho $ 不断迭代更新直到收敛。根据文献[18]，设置停止迭代的条件为：$$ {\epsilon}^{\left(t\right)}=\Vert {{{z}}}^{\left(t\right)}-{{{a}}}^{\left(t\right)}\odot {\rm{e}}^{{\rm{j}}{{{\varphi}} }^{\left(t\right)}}\Vert \leqslant \zeta $$ (49) IADPM算法求解问题
${{\cal{P}}_0}$ ，总结IADPM算法流程如下所示。输入：初始值
${{{z}}^{\left( 0 \right)}},{{{a}}^{\left( 0 \right)}},{{{\varphi }}^{\left( 0 \right)}},{{{\lambda }}^{\left( 0 \right)}},{\rho ^{\left( 0 \right)}}$ ，${{\rm{\delta }}_{{\rm{1c}}}},{{\rm{\delta }}_{{\rm{2c}}}}$ 和控制收敛的阈值${\rm{\zeta }}$ 输出：局部最优解
${{a}}$ 和${{\varphi }}$ ；while
$ {\epsilon}^{\left(t\right)}>\zeta $ 判断
$m{\gamma ^2}$ 是否大于$N$ ，根据式(33)和式(34)更新${{{z}}^{\left( t \right)}}$ ；for p=1:N
根据式(42)更新
$a_p^{\left( t \right)}$ ；end
for p= 1:N
根据式(48)更新
$\varphi _p^{\left( t \right)}$ ；end
分别根据式(29)和式(30)更新
${\rho ^{\left( t \right)}}$ 和${{{\lambda }}^{\left( t \right)}}$ ；根据式(49)计算
$ {\epsilon}^{\left(t\right)}$ ；end while
求出最优序列
${{a}}$ 和${{\varphi }}$ 后，根据${{a}}$ 和${{\varphi }}$ 的定义，即可计算发射脉冲串的幅度和相位序列。 -
基于IADPM算法求解问题
${{\cal{P}}_0}$ 的计算量主要与迭代次数和矢量${{z}},{{a}},{{\varphi }}$ 的更新有关。其中，更新矢量${{z}}$ 所用的Power Method计算量主要来自于二分法求$\beta $ 值，相应地计算量为$O\left( {I\left( N \right)} \right)$ ，其中$I$ 为二分法迭代次数；更新矢量${{a}}$ 中的计算量主要来自于计算$u$ 和$v$ ，其计算量为$O\left( N \right)$ ，遍历$N$ 个元素，其计算量为$O\left( {{N^2}} \right)$ ；更新矢量${{\varphi }}$ 中的计算量主要来自于构造$g$ ，其计算量为$O\left( N \right)$ ，遍历完$N$ 个元素，其计算量为$O\left( {{N^2}} \right)$ 。因此，在IADPM算法一次迭代中的计算量为$O\left( {I\left( N \right)} \right) + O\left( {{N^2}} \right)$ 。 -
本节主要通过仿真实验验证脉间幅度相位捷变设计抗速度欺骗干扰方法的有效性。
-
本小节首先对提出的IADPM算法的收敛性进行分析。考虑单目标与单干扰场景下，即
$P = 1$ ，$Q = 1$ ，假设雷达发射脉冲个数$N = 512$ ，干扰的归一化多普勒频率${f_J} = 0.4$ ，目标的归一化多普勒频率${f_T} = 0.47$ ，干扰信号滞后目标信号2个PRT，设置阻带范围为$\varOmega = \left[ {{f_{11}},{f_{12}}} \right] = \left[ {0.3,0.5} \right]$ 。随机初始化${{s}},{{\lambda }}$ ,${{z}}$ 的初始值与${{s}}$ 的初始值保持一样，${{\rm{\delta }}_{{\rm{1c}}}} = 0.999\;5$ ，${{\rm{\delta }}_{{\rm{2c}}}} = 1.000\;5, $ $ {\mu _1} = 0.9,\;M = 16$ 。将ADMM算法作为对比算法[18]，其中ADMM算法中惩罚因子$\rho $ 保持不变，其他参数迭代规则按照式(23)～式(30)所示。图2分别展示了IADPM与ADMM算法中的残差
$ \epsilon$ 随迭代次数变化曲线图。可以看到，对于IADPM算法，无论惩罚因子$\rho $ 的初始值取多少，残差$ \epsilon$ 总能不断保持下降趋势直至收敛，然而ADMM算法受惩罚因子$\rho $ 影响较大，并不能总是收敛，当$\rho $ 为1和100时，ADMM算法最后趋于收敛，然而当$\rho $ 取值较小为0.1时，ADMM算法无法收敛。
图 2
$ \epsilon$ 随迭代次数变化的曲线图 -
本小节针对设计的脉间幅度−相位波形序列的抗干扰效果进行仿真验证。首先考虑单目标与单干扰场景，参数设置与4.1节保持一致，设置噪声功率
${\sigma ^2} = 0$ dB，目标信号幅度${\left| {{\alpha _T}} \right|^2} = 10$ dB，干扰信号幅度${\left| {{\alpha _J}} \right|^2} = 20$ dB。图3为单目标单干扰场景下的抗干扰结果。可以看到，由于干扰的存在，利用传统脉间波形得到的慢时域多普勒谱上，目标和干扰同时存在，且干扰幅度比目标幅度高约10 dB。相反，利用优化后脉间波形得到的多普勒谱上，可以看到干扰尖峰消失，仅剩下目标尖峰，且在
$\left[ {0.3,0.5} \right]$ 阻带范围内存在凹槽。
图 3 单目标场景下抗干扰效果
考虑多目标与多干扰场景，设干扰个数
$P = 2$ ，目标个数$Q = 2$ ，干扰信号的归一化多普勒频率分别为${f_{J1}} = 0.35$ ，${f_{J2}} = 0.45 $ ，目标的归一化多普勒频率分别为${f_{T1}} = 0.42,\;{f_{T2}} = 0.44$ ，设置阻带范围为$\varOmega = \left[ {{f_{11}},{f_{12}}} \right] = \left[ {0.3,0.5} \right]$ ，其余参数保持不变。图4给出了多目标场景下抗干扰效果图。可以看到，利用传统脉冲序列得到的慢时域多普勒谱上，两个目标和两个干扰同时存在，无法分辨真实目标。通过优化设计的脉间波形的多普勒谱上，仅剩下两个真实目标，两个干扰尖峰消失，说明了该方法在多目标多干扰场景下依然适用。
图 4 多目标场景下抗干扰效果
-
针对速度欺骗干扰，本文通过建立最小化阻带内干扰能量与阻带内目标旁瓣能量的加权和的优化准则，考虑波形离散量化相位和PAR约束，构建了速度欺骗干扰下脉间波形幅相联合设计的优化问题，针对该优化问题，提出了IADPM算法求解，并理论分析了IADPM算法的计算复杂度。
仿真结果表明，IADPM算法在不同的惩罚因子下总能收敛，优于ADMM算法。离散相位情况下，设计的波形更加易于工程实现，同时可以有效抑制速度欺骗干扰，在多目标多干扰情况下，算法依然有效，普适性较高。
Optimized Amplitude-Phase Waveform Against Deceptive Jamming
-
摘要: 欺骗干扰作为一种典型的有源干扰，在距离维或多普勒维上产生虚假目标，使雷达发生虚警，严重影响雷达正常工作。针对速度欺骗干扰，该文提出了一种基于脉间波形幅相联合设计的抗速度欺骗干扰方法。首先，根据速度欺骗干扰下的雷达回波模型，建立最小化阻带内干扰能量与目标旁瓣能量的加权和的优化准则，同时考虑离散量化相位和峰均功率比(PAR)约束，形成脉间幅相联合设计抗速度欺骗干扰的优化问题；然后，提出非精确交替方向惩罚法(IADPM)算法求解该优化问题，并分析算法的计算复杂度；最后，仿真实验分析验证了该算法的有效性。Abstract: Deceptive jamming as a typical active jamming, produces false target on the range dimension or Doppler dimension, which increases probability of false alarm and significantly degrade the performance of radar. This paper is focused on the joint design of pulse amplitude and phase against velocity deceptive jamming. In particular, a design criterion incorporating a weighted sum of jamming energy and target sidelobe energy in the stopband is developed to minimize under phase quantized and peak-to-average power ratio (PAR) restriction. To tackle the resulting non-convex optimization problem, a new inexact alternating direction penalty method (IADPM) is proposed based on the ADPM framework and the computational complexity is analyzed. Finally, numerical results are provided to demonstrate the effectiveness of the proposed methodology.
-
Key words:
- jamming suppression /
- PAR constraint /
- velocity deceptive jamming /
- waveform design
-
[1] LI N, ZHANG Y. A survey of radar ECM and ECCM[J]. IEEE Transactions on Aerospace and Electronic Systems, 1995, 31(3): 1110-1120. doi: 10.1109/7.395232 [2] 赵国庆. 雷达对抗原理[M]. 第二版. 西安: 西安电子科技大学出版社, 2012. ZHAO Guo-qing. Principle of radar countermeasure[M]. 2nd ed. Xi’an: Xidian University Press, 2012. [3] 黄岩, 赵博, 陶明亮, 等. 合成孔径雷达抗干扰技术综述[J]. 雷达学报, 2020, 9(1): 86-106. HUANG Yan, ZHAO Bo, TAO Ming-liang, et al. Review of synthetic aperture radar interference suppression[J]. Journal of Radars, 2020, 9(1): 86-106. [4] 陈赓, 田波, 宫健, 等. 雷达有源干扰鉴别技术综述[J]. 现代防御技术, 2019, 47(5): 113-119. CHEN Geng, TIAN Bo, GONG Jian, et al. Overview of radar active jamming identification techniques[J]. Modern Defence Technology, 2019, 47(5): 113-119. [5] 马佳智, 施龙飞, 徐振海, 等. 单脉冲雷达多点源参数估计与抗干扰技术进展[J]. 雷达学报, 2019, 8(1): 125-139. MA Jia-zhi, SHI Long-fei, XU Zhen-hai, et al. Overview of multi-source parameter estimation and jamming mitigation for monopulse radars[J]. Journal of Radars, 2019, 8(1): 125-139. [6] AKHTAR J. Orthogonal block coded ECCM schemes against repeat radar jammers[J]. IEEE Transactions on Aerospace and Electronic Systems, 2009, 45(3): 1218-1226. doi: 10.1109/TAES.2009.5259195 [7] 罗美方, 郝志梅, 王强. LFM-频率编码复合低截获波形信号处理方法[J]. 电子科技大学学报, 2017, 46(2): 335-339. doi: 10.3969/j.issn.1001-0548.2017.02.003 LUO Mei-fang, HAO Zhi-mei, WANG Qiang. Signal processing technique for LFM/FSK hybrid LPI waveform[J]. Journal of University of Electronic Science and Technology of China, 2017, 46(2): 335-339. doi: 10.3969/j.issn.1001-0548.2017.02.003 [8] LU Shu-ping, CUI Guo-long, YI Wei, et al. Radar waveform design against signal-dependent jamming[C]// 2017 IEEE Radar Conference. Seattle, USA: IEEE, 2017: 1075-1080. [9] ZHOU Kai, LI De-xin, SU Yi, et al. Joint design of transmit waveform and mismatch filter in the presence of interrupted sampling repeater jamming[J]. IEEE Signal Processing Letters, 2020, 27: 1610-1614. doi: 10.1109/LSP.2020.3021667 [10] ZHANG Jin-dong, ZHU Dai-yin, ZHANG Gong. New anti-velocity deception jamming technique using pulses with adaptive initial phases[J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(2): 1290-1330. doi: 10.1109/TAES.2013.6494414 [11] YANG Ya, WU Jian, CUI Guo-long, et al. Optimized phase-coded waveform design against velocity deception[C]//2015 IEEE Radar Conference. Arlington, USA: IEEE, 2015: 1-4. [12] XIONG Wei, WANG Xin-hai, ZHANG Gong. Cognitive waveform design for anti-velocity deception jamming with adaptive initial phases[C]//2016 IEEE Radar Conference. Philadelphia, USA: IEEE, 2016: 1-5. [13] 吴健. 基于波形分集的雷达抗有源欺骗干扰技术研究[D]. 成都: 电子科技大学, 2015. WU Jian. Research of technology against radar active deception jamming based on waveform diversity[D]. Chengdu: University of Electronic Science and Technology of China, 2015. [14] 崔国龙, 余显祥, 杨婧, 等. 认知雷达波形优化设计综述[J]. 雷达学报, 2019, 8(5): 537-557. doi: 10.12000/JR19072 CUI Guo-long, YU Xian-xiang, YANG Jing, et al. An overview of waveform optimization methods for cognitive radar[J]. Journal of Radars, 2019, 8(5): 537-557. doi: 10.12000/JR19072 [15] YU Xian-xiang, CUI Guo-long, YANG Jing, et al. Quadratic optimization for unimodular sequence design via an ADPM framework[J]. IEEE Transactions on Signal Processing, 2020, 68: 3619-3634. doi: 10.1109/TSP.2020.2998637 [16] SOLTANALIAN M, STOICA P. Designing unimodular codes via quadratic optimization[J]. IEEE Transactions on Signal Processing, 2014, 62(5): 1221-1234. doi: 10.1109/TSP.2013.2296883 [17] YU Xian-xiang, CUI Guo-long, KONG Ling-jiang, et al. Constrained waveform design for colocated MIMO radar with uncertain steering matrices[J]. IEEE Transactions on Aerospace and Electronic Systems, 2019, 55(1): 356-370. doi: 10.1109/TAES.2018.2852200 [18] BOYD S, PARIKH N, CHU E. Distributed optimization and statistical learning via the alternating direction method of multipliers[M]. Hanover: Now Publishers Inc, 2011. -
点击查看大图
图(4)
计量
- 文章访问数: 4377
- HTML全文浏览量: 1623
- PDF下载量: 73
- 被引次数: 0
