间隙波导技术及其空间应用

陈翔; 孙冬全; 崔万照

doi:10.11999/JEIT211291

间隙波导技术及其空间应用

doi: 10.11999/JEIT211291

1.
中国空间技术研究院西安分院空间微波技术国家级重点实验室西安 710100
2.
西安电子科技大学物理学院西安 710071

基金项目: 国家自然科学基金(61901359, 61901320)

详细信息

作者简介:
陈翔：男，高级工程师，博士，研究方向为航天器微波毫米波技术

孙冬全：男，副教授，博士，研究方向为微波毫米波电路、天线理论及技术

崔万照：男，研究员，博士，研究方向为航天器微波技术

通讯作者:
孙冬全　dqsun87@163.com

中图分类号: TN81
计量
- 文章访问数: 3436
- HTML全文浏览量: 2401
- PDF下载量: 483
- 被引次数: 2
出版历程
- 收稿日期: 2021-11-18
- 修回日期: 2022-04-04
- 网络出版日期: 2022-04-12
- 刊出日期: 2023-01-17

Gap Waveguide Technology and its Space Applications

1.
National Key Laboratory of Science and Technology on Space Microwave, China Academy of Space Technology(Xi’an), Xi’an 710100, China
2.
School of Physics, Xidian University, Xi’an 710071, China

Funds: The National Natural Science Foundation of China (61901359, 61901320)

摘要

摘要: 间隙波导(GW)是一种基于非接触电磁带隙(EBG)结构的新型人工电磁(EM)材料，其独特的非接触结构和宽带电磁屏蔽特性在构建新型电磁传输线及屏蔽结构方面显示出极大的优势和灵活性，为微波毫米波部件、电路及天线等领域带来了新的研究视角和实现途径，近年来引起了广泛关注。该文首先简要介绍了间隙波导概念和原理，分析了其技术优势；进一步，根据不同的研究及应用领域分类，全方位地归纳总结了间隙波导技术相关的国内外研究进展情况；最后，结合空间技术背景和发展需求，探讨了间隙波导在空间微波毫米波技术中的应用前景，提出了基于间隙波导技术的非接触式无源互调干扰控制方法及堆叠集成毫米波电路系统两个重要的应用方向。该文工作可为间隙波导技术相关研究和应用提供一定的借鉴与参考。
- 微波毫米波 /
- 间隙波导 /
- 空间技术 /
- 无源互调 /
- 堆叠集成
Abstract: Gap Waveguide (GW) is a new kind of artificial ElectroMagnetic (EM) material based on contactless Electromagnetic Band Gap(EBG) structure. The unique contactless structure and wide EM forbidden band of GW show great advantages and flexibility in developing new EM transmission lines and shielding structures, providing a new research perspective and realization approach for microwave & millimeter-wave components, circuits and antennas, etc., and causing much attentions in recent years. Firstly, the concept and principle of GW are briefly introduced, and its technical advantages are analyzed. Then, the research progress of GW is comprehensively summarized, according to the classification of various research and application fields. In the end, the application prospect of GW in space microwave & millimeter-wave technology is discussed combining with space technology background and development needs, and the contactless suppression method of passive intermodulation and stacked-integrated technology of millimeter-wave systems are proposed. This work can provide some valuable references for the research and application of GW.
- Microwave & millimeter-wave /
- Gap Waveguide(GW) /
- Space technology /
- Passive intermodulation /
- Stacked-integrated

HTML全文

1. 引言

近年来，随着电子技术的不断发展，雷达所处的电磁工作环境正呈现复杂化和激烈化态势^[1]。尤其是数字射频存储器(Digital Radio Frequency Memory, DRFM)的出现和日臻完善，给雷达抗干扰带来了极大挑战^[2-4]。间歇采样转发干扰(Interrupted Sampling Repeater Jamming, ISRJ)作为DRFM衍生出的新型相参干扰样式，通过巧妙地对截获雷达信号进行分段转发，在雷达接收端脉压后形成与目标回波相似的密集假目标串^[5]，达到对干扰方实施压制和欺骗的双重效果。由于ISRJ具有响应速度快、工程实现简单的特点，传统方法往往难以有效抑制^[6,7]。因此，针对ISRJ的对抗技术成为当前抗干扰领域研究的热点问题。

目前，针对ISRJ的抑制方法可分为3个类别：第1类采用发射端波形设计的方法，具体做法是在发射端设计稀疏多普勒特性波形和正交LFM-相位编码波形，然后在接收端进行滑窗抽取检测和分段滤波来实现对ISRJ的识别和抑制^[8,9]。第2类采用接收端的信号处理方法，即基于干扰信号不同于目标回波的特征，在雷达接收端实现对ISRJ的滤除或抑制^[10-12]，该类方法大多依赖于对ISRJ的相关参数的估计^[13]；此外，也有学者通过深度学习对回波信号进行处理^[14,15]，在显著抑制干扰的同时摆脱了对参数估计的依赖。第3类方法综合采用发射端和接收端的联合设计^[6,7,16-19]，以提高抗干扰的自由度，简化接收端处理的复杂程度，代表性工作包括：文献[6]提出了一种抗ISRJ的联合设计算法，通过对非凸优化问题的求解实现了对ISRJ的显著抑制；文献[7]将模型扩展至多脉冲情形，并有效地解决了干扰机高速运动时的波形失配问题。

但需要指出的是，以文献[6]和文献[7]为代表的抗ISRJ联合设计方法中，性能准则“信干噪比”中的“信”并未考虑目标特性对发射信号的调制作用。但在宽带雷达条件下，目标常常不能再看作一个点目标，而是在距离维上具有多个散射中心的扩展目标^[19]。另外，值得注意的是，作为表征电磁波矢量性的参数，极化信息的利用可以显著提高雷达系统的目标检测^[19]、识别^[20]和抗干扰能力^[21,22]。在抗干扰方面，文献[21]在空域-极化域显著地抑制了多干扰源欺骗干扰；文献[22]利用正交极化高分辨距离像间的互相关系数作为鉴别量，达到了大于90%的正确鉴别概率。这充分说明了极化信息的引入对有效解决抗干扰问题具有重要价值。

针对现有技术的不足和短板，本文将在抗干扰设计中引入极化信息，并修正波形优化性能指标中信干噪比的数学形式，进而深入研究宽带全极化雷达的抗ISRJ发射波形和接收滤波器联合设计问题。实验结果表明，引入极化信息和考虑扩展目标对信号的调制作用显著提高了ISRJ的抑制性能。

2. 问题模型

假设雷达在一个相干处理间隔(Coherent Processing Interval, CPI)内发射 $K$ 个相位编码脉冲信号，发射波形序列的表达式为 ${\boldsymbol{s}} = [{\boldsymbol{s}}_1^{\text{T}},{\boldsymbol{s}}_2^{\text{T}},\cdots, {\boldsymbol{s}}_K^{\text{T}}]^{\text{T}} \in {\mathbb{C}^{NK \times 1}}$ ，其中 ${( \cdot )^{\text{T}}}$ 表示转置运算。不失一般性，序列中的第 $k$ 个脉冲表示为

${{\boldsymbol{s}}_k} = {[{s_k}(1),{s_k}(2), \cdots ,{s_k}(N)]^{\text{T}}} \in {\mathbb{C}^{N \times 1}},k = 1,2, \cdots ,K$

(1)

其中， $N$ 表示相位编码信号的码长。并且 $K$ 个相位编码脉冲信号满足恒模条件。

宽带雷达条件下，第 $k$ 个脉冲的目标冲激响应矩阵(Target Response Impulse Matrix, TIRM)表示为^[19]

$\begin{split} & {{\boldsymbol{H}}_k}(\theta ) = \left( {\begin{array}{*{20}{c}} {{H_0}(\theta )}&0& \cdots &0 \\ {{H_1}(\theta )}&{{H_0}(\theta )}& \cdots &0 \\ \vdots & \vdots & \ddots &0 \\ {{H_{Q - 1}}(\theta )}&{{H_{Q - 2}}(\theta )}& \cdots &0 \\ 0&{{H_{Q - 1}}(\theta )}& \cdots &0 \\ \vdots & \vdots & \ddots & \vdots \\ 0&0& \cdots &{{H_{Q - 1}}(\theta )} \end{array}} \right),\\ & Q = M - N + 1 \\[-10pt] \end{split}$

(2)

其中， $\theta$ 表示对应的目标方位角(Target Aspect Angle, TAA)。

假定目标在一个CPI内保持特性不变，则有 ${{\boldsymbol{H}}_1}(\theta ) = {{\boldsymbol{H}}_2}(\theta ) = \cdots = {{\boldsymbol{H}}_K}(\theta )$ ，则 $K$ 个目标回波可以表征为

${{\boldsymbol{s}}_{\text{T}}} = {\boldsymbol{H}}{\text{(}}\theta {\text{)}}{\boldsymbol{s}}$

(3)

其中， ${\boldsymbol{H}}{\text{(}}\theta {\text{)}} = {\text{diag(}}{{\boldsymbol{H}}_1}{\text{(}}\theta {\text{)}},\;{{\boldsymbol{H}}_2}{\text{(}}\theta {\text{)}}, \;\cdots ,\;{{\boldsymbol{H}}_K}{\text{(}}\theta {\text{))}} \in {\mathbb{C}^{MK \times NK}}$ 。

需要指出的是， ${\boldsymbol{H}}(\theta )$ 对方位角扰动十分敏感，故本文将给定方位角的先验范围 $[{\theta _{{\text{min}}}},{\theta _{{\text{max}}}}]$ ，并取该方位角范围内TIRM的均值作为目标特性，以此来增加系统设计的鲁棒性^[19]，即

$\overline {{\boldsymbol{H}}(\theta )} = \frac{1}{{{K_{\text{H}}}}}\sum\limits_{i = 1}^{{K_{\text{H}}}} {{\boldsymbol{H}}({\theta _i})}$

