王海洋,刘光斌,范志良,姚志成
(火箭军工程大学 导弹工程学院, 西安 710025)
摘要:
为解决干扰信号来向快速变化情况下,基于线性约束最小功率(linear constraint minimum power,LCMP)准则的传统抗干扰算法可能难以收敛,致使全球卫星导航系统(global navigation satellite system,GNSS)接收机失锁的问题,提出了一种基于协方差矩阵锥化(covariance matrix taper,CMT)的零陷加宽抗干扰算法,使变化的干扰信号来向能够始终保持在拓展的零陷之内.首先,在实际干扰信号来向变化不可预测的条件下,所提算法假设干扰信号来向变化服从三角形概率分布统计模型；然后,依据这一模型推导假设情况下的虚拟协方差矩阵与真实协方差矩阵的转换关系,再根据转换关系和真实采样协方差矩阵得到虚拟采样协方差矩阵.通过虚拟采样协方差矩阵求解阵列天线各通道权值；最后,将所提三角形分布算法与传统的均匀分布算法和拉普拉斯分布算法进行了对比分析.研究结果表明, 该算法可以有效加宽在干扰信号来向上的零陷并抑制干扰,而且相比于均匀分布算法和拉普拉斯分布算法,三角形分布算法在形成零陷的宽度和深度一致性之间更加平衡,阵列输出信干噪比也要更高,最后结合软件接收机验证了所提算法的有效性.
关键词: GNSS接收机 阵列天线 零陷加宽 协方差矩阵锥化 三角形分布
DOI：10.11918/j.issn.0367-6234.201806097
分类号:TN973.3+1
文献标识码:A
基金项目:国家自然科学基金(1,3)
A null widening anti-jamming algorithm for GNSS receivers
WANG Haiyang,LIU Guangbin,FAN Zhiliang,YAO Zhicheng
(School of Missile Engineering, Rocket Force University of Engineering, Xian 710025, China)
Abstract:
Since the directions of arrival (DOAs) of interferences change rapidly in high dynamic conditions, the traditional anti-jamming algorithm based on the linearly constrained minimum power (LCMP) criterion cannot suppress the interferences effectively, which may lead to the acquisition failure of GNSS receivers. To slove this problem, a covariance matrix taper (CMT) null widening algorithm was proposed, which can make the DOAs of interferences always in the widening nulls. Firstly, the DOAs of interferences which are unpredictable in practice were assumed to be triangular distributed. Then, the conversion relationship between the virtual covariance matrix and the real covariance matrix was obtained according to the assumption. Consequently, the virtual sample covariance matrix was derived, based on which the weight vector was acquired. Finally, the propsed algorithm was compared with the algorithms that adopt uniform distribution model and Laplace distribution model. The research results show that the proposed algorithm can widen the nulls effectively. Compared with uniform distribution model and Laplace distribution model, the proposed model achieved a superior balance between the null width and the null depth, and it had a higher array output signal to interference plus noise ratio (SINR). Moreover, software GNSS receiver was used to validate the effectiveness of the proposed algorithm.
Key words: GNSS receiver array antenna null widening covariance matrix taper triangular distribution
王海洋, 刘光斌, 范志良, 姚志成. 一种针对GNSS接收机的宽零陷抗干扰算法[J]. 哈尔滨工业大学学报, 2019, 51(4): 94-98,137. DOI: 10.11918/j.issn.0367-6234.201806097.
WANG Haiyang, LIU Guangbin, FAN Zhiliang, YAO Zhicheng. A null widening anti-jamming algorithm for GNSS receivers[J]. Journal of Harbin Institute of Technology, 2019, 51(4): 94-98,137. DOI: 10.11918/j.issn.0367-6234.201806097.
基金项目 国家自然科学基金(61501471, 61503393) 作者简介 王海洋(1992—)，男，博士研究生;
刘光斌(1963—)，男，教授，博士生导师 通信作者 王海洋，haiyang_wxdh@163.com 文章历史 收稿日期: 2018-06-14
Contents -->Abstract Full text Figures/Tables PDF
一种针对GNSS接收机的宽零陷抗干扰算法
王海洋
, 刘光斌, 范志良, 姚志成 火箭军工程大学 导弹工程学院, 西安 710025
收稿日期: 2018-06-14
基金项目: 国家自然科学基金(61501471, 61503393)
作者简介: 王海洋(1992—)，男，博士研究生;
刘光斌(1963—)，男，教授，博士生导师
通信作者: 王海洋，haiyang_wxdh@163.com
摘要: 为解决干扰信号来向快速变化情况下，基于线性约束最小功率(linear constraint minimum power, LCMP)准则的传统抗干扰算法可能难以收敛，致使全球卫星导航系统(global navigation satellite system, GNSS)接收机失锁的问题，提出了一种基于协方差矩阵锥化(covariance matrix taper, CMT)的零陷加宽抗干扰算法，使变化的干扰信号来向能够始终保持在拓展的零陷之内.首先，在实际干扰信号来向变化不可预测的条件下，所提算法假设干扰信号来向变化服从三角形概率分布统计模型；然后，依据这一模型推导假设情况下的虚拟协方差矩阵与真实协方差矩阵的转换关系，再根据转换关系和真实采样协方差矩阵得到虚拟采样协方差矩阵.通过虚拟采样协方差矩阵求解阵列天线各通道权值；最后，将所提三角形分布算法与传统的均匀分布算法和拉普拉斯分布算法进行了对比分析.研究结果表明, 该算法可以有效加宽在干扰信号来向上的零陷并抑制干扰，而且相比于均匀分布算法和拉普拉斯分布算法，三角形分布算法在形成零陷的宽度和深度一致性之间更加平衡，阵列输出信干噪比也要更高，最后结合软件接收机验证了所提算法的有效性.
关键词: GNSS接收机 阵列天线 零陷加宽 协方差矩阵锥化 三角形分布
A null widening anti-jamming algorithm for GNSS receivers
WANG Haiyang
, LIU Guangbin, FAN Zhiliang, YAO Zhicheng School of Missile Engineering, Rocket Force University of Engineering, Xi'an 710025, China
Abstract: Since the directions of arrival (DOAs) of interferences change rapidly in high dynamic conditions, the traditional anti-jamming algorithm based on the linearly constrained minimum power (LCMP) criterion cannot suppress the interferences effectively, which may lead to the acquisition failure of GNSS receivers. To slove this problem, a covariance matrix taper (CMT) null widening algorithm was proposed, which can make the DOAs of interferences always in the widening nulls. Firstly, the DOAs of interferences which are unpredictable in practice were assumed to be triangular distributed. Then, the conversion relationship between the virtual covariance matrix and the real covariance matrix was obtained according to the assumption. Consequently, the virtual sample covariance matrix was derived, based on which the weight vector was acquired. Finally, the propsed algorithm was compared with the algorithms that adopt uniform distribution model and Laplace distribution model. The research results show that the proposed algorithm can widen the nulls effectively. Compared with uniform distribution model and Laplace distribution model, the proposed model achieved a superior balance between the null width and the null depth, and it had a higher array output signal to interference plus noise ratio (SINR). Moreover, software GNSS receiver was used to validate the effectiveness of the proposed algorithm.
Keywords: GNSS receiver array antenna null widening covariance matrix taper triangular distribution
全球卫星导航系统因其导航精度高和不随时间积累误差的优点得到了广泛的应用，然而卫星导航信号到达GNSS接收机时的功率比环境基底噪声还要低20 dB，所以极易受到强干扰信号压制.此时GNSS接收机将无法正常捕获卫星导航信号，必然导致接收机搜星定位失败[1-2].因此目前普遍在接收机端采用阵列天线抑制强干扰信号.工程中常采用线性约束最小功率准则在强干扰信号来向上自适应生成零陷来抑制干扰[3]，并采用最小均方算法迭代求解阵列权值.但是在GNSS接收机载体处于高频振动或者高速运动的状态下，干扰信号来向在短时间内会发生快速变化，上述算法可能无法及时收敛，致使干扰信号不能被及时有效地抑制[4-6]，导致GNSS接收机失锁.
针对这一问题，通过采用基于协方差矩阵锥化的零陷加宽算法可以有效加宽干扰零陷，使得干扰信号变化方向始终保持在零陷之内，从而达到抑制干扰的目的.Mailloux[7]首先提出了协方差矩阵锥化零陷加宽算法，并假设在干扰源附近均匀分布着若干相同干扰源，求解假设情况下信号协方差矩阵与真实情况下协方差矩阵的关系，再利用真实协方差矩阵得到虚拟协方差矩阵，即可实现干扰零陷拓宽；Zatman[8]则假设干扰源带宽为真实带宽的若干倍，通过求解虚拟带宽情况下信号协方差矩阵与真实带宽下协方差矩阵的关系来拓宽干扰零陷.Guerci等[9-10]又系统总结了协方差矩阵锥化算法的一般表达式.在此基础上，李荣锋等[11]提出了新的零陷加宽抗干扰思路，即假设干扰信号来向变化服从高斯概率分布，实现了干扰零陷拓宽，并且证明了当干扰信号来向变化服从均匀分布时其算法与Zatman算法是等价的；随后，Lu等[12]提出了干扰信号来向变化服从拉普拉斯分布的零陷加宽方法；武思军等[13]提出了干扰信号来向变化服从伯努利分布的零陷加宽方法；张柏华等[14-15]研究了基于空时阵列处理的零陷加宽算法；王妙等[16]提出了基于降秩共轭梯度法的零陷加宽技术，能够有效降低计算量；柯熙政等[17]推导了基于均匀圆阵的零陷加宽算法.
需要说明的是，干扰信号来向变化实际上是未知的，并非一定是均匀分布，拉普拉斯分布，抑或是高斯分布.文献[11-13]中借用概率分布统计模型的目的是为了拓展干扰零陷宽度，与Mailloux虚拟干扰源的本质是相同的，只是采用不同概率分布统计模型所形成的零陷形状有所区别，性能也会有差异.正因为干扰来向变化是未知的，干扰零陷深度在零陷宽度范围内更应该保持一致.
本文提出了一种假设干扰信号来向变化服从三角形分布概率统计模型，基于线性约束最小功率准则的协方差矩阵锥化零陷加宽抗干扰算法，并将其与具有代表性的均匀分布、拉普拉斯分布算法进行比较分析，理论分析与仿真实验证明了所提算法拓宽的零陷更加平稳，输出信干噪比也更高，最后结合BD2软件接收机验证了该算法的有效性.
1 信号模型GNSS接收机阵列天线模型如图 1所示，阵列天线包含M个阵元，阵元之间间距为d，每个阵元对应一个独立的处理通道.远场信号源波面到达阵列天线前端可视为平面波(plane wave)，各阵元接收信号后首先在射频前端(radio frequency front end, RFFE)中对其进行放大去噪，再经模数转换器(analog to digital converter, ADC)抽样提取为数字信号，最后对每一路信号进行加权求和得到输出值.
Fig. 1
图 1 阵列天线模型 Fig. 1 Array antenna model 假设L个卫星导航信号和Q个窄带干扰信号从远场入射，那么阵列天线的接收信号模型可以表示为
$\mathit{\boldsymbol{x}}\left( t \right) = \sum\limits_{l = 1}^L {\mathit{\boldsymbol{a}}\left( {{\theta _l}} \right){s_l}\left( t \right)} + \sum\limits_{q = 1}^Q {\mathit{\boldsymbol{a}}\left( {{\theta _q}} \right){j_q}\left( t \right)} + \mathit{\boldsymbol{e}}\left( t \right).$
式中：θl为第l个卫星导航信号来向; sl(t)为第l个卫星导航信号包络；θq为第q个干扰信号来向; jq(t)为第q个干扰信号包络；e(t)为噪声向量；a(θ)为信号空域导向矢量，对于半波长均匀线阵，其表达式为[1, e-jπdsin θ, …, e-jπ(M-1)dsin θ].
线性约束最小功率准则是在线性约束条件下，通过使阵列输出功率最小来求解权矢量.工程中常设置约束条件为天线参考阵元无失真，其目标函数为
$\mathop {\min }\limits_w {\mathit{\boldsymbol{w}}^{\rm{H}}}{\mathit{\boldsymbol{R}}_x}\mathit{\boldsymbol{w}},{\rm{s}}.\;{\rm{t}}.\;{\mathit{\boldsymbol{w}}^{\rm{H}}}{\mathit{\boldsymbol{c}}_M} = 1.$ (1)
式中：上标“H”为Hermit转置；cM=[1, 0, …, 0]H为M维约束列向量；Rx为仅包含干扰信号与噪声的理论协方差矩阵.该算法也称功率倒置(power inversion, PI)算法，本文所提算法将基于该准则求解阵列权值.
利用拉格朗日乘数法，设置式(1)的最小化函数为
$J = \frac{1}{2}{\mathit{\boldsymbol{w}}^{\rm{H}}}{\mathit{\boldsymbol{R}}_x}\mathit{\boldsymbol{w}} + \lambda \left( {{\mathit{\boldsymbol{w}}^{\rm{H}}}{\mathit{\boldsymbol{c}}_M} - 1} \right),$ (2)
式中λ为拉格朗日乘数.为了求解式(2)的最优化问题，求解其关于w的梯度，并令其为0，可得
${\nabla _w}J\left( \lambda \right) = {\mathit{\boldsymbol{R}}_x}\mathit{\boldsymbol{w}} + \lambda {\mathit{\boldsymbol{c}}_M} = 0,$
进而可知
$\mathit{\boldsymbol{w}} = - \lambda \mathit{\boldsymbol{R}}_x^{ - 1}{\mathit{\boldsymbol{c}}_M}.$ (3)
为了求解拉格朗日参数λ，将式(3)代入式(1)中的约束条件，可得
${\left( { - \lambda \mathit{\boldsymbol{R}}_x^{ - 1}{\mathit{\boldsymbol{c}}_M}} \right)^{\rm{H}}}{\mathit{\boldsymbol{c}}_M} = 1,$
则拉格朗日参数λ表示如下:
$\lambda = - \frac{1}{{\mathit{\boldsymbol{c}}_M^{\rm{H}}\mathit{\boldsymbol{R}}_x^{ - 1}{\mathit{\boldsymbol{c}}_M}}}.$ (4)
将式(4)代入式(3)，即可得到权矢量的最优解wopt为
${\mathit{\boldsymbol{w}}_{{\rm{opt}}}} = \frac{{\mathit{\boldsymbol{R}}_x^{ - 1}{\mathit{\boldsymbol{c}}_M}}}{{\mathit{\boldsymbol{c}}_M^{\rm{H}}\mathit{\boldsymbol{R}}_x^{ - 1}{\mathit{\boldsymbol{c}}_M}}}.$ (5)
由于实际中难以得到理论协方差矩阵，一般利用快拍得到的样本协方差矩阵
${{\mathit{\boldsymbol{\tilde R}}}_x} = \frac{1}{N}\sum\limits_{n = 1}^N {\mathit{\boldsymbol{x}}\left( n \right){\mathit{\boldsymbol{x}}^{\rm{H}}}\left( n \right)} .$
式中：N为采样得到的快拍数；x(n)为阵列天线实际接收到的包含目标信号、干扰信号与噪声的采样序列.
那么最终得到基于功率倒置算法的阵列权矢量
$\mathit{\boldsymbol{\tilde w}} = \frac{{\mathit{\boldsymbol{\tilde R}}_x^{ - 1}{\delta _M}}}{{\delta _M^{\rm{H}}\mathit{\boldsymbol{\tilde R}}_x^{ - 1}{\delta _M}}}.$ (6)
2 Triangular-CMT-PI算法基于信号模型和阵列天线模型，假设干扰信号来向发生变化，即
$\tilde \theta = {\theta _0} + \Delta \theta .$
式中：θ0为干扰信号初始来向；Δθ为干扰信号来向变化幅度；
假设第q个干扰信号来向变化幅度Δθq在区间[-β, β]内符合三角形分布，那么其概率分布函数如下
$f\left( {\Delta {\theta _q}} \right) = \left\{ \begin{array}{l}\frac{{\Delta {\theta _q} + \beta }}{{{\beta ^2}}}, - \beta \le \Delta {\theta _q} \le 0;\\\frac{{\beta - \Delta {\theta _q}}}{{{\beta ^2}}},0 < \Delta {\theta _q} \le \beta .\end{array} \right.$
阵列天线接收的信号包括卫星导航信号、干扰信号和噪声，由于卫星导航信号的功率很弱，一般比噪声还要低20 dB以上，在计算协方差矩阵时可以将其忽略，所以有
$\mathit{\boldsymbol{\tilde R}}\left( x \right) \approx {{\mathit{\boldsymbol{\tilde R}}}_j} + {{\mathit{\boldsymbol{\tilde R}}}_{\rm{e}}} = \sum\limits_{q = 1}^Q {\sigma _q^2\mathit{\boldsymbol{a}}\left( {{\theta _q}} \right){\mathit{\boldsymbol{a}}^{\rm{H}}}\left( {{\theta _q}} \right)} + \sigma _{\rm{e}}^2\mathit{\boldsymbol{I}}.$
式中：σq2为第q个干扰信号的功率；a(θq)为第q个干扰信号的导向矢量；σe2为噪声信号的功率；I为单位矩阵.当干扰信号来向发生变化时，式(5)应该变化如下
${{\mathit{\boldsymbol{\tilde R}}}_L}\left( x \right) \approx \sum\limits_{q = 1}^Q {\sigma _q^2\int {f\left( {\Delta {\theta _q}} \right)\mathit{\boldsymbol{a}}\left( {{{\tilde \theta }_q}} \right){\mathit{\boldsymbol{a}}^{\rm{H}}}\left( {{{\tilde \theta }_q}} \right){\rm{d}}\Delta \theta } } + \sigma _{\rm{e}}^2\mathit{\boldsymbol{I}},$
则可得
$\begin{array}{l}{{\mathit{\boldsymbol{\tilde R}}}_L}\left( {m,n} \right) \approx \\\sigma _{\rm{e}}^2{\delta _{mn}} + \sum\limits_{q = 1}^Q {\sigma _q^2\int_{ - \beta }^\beta {f\left( {\Delta {\theta _q}} \right)} } \cdot \\\exp \left( { - j{\rm{ \mathsf{ π} }}\left( {m - n} \right)\sin {{\tilde \theta }_q}} \right){\rm{d}}\Delta {\theta _q} = \\\sigma _{\rm{e}}^2{\delta _{mn}} + \sum\limits_{q = 1}^Q {\sigma _q^2\int_{ - \beta }^\beta {f\left( {\Delta {\theta _q}} \right)} } \cdot \\\exp \left( { - j{\rm{ \mathsf{ π} }}\left( {m - n} \right)\sin \left( {{\theta _q} + \Delta {\theta _q}} \right)} \right){\rm{d}}\Delta {\theta _q} \approx \\\sigma _{\rm{e}}^2{\delta _{mn}} + \sum\limits_{q = 1}^Q {\sigma _q^2\exp \left( { - j{\rm{ \mathsf{ π} }}\left( {m - n} \right)\sin {\theta _q}} \right)} \cdot \\\int_{ - \beta }^\beta {f\left( {\Delta {\theta _q}} \right)\exp \left( { - j{\rm{ \mathsf{ π} }}\left( {m - n} \right)\Delta {\theta _q}\cos {\theta _q}} \right){\rm{d}}\Delta {\theta _q}} = \\\sigma _{\rm{e}}^2{\delta _{mn}} + \sum\limits_{q = 1}^Q {\sigma _q^2\exp \left( { - j{\rm{ \mathsf{ π} }}\left( {m - n} \right)\sin {\theta _q}} \right)} \cdot \\\frac{{2\left[ {1 - \cos \left( {{\rm{ \mathsf{ π} }}\left( {m - n} \right)\beta \cos {\theta _q}} \right)} \right]}}{{{{\left[ {{\rm{ \mathsf{ π} }}\left( {m - n} \right)\beta \cos {\theta _q}} \right]}^2}}}.\end{array}$ (7)
因此，扩张矩阵TL的第(m, n)个元素可以表示为
${\mathit{\boldsymbol{T}}_L}\left( {m,n} \right) = \left\{ \begin{array}{l}\frac{{2\left[ {1 - \cos \left( {{\rm{ \mathsf{ π} }}\left( {m - n} \right)\beta \cos {\theta _q}} \right)} \right]}}{{{{\left[ {{\rm{ \mathsf{ π} }}\left( {m - n} \right)\beta \cos {\theta _q}} \right]}^2}}},m \ne n;\\1,m = n.\end{array} \right.$
可以看出当m≠n时，TL(m, n)与每个具体干扰信号来向都有关系，在满足零陷宽度要求的前提下，为了避免求解干扰信号来向信息，假设每个干扰信号来向变化对应的零陷宽度均为最大值，即用βmax替代βmaxcos θq，这样新的扩张矩阵TL的第(m, n)个元素可以表示为
${\mathit{\boldsymbol{T}}_L}\left( {m,n} \right) = \left\{ \begin{array}{l}\frac{{2\left[ {1 - \cos \left( {{\rm{ \mathsf{ π} }}\left( {m - n} \right){\beta _{\max }}} \right)} \right]}}{{{{\left[ {{\rm{ \mathsf{ π} }}\left( {m - n} \right){\beta _{\max }}} \right]}^2}}},m \ne n;\\1,m = n.\end{array} \right.$
由于新的扩张矩阵不需要包含干扰信号来向信息，式(7)可以进一步简化如下
$\begin{array}{l}{{\mathit{\boldsymbol{\tilde R}}}_L}\left( {m,n} \right) \approx \sigma _{\rm{e}}^2{\delta _{mn}} + {\mathit{\boldsymbol{T}}_L}\left( {m,n} \right) \cdot \\\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\sum\limits_{q = 1}^Q {\sigma _q^2\exp \left( { - {\rm{j \mathsf{ π} }}\left( {m - n} \right)\sin {\theta _q}} \right)} = \\\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\mathit{\boldsymbol{T}}_L}\left( {m,n} \right)\mathit{\boldsymbol{\tilde R}}\left( {m,n} \right).\end{array}$
定义“
${{\mathit{\boldsymbol{\tilde R}}}_L}\left( x \right) \approx \mathit{\boldsymbol{\tilde R}}\left( x \right) \circ {\mathit{\boldsymbol{T}}_L}.$
再根据式(6)，可得到协方差矩阵锥化后的阵列权矢量wT为
${\mathit{\boldsymbol{w}}_{\rm{T}}} = \frac{{{{\left( {\mathit{\boldsymbol{\tilde R}}\left( x \right) \circ {\mathit{\boldsymbol{T}}_L}} \right)}^{ - 1}}{\mathit{\boldsymbol{c}}_M}}}{{\mathit{\boldsymbol{c}}_M^{\rm{H}}{{\left( {\mathit{\boldsymbol{\tilde R}}\left( x \right) \circ {\mathit{\boldsymbol{T}}_L}} \right)}^{ - 1}}{\mathit{\boldsymbol{c}}_M}}}.$
以上是基于半波长均匀线阵推导出的关系式，实际上该推导过程适用于任意形状的阵列结构.
3 仿真实验实验采用七阵元半波长均匀线阵，BD2 B1频点信号来向为0°，信噪比为-20 dB；窄带干扰信号来向为-20°，干噪比为40 dB；BD2数字中频信号频率为2.42 MHz，信号采样率为8.184 MHz，快拍数为1 024.根据文献[15]中关于高动态的阐述，假设权值更新期间干扰来向变化幅度在3°以内；为消除随机误差影响，进行100次蒙特卡洛仿真实验.
仿真1??比较干扰信号来向变化统计模型分别服从三角形分布、均匀分布、拉普拉斯分布时协方差矩阵锥化零陷加宽算法的阵列方向增益特性.为了便于识别图形，将三角形分布、均匀分布、拉普拉斯分布对应的曲线分别标注为Tri、Uni、Lap.
图 2是PI算法、均匀分布算法、拉普拉斯分布算法和三角形分布算法对应的阵列天线方向增益图.首先，相比较于传统PI算法，3种分布都能有效展宽零陷.其次，就零陷深度在零陷宽度范围内的一致性而言，三角形分布和均匀分布算法对应的零陷深度在整个零陷范围内都保持较好，且三角形分布算法对应的零陷深度在整个零陷宽度上基本是优于均匀分布的；另外，虽然三角形分布算法对应的零陷极值没有拉普拉斯分布的零陷极值大，但是拉普拉斯分布的零陷深度也仅在不到1°幅度范围内保持优于三角形分布.
Fig. 2
图 2 幅度为3°时不同算法的阵列方向增益图 Fig. 2 Array patterns of different algorithms when the DOA of the interference varies within 3° 造成上述特点的原因可以从3种分布的概率密度曲线上进行解释，如图 3所示，均匀分布对应的概率密度曲线在整个变化区间上是等概率的，那么产生的零陷深度一致性保持要相对最好；三角形分布和拉普拉斯分布对应的概率密度曲线则在整个变化区间上呈现出中间高、两边低的特点，因此其零陷极值要大于均匀分布，拉普拉斯分布的概率密度曲线更为陡峭，其极值也更大，但是其在整个区间的零陷深度一致性保持上相对较差.由此可以看出，三角形分布在零陷深度以及宽度一致性之间要更加平衡，由于干扰来向变化规律实际未知，在拓展零陷范围内更加平稳的零陷深度将有利于干扰抑制，这是所提算法的优势.
Fig. 3
图 3 3种分布的概率密度曲线 Fig. 3 Probability density curves of different distributions 图 4是PI算法、均匀分布算法、拉普拉斯分布算法和三角形分布算法对应的阵列天线输出信干噪比随输入信噪比变化的比较图.首先，3种分布算法对应的阵列输出信干噪比要明显高于PI算法，可知零陷展宽算法不仅展宽了零陷，也提高了信号的输出信干噪比；其次，三角形分布对应的阵列输出信干噪比又要好于均匀分布、拉普拉斯分布.
Fig. 4
图 4 不同算法的阵列输出信干噪比线 Fig. 4 Output SINRs of different algorithms 仿真2??结合BD2软件接收机，验证所提算法的可行性.首先利用PI算法、三角形分布算法、均匀分布算法和拉普拉斯分布算法分别对包含静态干扰的BD2 B1频点数据进行处理，分别得到权矢量；为考察展宽性能，以静态干扰信号来向为初始值，根据以上的仿真条件构造高动态干扰条件下的BD2数据，利用得到的初始权矢量分别处理该数据，最后送入BD2软件接收机对其进行捕获.
图 5(a)显示的是利用PI算法对处理后数据的捕获结果，可以看出PI算法并未有效抑制高动态干扰条件下的干扰信号，因此接收机无法有效捕获BD2信号.图 5(b)、(c)、(d)显示的是分别利用三角形分布、均匀分布和拉普拉斯分布算法得到的捕获结果，可以看出本文所提Triangular-CMT-PI算法能够有效抑制高动态干扰条件下的干扰信号，从而使接收机成功捕获到了BD2信号，验证了算法的可行性.需要注意的是均匀分布算法与拉普拉斯算法也是可行的，这在之前的文献中已经被证明过了.
Fig. 5
图 5 不同算法处理后的捕获图 Fig. 5 Receiver acquisition with different algorithms 4 结论1) 本文提出的一种假设干扰来向变化服从三角形概率分布的协方差矩阵锥化零陷加宽抗干扰算法，能够在干扰来向上形成展宽零陷，在干扰来向快速变化的情况下保证其始终处于零陷之内.
2) 相比较于均匀分布算法和拉普拉斯分布算法，三角形分布算法拓宽后的干扰零陷更加平稳，阵列输出信干噪比要更高.正是由于干扰信号来向变化情况实际未知，平稳的零陷深度更加有利于抑制干扰信号，凸显了该算法的优点.
3) 通过与GNSS软件接收机结合，证明所提算法可以有效抑制干扰来向变化条件下的干扰信号，能够使软件接收机成功捕获卫星导航信号，具有一定的工程应用价值.
参考文献
[1] 梅浩, 战兴群, 李源, 等. 城市道路GNSS脆弱性评估技术[J]. 哈尔滨工业大学学报, 2017, 49(4): 101.
MEI Hao, ZHAN Xingqun, LI Yuan, et al. A GNSS vulnerability assessment technique of urban roads[J]. Journal of Harbin Institute of Technology, 2017, 49(4): 101. DOI:10.11918/j.issn.0367-6234.201601111
[2] 赵思浩, 陆明泉, 冯振明, 等. 基于自适应卡尔曼滤波的GNSS矢量锁定环路[J]. 哈尔滨工业大学学报, 2012, 44(7): 139.
ZHAO Sihao, LU Mingquan, FENG Zhenming, et al. GNSS vector lock loop based on adaptive Kalman filter[J]. Journal of Harbin Institute of Technology, 2012, 44(7): 139. DOI:10.11918/j.issn.0367-6234.2012.07.027
[3] COMPTON R T. The power-inversion adaptive array: concept and performance[J]. IEEE Transaction on Aerospace and Electronic Systems, 1979, AES-15(6): 803. DOI:10.1109/taes.1979.308765
[4] ALSHRAFI W, ENGEL U, BERTUCH T. Compact controlled reception pattern antenna for interference mitigation tasks of global navigation satellite system receiver[J]. IET Microwaves, Antennas & Propagation, 2015, 9(6): 593. DOI:10.1049/iet-map.2014.0357
[5] LI Qiang, XU Dingjie, WANG Wei, et al. Anti-jamming scheme for GPS receiver with vector tracking loop and blind beamformer[J]. Electronics Letters, 2014, 50(19): 1386. DOI:10.1049/el.2014.2274
[6] SHIN C G, JU J H, KANG D G, et al. Implementation of an antenna array for satellite communications with the capability of cancelling jammers[J]. IEEE Antennas and Propagation Magazine, 2013, 55(1): 32. DOI:10.1109/map.2013.6474483
[7] MAILLOUX R J. Covariance matrix augmentation to produce adaptive array pattern troughs[J]. Electronics Letters, 1995, 31(10): 771. DOI:10.1049/el:19950537
[8] ZATMAN M. Production of adaptive array troughs by dispersion synthesis[J]. Electronics Letters, 1995, 31(25): 2141. DOI:10.1049/el:19951486
[9] GUERCI J R. Theory and application of covariance matrix tapers for robust adaptive beamforming[J]. IEEE Transaction on Signal Processing, 1999, 47(4): 977. DOI:10.1109/78.752596
[10] ZATMAN M, GUERCI J R. Comments on "Theory and application of covariance matrix tapers for robust adaptive beamforming"[with reply][J]. IEEE Transaction on Signal Processing, 2000, 48(6): 1796. DOI:10.1109/78.845937
[11] 李荣锋, 王永良, 万山虎. 自适应天线方向图干扰零陷加宽方法研究[J]. 现代雷达, 2003, 25(2): 42.
LI Rongfeng, WANG Yongliang, WAN Shanhu. Research on adapted pattern null widening techniques[J]. Modern Radar, 2003, 25(2): 42. DOI:10.3969/j.issn.1004-7859.2003.02.012
[12] LU Dan, WU Renbiao, WANG Wenyi. Robust widen antijamming algorithm for high dynamic GPS[C]//Proceedings of the 11th International Conference on Signal Processing. Beijing: IEEE Press, 2012: 378. DOI: 10.1109/cosp.2012.6491679 https://ieeexplore.ieee.org/document/6491679
[13] 武思军, 张锦中, 张曙. 阵列波束的零陷加宽算法研究[J]. 哈尔滨工程大学学报, 2004, 25(5): 658.
WU Sijun, ZHANG Jinzhong, ZHANG Shu. Research on beamforming of widen nulling algorithm[J]. Journal of Harbin Engineering University, 2004, 25(5): 658. DOI:10.3969/j.issn.1006-7043.2004.05.025
[14] 张柏华, 马红光, 孙新利, 等. 高动态条件下统计空时零陷加宽方法[J]. 电子与信息学报, 2016, 38(4): 913.
ZHANG Baihua, MA Hongguang, SUN Xinli, et al. Space time null widening method of navigation receiver in missile for high dynamic conditions[J]. Journal of Electronics & Information Technology, 2016, 38(4): 913. DOI:10.11999/JEIT150654
[15] ZHANG Baihua, MA Hongguang, SUN Xinli, et al. Robust antijamming method for high dynamic global positioning system receiver[J]. IET Signal Processing, 2016, 10(4): 342. DOI:10.1049/iet-spr.2015.0122
[16] 王妙, 方明. 一种基于降秩共轭梯度法的零陷加宽技术[J]. 电子信息对抗技术, 2011, 26(4): 14.
WANG Miao, FANG Ming. A novel null broadening technique based on reducedrank conjugate gradient algorithm[J]. Electronic Information Warfare Technology, 2011, 26(4): 14. DOI:10.3969/j.issn.1674-2230.2011.04.004
[17] 柯熙政, 刘娜. 卫星导航自适应抗干扰算法及零陷加宽研究[J]. 电子测量与仪器学报, 2010, 24(12): 1082.
KE Xizheng, LIU Na. Research on adaptive antijamming and null widening algorithm for satellite navigation[J]. Journal of Electronic Measurement and Instrument, 2010, 24(12): 1082. DOI:10.3724/SP.J.1187.2010.01082
