CN104049262A - 一种基于矢量跟踪环路的波束形成抗干扰方法 - Google Patents
一种基于矢量跟踪环路的波束形成抗干扰方法 Download PDFInfo
- Publication number
- CN104049262A CN104049262A CN201410298930.XA CN201410298930A CN104049262A CN 104049262 A CN104049262 A CN 104049262A CN 201410298930 A CN201410298930 A CN 201410298930A CN 104049262 A CN104049262 A CN 104049262A
- Authority
- CN
- China
- Prior art keywords
- signal
- carrier
- vector
- satellite
- theta
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 5
- 238000005311 autocorrelation function Methods 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000005764 inhibitory process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000116 mitigating effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及阵列天线和全球定位***接收机矢量跟踪环路领域,具体的涉及的是一种为基于阵列天线的波束形成器提供本地参考信号,保证矢量跟踪环路正常工作的基于矢量跟踪环路的波束形成抗干扰方法。本发明包括;阵列天线接收信号;剥离卫星信号中载波;求取阵列最优权值;计算当前各跟踪通道伪距测量值,获得载***置;构建矢量跟踪环路。本发明够抑制与卫星信号不相干的各种干扰,属于盲波束形成的一种。同时,矢量跟踪环路能够容易跟踪抑制干扰后的卫星信号,不易失锁,从而保证GPS接收机的在干扰环境中的稳健性。
Description
技术领域
本发明涉及阵列天线和全球定位***接收机矢量跟踪环路领域,具体的涉及的是一种为基于阵列天线的波束形成器提供本地参考信号,保证矢量跟踪环路正常工作的基于矢量跟踪环路的波束形成抗干扰方法。
背景技术
GPS卫星信号到达地球表面时的信号功率约为-160dBW,因此GPS极易受到干扰信号的影响。GPS信号依靠其扩频体制能够为其获得43dB的处理增益,然而随着人为干扰技术的提高及干扰设备成本的降低,卫星信号只依靠其扩频体制进行抗干扰已经不能满足用户需求。另外,现有GPS信号是在众所周知的频率上发射的,其调制特征广为人知,信噪比又比较低,因而很容易进行干扰或欺骗。因此,研究GPS抗干扰技术已成为热点和重点。
目前,国内外抗干扰抑制技术主要分为两类:第一类是前相关抗干扰技术
(Pre-correlation mitigation techniques),主要是在信号进入接收机通道之前完成
干扰信号的监测与抑制,如基于阵列天线的空间调零和波束形成技术。第二类就是后相关抗干扰技术(Post-correlation mitigation techniques),主要针对接收机通道中的码环和载波环,具体方法有自适应跟踪环路技术、矢量跟踪环路技术和外部惯性测量单元(IMU)辅助环路技术。
自适应天线是将天线阵列与信号处理相结合,用信号处理的理论和方法、自动控制的技术调准天线阵的方向图,使天线方向图主波束对准有用信号方向,零点实时指向干扰方向。从而使它在干扰方向形成零陷,将干扰信号抵消,还可以使有用信号得到加强,从而达到抗干扰的目的。信号和干扰变化时,信号处理装置按照设计的自适应准则和算法及时调整权集,使波束自动、实时地跟踪这种变化。天线阵在接收到信号后,通过对天线阵的每个阵元上引入可以调整的权系数,来实现对天线阵主波束方向的灵活控制,达到对整个空域信号的接收,这个天线阵元上的权系数的调整也即实现了“空域滤波”功能。用这种方法产生的天线阵主波束明显比全向天线的主波束窄,并且主波束的宽度与天线阵元的数目有关,天线阵元数目越多,则形成的主波束也就越窄,反之越宽。在形成窄的主波束同时也最大限度的抑制了与来波方向不一致的干扰,使得接收机中信号的信干比(信号干扰功率比)有明显的改善,从而增强了抗干扰性能。
在矢量跟踪环路中,接收机各个通道的信息通过滤波器融合在一起,不用考虑每个跟踪通道滤波器设计。传统GPS接收机利用接收信号的相位和频率估计出接收机的位置和速度,矢量跟踪环路接收机逆向使用了传统接收机的基本工作原理:从载体的位置、速度等导航解算信息估计出接收信号的C/A码码相位、载波相位、频率等参数信息,这些预测信息与接收机接收到的下一时刻信息在每个通道产生差值,利用这些差值去更新校正下一时刻载体的伪距、伪距率。另外,当存在干扰时,卫星信号的等效载噪比将会下降。而在矢量跟踪环路中,载噪比较强的跟踪通道通过导航滤波器可以辅助载噪比较弱的通道对卫星信号进行跟踪,从而增强接收机在干扰情况下的稳健性。
然而,自适应阵列天线抗干扰方法和矢量跟踪环路抗干扰方法都存在一定的缺陷。在基于阵列天线的波束形成中,常用的求取最优权矢量准则有最小方差无失真响应准则(MVDR)和最小均方误差准则(MMSE)。MVDR方法需要知道期望信号的空间角度,在实际中,由于GPS信号功率很弱,其到达阵列天线处的角度不容易估计出来。而MMSE方法则不需要知道卫星信号的来波方向,但是需要知道卫星信号的本地参考信号。本发明采用基于最小均方误差准则的自相干恢复(self-coherence restoral,SCORE)算法来求取阵列天线的自适应权矢量,由接收机矢量跟踪环路为SCORE波束形成器提供本地参考信号。同时波束形成抑制干扰后的卫星信号进入矢量跟踪环路进行信号处理,从而保证矢量跟踪环路能够正常工作。
发明内容
本发明的目的在于提供接收机在干扰环境中的稳健性的一种基于矢量跟踪环路的波束形成抗干扰方法。
本发明的目的是这样实现的:
(1)阵列天线接收信号:
天线阵元数M,相邻阵元的间隔为λ/2,λ是GPS信号的波长,期望信号和干扰信号的到达角度分别为θ0和θk,k=1、2,…P,则阵列天线接收信号为:
s0(t)为卫星直达信号,sk(t)为第k个干扰信号,n(t)为空间噪声为相互独立的零均值高斯白噪声,n(t)=[n0(t),n1(t),…,np(t)]T,a(θ0)为卫星信号的导向矢量,
a(θk)为第k个干扰信号的导向矢量,
GPS信号s0(t)为由伪随机码、L1载波和导航数据电文组成,
s0(t)=AC(t)D(t)cos(ωt+φ0)
A为信号幅值,C(t)为C/A码,D(t)为导航电文数据,ω为L1载波角频率,φ0为L1载波频段初始载波相位;
(2)剥离卫星信号中载波:
在跟踪环路的混频器中,GPS信号分别与本地载波数控振荡器产生的同相信号cos(ωt+φ)、正交信号sin(ωt+φ)相乘得到同向支路信号I(t)和正交支路信号Q(t),φ为本地复制载波初始相位,I(t)和Q(t)两路输出的信号经过低通滤波器后,信号中的高频分量被滤除,当跟踪环路跟踪上卫星信号的载波初始相位时,即φ0=φ,卫星信号中的载波被完全剥离;
(3)求取阵列最优权值:
干扰和噪声与期望信号不相关,干扰信号为独立的零均值高斯白噪声,阵列天线接收信号为
x(t)=a(θ0)s0(t)+v(t)
v(t)为干扰信号加噪声,
GPS信号s0(t)采用BPSK调制方式,循环自相干函数在频移α处不为零,即v(t)为独立的零均值高斯白噪声,所以循环自相干函数在频移α处为零,对于时延τ,天线阵列接收信号x(t)的循环自相关函数为
τ为本地复制C/A码相位的延迟,本地载波已经与卫星信号中的载波同步,频移α=0;
d(t)为矢量跟踪环路需要提供的本地参考信号,误差信号ε(t)为
ε(t)=d(t)-wHx(t)
SCORE波束形成器的代价函数为
Rxd=E{x(t)d*(t)}表示阵列天线接收数据与参考信号的互相关矩阵,Rxx=E{x(t)x(t)H}表示阵列接收数据的协方差矩阵,f为控制矢量,Rdd为参考信号的协方差矩阵,
最优权矢量w为最大特征值对应的特征矢量:
λ为最大特征值,最优权矢量w为与λ相对应的特征向量;
(4)计算当前各跟踪通道伪距测量值,获得载***置:
第n个跟踪通道量测的伪距ρn为
ρn=(τ0+τ)·c/fc/a
τ0为本地NCO中的伪码相位基准值,τ仍为本地复制C/A码相位的延迟,c为光速,fc/a为码片速率1.023MHz;
在k时刻,从跟踪通道内获得n个伪距测量值,第n颗卫星的位置坐标可由星历计算得到,计算出载体k时刻的位置坐标Pk=[x,y,z]T;
(5)构建矢量跟踪环路
采用矢量跟踪环路为波束形成提供参考信号,获得k时刻的载***置Pk后,估计出载体在k+1时刻的位置Pk+1,利用卫星星历预报出k+1时刻第n颗卫星的位置由Pk+1和计算出K+1时刻载体与第n颗卫星间的伪距预测的伪距反馈到本地NCO产生本地复制信号,反推出下一时刻码相位值,码相位预测值经过码相关器为波束形成器提供参考信号,码相关器获得预测码相位误差后,获得K+1时刻伪距误差与预测的伪距值叠加后得到K+1时刻准确的伪距值
用更新载***置,完成闭合跟踪环路。
本发明的有益效果在于:
在许多阵列天线抗干扰方法中,往往需要知道期望信号的导向矢量,即需要知道期望信号到达阵列天线的空间角度。而在基于阵列天线的GPS抗干扰中,由于GPS信号功率很弱,GPS信号的空间角度不容易被估计出来。另一方面,矢量跟踪环路也有着一定的抗干扰能力,但是抗干扰能力较弱。
本发明采用阵列天线与矢量跟踪环路相结合的方法来提高GPS接收机的抗干扰能力,本发明采用基于最小均方误差准则的SCORE算法来求取阵列天线的自适应权矢量。该方法不需要知道卫星信号的来波方向,而是根据卫星信号的特点,由矢量跟踪环路为基于阵列天线的SCORE波束形成器提供本地参考信号。该方法能够抑制与卫星信号不相干的各种干扰,属于盲波束形成的一种。同时,矢量跟踪环路能够容易跟踪抑制干扰后的卫星信号,不易失锁,从而保证GPS接收机的在干扰环境中的稳健性。
附图说明
图1是矢量跟踪环路结构框图;
图2是基于最小均方误差准则的自适应波束形成示意图;
图3是基于矢量跟踪环路的波束形成流程图。
具体实施方式
下面结合说明书附图对本发明做进一步的阐述:
主要利用矢量跟踪环路为基于阵列天线的波束形成器提供本地参考信号,然后求取阵列的自适应权值来对干扰信号进行抑制。抑制干扰后的卫星信号进入矢量跟踪环路来对伪码进行实时跟踪,从而保证接收机在干扰环境中的稳健性。
本发明提出的一种基于矢量跟踪环路的波束形成抗干扰方法,具体步骤如下:
步骤1:建立阵列天线接收信号模型
假设有M个天线阵元,相邻阵元的间隔为λ/2,λ是GPS信号的波长。假设期望信号和干扰信号的到达角度分别为θ0和θk(k=1、2,…P),则阵列天线接收信号可以表示为
其中s0(t)代表卫星直达信号,sk(t)为第k个干扰信号,n(t)代表空间噪声为相互独立的零均值高斯白噪声,n(t)可以表示为n(t)=[n0(t),n1(t),…,np(t)]T。a(θ0)为卫星信号的导向矢量,表示为a(θk)为第k个干扰信号的导向矢量,表示为
在公式(1)中,GPS信号s0(t)为由伪随机码(C/A码)、L1载波(载波频率为1575.42MHz得余弦波)和导航数据电文组成,可以表示为
s0(t)=AC(t)D(t)cos(ωt+φ0) (2)
式(2)中,A为信号幅值,C(t)表示C/A码,D(t)表示导航电文数据,ω为L1载波角频率,φ0为L1频段初始载波相位。
步骤2:剥离卫星信号中载波
在跟踪环路的混频器中,式(2)中的GPS信号分别与本地载波数控振荡器(NCO)产生的同相信号cos(ωt+φ)、正交信号sin(ωt+φ)相乘,φ为本地复制载波初始相位,得到的同向支路I(t)和正交支路Q(t)两路输出信号。I(t)和Q(t)两路输出的信号经过低通滤波器后,信号中的高频分量被滤除。当跟踪环路跟踪上卫星信号的载波初始相位时,即φ0=φ,卫星信号中的载波被完全剥离。
步骤3:建立SCORE波束形成器算法模型,求取阵列最优权值
本发明将基于MMSE准则的波束形成器与矢量跟踪环路相结合,其中,矢量跟踪环路指的是伪码矢量跟踪环路。由于卫星信号中的载波已经被载波混频器所剥离,因此,本地参考信号只需提供伪码参考信号。
本发明所提出的抗干扰方法是在自相干恢复(self-coherence restoral,SCORE)算法的基础上改进得到的。对于一个信号s(t)来说,如果对于某个时延τ,s(t)与其频移α后所得信号的相关值不为0,则称s(t)具有谱相干特性,即
式(3)中,(·)*表示为复共轭,<·>∞表示无限时间平均运算,为自相干函数,Rss(0)为s(t)的平均功率,为s(t)的循环自相干函数。和Rss(0)分别为
在本发明中,假设干扰和噪声与期望信号不相关。为简化方便,将干扰信号也假设为独立的零均值高斯白噪声。因此,公式(1)可以简写为
x(t)=a(θ0)s0(t)+v(t) (6)
式(6)中,v(t)为干扰信号加噪声,仍未独立的零均值高斯白噪声。
由于GPS信号s0(t)采用的是BPSK调制方式,其循环自相干函数在频移α处不为零,即由于v(t)为独立的零均值高斯白噪声,所以其循环自相干函数在频移α处为零。因此,对于某个时延τ,天线阵列接收信号x(t)的循环自相关函数可以表示为
本发明中,τ为本地复制C/A码相位的延迟。由于波束形成器位于跟踪环路中载波解调过程之后,这里认为本地载波已经与卫星信号中的载波同步,所以步骤3中频移α=0。
本发明所提出的SCORE波束形成器算法是基于最小均方误差准则(MMSE)的。MMSE使参考信号d(t)与阵列输出y(t)=wHx(t)的均方误差最小,利用参考信号求解自适应权矢量w。在本发明中,d(t)即为矢量跟踪环路需要提供的本地参考信号。误差信号ε(t)可以表示为
ε(t)=d(t)-wHx(t) (8)
SCORE波束形成器的代价函数可以写为
式(9)中,Rxd=E{x(t)d*(t)}表示阵列天线接收数据与参考信号的互相关矩阵,Rxx=E{x(t)x(t)H}表示阵列接收数据的协方差矩阵,f为控制矢量,Rdd为参考信号的协方差矩阵。公式(9)的目的即求取自适应权值w使得该代价函数取得最大值。
通过求解,最优权矢量w为下面广义特征值问题中最大特征值对应的特征矢量:
式(10)中,λ为最大特征值,最优权矢量w为与λ相对应的特征向量。
步骤4:计算当前各跟踪通道伪距测量值,获得载***置
本发明是根据MMSE准则,通过使参考信号d(t)与阵列输出y(t)的均方误差最小来抑制干扰的。当抑制干扰后,矢量跟踪环路即可对卫星信号中的码相位进行跟踪。第n个跟踪通道量测的伪距ρn可以表示为
ρn=(τ0+τ)·c/fc/a (11)
式(11)中,τ0为本地NCO中的伪码相位基准值,τ仍为本地复制C/A码相位的延迟,c为光速,fc/a为码片速率1.023MHz。
在k时刻,若矢量跟踪环路跟踪上n(n≥4)颗卫星,即可从跟踪通道内获得n个伪距测量值。另外,第n颗卫星的位置坐标X(n)=[x(n),y(n),z(n)]T可由星历计算得到。然后根据GPS定位原理,可以计算出载体k时刻的位置坐标Pk=[x,y,z]T
步骤5:构建矢量跟踪环路
而在本发明中,采用矢量跟踪环路为波束形成提供参考信号,附图3为矢量跟踪环路结构框图。在步骤4中,获得k时刻的载***置Pk后,根据载体运动模型,可以估计出载体在k+1时刻的位置Pk+1。同时,利用卫星星历预报出k+1时刻第n颗卫星的位置由Pk+1和可以计算出K+1时刻载体与第n颗卫星间的伪距这些预测的伪距反馈到本地NCO产生本地复制信号,即再次利用公式(11)反推出下一时刻码相位值。这些码相位预测值经过码相关器为波束形成器提供参考信号。另外,码相关器获得预测码相位误差后,获得K+1时刻伪距误差然后与预测的伪距值叠加后得到K+1时刻准确的伪距值
然后,用来更新载***置,从而完成闭合跟踪环路
步骤1:建立阵列天线接收信号模型
本发明考虑均匀线性阵列天线结构,假设有M个天线阵元,相邻阵元的间隔为λ/2,λ是GPS信号的波长。为表示方便,这里假设存在一个卫星信号,P个干扰信号,且干扰信号、噪声均与卫星信号不相干。
矢量跟踪环路为基于阵列天线的波束形成器提供本地参考信号,然后采用采用基于最小均方误差准则的自相干恢复(self-coherence restoral,SCORE)算法来求取阵列天线的自适应权矢量。抑制干扰后的卫星信号经过积分清除器和鉴别器得到各通道的伪距测量值,利用伪距计算载体当前位置,结合当前卫星位置预测下一时刻载体与卫星间的伪距。预测的伪距反馈到伪码数控振荡器(NCO)产生本地复制信号,经过相关器为波束形成器提供参考信号,同时计算下一时刻伪距误差,然后与预测的伪距值叠加后得到下一时刻准确的伪距值,完成闭合跟踪环路。
阵列天线接收到的信号经过载波混频器,卫星信号中的载波已经被剥离。矢量跟踪环路只考虑对卫星信号中的伪码进行跟踪和解算,因此,本地参考信号也只含伪码信息。
假设期望信号和干扰信号的到达角度分别为θ0和θk(k=1、2,…P),则阵列天线接收信号可以表示为
其中s0(t)代表GPS直达信号,sk(t)为第k个干扰信号,n(t)代表空间噪声为相互独立的零均值高斯白噪声,n(t)可以表示为n(t)=[n0(t),n1(t),…,np(t)]T。a(θ0)为卫星信号的导向矢量,表示为a(θk)为第k个干扰信号的导向矢量,表示为
在公式(1)中,GPS信号s0(t)为由伪随机码(C/A码)、L1载波(载波频率为1575.42MHz得余弦波)和导航数据电文组成,可以表示为
s0(t)=AC(t)D(t)cos(ωt+φ0) (2)
式(2)中,A为信号幅值,C(t)表示C/A码,D(t)表示导航电文数据,ω为L1载波角频率,φ0为L1频段初始载波相位。
步骤2:卫星信号中载波被剥离
在载波跟踪过程中,可将式(2)中的GPS信号分别与本地载波数控振荡器(NCO)产生的同相信号cos(ωt+φ)、正交信号sin(ωt+φ)相乘(即混频器中的过程),φ为本地复制载波初始相位,得到的同向支路I(t)和正交支路Q(t)两路输出信号,分别可表示为:
经过低通滤波后,式(3)和(4)中的高频分量被滤除,则I(t)和Q(t)两路的输出可表示为:
当跟踪环路跟踪上卫星信号的载波初始相位时,即φ0=φ,I(t)和Q(t)两路的输出为
Q(t)=0 (8)
此时,卫星信号中的载波被完全剥离。
步骤3:建立SCORE波束形成器算法模型,求取阵列最优权值
本发明将基于MMSE准则的波束形成器与矢量跟踪环路相结合,其中,矢量跟踪环路指的是伪码矢量跟踪环路。由于卫星信号中的载波已经被载波混频器所剥离,因此,本地参考信号只需提供伪码参考信号。附图1即为基于矢量跟踪环路的波束形成流程图。
本发明所提出的抗干扰方法是在自相干恢复(self-coherence restoral,SCORE)算法的基础上改进得到的。对于一个信号s(t)来说,如果对于某个时延τ,s(t)与其频移α后所得信号的相关值不为0,则称s(t)具有谱相干特性,即
式(3)中,(·)*表示为复共轭,<·>∞表示无限时间平均运算,为自相干函数,Rss(0)为s(t)的平均功率,为s(t)的循环自相干函数。和Rss(0)分别为
在本发明中,假设干扰和噪声与期望信号不相关。为简化方便,将干扰信号也假设为独立的零均值高斯白噪声。因此,公式(1)可以简写为
x(t)=a(θ0)s0(t)+v(t) (6)
式(6)中,v(t)为干扰信号加噪声,仍未独立的零均值高斯白噪声。
由于GPS信号s0(t)采用的是BPSK调制方式,其循环自相干函数在频移α处不为零,即由于v(t)为独立的零均值高斯白噪声,所以其循环自相干函数在频移α处为零。因此,对于某个时延τ,天线阵列接收信号x(t)的循环自相关函数可以表示为
本发明中,τ为本地复制C/A码相位的延迟。由于波束形成器位于跟踪环路中载波解调过程之后,这里认为本地载波已经与卫星信号中的载波同步,所以步骤3中频移α=0。
本发明所提出的SCORE波束形成器算法是基于最小均方误差准则(MMSE)的。MMSE使参考信号d(t)与阵列输出y(t)=wHx(t)的均方误差最小,利用参考信号求解自适应权矢量w,如附图2所示。在本发明中,d(t)即为矢量跟踪环路需要提供的本地参考信号。误差信号ε(t)可以表示为
ε(t)=d(t)-wHx(t) (8)
对上式两边取模的平方并且求数学期望,可得其均方误差为
E{|ε(t)|2}=E{[d(t)-wHx(t)][d*(t)-x(t)Hw]}
(9)
=E{d2(t)}-2wHRxd+wHRxxw
其中,表示阵列天线接收数据与参考信号的互相关矩阵,Rxx=E{x(t)x(t)H}表示阵列接收数据的协方差矩阵。为使均方误差最小,将式(9)对w求导并令其为0,可得
2Rxxw-2Rxd=0 (10)
则MMSE的最优权矢量为
此最优权矢量被称为维纳解。
SCORE波束形成器的代价函数可以写为
式(12)中,Rxd和Rxx已经在公式(9)下方定义,f为控制矢量,Rdd为参考信号的协方差矩阵。公式(12)的目的即求取自适应权值w使得该代价函数取得最大值。
通过求解,最优权矢量w为下面广义特征值问题中最大特征值对应的特征矢量:
式(13)中,λ为最大特征值,最优权矢量w为与λ相对应的特征向量。
步骤4:计算当前各跟踪通道伪距测量值,获得载***置
本发明是根据MMSE准则,通过使参考信号d(t)与阵列输出y(t)的均方误差最小来抑制干扰的。当抑制干扰后,矢量跟踪环路即可对卫星信号中的码相位进行跟踪。第n个跟踪通道量测的伪距ρn可以表示为
ρn=(τ0+τ)·c/fc/a (14)
式(14)中,τ0为本地NCO中的伪码相位基准值,τ仍为本地复制C/A码相位的延迟,c为光速,fc/a为码片速率1.023MHz。
在k时刻,若矢量跟踪环路跟踪上n颗卫星,第n颗卫星的位置坐标X(n)=[x(n),y(n),z(n)]T可由星历计算得到。因此,载体在k时刻的位置坐标Pk=[x,y,z]T可以通过式(15)计算得到
式(15)中,δtu为时钟误差。从式(15)可以看出,如果接收机跟踪上4颗卫星即可解算出载***置。
步骤5:构建矢量跟踪环路
在传统的GPS接收机伪码跟踪环路中,GPS信号经过射频前段下变频后变为中频信号,各个跟踪通道的基带信号模块对中频信号进行处理,得到各通道的伪距误差。这些伪距误差一方面作为误差校正量反馈给本地信号发生器,对复制信号的码相位进行校正,一方面输入到接收机的导航解算模块进行位置解算。
而在本发明中,采用矢量跟踪环路为波束形成提供参考信号,附图3为矢量跟踪环路结构框图。在步骤4中,通过式(15)获得k时刻的载***置Pk后,根据载体运动模型,可以估计出载体在k+1时刻的位置Pk+1。
以载体做匀速直线运动为例,载体运动速度为V,Pk+1可以表示为
Pk+1=Pk+VT (16)
式中,T为导航滤波器更新周期。同时,利用卫星星历中的轨道升交点赤经、轨道倾角、近地点角距、长半径、偏心率和卫星的真近点角等参数预报出k+1时刻第n颗卫星的位置由Pk+1和可以计算出K+1时刻载体与第n颗卫星间的伪距这些预测的伪距反馈到本地NCO产生本地复制信号,即再次利用公式(14)反推出下一时刻码相位值。这些码相位预测值经过码相关器为波束形成器提供参考信号。另外,码相关器获得预测码相位误差后,获得K+1时刻伪距误差然后与预测的伪距值叠加后得到K+1时刻准确的伪距值
然后,再次代入公式(15)来更新载***置,从而完成闭合跟踪环路。
Claims (1)
1.一种基于矢量跟踪环路的波束形成抗干扰方法,其特征在于:
(1)阵列天线接收信号:
天线阵元数M,相邻阵元的间隔为λ/2,λ是GPS信号的波长,期望信号和干扰信号的到达角度分别为θ0和θk,k=1、2,…P,则阵列天线接收信号为:
s0(t)为卫星直达信号,sk(t)为第k个干扰信号,n(t)为空间噪声为相互独立的零均值高斯白噪声,n(t)=[n0(t),n1(t),…,np(t)]T,a(θ0)为卫星信号的导向矢量,
a(θk)为第k个干扰信号的导向矢量,
GPS信号s0(t)为由伪随机码、L1载波和导航数据电文组成,
s0(t)=AC(t)D(t)cos(ωt+φ0)
A为信号幅值,C(t)为C/A码,D(t)为导航电文数据,ω为L1载波角频率,φ0为L1载波频段初始载波相位;
(2)剥离卫星信号中载波:
在跟踪环路的混频器中,GPS信号分别与本地载波数控振荡器产生的同相信号cos(ωt+φ)、正交信号sin(ωt+φ)相乘得到同向支路信号I(t)和正交支路信号Q(t),φ为本地复制载波初始相位,I(t)和Q(t)两路输出的信号经过低通滤波器后,信号中的高频分量被滤除,当跟踪环路跟踪上卫星信号的载波初始相位时,即φ0=φ,卫星信号中的载波被完全剥离;
(3)求取阵列最优权值:
干扰和噪声与期望信号不相关,干扰信号为独立的零均值高斯白噪声,阵列天线接收信号为
x(t)=a(θ0)s0(t)+v(t)
v(t)为干扰信号加噪声,
GPS信号s0(t)采用BPSK调制方式,循环自相干函数在频移α处不为零,即v(t)为独立的零均值高斯白噪声,所以循环自相干函数在频移α处为零,对于时延τ,天线阵列接收信号x(t)的循环自相关函数为
τ为本地复制C/A码相位的延迟,本地载波已经与卫星信号中的载波同步,频移α=0;
d(t)为矢量跟踪环路需要提供的本地参考信号,误差信号ε(t)为
ε(t)=d(t)-wHx(t)
SCORE波束形成器的代价函数为
Rxd=E{x(t)d*(t)}表示阵列天线接收数据与参考信号的互相关矩阵,Rxx=E{x(t)x(t)H}表示阵列接收数据的协方差矩阵,f为控制矢量,Rdd为参考信号的协方差矩阵,
最优权矢量w为最大特征值对应的特征矢量:
λ为最大特征值,最优权矢量w为与λ相对应的特征向量;
(4)计算当前各跟踪通道伪距测量值,获得载***置:
第n个跟踪通道量测的伪距ρn为
ρn=(τ0+τ)·c/fc/a
τ0为本地NCO中的伪码相位基准值,τ仍为本地复制C/A码相位的延迟,c为光速,fc/a为码片速率1.023MHz;
在k时刻,从跟踪通道内获得n个伪距测量值,第n颗卫星的位置坐标可由星历计算得到,计算出载体k时刻的位置坐标Pk=[x,y,z]T;
(5)构建矢量跟踪环路:
采用矢量跟踪环路为波束形成提供参考信号,获得k时刻的载***置Pk后,估计出载体在k+1时刻的位置Pk+1,利用卫星星历预报出k+1时刻第n颗卫星的位置由Pk+1和计算出K+1时刻载体与第n颗卫星间的伪距预测的伪距反馈到本地NCO产生本地复制信号,反推出下一时刻码相位值,码相位预测值经过码相关器为波束形成器提供参考信号,码相关器获得预测码相位误差后,获得K+1时刻伪距误差与预测的伪距值叠加后得到K+1时刻准确的伪距值
用更新载***置,完成闭合跟踪环路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298930.XA CN104049262B (zh) | 2014-06-26 | 2014-06-26 | 一种基于矢量跟踪环路的波束形成抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298930.XA CN104049262B (zh) | 2014-06-26 | 2014-06-26 | 一种基于矢量跟踪环路的波束形成抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104049262A true CN104049262A (zh) | 2014-09-17 |
CN104049262B CN104049262B (zh) | 2017-01-25 |
Family
ID=51502348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410298930.XA Expired - Fee Related CN104049262B (zh) | 2014-06-26 | 2014-06-26 | 一种基于矢量跟踪环路的波束形成抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104049262B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536017A (zh) * | 2015-01-06 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种先子空间投影后波束合成的导航接收机stap算法 |
CN105954772A (zh) * | 2016-05-16 | 2016-09-21 | 中国人民解放军国防科学技术大学 | 一种稳健无偏的导航信号矢量跟踪方法 |
CN106226753A (zh) * | 2016-08-22 | 2016-12-14 | 西安电子科技大学 | 基于最小方差法谱函数二阶导数的波达方向估计方法 |
CN106483536A (zh) * | 2015-08-26 | 2017-03-08 | 罗克韦尔柯林斯公司 | 利用具备波束赋形能力的天线的卫星信号采集 |
CN107045131A (zh) * | 2017-05-24 | 2017-08-15 | 西北工业大学 | 张量分解卫星导航抗干扰方法 |
CN107179540A (zh) * | 2017-05-26 | 2017-09-19 | 中国民航大学 | 基于解重扩算法的gnss矢量接收机抗干扰方法 |
CN108811078A (zh) * | 2017-05-04 | 2018-11-13 | 三星电子株式会社 | 用于5g-nr***的多波束初始同步的装置和方法 |
CN109298435A (zh) * | 2017-11-02 | 2019-02-01 | 天津讯联科技有限公司 | 卫星反射信号的跟踪方法、装置及接收机 |
CN113281701A (zh) * | 2021-04-28 | 2021-08-20 | 中国人民解放军战略支援部队信息工程大学 | 协同短波多站角度与三星时差的超视距目标直接定位方法 |
CN113534198A (zh) * | 2021-06-16 | 2021-10-22 | 北京遥感设备研究所 | 基于协方差矩阵重构的卫星导航动态抗干扰方法及其*** |
CN113866796A (zh) * | 2021-09-10 | 2021-12-31 | 西安电子科技大学 | 基于空域滤波与波束成形的导航接收机抗干扰方法 |
CN115037331A (zh) * | 2022-08-10 | 2022-09-09 | 中国电子科技集团公司第十研究所 | 一种基于逆向外推的异步突发信号定时同步方法 |
CN116035558A (zh) * | 2023-03-02 | 2023-05-02 | 中国科学技术大学 | 基于波束形成的抗干扰呼吸检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120326926A1 (en) * | 2011-06-24 | 2012-12-27 | Mayflower Communications Company, Inc. | High sensitivity gps/gnss receiver |
CN103116169A (zh) * | 2013-01-20 | 2013-05-22 | 哈尔滨工程大学 | 一种基于矢量跟踪环路的抗干扰方法 |
-
2014
- 2014-06-26 CN CN201410298930.XA patent/CN104049262B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120326926A1 (en) * | 2011-06-24 | 2012-12-27 | Mayflower Communications Company, Inc. | High sensitivity gps/gnss receiver |
CN103116169A (zh) * | 2013-01-20 | 2013-05-22 | 哈尔滨工程大学 | 一种基于矢量跟踪环路的抗干扰方法 |
Non-Patent Citations (3)
Title |
---|
卢丹: "稳健的全球卫星导航***抗干扰技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
李平: "基于周期平稳特性的GPS自适应抗干扰技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
***: "GNSS/INS深组合跟踪环路设计及滤波方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536017A (zh) * | 2015-01-06 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种先子空间投影后波束合成的导航接收机stap算法 |
CN106483536B (zh) * | 2015-08-26 | 2019-04-02 | 罗克韦尔柯林斯公司 | 利用具备波束赋形能力的天线的卫星信号采集 |
CN106483536A (zh) * | 2015-08-26 | 2017-03-08 | 罗克韦尔柯林斯公司 | 利用具备波束赋形能力的天线的卫星信号采集 |
CN105954772A (zh) * | 2016-05-16 | 2016-09-21 | 中国人民解放军国防科学技术大学 | 一种稳健无偏的导航信号矢量跟踪方法 |
CN105954772B (zh) * | 2016-05-16 | 2018-04-24 | 中国人民解放军国防科学技术大学 | 一种稳健无偏的导航信号矢量跟踪方法 |
CN106226753A (zh) * | 2016-08-22 | 2016-12-14 | 西安电子科技大学 | 基于最小方差法谱函数二阶导数的波达方向估计方法 |
CN108811078B (zh) * | 2017-05-04 | 2021-11-09 | 三星电子株式会社 | 用于5g-nr***的多波束初始同步的装置和方法 |
CN108811078A (zh) * | 2017-05-04 | 2018-11-13 | 三星电子株式会社 | 用于5g-nr***的多波束初始同步的装置和方法 |
CN107045131A (zh) * | 2017-05-24 | 2017-08-15 | 西北工业大学 | 张量分解卫星导航抗干扰方法 |
CN107179540B (zh) * | 2017-05-26 | 2020-02-18 | 中国民航大学 | 基于解重扩算法的gnss矢量接收机抗干扰方法 |
CN107179540A (zh) * | 2017-05-26 | 2017-09-19 | 中国民航大学 | 基于解重扩算法的gnss矢量接收机抗干扰方法 |
CN109298435A (zh) * | 2017-11-02 | 2019-02-01 | 天津讯联科技有限公司 | 卫星反射信号的跟踪方法、装置及接收机 |
CN109298435B (zh) * | 2017-11-02 | 2023-05-23 | 天津讯联科技有限公司 | 卫星反射信号的跟踪方法、装置及接收机 |
CN113281701A (zh) * | 2021-04-28 | 2021-08-20 | 中国人民解放军战略支援部队信息工程大学 | 协同短波多站角度与三星时差的超视距目标直接定位方法 |
CN113281701B (zh) * | 2021-04-28 | 2024-03-12 | 中国人民解放军战略支援部队信息工程大学 | 协同短波多站角度与三星时差的超视距目标直接定位方法 |
CN113534198A (zh) * | 2021-06-16 | 2021-10-22 | 北京遥感设备研究所 | 基于协方差矩阵重构的卫星导航动态抗干扰方法及其*** |
CN113534198B (zh) * | 2021-06-16 | 2023-05-23 | 北京遥感设备研究所 | 基于协方差矩阵重构的卫星导航动态抗干扰方法及其*** |
CN113866796A (zh) * | 2021-09-10 | 2021-12-31 | 西安电子科技大学 | 基于空域滤波与波束成形的导航接收机抗干扰方法 |
CN113866796B (zh) * | 2021-09-10 | 2024-04-19 | 西安电子科技大学 | 基于空域滤波与波束成形的导航接收机抗干扰方法 |
CN115037331A (zh) * | 2022-08-10 | 2022-09-09 | 中国电子科技集团公司第十研究所 | 一种基于逆向外推的异步突发信号定时同步方法 |
CN116035558A (zh) * | 2023-03-02 | 2023-05-02 | 中国科学技术大学 | 基于波束形成的抗干扰呼吸检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104049262B (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104049262B (zh) | 一种基于矢量跟踪环路的波束形成抗干扰方法 | |
US10330791B2 (en) | Method for locating a jamming source jamming signals of a satellite navigation system and associated system | |
US7471241B1 (en) | Global navigation satellite system (GNSS) receivers based on satellite signal channel impulse response | |
EP2049915B1 (en) | A global position system ( gps) user receiver and geometric surface processing for all-in-view coherent gps signal pseudo-random noise (prn) codes acquisition and navigation solution determination | |
Li et al. | A robust anti-jamming navigation receiver with antenna array and GPS/SINS | |
CN106908817B (zh) | 辅助导航定位方法及*** | |
CN102353970A (zh) | 一种高抗干扰性能gps/sins组合导航***及实现方法 | |
CN104199052A (zh) | 一种基于范数约束的波束旁瓣抑制方法 | |
Abdallah et al. | Multipath mitigation via synthetic aperture beamforming for indoor and deep urban navigation | |
CN104155662A (zh) | 基于gnss相关峰值探测器的自适应互干扰抑制方法 | |
CN103926603A (zh) | Gnss接收机极弱信号的跟踪方法 | |
CN106019333A (zh) | 一种基于非相干鉴别器的北斗导航信号矢量跟踪方法 | |
Lohan et al. | Accuracy limits in multi-GNSS | |
CN107607971A (zh) | 基于gnss共视时间比对算法的时间频率传递方法及接收机 | |
CN102565825B (zh) | 接收信号可靠度判定装置、方法及码相位误差算出方法 | |
CN104199067A (zh) | 多径环境下全球卫星导航***接收机无模糊处理方法 | |
CN105954772B (zh) | 一种稳健无偏的导航信号矢量跟踪方法 | |
Miller et al. | Multipath effects in GPS receivers: a primer | |
Musumeci et al. | Design of a very high sensitivity acquisition system for a space GNSS receiver | |
Hsu | Integration of vector tracking loop and multipath mitigation technique and its assessment | |
Sadrieh et al. | Doppler characterization of a mobile GNSS receiver in multipath fading channels | |
CN102508265B (zh) | 基于信号分离估计理论的卫星导航信号多径干扰抑制方法 | |
Chang et al. | An adaptive multipath mitigation filter for GNSS applications | |
Fohlmeister et al. | Dual Kalman filtering based GNSS phase tracking for scintillation mitigation | |
Tamazin | High resolution signal processing techniques for enhancing GPS receiver performance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 |