(4)

其中， ${K_{\text{H}}}$ 为对方位角区间 $[{\theta _{{\text{min}}}},{\theta _{{\text{max}}}}]$ 离散化的份数。

令 ${\boldsymbol{\varepsilon}} = {[{{\boldsymbol{\varepsilon}} _{\text{H}}},{{\boldsymbol{\varepsilon}} _{\text{V}}}]^{\text{T}}}$ 为全极化雷达发射天线的Jones矢量，则 ${\boldsymbol{\varepsilon}}$ 可以写为

${\boldsymbol{\varepsilon}} = \left[ {\begin{array}{*{20}{c}} {\cos \alpha }&{\sin \alpha } \\ { - \sin \alpha }&{\cos \alpha } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\cos \beta } \\ {{\text{j}}\sin \beta } \end{array}} \right]$

(5)

其中， $\alpha \in [ - {\pi}/2,{\pi}/2]$ 为极化旋转角， $\beta \in [ - {\pi}/4,{\pi}/4]$ 为极化椭圆角。假设雷达发射和接收信号共用一副天线，且一个CPI内不改变天线的极化方式，则有接收天线极化 ${\boldsymbol{\eta}} = {\boldsymbol{\varepsilon}}$ 。则全极化雷达接收到的回波序列表示为

${{\boldsymbol{x}}_{\text{T}}} = ({{\boldsymbol{\eta}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} )\overline {{\boldsymbol{H}}{\text{(}}\theta {\text{)}}} {\boldsymbol{s}} = ({{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T}}{\boldsymbol{\varepsilon}} )\overline {{\boldsymbol{H}}{\text{(}}\theta {\text{)}}} {\boldsymbol{s}} \in {\mathbb{C}^{MK \times 1}}$

(6)

其中， ${\boldsymbol{T}}$ 为目标的极化散射矩阵。

另外，ISRJ可用矩阵 ${\boldsymbol{J}} = {\text{Diag(}}{\boldsymbol{j}}{\text{)}} \in {\mathbb{C}^{N \times N}}$ 表示^[7]，则对应于第 $k$ 个发射脉冲的ISRJ信号可表示为

${{\boldsymbol{s}}_{{\text{J}}k}} = ({{\boldsymbol{W}}_{\text{J}}}{{\boldsymbol{J}}_k}){{\boldsymbol{s}}_k} = {{\boldsymbol{J}}_M}{{\boldsymbol{s}}_k}$

(7)

其中， ${{\boldsymbol{W}}}_{\text{J}}\in {ℂ}^{M\times N}$ 且有

${{\boldsymbol{W}}_{\rm{J}}}\left[ {p,q} \right] = \left\{ {\begin{array}{*{20}{c}} {1{\text{ , }}p = q} \\ {0{\text{ , }}p \ne q} \end{array}} \right.$

(8)

其中， $p = 1,2, \cdots ,M$ , $q = 1,2, \cdots ,N$ ，等价于在干扰矢量 ${\boldsymbol{j}}$ 后补 $M - N$ 个 $0$ 。同目标特性相似，假定干扰特性在一个CPI内也保持不变，即有 ${{\boldsymbol{J}}_1} = {{\boldsymbol{J}}_2} = \cdots ={{\boldsymbol{J}}_K} = {\boldsymbol{J}}$ ，则雷达在一个CPI内接收到的 $K$ 个ISRJ信号可综合表示为

${{\boldsymbol{s}}_{\text{J}}} = {[{\boldsymbol{s}}_{{\text{J}}1}^{\text{T}},{\boldsymbol{s}}_{{\text{J}}2}^{\text{T}},\cdots,{\boldsymbol{s}}_{{\text{J}}K}^{\text{T}}]^{\text{T}}} = {\boldsymbol{C}}\cdot {\boldsymbol{s}}$

(9)

其中， ${\boldsymbol{C}}$ 是 $K$ 个ISRJ信号对应的干扰特性矩阵，有 ${\boldsymbol{C}} = {\text{diag}}({{\boldsymbol{W}}_{\text{J}}}J,{{\boldsymbol{W}}_{\text{J}}}{\boldsymbol{J}}, \cdots ,{{\boldsymbol{W}}_{\text{J}}}{\boldsymbol{J}}) \in {\mathbb{C}^{MK \times NK}}$ 。可以看出， ${\boldsymbol{C}}$ 是由K个发射脉冲干扰特性矩阵拼接而成的分块对角矩阵。则全极化雷达接收到的干扰回波序列为

${{\boldsymbol{x}}_{\text{J}}} = ({{\boldsymbol{\eta}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}){{\boldsymbol{s}}_{\text{J}}} = ({{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}){{\boldsymbol{s}}_{\text{J}}} \in {\mathbb{C}^{MK \times 1}}$

(10)

其中， ${{\boldsymbol{\varepsilon}} _{\text{J}}}$ 为干扰机天线极化方式的Jones矢量。

设接收滤波器为 ${\boldsymbol{w}} = {[{\boldsymbol{w}}_1^{\text{T}},{\boldsymbol{w}}_2^{\text{T}},\cdots,{\boldsymbol{w}}_K^{\text{T}}]^{\text{T}}} \in {\mathbb{C}^{MK \times 1}}$ ，则不失一般性，雷达回波经接收滤波器处理后的输出结果为

${\boldsymbol{r}} = ({{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T}}{\boldsymbol{\varepsilon}} )\cdot {{\boldsymbol{w}}^{\text{H}}}\overline {{\boldsymbol{H}}{\text{(}}\theta {\text{)}}} {\boldsymbol{s}} + ({{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}})\cdot {{\boldsymbol{w}}^{\text{H}}}{{\boldsymbol{s}}_{\text{J}}} + {{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{v}}$

(11)

其中， ${\boldsymbol{v}} \in {\mathbb{C}^{MK \times 1}}$ 为复加性高斯白噪声序列。

设 $r_{{x_{\text{T}}}h}^k(n)$ 为第 $k$ 个目标回波 ${{\boldsymbol{x}}_{{\text{T}}k}}$ 和相应接收滤波器 ${{\boldsymbol{h}}_k}$ 的互相关函数，即

$\begin{split} & r_{{{\boldsymbol{x}}_{\text{T}}}h}^k(n) = \sum\limits_{m = 1}^M ({{\boldsymbol{x}}_{{\text{TH}}k}}(m)h_{{\text{H}}k}^*(m + n) \\ & \quad + {{\boldsymbol{x}}_{{\text{TV}}k}}(m)h_{{\text{V}}k}^*(m + n)) ,{\text{ }}\\ & n = 1 - M,2 - M,\cdots,0,\cdots,M - 2,M - 1 \end{split}$

(12)

则目标回波与接收滤波器的互模糊函数为

${\chi _{{{\boldsymbol{x}}_{\text{T}}}h}}(n,f) = \sum\limits_{k = 1}^K {r_{{{\boldsymbol{x}}_{\text{T}}}h}^k(n)\exp ({\text{j}}kf)}$

(13)

其中， $f = 2{\pi}{f_{\text{d}}}{T_{\text{p}}}$ , ${f_{\text{d}}}$ 为多普勒频移， ${T_{\text{p}}}$ 为雷达的脉冲重复周期。干扰信号与接收滤波器的互模糊函数为

${\chi _{{{\boldsymbol{x}}_{\text{J}}}h}}(n,f) = \sum\limits_{k = 1}^K {r_{{{\boldsymbol{x}}_{\text{J}}}h}^k(n)\exp ({\text{j}}kf)}$

(14)

联合设计的目的，一方面在于在干扰机高速运动的情况下，目标回波能够最大限度地匹配目标特性，而在一定的多普勒频移区间内保持较高的脉压增益。另一方面，尽可能地让多普勒频移区间内干扰信号与接收滤波器之间的互模糊函数能量保持为零，从而对ISRJ的形成抑制。综合以上要求，对多普勒频移区间(f₁, f₂)进行离散化处理后，构造如式(15)的目标函数分量：

$\begin{split} \xi = & \sum\nolimits^{L}_{l=1} \left( \sum\nolimits_{n = 1 - N,n \ne 0}^{N - 1} \left| \sum\nolimits_{k = 1}^K (r_{{{\boldsymbol{x}}_{{\text{TH}}}}{{\boldsymbol{h}}_{\text{H}}}}^k(n) \right.\right.\\ & + r_{{{\boldsymbol{x}}_{{\text{TV}}}}{{\boldsymbol{h}}_{\text{V}}}}^k(n))\exp ({\text{j}}kf) |^2 + \sum\nolimits_{n = 1 - M}^{M - 1} \\ & \cdot\left. {{{\left| {\sum\limits_{k = 1}^K {(r_{{{\boldsymbol{x}}_{{\text{JH}}}}{{\boldsymbol{h}}_{\text{H}}}}^k(n) + r_{{{\boldsymbol{x}}_{{\text{JV}}}}{{\boldsymbol{h}}_{\text{V}}}}^k(n))\exp ({\text{j}}kf)} } \right|}^2}}\right) \end{split}$

(15)

其中， ${f_l} = {f_1} + (l - 1){{\Delta }}f,l = 1,2,\cdots,L$ ， $L$ 为离散化份数， ${{\Delta }}f = ({f_2} - {f_1})/(L - 1)$ 为等份划分后的频率间隔。

另外，由于非匹配滤波体制的使用，在对目标回波进行接收滤波处理时，脉压峰值将存在峰值损耗，即存在信噪比损失(Signal-to-Noise Ratio Loss, SNRL)。由矩阵的相容性，有

$\left\| {({{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} )\cdot{{\boldsymbol{w}}^{\text{H}}}\overline {{\boldsymbol{H}}(\theta )} {\boldsymbol{s}}} \right\| \le \left\| {{{\boldsymbol{\varepsilon}} ^2}} \right\|\left\| T \right\|\left\| {\boldsymbol{w}} \right\|\left\| {\overline {{\boldsymbol{H}}(\theta )} } \right\|\left\| {\boldsymbol{s}} \right\|$

(16)

其中， $\left\| \cdot \right\|$ 表示向量或矩阵的 $2$ -范数。由此可定义SNRL的表达式为

${\text{SNRL}} = 20\lg \frac{{\left\| {\boldsymbol{T}} \right\|\left\| {\boldsymbol{w}} \right\|\left\| {\overline {{\boldsymbol{H}}(\theta )} } \right\|\left\| {\boldsymbol{s}} \right\|}}{{\left| {({{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} )\cdot {{\boldsymbol{w}}^{\text{H}}}\overline {{\boldsymbol{H}}(\theta )} {\boldsymbol{s}}} \right|}}$

(17)

故可给目标回波施加如式(18)的峰值约束条件，保证脉压峰值具备一定的稳定性，即

$\begin{split} {g_1}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) =& {\left| {({\boldsymbol{h}}_{\text{H}}^{\text{H}}{{\boldsymbol{x}}_{{\text{TH}}}} + {\boldsymbol{h}}_{\text{V}}^{\text{H}}{{\boldsymbol{x}}_{{\text{TV}}}}) - {a_{{\text{max}}}}} \right|^2} \\ & = {\left| {({{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} )\cdot {{\boldsymbol{w}}^{\text{H}}}\overline {{\boldsymbol{H}}(\theta )} {\boldsymbol{s}} - {a_{{\text{max}}}}} \right|^2} \end{split}$

(18)

代入恒模约束条件 ${\boldsymbol{s}} = NK$ ，以及 ${{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{w}} = {M_{\boldsymbol{w}}}$ ，其中 ${M_{\boldsymbol{w}}}$ 为接收滤波器能量，则有

${a_{{\text{max}}}} = \sqrt {NK{M_{\boldsymbol{w}}}} \left\| {\boldsymbol{T}} \right\|\left\| {\overline {{\boldsymbol{H}}(\theta )} } \right\|{10^{ - \frac{\mu }{{20}}}}$

(19)

然后，为进一步抑制干扰信号，同样对干扰施加峰值约束，以将干扰信号的脉压峰值尽可能降低，即

$\begin{split} {g_2}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) & = {\left| {({\boldsymbol{h}}_{\text{H}}^{\text{H}}{{\boldsymbol{x}}_{{\text{JH}}}} + {\boldsymbol{h}}_{\text{V}}^{\text{H}}{{\boldsymbol{x}}_{{\text{JV}}}}) - {a_{{\text{min}}}}} \right|^2} \\ & = {\left| {({{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}})\cdot {{\boldsymbol{w}}^{\text{H}}}{{\boldsymbol{s}}_{\text{J}}} - {a_{{\text{min}}}}} \right|^2} \end{split}$

(20)

综上，能量约束下天线极化方式、发射波形序列和接收滤波器序列的联合优化问题可以表示为

$\left. \begin{split} & \mathop {{\text{min }}}\limits_{{\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} } \varGamma ({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) = \xi ({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) + {\gamma _1}{g_1}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) \\ & \quad + {\gamma _2}{g_2}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} )\\ & \text{s}\text{.t}\text{.}\;\;{{\boldsymbol{w}}}^{\text{H}}{\boldsymbol{w}}={M}_{w}\text{,}\left|{s}_{k}(n)\right|=1,k=1,2,\cdots, K,\\ & \quad n=1,2,\cdots ,N \end{split}\right\}$

(21)

其中， ${\gamma _1}$ 和 ${\gamma _2}$ 表示对应的峰值约束权值，需要指出的是，为保证惩罚力度， ${\gamma _1}$ 和 ${\gamma _2}$ 通常取较大数值^[23]，这样才能有效地滤波器的输出，同时抑制干扰信号。至此，本节完成了宽带全极化雷达抗ISRJ问题的数学建模。

3. 问题求解

对于式 $\left(21\right)$ ，本文将采用交替迭代的方法进行分解，具体做法为：在给定发射波形序列 ${{\boldsymbol{s}}^{(0)}}$ 和Jones矢量 ${{\boldsymbol{\varepsilon}} ^{(0)}}$ 初始值的基础上，先固定 ${{\boldsymbol{s}}^{(i)}}$ 和 ${{\boldsymbol{\varepsilon}} ^{(i)}}$ ，可求解 ${{\boldsymbol{w}}^{(i + 1)}}$ ；再固定 ${{\boldsymbol{w}}^{(i + 1)}}$ 和 ${{\boldsymbol{\varepsilon}} ^{(i)}}$ ，可求解 ${{\boldsymbol{s}}^{(i + 1)}}$ ；最后固定 ${{\boldsymbol{s}}^{(i + 1)}}$ 和 ${{\boldsymbol{w}}^{(i + 1)}}$ ，求解 ${{\boldsymbol{\varepsilon}} ^{(i + 1)}}$ ，重复上述流程，直至达到设定的收敛条件，上述迭代过程的数学形式依次可表示为

${{\boldsymbol{w}}^{(i + 1)}} = {\text{arg min }}\varGamma \left( {{{\boldsymbol{s}}^{(i)}},{\boldsymbol{w}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right)$

(22)

${{\boldsymbol{s}}^{(i + 1)}} = {\text{arg min }}\varGamma \left( {{\boldsymbol{s}},{{\boldsymbol{w}}^{(i + 1)}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right)$

(23)

${{\boldsymbol{\varepsilon}} ^{(i + 1)}} = {\text{arg min }}\varGamma \left( {{{\boldsymbol{s}}^{(i + 1)}},{{\boldsymbol{w}}^{(i + 1)}},{\boldsymbol{\varepsilon}} } \right)$

(24)

其中， ${{\boldsymbol{s}}^{(i)}}$ , ${{\boldsymbol{w}}^{(i)}}$ 和 ${{\boldsymbol{\varepsilon }}^{(i)}}$ 分别表示第 $i$ 次迭代后发射波形序列 ${\boldsymbol{s}}$ 、接收滤波器序列 ${\boldsymbol{w}}$ 和Jones矢量 ${\boldsymbol{\varepsilon}}$ 的最优解。收敛条件方面，可设置 $\left| {{\varGamma ^{(i + 1)}}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} ) - {\varGamma ^{(i)}}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon}} )} \right|/ \left| {{\varGamma ^{(i)}}({\boldsymbol{s}},{\boldsymbol{w}},{\boldsymbol{\varepsilon }})} \right| < {\eta _\varGamma }$ 作为算法的收敛条件^[24,25]， ${\eta _\varGamma }$ 为预设值。下面将分别求解迭代过程中接收滤波器 ${{\boldsymbol{s}}^{(i + 1)}}$ 、发射波形 ${{\boldsymbol{w}}^{(i + 1)}}$ 和Jones矢量 ${{\boldsymbol{\varepsilon}} ^{(i + 1)}}$ 对应的子问题。

根据式 $\left(21\right)$ 和式 $\left(22\right)$ ，在迭代求解的第i步，固定发射波形 ${{\boldsymbol{s}}^{(i)}}$ 和Jones矢量 ${{\boldsymbol{\varepsilon}} ^{{\text{(}}i{\text{)}}}}$ ，求解 ${{\boldsymbol{w}}^{(i + 1)}}$ 的问题为

$\begin{split} & \mathop {{\text{min }}}\limits_{\boldsymbol{w}} \varGamma ({\boldsymbol{w}}) = \xi \left( {{{\boldsymbol{s}}^{(i)}},{\boldsymbol{w}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right) + {\gamma _1}{g_1}\left( {{{\boldsymbol{s}}^{(i)}},{\boldsymbol{w}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right) \\ & \quad + {\gamma _2}{g_2}\left( {{{\boldsymbol{s}}^{(i)}},{\boldsymbol{w}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right)\\ & {\text{s}}{\text{.t}}{\text{. }}{{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{w}} = {M_{\boldsymbol{w}}}\\[-10pt] \end{split}$

(25)

为了便于对目标函数进行优化，分别定义目标回波矩阵 ${{\boldsymbol{S}}_{{\text{T}}l}}$ 和ISRJ矩阵 ${{\boldsymbol{S}}_{{\text{J}}l}}$ 为

$\begin{split} & {{\boldsymbol{S}}}_{\text{T}l}\left[p,q\right]\\ & =\left\{ \begin{aligned} & {s}_{\text{T}k+1}\left(M - n\right)\text{exp}\left({\rm{j}}\left(k + 1\right){f}_{l}\right),q - \left(p - Mk\right)=n\\ & 0,\;{\text{其他}} \end{aligned}\right. \end{split}$

(26)

$\begin{split} & {{\boldsymbol{S}}}_{\text{J}l}\left[p,q\right]\\ & =\left\{\begin{aligned} & {s}_{\text{J}k+1}\left(M - n\right)\text{exp}\left({\rm{j}}\left(k + 1\right){f}_{l}\right),q - \left(p - Mk\right)=n\\ & 0,\; 其他\end{aligned}\right. \end{split}$

(27)

其中， $l = 1,2, \cdots ,L$ , $n = 0,1, \cdots ,M - 1$ , $p = 1, 2, \cdots ,MK$ , $q = 1,2, \cdots ,2M - 1$ , $k = {\text{quo}}(p - 1,M)$ 表示 $\left\lfloor {\dfrac{{p - 1}}{M}} \right\rfloor$ , $\left\lfloor \cdot \right\rfloor$ 表示向下取整操作，于是，优化问题(25)可以等效写为

$\begin{split} & \mathop {{\text{min }}}\limits_{\boldsymbol{w}} \varGamma ({\boldsymbol{w}}) = {{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{Pw}} - 2{\text{Re}}\left( {{{\boldsymbol{p}}^{\text{H}}}{\boldsymbol{w}}} \right)\\ & {\text{s}}{\text{.t}}{\text{. }}\;{{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{w}} = {M_{\boldsymbol{w}}} \end{split}$

(28)

其中

$\begin{split} {\boldsymbol{P}} =& \sum\limits_{l = 1}^L {\left( {{{\left| {{\boldsymbol{{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} } \right|}^2}{{\boldsymbol{S}}_{{\text{T}}l}}{\boldsymbol{W}}{\boldsymbol{S}}_{{\text{T}}l}^{\text{H}} + {{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}} \right|}^2}{{\boldsymbol{S}}_{{\text{J}}l}}{{\boldsymbol{S}}_{{\text{J}}l}}^{\text{H}}} \right)} \\ & + {\gamma _1}{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} } \right|^2}{{\boldsymbol{s}}_{\text{T}}}{{\boldsymbol{s}}_{\text{T}}}^{\text{H}} + {\gamma _2}{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}} \right|^2}{{\boldsymbol{s}}_{\text{J}}}{{\boldsymbol{s}}_{\text{J}}}^{\text{H}}\\[-1pt] \end{split}$

(29)

${\boldsymbol{p}} = {\gamma _1}{a_{{\text{max}}}}{{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T\varepsilon}} \cdot {{\boldsymbol{s}}_{\text{T}}} + {\gamma _2}{a_{{\text{min}}}}{{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}\cdot {{\boldsymbol{s}}_{\text{J}}}$

(30)

${\text{Re(}} \cdot {\text{)}}$ 表示取实部操作， ${\boldsymbol{W}}$ 表示维度为 $2M - 1$ 的对角方阵(除了第 $M$ 个对角元素等于0之外，其余对角元素为1，且对角元素以外的元素均为0)。

根据文献[]，式 $\left(28\right)$ 中的最优解为

${{\boldsymbol{w}}^{(i + 1)}} = - \sqrt {\frac{{{M_{\boldsymbol{w}}}}}{{{{\left\|u\left( {{{\boldsymbol{w}}^{(i)}}} \right)\right\|}^2}}}} u\left( {{{\boldsymbol{w}}^{(i)}}} \right)$

(31)

其中

$u\left( {{{\boldsymbol{w}}^{(i)}}} \right) = \left( {{{\boldsymbol{P}}^{(i)}} - {\text{tr}}\left( {{{\boldsymbol{P}}^{(i)}}} \right){{\boldsymbol{I}}_{MK}}} \right){{\boldsymbol{w}}^{(i)}} - {{\boldsymbol{p}}^{(i)}}$

(32)

根据式 $\left(21\right)$ 和式 $\left(23\right)$ ，固定接收滤波器序列 ${{\boldsymbol{w}}^{(i + 1)}}$ 和Jones矢量 ${{\boldsymbol{\varepsilon}} ^{(i)}}$ ，求解 ${{\boldsymbol{s}}^{(i + 1)}}$ 的问题形式可以表示为

$\begin{split} & \mathop {{\text{min }}}\limits_s \varGamma ({\boldsymbol{s}}) = \xi \left( {{\boldsymbol{s}},{{\boldsymbol{w}}^{(i + 1)}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right) + {\gamma _1}{g_1}\left( {{\boldsymbol{s}},{{\boldsymbol{w}}^{(i + 1)}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right) \\ & \quad + {\gamma _2}{g_2}\left( {{\boldsymbol{s}},{{\boldsymbol{w}}^{(i + 1)}},{{\boldsymbol{\varepsilon}} ^{(i)}}} \right)\\ & {\text{s}}{\text{.t}}{\text{. }}\;\left| {{s_k}(n)} \right| = 1,k = 1,2,\cdots,K,n = 1,2, \cdots ,N \\[-10pt] \end{split}$

(33)

为求解式(33)，定义接收滤波器矩阵 ${{\boldsymbol{R}}_l}$

$\begin{split} & {{\boldsymbol{R}}}_{l}\left[p,q\right]\\ & =\left\{ \begin{aligned} & {w}_{k+1}\left(n\right)\text{exp}\left(j\left(k + 1\right){f}_{l}\right), q + \left(p - Mk\right) = M + n\\ & 0, {\text{其他}}\end{aligned}\right. \end{split}$

(34)

其中， $l = 1,2, \cdots ,L$ , $n = 1,2, \cdots ,M$ , $p = 1,2, \cdots ,MK$ , $q=1,2,\cdots ,2M-1,\;k=\text{quo}(p-1,M)$ 。

于是，式 $\left(33\right)$ 可以等效为

$\begin{split} & \mathop {{\text{min }}}\limits_{\boldsymbol{s}} \varGamma (s) = {s^{\text{H}}}{\boldsymbol{Qs}} - 2{\text{Re}}\left( {{{\boldsymbol{q}}^{\text{H}}}{\boldsymbol{s}}} \right)\\ & \text{ s}\text{.t}\text{.}\;\left|{s}_{k}(n)\right|=1,k=1,2,\cdots,K,n=1,2,\cdots ,N \end{split}$

(35)

其中

$\begin{split} {\boldsymbol{Q}} =& \sum\limits_{l = 1}^L \left( {{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T}}{\boldsymbol{\varepsilon}} } \right|}^2}{{\overline {{\boldsymbol{H}}(\theta )} }^{\text{H}}}{{\boldsymbol{R}}_l}{\boldsymbol{W}}{{\boldsymbol{R}}_l}^{\text{H}}\overline {{\boldsymbol{H}}(\theta )} \right. \\ & \left.+ {{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}} \right|}^2}{{\boldsymbol{C}}^{\text{H}}}{{\boldsymbol{R}}_l}{{\boldsymbol{R}}_l}^{\text{H}}{\boldsymbol{C}} \right) \\ & + {\gamma _1}{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{\boldsymbol{T}}{\boldsymbol{\varepsilon}} } \right|^2}{\overline {{\boldsymbol{H}}(\theta )} ^{\text{H}}}{\boldsymbol{w}}{{\boldsymbol{w}}^{\text{H}}}\overline {{\boldsymbol{H}}(\theta )} \\ & + {\gamma _2}{\left| {{{\boldsymbol{\varepsilon}} ^{\text{H}}}{{\boldsymbol{\varepsilon}} _{\text{J}}}} \right|^2}{{\boldsymbol{C}}^{\text{H}}}{\boldsymbol{w}}{{\boldsymbol{w}}^{\text{H}}}{\boldsymbol{C}} \end{split}$

(36)

${\boldsymbol{q}}={\gamma }_{1}{a}_{\text{max}}{({{\boldsymbol{\varepsilon}} }^{\text{H}}{\boldsymbol{T}}{\boldsymbol{\varepsilon}} )}^{\text{H}}{\overline{{\boldsymbol{H}}(\theta )}}^{\text{H}}{\boldsymbol{w}}+{\gamma }_{2}{a}_{\text{min}}{({{\boldsymbol{\varepsilon}} }^{\text{H}}{{\boldsymbol{\varepsilon}} }_{\text{J}})}^{\text{H}}{{\boldsymbol{C}}}^{\text{H}}{\boldsymbol{w}}$

(37)

对于问题式(35)，有最优解为^[7]

${{\boldsymbol{s}}^{(i + 1)}} = - {\text{exp}}\left( {{\text{j}} \cdot {\text{arg}}\left( {v\left( {{{\boldsymbol{s}}^{(i)}}} \right)} \right)} \right)$

(38)

其中

$v\left( {{{\boldsymbol{s}}^{(i)}}} \right) = \left( {{{\boldsymbol{Q}}^{(i)}} - {\text{tr}}\left( {{{\boldsymbol{Q}}^{(i)}}} \right){{\boldsymbol{I}}_{NK}}} \right){{\boldsymbol{s}}^{(i)}} - {{\boldsymbol{q}}^{(i)}}$

(39)

根据式 $\left(21\right)$ 和式 $\left(24\right)$ ，固定接收滤波器序列 ${{\boldsymbol{w}}^{(i + 1)}}$ 和发射波形序列 ${{\boldsymbol{s}}^{(i + 1)}}$ ，求解Jones矢量 ${{\boldsymbol{\varepsilon}} ^{(i + 1)}}$ 的问题为

$\begin{split} \mathop {{\text{min }}}\limits_{\boldsymbol{\varepsilon}} \varGamma {\text{(}}{\boldsymbol{\varepsilon}} {\text{)}} = & \xi \left( {{{\boldsymbol{s}}^{(i + 1)}},{{\boldsymbol{w}}^{(i + 1)}},{\boldsymbol{\varepsilon}} } \right) \\ & + {\gamma _1}{g_1}\left( {{{\boldsymbol{s}}^{(i + 1)}},{{\boldsymbol{w}}^{(i + 1)}},{\boldsymbol{\varepsilon}} } \right) \\ & + {\gamma _2}{g_2}\left( {{{\boldsymbol{s}}^{(i + 1)}},{{\boldsymbol{w}}^{(i + 1)}},{\boldsymbol{\varepsilon}} } \right) \\ \text{ s}\text{.t}\text{. }{{\boldsymbol{\varepsilon}} }^{\text{H}}{\boldsymbol{\varepsilon}} =1& \end{split}$

(40)

由于Jones矢量 ${\boldsymbol{\varepsilon}}$ 的未知量个数恒定为2，可采用网格法求解。具体地，将 $\alpha \in [ - {\pi}/2,{\pi}/2]$ 等分为 ${L_\alpha }$ 份，将 $\beta \in [ - {\pi}/4,{\pi}/4]$ 等分为 ${L_\beta }$ 份，算法遍历 ${L_\alpha }\cdot {L_\beta }$ 次后求得最优解 ${{\boldsymbol{\varepsilon}} ^{(i + 1)}}$ 。

综合以上发射波形序列 ${\boldsymbol{s}}$ 、接收滤波器序列 ${\boldsymbol{w}}$ 和Jones矢量 ${\boldsymbol{\varepsilon}}$ 的交替迭代求解方法，即可完成抗ISRJ的恒模互补波形最优化设计。为进一步缩短算法收敛所需的时长，可在不牺牲算法性能的前提下，引入平方迭代框架(SQUAREd iterative Method, SQUAREM)实现算法收敛速度的显著提升^[25]。使用平方迭代框架后，宽带全极化雷达的恒模互补波形设计算法流程如算法1所示。

算法1　宽带全极化雷达抗ISRJ恒模互补波形设计流程
输入：目标脉冲响应矩阵 $\overline {{\boldsymbol{H}}(\theta )}$ 、干扰特性矩阵 ${\boldsymbol{C}}$ 、初始 ${\boldsymbol{s}}$ , ${\boldsymbol{w}}$ 和　　　　 ${\boldsymbol{\varepsilon}}$
输出：最优发射波形序列 ${{\boldsymbol{s}}^\nabla }$ 、接收滤波器序列 ${{\boldsymbol{w}}^\nabla }$ 和Jones矢量 ${{\boldsymbol{\varepsilon}} ^\nabla }$
1 重复
2 将 ${{\boldsymbol{w}}^{(i)} }$ 代入式(31)计算 ${{\boldsymbol{w}}_1}$ ，将 ${{\boldsymbol{w}}_1}$ 代入式(31)计算 ${{\boldsymbol{w}}_2}$
3 　　 ${{\boldsymbol{r}}_1} = {{\boldsymbol{w}}_1} - {{\boldsymbol{w}}^{(i)} }$ , ${{\boldsymbol{v}}_1} = {{\boldsymbol{w}}_2} - {{\boldsymbol{w}}_1} - {{\boldsymbol{r}}_1}$ $,$ ${\alpha _1} = - {{\boldsymbol{r}}_1}/{{\boldsymbol{v}}_1}$
4 　　将 ${ {\boldsymbol{w} }^{(i)} } - 2{\alpha _1}{ {\boldsymbol{r} }_1} + {\alpha _1}^2{{\boldsymbol{v}}_1}$ 代入式(31)计算 ${{\boldsymbol{w}}^{(i + 1)}}$
5 　　若 $\varGamma \left( {{{\boldsymbol{s}}^{{\text{(}}i{\text{)}}}},{{\boldsymbol{w}}^{(i + 1)}},{{\boldsymbol{\varepsilon }}^{{\text{(}}i{\text{)}}}}} \right) > \varGamma \left( {{{\boldsymbol{s}}^{{\text{(}}i{\text{)}}}},{{\boldsymbol{w}}^{{\text{(}}i{\text{)}}}},{{\boldsymbol{\varepsilon}} ^{{\text{(}}i{\text{)}}}}} \right)$ ，循环
6 　　　　 ${\alpha _1} \leftarrow \left( {{\alpha _1} - 1} \right)/2$
7 　　　　将 ${{\boldsymbol{w}} ^{{\text{(}}i{\text{)}}}} - 2{\alpha _1}{{\boldsymbol{r}}_1} + {\alpha _1}^2{{\boldsymbol{v}}_1}$ 代入式(31)计算 ${{\boldsymbol{w}}^{{\text{(}}i + 1{\text{)}}}}$
8 　　结束
9 　　将 ${{\boldsymbol{s}}^{{\text{(}}i{\text{)}}}}$ 代入式(38)计算 ${{\boldsymbol{s}}_1}$ ，将 ${{\boldsymbol{s}}_1}$ 代入式(38)计算 ${{\boldsymbol{s}}_2}$
10 　 ${{\boldsymbol{r}}_2} = {{\boldsymbol{s}}_1} - {{\boldsymbol{s}}^{{\text{(}}i{\text{)}}}}$ , ${{\boldsymbol{v}}_2} = {{\boldsymbol{s}}_2} - {{\boldsymbol{s}}_1} - {{\boldsymbol{r}}_1},{\text{ }}{\alpha _2} = - {{\boldsymbol{r}}_2}/{{\boldsymbol{v}}_2}$
11 　将 ${{\boldsymbol{s}}^{ {\text{(} }i{\text{)} } } } - 2{\alpha _2}{{\boldsymbol{r}}_2} + {\alpha _2}^2{{\boldsymbol{v}}_2}$ 代入式(38)计算 ${s^{{\text{(}}i + 1{\text{)}}}}$
12 　若 $\varGamma \left( {{{\boldsymbol{s}}^{{\text{(}}i + 1{\text{)}}}},{{\boldsymbol{w}}^{{\text{(}}i + 1{\text{)}}}},{{\boldsymbol{\varepsilon}} ^{{\text{(}}i{\text{)}}}}} \right) > \varGamma \left( {{{\boldsymbol{s}}^{{\text{(}}i{\text{)}}}},{{\boldsymbol{w}}^{{\text{(}}i + 1{\text{)}}}},{{\boldsymbol{\varepsilon}} ^{{\text{(}}i{\text{)}}}}} \right)$ ，循环
13 　　　　 ${\alpha _2} \leftarrow \left( {{\alpha _2} - 1} \right)/2$
14 　　　　将 ${ {\boldsymbol{s} }^{ {\text{(} }i{\text{)} } } } - 2{\alpha _2}{{\boldsymbol{r}}_2} + {\alpha _2}^2{ {\boldsymbol{v} }_2}$ 代入式(38)计算 ${{\boldsymbol{s}}^{{\text{(}}i + 1{\text{)}}}}$
15 　结束
16 　使用网格法迭代求解 ${{\boldsymbol{\varepsilon}} ^{{\text{(}}i + 1{\text{)}}}}$
17 　 $i \leftarrow i + 1$
18 达到收敛条件，结束

下载: 导出CSV

| 显示表格

计算复杂度方面，第 $i$ 次迭代求解 ${{\boldsymbol{w}}^{{\text{(}}i + 1{\text{)}}}}$ 和 ${{\boldsymbol{s}}^{{\text{(}}i + 1{\text{)}}}}$ 的运算量主要来自对矩阵 ${\boldsymbol{P}}$ 和 ${\boldsymbol{Q}}$ 的计算，求解 ${{\boldsymbol{\varepsilon}} ^{{\text{(}}i + 1{\text{)}}}}$ 的运算量则可忽略不计。使用式(29)、式(36)给出的定义式运算的计算复杂度分别为 ${{O}}(LK{M^3})$ 和 ${{O}}(LKN{M^2})$ 。若直接计算式(32)、式(39)中的第1项和第2项，可大幅度降低算法复杂度。经由证明，求解 ${\text{tr}}({\boldsymbol{P}})$ 和 ${\text{tr}}({\boldsymbol{Q}})$ 的复杂度为 ${{O}}(LK{M^2})$ ，故本文算法经简化后的计算复杂度为 ${{O}}(LK{M^2})$ ，而引入平方迭代框架不改变计算复杂度。

4. 实验验证

本节设计了基于实测目标数据的实验验证。目标数据采用佐治亚理工学院T-72坦克全极化雷达 $\mathrm{X}$ 波段实测数据^[26]，测量信号为中心频率 ${f_0} = 9.6{\text{ GHz}}$ 的步进频脉冲信号，覆盖带宽为660 MHz，雷达俯仰角固定为30.42°，方位角范围为 $[0^\circ , 360^\circ ]$ 。具体来说，数据库中，转台方位角的变动步长为0.05°，每个数据文件中宽带雷达共进行 79 次扫描，因此，单个文件的方位角变化范围为3.90°，同时，相邻数据文件中心方位角间相差4.25°。这表明在 79 次扫描之后，存在 0.35°的测量盲区。高分辨测量中，单个 $\mathrm{T}\mathrm{I}\mathrm{R}\mathrm{M}$ 的距离支撑长度为 $Q = 37$ ^[19]。非特别说明，仿真中求平均采用的方位角范围 $[{\theta _{{\text{min}}}},{\theta _{{\text{max}}}}]$ 取 $[161.55^\circ ,162.75^\circ ]$ 。

4.1 天线极化矢量固定条件下算法的有效性验证

本小节固定雷达的天线极化方式，分析本文方法的有效性。仿真参数设置为：雷达和干扰机均采用垂直极化 ${\boldsymbol{\varepsilon}} = {[0,1]^{\text{T}}}$ ，雷达参数方面，发射波形时宽为 $T = 81.28{\text{ ns}}$ ，发射脉冲个数为 $K = 4$ ，发射信号的码长为 $N = 40$ ，回波点数为 $M = 76$ ，对应信号带宽 $B = 0.5{\text{ GHz}}$ ，占空比 ${D_{\text{t}}} = 1\%$ ；干扰机参数方面，设置ISRJ占空比为 ${D_{\text{J}}} = 20\%$ ，ISRJ转发周期为 ${T_{\text{J}}} = 20{\text{ ns}}$ ；目标方面，极化散射矩阵设置为 $[1,{\text{0}}{\text{.1}}{{\text{e}}^{{{ - {\rm{j}}\pi }}/3}};{\text{0}}{\text{.1}}{{\text{e}}^{{{ - {\rm{j}}\pi }}/3}},{{\text{e}}^{{{{\rm{j}}\pi }}/3}}]$ ，信噪比损失因子为 $\mu = 1{\text{ dB}}$ ，权重因子 ${\gamma _1} = {\gamma _2} = 100$ ，峰值约束关系为 ${a_{{\text{min}}}} = {a_{{\text{max}}}} \times {10^{ - 4}}$ ，并将 $[{f_1},{f_2}]$ = $[ - 0.3,0.3]$ 离散化为 $L = 10$ 份，预设收敛精度为 ${\eta _\varGamma }$ = ${10^{ - 7}}$ ，设置信噪比为 $15{\text{ dB}}$ ，信干比为 $- 15{\text{ dB}}$ 。

、分别给出了目标函数 $\varGamma$ 、信噪比损失 $\mu$ 随迭代次数的取值变化情况。可以看出，本文算法具有收敛性。给出了优化后的目标回波、 $\mathrm{I}\mathrm{S}\mathrm{R}\mathrm{J}$ 的脉压结果及归一化模糊函数能量图。从可以看出，本文方法的目标回波主瓣峰值约为 $- 7{\text{ dB}}$ ，文献[]方法的目标回波主瓣峰值约为 $- 12{\text{ dB}}$ ；由可见，本文方法的干扰峰值约为 $- 29{\text{ dB}}$ ，文献[]方法的干扰峰值约为 $- 24{\text{ dB}}$ 。这表明，本文方法较好地实现了波形对目标特性的匹配和对ISRJ的抑制。对比、可以看出，目标回波的模糊函数在给定的频移区间内保持了良好的多普勒容忍性；而ISRJ的模糊函数也实现了对给定频移区间的干扰抑制。针对算法的鲁棒性，给出了不同方位角区间 $[{\theta _{{\text{min}}}},{\theta _{{\text{max}}}}]$ 对目标回波脉压峰值的影响，对比TIRM精确已知的情况，跨度1.2°的方位角区间 $[161.55^\circ ,162.75^\circ ]$ 性能下降约2 dB、跨度3.9°的方位角区间 $[161.55^\circ ,165.45^\circ ]$ 性能下降约4 dB、跨度10.0°的方位角区间 $[161.55^\circ , 171.55^\circ ]$ 性能下降约4.5 dB。由此表明：先验方位角范围越大，本文所提算法的脉压增益下降越大，但鲁棒性随之增强。给出了不同脉冲数 $K$ 和信噪比损失 $\mu$ 时，本文算法重复1000次独立实验后旁瓣峰值随信干比的变化情况。图5(a)表明，随信干比增加，旁瓣峰值呈现下降-稳中略升-下降-稳定的趋势，当信干比较大时，回波几乎仅含有目标回波。图5(b)表明，旁瓣峰值呈现单调下降至稳定的趋势，当信噪比很大时，回波几乎不含噪声，而信干比固定，因而趋于稳定。

图 1 目标函数和信噪比损失的收敛情况

下载: 全尺寸图片幻灯片

图 2 目标回波和干扰的脉压性能比较

下载: 全尺寸图片幻灯片

图 3 目标回波和干扰的互模糊函数性能比较

下载: 全尺寸图片幻灯片

图 4 不同方位角范围的鲁棒性算法的归一化脉压峰值

下载: 全尺寸图片幻灯片

图 5 旁瓣峰值随信干比和信噪比的变化曲线

下载: 全尺寸图片幻灯片

4.2 天线极化矢量优化条件下算法的有效性验证

本小节引入天线极化矢量优化，分析极化优化对目标回波和ISRJ脉压性能的影响。仿真过程中设置 $N = M = 80$ ，将极化旋转角 $\alpha \in [ - {\pi}/2,{\pi}/2]$ 、极化椭圆角 $\beta \in [ - {\pi}/4,{\pi}/4]$ 等分为50份，步长分别为 $0.02{\pi}$ , $0.01{\pi}$ ，其他参数设置同4.1节。干扰机天线极化方式依次考虑以下几种设置：60°线极化 ${{\boldsymbol{\varepsilon}} _{\text{J}}} = {\left[ {\dfrac{{\sqrt 2 }}{2}{{\text{e}}^{{\text{j}}\frac{{\pi}}{3}}},\dfrac{{\sqrt 2 }}{2}{{\text{e}}^{{\text{j}}\frac{{\pi}}{3}}}} \right]^{\text{T}}}$ 、左旋圆极化 ${{\boldsymbol{\varepsilon}} _{\text{J}}} = \left[ \dfrac{{\sqrt 2 }}{2}{{\text{e}}^{{\text{j}}\frac{{5{\pi}}}{6}}}, \dfrac{{\sqrt 2 }}{2}{{\text{e}}^{{\text{j}}\frac{{\pi}}{3}}} \right]^{\text{T}}$ 、右旋圆极化 ${{\boldsymbol{\varepsilon}} _{\text{J}}} = {\left[ {\dfrac{{\sqrt 2 }}{2}{{\text{e}}^{ - {\text{j}}\frac{{\pi}}{6}}},\dfrac{{\sqrt 2 }}{2}{{\text{e}}^{{\text{j}}\frac{{\pi}}{3}}}} \right]^{\text{T}}}$ 。

图6依次给出了干扰机在3种不同极化工作方式下，进行天线极化矢量优化前后的目标回波脉压归一化幅度图。可以看出，经极化矢量优化后目标回波脉压结果对比优化前旁瓣平均下降了13～19 dB。图7则给出了干扰机在不同极化工作方式下，进行天线极化矢量优化前后的ISRJ信号脉压归一化幅度图。可以看出，经极化矢量优化后目标的ISRJ信号脉压峰值对比优化前有不同程度的下降，其中在干扰机采用60°线极化时，干扰脉压峰值近似相同，整体水平却有所下降；在干扰机采用左旋圆极化和右旋圆极化时，干扰脉压后整体下降至非常低的水平(此时近似被完全滤除)。综上说明了经优化后的天线极化矢量在获得更好目标位置检测能力的同时，显著抑制了间歇采样干扰信号，进一步增强了宽带场景下的雷达探测能力。

图 6 不同干扰机天线极化方式设置条件下目标回波脉压幅度图

下载: 全尺寸图片幻灯片

图 7 不同干扰机天线极化方式设置条件下ISRJ回波脉压幅度图

下载: 全尺寸图片幻灯片

4.3 本文方法的抗ISRJ性能分析

为进一步验证本文方法在宽带全极化雷达条件下抗ISRJ的有效性，本节考虑干扰机高速运动，即设置干扰机与雷达间距离为 ${d_{\text{J}}} = 1100{\text{ m}}$ ，干扰机速度 ${v_{\text{J}}} = 30{\text{ m}}/{\text{s}}$ ；目标距离雷达 $d = 600{\text{ m}}$ ，目标速度 $v = 40{\text{ m}}/{\text{s}}$ ，其余参数同4.1节，对比本文方法与文献[7]方法在目标探测和干扰抑制方面的差异。

给出了不同脉冲数 $K$ 对应的回波脉压效果。根据可以看出，在 $K = 2$ 时，本文方法对目标回波有约 $- 9{\text{ dB}}$ 的脉压增益(较文献[]方法高出约 $8{\text{ dB}}$ )；从可看出在 $K = 4$ 时目标回波增益较文献[]方法高出约6 dB；从则可看出在 $K = 6$ 时目标回波增益较文献[7]方法高出约4 dB。而对于ISRJ信号，则在3种场景下都近似被完全滤除，图8结果进一步表明，本文方法相比文献[7]的目标探测性能更好，同时对ISRJ信号起到了显著的抑制效果。

图 8 不同脉冲数条件下本文方法与文献[7]方法的脉压结果对比

下载: 全尺寸图片幻灯片

5. 结语

本文对宽带全极化雷达发射波形、接收滤波器和天线极化进行联合优化以实现间歇采样转发干扰的抑制。对于多分量非凸优化问题，基于交替迭代求解方法，依次对最优发射波形、接收滤波器和天线极化进行了有效求解。对比现有方法，本文方法不仅实现了发射波形-目标的相互匹配，还充分利用了目标和干扰机的极化域信息，使得脉压后显著抑制了ISRJ并保证良好的目标回波增益，同时其多普勒容限显著增强了算法对高速干扰机的适用性，对解决高分辨场景下抗间歇采样转发干扰问题具有一定的实际意义和应用价值。而下一步将重点研究 $\mathrm{I}\mathrm{S}\mathrm{R}\mathrm{J}$ 参数不确定条件下的波形鲁棒性设计问题，以降低对先验信息的精度要求。

图 1 间隙波导理论和实际模型

下载: 全尺寸图片幻灯片

图 2 典型的间隙波导传输线

下载: 全尺寸图片幻灯片

图 3 近年来所提出的几种重要的新型间隙波导非接触EBG结构

下载: 全尺寸图片幻灯片

图 4 基于间隙波导技术的新型毫米波传输线

下载: 全尺寸图片幻灯片

图 5 基于间隙波导的非接触法兰及扩展应用

下载: 全尺寸图片幻灯片

图 6 基于间隙波导技术的新型无源器件

下载: 全尺寸图片幻灯片

图 7 基于脊间隙波导传输线的无基片功率放大器

下载: 全尺寸图片幻灯片

图 8 间隙波导阵列天线

下载: 全尺寸图片幻灯片

图 9 S频段折叠小型化非接触式低PIM法兰及PIM实测结果

下载: 全尺寸图片幻灯片

图 10 基于间隙波导技术的堆叠式毫米波电路系统

下载: 全尺寸图片幻灯片

表 1 间隙波导国内外研究领域分类

研究领域	细节分类
非接触EBG结构	理论分析；新型非接触EBG；其他
GW传输线	经典GW传输线及应用；新型GW传输线及应用；过渡、互联方法
GW的微波毫米波技术应用	无源器件：非接触法兰、滤波器、耦合器、移相器、旋转关节等；有源电路：无基片有源电路；天线：缝隙阵列天线、GW喇叭天线、其他新型天线等；其他：电路封装、快速测试、3D打印、先进制造等

下载: 导出CSV

参考文献(57)

[1]	徐常志, 靳一, 李立, 等. 面向6G的星地融合无线传输技术[J]. 电子与信息学报, 2021, 43(1): 28–36. doi: 10.11999/JEIT200363 XU Changzhi, JIN Yi, LI Li, et al. Wireless transmission technology of satellite-terrestrial integration for 6G mobile communication[J]. Journal of Electronics &Information Technology, 2021, 43(1): 28–36. doi: 10.11999/JEIT200363
[2]	KILDAL P S, ALFONSO E, VALERO-NOGUEIRA A, et al. Local metamaterial-based waveguides in gaps between parallel metal plates[J]. IEEE Antennas and Wireless Propagation Letters, 2009, 8: 84–87. doi: 10.1109/LAWP.2008.2011147
[3]	MOROZOV G V and SPRUNG D W L. Floquet-Bloch waves in one-dimensional photonic crystals[J]. Europhysics Letters, 2011, 96(5): 54005. doi: 10.1209/0295-5075/96/54005
[4]	BERENGUER A, FUSCO V, ZELENCHUK D E, et al. Propagation characteristics of groove gap waveguide below and above cutoff[J]. IEEE Transactions on Microwave Theory and Techniques, 2016, 64(1): 27–36. doi: 10.1109/TMTT.2015.2504501
[5]	KILDAL P S, ZAMAN A U, RAJO-IGLESIAS E, et al. Design and experimental verification of ridge gap waveguide in bed of nails for parallel-plate mode suppression[J]. IET Microwaves, Antennas & Propagation, 2011, 5(3): 262–270. doi: 10.1049/iet-map.2010.0089
[6]	RAZA H, YANG Jian, KILDAL P S, et al. Microstrip-ridge gap waveguide–study of losses, bends, and transition to WR-15[J]. IEEE Transactions on Microwave Theory and Techniques, 2014, 62(9): 1943–1952. doi: 10.1109/TMTT.2014.2327199
[7]	BRAZÁLEZ A A, RAJO-IGLESIAS E, VÁZQUEZ-ROY J L, et al. Design and validation of microstrip gap waveguides and their transitions to rectangular waveguide, for millimeter-wave applications[J]. IEEE Transactions on Microwave Theory and Techniques, 2015, 63(12): 4035–4050. doi: 10.1109/TMTT.2015.2495141
[8]	RAZA H, YANG Jian, KILDAL P S, et al. Resemblance between gap waveguides and hollow waveguides[J]. IET Microwaves, Antennas & Propagation, 2013, 7(15): 1221–1227. doi: 10.1049/iet-map.2013.0178
[9]	NASR M A and KISHK A A. Vertical coaxial-to-ridge waveguide transitions for ridge and ridge gap waveguides with 4: 1 bandwidth[J]. IEEE Transactions on Microwave Theory and Techniques, 2019, 67(1): 86–93. doi: 10.1109/TMTT.2018.2873312
[10]	BIRGERMAJER S, JANKOVIĆ N, RADONIĆ V, et al. Microstrip-ridge gap waveguide filter based on cavity resonators with mushroom inclusions[J]. IEEE Transactions on Microwave Theory and Techniques, 2018, 66(1): 136–146. doi: 10.1109/TMTT.2017.2750149
[11]	SHI Yongrong, ZHANG Junzhi, ZENG Sheng, et al. Novel W-band millimeter-wave transition from microstrip line to groove gap waveguide for MMIC integration and antenna application[J]. IEEE Transactions on Antennas and Propagation, 2018, 66(6): 3172–3176. doi: 10.1109/TAP.2018.2819902
[12]	PENG Songtao, PU Youlei, WU Zewei, et al. A broadband transition from ridge gap waveguide to microstrip using suspended line coupling[J]. IEEE Microwave and Wireless Components Letters, 2021, 31(3): 253–256. doi: 10.1109/LMWC.2020.3040651
[13]	SHI Yongrong, FENG Wenjie, WANG Hao, et al. Novel W-band LTCC transition from microstrip line to ridge gap waveguide and its application in 77/79 GHz antenna array[J]. IEEE Transactions on Antennas and Propagation, 2019, 67(2): 915–924. doi: 10.1109/TAP.2018.2882625
[14]	MORALES-HERNÁNDEZ A, FERRANDO-ROCHER M, SÁNCHEZ-SORIANO M Á, et al. Design strategy and considerations to improve corona discharge breakdown in groove gap waveguides[C]. IEEE 15th European Conference on Antennas and Propagation (EuCAP), Dusseldorf, Germany, 2021: 1–5.
[15]	FAN Fangfang, YANG Jian, VASSILEV V, et al. Bandwidth investigation on half-height pin in ridge gap waveguide[J]. IEEE Transactions on Microwave Theory and Techniques, 2018, 66(1): 100–108. doi: 10.1109/TMTT.2017.2732983
[16]	EBRAHIMPOURI M, RAJO-IGLESIAS E, SIPUS Z, et al. Cost-effective gap waveguide technology based on glide-symmetric holey EBG structures[J]. IEEE Transactions on Microwave Theory and Techniques, 2018, 66(2): 927–934. doi: 10.1109/TMTT.2017.2764091
[17]	SUN Dongquan, CHEN Xiang, DENG Jingya, et al. Gap waveguide with interdigital-pin bed of nails for high-frequency applications[J]. IEEE Transactions on Microwave Theory and Techniques, 2019, 67(7): 2640–2648. doi: 10.1109/TMTT.2019.2914907
[18]	SUN Dongquan, CHEN Xiang, and GUO Lixin. Compact corrugated plate for double-sided contactless waveguide flange[J]. IEEE Microwave and Wireless Components Letters, 2021, 31(2): 129–132. doi: 10.1109/LMWC.2020.3042279
[19]	SUN Dongquan, CHEN Xiang, GUO Lixin, et al. Hard–soft groove gap waveguide based on perpendicularly stacked corrugated metal plates[J]. IEEE Transactions on Microwave Theory and Techniques, 2021, 69(8): 3684–3692. doi: 10.1109/TMTT.2021.3086497
[20]	PENG Songtao, PU Youlei, WU Zewei, et al. Embedded bed of nails with robustness suitable for broadband gap waveguide technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2021, 69(12): 5317–5326. doi: 10.1109/TMTT.2021.3116178
[21]	BAYAT-MAKOU N and KISHK A A. Contactless air-filled substrate integrated waveguide[J]. IEEE Transactions on Microwave Theory and Techniques, 2018, 66(6): 2928–2935. doi: 10.1109/TMTT.2018.2818137
[22]	ZHAO Xiaofei, DENG Jingya, YIN Jiayuan, et al. Novel suspended-line gap waveguide packaged with stacked-mushroom EBG structures[J]. IEEE Transactions on Microwave Theory and Techniques, 2021, 69(5): 2447–2457. doi: 10.1109/TMTT.2021.3068260
[23]	PUCCI E and KILDAL P S. Contactless non-leaking waveguide flange realized by bed of nails for millimeter wave applications[C]. IEEE 6th European Conference on Antennas and Propagation (EUCAP), Prague, Czech Republic, 2012: 3533–3536.
[24]	孙冬全. 毫米波间隙波导技术及FMCW反射功率对消系统应用研究[D]. [博士论文], 东南大学, 2017. SUN Dongquan. Research on millimeter-wave gap waveguide technology and FMCW reflected power cancelation systems[D]. [Ph. D. dissertation], Southeast University, 2017.
[25]	RAHIMINEJAD S, PUCCI E, VASSILEV V, et al. Polymer gap adapter for contactless, robust, and fast measurements at 220–325 GHz[J]. Journal of Microelectromechanical Systems, 2016, 25(1): 160–169. doi: 10.1109/JMEMS.2015.2500277
[26]	SUN Dongquan, CHEN Zhenhua, YAO Changfei, et al. Flexible rectangular waveguide based on cylindrical contactless flange[J]. Electronics Letters, 2016, 52(25): 2042–2044. doi: 10.1049/el.2016.3536
[27]	SUN Dongquan and XU Jinping. Real time rotatable waveguide twist using contactless stacked air-gapped waveguides[J]. IEEE Microwave and Wireless Components Letters, 2017, 27(3): 215–217. doi: 10.1109/LMWC.2017.2661881
[28]	AHMADI B and BANAI A. Direct coupled resonator filters realized by gap waveguide technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2015, 63(10): 3445–3452. doi: 10.1109/TMTT.2015.2457916
[29]	XIU Tao, YAO Yuan, JIANG Hang, et al. Design of a compact and low-loss E-band filter based on multilayer groove gap waveguide[J]. IEEE Microwave and Wireless Components Letters, 2021, 31(11): 1211–1214. doi: 10.1109/LMWC.2021.3111955
[30]	SUN Dongquan and XU Jinping. A novel iris waveguide bandpass filter using air gapped waveguide technology[J]. IEEE Microwave and Wireless Components Letters, 2016, 26(7): 475–477. doi: 10.1109/LMWC.2016.2574822
[31]	REZAEE M and ZAMAN A U. Groove gap waveguide filter based on horizontally polarized resonators for V-band applications[J]. IEEE Transactions on Microwave Theory and Techniques, 2020, 68(7): 2601–2609. doi: 10.1109/TMTT.2020.2986111
[32]	SUN Dongquan and XU Jinping. Rectangular waveguide coupler with adjustable coupling coefficient using gap waveguide technology[J]. Electronics Letters, 2017, 53(3): 167–169. doi: 10.1049/el.2016.4039
[33]	PALOMARES-CABALLERO Á, ALEX-AMOR A, ESCOBEDO P, et al. Low-loss reconfigurable phase shifter in gap-waveguide technology for mm-wave applications[J]. IEEE Transactions on Circuits and Systems II:Express Briefs, 2020, 67(12): 3058–3062. doi: 10.1109/TCSII.2020.3000058
[34]	SÁNCHEZ-ESCUDEROS D, HERRANZ-HERRUZO J I, FERRANDO-ROCHER M, et al. True-time-delay mechanical phase shifter in gap waveguide technology for slotted waveguide arrays in Ka-band[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(5): 2727–2740. doi: 10.1109/TAP.2020.3030993
[35]	陈翔, 孙冬全, 李小军, 等. U型超宽带非接触式波导旋转关节、控制系统、方法及应用[P]. 中国专利, 111934062A, 2020. CHEN Xiang, SUN Dongquan, LI Xiaojun, et al. U-shaped ultra-wideband non-contact waveguide rotary joint, control system, method and application[P]. China Patent, 111934062A, 2020.
[36]	FARAHBAKHSH A. Wideband rotary joint based on gap waveguide technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2021, 69(10): 4385–4391. doi: 10.1109/TMTT.2021.3090988
[37]	HORESTANI A K, SHATERIAN Z, and MROZOWSKI M. High dynamic range microwave displacement and rotation sensors based on the phase of transmission in groove gap waveguide technology[J]. IEEE Sensors Journal, 2022, 22(1): 182–189. doi: 10.1109/JSEN.2021.3130658
[38]	WANG Minxing, WU Zewei, LIAO Xiaoyi, et al. Exciting circular TM₁₁ mode using symmetric probes based on ridge gap waveguide[J]. IEEE Transactions on Microwave Theory and Techniques, 2022, 70(1): 334–342. doi: 10.1109/TMTT.2021.3121324
[39]	TAMAYO-DOMÍNGUEZ A, FERNÁNDEZ-GONZÁLEZ J M, and SIERRA-CASTAÑER M. 3-D-printed modified butler matrix based on gap waveguide at W-band for monopulse radar[J]. IEEE Transactions on Microwave Theory and Techniques, 2020, 68(3): 926–938. doi: 10.1109/TMTT.2019.2953164
[40]	SHI Yongrong, ZHANG Junzhi, ZHOU Ming, et al. Miniaturized W-band gap waveguide bandpass filter using the MEMS technique for both waveguide and surface mounted packaging[J]. IEEE Transactions on Circuits and Systems II:Express Briefs, 2019, 66(6): 938–942. doi: 10.1109/TCSII.2018.2873236
[41]	FARJANA S, GHADERI M, ZAMAN A U, et al. Low-loss gap waveguide transmission line and transitions at 220–320 GHz using dry film micromachining[J]. IEEE Transactions on Components, Packaging and Manufacturing Technology, 2021, 11(11): 2012–2021. doi: 10.1109/TCPMT.2021.3111137
[42]	AHMADI B and BANAI A. Substrateless amplifier module realized by ridge gap waveguide technology for millimeter-wave applications[J]. IEEE Transactions on Microwave Theory and Techniques, 2016, 64(11): 3623–3630. doi: 10.1109/TMTT.2016.2607177
[43]	FERRANDO-ROCHER M, HERRANZ-HERRUZO J I, VALERO-NOGUEIRA A, et al. 8×8 Ka-band dual-polarized array antenna based on gap waveguide technology[J]. IEEE Transactions on Antennas and Propagation, 2019, 67(7): 4579–4588. doi: 10.1109/TAP.2019.2908109
[44]	SHI Yongrong, FENG Wenjie, and CAO Baolin. W-band gap waveguide antenna array: Passive/active component gap waveguide transition interface for system integration[J]. IEEE Antennas and Propagation Magazine, 2021, 63(2): 40–49. doi: 10.1109/MAP.2019.2943306
[45]	LIU Ying, YUE Zhenzhen, JIA Yongtao, et al. Dual-band dual-circularly polarized antenna array with printed ridge gap waveguide[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(8): 5118–5123. doi: 10.1109/TAP.2020.3048504
[46]	VILENSKIY A R, MAKURIN M N, LEE C, et al. Reconfigurable transmitarray with near-field coupling to gap waveguide array antenna for efficient 2-D beam steering[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(12): 7854–7865. doi: 10.1109/TAP.2020.2998904
[47]	TAMAYO-DOMÍNGUEZ A, FERNÁNDEZ-GONZÁLEZ J M, and SIERRA-CASTAÑER M. Monopulse radial line slot array antenna fed by a 3-D-printed cavity-ended modified butler matrix based on gap waveguide at 94 GHz[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(8): 4558–4568. doi: 10.1109/TAP.2021.3060045
[48]	LIU Jinlin, YANG Fei, FAN Kuikui, et al. Unequal power divider based on inverted microstrip gap waveguide and its application for low sidelobe slot array antenna at 39 GHz[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(12): 8415–8425. doi: 10.1109/TAP.2021.3096981
[49]	QUAN Yu, WANG Hao, TAO Shifei, et al. A double-layer multibeam antenna with 45° linear polarization based on gap waveguide technology[J]. IEEE Transactions on Antennas and Propagation, 2022, 70(1): 56–66. doi: 10.1109/TAP.2021.3090507
[50]	SUN Dongquan and XU Jinping. Compact phase corrected H-plane horn antenna using slow-wave structures[J]. IEEE Antennas and Wireless Propagation Letters, 2017, 16: 1032–1035. doi: 10.1109/LAWP.2016.2618843
[51]	MOHAMMADPOUR M, MOHAJERI F, and RAZAVI S A. A new wide band and compact H-plane horn antenna based on groove gap waveguide technology[J]. IEEE Transactions on Antennas and Propagation, 2022, 70(1): 221–228. doi: 10.1109/TAP.2021.3111342
[52]	YUAN Wei, CHEN Jianfeng, ZHANG Cheng, et al. Glide-symmetric lens antenna in gap waveguide technology[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(4): 2612–2620. doi: 10.1109/TAP.2019.2955919
[53]	BRAZALEZ A A, ZAMAN A U, and KILDAL P S. Improved microstrip filters using PMC packaging by lid of nails[J]. IEEE Transactions on Components, Packaging and Manufacturing Technology, 2012, 2(7): 1075–1084. doi: 10.1109/TCPMT.2012.2190931
[54]	ASHRAF N, SEBAK A R, and KISHK A A. PMC packaged single-substrate 4 × 4 butler matrix and double-ridge gap waveguide horn antenna array for multibeam applications[J]. IEEE Transactions on Microwave Theory and Techniques, 2021, 69(1): 248–261. doi: 10.1109/TMTT.2020.3022092
[55]	CHEN Xiang, SUN Dongquan, CUI Wanzhao, et al. A folded contactless waveguide flange for low passive-intermodulation applications[J]. IEEE Microwave and Wireless Components Letters, 2018, 28(10): 864–866. doi: 10.1109/LMWC.2018.2865506
[56]	JIANG Xun, JIA Fangxiu, SHI Yongrong, et al. Ridge gap waveguide layer transition for compact 3-D waveguide packaging application[J]. IEEE Transactions on Components, Packaging and Manufacturing Technology, 2019, 9(10): 2136–2139. doi: 10.1109/TCPMT.2019.2937436
[57]	WANG Lin, DING Dazhi, CHEN Rushan, et al. Transient analysis of high-power microwave air breakdown under external DC magnetic field[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(6): 4894–4903. doi: 10.1109/TAP.2020.2969894

施引文献

期刊类型引用(1)
1. 王颖，郭睿，梁毅. 实测双极化雷达压制干扰的特性分析与抑制. 海军航空大学学报. 2025(01): 163-170+196 . 百度学术
其他类型引用(1)

资源附件(0)

访问统计

图(10) / 表(1)

计量

文章访问数: 3436
HTML全文浏览量: 2401
PDF下载量: 483
被引次数: 2

1. 引言
2. 问题模型
3. 问题求解
4. 实验验证
4.1 天线极化矢量固定条件下算法的有效性验证
4.2 天线极化矢量优化条件下算法的有效性验证
4.3 本文方法的抗ISRJ性能分析
5. 结语

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

间隙波导技术及其空间应用

doi: 10.11999/JEIT211291

作者简介:
陈翔：男，高级工程师，博士，研究方向为航天器微波毫米波技术

孙冬全：男，副教授，博士，研究方向为微波毫米波电路、天线理论及技术

崔万照：男，研究员，博士，研究方向为航天器微波技术

通讯作者:
孙冬全　dqsun87@163.com

计量

Gap Waveguide Technology and its Space Applications

1. 引言

2. 问题模型

3. 问题求解

4. 实验验证

4.1 天线极化矢量固定条件下算法的有效性验证

4.2 天线极化矢量优化条件下算法的有效性验证

4.3 本文方法的抗ISRJ性能分析

5. 结语

期刊类型引用(1)

其他类型引用(1)

计量

目录

1. 引言

2. 问题模型

3. 问题求解

4. 实验验证

4.1 天线极化矢量固定条件下算法的有效性验证

4.2 天线极化矢量优化条件下算法的有效性验证

4.3 本文方法的抗ISRJ性能分析

5. 结语

留言板

间隙波导技术及其空间应用

doi: 10.11999/JEIT211291

作者简介: 陈翔：男，高级工程师，博士，研究方向为航天器微波毫米波技术 孙冬全：男，副教授，博士，研究方向为微波毫米波电路、天线理论及技术 崔万照：男，研究员，博士，研究方向为航天器微波技术

通讯作者: 孙冬全 dqsun87@163.com

计量

出版历程

Gap Waveguide Technology and its Space Applications

1. 引言

2. 问题模型

3. 问题求解

4. 实验验证

4.1 天线极化矢量固定条件下算法的有效性验证

4.2 天线极化矢量优化条件下算法的有效性验证

4.3 本文方法的抗ISRJ性能分析

5. 结语

期刊类型引用(1)

其他类型引用(1)

计量

出版历程

目录

1. 引言

2. 问题模型

3. 问题求解

4. 实验验证

4.1 天线极化矢量固定条件下算法的有效性验证

4.2 天线极化矢量优化条件下算法的有效性验证

4.3 本文方法的抗ISRJ性能分析

5. 结语

作者简介:
陈翔：男，高级工程师，博士，研究方向为航天器微波毫米波技术

孙冬全：男，副教授，博士，研究方向为微波毫米波电路、天线理论及技术

崔万照：男，研究员，博士，研究方向为航天器微波技术

通讯作者:
孙冬全　dqsun87@163.com