DOA Estimation Under Active Deception Jamming Environment
-
摘要:
针对有源欺骗干扰环境下基于小样本的DOA估计问题，该文提出自适应极化滤波(APF)联合块稀疏贝叶斯学习(BSBL)算法的DOA估计方法。首先，通过APF抑制干扰能量，提高信干比。然后，建立有源欺骗干扰环境下的稀疏贝叶斯模型，基于相邻快拍相关性，利用BSBL算法进行DOA估计。仿真和实测数据处理结果表明，所提方法降低了干扰对BSBL算法的影响，且与APF联合子空间类算法或最大似然算法(ML)相比，具有更高的空间分辨率和DOA估计精度。
Abstract:For the target DOA estimation under active deception jamming environment with limited samples, a novel DOA estimation method based on the combination of Adaptive Polarization Filter(APF) and Block Sparse Bayesian Learning(BSBL) algorithm is proposed. First, the interference energy is suppressed using APF. Then, the proposed method constructs a sparse Bayesian model under active deception jamming environment. The target DOA is estimated using the BSBL algorithm based on the neighbor time sampling correlation. Simulated and measured data processing results prove that the proposed method reduces the influence of interference on the BSBL algorithm, and has higher spatial resolution and higher angle measurement accuracy, comparing with the method based on the combination of APF and subspace-based DOA algorithms or maximum likelihood DOA algorithm.
-
Key words:
- DOA estimation /
- Anti-jamming /
- Adaptive polarization filter /
- Sparse Bayesian learning
-
1. 引言
现代战场环境复杂，电子战装备呈现高度集成化、综合化发展趋势，释放干扰形式多样，特别是数字射频存储器(Digital Radio Frequency Memory, DRFM)技术在雷达对抗领域中的广泛应用[1,2]，使得干扰信号与雷达信号在多个参数域重叠，导致常用时频域抗干扰方法失效。在如此严峻的干扰环境下，目标源DOA估计方法的性能将大打折扣。此外，目标源的非合作性导致有效样本数少。所以，研究有源欺骗干扰环境下，基于小样本的DOA估计方法具有重要意义。
常规的DOA估计算法有两类，第1类是最大似然DOA估计算法，在高信噪比情况下可接近克拉美罗界，但运算量非常大。虽然存在很多优化算法应用到最大似然DOA估计中以降低运算量[3–5]，但还是面临寻优速度慢、精度不够的问题；第2类是子空间类算法，这类算法基于空间搜索，运算量大，而且不适用于相干信号源DOA估计。虽然存在很多算法用于改善相干源下子空间类算法的DOA估计性能[6,7]，但多数算法基于特定的阵列结构，限制了算法的实用性。近几年兴起的稀疏重构类算法为DOA估计开辟了新思路[8–11]。相比于其他稀疏重构类算法，稀疏贝叶斯学习算法具备较大优势。首先，稀疏贝叶斯学习算法的全局最小点通常是稀疏解，而广泛使用的基于L1惩罚项的算法在限定条件下才可以满足[12]；其次，当感知矩阵列相关性很强时，大多数稀疏重构类算法失效，但稀疏贝叶斯学习算法仍有较好性能[13]。稀疏贝叶斯DOA估计通常将DOA估计问题转化为基于多测量矢量(Multiple Measurement Vectors, MMV)模型的稀疏重构问题，根据贝叶斯准则获取后验概率函数，利用EM(Evidence Maximization)等算法得到超参数估计，进而估计来波方向。基于块稀疏贝叶斯学习(Block Sparse Bayesian Learning, BSBL)[13–15]的DOA估计利用Kronecker积将MMV模型转化为具有块结构的SMV(Single Measurement Vector, SMV)模型，MMV模型下的非零行对应SMV模型下的非零块，再利用稀疏贝叶斯学习得到DOA估计。因为，BSBL算法既利用了统计信息又利用了结构信息，所以DOA估计的稳健性得以提高。
以上DOA估计算法工作在无干扰环境时性能较好，但有源欺骗干扰的存在会降低算法的DOA估计性能[16,17]。极化域作为继时域、频域、空域后的又一重要信号域，若可以充分利用雷达信号和干扰信号在该域的差异，抗干扰性能将得以提升。Nathanson[18]早先提出自适应极化对消器，用于研究雨杂波对消问题，之后，出现实时性较好的自适应极化滤波(Adaptive Polarization Filter, APF)包括极化状态估计和最佳极化接收两部分，通过调整接收天线极化形式，使其与干扰极化形式正交，从而抑制干扰信号，提高信干比[19,20]。
针对有源欺骗干扰环境下基于小样本的DOA估计问题，本文将APF与BSBL算法联合，实现了干扰环境下的高精度DOA估计。首先通过APF抑制干扰，提高信干比。之后，建立有源欺骗干扰环境下的稀疏贝叶斯模型，利用相邻采样时间相关性，将MMV模型转化为具有块稀疏特性的SMV模型，采用BSBL算法进行DOA估计。在仿真中，先在不同信噪比条件下，分析了所提算法的DOA估计性能与干扰和目标的极化角度差以及算法迭代次数的关系。然后，在有源欺骗干扰、小样本环境下，将所提方法与BSBL算法、APF联合子空间类算法和APF联合最大似然DOA估计算法进行对比，仿真和实测数据处理结果验证了所提方法具有较高的空间分辨率和测角精度。
2. 有源欺骗干扰环境下信号模型
正交极化双通道测量是获取入射电磁波极化信息的物理前提，到达接收天线的信号极化矢量表示为
h=[hh,hv]=[cos(α)sin(α)e(jϕ)] (1) 其中，h和v表示接收天线正交极化双通道，
α 和ϕ 决定信号的极化状态。雷达信号的极化矢量表示为
hq=[cos(αq)sin(αq)e(jϕq)],q=1,2,···,Q (2) 其中，
Q 为目标数目。干扰信号极化矢量表示为h1k=[cos(α1k)sin(α1k)e(jϕk)],k=1,2,···,K (3) 其中，
K 为干扰数目。弹载雷达接收信号由3部分构成：目标信号、干扰信号和噪声信号。目标信号表示为
Sq=[e[j2πλdsin(θq)(0:M−1)T]gq(t)]⊗hq,q=1,2,···,Q (4) 其中，
θq 表示目标方位角，d为均匀线阵的阵元间距，M 为阵元数。gq(t) 表示第q 个目标的回波复包络，⊗ 表示Kronecker积运算。雷达接收到的干扰信号表示为Jk=[e[j2πλdsin(θ1k)(0:M−1)T]jk(t)]⊗h1k,k=1,2,···,K (5) 其中，
θ1k 表示干扰方位角，jk(t) 表示第k 个干扰信号的复包络。天线接收到的总信号X 表示为X=Q∑q=1Sq+K∑k=1Jk+Ns (6) 其中，
Ns 表示互不相关的加性高斯白噪声。3. APF联合BSBL算法的DOA估计方法
弹载雷达接收信号
X 中包含多个有源欺骗干扰信号，在有限信噪比和迭代次数下，影响的BSBL算法DOA估计性能。将APF联合BSBL算法可提高BSBL算法在有源欺骗干扰环境下DOA估计性能。APF抑制干扰的前提是知道干扰信号的极化信息。记接收电场的极化相干矩阵为
C ，表示如式(7)C=[ChhChvCvhCvv] (7) 极化信息可完全由
C 决定，实际应用中不可先验得到，它的最大似然估计为ˆC=1NN∑n=1XHnXn (8) 其中，
N 为回波样本数，Xn=[Xnh,Xnv] ，Xnh 和Xnv 分别表示Xh 和Xv 的第n 列。假设n 时刻极化相干矩阵估计为ˆCn ,n+1 时刻的极化相干矩阵估计ˆCn+1 通过迭代得到，具体步骤如下：(1)初始值
ˆC0=[1j−j1] 。(2)
ˆCn+1 的迭代公式如式(9)ˆCn+1=(1−ζ)ˆCn+ζXHn+1Xn+1 (9) 其中，
n=0,1,⋯,N−1 。0≤ζ≤1 为新数据权重因子，反映新数据对相干矩阵估计的重要性。接收回波样本数据量越大，每个新样本对ˆCn+1 的作用越小。所以，可以认为与样本数N 具有反比关系，取为1/N 。极化相干矩阵
ˆC 迭代完成后，对其进行特征分解，小特征值对应的特征矢量记为u2 。天线最佳接收极化矢量表示为hopt=u2T‖ (10) 获得天线最佳接收极化矢量
{{{h}}_{{\rm{opt}}}} 后，对雷达接收信号{{X}} 进行干扰滤除得{{{X}}^{\rm{A}}} \begin{array}{*{20}{c}} {{{X}}_i^{\rm{A}} = {{{X}}_i}{{h}}_{{\rm{opt}}}^{\rm{T}}, }&{i = 1, 2, ·\!·\!· , N} \end{array} (11) 其中，
{{X}}_i^{\rm{A}} 表示{{{X}}^{\rm{A}}} 的第i 列。利用APF提高信干比得到
{{{X}}^{\rm{A}}} 后，通过BSBL算法进行DOA估计。BSBL算法的数学模型是将MMV模型转化为具有块稀疏特性的SMV模型。{{{X}}^{\rm{A}}} 在基于MMV的稀疏贝叶斯模型中表示为{{{X}}^{\rm{A}}} = {{fs}} + {{V}} (12) 其中，
{{f}} 是M \times L 的学习字典，L 是空间划分格数。{{s}} 是L \times N 维的待恢复稀疏矩阵，为了得到唯一的全局最优解，应保证{{s}} 是行稀疏的。{{V}} 是M \times N 维的噪声矢量。无角度欺骗干扰时，{{s}} 的元素满足{{s}}\left( {l, :} \right) = \left\{ {\begin{array}{*{20}{l}} {{{{z}}_q}, } & {{\theta _l} = {\theta _q}}\\ {0, } & {{其它}} \end{array}} \right. , {l = 1, 2, ·\!·\!· , L} (13) 其中，
{{{z}}_q} 表示目标源信号。有角度欺骗干扰时，{{s}} 的元素满足{{s}}\left( {l, :} \right) = \left\{ {\begin{array}{*{20}{l}} {{{{z}}_q}, } & {{\theta _l} = {\theta _q}} \\ {{{{z}}_k}, } & {{\theta _l} = \theta _k^1} \\ {0, } & {{其它}} \end{array}} \right., {l = 1, 2, ·\!·\!· , L} (14) 其中，
{{{z}}_k} 表示干扰转发信号。令
{{y}} = {\rm{vec}}({({{{X}}^{\rm{A}}})^{\rm{T}}}) (15) {{D}} = {{f}} \otimes {{{I}}_N} \hspace{20pt} (16) {{x}} = {\rm{vec}}({{{s}}^{\rm{T}}}) \hspace{20pt} (17) {{v}} = {\rm{vec}}({{{V}}^{\rm{T}}}) \hspace{18pt} (18) 其中，
{\rm{vec}}( \bullet ) 代表矩阵矢量化，{{{I}}_N} 表示N 维单位阵。利用式(15)到式(18)，MMV模型转化为SMV模型，表示为
{{y}} = {{Dx}} + {{v}} (19) 式(19)的等价表示为
\begin{aligned} {{y}} =& \left[ {\begin{aligned} {{{{f}}_1} \otimes {{{I}}_N}}\;\; {{{{f}}_2} \otimes {{{I}}_N}}\;\;·\!·· &{{{{f}}_L} \otimes {{{I}}_N}} \end{aligned}}\right]\\ &\cdot{\left[ {\begin{aligned} {{{x}}_1^{\rm{T}}}\;\; {{{x}}_2^{\rm{T}}}\;\;·\!·\!· \;\;{{{x}}_L^{\rm{T}}} \end{aligned}} \right]^{\rm{T}}} + {{v}} \end{aligned} (20) 其中，
{{x}}_i^{\rm{T} } 表示{{s}} 的第i 行，{{{f}}_i} 表示{{f}} 的第i 列，i = 1, 2, ·\!·\!· , L 。由展开式可以看出，{{x}} 由L 个块组成，每个块中包含N 个元素，由于接收信号采样频率较高，N 个元素间通常具有时间相关性。假设噪声矢量中元素独立且同服从参数为\partial 的高斯分布，{{y}} 的似然函数为p({{y}}|{{x}};\partial )\sim {{\rm{N}}_{{{y}}|{{x}}}}({{Dx}}, \partial {{I}}) (21) 似然函数的概率密度函数形式确定后，未知参数为期望和方差，需要用样本对参数进行估计。
假定待恢复稀疏矩阵
{{s}} 中的每行之间独立，且服从正态分布\begin{array}{*{20}{c}} {p({{{s}}_i};{\gamma _i}, {{{B}}_i})\sim{\rm{N}}({{0}}, {\gamma _i}{{{B}}_i}), }&{i = 1, 2, ·\!·\!· , L} \end{array} (22) 其中，
{\gamma _i} 控制{{s}} 的行稀疏性，{{{B}}_i} 为正定矩阵，用于描述{{{s}}_i} 的相关结构。则
{{x}} 的先验概率满足p({{x}};{\gamma _i}, {{{B}}_i}, \forall i)\sim{{\rm{N}}_{{x}}}({{0}}, {{{Σ}} _0}) (23) 其中，
{{{Σ}} _0} 表示为{{{Σ}} _0} = \left[ {\begin{array}{*{20}{c}} {{\gamma _1}{{{B}}_1}}&{}&{}&{} \\ {}&{{\gamma _2}{{{B}}_2}}&{}&{} \\ {}&{}& \ddots &{} \\ {}&{}&{}&{{\gamma _L}{{{B}}_L}} \end{array}} \right] (24) 为了防止过拟合，
L 个块采用同一个{{B}} 描述相关结构[13]，得{{{Σ}} _0} = {{Γ}} \otimes {{B}} (25) 其中，
{{Γ}} = {\rm{diag}}({\gamma _1}, {\gamma _2}, ·\!·\!· , {\gamma _L}) ，根据贝叶斯公式，后验概率p({{x}}|{{y}}) 由先验概率和似然函数表示为p({{x}}|{{y}}) = \frac{{p({{y}}|{{x}}){{p}}({{x}})}}{{p({{y}})}} 利用贝叶斯准则得后验概率分布
p({{x}}|{{y}};\partial , {\gamma _i}, {{B}}, \forall i)\sim{{\rm{N}}_{{{x}}|{{y}}}}({{{μ}} _{{x}}}, {{{Σ}} _{{x}}}) (26) 其中，
\begin{aligned} {{{Σ}} _{{x}}} =& {\left( {{{Σ}} _0^{ - 1} + \frac{1}{\partial }{{{D}}^{\rm{T}}}{{D}}} \right)^{ - 1}}\\ =& {{{Σ}} _0} - {{{Σ}} _0}{{{D}}^{\rm{T}}}{\left( {\partial {{I}} + {{D}}{{{Σ}} _0}{{{D}}^{\rm{T}}}} \right)^{ - 1}}{{D}}{{{Σ}} _0} \end{aligned} (27) $ {{{μ}} _{{x}}} = \frac{1}{\partial }{{{Σ}} _{{x}}}{{{D}}^{\rm{T}}}{{y}} \hspace{115pt} (28) {{x}} 的MAP(Maximum-A-Posterior)估计为\begin{aligned} {{{x}}^*} \triangleq {{{μ}} _x} =& {\left( {\partial {{{Σ}} _0}^{ - 1} + {{{D}}^{\rm{T}}}{{D}}} \right)^{ - 1}}{{{D}}^{\rm{T}}}{{y}} \\ =& {{{Σ}} _0}{{{D}}^{\rm{T}}}{\left( {\partial {{I}} + {{D}}{{{Σ}} _0}{{{D}}^{\rm{T}}}} \right)^{ - 1}}{{y}} \end{aligned} (29) 由此可知，对参数集
{{Θ}} = \left\{ {{\gamma _1}, {\gamma _2}, ·\!·\!· , {\gamma _L}, {{B}}, \partial } \right\} 进行估计，即求得稀疏解。最大化
p({{y}};{{Θ}} ) 可得参数估计，相当于最小化- \lg p({{y}};{{Θ}} ) ，目标函数为L\left( {{Θ}} \right) = {{{y}}^{\rm{T}}}{{Σ}} _{{y}}^{ - 1}{{y}} + \lg \left| {{{{Σ}} _{{y}}}} \right| (30) 通过最大期望(Expectation-Maximization, EM)算法，得到参数更新[13]
\!\!\!\begin{array}{*{20}{c}} {{\gamma _i}\! =\! \frac{{{{\rm{t}}{\rm{r}}}\left[ {{{{B}}^{ - 1}}({{Σ}} _x^i \!+\! {{μ}} _x^i{{({{μ}} _x^i)}^{\rm{T}}})} \right]}}{N}, }&{i \!=\! 1, 2, ·\!·\!· , L} \end{array} (31) 其中，
{{\rm{tr}}}\left[ \bullet \right] 表示矩阵求迹，{{{μ}} _x^i}{\underline{\underline \Delta } } {{{μ}} _x}\left(\! {\left( {i \!-\! 1} \right)\!N \!+\! 1\!:iN} \right) ,{{{Σ}} _x^i}{\underline{\underline \Delta } } {{{Σ}} _x}\left( {\left( {i - 1} \right)N + 1:iN, \left( {i - 1} \right)N + 1:iN} \right) 。{{B}} = \frac{1}{L}\sum\limits_{i = 1}^L {\frac{{{{Σ}} _x^i + {{μ}} _x^i{{({{μ}} _x^i)}^{\rm{T}}}}}{{{\gamma _i}}}} (32) 为了增加鲁棒性，按以下规则对
{{B}} 进行估计{{B}} = {{\mathop {{B}}\limits^ \approx } /{||\mathop {{B}}\limits^ \approx |{|_{\rm{F}}}}} (33) 其中，
\mathop {{B}}\limits^ \approx = \displaystyle\sum\nolimits_{i = 1}^L {\frac{{{{μ}} _x^i{{({{μ}} _x^i)}^{\rm{T}}}}}{{{\gamma _i}}}} + h{{I}} ,h 为一正的常量，||\mathop {{B}}\limits^ \approx |{|_{\rm{F}}} 表示\mathop {{B}}\limits^ \approx 的Frobenius范数。\partial 的迭代公式为\partial = \frac{{\left\| {{{y}} - {{D}}{{{μ}} _x}} \right\|_2^2 + \partial \left[ {LN - {{\rm{tr}}}({{{Σ}} _x}{{Σ}} _0^{ - 1})} \right]}}{{MN}} (34) 当达到设定的迭代次数小于门限时，停止迭代，此时，
{{\gamma _i}, }i = 1, 2, ·\!·\!· , L 中最大值对应目标方位。4. 仿真与实测分析
4.1 仿真实验
仿真条件：仿真采用正交极化双通道阵列天线接收信号，基于7阵元均匀线阵，阵元间距取波长的一半。载波波段为C波段，信号形式为线性调频信号，带宽为4 MHz，采样频率为5 MHz。
(1)为验证干扰信号和雷达信号的极化角度差以及迭代次数对本文方法测角性能的影响，下面对本文方法测角误差进行仿真分析。目标雷达数目为1，方位角为10°，极化矢量固定为
{{{h}}_1} \!=\!\! \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}0^\circ)}}} \end{array}} \!\!\!\!\right] 。干扰数目为1，方位角为–10°，极化矢量为{{h}}_1^1 = \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}\Delta \phi )}}} \end{array}} \!\!\!\!\right] ，其中，\Delta \phi 的变化间隔为0.5°，变化范围为[12.0° 24.5°]。接收信号前面采样数据信干比为–3 dB，用于干扰极化状态估计，后面采样数据信干比为0 dB。固定BSBL算法迭代次数为100次，本文方法的测角误差随目标和干扰的极化角度差变化曲线如图1(a)所示。然后，固定干扰极化矢量{{h}}_1^1 = \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}40^\circ)}}} \end{array}} \!\!\!\!\right] ，迭代次数的变化间隔为2次，变化范围为[108次 160次]，本文方法的测角误差随迭代次数变化的曲线如图1(b)所示。由图1(a)可知，信噪比越高，测角误差越小，同一信噪比下，目标和干扰极化角度差越大，测角误差越小，并逐步趋于稳定；由图1(b)可知，信噪比越高，对迭代次数的要求越小，同一信噪比下，迭代次数越大，测角误差越小，并逐步趋于稳定。
(2)为验证有源欺骗干扰环境和小样本条件下所提方法DOA估计的性能，将本文方法与BSBL算法、APF联合子空间类算法、APF联合最大似然(Maximum Likelihood, ML)算法进行对比。信噪比取为20 dB，样本数为5，干扰数目为1，方位为10°，目标数目为2，方位分别为18°和25°。干扰极化矢量为
{{h}}_1^1 = \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}90^\circ)}}} \end{array}} \!\!\!\!\right] ，目标1极化矢量为{{{h}}_1} = \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}0^\circ)}}} \end{array}} \!\!\!\!\right] ，目标2极化矢量为{{{h}}_2} = \left[\!\!\!\! {\begin{array}{*{20}{c}} {\cos(45^\circ)}\;{\sin(45^\circ){{\rm{e}}^{({\rm{j}}20^\circ)}}} \end{array}} \!\!\!\!\right] 。不同方法DOA估计结果如图2所示。由图2可知，BSBL算法的估计结果偏离正确的方位角，干扰信号影响BSBL算法的稀疏恢复性能，这是因为待测源个数增多要求信噪比和BSBL算法迭代次数增加，导致BSBL算法测角性能下降；通过APF联合子空间类算法进行DOA估计，由于阵元数目较少、样本数少，子空间类算法空间分辨率有限，不可以同时对两个目标源进行准确的DOA估计；APF联合最大似然算法也分辨率不足，无法对两个目标进行准确的DOA估计；本文方法先通过APF滤除干扰，再通过BSBL进行DOA估计，可以分辨两个目标源，测得两个目标源方位角分别为18°和26°。综上所述，在有源欺骗干扰环境、阵元数少和小样本条件下，本文方法具有较高的空间分辨率和测角精度。
4.2 实测实验
实测场景和条件：实验采用基于弹载的正交极化双通道阵列天线接收信号，基于7阵元非均匀圆阵。为有效掩护雷达信号，干扰辐射功率高于雷达平均副瓣1～3 dB，并低于雷达主瓣约20～25 dB。各干扰源的辐射强度基本相同，天线指向固定，因此，弹载雷达收到的干扰信号功率相对稳定。典型配置为3个干扰，为了更好的干扰弹载雷达系统，使其不能准确测量雷达信号源方向，一般沿雷达两侧及径向布置有源欺骗干扰，布设距离100～300 m。雷达主瓣极化形式多为线极化，而雷达副瓣的极化形式较为杂乱，但此时雷达功率低于干扰功率较多。干扰场景如图3所示。
为在时-频-空域保护目标，干扰信号脉冲覆盖雷达信号脉冲的前后沿。干扰信号与雷达信号的载频、重频等参数相同或接近，3个干扰分布在雷达主瓣内，时-频-空域抑制干扰手段失效。导弹对雷达进行初始探测时，接收到的信号来自雷达副瓣和干扰。所以，初始探测采样数据可以近似为干扰数据，用于极化状态估计，进而在极化域实现抗干扰。在导弹攻击过程中，初始阶段干扰与弹载雷达的空间角度差较小，随着弹目不断接近，空间角度差逐渐变大。参考方向如图4所示。
因为仰角和方位角的2维测量导致冗余字典长度增加，所以利用本文方法进行DOA估计时，先进行大范围粗略估计，再进行小范围精确估计。仿真实验中用到的ESPRIT算法利用自相关矩阵信号子空间旋转不变性，从子阵固定关系中求得来波方向，在测角前需要先将总阵分为几个布局一致的子阵。由此可知，ESPRIT算法不适用于基于非均匀圆阵的DOA估计。所以实测数据处理时仅将所提方法与APF联合MUSIC算法和APF联合ML算法进行对比。为了定量评估DOA估计精度，测角误差表示为方位和仰角测角误差的均方根。基于接收到的20个脉冲信号，3种方法的测角均方根误差如图5所示。
由图5可知，本文方法同样适用于干扰环境下基于非均匀圆阵的2维DOA估计。在有源欺骗干扰和小样本环境下，相比APF联合MUSIC算法或最大似然算法，本文方法具有更高的测角精度。
5. 结束语
有源欺骗干扰的存在使得待测源个数增多，影响BSBL算法测向性能。因此，本文提出APF联合BSBL算法，用于有源欺骗干扰环境下基于小样本的DOA估计。本文方法先基于接收数据估计干扰极化状态，进而在极化域抑制干扰信号。然后，通过BSBL算法估计目标方向，有效提高了有源欺骗干扰环境下BSBL算法的测向性能。仿真和实测实验处理结果表明本文方法在有源欺骗干扰环境下具有较高的测角精度和空间分辨率。如何进一步降低本文方法运算量，使其可以应用于更为复杂的战场环境，将是今后的研究重点。
-
ZHOU Chao, LIU Quanhua, and CHEN Xinliang. Parameter estimation and suppression for DRFM-based interrupted sampling repeater jammer[J]. IET Radar, Sonar & Navigation, 2018, 12(1): 56–63. doi: 10.1049/iet-rsn.2017.0114 LI Zheng, TAI Ning, WANG Chao, et al. A study on blanket noise jamming to LFM pulse compression radar[C]. 2017 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Xiamen, China, 2017: 1–5. doi: 10.1109/ICSPCC.2017.8242411. WU Xiaohuan, ZHU Weiping, and YAN Jun. A fast gridless covariance matrix reconstruction method for one- and two-dimensional direction-of-arrival estimation[J]. IEEE Sensors Journal, 2017, 17(15): 4916–4927. doi: 10.1109/JSEN.2017.2709329 贾伟娜, 刘顺兰. 模拟退火遗传算法在DOA估计技术中的应用[J]. 计算机工程与应用, 2014, 50(12): 266–270. doi: 10.3778/j.issn.1002-8331.1206-0247JIA Weina and LIU Shunlan. Application of simulated annealing genetic algorithm in DOA estimation technique[J]. Computer Engineering and Applications, 2014, 50(12): 266–270. doi: 10.3778/j.issn.1002-8331.1206-0247 ZHANG T T, LU Y L, and HUI H T. Compensation for the mutual coupling effect in uniform circular arrays for 2D DOA estimations employing the maximum likelihood technique[J]. IEEE Transactions on Aerospace and Electronic Systems, 2008, 44(3): 1215–1221. doi: 10.1109/TAES.2008.4655375 WAN Liangtian, HAN Guangjie, JIANG Jinfang, et al. DOA estimation for coherently distributed sources considering circular and noncircular signals in massive MIMO systems[J]. IEEE Systems Journal, 2017, 11(1): 41–49. doi: 10.1109/JSYST.2015.2445052 蔡晶晶, 宗汝, 蔡辉. 基于空域平滑稀疏重构的DOA估计算法[J]. 电子与信息学报, 2016, 38(1): 168–173. doi: 10.11999/JEIT150538CAI Jingjing, ZONG Ru, and CAI Hui. DOA estimation via sparse representation of the smoothed array covariance matrix[J]. Journal of Electronics &Information Technology, 2016, 38(1): 168–173. doi: 10.11999/JEIT150538 AL-SHOUKAIRI M, SCHNITER P, and RAO B D. A GAMP-based low complexity sparse Bayesian learning algorithm[J]. IEEE Transactions on Signal Processing, 2018, 66(2): 294–308. doi: 10.1109/TSP.2017.2764855 HOU Huijun and MAO Xingpeng. Oblique projection and sparse reconstruction based DOA estimation of hybrid completely and partially polarized signals with arbitrary polarimetric arrays[C]. Proceedings of the 30th Canadian Conference on Electrical and Computer Engineering (CCECE), Windsor, Canada, 2017: 1–4. doi: 10.1109/CCECE.2017.7946799. ZHEN Jiaqi and WANG Zhifang. DOA estimation method for wideband signals by sparse recovery in frequency domain[J]. Journal of Systems Engineering and Electronics, 2017, 28(5): 871–878. doi: 10.21629/JSEE.2017.05.06 WANG Yi, CHEN Baixiao, ZHENG Yisong, et al. Joint power distribution and direction of arrival estimation for wideband signals using sparse Bayesian learning[J]. IET Radar, Sonar & Navigation, 2017, 11(1): 52–59. doi: 10.1049/iet-rsn.2015.0610 王洪雁, 房云飞, 裴炳南. 基于矩阵补全的二阶统计量重构DOA估计方法[J]. 电子与信息学报, 2018, 40(6): 1383–1389. doi: 10.11999/JEIT170826WANG Hongyan, FANG Yunfei, and PEI Bingnan. Matrix completion based second order statistic reconstruction DOA estimation method[J]. Journal of Electronics &Information Technology, 2018, 40(6): 1383–1389. doi: 10.11999/JEIT170826 ZHANG Zhilin and RAO B D. Sparse signal recovery with temporally correlated source vectors using sparse Bayesian learning[J]. IEEE Journal of Selected Topics in Signal Processing, 2011, 5(5): 912–926. doi: 10.1109/JSTSP.2011.2159773 WANG Lu, ZHAO Lifan, RAHARDJA S, et al. Alternative to extended block sparse Bayesian learning and its relation to pattern-coupled sparse Bayesian learning[J]. IEEE Transactions on Signal Processing, 2018, 66(10): 2759–2771. doi: 10.1109/TSP.2018.2816574 HUANG Qinghua, ZHANG Guangfei, and FANG Yong. DOA estimation using block variational sparse Bayesian learning[J]. Chinese Journal of Electronics, 2017, 26(4): 768–772. doi: 10.1049/cje.2017.04.004 宫健, 楼顺天, 张伟涛. 一种强干扰条件下阵列天线波达方向估计方法[J]. 西安电子科技大学学报: 自然科学版, 2018, 45(1): 168–172. doi: 10.3969/j.issn.1001-2400.2018.01.030GONG Jian, LOU Shuntian, and ZAHNG Weitao. Method of array antenna DOA under strong interference presence[J]. Journal of Xidian University, 2018, 45(1): 168–172. doi: 10.3969/j.issn.1001-2400.2018.01.030 MA Jiazhi, SHI Longfei, LI Yongzhen, et al. Angle estimation of extended targets in main-lobe interference with polarization filtering[J]. IEEE Transactions on Aerospace and Electronic Systems, 2017, 53(1): 169–189. doi: 10.1109/TAES.2017.2649783 NATHANSON F E. Adaptive circular polarization[C]. IEEE International Radar Conference, Arlington, USA, 1975: 221–225. 王雪松, 汪连栋, 肖顺平, 等. 自适应极化滤波器的理论性能分析[J]. 电子学报, 2004, 32(8): 1326–1329. doi: 10.3321/j.issn:0372-2112.2004.08.023WANG Xuesong, WANG Liandong, XIAO Shunping, et al. Theoretical performance analysis of adaptive polarization filters[J]. Acta Electronica Sinica, 2004, 32(8): 1326–1329. doi: 10.3321/j.issn:0372-2112.2004.08.023 任博, 罗笑冰, 邓方刚, 等. 应用极化聚类中心设计快速自适应极化滤波器[J]. 国防科技大学学报, 2015, 37(4): 87–92. doi: 10.11887/j.cn.201504015REN Bo, LUO Xiaobing, DENG Fanggang, et al. Design of fast adaptive polarization filters utilizing polarizing cluster center[J]. Journal of National University of Defense Technology, 2015, 37(4): 87–92. doi: 10.11887/j.cn.201504015 期刊类型引用(7)
1. 赵书红，董绍武，白杉杉，高喆. 一种优化的频率驾驭算法研究. 电子与信息学报. 2021(05): 1457-1464 . 
本站查看2. 陶贵丽，刘文强，张兴华，牛晓霞. 带丢包不确定广义系统鲁棒Kalman预报器. 系统科学与数学. 2021(05): 1215-1232 . 
百度学术3. 谢卫，王前东. 一种基于自适应网格剖分的协方差交集融合新算法. 电讯技术. 2019(09): 1067-1074 . 百度学术4. 刘振亚，高敏，程呈. 基于理想弹道鲁棒容积卡尔曼滤波视线角估计. 系统工程与电子技术. 2018(02): 409-416 . 百度学术5. 秦文利，胡捍英，陈松. 基于带势概率假设密度粒子滤波的MIMO雷达检测前跟踪算法. 信息工程大学学报. 2018(02): 140-145 . 百度学术6. 王雪梅，刘文强，邓自立. 带丢失观测和不确定噪声方差系统改进的鲁棒协方差交叉融合稳态Kalman滤波器. 控制理论与应用. 2016(07): 973-979 . 百度学术7. 王雪梅，刘文强，邓自立. 带不确定协方差线性相关白噪声系统改进的鲁棒协方差交叉融合稳态Kalman估值器. 控制与决策. 2016(10): 1749-1756 . 百度学术其他类型引用(5)
-
下载:
下载:
百度学术
计量
- 文章访问数: 2582
- HTML全文浏览量: 773
- PDF下载量: 135
- 被引次数: 12
下载:
