CN107728166A - 一种基于时域分组处理的卫星导航接收机多干扰抑制方法 - Google Patents
一种基于时域分组处理的卫星导航接收机多干扰抑制方法 Download PDFInfo
- Publication number
- CN107728166A CN107728166A CN201710806427.4A CN201710806427A CN107728166A CN 107728166 A CN107728166 A CN 107728166A CN 201710806427 A CN201710806427 A CN 201710806427A CN 107728166 A CN107728166 A CN 107728166A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msup
- 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.)
- Granted
Links
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
-
- 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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Noise Elimination (AREA)
Abstract
本发明公开了一种基于时域分组处理的卫星导航接收机多干扰抑制方法,属于射频干扰抑制领域。所述方法包括:利用N阵元天线阵获得N通道的接收信号数据;采用改进的奇异值比谱方法估计接收信号中各周期分量的周期,进而计算所有周期分量的最小公周期;根据估计的周期信息将各通道接收信号中具有相同频率特性的数据点重组得到新的子空时数据矩阵;利用基于最小功率无畸变准则的空时处理算法完成各组信号的干扰消除;将消除干扰后的数据进行重构以得到后续处理所需信号。本申请发明的方法能够根据周期调频信号的时频稀疏性与广义周期特性,在不增加接收机阵元个数的前提下,提升接收机干对抗宽带调频干扰的个数。
Description
技术领域
本发明属于射频干扰抑制领域,具体涉及一种基于时域分组处理的卫星导航接收机多干扰抑制方法。
背景技术
无论是GPS,GLONASS,GALILEO还是“北斗”卫星导航***,都将存在一个共同的缺点——容易受到多种形式的有意或无意的干扰,导致接收机定位、导航性能下降,甚至无法正常工作。而且,随着卫星导航***的在军事、民用领域的广泛应用,接收机所处的工作环境日益复杂,卫星导航接收机的抗干扰能力已成为其能否正常工作的关键问题。研究卫星导航***接收机抗干扰技术,对于提高卫星导航***在复杂电磁环境下的工作性能、增强其在各种环境中的适用性具有重要意义,对卫星导航***用户端接收机的研究已成为各国学者关注的热点之一。
由于卫星导航***最初的研制目标是为动态用户提供全天候、实时、高精度的位置信息,没有充分考虑用户端接收机在各种复杂电磁环境下的工作问题,而且由于卫星信号自身所具有的特性,致使用户端对卫星信号进行接收时存在以下问题:
卫星所发射信号的载波频率固定、调制方式已知、数据传输格式不变,致使用户端接收机很容易受到人为干扰的影响。如果将干扰的频率集中对准卫星信号的载波频率,将导致接收机失锁。而且,一般的用户端接收机只辨别卫星信号的结构,无法识别卫星信号的真伪,当干扰与卫星信号结构相同时,干扰与卫星信号会同时进入用户端接收机,导致接收机误捕获。卫星距离地球遥远,到达用户端接收机的卫星信号十分微弱。根据信号传播损耗理论,信号在空间传播的功率与其传播距离的平方成反比,传播距离越远,信号功率衰减越大。导航卫星一般位于距离地球表面20000km乃至更高的轨道平面上,而且卫星发射机功率不可能很大,信号到达用户端接收机时已相当微弱。卫星信号本身抗干扰裕度不大。卫星信号多为扩频信号,具有一定抗干扰能力,但当干扰功率大于抗干扰裕度时,用户端接收机的性能就会受到显著影响。
随着全球卫星导航***(GNSS)的蓬勃发展以及其在国防和国民经济各个领域的广泛应用,卫星导航信号的易损性和脆弱性所带来的问题也日益凸显,为了保障GNSS服务的可靠性、精确性和连续性,大量的学者对该领域进行了研究并取得了一系列成果。
根据干扰算法所需天线阵元数,可以将现有的干扰抑制算法分为:单天线干扰抑制技术、基于阵列天线的干扰抑制技术。单天线干扰抑制技术(时域/频域、时频域)干扰抑制方法,具有体积小、硬件复杂度低的优点;但他们只能处理具有时频稀疏性的干扰,且对多个干扰处理能不强。基于阵列天线的空域干扰抑制技术,可以在不考虑干扰信号时频特性的条件下,利用空域滤波方法对抗宽带或者窄带干扰,但其可对抗干扰的个数受限于阵元数目。为了在不增加天线阵元个数的前提下提高干扰抑制自由度,基于空时处理的干扰抑制算法被引入GNSS干扰抑制领域并得到了广泛研究。但是该类方法只增加了抑制窄带干扰的个数,无法处理过多的宽带干扰。
发明内容
本发明的目的在于提供在不增加接收机阵元个数的前提下,提升接收机干对抗宽带调频干扰的个数的一种基于时域分组处理的卫星导航接收机多干扰抑制方法。
本发明的目的通过以下技术方案来实现:
一种基于时域分组处理的卫星导航接收机多干扰抑制方法,包括以下步骤:
(1.1)利用阵元数为N的阵列天线,经过射频处理,获得N通道M采样点的数字接收信号x(m)=[x1(m) x2(m) L xN(m)]T,其中m=1,2,L M;
(1.2)选取某一通道信号xn,计算接收信号中各周期分量的最小公周期Tc;
(1.3)根据估计的小公周期Tc,将具有相似时频特性的信号数据进行重组,获得Q组空时域数据,X1,X2,…XQ;
(1.4)对各组信号数据(X1,X2,…XQ)进行空时滤波处理,获得消除干扰后信号数据Y1,Y2,…YQ。
(1.5)对各组干扰抑制后数据进行重构,得到后续处理所需要的数据Y。
本发明的核心技术内容在于将利用各周期干扰信号的周期特性将具有相似时频特性的信号数据进行重组,在分组进行干扰消除,从而提出了一种基于时域分组处理的卫星导航接收机多干扰抑制方法。
所述的选取某一通道信号xn,计算接收信号中各周期分量的最小公周期Tc,主要内容为:利用奇异值比谱方法估计得到接收信号xn中各周期分量的周期Tk,其中k=1,2,…,K为周期分量的序号;进而计算所有周期分量周期的最小公倍数,获得最小公周期Tc;
所述的利用奇异值比谱方法估计得到接收信号xn中各周期分量的周期Tk,其主要内容为:构建信号xn的奇异值比谱,然后提取大于门限值峰值序列并所得峰值进行周期序列检测,每个序列的时间间隔对应于一个周期调频信号的调频周期Tk。该方法避免了传统基于奇异值比谱周期估计方法所估计的周期可能是真实周期的整数倍的问题。
所述的据估计的小公周期Tc,将具有相似时频特性的信号数据进行重组,其主要内容为:根据所估计的最小公周期Tc,将各通道信号们可以将时域数据进行分割,再把具有相似时频特性的信号数据进行重组。该方法可以获得足够具有同时频特性的采样数据,以利于获得更精确的相应时频信号的协方差矩阵信息。
所述的对各组信号数据进行空时滤波处理,获得消除干扰后信号数据,其主要内容为:对每组数据利用基于最小输出功率准则的空时滤波器进行干扰抑制处理。
本发明的有益效果在于:
将时域分组处理方式引入干扰抑制领域,利用周期调频信号的时频稀疏性与广义周期特性,通过将同时频特性的数据点重组得到新的子空时数据矩阵,然后对分组后信号的进行干扰消除;从而形成一种基于时域分组处理的卫星导航接收机多干扰抑制方法。
能够充分利用周期调频信号在时频域的聚集性和周期性,构建空时频矩阵以避免对空域自由度的重复利用,然后再利用干扰信号在空域的稀疏性进行干扰消除。所以该方法能够在不增加阵列天线阵元个数的前提下,增加接收机抗干扰个数。
附图说明
表1为干扰信号参数;
图1为基于时域分组处理的卫星导航接收机多干扰抑制方法流程图;
图2为利用不同算法干扰抑制后C/A码捕获结果对比。
具体实施方式
下面结合附图对本发明的具体实施内容作进一步说明:
本申请实施例在分析周期调频信号的时频稀疏性与广义周期特性的基础上,提出一种基于时域数据重组的空时干扰抑制算法。所申请方法能够根据周期调频信号的时频稀疏性与广义周期特性,在不增加接收机阵元个数的前提下,提升接收机干对抗宽带调频干扰的个数。
为更加清晰的说明所申请方法,本申请实施例通过仿真实验进行流程说明以及效果展示,但不限制本申请实施例的范围。实验条件为:天线阵选取间隔为半波长的均匀线阵,天线单元数N=4。期望信号为全球定位***(GPS)的C/A码,中频中心频为1.023MHz,采样频率为4.096MHz;接收信号的入射方向为80°;信噪比为-15dB。干扰信号的干噪比为60dB,其余参数如表1所示。
图1是基于时域数据重组的空时GNSS多干扰干扰抑制算法的流程图,该方法包括:
S110获取N(N=4)个通道的信号数据:
x(m)=[x1(m) x2(m) L xN(m)]T
式中,m代表第m个采样时刻,T代表转置运算。
S120选取任意一个通道的信号数据,估计接收信号中周期干扰信号的最小公周期Tc;
步骤S120可进一步包括:
S121将所选第n(n=1)通道信号数据按照长度Ω连续截取2段构造矩
对O进行奇异值分解,构建奇异值比谱。传统周期分量估计方法是以最大峰值的对应的时间点作为估计的周期,但是每个周期信号对应的奇异值比谱峰值周期性出现,以峰值大小作为检测准则,可能会导致所估计的周期是真实周期的整数倍。针对该问题,提出峰值序列检测法以估计各周期分量信号的周期。首先,提取大于门限值峰值序列,
其中,ρ为门限值,其定义为:
S122对所得峰值进行周期序列检测,搜索具有相同时间间隔的峰值序列,每个序列的时间间隔对应一个周期调频信号的广义周期Tk。进而根据所得的周期计算最小公周期:
Tc=T1,T2,…,TK
●表示求一列数的公倍数。
S130利用估计出公周期Tc,构建空时数据矩阵;
将数据分为Q组空域数据,连续数据的个数为p(为正整数),则第n通道第q组空域数据为:
其中,TcG为所截取的信号长度;
则每个通道有Mτ-1个时域抽头,则第q组空时空时矩阵可以表述为
S140对每组数据可以利用最小输出功率准则准则进行空域滤波处理,最小输出功率准则准则为:
其中,wq为第q组数据对应的空域权矢量,Rq为根据最大似然准则估计获得的第q组数据的协方差矩阵:
c为空时约束向量,定义为:
其中,L为数据长度,最优权值为:
干扰抑制后数据输出为:
S150最后对各组干扰抑制后数据进行重组,得到后续处理所需要的数据:
Y=[Y1,1Y1,2…YQ,1 Y2,1Y2,2…YQ,2…Y1,GY1,G…YQ,G]
其中
干扰抑制后,导航信号捕获结果由图2给出,为了说明所申请算法的有效性,将所申请的基于时域数据重组的空时干扰抑制算法与经典空域算法和空时算法进行比较,可以发现本发明所提出方法可以不增加阵列天线阵元个数的前提下,增加接收机抗干扰个数。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于时域分组处理的卫星导航接收机多干扰抑制方法,其特征在于,包括以下步骤:
(1)利用阵元数为N的阵列天线,经过射频处理,获得N通道M采样点的数字接收信号数据矩阵x;x(m)=[x1(m) x2(m) L xN(m)]T,其中xn表示第n通道信号数据向量;m=1,2,L M,表示采样时刻;
(2)选取某一通道信号xn,计算接收信号中各周期分量的最小公周期Tc;
(3)根据估计的小公周期Tc,将具有相似时频特性的信号数据进行重组,获得Q组空时域数据,X1,X2,…XQ;
(4)对各组信号数据(X1,X2,…XQ)进行空时滤波处理,获得消除干扰后信号数据Y1,Y2,…YQ;
(5)对各组干扰抑制后数据进行重构,得到后续处理所需要的数据Y。
2.根据权利要求1所述一种基于时域分组处理的卫星导航接收机多干扰抑制方法,其特征在于,所述的步骤(2)具体为:
(2.1)将所选第n通道信号数据按照长度Ω连续截取2段构造矩
<mrow>
<mi>O</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Omega;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Omega;</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&Omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
对O进行奇异值分解,构建奇异值比谱;
(2.2)提取大于门限值峰值序列,
<mrow>
<msub>
<mi>R</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>&rho;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>></mo>
<mi>&rho;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,ρ为门限值,其定义为:
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>R</mi>
<mi>&sigma;</mi>
</msub>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
(2.3)对所得峰值进行周期序列检测,搜索具有相同时间间隔的峰值序列,每个序列的时间间隔对应一个周期调频信号的广义周期Tk,
进而根据所得的周期计算最小公周期:
Tc=T1,T2,…,TK
●表示求一列数的公倍数。
3.根据权利要求1所述一种基于时域分组处理的卫星导航接收机多干扰抑制方法,其特征在于,所述的步骤(3)具体为:
(3.1)将数据分为Q组空域数据,连续数据的个数为p(为正整数),则第n通道第q组空域数据为:
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>n</mi>
<mi>q</mi>
</msub>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<mi>G</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
其中,TcG为所截取的信号长度;
<mrow>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>T</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
(3.2)每个通道有Mτ-1个时域抽头,则第q组空时空时矩阵可以表述为
<mrow>
<msub>
<mi>X</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>1</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>1</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<mn>...</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>1</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<msub>
<mi>M</mi>
<mi>&tau;</mi>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>2</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>2</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mn>2</mn>
<mi>q</mi>
</msub>
<mo>,</mo>
<msub>
<mi>M</mi>
<mi>&tau;</mi>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>q</mi>
</msub>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>n</mi>
<mi>q</mi>
</msub>
<mo>,</mo>
<msub>
<mi>M</mi>
<mi>&tau;</mi>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>.</mo>
</mrow>
4.根据权利要求1所述一种基于时域分组处理的卫星导航接收机多干扰抑制方法,其特征在于,所述的步骤(4)具体为:
(4.1)根据最大似然准则估计获得的第q组数据的协方差矩阵Rq:
<mrow>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<msub>
<mi>X</mi>
<mi>q</mi>
</msub>
<msub>
<msup>
<mi>X</mi>
<mi>H</mi>
</msup>
<mi>q</mi>
</msub>
<mo>;</mo>
</mrow>
(4.2)利用最小输出功率准则进行空域滤波处理,最小输出功率准则准则为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>w</mi>
<mi>q</mi>
</msub>
<mi>H</mi>
</msup>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
<msub>
<mi>w</mi>
<mi>q</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>w</mi>
<mi>q</mi>
</msub>
<mi>H</mi>
</msup>
<mi>c</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中c为空时约束向量,定义为:
<mrow>
<mi>c</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>a</mi>
<mi>T</mi>
</msup>
</mtd>
<mtd>
<msup>
<munder>
<mn>0</mn>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</munder>
<mi>T</mi>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<munder>
<mn>0</mn>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</munder>
<mi>T</mi>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
<mi>a</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<munder>
<mn>0</mn>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</munder>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
利用拉格朗日法,可求解最优解为:
<mrow>
<msub>
<mi>w</mi>
<msub>
<mi>q</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>c</mi>
</mrow>
<mrow>
<msup>
<mi>c</mi>
<mi>H</mi>
</msup>
<msup>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>c</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
(4.3)干扰抑制后数据输出为:
<mrow>
<msub>
<mi>Y</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>w</mi>
<msub>
<mi>q</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
</msub>
<mi>H</mi>
</msup>
<msub>
<mi>X</mi>
<mi>q</mi>
</msub>
<mo>.</mo>
</mrow>
5.根据权利要求1,2所述一种基于时域分组处理的卫星导航接收机多干扰抑制方法,其特征在于,所述各组干扰抑制后数据进行重构中:
对各组干扰抑制后数据进行重组,得到后续处理所需要的数据:
Y=[Y1,1Y1,2…YQ,1Y2,1Y2,2…YQ,2…Y1,GY1,G…YQ,G]
其中
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>q</mi>
<mo>,</mo>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mi>g</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>p</mi>
<mo>+</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710806427.4A CN107728166B (zh) | 2017-09-08 | 2017-09-08 | 一种基于时域分组处理的卫星导航接收机多干扰抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710806427.4A CN107728166B (zh) | 2017-09-08 | 2017-09-08 | 一种基于时域分组处理的卫星导航接收机多干扰抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107728166A true CN107728166A (zh) | 2018-02-23 |
CN107728166B CN107728166B (zh) | 2020-12-22 |
Family
ID=61205019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710806427.4A Active CN107728166B (zh) | 2017-09-08 | 2017-09-08 | 一种基于时域分组处理的卫星导航接收机多干扰抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107728166B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110320535A (zh) * | 2019-07-03 | 2019-10-11 | 哈尔滨工程大学 | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 |
CN112083448A (zh) * | 2020-09-04 | 2020-12-15 | 哈尔滨工程大学 | 面向卫星导航***的干扰信号分类识别特征提取方法及*** |
CN114675307A (zh) * | 2022-03-14 | 2022-06-28 | 合肥工业大学 | 基于分数阶伪Wigner-Ville分布的卫星导航接收机干扰检测方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1463105A (zh) * | 2003-06-18 | 2003-12-24 | 清华大学 | 用于异步码分多址***的分组单次干扰抵消方法 |
CN101114863A (zh) * | 2006-07-28 | 2008-01-30 | 美国博通公司 | 通信***内处理信号的方法和*** |
CN101615982A (zh) * | 2009-07-17 | 2009-12-30 | 中兴通讯股份有限公司 | 分组干扰抑制方法及装置 |
CN101958729A (zh) * | 2010-09-30 | 2011-01-26 | 哈尔滨工程大学 | 直扩通信***中快速线性调频干扰检测和抑制方法 |
CN103929396A (zh) * | 2014-05-08 | 2014-07-16 | 西安电子科技大学 | Mimo-ofdm***下行信息数据的处理方法 |
US20140250972A1 (en) * | 2012-09-08 | 2014-09-11 | Schweitzer Engineering Laboratories, Inc. | Quality of Precision Time Sources |
CN105223590A (zh) * | 2015-10-27 | 2016-01-06 | 国信军创(岳阳)六九零六科技有限公司 | 一种混合卫星导航抗干扰方法及装置 |
CN105634543A (zh) * | 2015-12-30 | 2016-06-01 | 航天恒星科技有限公司 | 一种抗窄带干扰方法及*** |
CN105699948A (zh) * | 2015-11-27 | 2016-06-22 | 中国人民解放军理工大学 | 基于支持向量机并改进均方误差性能的波束形成方法及*** |
US20160266259A1 (en) * | 2014-08-06 | 2016-09-15 | Raytheon Company | Position And Motion Informed Navigation System |
US20170227650A1 (en) * | 2016-02-10 | 2017-08-10 | Raytheon Company | Mitigation of Spoofer Satellite Signals |
-
2017
- 2017-09-08 CN CN201710806427.4A patent/CN107728166B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1463105A (zh) * | 2003-06-18 | 2003-12-24 | 清华大学 | 用于异步码分多址***的分组单次干扰抵消方法 |
CN101114863A (zh) * | 2006-07-28 | 2008-01-30 | 美国博通公司 | 通信***内处理信号的方法和*** |
CN101615982A (zh) * | 2009-07-17 | 2009-12-30 | 中兴通讯股份有限公司 | 分组干扰抑制方法及装置 |
CN101958729A (zh) * | 2010-09-30 | 2011-01-26 | 哈尔滨工程大学 | 直扩通信***中快速线性调频干扰检测和抑制方法 |
US20140250972A1 (en) * | 2012-09-08 | 2014-09-11 | Schweitzer Engineering Laboratories, Inc. | Quality of Precision Time Sources |
CN103929396A (zh) * | 2014-05-08 | 2014-07-16 | 西安电子科技大学 | Mimo-ofdm***下行信息数据的处理方法 |
US20160266259A1 (en) * | 2014-08-06 | 2016-09-15 | Raytheon Company | Position And Motion Informed Navigation System |
CN105223590A (zh) * | 2015-10-27 | 2016-01-06 | 国信军创(岳阳)六九零六科技有限公司 | 一种混合卫星导航抗干扰方法及装置 |
CN105699948A (zh) * | 2015-11-27 | 2016-06-22 | 中国人民解放军理工大学 | 基于支持向量机并改进均方误差性能的波束形成方法及*** |
CN105634543A (zh) * | 2015-12-30 | 2016-06-01 | 航天恒星科技有限公司 | 一种抗窄带干扰方法及*** |
US20170227650A1 (en) * | 2016-02-10 | 2017-08-10 | Raytheon Company | Mitigation of Spoofer Satellite Signals |
Non-Patent Citations (3)
Title |
---|
DANESHMAND, S等: "GNSS space-time interference mitigation and attitude determination in the presence of interference signals", 《SENSORS》 * |
QIANG GUO等: "Combining DCQGMP-Based Sparse Decomposition and MPDR Beamformer for Multi-Type Interferences Mitigation for GNSS Receivers", 《SENSORS》 * |
周浩等: "应用时频分析进行高频雷达射频干扰抑制", 《电子学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110320535A (zh) * | 2019-07-03 | 2019-10-11 | 哈尔滨工程大学 | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 |
CN110320535B (zh) * | 2019-07-03 | 2023-05-05 | 哈尔滨工程大学 | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 |
CN112083448A (zh) * | 2020-09-04 | 2020-12-15 | 哈尔滨工程大学 | 面向卫星导航***的干扰信号分类识别特征提取方法及*** |
CN112083448B (zh) * | 2020-09-04 | 2023-08-15 | 哈尔滨工程大学 | 面向卫星导航***的干扰信号分类识别特征提取方法及*** |
CN114675307A (zh) * | 2022-03-14 | 2022-06-28 | 合肥工业大学 | 基于分数阶伪Wigner-Ville分布的卫星导航接收机干扰检测方法 |
CN114675307B (zh) * | 2022-03-14 | 2024-04-30 | 合肥工业大学 | 基于分数阶伪Wigner-Ville分布的卫星导航接收机干扰检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107728166B (zh) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104536017B (zh) | 一种先子空间投影后波束合成的导航接收机stap方法 | |
CN101821646B (zh) | 判断涉及固定或移动发射机、反射器和传感器的传播信道的脉冲响应的方法和装置 | |
CN101807977B (zh) | 基于波形特征的空时盲自适应抗干扰方法 | |
CN101718873A (zh) | 导航信号空时联合抗干扰数字信号处理器 | |
CN101799551B (zh) | 基于解重扩技术的空时盲自适应gps干扰抑制方法 | |
CN104865586A (zh) | 一种阵列天线导航接收机抗干扰及多径抑制方法 | |
CN104536018B (zh) | 一种使用阵列天线抗干扰技术的gnss多星联合捕获方法 | |
Wu et al. | Adaptive interference mitigation in GNSS | |
CN105866750A (zh) | 一种多个gps卫星微弱回波信号检测的方法 | |
CN102798870B (zh) | 一种卫星导航***兼容性分析方法 | |
CN107272030B (zh) | 一种基于多个gps卫星信息融合的目标探测方法 | |
CN107728166A (zh) | 一种基于时域分组处理的卫星导航接收机多干扰抑制方法 | |
CN103105614B (zh) | 基于惯导辅助的空时域联合抗干扰方法 | |
CN104155662A (zh) | 基于gnss相关峰值探测器的自适应互干扰抑制方法 | |
Bao et al. | Spoofing mitigation in Global Positioning System based on C/A code self-coherence with array signal processing | |
CN111198387A (zh) | 一种抗欺骗干扰的空时采样导航定位方法 | |
Qiao et al. | A survey of GNSS interference monitoring technologies | |
Öktem et al. | Power delay doppler profile fingerprinting for mobile localization in NLOS | |
Lv et al. | General method to mitigate the continuous wave interference and narrowband interference for GNSS receivers | |
Enneking et al. | Exploiting WSSUS multipath for GNSS ranging | |
Anantharamu et al. | Space-time equalization techniques for new GNSS signals | |
CN106338742A (zh) | 基于互谱准则的降维自适应多波束gps信号抗干扰方法 | |
Karasek et al. | The DVB-T-based positioning system and single frequency network offset estimation | |
Brahim et al. | Precise acquisition of global navigation satellite system signals in the presence of multipath and influence on tracking capability | |
Wen et al. | A flexible GNSS spoofer localization system: Spoofing discrimination and localization method |
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 |