首页> 中国专利> 一种GPS天线相位中心变化对相对定位精度的修正方法

一种GPS天线相位中心变化对相对定位精度的修正方法

摘要

本发明公开了一种GPS天线相位中心变化对相对定位精度的修正方法,包括以下步骤:两个目标均载有GPS接收机及GPS天线,利用GPS卫星信号对两个目标进行相对定位;确定两个目标的GPS天线相位中心变化;根据所述天线相位中心变化确定与目标相对定位精度修正量的关系;利用所述修正量对目标相对定位精度进行修正。采用本发明实现了根据天线相位中心变化修正目标相对定位结果,提高了目标相对定位精度。

著录项

  • 公开/公告号CN102981174A

    专利类型发明专利

  • 公开/公告日2013-03-20

    原文格式PDF

  • 申请/专利权人 航天东方红卫星有限公司;

    申请/专利号CN201210334296.1

  • 申请日2012-09-12

  • 分类号G01S19/40;

  • 代理机构中国航天科技专利中心;

  • 代理人安丽

  • 地址 100094 北京市5616信箱

  • 入库时间 2024-02-19 17:37:56

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2014-10-08

    授权

    授权

  • 2013-04-17

    实质审查的生效 IPC(主分类):G01S19/40 申请日:20120912

    实质审查的生效

  • 2013-03-20

    公开

    公开

说明书

技术领域

本发明涉及一种GPS相对定位精度修正方法,特别涉及一种利用GPS天 线相位中心变化对相对定位精度的修正方法。

背景技术

GPS(Global Position System,全球定位系统)天线与GPS接收机共同 完成对全球导航卫星(GPS卫星)信号的接收、处理及解算工作。GPS接收 机通过天线接收GPS导航卫星信号及相关数据信息。由于GPS天线本身的特 性,天线几何中心与相位中心不可能完全一致,两者在水平方向和高程方向均 存在偏差,并且随着接收到的GPS卫星的方位和高度的变化而变化。天线相 位中心的误差,必将对测量结果带来影响,尤其是在高精度的GPS测量中影 响不可忽视。

天线相位中心是指天线的辐射场分量等效的由这点辐射出去,辐射波为一 球面。理想天线存在唯一的相位中心,其等相面为球面,因此接收来自不同方 向的卫星信号时,天线因自身产生额外的相位差而造成定位测量结果的偏差。 然而在整个波束空间存在唯一相位中心的天线几乎是不存在的,只可能在主瓣 某一范围内相位保持相对恒定,因此接收天线在接收不同方向的卫星信号时会 引入额外的相位差异,而引入定位测量结果的误差。显然,天线类型不同,其 相位特性不同,对接收机测量精度的影响也不同。天线相位中心与几何中心以 及等相面关系示意图见图2,ARP为天线几何中心(天线参考点)一般定义为 天线上表面的中心点,E为天线的平均相位中心。可以看出,正常情况下相位 中心存在两个偏差。天线相位中心E与天线参考点ARP的偏差称为天线相位 中心偏移PCO(Phase Center Offset),实际相位中心P点与平均相位中心E点 的偏差称为天线相位中心变化PCV(Phase Center Variations),其数值等效于 实际等相位面与拟合等相位球面的偏差。

发明内容

本发明的技术解决问题是:提供了一种利用GPS天线相位中心变化对相 对定位精度的修正方法,实现了根据天线相位中心变化修正目标相对定位结果, 提高了目标相对定位精度。

本发明的技术解决方案是:

一种GPS天线相位中心变化对相对定位精度的修正方法,包括以下步骤:

1)两个目标均载有GPS接收机及GPS天线,利用GPS卫星信号对两个 目标进行相对定位;

2)确定2个目标的GPS天线相位中心变化;

3)根据GPS天线相位中心变化确定目标相对定位精度修正量:

单差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目 标相对定位精度修正量的关系:

σδX=PDOP·σ=PDOP·2σPCV

其中,σδX为用于表示目标相对定位精度修正量的标准差;PDOP为空间 位置精度因子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差 的标准差;σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。

双差方式进行相对定位时,根据所述天线相位中心变化通过下式确定与目 标相对定位精度修正量的关系:

σδX=DOP·σ=DOP·2σPCV

其中,σδX为用于表示目标相对定位精度修正量的标准差;DOP为精度因 子;σ为用于表示GPS天线相位中心变化在均值为0时的测量误差的标准差; σpcv为用于表示GPS天线相位中心变化在均值为0时的标准差。

4)利用目标相对定位精度修正量对目标相对定位精度进行修正。

单差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对 定位修正量的关系通过下式在GPS天线到各卫星产生的天线相位变化呈正态 分布,且不同卫星到接收机之间的相位中心变化互不相关时获得:

δXP-t12=-(ABTAB)-1ABTD12

其中,δXP为只考虑天线相位中心变化影响的δX2

δX2=δx2δy2δz2;[δx2,δy2,δz2]为接收机坐标偏差;

δt12=δt2-δt1

δt1、δt2分别为两个接收机的钟差,c为光速;

A为δX2的系数矩阵,A=k21(t1)l21(t1)m21(t1)k22(t1)l22(t1)m22(t1)·········k2n(t1)l2n(t1)m2n(t1);B=11···1;k2p(t1)=xp(t1)-x2(t1)0ρ2p(t1)0,l1p(t1)=yp(t1)-y2(t1)0ρ2p(t1)0,m2p(t1)=zp(t1)-z2(t1)0ρ2p(t1)0ρ2p(t1)0=(x2(t1)0-xp(t1))2+(y2(t1)0-yp(t1))2+(z2(t1)0-zp(t1)0)2,

其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。 D12为只含有天线相位中心变化信息,D12=Δdr121Δdr122···Δdr12n

分别为卫星p信号相对于接收机i的 方位角和仰角,p=1,2,…n。

双差方式进行相对定位时,所述GPS天线相位中心变化与所述目标相对定 位修正量的关系通过下式在GPS天线到各GPS卫星产生的天线相位变化呈正 态分布,且不同卫星到GPS接收机之间的相位中心变化互不相关时获得:

δXP=-(ATA)-1AT·ΔD

其中,δXP为只考虑天线相位中心变化影响的δX2

δX2=δx2δy2δz2;[δx2,δy2,δz2]为接收机坐标偏差;

A为δX2的系数矩阵,A=k21(t1)l21(t1)m21(t1)k22(t1)l22(t1)m22(t1)·········k2n(t1)l2n(t1)m2n(t1);

k2p(t1)=xp(t1)-x2(t1)0ρ2p(t1)0,l1p(t1)=yp(t1)-y2(t1)0ρ2p(t1)0,m2p(t1)=zp(t1)-z2(t1)0ρ2p(t1)0,

ρ2p(t1)0=(x2(t1)0-xp(t1))2+(y2(t1)0-yp(t1))2+(z2(t1)0-zp(t1)0)2,

其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。

ΔD为只含有天线相位中心变化信息,ΔD=Δdr1212Δdr1213···Δdr121n

Δdr12pq=Δdr12q-Δdr12p,

分别为卫星p信号相对于接收机i的 方位角和仰角,p=1,2,…n。

本发明与现有技术相比的优点在于:

1)采用本发明获得了利用GPS天线相位中心变化对GPS相对定位精度 的影响,实现了对目标相对定位精度的修正。

2)本发明给出了星载GPS天线相位中心变化与GPS相对定位精度之间的 定量关系,给出了精确的数学模型。现有的技术虽然分析到GPS天线相位中心 变化的存在,对GPS测量的水平位置和高程位置精度都会产生影响。但从来没 有提出过明确的定量关系和数学模型。

附图说明

图1为本发明流程图。

图2为本发明天线相位中心与几何中心以及等相面关系示意图

图3为本发明卫星信号的方位角和仰角示意图

图4为本发明天线相位中心变化导致的径向距离变化

图5为GPS单差相对定位示意图

图6为GPS双差相对定位示意图

图7为本发明天线相位中心变化对单差相对定位影响总误差实例

图8为本发明天线相位中心变化对单差相对定位各个方向上误差实例

图9为本发明天线相位中心变化对双差相对定位影响总误差实例

图10为本发明天线相位中心变化对双差相对定位各个方向上误差实例

图11为不同PDOP值下相对定位精度统计曲线(单差方式)

图12为不同DOP值下相对定位精度统计曲线(双差方式)

具体实施方式

下面就结合附图对本发明做进一步介绍。

由于星载GPS天线几何中心与相位中心不一致,随着接收到的GPS卫星 的高度和方位的变化而变化,天线相位中心的误差直接对测量结果(目标相对 定位精度)带来影响,尤其是在高精度的GPS测量中需要分析并消除部分影 响。

天线相位中心偏移PCO是一个固定不变的矢量,而天线相位中心变化PCV 则随着卫星入射信号高度角和方位角的变化而变化,因此会给目标相对定位精 度引入额外的相位偏差。如图1所示为本发明流程图,下面就针对流程图对本 发明进行介绍。

天线平均相位中心E点,PCO、PCV、天线参考点APR的关系如图2所 示。对于每个天线,向量表示天线相位中心偏移PCO,其长度为一个固定的 值。向量表示卫星信号方向上的单位向量,可以用卫星信号相对于大地局部 坐标系的方位角θ、仰角表示,如图3所示。在图3中,S为卫星位置,O为 接收机位置,向量ON为大地正北方向,OE为大地正东方向,OU为大地天顶 方向。向量OS在大地平面EON上的投影为向量OA,那么方位角θ为向量ON 到OA的角度,范围为0~360度;仰角为向量OA到向量OS的角度,范围 为0~90度。

天线相位中心变化量PCV与卫星信号相对于接收机天线盘面坐标系的方 位角θr和仰角有关,其值用表示。则在方向上实际相位中心相对于 天线参考点的距离变化如图4所示,可以表示为:

公式(1)左边表示距离观测量的偏差,右边的PCO向量与PCV值可通过天线相位特性实际测量数据给出。

GPS相对定位的基本原理是:用2台接收机分别安置在基线的2个端点, 同步观测相同的4颗以上GPS卫星,确定基线2个端点在地球坐标系中的相 对位置。

在2个观测站或多个观测站同步观测相同卫星的情况下,卫星的轨道误差、 卫星钟差、接收机钟差以及电离层和对流层的折射误差等对观测值的影响具有 一定的相关性。利用这些观测量的不同组合(求差)进行相对定位,可有效削弱 相关误差的影响,从而提高相对定位的精度。对GPS载波相位观测值的求差, 可以分别在卫星之间、接收机之间和不同历元之间进行,根据求差次数的不同 又可分为单差,双差等方式。实际工程应用中,可以采用单差或者双差方式。

单差(Single Different,SD)是指不同观测站同步观测相同卫星所得观测量 之差。

假设观测站1(接收机1)和观测站2(接收机2)分别在t1时刻(历元)对卫星p 进行了同步观测,如图5所示,则可得载波相位观测量:则单 差方程为:

在图5中,若在t1时刻在观测站1、观测站2同时对卫星p进行了载波相位 测量,得观测方程:

公式(3)和公式(4)中,表示t1时刻接收机i到卫星j的相位观测量, 表示t1时刻接收机i到卫星j的几何距离,分别表示t1时刻 接收机i到卫星j的电离层延时距离误差和对流层延时距离误差,δti表示接收 机i的钟差,δtj表示卫星j的钟差,表示观测初始时刻t0接收机i到卫星 j的整周未知数,f为载波频率,c为光速。

将公式(3)和公式(4)代入公式(2)中,可处理得:

δXP-t12=-(ABTAB)-1ABT(L-λCN)---(5)

其中,A=k21(t1)l21(t1)m21(t1)k22(t1)l22(t1)m22(t1)·········k2n(t1)l2n(t1)m2n(t1),δX2=δx2δy2δz2,B=11···1,k2p(t1)=xp(t1)-x2(t1)0ρ2p(t1)0,l2p(t1)=yp(t1)-y2(t1)0ρ2p(t1)0,m2p(t1)=zp(t1)-z2(t1)0ρ2p(t1)0ρ2p(t1)0=(x2(t1)0-xp(t1))2+(y2(t1)0-yp(t1))2+(z2(t1)0-zp(t1)0)2,

其中[xp(t1),yp(t1),zp(t1)]为t1时刻卫星p的位置坐标,p=1,2,…n。

