CN115639582B - 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 - Google Patents
一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 Download PDFInfo
- Publication number
- CN115639582B CN115639582B CN202211266864.9A CN202211266864A CN115639582B CN 115639582 B CN115639582 B CN 115639582B CN 202211266864 A CN202211266864 A CN 202211266864A CN 115639582 B CN115639582 B CN 115639582B
- Authority
- CN
- China
- Prior art keywords
- satellite
- station
- orbit determination
- central reference
- reference station
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 15
- 230000001360 synchronised effect Effects 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 16
- 239000005433 ionosphere Substances 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000007774 longterm Effects 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000005436 troposphere Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法。该方法包括:构建由中心基准站、若干定轨站和GEO卫星组成的GEO卫星运动学定轨系统,中心基准站包括时频基准系统、星地时间同步天线和基准接收机,定轨站包括卫星导航接收机。通过多星共视法实时测量各个定轨站与中心基准站之间的钟差,利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差,根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨。本发明利用共视法授时获得地面站的精确时间,通过多个定轨站同时测距,进行运动学定轨获得为卫星精确位置,为地面天线对卫星连续跟踪提供目标指引。
Description
技术领域
本发明涉及导航卫星技术领域,尤其涉及一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法。
背景技术
GEO(Geosynchronous Earth Orbit,地球同步轨道)卫星在航天中具有重要地位,也是卫星导航系统的重要的信息转发器,承载了星基增强服务、精密单点定位服务和短报文服务等高实时性服务的广播信息播发任务,为保证上述服务的连续性,需精确地卫星位置信息,保证地面天线对GEO卫星的连续跟踪。但GEO卫星轨道资源极其宝贵,按照东西位置保持精度优于0.1°要求,平均每二十多天就要进行一次轨道机动,GEO卫星连续使用需求与频繁的轨道机动之间的矛盾,给卫星位置精确确定带来了一定压力。
由于轨位资源宝贵,GEO卫星需要通过频繁轨道机动来维持位置,轨道机动期间受发动机推力的影响,GEO卫星轨道不再符合原来的动力学轨道,卫星动力学定轨不再适用。目前,现有技术中的卫星定轨方法包括地基测量定轨和天基测量定轨两大类,地基定轨技术包括测距跟踪定轨和干涉测角跟踪定轨,距离测量的技术手段包括雷达测距、激光测距、伪距和载波测距,角度测量的技术手段包括光学测角和射电干涉测角。
鉴于伪距和载波测距定轨最为经济、便捷和高效,卫星导航多采用地面定轨站的伪距和载波测距观测量进行卫星轨道测定。正常情况下,采用多星多站动力学定轨,综合利用多个定轨站对多个卫星的长弧段伪距和载波观测数据以及卫星动力学信息进行动力学定轨(下文简称多星定轨),得到高精度卫星位置。在轨道机动时,由于受发动机推力的影响,GEO卫星运动轨迹不再符合原有动力学模型,因此采用多个定轨站伪距观测数据进行不依靠卫星动力学模型的运动学定轨,解算GEO卫星位置。
GEO卫星轨道机动期间采用基于多星定轨钟差预报修正的运动学定轨方法。该方法以伪距作为观测量,伪距中包含卫星位置至接收机之间的几何距离、卫星钟差、定轨站钟差以及各种误差,定轨之前需要从伪距观测量中扣除定轨站钟差和卫星钟差的影响。卫星钟差通过星地双向时间同步获得,定轨站钟差则通过多星动力学定轨(定轨解算中不包含机动卫星)解算的钟差结果预报得到。以修正后的伪距作为输入,采用前方交会法,计算卫星位置。但是该方法存在计算量大、实时性低、定轨精度低、可选定轨站少等问题,影响天线对GEO卫星的跟踪。
上述现有技术中的卫星定轨方法的缺点包括:基于多星定轨钟差预报修正的运动学定轨,有计算量大、实时性低、定轨精度低和可选定轨站少等问题:
一是计算量大。多星定轨利用多个均匀分布的定轨站对卫星的长弧段观测数据,同时估计卫星轨道钟差、接收机坐标钟差、整周模糊度、对流层参数、电离层延迟、各类时延,以期获得较高的精度,由于待估参数数量大,法方程阶数高,导致计算量巨大,即使采用了消参数法、分区法减小计算量,但是每次定轨都要耗费十几甚至几十分钟。
二是实时性低。由于多星定轨计算量大、耗时较长,通常多星定轨采用定时启动策略,而轨道机动期间,卫星运动学定轨需要实时解算卫星位置,因此只能采用钟差预报的方法,通过对多星定轨计算得到的钟差进行拟合预报,获取定轨站钟差预报值。
三是定轨精度低。由于GEO卫星轨道高度较高,相对与我国国土位置偏置且定轨站分布范围较小,导致GEO卫星定轨的观测几何结构差,定轨站钟差预报误差导致的测距修正误差,将数倍放大于卫星定轨结果中,若想得到优于20米的卫星位置精度,要求定轨站钟差精度优于20ns。多星动力学定轨每小时运行一次,加之定轨解算所需的时间,需要对多星多站钟差进行最长1.5小时的预报,然而地面定轨站时频基准本身没有长期预报需求,通常选用晶振或性能较差的原子钟,钟差长期预报误差往往可达数百甚至数千ns,最大可导致数百米的卫星位置误差。采用性能更好的氢原子钟也可实现定轨站钟差的高精度长期预报,但是氢原子钟价格过高,经济代价大。
四是可选定轨站少。由于多星动力学定轨计算量会随着定轨站数量增加呈指数级上升,因此考虑到实时性要求,多星动力学定轨通常仅选取少量定轨站。多站支持条件下的动力学定轨理论上,最少需要4个定轨站,一旦个别定轨站产生故障,则可能无法找到精度符合要求的定轨站进行替代,导致定轨程序中断。
发明内容
本发明的实施例提供了一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法,以实现对GEO卫星进行有效地定位。
为了实现上述目的,本发明采取了如下技术方案。
一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法,构建由中心基准站、若干定轨站和GEO卫星组成的GEO卫星运动学定轨系统,中心基准站包括时频基准系统、星地时间同步天线和基准接收机,定轨站包括卫星导航接收机,所述方法包括:
通过多星共视法实时测量各个定轨站与中心基准站之间的钟差;
利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差;
根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨。
优选地,所述的通过多星共视法实时测量各个定轨站与中心基准站之间的钟差,包括:
在中心基准站和若干定轨站部署相同型号的接收机,中心基准站和定轨站在异地同时观测同一组卫星,通过站间单差法将定轨站接收机的观测数据与中心基准站接收机的观测数据做差,计算得到定轨站相对于中心基准站的钟差,根据为已知值的中心基准站的钟差和定轨站相对于中心基准站的钟差得到定轨站的钟差。
优选地,所述的通过多星共视法实时测量各个定轨站与中心基准站之间的钟差,包括:
中心基准站和定轨站的卫星导航伪距观测方程为:
其中,为伪距观测量,s代表卫星号,r代表定轨站号,j代表频率号,/>为星地距离,根据卫星和接收机位置得到,c为光速,dtr和dts为中心基准站或定轨站接收机和卫星钟差,δtr为对流层延迟,δion为电离层延迟,δph为卫星相位中心影响,δrel为相对论效应影响,δecc接收机相位中心影响,δtide为潮汐影响,dr,j和/>分别为接收机端和卫星端伪距硬件延迟偏差,εp为未模型化误差;
利用同一时刻,同一接收机对同一卫星的1号频率和3号频率的观测数据和/>构建消电离层组合,利用电离层延迟量δion,j与观测频率fj的关系,消除一阶电离层延迟:
其中,IF标识为消电离层组合,消电离层组合的观测量,dr,IF和/>分别为接收机端和卫星端伪距消电离层组合硬件延迟偏差,fj为频率值;
在某一时刻,中心基准站r1与定轨站r2对卫星s形成共视,得到两站的双频消电离层组合和/>
对式(3)和式(4)做差,得到站间单差
其中,dts、δph、被消除;中心基准站和定轨站位置为已知值,中心基准站r1的钟差/>为已知值;卫星位置通过广播星历参数计算得到,经过站间单差消除/>中由广播星历引起的误差;
在同一时刻,中心基准站r1与定轨站r2有n个共视卫星,得到至
共n个站间单差观测值,将观测方程改写矩阵的形式:
V=AX+l (8)
X=dtr2 (9)
其中,i为卫星号,A参数矩阵,X为待估参数定轨站r2的接收机钟差dtr2,V为观测误差;
解算法方程得到定轨站接收机与中心基准站之间的钟差的估计值
V=l-A(ATPA)-1ATPl (14)。
优选地,所述的利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差,包括:
GEO卫星和中心基准站同时分别向对方发射伪距测距信号,GEO卫星和中心基准站接收对方发送的测距信号,并记录观测值,卫星通过数据链路将观测值下传,通过对同一时刻的星地双向观测数据进行做差,计算出GEO卫星与中心基准站之间的钟差值;
假设中心基准站r1和GEO卫星s进行星地双向观测,中心基准站r1在t0时刻测量得到下行观测值GEO卫星s在同一时刻测量得到上行观测值/>则有:
其中,为t0时刻下行观测值,/>为t0时刻上行观测值,/>是信号在空间中传播的真实距离,ε为未模型化误差和观测噪声;
将式(16)减去式(15)除以2,得到GEO卫星s的钟差dts:
优选地,所述的根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨,包括:
设在全球导航卫星系统GNSS数据处理中有多个地面站伪距信号数据,利用多站交会得到卫星位置信息;
将公式(2)改写为以下格式:
其中,卫星位置坐标(Xs,Ys,Zs)为待估参数,定轨站坐标(xr,yr,zr)为已知值,轨道机动卫星钟差dts由星地双向时间同步得到,定轨站钟差dtr由共视法钟差测量得到;
设在定轨时,同时有n个测站r1至rn,跟踪观测轨道机动卫星,得到至/>共n个观测量,对伪距观测方程进行线性化,得到:
其中,i为定轨站号,Gi(·)为待估参数与观测数据的函数,为卫星位置向量(Xs,Ys,Zs),/>为初始参考位置,改写为矩阵形式:
V=AδX+l (20)
n为定轨站数量,解算方程(21)得到实时卫星位置
由上述本发明的实施例提供的技术方案可以看出,本发明方法利用共视法授时获得地面站的精确时间,通过多个定轨站同时测距,进行运动学定轨,获得为卫星精确位置,为地面天线对卫星连续跟踪提供目标指引,确保GEO卫星上承载的星基增强服务、精密单点定位服务和短报文服务连续稳定。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种GEO卫星运动学定轨系统的结构图;
图2为本发明实施例提供的一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法的处理流程图;
图3为本发明实施例提供的一种多星共视法钟差测量工作场景示意图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
针对基于多星动力学定轨钟差预报修正伪距进行运动学定轨方法定轨精度低、可选定轨站少的问题,本发明实施例提出了基于共视法授时的GEO卫星轨道机动期间位置测定方法,利用共视法实时测定地面定轨站与主控站BDT钟差,利用实时解算钟差修正伪距观测量,进而实现GEO卫星轨道机动期间的运动学定轨,解决基于多星动力学定轨钟差预报计算量大、实时性差、精度低、可选定轨站少的问题,提高GEO卫星轨道机动期间定轨算法精度、可靠性和稳定性。
本发明实施例提供的一种GEO卫星运动学定轨系统的结构图如图1所示,该系统由一个中心基准站、若干定轨站和轨道机动卫星组成。中心基准站包括时频基准系统、星地时间同步天线和基准接收机,定轨站由卫星导航接收机组成。
基于图1所示的系统,本发明实施例提供的一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法的处理流程如图2所示,包括如下的处理步骤:
步骤S10、通过多星共视法实时测量各个定轨站与中心基准站之间的钟差;
步骤S20、利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差;
步骤S30、根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨。
具体的,上述步骤S10包括:图3为本发明实施例提供的一种多星共视法钟差测量工作场景示意图。共视法钟差测量的本质是在中心基准站和若干定轨站部署相同型号的接收机,中心基准站和定轨站在异地同时观测同一组卫星,通过站间单差法,将定轨站接收机的观测数据与中心基准站接收机的观测数据做差,计算得到定轨站相对于中心基准站的钟差,中心基准站的钟差为已知值,进而得到定轨站的钟差。该方法利用卫星导航测量误差的强时空相关性,消除共同误差,得到较高的钟差估计精度。
中心基准站和定轨站的卫星导航伪距观测方程为:
其中,为伪距观测量,s代表卫星号,r代表定轨站号,j代表频率号,/>为星地距离,根据卫星和接收机位置得到,c为光速,dtr和dts为中心基准站或定轨站接收机和卫星钟差,δtr为对流层延迟,δion为电离层延迟,δph为卫星相位中心影响,δrel为相对论效应影响,δecc接收机相位中心影响,δtide为潮汐影响,dr,j和/>分别为接收机端和卫星端伪距硬件延迟偏差,εp为多路径、观测噪声等未模型化误差。
利用同一时刻,同一接收机对同一卫星的1号频率和3号频率的观测数据和构建消电离层组合,利用电离层延迟量δion,j与观测频率fj的关系,消除一阶电离层延迟:
其中,IF标识为消电离层组合,消电离层组合的观测量,dr,IF和/>分别为接收机端和卫星端伪距消电离层组合硬件延迟偏差,fj为频率值。
在某一时刻,中心基准站r1与定轨站r2对卫星s形成共视,可同时得到两站的双频消电离层组合和/>
对式(3)和式(4)做差,可得到站间单差
其中,dts、δph、被消除;中心基准站和定轨站位置为已知值,中心基准站r1的钟差/>为已知值;卫星位置可以通过广播星历参数计算得到,经过站间单差,可基本消除中由广播星历引起的误差;站间单差得到的对流层延迟Δδtr残差、接收机中心影响Δδecc残差、相对论效应影响Δδrel残差和潮汐影响Δδtide残差,先利用各自模型进行误差修正,再通过做差消除公共部分,剩余影响可以忽略不计。
在同一时刻,中心基准站r1与定轨站r2有n个共视卫星,可得到至/>共n个站间单差观测值,将观测方程改写矩阵的形式:
V=AX+l (8)
X=dtr2 (9)
其中,i为卫星号,A参数矩阵,X为待估参数定轨站r2的接收机钟差dtr2,V为观测误差。
解算法方程,可得定轨站接收机钟差估计值
V=l-A(ATPA)-1ATPl (14)
具体的,上述步骤S20包括:利用星地双向时间同步法测量卫星钟差,其基本原理为:卫星和中心基准站同时分别向对方发射伪距测距信号,与此同时,卫星和中心基准站接收对方发送的测距信号,并记录观测值,卫星通过数据链路将观测值下传,通过对同一时刻的星地双向观测数据进行做差,可以计算星地钟差值。
假设中心基准站r1和卫星s进行星地双向观测,中心基准站r1在t0时刻,测量得到下行观测值卫星s在同一时刻,测量得到上行观测值/>则有:
其中,为t0时刻下行观测值,/>为t0时刻上行观测值,/>是信号在空间中传播的真实距离,ε为未模型化误差和观测噪声,为了简洁,省去了电离层延迟、对流层延迟等误差,其与公式(1)相同,可用模型改正,且上下行基本相同,做差后可以消除。
将式(16)减去式(15)除以2,可得GEO卫星钟差dts
公式(17)消除了对流层延迟、卫星星历误差和中心基准站坐标误差等误差的影响,与信号频率有关的电离层延迟也被很大程度地削弱。
具体的,上述步骤S30包括:卫星运动学定轨。在GNSS(Global NavigationSatellite System,全球导航卫星系统)数据处理中,有多个地面站伪距信号数据,因此,可以直接利用多站交会得到卫星位置信息,也可进一步利用位置信息差分得到卫星速度。
将公式(2)改写为以下格式:
其中,卫星位置坐标(Xs,Ys,Zs)为待估参数,定轨站坐标(xr,yr,zr)为已知值,轨道机动卫星钟差dts由星地双向时间同步得到,定轨站钟差dtr由共视法钟差测量得到。
在定轨时,同时有n个测站(r1至rn),跟踪观测轨道机动卫星,同时可以得到至/>共n个观测量,对伪距观测方程进行线性化。
其中,i为定轨站号,Gi(·)为待估参数与观测数据的函数,为卫星位置向量(Xs,Ys,Zs),/>为初始参考位置,改写为矩阵形式:
V=AδX+l (20)
其中,
n为定轨站数量,解算方程可得,
得到实时卫星位置。
综上所述,本发明实施例利用共视法计算定轨站钟差,进而进行几何法定轨,本方法有以下优点:
第一,避免了进行站钟差预报,从根本上解决了定轨站原子钟长期稳定性较差,导致钟差预报精度低的问题,有效提高了钟差计算的精度。
第二,通过定轨算法的优化,避免了对于硬件的大规模更新,性价比较高。
第三,通过将定轨站钟差计算与定轨的解绑,大大扩展了几何法定轨定轨站的可选范围,避免了个别定轨站故障对于定轨稳定性的影响,提高了定轨的稳定性和可靠性。
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或系统实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及系统实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种基于共视法授时的GEO卫星轨道机动期间轨道测定方法,其特征在于,构建由中心基准站、若干定轨站和GEO卫星组成的GEO卫星运动学定轨系统,中心基准站包括时频基准系统、星地时间同步天线和基准接收机,定轨站包括卫星导航接收机,所述方法包括:
通过多星共视法实时测量各个定轨站与中心基准站之间的钟差;
利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差;
根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨;
所述的通过多星共视法实时测量各个定轨站与中心基准站之间的钟差,包括:
中心基准站和定轨站的卫星导航伪距观测方程为:
其中,为伪距观测量,s代表卫星号,r代表定轨站号,j代表频率号,/>为星地距离,根据卫星和接收机位置得到,c为光速,dtr和dts为中心基准站或定轨站接收机和卫星钟差,δtr为对流层延迟,δion,j为电离层延迟量,δph为卫星相位中心影响,δrel为相对论效应影响,δecc接收机相位中心影响,δtide为潮汐影响,dr,j和/>分别为接收机端和卫星端伪距硬件延迟偏差,εp为未模型化误差;
利用同一时刻,同一接收机对同一卫星的1号频率和2号频率的观测数据和/>构建消电离层组合,利用电离层延迟量δion,j与观测频率fj的关系,消除一阶电离层延迟:
其中,IF标识为消电离层组合,消电离层组合的观测量,dr,IF和/>分别为接收机端和卫星端伪距消电离层组合硬件延迟偏差,fj为频率值;
在某一时刻,中心基准站r1与定轨站r2对卫星s形成共视,得到两站的双频消电离层组合和/>
对式(3)和式(4)做差,得到站间单差
其中,dts、δph、被消除;中心基准站和定轨站位置为已知值,中心基准站r1的钟差/>为已知值;卫星位置通过广播星历参数计算得到,经过站间单差消除/>中由广播星历引起的误差;
在同一时刻,中心基准站r1与定轨站r2有n个共视卫星,得到至/>共n个站间单差观测值,将所述卫星导航伪距观测方程改写矩阵的形式:
V=AX+l (8)
X=dtr2 (9)
其中,i为卫星号,A参数矩阵,X为待估参数定轨站r2的接收机钟差dtr2,V为观测误差;
求解方程(8),得到定轨站接收机与中心基准站之间的钟差的估计值
V=l-A(ATPA)-1ATPl (14)
所述的根据定轨站的位置信息、各个定轨站与中心基准站之间的钟差和GEO卫星与中心基准站之间的钟差,进行GEO卫星的运动学定轨,包括:
设在全球导航卫星系统GNSS数据处理中有多个地面站伪距信号数据,利用多站交会得到卫星位置信息;
将公式(2)改写为以下格式:
其中,卫星位置坐标(Xs,Ys,Zs)为待估参数,定轨站坐标(xr,yr,zr)为已知值,轨道机动卫星钟差dts由星地双向时间同步得到,定轨站钟差dtr由共视法钟差测量得到;
设在定轨时,同时有n个测站r1至rn,跟踪观测轨道机动卫星,得到至/>共n个观测量,对伪距观测方程进行线性化,得到:
其中,i为定轨站号,Gi(·)为待估参数与观测数据的函数,为卫星位置向量(Xs,Ys,Zs),/>为初始参考位置,改写为矩阵形式:
V=AδX+l (20)
n为定轨站数量,解算方程(21)得到实时卫星位置
2.根据权利要求1所述的方法,其特征在于,所述的通过多星共视法实时测量各个定轨站与中心基准站之间的钟差,包括:
在中心基准站和若干定轨站部署相同型号的接收机,中心基准站和定轨站在异地同时观测同一组卫星,通过站间单差法将定轨站接收机的观测数据与中心基准站接收机的观测数据做差,计算得到定轨站相对于中心基准站的钟差,根据为已知值的中心基准站的钟差和定轨站相对于中心基准站的钟差得到定轨站的钟差。
3.根据权利要求2所述的方法,其特征在于,所述的利用双向时间同步天线实时测量GEO卫星与中心基准站之间的钟差,包括:
GEO卫星和中心基准站同时分别向对方发射伪距测距信号,GEO卫星和中心基准站接收对方发送的测距信号,并记录观测值,卫星通过数据链路将观测值下传,通过对同一时刻的星地双向观测数据进行做差,计算出GEO卫星与中心基准站之间的钟差值;
假设中心基准站r1和GEO卫星s进行星地双向观测,中心基准站r1在t0时刻测量得到下行观测值GEO卫星s在同一时刻测量得到上行观测值/>则有:
其中,为t0时刻下行观测值,/>为t0时刻上行观测值,/>是信号在空间中传播的真实距离,ε为未模型化误差和观测噪声;
将式(16)减去式(15)除以2,得到GEO卫星s的钟差dts:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211266864.9A CN115639582B (zh) | 2022-10-17 | 2022-10-17 | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211266864.9A CN115639582B (zh) | 2022-10-17 | 2022-10-17 | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115639582A CN115639582A (zh) | 2023-01-24 |
CN115639582B true CN115639582B (zh) | 2023-11-17 |
Family
ID=84944097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211266864.9A Active CN115639582B (zh) | 2022-10-17 | 2022-10-17 | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115639582B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116931027B (zh) * | 2023-09-15 | 2023-11-28 | 成都子辰时频科技有限公司 | 一种基于sbas的动态高精度卫星共视算法系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN108931915A (zh) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | 利用导航卫星的授时方法和装置、计算机可读存储介质 |
EP2741108B1 (en) * | 2012-12-07 | 2019-08-28 | Technische Universität München | Method for determining the position of a high orbit satellite |
WO2020133711A1 (zh) * | 2018-12-28 | 2020-07-02 | 上海海积信息科技股份有限公司 | 一种卫星定轨方法、装置及电子设备 |
CN111505677A (zh) * | 2020-04-15 | 2020-08-07 | 中国科学院国家授时中心 | 一种基于地面参考站观测的geo卫星轨道机动修复方法 |
CN112924992A (zh) * | 2021-01-25 | 2021-06-08 | 中国科学院国家授时中心 | 一种geo轨道精度评估方法、装置、电子设备及存储介质 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
CN113885055A (zh) * | 2021-09-07 | 2022-01-04 | 中国科学院上海天文台 | 星地双向钟差约束的导航卫星多星定轨系统及方法 |
-
2022
- 2022-10-17 CN CN202211266864.9A patent/CN115639582B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2741108B1 (en) * | 2012-12-07 | 2019-08-28 | Technische Universität München | Method for determining the position of a high orbit satellite |
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN108931915A (zh) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | 利用导航卫星的授时方法和装置、计算机可读存储介质 |
WO2020133711A1 (zh) * | 2018-12-28 | 2020-07-02 | 上海海积信息科技股份有限公司 | 一种卫星定轨方法、装置及电子设备 |
CN111505677A (zh) * | 2020-04-15 | 2020-08-07 | 中国科学院国家授时中心 | 一种基于地面参考站观测的geo卫星轨道机动修复方法 |
CN112924992A (zh) * | 2021-01-25 | 2021-06-08 | 中国科学院国家授时中心 | 一种geo轨道精度评估方法、装置、电子设备及存储介质 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
CN113885055A (zh) * | 2021-09-07 | 2022-01-04 | 中国科学院上海天文台 | 星地双向钟差约束的导航卫星多星定轨系统及方法 |
Non-Patent Citations (7)
Title |
---|
Common-View Time Transfer Using Geostationary Satellite;Wei Pei等;IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL;全文 * |
GEO导航卫星定轨观测模型;唐波;张军;陈刘成;;测绘科学技术学报(第05期);全文 * |
Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver;Fen Cao等;Advances in Space Research;全文 * |
北斗卫星超短弧运动学定轨方法优化与试验分析;郭睿;陈金平;朱陵凤;唐桂芬;吴晓莉;;测绘学报(第04期);全文 * |
基于自发自收测距的GEO卫星精密定轨;郭睿;刘雁雨;谭红力;何峰;董恩强;;测绘科学技术学报(第05期);全文 * |
监测站时间同步条件下的单星定轨;宋小勇;贾小林;毛悦;吴显兵;;测绘科学(第05期);全文 * |
钟差和动力学参数联合约束下的北斗卫星轨道快速确定;郭睿;陈金平;李晓杰;常志巧;;空间科学学报(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115639582A (zh) | 2023-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110231037B (zh) | 一种附加钟差模型约束的gnss机动卫星定轨方法 | |
Li et al. | LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services | |
CN109061696B (zh) | 一种确定导航卫星轨道和钟差的方法 | |
Yang et al. | Inter-satellite link enhanced orbit determination for BeiDou-3 | |
CN103033188B (zh) | 基于综合孔径观测的导航卫星自主时间同步方法 | |
CN108061908B (zh) | 基于低轨卫星星载gnss技术的slr台站三维坐标几何解算方法 | |
EP2406652A1 (en) | Removing biases in dual frequency gnss receivers using sbas | |
Li et al. | Review of PPP–RTK: Achievements, challenges, and opportunities | |
CN113885055A (zh) | 星地双向钟差约束的导航卫星多星定轨系统及方法 | |
Sibthorpe et al. | An evaluation of solar radiation pressure strategies for the GPS constellation | |
Li et al. | Precise orbit determination for the FY-3C satellite using onboard BDS and GPS observations from 2013, 2015, and 2017 | |
Wang et al. | Evaluating the impact of CNES real-time ionospheric products on multi-GNSS single-frequency positioning using the IGS real-time service | |
CN104316943A (zh) | 一种伪距离和多普勒组合差分定位系统及方法 | |
Guo et al. | Precise orbit determination and rapid orbit recovery supported by time synchronization | |
CN115639582B (zh) | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 | |
Li et al. | Orbit determination and prediction for Beidou GEO satellites at the time of the spring/autumn equinox | |
Qing et al. | Benefits of satellite clock modeling in BDS and Galileo orbit determination | |
CN113608247A (zh) | 一种卫星的定轨方法 | |
Ye et al. | Analysis of BDS-3 distributed autonomous navigation based on BeiDou system simulation platform | |
CN104309817B (zh) | 基于多台并址接收机的北斗导航卫星区域定轨方法 | |
Yu et al. | Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products | |
Liu et al. | Beidou wide-area augmentation system clock error correction and performance verification | |
CN112731504B (zh) | 对月球探测器自主定轨的方法及装置 | |
Zhang et al. | Precise orbit determination for TH02-02 satellites based on BDS3 and GPS observations | |
Li et al. | Performance assessment of real-time orbit determination for the Haiyang-2D using onboard BDS-3/GPS observations |
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 |