CN113624243B - 一种近地轨道卫星的星上实时轨道预报方法 - Google Patents

一种近地轨道卫星的星上实时轨道预报方法 Download PDF

Info

Publication number
CN113624243B
CN113624243B CN202110941066.0A CN202110941066A CN113624243B CN 113624243 B CN113624243 B CN 113624243B CN 202110941066 A CN202110941066 A CN 202110941066A CN 113624243 B CN113624243 B CN 113624243B
Authority
CN
China
Prior art keywords
satellite
orbit
earth
low earth
data
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
Application number
CN202110941066.0A
Other languages
English (en)
Other versions
CN113624243A (zh
Inventor
袁俊军
唐成盼
李凯
巩秀强
周善石
曹月玲
胡小工
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Astronomical Observatory of CAS
Original Assignee
Shanghai Astronomical Observatory of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Astronomical Observatory of CAS filed Critical Shanghai Astronomical Observatory of CAS
Priority to CN202110941066.0A priority Critical patent/CN113624243B/zh
Publication of CN113624243A publication Critical patent/CN113624243A/zh
Application granted granted Critical
Publication of CN113624243B publication Critical patent/CN113624243B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • 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/14—Receivers specially adapted for specific applications
    • 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

Abstract

本发明涉及卫星轨道预报技术领域，提出一种近地轨道卫星的星上实时轨道预报方法，包括由近地轨道卫星接收数据，并且对所述数据进行内容检验以及格式转换；由近地轨道卫星对所述数据进行预处理；以及由近地轨道卫星根据预处理后的数据确定预报轨道，其中包括由近地轨道卫星执行下列动作：根据预处理后的数据进行精密单点定位解算以确定离散轨道点；对所述离散轨道点进行拟合以确定初始轨道；计算地球中心引力以及非球形地球重力场摄动；以及根据所述初始轨道、地球中心引力以及非球形地球重力场摄动进行轨道积分外推以确定预报轨道。

Description

一种近地轨道卫星的星上实时轨道预报方法
技术领域
本发明总的来说涉及卫星轨道预报技术领域。具体而言，本发明涉及一种近地轨道卫星的星上实时轨道预报方法。
背景技术
近地轨道(low earth orbit LEO)卫星通常运行于300-2000km的近地球轨道，其飞行速度快，相对于地面监测站的图形几何构型变化快，是天然的全球卫星导航系统(Global Navigation Sate ll i te System GNSS)监测站和导航信号播发平台。相比高地球轨道的GNSS卫星，LEO卫星所播发的导航信号传播过程中具有损失更小、落地功率大、并且拥有较快的图形变化结构的优点，有利于地面精密单点定位的模糊度快速收敛。因此，基于大规模LEO星座构建低轨导航增强系统成为导航增强领域的研究热点。
为了保证LEO卫星提供实时高精度的定位、导航、授时(Positioning、Navigation、Timing PNT)服务，必须实现LEO卫星星上实时高精度轨道预报。目前通常使用的LEO卫星轨道预报方法包括多项式拟合法、轨道根数分析法以及动力学拟合法。然而，多项式拟合法及轨道根数分析法精度过低，无法满足低轨导航增强系统的轨道预报精度需求。而动力学拟合法尽管精度有所保障，但需依赖强大的计算能力，不适用于处理能力十分有限的LEO卫星星载处理器，只能将星载观测数据下传到地面再进行计算，无法满足实时性要求。
发明内容
为至少部分解决现有技术中的上述问题，本发明提出一种近地轨道卫星的星上实时轨道预报方法，包括下列步骤：
由近地轨道卫星接收数据，并且对所述数据进行内容检验以及格式转换；
由近地轨道卫星对所述数据进行预处理；以及
由近地轨道卫星根据预处理后的数据确定预报轨道，其中包括由近地轨道卫星执行下列动作：
根据预处理后的数据进行精密单点定位解算以确定离散轨道点；
对所述离散轨道点进行拟合以确定初始轨道；
计算地球中心引力以及非球形地球重力场摄动；以及
根据所述初始轨道、地球中心引力以及非球形地球重力场摄动进行轨道积分外推以确定预报轨道。
在本发明一个实施例中规定，所述近地轨道卫星具有星载全球卫星导航系统接收机，其中所述星载全球卫星导航系统接收机接收全球卫星导航系统卫星发出的数据。
在本发明一个实施例中规定，其特征在于，近地轨道卫星接收的数据包括：近地轨道星载相位、近地轨道伪距、全球卫星导航系统轨道数据、全球卫星导航系统钟差、地球定向参数以及卫星天线相位中心改正值。
在本发明一个实施例中规定，所述数据内容检验包括数据容量检验、数据日期检验以及数据来源检验。
在本发明一个实施例中规定，由近地轨道卫星对所述数据进行预处理包括：
进行电离层延迟修正以确定电离层延迟误差；
进行全球卫星导航系统卫星/近地轨道卫星的天线相位中心修正以确定全球卫星导航系统卫星/近地轨道卫星的天线相位中心偏差；以及
进行相对论修正以确定相对论延迟误差。
在本发明一个实施例中规定，所述全球卫星导航系统卫星/近地轨道卫星的天线相位中心修正包括将全球卫星导航系统卫星/近地轨道卫星的天线相位中心与全球卫星导航系统卫星/近地轨道卫星的质心之间的偏差α从卫星星固坐标系转换到J2000惯性坐标系中，表示为下式：
α′＝R(t)α
R(t)＝[g1，g2，g3]T
Figure BDA0003213595240000021
g3＝g1×g2
其中，R(t)表示在观测历元t时寸刻卫星星固坐标系与J2000惯性坐标系之间的转换矩阵，α′表示J2000惯性坐标系下全球卫星导航系统卫星/近地轨道卫星的天线相位中心与全球卫星导航系统卫星/近地轨道卫星的质心之间的偏差，r(t)、r′(t)分别表示全球卫星导航系统卫星/近地轨道卫星在J2000惯性坐标系下的位置和速度。
在本发明一个实施例中规定，所述相对论修正包括计算狭义相对论效应周期性误差并且进行修正，表示为下式：
Figure BDA0003213595240000031
其中，Δρrel表示狭义相对论效应周期性误差，
Figure BDA0003213595240000032
及
Figure BDA0003213595240000033
表示近地轨道卫星的坐标及速度，c表示光速，
Figure BDA0003213595240000034
及
Figure BDA0003213595240000035
表示全球卫星导航系统卫星的坐标及速度。
在本发明一个实施例中规定，由近地轨道卫星根据预处理后的数据进行精密单点定位解算包括下列步骤：
由近地轨道卫星根据预处理后的数据构建双频消电离层方程，表示为下式：
Figure BDA0003213595240000036
Figure BDA0003213595240000037
其中，角标r、s分别表示星载全球卫星导航系统接收机r以及全球卫星导航系统卫星s，
Figure BDA0003213595240000038
伪距观测值，
Figure BDA0003213595240000039
表示星载相位观测值，
Figure BDA00032135952400000310
表示全球卫星导航系统卫星s至星载全球卫星导航系统接收机r的几何距离，c表示光速，δtr表示星载全球卫星导航系统接收机钟差，δts表示全球卫星导航系统卫星钟差，δρion表示电离层延迟误差，δρrel表示相对论延迟误差，δρpco表示星载全球卫星导航系统接收机天线相位中心偏差，
Figure BDA00032135952400000311
表示全球卫星导航系统卫星天线相位中心偏差，λ表示波长，N表示相位整周模糊度，ε表示观测噪声；
将所述双频消电离层方程在近似坐标(x0，y0，z0)处展开，表示为下式：
Figure BDA00032135952400000312
其中，ρ表示近似展开项，xsat、ysat、zsat表示GNSS卫星坐标，ρ0表示GNSS卫星与LEO卫星的近似几何距离；
基于伪距观测值对同一历元下所有观测到的GNSS卫星构建伪距定位方程，表示为下式：
Figure BDA0003213595240000041
基于星载相位观测值，并且添加模糊度参数以构建相位定位方程；以及
通过卡尔曼滤波参数估计方法解算所述伪距定位方程和相位定位方程以确定离散轨道点。
在本发明一个实施例中规定，由近地轨道卫星计算地球中心引力表示为下式：
Figure BDA0003213595240000042
该式中，
Figure BDA0003213595240000043
表示地球中心引力，G表示万有引力常数，ME表示地球质量，
Figure BDA0003213595240000044
表示近地轨道卫星与地心的距离矢量，r表示近地轨道卫星与地心的距离数值，m表示近地轨道卫星质量；以及
由近地轨道卫星计算非球形地球重力场摄动，其中包括将所述非球形地球重力场摄动进行归一化处理，归一化处理后的摄动函数T表示为下式：
Figure BDA0003213595240000045
该式中，GM表示地球引力常数，r表示地心向经，φ表示地心纬度，λ表示地心经度，a表示地球赤道平均半径，
Figure BDA0003213595240000046
表示规格化的位系数，
Figure BDA0003213595240000047
表示n阶m次的规格化勒让德函数，N表示地球引力位的最大阶次。
在本发明一个实施例中规定，本发明方法还包括对所述预报轨道进行精度检核，其中当所述预报轨道符合精度要求时由近地轨道卫星输出所述预报轨道，当所述预报轨道不符合精度要求时由近地轨道卫星重复根据预处理后的数据确定预报轨道的步骤。
本发明至少具有如下有益效果：本发明提出一种算法简洁高效，并且能够满足高精度需求的LE0星上实时轨道预报方法，避免了复杂的动力学信息计算负担，相对于传统动力学拟合法计算所需资源更少，十分适合用于星上实时计算；本发明保留部分重要动力学信息，相对于多项式拟合法、轨道根数分析法等精度更高；本发明的计算效率高，利用本发明方法预报轨道平均速度快；以及本发明的计算可靠性高。
附图说明
为进一步阐明本发明的各实施例中具有的及其它的优点和特征，将参考附图来呈现本发明的各实施例的更具体的描述。可以理解，这些附图只描绘本发明的典型实施例，因此将不被认为是对其范围的限制。在附图中，为了清楚明了，相同或相应的部件将用相同或类似的标记表示。
图1示出了本发明一个实施例中近地轨道卫星进行星上实时轨道预报的流程示意图。
图2示出了本发明一个实施例中近地轨道卫星进行精密单点定位解算的定位结果精度示意图。
图3a-c示出了本发明一个实施例中GRACE-C卫星进行轨道预报的结果示意图。
图4a-c示出了本发明一个实施例中HY2A卫星进行轨道预报的结果示意图。
具体实施方式
应当指出，各附图中的各组件可能为了图解说明而被夸大地示出，而不一定是比例正确的。在各附图中，给相同或功能相同的组件配备了相同的附图标记。
在本发明中，除非特别指出，“布置在…上”、“布置在…上方”以及“布置在…之上”并未排除二者之间存在中间物的情况。此外，“布置在…上或上方”仅仅表示两个部件之间的相对位置关系，而在一定情况下、如在颠倒产品方向后，也可以转换为“布置在…下或下方”，反之亦然。
在本发明中，各实施例仅仅旨在说明本发明的方案，而不应被理解为限制性的。
在本发明中，除非特别指出，量词“一个”、“一”并未排除多个元素的场景。
在此还应当指出，在本发明的实施例中，为清楚、简单起见，可能示出了仅仅一部分部件或组件，但是本领域的普通技术人员能够理解，在本发明的教导下，可根据具体场景需要添加所需的部件或组件。另外，除非另行说明，本发明的不同实施例中的特征可以相互组合。例如，可以用第二实施例中的某特征替换第一实施例中相对应或功能相同或相似的特征，所得到的实施例同样落入本申请的公开范围或记载范围。
在此还应当指出，在本发明的范围内，“相同”、“相等”、“等于”等措辞并不意味着二者数值绝对相等，而是允许一定的合理误差，也就是说，所述措辞也涵盖了“基本上相同”、“基本上相等”、“基本上等于”。以此类推，在本发明中，表方向的术语“垂直于”、“平行于”等等同样涵盖了“基本上垂直于”、“基本上平行于”的含义。
另外，本发明的各方法的步骤的编号并未限定所述方法步骤的执行顺序。除非特别指出，各方法步骤可以以不同顺序执行。
下面结合具体实施方式参考附图进一步阐述本发明。
如图1所示，在本发明一个实施例中提出一个针对低轨导航增强系统的LEO卫星星上实时轨道预报方法，其中可以包括下列步骤：
步骤一、LEO卫星接收数据、并且进行内容检验与格式转换。
LEO卫星进行星上轨道预报所需接收的数据主要包括LEO星载相位、LEO伪距数据、GNSS轨道、GNSS钟差数据、以及地球定向参数(Earth 0rientation Parameters，EOP)和卫星天线相位中心改正值等辅助数据。
在数据接收后LEO卫星需要进行数据内容检验及格式转换。内容检验例如可以包括数据容量检验、数据日期检验以及数据来源检验，其中通过所述内容检验可以保证LEO卫星输入数据合法有效。格式转换则是将原始的数据格式转换为处理内部可识别的格式，并根据不同的数据类型分别将转换后的数据送至LEO卫星的预处理模块。
步骤二、LEO卫星进行数据预处理。
LEO卫星进行精密单点定位解算(LEO-PPP)需要基于步骤一中所接收的数据，然而这些数据会受到电离层延迟、多路径等误差的影响，因此必须对其进行数据预处理以修正公共误差，其中所述数据预处理例如可以包括电离层延迟修正、GNSS/LE0卫星的天线相位中心修正以及相对论修正。
需要进行GNSS/LE0卫星的天线相位中心修正是由于GNSS/LE0卫星天线相位中心与GNSS/LE0卫星质心之间存在偏差。GNSS/LE0卫星的天线相位中心修正包括将GNSS/LE0卫星天线相位中心与GNSS/LE0卫星质心之间的偏差α从卫星星固坐标系转换到J2000惯性坐标系中，表示为下式：
α′＝R(t)α
R(t)＝[g1，g2，g3]T
Figure BDA0003213595240000071
g3＝g1×g2
其中，R(t)表示在观测历元t时寸刻卫星星固坐标系与J2000惯性坐标系之间的转换矩阵，α′表示J2000惯性坐标系下GNSS/LE0卫星天线相位中心与GNSS/LE0卫星质心之间的偏差，r(t)、r′(t)表示GNSS/LEO卫星在J2000惯性坐标系下的位置和速度。
需要进行相对论修正是由于GNSS卫星以及LEO卫星轨道高度差异较大，因此必须考虑相对论效应对GNSS/LEO钟差的影响。相对论效应误差包括相对论效应常量误差和周期误差，其中，常量部分通常在卫星发射前通过人工干预消除，而相对论效应周期性误差可以通过下式计算：
Figure BDA0003213595240000072
其中，Δρrel表示狭义相对论效应周期性误差，
Figure BDA0003213595240000073
及
Figure BDA0003213595240000074
表示LEO卫星的坐标及速度，c表示光速，
Figure BDA0003213595240000075
及
Figure BDA0003213595240000076
表示GNSS卫星的坐标及速度。
步骤三、LEO卫星根据预处理后的数据进行LEO-PPP解算以获得离散轨道点，其中包括下列步骤：
利用预处理后的星载GNSS接收机r接收到GNSS卫星s的伪距、相位观测量，LEO卫星可以构建双频消电离层方程，表示为下式：
Figure BDA0003213595240000077
Figure BDA0003213595240000078
其中，
Figure BDA0003213595240000079
表示伪距观测值，
Figure BDA00032135952400000710
表示相位观测值，
Figure BDA00032135952400000711
表示GNSS卫星s至星载接收机r的几何距离，c表示光速，δtr表示星载GNSS接收机钟差，δts表示GNSS卫星钟差，δρion表示电离层延迟误差，δρrel表示相对论延迟误差，δρpco表示星载GNSS接收机天线相位中心偏差，
Figure BDA00032135952400000712
表示GNSS卫星天线相位中心偏差，λ表示波长，N表示相位整周模糊度，ε表示观测噪声。
将所述双频消电离层方程在近似坐标(x0，y0，z0)处展开，表示为下式：
Figure BDA0003213595240000081
其中，xsat、ysat、zsat表示GNSS卫星坐标，ρ0表示GNSS卫星与LEO卫星的近似几何距离。
使用伪距观测数据可以对同一历元下所有观测到的GNSS卫星建立伪距定位方程，表示为下式：
Figure BDA0003213595240000082
所述伪距定位方程仅使用伪距观测数据，使用相位数据时需要添加一个模糊度参数即可以相位定位方程。
以及在建立伪距、相位定位方程后，可以通过卡尔曼滤波参数估计方法进行解算，便可得到LEO-PPP定位结果，其中所获得的定位结果为离散轨道点。本发明一个实施例中LEO卫星在一天内的LEO-PPP定位结果精度如图2所示。
步骤四、LEO卫星对所述离散轨道点进行拟合以获取初始轨道。
为获取指定时刻的初始轨道对所述离散轨道点进行拟合，其中可以采用最小二乘拟合方法，以获取所述指定时刻的坐标及速度信息。
步骤五、LEO卫星计算地球中心引力以及非球形摄动力。
读取重力场文件，计算摄动力量级最大的地球中心引力及非球形摄动力，其中，地球和卫星之间的万有引力称为二体问题。引力
Figure BDA0003213595240000083
的表达式如下：
Figure BDA0003213595240000084
因此，卫星的加速度为：
Figure BDA0003213595240000085
该式中，
Figure BDA0003213595240000086
表示地球中心引力，G表示万有引力常数，ME表示地球质量，
Figure BDA0003213595240000087
表示近地轨道卫星与地心的距离矢量，r表示近地轨道卫星与地心的距离数值，m表示近地轨道卫星质量。
真实地球实际是不规则椭球体，且地球内部质量分布不均匀，因此地球非球形摄动对LEO卫星的影响不可忽视，它的量级也是除地球中心引力外最大的。由LEO卫星计算非球形地球重力场摄动，其中包括将所述非球形地球重力场摄动进行归一化处理，归一化处理后的摄动函数T表示为下式
Figure BDA0003213595240000091
该式中，，GM表示地球引力常数，r表示地心向经，φ表示地心纬度，λ表示地心经度，a表示地球赤道平均半径，
Figure BDA0003213595240000092
表示规格化的位系数，
Figure BDA0003213595240000093
表示n阶m次的规格化勒让德函数，N表示地球引力位的最大阶次。
步骤六、LEO卫星根据步骤四中所述初始轨道以及步骤五中所述地球中心引力和非球形地球重力场摄动进行轨道积分外推以确定预报轨道。
计算过程中不考虑N体摄动、大气阻力、太阳光压、经验力等复杂摄动力，在本发明一个实施例中通过考虑多种地球重力场阶数，分析采用本发明方法基于不同策略对于GRACE-C卫星及HY2A卫星的轨道预报精度的影响。
策略1：考虑10×10阶的重力场摄动，海潮、固体潮对重力场的影响不予考虑，其余N体摄动、相对论、大气阻力、太阳光压、经验力等全部不予考虑。轨道预报长度为1min、5min、10min；
策略2：考虑30×30阶的重力场摄动，海潮、固体潮对重力场的影响不予考虑，其余N体摄动、相对论、大气阻力、太阳光压、经验力等全部不予考虑。轨道预报长度为1min、5min、10min；
策略3：考虑60×60阶的重力场摄动，海潮、固体潮对重力场的影响不予考虑，其余N体摄动、相对论、大气阻力、太阳光压、经验力等全部不予考虑。轨道预报长度为1min、5min、10min；
图3a-c和图4a-c分别示出了GRACE-C卫星以及HY2A卫星基于上述三种策略下三维位置(3D)方向预报精度的详细变化。其中，图3a示出了使用策略1预报GRACE-C卫星轨道的3D方向误差，图3b示出了使用策略2预报GRACE-C卫星轨道的3D方向误差，图3c示出了使用策略3预报GRACE-C卫星轨道的3D方向误差；图4a示出了使用策略1预报HY2A卫星轨道的3D方向误差，图4b示出了使用策略2预报HY2A卫星轨道的3D方向误差，图4c示出了使用策略3预报HY2A卫星轨道的3D方向误差。另外，上述轨道预报在R、T、N以及3D方向上的平均精度统计结果如表1所示
表1LEO卫星轨道预报精度统计单位：cm
Figure BDA0003213595240000101
此外，LEO卫星在完成上述步骤一至步骤六的实时轨道预报后，可以对轨道预报结果进行精度检核，若精度检核不合格则重复步骤三至步骤六，直至精度检核合格生成并且输出预报轨道。
尽管上文描述了本发明的各实施例，但是，应该理解，它们只是作为示例来呈现的，而不作为限制。对于相关领域的技术人员显而易见的是，可以对其做出各种组合、变型和改变而不背离本发明的精神和范围。因此，此处所公开的本发明的宽度和范围不应被上述所公开的示例性实施例所限制，而应当仅根据所附权利要求书及其等同替换来定义。

Claims (10)

1.一种近地轨道卫星的星上实时轨道预报方法，其特征在于，包括下列步骤：
由近地轨道卫星接收数据，并且对所述数据进行内容检验以及格式转换；
由近地轨道卫星对所述数据进行预处理；以及
由近地轨道卫星根据预处理后的数据确定预报轨道，其中包括由近地轨道卫星执行下列动作：
根据预处理后的数据进行精密单点定位解算以确定离散轨道点；
对所述离散轨道点进行拟合以确定初始轨道；
计算地球中心引力以及非球形地球重力场摄动；以及
根据所述初始轨道、地球中心引力以及非球形地球重力场摄动进行轨道积分外推以确定预报轨道。
2.根据权利要求1所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，所述近地轨道卫星具有星载全球卫星导航系统接收机，其中所述星载全球卫星导航系统接收机接收全球卫星导航系统卫星发出的数据。
3.根据权利要求2所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，近地轨道卫星接收的数据包括：近地轨道星载相位、近地轨道伪距、全球卫星导航系统轨道数据、全球卫星导航系统钟差、地球定向参数以及卫星天线相位中心改正值。
4.根据权利要求1所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，所述数据内容检验包括数据容量检验、数据日期检验以及数据来源检验。
5.根据权利要求3所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，由近地轨道卫星对所述数据进行预处理包括：
进行电离层延迟修正以确定电离层延迟误差；
进行全球卫星导航系统卫星/近地轨道卫星的天线相位中心修正以确定全球卫星导航系统卫星/近地轨道卫星的天线相位中心偏差；以及
进行相对论修正以确定相对论延迟误差。
6.根据权利要求5所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，所述全球卫星导航系统卫星/近地轨道卫星的天线相位中心修正包括将全球卫星导航系统卫星/近地轨道卫星的天线相位中心与全球卫星导航系统卫星/近地轨道卫星的质心之间的偏差α从卫星星固坐标系转换到J2000惯性坐标系中，表示为下式：
α′＝R(t)α
R(t)＝[g1，g2，g3]T
Figure FDA0003213595230000021
g3＝g1×g2
其中，R(t)表示在观测历元t时刻卫星星固坐标系与J2000惯性坐标系之间的转换矩阵，α′表示J2000惯性坐标系下全球卫星导航系统卫星/近地轨道卫星的天线相位中心与全球卫星导航系统卫星/近地轨道卫星的质心之间的偏差，r(t)、r′(t)分别表示全球卫星导航系统卫星/近地轨道卫星在J2000惯性坐标系下的位置和速度。
7.根据权利要求5所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，所述相对论修正包括计算狭义相对论效应周期性误差并且进行修正，表示为下式：
Figure FDA0003213595230000022
其中，Δρrel表示狭义相对论效应周期性误差，
Figure FDA0003213595230000023
及
Figure FDA0003213595230000024
表示近地轨道卫星的坐标及速度，c表示光速，
Figure FDA0003213595230000025
及
Figure FDA0003213595230000026
表示全球卫星导航系统卫星的坐标及速度。
8.根据权利要求5所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，由近地轨道卫星根据预处理后的数据进行精密单点定位解算包括下列步骤：
由近地轨道卫星根据预处理后的数据构建双频消电离层方程，表示为下式：
Figure FDA0003213595230000027
Figure FDA0003213595230000028
其中，角标r、s分别表示星载全球卫星导航系统接收机r以及全球卫星导航系统卫星s，
Figure FDA0003213595230000029
伪距观测值，
Figure FDA00032135952300000210
表示星载相位观测值，
Figure FDA00032135952300000211
表示全球卫星导航系统卫星s至星载全球卫星导航系统接收机r的几何距离，c表示光速，δtr表示星载全球卫星导航系统接收机钟差，δts表示全球卫星导航系统卫星钟差，δρion表示电离层延迟误差，δρrel表示相对论延迟误差，δρpco表示星载全球卫星导航系统接收机天线相位中心偏差，
Figure FDA0003213595230000031
表示全球卫星导航系统卫星天线相位中心偏差，λ表示波长，N表示相位整周模糊度，ε表示观测噪声；
将所述双频消电离层方程在近似坐标(x0，y0，z0)处展开，表示为下式：
Figure FDA0003213595230000032
其中，p表示近似展开项，xsat、ysat、zsat表示GNSS卫星坐标，ρ0表示GNSS卫星与LEO卫星的近似几何距离；
基于伪距观测值对同一历元下所有观测到的GNSS卫星构建伪距定位方程，表示为下式：
Figure FDA0003213595230000033
基于星载相位观测值，并且添加模糊度参数以构建相位定位方程；以及
通过卡尔曼滤波参数估计方法解算所述伪距定位方程和相位定位方程以确定离散轨道点。
9.根据权利要求1所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，由近地轨道卫星计算地球中心引力表示为下式：
Figure FDA0003213595230000034
其中，
Figure FDA0003213595230000035
表示地球中心引力，G表示万有引力常数，ME表示地球质量，
Figure FDA0003213595230000036
表示近地轨道卫星与地心的距离矢量，r表示近地轨道卫星与地心的距离数值，m表示近地轨道卫星质量；以及
由近地轨道卫星计算非球形地球重力场摄动，其中包括将所述非球形地球重力场摄动进行归一化处理，归一化处理后的摄动函数T表示为下式：
Figure FDA0003213595230000037
其中，GM表示地球引力常数，r表示地心向经，φ表示地心纬度，λ表示地心经度，a表示地球赤道平均半径，
Figure FDA0003213595230000038
表示规格化的位系数，
Figure FDA0003213595230000039
表示n阶m次的规格化勒让德函数，N表示地球引力位的最大阶次。
10.根据权利要求1所述的近地轨道卫星的星上实时轨道预报方法，其特征在于，还包括对所述预报轨道进行精度检核，其中当所述预报轨道符合精度要求时由近地轨道卫星输出所述预报轨道，当所述预报轨道不符合精度要求时由近地轨道卫星重复根据预处理后的数据确定预报轨道的步骤。
CN202110941066.0A 2021-08-16 2021-08-16 一种近地轨道卫星的星上实时轨道预报方法 Active CN113624243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110941066.0A CN113624243B (zh) 2021-08-16 2021-08-16 一种近地轨道卫星的星上实时轨道预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110941066.0A CN113624243B (zh) 2021-08-16 2021-08-16 一种近地轨道卫星的星上实时轨道预报方法

Publications (2)

Publication Number Publication Date
CN113624243A CN113624243A (zh) 2021-11-09
CN113624243B true CN113624243B (zh) 2022-08-26

Family

ID=78385966

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110941066.0A Active CN113624243B (zh) 2021-08-16 2021-08-16 一种近地轨道卫星的星上实时轨道预报方法

Country Status (1)

Country Link
CN (1) CN113624243B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114383619B (zh) * 2021-12-07 2023-09-05 上海航天控制技术研究所 一种高精度轨道计算方法
CN116208221B (zh) * 2022-09-07 2023-11-21 北京航天驭星科技有限公司 一种超低轨卫星地面站数据传输跟踪方法及相关设备
CN117194869B (zh) * 2023-11-07 2024-03-19 中国科学院国家授时中心 顾及姿态的低轨卫星天线相位中心预报及拟合方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101435863A (zh) * 2008-12-25 2009-05-20 武汉大学 一种导航卫星实时精密定轨的方法
CN109946727A (zh) * 2019-03-21 2019-06-28 武汉大学 一种低轨导航星座增强的网络rtk方法
CN110837094A (zh) * 2019-11-21 2020-02-25 中国人民解放军军事科学院国防科技创新研究院 基于低轨卫星的无奇点20轨道根数拟合方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100809425B1 (ko) * 2006-09-29 2008-03-05 한국전자통신연구원 ＧＰＳ와 Ｇａｌｉｌｅｏ 데이터를 이용한 정밀 궤도결정시스템 및 그 방법
US9618625B2 (en) * 2012-07-06 2017-04-11 Apple Inc. System and method for continuous carrier wave reconstruction
CN107229061B (zh) * 2017-07-18 2019-09-03 武汉大学 一种基于低轨卫星的星地差分实时精密定位方法
CN111505679B (zh) * 2020-04-20 2022-05-03 中国科学院国家空间科学中心 一种基于星载gnss的leo初轨确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101435863A (zh) * 2008-12-25 2009-05-20 武汉大学 一种导航卫星实时精密定轨的方法
CN109946727A (zh) * 2019-03-21 2019-06-28 武汉大学 一种低轨导航星座增强的网络rtk方法
CN110837094A (zh) * 2019-11-21 2020-02-25 中国人民解放军军事科学院国防科技创新研究院 基于低轨卫星的无奇点20轨道根数拟合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
LEO星座增强GNSS的精密单点定位初步分析;赵兴隆等;《大地测量与地球动力学》;20200630;第40卷(第6期);第624-628页 *

Also Published As

Publication number Publication date
CN113624243A (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
CN113624243B (zh) 一种近地轨道卫星的星上实时轨道预报方法
Montenbruck et al. Precision real-time navigation of LEO satellites using global positioning system measurements
CN111947667B (zh) 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法
CN108120994B (zh) 一种基于星载gnss的geo卫星实时定轨方法
Jäggi et al. Phase center modeling for LEO GPS receiver antennas and its impact on precise orbit determination
Yang et al. Chinese navigation satellite systems
CN110764127B (zh) 易于星载在轨实时处理的编队卫星相对定轨方法
Hwang et al. GPS‐Based Orbit Determination for KOMPSAT‐5 Satellite
Bock et al. GPS single-frequency orbit determination for low Earth orbiting satellites
Li et al. Improving BDS-3 precise orbit determination for medium earth orbit satellites
CN113219815B (zh) 一种基于x射线脉冲星的深空授时方法
Li et al. Precise orbit determination for the FY-3C satellite using onboard BDS and GPS observations from 2013, 2015, and 2017
Beran Single-frequency, single receiver terrestrial and spaceborne point positioning
Li et al. LEO–BDS–GPS integrated precise orbit modeling using FengYun-3D, FengYun-3C onboard and ground observations
CN113253314A (zh) 一种低轨卫星间时间同步方法及系统
Zehentner Kinematic orbit positioning applying the raw observation approach to observe time variable gravity
CN113581501B (zh) 一种适用于组网低轨卫星联合定轨的系统及方法
Lu et al. Improving LEO precise orbit determination with BDS PCV calibration
Gu et al. Real-time precise point positioning based on BDS-3 global short message communication
Iiyama et al. Precise Positioning and Timekeeping in a Lunar Orbit via Terrestrial GPS Time-Differenced Carrier-Phase Measurements
CN110988932B (zh) 一种提高星载gps接收机实时钟差解算精度的方法
Schreiner et al. Precise GPS data processing for the GPS/MET radio occultation mission at UCAR
Gong et al. Precise real-time navigation of the small TJU-1 satellite using GPS, GLONASS and BDS
CN116859420A (zh) 一种低轨卫星增强全球卫星导航系统(LeGNSS)数据仿真的方法
Mikrin et al. Circumlunar spacecraft navigation using the measurements from global navigation satellite systems glonass, gps, galileo and beidou

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