L=ΔL121(t1)ΔL122(t1)···ΔL12n(t1),N=N121(t0)N122(t0)···N12n(t0).

λ为波长

表示t1时刻接收机i到卫星p的几何距离,i=1,2。

N12p(t0)=N2p(t0)-N1p(t0),

表示观测初始时刻t0接收机i到卫星p的整周未知数,i=1,2。

假设接收机i的天线相位中心变化在卫星j的信号方向上引起的观测距离 变化为分别为卫星j信号相对于接收机 i的方位角和仰角;

对于观测站1和观测站2,对其相位观测值进行修正,并经过一系列的处 理,可得

L=Ln-D12                        (6)

L=ΔL121(t1)ΔL122(t1)···ΔL12n(t1),D12=Δdr121Δdr122···Δdr12n.

则有,δXP-t12=-(ABTAB)-1ABT(Ln-D12-λCN)---(7)

其中,只有D12矩阵包含天线相位中心变化因素,因此,天线相位中心变化 引起的相对定位单差观测方式的误差为:

δXP-t12=-(ABTAB)-1ABTD12---(8)

双差(Double Different,DD)是指不同观测站同步观测同一组卫星,所得单 差观测值之差。

如图6所示,观测站1、观测站2在t1时刻同时观测了2颗卫星p、q,那 么对p、q卫星分别有如公式(2)所示的单差方程。若忽略大气折射残差,双差 观测方程为

经过一系列处理,可得

δX2=λ(ATA)-1AT(ΔL-CΔN)                    (10)

其中,A=Δk212(t1)Δl212(t1)Δm212(t1)Δk213(t1)Δl213(t1)Δm213(t1)·········Δk21n(t1)Δl21n(t1)Δm12n(t1),δX2=δx2δy2δz2,

ΔN=ΔN1212(t0)ΔN1213(t0)···ΔN121n(t0),ΔL=ΔL1212(t1)ΔL1213(t1)···ΔL121n(t1)

接收机i的天线相位中心变化在卫星j的信号方向上引起的观测距离变化为 分别为卫星j信号相对于接收机i的方位 角和仰角;

对于观测站1和观测站2,对其相位观测值进行修正,并经过一系列的处 理,可得

λΔL=λΔLn-ΔD                        (11)

其中,ΔD=Δdr1212Δdr1213···Δdr121n.

Δdr12pq=Δdr12q-Δdr12p,

参考双差观测的解算公式(9),得到

δX2=λ(ATA)-1AT(ΔLn-1λΔD-CΔN)---(12)

其中,只有ΔD矩阵包含天线相位中心变化因素,因此,天线相位中心变化 引起的双差观测方式的相对定位精度为

δXP=-(ATA)-1AT·ΔD                    (13)

利用载波相位测量的星载GPS天线相位中心变化与GPS接收机相对定位 精度的数学关系公式(8)和(13),分别建立星载GPS天线相位中心变化对相对 定位精度影响的单差、双差数学模型,并设计仿真软件,定量的分析天线相位 中心变化对相对定位精度的影响。从而实现接收机天线相位中心变化对相对定 位精度的影响数学模型的验证。该软件要求输入天线相位中心变化数据(自定义 格式或者标准格式),卫星位置数据(自定义格式或者标准格式),以及模型验证 数据。用于数学模型的仿真验证。

GPS天线相位中心变化引起的单差相对定位精度公式(8)中,假设GPS接 收机天线到各个卫星产生的相位中心变化均呈相同的正态分布,其均值为0, 并且不同卫星到接收机之间的相位中心变化互不相关。则三维空间单差相对定 位精度的标准差:

σδX=σx2+σy2+σz2=h11+h22+h33·2σPCV=PDOP·2σPCV---(14)

其中,PDOP定义为空间位置精度因子,它的值为

PDOP=h11+h22+h33---(15)

如果用hii代表权系数矩阵H的对角元素,其中i=1,2,...,4,H矩阵为 H=(CTC)-1,H通常称为权系数阵,他是一个4*4的对称矩阵。

