CN106842263A - 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 - Google Patents
基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 Download PDFInfo
- Publication number
- CN106842263A CN106842263A CN201710046067.2A CN201710046067A CN106842263A CN 106842263 A CN106842263 A CN 106842263A CN 201710046067 A CN201710046067 A CN 201710046067A CN 106842263 A CN106842263 A CN 106842263A
- Authority
- CN
- China
- Prior art keywords
- signal
- satellite
- navigation
- array
- carrier
- 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.)
- Pending
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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
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
本发明公开了一种基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,包括以下步骤:坐标系的建立;接收机对卫星导航信号进行接收;将接收到的卫星信号进行空时处理;空时相干约束,利用相干的本地合成载波和扩频信号,实现阵列的相干调向约束,相干调向约束后的数据送入滤波器组,最后提取观测量送至后端实现导航解算;卫星导航信号的捕获跟踪和矢量跟踪,利用空时联合相干约束,得到相对于载体坐标系原点的卫星信号,对该信号进行载波和码相位鉴相,则可获取该卫星信号的星历数据,利用该通道的EPL数据进一步获取载波和码相位,形成观测量,并可实现定位运算。本发明利用空时相干约束,解决由于阵列信号的空时处理造成的精度下降问题。
Description
技术领域:
本发明涉及一种卫星导航方法,特别是涉及一种基于导航信号空时处理与矢量跟踪相结合的卫星导航方法。
背景技术:
近年来卫星导航在国民经济和军事国防中得到广泛的应用,中美欧俄日等国家和地区都在积极发展自主的卫星导航系统及终端应用系统,以期提供更优的服务应用能力,特别是随着卫星导航的深入应用,一些新的导航的需求也浮现出来,迫切需要提升卫星导航的定位、定姿、抗干扰、弱信号跟踪和多系统融合等方面的能力。
卫星导航能力提升主要体现在系统能力和终端应用能力的建设发展上,其中终端朝着低成本和高性能两个方向发展,特别是在导航接收机方面的抗干扰、多系统融合、组合导航等技术,是目前的研究热点,也是提高导航抗干扰能力和提升导航定位精度、可靠性的关键技术,国内外已相继提出和开展了导航信号的矢量跟踪与导航解算、空时抗干扰等新技术的研究,并提出了信号的空时处理(Space-Time Processing)、矢量跟踪(VectorTracking)等新的导航信号处理与导航解算方法,通过充分发掘导航信号的时间空间结构和相干关系,融合信号处理与导航解算过程,提升导航性能,取得了较好的效果。
发明内容:
本发明所要解决的技术问题是:克服现有技术的不足,提供一种能够提高精度的基于导航信号空时处理与矢量跟踪相结合的卫星导航方法。
本发明为解决技术问题所采取的技术方案是:
一种基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,包括以下步骤:
步骤一:坐标系的建立:设载体安装阵列P,共有N个理想全向点状天线阵元,建立载体坐标系b,则阵元n在载体坐标系中的坐标则阵列p在载体坐标系中的坐标记为:
若某卫星信号si,在载体坐标系的入射方向为其中φ、θ分别为入射方位角和极距角,则其方向余弦为:
则地心地固坐标系e中,阵列p的坐标为:
在坐标系e中,载体坐标b的原点坐标为坐标系e旋转到载体坐标b的欧拉角为
则:
则坐标系e转换到坐标系b的转换矩阵为:
则坐标系b转换到坐标系e的转换矩阵为:
信号si若在e坐标系的方向余弦为则有:
步骤二:接收机对卫星导航信号进行接收;
对于卫星i发射的某导航信号si,设共有M个卫星,则对于卫星i,其发射的信号si的一般形式为:
其中,Ai为信号幅度,di为导航数据,Ci为该信号的PRN码,ωi表示该信号的中心频率或载波频率;
对于信号si,有一理想阵列p接收,则阵元pn接收到的卫星信号为延时后的si:
为阵元pn接收到信号si幅度,为信号si的传输时延;
对于阵元pn接收到的全部卫星信号及噪声可写为:
步骤三:将接收到的卫星信号进行载波信号空时处理和PRN码信号空时处理;
1)载波信号空时处理
若为信号si的载波信号,为波长,入射方向为ai(φi,θi),则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于卫星导航信号,若剥离掉调制信号,则余留单频信号其带宽bs≈0,属于窄带信号,对于窄带信号,为了获取阵列处理增益,应将每个阵元输出补偿一个时移τi,使各个阵元输出的信号恢复到共同相位上,再进行信号相加,利用信号的空间相干特性实现阵列接收的调向,在期望信号方向上形成波峰,信号的延时补偿τn可用中心频率相移来实现。
定义:
其中,ks为期望的平面波波数,且若忽略各类噪声,阵列调向后接收到卫星信号的响应为:
由公式可得总的频率-波数响应为
2)PRN码信号空时处理
对于卫星导航信号si,若剥离掉载波信号则得到某种调制信号该信号由导航数据和PRN码共同构成,剥离掉导航数据,该信号具有PRN码的周期特性和相关特性,可以在更长时间和空间尺度上,实现相干处理。
定义分别为已调制信号的等效频率、波长,和si具有相同信号到达方向,则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于PRN码信号是通过对每个阵元输出补偿一个时移,使各个阵元输出信号恢复到共同相位上,再进行信号相加,从而实现阵列接收的调向,在信号到达方向上形成一个波峰,对于卫星导航信号,PRN码周期一般为1ms,约为300km,所以不能将其作为窄带信号处理,不能将时域内的时移转化为频域的相移来实现。
若为期望的卫星信号方向或k空间的主响应轴,根据公式可得总的频率-波数响应为:
步骤四:空时相干约束
通过剥离掉有用信号的载波和扩频码,使有用信号带宽从扩频码带宽降低到数据带宽,从而在后继的低通滤波和波数成形运算中,降低对滤波器的相位处理要求,抑制检测噪声,提高观测量质量和导航解算精度。
在导航解算建立后导航信号到达方向已知的情况下,利用相干的本地合成载波和扩频信号,按照某个待接收卫星导航信号到达方向对阵列进行空时相干运算,实现阵列的相干调向约束,相干调向约束后的数据送入滤波器组,进行积分清洗或波束成形运算,最后提取观测量送至后端实现导航解算。
相关调向实现相干信号的解调解扩和阵列调向,通过载波的相关运算,实现了下变频和单频信号调向,而不改变扩频信号的相位突变时刻,利用扩频信号的相关运算,提取出扩频序列时刻,从而保证了观测量不包含由于波束成形和滤波造成的噪声,滤波器组实现积分清洗滤波或波束成形,由于在相关调向环节已经形成了基于载波的调相,因此此处的波束成形可在复符号基础上组阵,以获得更高的多径干扰抑制能力。
为了便于准确分析相关调向约束的性能,特作如下假设:接收机已经完全与卫星信号同步,忽略各类噪声和频率漂移的影响且以积分清洗即均匀加权的方法统一分析;以单个卫星信号为例分析,若si的入射方位角和极距角为ai(φi,θi):
1)依据公式计算阵元n延时
2)解调和解扩
3)空域滤波(均匀加权)
4)时域滤波(均匀加权)
按照数据d的周期Td内,积分得到:
步骤五:卫星导航信号的捕获跟踪和矢量跟踪
接收机初始工作时,由于接收信号视距动态造成多普勒频移,而载波环路检测带宽较窄,需要遍历所有可能的多普勒频移带宽即频域搜索,由于接收机钟差未知,而码鉴相范围只有±0.5个码片,因此需要遍历所有码状态即时域搜索,实现对卫星信号的捕获,之后通过锁相环路实现信号的跟踪,这个过程包括载波和码相位的检测、环路滤波器、维特比算法等。
对于空时联合矢量跟踪接收机,不仅要获取接收机的位置信息,还要获取其姿态信息,分别利用正交检测和迟早门检测,形成每个阵元接收信号的IQ支路采样值和阵列中心等效接收信号IQ支路采样值
在积分和清洗时刻,阵元n形成IQ支路采样值下面以阵元n为例,分别介绍单个阵元载波和码相位算法:
对于载波相位检测采用即时门IQ支路采样值计算,则在t时刻,载波相位差为:
假设在时刻t1、t2获得的IQ支路采样值分别为则该观测时刻频率差为:
其中ATAN2为四象限反正切,在高信噪比和低信噪比时最佳(最大似然估计),鉴相斜率与信号幅度无关。
对于码相位检测采用迟门和早门IQ支路采样值计算,通常利用归一化的超前滞后包络鉴别:
其中,当输入误差在±0.5个码片以内时,能够输出接近真实的误差。
对于最大似然估计采用维特比算法实现其最大似然的符号位同步,可采用如下算法:
(1)在第i个PRN码周期内,I支路输出信号模型为:
Ii=Adi+nIi
其中,di为数据位,nIi为第i个周期内的噪声,对进行20个PRN周期积累,若没有跨数据位,即di在累积期间保持不变,则获得累积最大值。
(2)由于存在20个不确定的数据位边,故构成20个不同的观测数据组,用下标φ=1,2,…,20表示,定义观测数据为sm,φ:
其中m表示以20个PRN周期为单位的观测数,则有:
sm,φ=Am,φdm,φ+nm,φ
(3)对于给定的序列:
采用VA算法求解,并设定每一个φ由一个单独的状态图来表示,每个状态图具有2个状态G和Κ,分别对应数据位值﹢1和﹣1。
(4)每1ms更新状态图,并将状态图φ加1,根据状态图状态G或Κ,分别对状态累加值加或减去Ι支路数据,当φ大于20时,φ设为1并输出状态累加值。
(5)对每个状态图进行比较,选取最大值作为输出,该最大值所对应的φ,代表了数据位。
利用空时联合相干约束,得到相对于载体坐标系原点的卫星信号si(t),对该信号进行载波和码相位鉴相,得到IQ支路采样值对该数据进行VA符号同步判别和数据剥离,则可获取该卫星信号的星历数据,利用该通道的EPL数据进一步获取载波和码相位,形成观测量,并可实现定位运算。
若实现位置测量,仍需对实际的每个阵元进行观测,由于接收的同一个信号,故可利用上述si(t)的观测值和同步值。
对于阵元n的IQ支路采样值利用si(t)同步指示,实现数据剥离,形成通道的EPL数据,经载波和码相位鉴相,形成阵元n的观测量。
在步骤五中,载波的相位检测由锁相环实现,即将数字中频信号与接收机复现的正交载波信号分别相乘累积,得到IQ支路信号;锁相环鉴相即通过IQ支路信号,利用相干原理,鉴别出本地复现信号与接收信号的相位差;码相位的检测一般采用延时锁定环实现,即通过比复现PRN码相位早半个码片和迟半个码片的信号,与剥离载波后的IQ支路分别相关,利用迟早门数值鉴别复现PRN 码与接收到的卫星信号码相位差,环路滤波通过对鉴相函数输出进行滤波处理。
在步骤五中,矢量跟踪是将信号跟踪和位置速度估计结合在一起,对于矢量跟踪典型的状态向量通常选取位置、速度、钟差和钟漂,如下式所示:
对于通道j的残差,其EKF的测量方程为:
在矢量跟踪中,对于N个卫星信号的伪距是利用N个伪距残差估计四个状态(位置和钟差),对于N个信号的伪距残差估计值,可认为是真实伪距残差加上白噪声:
在矢量跟踪中,对于N个伪距残差,估计位置误差和钟差,有:
则上式的最小二乘估计及其协方差为:
则对于矢量跟踪,其伪距协方差为:
获取各阵元的卫星信号观测值后,期望建立观测值与阵元姿态直接的关系,以满足导航递推解算要求,设阵列P具有N个阵元,在载体坐标系b中,阵元位置分别为由于运动,该载体坐标系移动到新位置,并形成新的载体坐标系c。从坐标系b到c,以b坐标系为参考,则坐标系的运动可表示为平移分量εp和旋转分量εψ,其中,旋转分量按照欧拉转动角旋转:
若载体和阵元是刚体,则在新的坐标系c中,阵元的坐标数值不变。则阵元pn新位置p′n在载体坐标系中坐标为:
若旋转角度εψ较小,则有:
对于导航信号si,ui为其方向余弦,则由于空间位置变化δpn引起的阵元pn通道处理延时变化为:
其中,可通过对载波或码相位检测值获得:
为载波波长的整数倍,当载体坐标移动较小时,有则可变化为:
对于载体坐标原点的等效接收信号,有
对于所有的M个卫星信号:
解上方程,可求得εp。
则最小二乘估计为:
有上式可得,在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εp。
定义Rp为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
对应其他N个阵元,写成εψ的形式:
对于卫星信号i的N个阵元的观测量,写成矩阵形式:
若:
组合M个卫星信号:
则最小二乘估计为:
有上式可得,在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εψ。
定义Rψ为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
本发明的积极有益效果是:
1、本发明将导航信号的空时特性引入矢量跟踪环路中,充分利用卫星导航信号的空域和时频域相干特性,融合卫星信号的空时处理、矢量跟踪和导航解算过程,进一步提高导航接收机的综合性能。
2、本发明具有广泛的应用前景,可用于卫星导航的多系统融合、微弱信号跟踪、抗干扰与多径抑制、精密定位和组合导航等应用场合;对研究解决导航信号接收处理技术发展所必然面对的空时处理问题,探讨下一代卫星导航信号的空时联合矢量跟踪架构具有积极的作用和意义。
3、本发明利用空时相干约束,解决由于阵列信号的空时处理造成的精度下降问题;建立了天线阵列位置姿态偏差与导航信号相位检测的解析关系,获得了具有明确空间几何意义的观测量,对实现基于空间几何测量的异种导航定位传感器组合具有重要意义。
附图说明:
图1是本发明的流程示意图;
图2是本发明中导航信号阵列的相移调向示意图;
图3是本发明中导航信号阵列的时移调向示意图;
图4是本发明中导航信号空时相干约束示意图;
图5是本发明中阵列信号鉴相示意图;
图6是本发明中VA算法结构图和观测量形成示意图;
图7是本发明算法模型示意图。
具体实施方式:
下面结合附图和具体实施例对本发明作进一步的解释和说明:
参见图1、图2、图3、图4、图5、图6和图7。
实施例:一种基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,包括以下步骤:
步骤一:坐标系的建立:设载体安装阵列P,共有N个理想全向点状天线阵元,建立载体坐标系b,则阵元n在载体坐标系中的坐标则阵列p在载体坐标系中的坐标记为:
若某卫星信号si,在载体坐标系的入射方向为其中φ、θ分别为入射方位角和极距角,则其方向余弦为:
则地心地固坐标系e中,阵列p的坐标为:
在坐标系e中,载体坐标b的原点坐标为坐标系e旋转到载体坐标b的欧拉角为
则:
则坐标系e转换到坐标系b的转换矩阵为:
则坐标系b转换到坐标系e的转换矩阵为:
信号si若在e坐标系的方向余弦为则有:
步骤二:接收机对卫星导航信号进行接收;
对于卫星i发射的某导航信号si,设共有M个卫星,则对于卫星i,其发射的信号si的一般形式为:
其中,Ai为信号幅度,di为导航数据,Ci为该信号的PRN码,ωi表示该信号的中心频率或载波频率;
对于信号si,有一理想阵列p接收,则阵元pn接收到的卫星信号为延时后的si:
为阵元pn接收到信号si幅度,为信号si的传输时延;
对于阵元pn接收到的全部卫星信号及噪声可写为:
步骤三:将接收到的卫星信号进行载波信号空时处理和PRN码信号空时处理;
1)载波信号空时处理
若为信号si的载波信号,为波长,入射方向为ai(φi,θi),则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于卫星导航信号,若剥离掉调制信号,则余留单频信号其带宽bs≈0,属于窄带信号,对于窄带信号,为了获取阵列处理增益,应将每个阵元输出补偿一个时移τi,使各个阵元输出的信号恢复到共同相位上,再进行信号相加,利用信号的空间相干特性实现阵列接收的调向,在期望信号方向上形成波峰,信号的延时补偿τn可用中心频率相移来实现。
定义:
其中,ks为期望的平面波波数,且若忽略各类噪声,阵列调向后接收到卫星信号的响应为:
由公式可得总的频率-波数响应为
2)PRN码信号空时处理
对于卫星导航信号si,若剥离掉载波信号则得到某种调制信号该信号由导航数据和PRN码共同构成,剥离掉导航数据,该信号具有PRN码的周期特性和相关特性,可以在更长时间和空间尺度上,实现相干处理。
定义分别为已调制信号的等效频率、波长,和si具有相同信号到达方向,则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于PRN码信号是通过对每个阵元输出补偿一个时移,使各个阵元输出信号恢复到共同相位上,再进行信号相加,从而实现阵列接收的调向,在信号到达方向上形成一个波峰,对于卫星导航信号,PRN码周期一般为1ms,约为300km,所以不能将其作为窄带信号处理,不能将时域内的时移转化为频域的相移来实现。
若为期望的卫星信号方向或k空间的主响应轴,根据公式可得总的频率-波数响应为:
步骤四:空时相干约束
通过剥离掉有用信号的载波和扩频码,使有用信号带宽从扩频码带宽降低到数据带宽,从而在后继的低通滤波和波数成形运算中,降低对滤波器的相位处理要求,抑制检测噪声,提高观测量质量和导航解算精度。
在导航解算建立后导航信号到达方向已知的情况下,利用相干的本地合成载波和扩频信号,按照某个待接收卫星导航信号到达方向对阵列进行空时相干运算,实现阵列的相干调向约束,相干调向约束后的数据送入滤波器组,进行积分清洗或波束成形运算,最后提取观测量送至后端实现导航解算;
相关调向实现相干信号的解调解扩和阵列调向,通过载波的相关运算,实现了下变频和单频信号调向,而不改变扩频信号的相位突变时刻,利用扩频信号的相关运算,提取出扩频序列时刻,从而保证了观测量不包含由于波束成形和滤波造成的噪声,滤波器组实现积分清洗滤波或波束成形,由于在相关调向环节已经形成了基于载波的调相,因此此处的波束成形可在复符号基础上组阵,以获得更高的多径干扰抑制能力。
为了便于准确分析相关调向约束的性能,特作如下假设:接收机已经完全与卫星信号同步,忽略各类噪声和频率漂移的影响且以积分清洗即均匀加权的方法统一分析;以单个卫星信号为例分析,若si的入射方位角和极距角为ai(φi,θi):
1)依据公式计算阵元n延时
2)解调和解扩
3)空域滤波(均匀加权)
4)时域滤波(均匀加权)
按照数据d的周期Td内,积分得到:
步骤五:卫星导航信号的捕获跟踪和矢量跟踪
接收机初始工作时,由于接收信号视距动态造成多普勒频移,而载波环路检测带宽较窄,需要遍历所有可能的多普勒频移带宽即频域搜索,由于接收机钟差未知,而码鉴相范围只有±0.5个码片,因此需要遍历所有码状态即时域搜索,实现对卫星信号的捕获,之后通过锁相环路实现信号的跟踪,这个过程包括载波和码相位的检测、环路滤波器、维特比算法等。
对于空时联合矢量跟踪接收机,不仅要获取接收机的位置信息,还要获取其姿态信息,分别利用正交检测和迟早门检测,形成每个阵元接收信号的IQ支路采样值和阵列中心等效接收信号IQ支路采样值
在积分和清洗时刻,阵元n形成IQ支路采样值下面以阵元n为例,分别介绍单个阵元载波和码相位算法:
对于载波相位检测采用即时门IQ支路采样值计算,则在t时刻,载波相位差为:
假设在时刻t1、t2获得的IQ支路采样值分别为则该观测时刻频率差为:
其中ATAN2为四象限反正切,在高信噪比和低信噪比时最佳(最大似然估计),鉴相斜率与信号幅度无关;
对于码相位检测采用迟门和早门IQ支路采样值计算,通常利用归一化的超前滞后包络鉴别:
其中,当输入误差在±0.5个码片以内时,能够输出接近真实的误差。
对于最大似然估计采用维特比算法实现其最大似然的符号位同步,可采用如下算法:
(1)在第i个PRN码周期内,I支路输出信号模型为:
Ii=Adi+nIi
其中,di为数据位,nIi为第i个周期内的噪声,对进行20个PRN周期积累,若没有跨数据位,即di在累积期间保持不变,则获得累积最大值。
(2)由于存在20个不确定的数据位边,故构成20个不同的观测数据组,用下标φ=1,2,…,20表示,定义观测数据为sm,φ:
其中m表示以20个PRN周期为单位的观测数,则有:
sm,φ=Am,φdm,φ+nm,φ
(3)对于给定的序列:
采用VA算法求解,并设定每一个φ由一个单独的状态图来表示,每个状态图具有2个状态G和Κ,分别对应数据位值﹢1和﹣1。
(4)每1ms更新状态图,并将状态图φ加1,根据状态图状态G或Κ,分别对状态累加值加或减去Ι支路数据,当φ大于20时,φ设为1并输出状态累加值。
(5)对每个状态图进行比较,选取最大值作为输出,该最大值所对应的φ,代表了数据位。
利用空时联合相干约束,得到相对于载体坐标系原点的卫星信号si(t),对该信号进行载波和码相位鉴相,得到IQ支路采样值对该数据进行VA符号同步判别和数据剥离,则可获取该卫星信号的星历数据,利用该通道的EPL数据进一步获取载波和码相位,形成观测量,并可实现定位运算。
若实现位置测量,仍需对实际的每个阵元进行观测,由于接收的同一个信号,故可利用上述si(t)的观测值和同步值。
对于阵元n的IQ支路采样值利用si(t)同步指示,实现数据剥离,形成通道的EPL数据,经载波和码相位鉴相,形成阵元n的观测量。
在步骤五中,载波的相位检测由锁相环实现,即将数字中频信号与接收机复现的正交载波信号分别相乘累积,得到IQ支路信号;锁相环鉴相即通过IQ支路信号,利用相干原理,鉴别出本地复现信号与接收信号的相位差;码相位的检测一般采用延时锁定环实现,即通过比复现PRN码相位早半个码片和迟半个码片的信号,与剥离载波后的IQ支路分别相关,利用迟早门数值鉴别复现PRN码与接收到的卫星信号码相位差,环路滤波通过对鉴相函数输出进行滤波处理。
在步骤五中,矢量跟踪是将信号跟踪和位置速度估计结合在一起,对于矢量跟踪典型的状态向量通常选取位置、速度、钟差和钟漂,如下式所示:
对于通道j的残差,其EKF的测量方程为:
在矢量跟踪中,对于N个卫星信号的伪距是利用N个伪距残差估计四个状态(位置和钟差),对于N个信号的伪距残差估计值,可认为是真实伪距残差加上白噪声:
在矢量跟踪中,对于N个伪距残差,估计位置误差和钟差,有:
则上式的最小二乘估计及其协方差为:
则对于矢量跟踪,其伪距协方差为:
获取各阵元的卫星信号观测值后,期望建立观测值与阵元姿态直接的关系,以满足导航递推解算要求,设阵列P具有N个阵元,在载体坐标系b中,阵元位置分别为由于运动,该载体坐标系移动到新位置,并形成新的载体坐标系c。从坐标系b到c,以b坐标系为参考,则坐标系的运动可表示为平移分量εp和旋转分量εψ,其中,旋转分量按照欧拉转动角旋转:
若载体和阵元是刚体,则在新的坐标系c中,阵元的坐标数值不变。则阵元pn新位置p′n在载体坐标系中坐标为:
若旋转角度εψ较小,则有:
对于导航信号si,ui为其方向余弦,则由于空间位置变化δpn引起的阵元pn通道处理延时变化为:
其中,可通过对载波或码相位检测值获得:
为载波波长的整数倍,当载体坐标移动较小时,有则可变化为:
对于载体坐标原点的等效接收信号,有
对于所有的M个卫星信号:
解上方程,可求得εp。
则最小二乘估计为:
有上式可得,在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εp。
定义Rp为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
对应其他N个阵元,写成εψ的形式:
对于卫星信号i的N个阵元的观测量,写成矩阵形式:
若:
组合M个卫星信号:
则最小二乘估计为:
有上式可得,在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εψ。
定义Rψ为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
以上所述,仅是本发明的优先实施例而已,并未对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的保护范围内。
Claims (4)
1.一种基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,其特征在于:包括以下步骤:
步骤一:坐标系的建立:设载体安装阵列P,共有N个理想全向点状天线阵元,建立载体坐标系b,则阵元n在载体坐标系中的坐标则阵列p在载体坐标系中的坐标记为:
若某卫星信号si,在载体坐标系的入射方向为其中φ、θ分别为入射方位角和极距角,则其方向余弦为:
则地心地固坐标系e中,阵列p的坐标为:
在坐标系e中,载体坐标b的原点坐标为坐标系e旋转到载体坐标b的欧拉角为
则:
则坐标系e转换到坐标系b的转换矩阵为:
则坐标系b转换到坐标系e的转换矩阵为:
信号si若在e坐标系的方向余弦为则有:
步骤二:接收机对卫星导航信号进行接收;
对于卫星i发射的某导航信号si,设共有M个卫星,则对于卫星i,其发射的信号si的一般形式为:
其中,Ai为信号幅度,di为导航数据,Ci为该信号的PRN码,ωi表示该信号的中心频率或载波频率;
对于信号si,有一理想阵列p接收,则阵元pn接收到的卫星信号为延时后的si:
为阵元pn接收到信号si幅度,为信号si的传输时延;
对于阵元pn接收到的全部卫星信号及噪声可写为:
步骤三:将接收到的卫星信号进行空时处理;
1)载波信号空时处理
若为信号si的载波信号,为波长,入射方向为ai(φi,θi),则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于卫星导航信号,若剥离掉调制信号,则余留单频信号其带宽bs≈0,属于窄带信号,对于窄带信号,为了获取阵列处理增益,应将每个阵元输出补偿一个时移τi,使各个阵元输出的信号恢复到共同相位上,再进行信号相加,利用信号的空间相干特性实现阵列接收的调向,在期望信号方向上形成波峰,信号的延时补偿τn可用中心频率相移来实现;
定义:
其中,ks为期望的平面波波数,且若忽略各类噪声,阵列调向后接收到卫星信号的响应为:
由公式可得总的频率-波数响应为
2)PRN码信号空时处理
对于卫星导航信号si,若剥离掉载波信号则得到某种调制信号该信号由导航数据和PRN码共同构成,剥离掉导航数据,该信号具有PRN码的周期特性和相关特性,可以在更长时间和空间尺度上,实现相干处理;
定义分别为已调制信号的等效频率、波长,和si具有相同信号到达方向,则信号的波数为:
其中ui为信号si的方向余弦,则该信号的阵列流行矢量:
对于PRN码信号是通过对每个阵元输出补偿一个时移,使各个阵元输出信号恢复到共同相位上,再进行信号相加,从而实现阵列接收的调向,在信号到达方向上形成一个波峰,若为期望的卫星信号方向或k空间的主响应轴,根据公式可得总的频率-波数响应为:
步骤四:空时相干约束
通过剥离掉有用信号的载波和扩频码,使有用信号带宽从扩频码带宽降低到数据带宽,从而在后继的低通滤波和波数成形运算中,降低对滤波器的相位处理要求,抑制检测噪声,提高观测量质量和导航解算精度;
在导航解算建立后导航信号到达方向已知的情况下,利用相干的本地合成载波和扩频信号,按照某个待接收卫星导航信号到达方向对阵列进行空时相干运算,实现阵列的相干调向约束,相干调向约束后的数据送入滤波器组,进行积分清洗或波束成形运算,最后提取观测量送至后端实现导航解算;
步骤五:卫星导航信号的捕获跟踪和矢量跟踪
接收机初始工作时,由于接收信号视距动态造成多普勒频移,而载波环路检测带宽较窄,需要遍历所有可能的多普勒频移带宽即频域搜索,由于接收机钟差未知,而码鉴相范围只有±0.5个码片,因此需要遍历所有码状态即时域搜索,实现对卫星信号的捕获,之后通过锁相环路实现信号的跟踪,这个过程包括载波和码相位的检测、环路滤波器、维特比算法等;
利用空时联合相干约束,得到相对于载体坐标系原点的卫星信号si(t),对该信号进行载波和码相位鉴相,得到IQ支路采样值对该数据进行VA符号同步判别和数据剥离,则可获取该卫星信号的星历数据,利用该通道的EPL数据进一步获取载波和码相位,形成观测量,并可实现定位运算;获取各阵元的卫星信号观测值后,期望建立观测值与阵元姿态直接的关系,以满足导航递推解算要求。
2.根据权利要求1所述的基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,其特征是:在步骤五中,载波的相位检测由锁相环实现,即将数字中频信号与接收机复现的正交载波信号分别相乘累积,得到IQ支路信号;锁相环鉴相即通过IQ支路信号,利用相干原理,鉴别出本地复现信号与接收信号的相位差;码相位的检测一般采用延时锁定环实现,即通过比复现PRN码相位早半个码片和迟半个码片的信号,与剥离载波后的IQ支路分别相关,利用迟早门数值鉴别复现PRN码与接收到的卫星信号码相位差,环路滤波通过对鉴相函数输出进行滤波处理。
3.根据权利要求1所述的基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,其特征是:在步骤五中,矢量跟踪是将信号跟踪和位置速度估计结合在一起,对于矢量跟踪典型的状态向量通常选取位置、速度、钟差和钟漂,如下式所示:
对于通道j的残差,其EKF的测量方程为:
在矢量跟踪中,对于N个卫星信号的伪距是利用N个伪距残差估计四个状态(位置和钟差),对于N个信号的伪距残差估计值,可认为是真实伪距残差加上白噪声:
在矢量跟踪中,对于N个伪距残差,估计位置误差和钟差,有:
则上式的最小二乘估计及其协方差为:
则对于矢量跟踪,其伪距协方差为:
4.根据权利要求1所述的基于导航信号空时处理与矢量跟踪相结合的卫星导航方法,其特征是:在步骤五中,获取各阵元的卫星信号观测值后,期望建立观测值与阵元姿态直接的关系,以满足导航递推解算要求,其中,旋转分量按照欧拉转动角旋转:
在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εp。
定义Rp为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
在已知阵元坐标和卫星信号的方向余弦情况下,通过测量载波或码相位可获得平移分量εψ。
定义Rψ为M个卫星信号所对应的权阵,则加权的最小二乘估计为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710046067.2A CN106842263A (zh) | 2017-01-21 | 2017-01-21 | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710046067.2A CN106842263A (zh) | 2017-01-21 | 2017-01-21 | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106842263A true CN106842263A (zh) | 2017-06-13 |
Family
ID=59119406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710046067.2A Pending CN106842263A (zh) | 2017-01-21 | 2017-01-21 | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842263A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782219A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于相关性的信号累积测相方法 |
CN111025367A (zh) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | 基于导航卫星信号双短基线的定姿方法 |
CN111025358A (zh) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | 基于导航卫星信号单短基线的定向方法 |
CN111147415A (zh) * | 2019-12-23 | 2020-05-12 | 东方红卫星移动通信有限公司 | 一种低轨卫星mapsk通信系统的相位跟踪方法 |
CN112014862A (zh) * | 2019-05-30 | 2020-12-01 | 上海海积信息科技股份有限公司 | 一种载波相位观测数据生成方法及装置 |
CN112433238A (zh) * | 2020-11-26 | 2021-03-02 | 中国人民解放军战略支援部队信息工程大学 | 基于大数据分析的可信gnss定位方法 |
CN112731476A (zh) * | 2020-10-23 | 2021-04-30 | 中国人民解放军63891部队 | 基于短码解扩和定时恢复联合的gps卫星m码信号解调方法 |
CN112763822A (zh) * | 2020-12-23 | 2021-05-07 | 北京无线电计量测试研究所 | 一种天线抗干扰测量装置及测量系统 |
CN116482724A (zh) * | 2023-03-09 | 2023-07-25 | 北京理工大学 | 一种导航接收机高精度观测量计算方法 |
CN117233807A (zh) * | 2023-11-13 | 2023-12-15 | 中国船舶集团有限公司第七〇七研究所 | 一种适应高海况的北斗三号捕获跟踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101900819A (zh) * | 2010-07-15 | 2010-12-01 | 哈尔滨工业大学 | 导航接收机盲自适应空时阵列抗干扰方法 |
CN102879790A (zh) * | 2011-07-13 | 2013-01-16 | 北京泰豪联星技术有限公司 | 基于数字波束形成与空时调零级联的抗干扰系统和方法 |
CN105044737A (zh) * | 2015-06-30 | 2015-11-11 | 北京时代民芯科技有限公司 | 一种导航接收机中检测观测量有效性的方法 |
-
2017
- 2017-01-21 CN CN201710046067.2A patent/CN106842263A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101900819A (zh) * | 2010-07-15 | 2010-12-01 | 哈尔滨工业大学 | 导航接收机盲自适应空时阵列抗干扰方法 |
CN102879790A (zh) * | 2011-07-13 | 2013-01-16 | 北京泰豪联星技术有限公司 | 基于数字波束形成与空时调零级联的抗干扰系统和方法 |
CN105044737A (zh) * | 2015-06-30 | 2015-11-11 | 北京时代民芯科技有限公司 | 一种导航接收机中检测观测量有效性的方法 |
Non-Patent Citations (1)
Title |
---|
贾学东: "卫星导航信号空时联合矢量跟踪技术研究", 《万方数据库》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782219A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于相关性的信号累积测相方法 |
CN112014862A (zh) * | 2019-05-30 | 2020-12-01 | 上海海积信息科技股份有限公司 | 一种载波相位观测数据生成方法及装置 |
CN112014862B (zh) * | 2019-05-30 | 2024-03-29 | 上海海积信息科技股份有限公司 | 一种载波相位观测数据生成方法及装置 |
CN111025367B (zh) * | 2019-12-05 | 2021-11-16 | 中国人民解放军63620部队 | 基于导航卫星信号双短基线的定姿方法 |
CN111025358A (zh) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | 基于导航卫星信号单短基线的定向方法 |
CN111025358B (zh) * | 2019-12-05 | 2021-11-16 | 中国人民解放军63620部队 | 基于导航卫星信号单短基线的定向方法 |
CN111025367A (zh) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | 基于导航卫星信号双短基线的定姿方法 |
CN111147415A (zh) * | 2019-12-23 | 2020-05-12 | 东方红卫星移动通信有限公司 | 一种低轨卫星mapsk通信系统的相位跟踪方法 |
CN111147415B (zh) * | 2019-12-23 | 2022-05-31 | 东方红卫星移动通信有限公司 | 一种低轨卫星mapsk通信系统的相位跟踪方法 |
CN112731476B (zh) * | 2020-10-23 | 2023-09-08 | 中国人民解放军63891部队 | 基于短码解扩和定时恢复联合的gps卫星m码信号解调方法 |
CN112731476A (zh) * | 2020-10-23 | 2021-04-30 | 中国人民解放军63891部队 | 基于短码解扩和定时恢复联合的gps卫星m码信号解调方法 |
CN112433238A (zh) * | 2020-11-26 | 2021-03-02 | 中国人民解放军战略支援部队信息工程大学 | 基于大数据分析的可信gnss定位方法 |
CN112433238B (zh) * | 2020-11-26 | 2023-07-21 | 中国人民解放军战略支援部队信息工程大学 | 基于大数据分析的可信gnss定位方法 |
CN112763822A (zh) * | 2020-12-23 | 2021-05-07 | 北京无线电计量测试研究所 | 一种天线抗干扰测量装置及测量系统 |
CN112763822B (zh) * | 2020-12-23 | 2023-01-06 | 北京无线电计量测试研究所 | 一种天线抗干扰测量装置及测量系统 |
CN116482724A (zh) * | 2023-03-09 | 2023-07-25 | 北京理工大学 | 一种导航接收机高精度观测量计算方法 |
CN116482724B (zh) * | 2023-03-09 | 2024-05-03 | 北京理工大学 | 一种导航接收机高精度观测量计算方法 |
CN117233807A (zh) * | 2023-11-13 | 2023-12-15 | 中国船舶集团有限公司第七〇七研究所 | 一种适应高海况的北斗三号捕获跟踪方法 |
CN117233807B (zh) * | 2023-11-13 | 2024-02-02 | 中国船舶集团有限公司第七〇七研究所 | 一种适应高海况的北斗三号捕获跟踪方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842263A (zh) | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 | |
JP3548853B2 (ja) | マルチビット相関器を備えたスペクトル拡散受信機 | |
CN103033828B (zh) | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 | |
CN101846746B (zh) | 一种基于gnss-r技术的载波相位测高装置及其方法 | |
CN102033236B (zh) | 一种卫星导航位置速度联合估计方法 | |
US7839916B1 (en) | Systems and methods for communication in a global positioning system (GPS) device | |
CN103777218B (zh) | Gnss/ins超紧组合导航系统的性能评估系统及方法 | |
US20120050103A1 (en) | Synthetic aperture device for receiving signals of a system comprising a carrier and means for determining its trajectory | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
CN104280746A (zh) | 一种惯性辅助gps的深组合半实物仿真系统 | |
CN107255822A (zh) | 多径环境下gnss接收机信号参数估计方法 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN104049262A (zh) | 一种基于矢量跟踪环路的波束形成抗干扰方法 | |
CN101246210A (zh) | 一种gps定位方法及基于双核处理器的gps软件接收机 | |
JPH11183586A (ja) | 特にgps型の、完全高周波航行用受信装置 | |
CN104459731A (zh) | 一种接收极弱gnss信号的高轨卫星定轨方法 | |
CN102253396A (zh) | 一种高动态gps载波环跟踪方法 | |
CN106802153B (zh) | 基于单频导航原始观测量地面处理的高精度测定轨方法 | |
CN106597492A (zh) | 卫星导航接收机及其抗远近效应的方法和室内定位方法 | |
Yang et al. | Post-correlation semi-coherent integration for high-dynamic and weak GPS signal acquisition | |
CN102023299B (zh) | 一种抑制多径的相关器方法 | |
CN109471137A (zh) | 一种自适应矢量跟踪环路的实现方法 | |
CN104931980B (zh) | 载波相位测量半周模糊解除方法 | |
ZhG | Entanglement difference of GNSS carrier phase for vehicle attitude determination | |
Pinell et al. | Receiver architectures for positioning with low earth orbit satellite signals: a survey |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |