CN106680837B - 一种卫星导航干扰抑制算法 - Google Patents
一种卫星导航干扰抑制算法 Download PDFInfo
- Publication number
- CN106680837B CN106680837B CN201611149124.1A CN201611149124A CN106680837B CN 106680837 B CN106680837 B CN 106680837B CN 201611149124 A CN201611149124 A CN 201611149124A CN 106680837 B CN106680837 B CN 106680837B
- Authority
- CN
- China
- Prior art keywords
- weight vector
- signal
- covariance matrix
- interference
- estimated value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001629 suppression Effects 0.000 title claims abstract description 15
- 239000013598 vector Substances 0.000 claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 6
- 238000010606 normalization Methods 0.000 claims description 4
- 238000000034 method Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 4
- 241001463014 Chazara briseis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种卫星导航干扰抑制算法,首先计算接收数据协方差矩阵的估计值,然后对接收数据协方差矩阵的估计值进行特征分解运算,获得特征值和相应的特征向量,将最小特征值作求倒数运算并与其对应的特征向量相乘得到未归一化的权向量,并进行归一化处理,采用权向量对接收到的中频数字复数信号进行加权叠加处理,得到干扰抑制后的输出信号。本发明可以同时对强干扰信号和弱干扰信号进行有效地抑制。
Description
技术领域
本发明属于卫星导航抗干扰领域,是一种对压制式卫星导航干扰信号实现抑制的算法。
背景技术
在卫星导航抗干扰天线***的信号处理过程中,干扰抑制算法是整个数字信号处理中核心的部分。干扰抑制算法可以分为自适应调零类算法和最优数字多波束类算法两大类。
自适应调零类算法的典型代表为功率倒置(Power Inversion,PI)算法,该算法不需要卫星方位信息等先验信息,零陷深度随干扰功率的强度自动调节,对于强干扰信号有较好的抑制效果,因而在卫星导航抗干扰领域的工程应用中得到了广泛地使用。但是,其缺点是对于弱干扰信号的抑制效果较差,使得其在弱干扰环境或同时存在强干扰和弱干扰的环境中的应用受到限制。
最优数字多波束算法采用数字复用技术,同时在视野空间内形成多个最优数字接收波束。典型的最优数字波束为最小方差无失真响应响应(Minimum Variance DistortedResponse,MVDR)波束。每个最优数字波束的主瓣指向一颗导航卫星,同时自适应地在干扰方向形成零陷。最优数字多波束算法在对干扰抑制的同时可以提升***输出信号的信噪比。但是该算法需要卫星方位、阵列姿态等先验信息辅助,算法结构复杂、计算量较大,很难在工程应用中得到广泛使用。
发明内容
为了克服现有技术的不足,本发明提供一种卫星导航干扰抑制算法,采用归一化后的接收数据协方差矩阵的噪声子空间特征向量作为自适应调零权向量对接收信号中的干扰分量进行抑制,可以同时对接收信号中的弱功率干扰信号分量和强功率干扰进行有效地抑制,易于工程实现。
本发明解决其技术问题所采用的技术方案包括以下步骤:
第一步,计算接收数据协方差矩阵Rx的估计值其中,x(n)为由M个阵元的阵列天线在n时刻接收到的中频数字复数信号向量,维数为M×1,M≥2;的维数大小为M×M;上标H为共轭转置运算符;N为计算接收数据协方差矩阵的估计值需要的信号快拍数,N≥4M;
第二步,对接收数据协方差矩阵的估计值采用Jacobi算法进行特征分解运算,获得特征值λm和相应的特征向量qm,m=1,2,…,M;
第三步:将最小特征值λmin作求倒数运算并与其对应的特征向量相乘得到未归一化的权向量 的维数大小为M×1;
第四步,对权向量进行归一化处理,得到归一化的权向量其中,为权向量的第一个元素;
第五步,采用权向量wopt对第一步中接收到的中频数字复数信号x(n)进行加权叠加处理,干扰抑制后的输出信号
本发明的有益效果是:采用本发明得到的权向量wopt省去对应于干扰特征向量的权分量,使得加权后输出的信号功率中只包含噪声信号分量,可以同时对强干扰信号和弱干扰信号进行有效地抑制。
附图说明
图1是本发明的算法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明算法实现步骤如下:
第一步:计算接收数据协方差矩阵Rx的估计值计算表达式如下:
其中,x(n)为由M(M≥2)个阵元的阵列天线在n时刻接收到的中频数字复数信号向量,维数为M×1;的维数大小为M×M;上标H为共轭转置运算符;N为计算接收数据协方差矩阵的估计值需要的信号快拍数,需满足N≥4M。
第二步:对接收数据协方差矩阵的估计值进行特征分解运算,特征分解采用Jacobi算法实现,获得特征值λm(m=1,2,…,M)和相应的特征向量qm(m=1,2,…,M)。
第三步:将最小特征值λmin作求倒数运算并与其对应的特征向量相乘得到未归一化的权向量
的维数大小为M×1。
第四步:对第三步中得到的未归一化的权向量进行归一化处理,得到归一化的权向量wopt,计算过程如下式所示:
其中,为权向量的第一个元素。
第五步:采用第四步中计算得到的权向量wopt对第一步中接收到的中频数字复数信号x(n)进行加权叠加处理,计算过程如下式所示:
其中,上标H为共轭转置运算符,y为干扰抑制后的输出信号。
对于功率倒置算法,其权向量计算如下:
其中,wPI为功率倒置算法计算得到的权向量,维数大小为M×1;为收数据协方差矩阵估计值的逆矩阵,维数大小为M×M;a0=[1,0,…,0]T为约束向量,维数大小为M×1;K为干扰源数目,且K<M;δ2为***噪声功率;wjam和wnoise分为干扰信号和噪声信号两部分分别对应的权向量,维数大小为M×1。
当外界存在强干扰信号时(干扰信号的功率远大于噪声信号功率),干扰信号对用的特征值将远大于噪声信号对应的特征值,即
由式(6)可知,随着干扰功率的增大,最优权值中对应于干扰特征向量的权分量将越来越小。因此,功率倒置法可以在强干扰方向上形成较深的零陷,抗干扰能力较强。反之,若干扰信号的功率较小(略大于噪声信号)时,最优权值中对应于干扰特征向量的权分量(式(5)中第二个等号右侧第一项)将略小于对应于噪声特征分量的权向量(式(5)中第二个等号右侧第二项),形成的零陷变浅,抗干扰能力较差。而采用本发明得到的权向量wopt省去对应于干扰特征向量的权分量,使得加权后输出的信号功率中只包含噪声信号分量,可以同时对强干扰信号和弱干扰信号进行有效地抑制。
本发明适用于GPS、BDS和GLONASS卫星导航***干扰抑制处理。以4单元BD2-B3频点抗干扰天线抗两个宽带干扰为例阐述本发明的具体实施方式。
步骤1:计算4元阵列接收数据协方差矩阵的估计值数据采样快拍数N=256,计算得到的是维数为4×4的Hermit矩阵,计算过程如下:
步骤2:采用Jacobi算法对接收数据协方差矩阵的估计值进行特征分解运算,获得特征值λm(m=1,2,…,4)和相应的特征向量qm(m=1,2,…,4)。
步骤3:将步骤2计算得到的最小特征值λmin作求倒数运算并与其对应的特征向量相乘得到未归一化的权向量计算过程如下:
步骤4:对步骤3中得到的未归一化的权向量进行归一化处理,得到归一化的权向量wopt,计算过程如式(3)所示。
步骤5:采用步骤4中计算得到的权向量wopt对第一步中接收到的中频数字复数信号x(n)进行加权叠加处理并输出,计算过程如式(4)所示。
至此,完成对4单元BD2-B3频点接收信号的干扰抑制处理。
Claims (1)
1.一种卫星导航干扰抑制算法,其特征在于包括下述步骤:
第一步,计算接收数据协方差矩阵Rx的估计值其中,x(n)为由M个阵元的阵列天线在n时刻接收到的中频数字复数信号向量,维数为M×1,M≥2;的维数大小为M×M;上标H为共轭转置运算符;N为计算接收数据协方差矩阵的估计值需要的信号快拍数,N≥4M;
第二步,对接收数据协方差矩阵的估计值采用Jacobi算法进行特征分解运算,获得特征值λm和相应的特征向量qm,m=1,2,…,M;
第三步:将最小特征值λmin作求倒数运算并与其对应的特征向量相乘得到未归一化的权向量 的维数大小为M×1;
第四步,对权向量进行归一化处理,得到归一化的权向量其中,为权向量的第一个元素;
第五步,采用权向量wopt对第一步中接收到的中频数字复数信号x(n)进行加权叠加处理,干扰抑制后的输出信号
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611149124.1A CN106680837B (zh) | 2016-12-14 | 2016-12-14 | 一种卫星导航干扰抑制算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611149124.1A CN106680837B (zh) | 2016-12-14 | 2016-12-14 | 一种卫星导航干扰抑制算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106680837A CN106680837A (zh) | 2017-05-17 |
CN106680837B true CN106680837B (zh) | 2019-04-19 |
Family
ID=58868312
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611149124.1A Active CN106680837B (zh) | 2016-12-14 | 2016-12-14 | 一种卫星导航干扰抑制算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106680837B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108089162A (zh) * | 2017-12-29 | 2018-05-29 | 中国电子科技集团公司第二十研究所 | 一种脉冲干扰信号的检测与抑制方法 |
CN108241160A (zh) * | 2017-12-29 | 2018-07-03 | 中国电子科技集团公司第二十研究所 | 一种卫星导航干扰源空间谱峰估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630910A (zh) * | 2013-12-13 | 2014-03-12 | 武汉大学 | 一种gnss接收机设备的抗干扰方法 |
CN104536018A (zh) * | 2015-01-06 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种使用阵列天线抗干扰技术的gnss多星联合捕获方法 |
CN105137454A (zh) * | 2015-07-22 | 2015-12-09 | 北京航空航天大学 | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 |
-
2016
- 2016-12-14 CN CN201611149124.1A patent/CN106680837B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630910A (zh) * | 2013-12-13 | 2014-03-12 | 武汉大学 | 一种gnss接收机设备的抗干扰方法 |
CN104536018A (zh) * | 2015-01-06 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种使用阵列天线抗干扰技术的gnss多星联合捕获方法 |
CN105137454A (zh) * | 2015-07-22 | 2015-12-09 | 北京航空航天大学 | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 |
Non-Patent Citations (1)
Title |
---|
"一种改进的功率倒置卫星导航抗干扰方法";江城 等;《现代导航》;20160831(第4期);第1-2节 |
Also Published As
Publication number | Publication date |
---|---|
CN106680837A (zh) | 2017-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109254261B (zh) | 基于均匀圆阵epuma的相干信号零陷加深方法 | |
CN105302936B (zh) | 基于相关计算和协方差矩阵重构的自适应波束形成方法 | |
Zhang et al. | Robust beamforming for coherent signals based on the spatial-smoothing technique | |
CN105204008A (zh) | 一种基于协方差矩阵扩展的自适应天线波束形成零陷展宽方法 | |
Wang et al. | Interference cancellation in GPS receiver using noise subspace tracking algorithm | |
Li et al. | On adaptive beamforming for coherent interference suppression via virtual antenna array | |
CN106680837B (zh) | 一种卫星导航干扰抑制算法 | |
CN113504549B (zh) | 基于广义旁瓣相消器的导航空时抗干扰方法 | |
CN108415040B (zh) | 一种基于子空间投影的csmg波束形成方法 | |
Wang et al. | Eigenspace-based beamforming technique for multipath coherent signals reception | |
CN107332601B (zh) | 一种基于波束赋形的自适应干扰抑制方法 | |
Redif et al. | Robust broadband adaptive beamforming via polynomial eigenvalues | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN104880714B (zh) | 一种基于两级滤波结构的gnss天线阵抗干扰方法 | |
CN111257863A (zh) | 一种高精度多点线性约束的自适应单脉冲测向方法 | |
Zhang et al. | Adaptive multiple-beamformers for reception of coherent signals with known directions in the presence of uncorrelated interferences | |
Miranda et al. | High accuracy and low complexity adaptive Generalized Sidelobe Cancelers for colored noise scenarios | |
Barnov et al. | Spatially robust GSC beamforming with controlled white noise gain | |
Liu et al. | A virtual space-time adaptive beamforming method for space-time antijamming | |
CN108717196B (zh) | 一种阵列天线接收信号的去干扰方法及*** | |
Cai et al. | Low-complexity reduced-dimension space–time adaptive processing for navigation receivers | |
Dungriyal et al. | Performance of MIR-LMS algorithm for adaptive beam forming in smart antenna | |
Bucris et al. | Robust focusing for wideband MVDR beamforming | |
KR101206297B1 (ko) | 간섭 널 투사에 의한 신호차단 방식의 적응 빔 형성 장치 | |
Luo et al. | Researchment of mainlobe anti-jamming algorithms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |