删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

一种针对GNSS接收机的宽零陷抗干扰算法

本站小编 哈尔滨工业大学/2019-10-24

一种针对GNSS接收机的宽零陷抗干扰算法

王海洋,刘光斌,范志良,姚志成

(火箭军工程大学 导弹工程学院, 西安 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, 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.

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]HM维约束列向量;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}$进行替代,样本协方差矩阵的表达式如下

${{\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}}}$

$\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为干扰信号初始来向;Δθ为干扰信号来向变化幅度;${\tilde \theta }$为干扰信号来向变化后的真实值,其规律是未知的.

假设第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}},$

则可得${{\mathit{\boldsymbol{\tilde R}}}_L}\left( x \right)$第(m, n)个元素的表达式为

$\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.$

可以看出当mn时,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}$

定义“$ \circ $”表示Hardmard积,则采样协方差矩阵${{\mathit{\boldsymbol{\tilde R}}}}\left( x \right)$与锥化后的采样协方差矩阵${{\mathit{\boldsymbol{\tilde R}}}_L}\left( x \right)$满足下述关系:

${{\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



相关话题/干扰 信号 统计 软件 概率

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 合成生物学基因设计软件:iGEM设计综述
    合成生物学基因设计软件:iGEM设计综述伍克煜1,刘峰江1,许浩1,张浩天1,王贝贝1,2(1.电子科技大学生命科学与技术学院,成都611731;2.电子科技大学信息生物学研究中心,成都611731)摘要:随着基因回路规模的扩大,和应用范围的拓展,传统的合成基因回路的设计思路面临着新的挑战。新合成基 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • Hsa-miR-210-5p靶基因预测及其相关信号通路的生物信息学分析
    Hsa-miR-210-5p靶基因预测及其相关信号通路的生物信息学分析蔡丹平,龙鼎新(南华大学船山学院,湖南衡阳421001)摘要:为深入研究miR-210-5p的调控机制及生物学功能提供理论机制,应用生物信息学方法分析miR-210-5p序列,预测其靶基因,用Veney2.1.0绘制韦恩图得到靶基 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 肺泡发育及其相关信号通路
    肺泡发育及其相关信号通路陈俏媛,林万华(广西高校干细胞与医药生物技术重点实验室(广西师范大学生命科学学院),广西桂林541004)摘要:肺是由分支的气道和血管组成的复杂结构,它们结合在最远端的肺泡进行气体交换。而肺泡则是由鳞状肺泡I型(AT1)细胞和立方状肺泡II型(AT2)细胞组成的气体交换囊,其 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 桥梁颤振临界风速的概率密度演化计算
    桥梁颤振临界风速的概率密度演化计算姜保宋,周志勇,唐峰(土木工程防灾国家重点实验室(同济大学),上海200092)摘要:针对桥梁结构自身特性以及外部环境的随机性(如刚度、质量、阻尼比、气动导数等因素)所造成的桥梁的颤振临界风速不确定,难以衡量桥梁颤振稳定性问题.将概率密度演化方法与桥梁颤振多模态耦合 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 城市交通信号自组织控制规则的邻域重构
    城市交通信号自组织控制规则的邻域重构钟馥声,王安麟,姜涛,花彬(同济大学机械与能源工程学院,上海201804)摘要:为解决城市交通信号自组织控制中,交通流相变频繁及路网拓扑结构复杂等时空条件带来的固定局部规则控制精度的问题,提出一种城市交通信号自组织控制规则邻域的重构方法.首先,定义邻域为当前路口自 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 陀螺飞轮信号的EMD/LPF混合去噪方法
    陀螺飞轮信号的EMD/LPF混合去噪方法赵昱宇,赵辉,霍鑫,姚郁(哈尔滨工业大学航天学院,哈尔滨150001)摘要:为提高陀螺飞轮系统的标定与辨识精度进而保证姿态测量性能,对其信号去噪方法进行研究以便从复杂噪声中提取标定辨识所需的有用信息.首先基于噪声产生机理对陀螺飞轮信号的噪声特性进行分析.其次, ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 概率语言环境下考虑专家心理行为的QFD方法
    概率语言环境下考虑专家心理行为的QFD方法鞠萍华,陈资,冉琰,涂顺泽(机械传动国家重点实验室(重庆大学),重庆400044)摘要:为解决传统质量功能展开(QFD)在实际运用过程中存在关于顾客需求和工程技术之间的关系评估,顾客需求权重的确定和工程技术的优先级排序等方面的固有缺陷,提出一种概率语言环境下 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 一种故障干扰解耦的航天器主动容错控制方法
    一种故障干扰解耦的航天器主动容错控制方法宗群1,杨希成1,张秀云1,刘文静2(1.天津大学电气自动化与信息工程学院,天津300072;2.空间智能控制技术重点实验室(北京控制工程研究所),北京100190)摘要:为提高航天器系统的可靠性,研究一种干扰影响下的航天器主动容错控制技术.首先,为实现干扰下 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 声振信号联合1D-CNN的大型电机故障诊断方法
    声振信号联合1D-CNN的大型电机故障诊断方法赵书涛1,王二旭1,陈秀新2,王科登1,李小双1(1.华北电力大学电气与电子工程学院,河北保定071003;2.华北电力大学控制与计算机工程学院,河北保定071003)摘要:针对复杂运行环境下大功率电动机故障诊断准确率不高、算法泛化能力差的问题,提出一种 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 立交区域交叉口交通信息识别概率预测
    立交区域交叉口交通信息识别概率预测沈强儒1,杨少伟2,曹慧1,顾镇媛1,葛婷3(1.南通大学交通与土木工程学院,江苏南通,226019;2.长安大学公路学院,西安710064;3.苏州科技大学土木工程学院,江苏苏州215011)摘要:为预测立交区域交叉口交通信息识别概率,运用汽车动力学理论、驾驶员特 ...
    本站小编 哈尔滨工业大学 2020-12-05