因此,在样本数足够多的情况下,对于单差相对定位方式,三维空间相对定 位精度的标准差σδX等于天线相位中心变化产生的测量误差的标准差σ乘以空 间位置精度因子(PDOP),理论上σ为倍σpcv,σpcv为天线相位中心变化标 准差。

在GPS星座中,已知在遮蔽角为5度时,PDOP值小于3的概率约为99%。 进一步分析相对定位精度的分布区间,从图11可以看出,在单差相对定位时, 定位误差的标准差小于3倍天线相位中心变化产生的测量误差的标准差的置信 度为99%。因此,可以认为:在不考虑其他误差影响的情况下,对于单差相对定位 方式,相对定位精度的标准差小于3倍天线相位中心变化引起的测量误差的标 准差。

GPS天线相位中心变化引起的双差相对定位精度公式(13)中,假设GPS天 线到各个卫星产生的相位中心变化均呈相同的正态分布,其均值为0,并且不 同卫星到接收机之间的相位中心变化互不相关。则三维空间双差相对定位精度 的标准差:

σδX=σx2+σy2+σz2=h11+h22+h33·2σPCV=DOP·2σPCV---(16)

其中,DOP定义为精度因子,它的值为

DOP=h11+h22+h33---(17)

如果用hii代表权系数矩阵H的对角元素,其中i=1,2,3,H矩阵为 H=(ATA)-1,H通常称为权系数阵,它是一个3*3的对称矩阵。

因此,在样本数足够多的情况下,对于双差相对定位方式,三维空间相对定 位精度的标准差σδX等于天线相位中心变化产生的测量误差的标准差σ乘以精 度因子(DOP),理论上σ为倍σpcv,σpcv为天线相位中心变化标准差。

通常双差定位的精度因子DOP小于单差定位的精度因子PDOP。由于在 遮蔽角为5度时,单差定位的精度因子PDOP小于3的概率为99%。我们可 以得到DOP通常小于3。

进一步分析相对定位精度的分布区间,从图12可以看出,在双差相对定 位时,相对定位精度的标准差小于3倍天线相位中心变化产生的测量误差的标 准差的置信度为99%。因此,可以认为:在不考虑其他误差影响的情况下,对于双 差相对定位方式,相对定位精度的标准差小于3倍天线相位中心变化引起的测 量误差的标准差。

对于该统计结果,现有技术虽然分析了GPS天线几何中心与其电学相位中 心不一致,不仅在水平、高程方向存在一个固定的偏差,而且随着接收到的GPS 卫星的高度和方位的变化而变化。但没有做出过详细的统计特性估计。

实施例

观测站、参考站接收机均静止,其中观测站的空间坐标为[-113402.1865, -5504362.8394,3209404.3787](单位:米),参考站的空间坐标为 [-112409.5336,-5504317.7352,3209516.6545](单位:米)。接收机天线类型 为:XX-XX,PCV的三个分量分别表示North、East、Up方向上的长度。选择 的观测历元区间为2011-7-1 0:0:00至2011-7-1 1:0:00,历元采样间隔为2分 钟。考查天线相位中心变化PCV对相对定位精度的影响,PCV的数据采取标 准格式输入,单差方式的仿真结果如图7、图8所示,双差方式的仿真结果如 图9、图10所示。

图7和图9是相对定位精度总的长度,表示定位位置与实际位置之间的距 离误差;图8和图10中,“定位精度(东方向)”表示定位精度在大地正东方向上 的分量,“定位精度(北方向)”表示定位精度在大地正北方向上的分量,“定位精 度(水平方向)”表示定位精度在大地水平面上的分量,是东方向和北方向分量的 合成,“定位精度(高程方向)”表示定位精度在大地垂直方向上的分量,即天顶方 向上的分量。

在图7、图8、图9和图10中,星号表示数学模型分析得到的误差,圆圈 表示实际定位程序分析得到的误差。从仿真结果中可以看出,在仅讨论天线相 位中心变化对GPS相对定位精度的影响时,单差方式和双差方式,其相对定 位精度很相近,差别不大。采用单差方式或者双差方式,由实际工程应用的具 体情况决定,这里不作讨论。可以看出,相对定位精度修正量很小,约为0.0147 毫米,基本可以忽略,因此相对定位可以大大提高定位精度。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号