CN103091723A - 降低重力卫星质心调整误差对地球重力场精度影响的方法 - Google Patents
降低重力卫星质心调整误差对地球重力场精度影响的方法 Download PDFInfo
- Publication number
- CN103091723A CN103091723A CN201310047888XA CN201310047888A CN103091723A CN 103091723 A CN103091723 A CN 103091723A CN 201310047888X A CN201310047888X A CN 201310047888XA CN 201310047888 A CN201310047888 A CN 201310047888A CN 103091723 A CN103091723 A CN 103091723A
- Authority
- CN
- China
- Prior art keywords
- partiald
- centerdot
- double star
- expression
- star
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种降低重力卫星质心调整误差对地球重力场精度影响的方法;基于新型扰动能量原理,通过建立的重力卫星质心调整观测方程,有效降低重力卫星质心调整误差对地球引力位系数精度影响,进而高精度和高空间分辨率反演地球重力场;该方法卫星重力反演精度高,质心调整误差对重力场精度影响小,易于卫星重力***需求分析,卫星观测方程物理含义明确,计算机性能要求低。因此,扰动能量法是建立高精度和高空间分辨率地球重力场模型的有效方法。
Description
一、技术领域
本发明涉及卫星大地测量学、卫星重力学、航空航天等交叉技术领域,特别是涉及一种基于新型扰动能量原理,通过建立重力卫星质心调整观测方程,有效降低重力卫星质心调整误差对地球引力位系数精度影响的方法,进而高精度和高空间分辨率反演地球重力场。
二、背景技术
自1957年10月4日成功发射第一颗人造卫星Sputnik-1以来,国内外的许多学者在利用卫星技术精密探测地球重力场方面取得了辉煌成就。21世纪是人类利用SST(Satellite-to-Satellite Tracking)和SGG(Satellite Gravity Gradiometry)技术提升对“数字地球”认知能力的新纪元。地球重力场及其时变特性反映地球表层及内部物质的空间分布、运动和变化,同时决定着大地水准面的起伏和变化。因此,确定地球重力场的精细结构及其时变不仅是大地测量学、地震学、海洋学、空间科学、国防建设等的需求,同时也将为寻求资源、保护环境和预测灾害提供重要的信息资源。
GRACE由美国宇航局(NASA)和德国航天局(DLR)共同研制开发,已于2002年3月17日发射升空。如图1所示,GRACE初始轨道高度为500km,采用近圆和极地轨道设计,轨道倾角89°,轨道离心率e<0.004,在10年的飞行使命中轨道高度由500km降到300km。GRACE-A/B双星采用卫星跟踪卫星高低/低低相结合的飞行模式(SST-HL/LL),除利用高轨GPS卫星对低轨双星精密跟踪定位,同时两颗低轨卫星在同一轨道平面内前后相互跟踪(星间距离220±50km)编队飞行,并利用共轨双星轨道摄动之差高精度测量地球重力场。GRACE卫星利用冷气微推进器和磁力矩器辅助双频GPS接收机精密定轨,利用K波段高频链路高精度测量星间距离和星间速度,利用高精度SuperSTAR静电悬浮加速度计测量作用于卫星的非保守力。
在地心惯性系中研究卫星绕地球的运动规律,通常将卫星视为质点。因此,在卫星飞行中作用于卫星的非保守力可以等效为作用于卫星的质点处。在卫星重力测量中,为了将地球引力从卫星受到的合外力中有效分离,作用于卫星的非保守力的精确扣除是能否反演高精度和高空间分辨率地球重力场的关键技术,因此GRACE星载加速度计检验质量的质心要求精确定位于卫星体的质心处。由于卫星在实际飞行中,卫星体的质心和星载加速度计检验质量的质心实时存在偏移,因此质心偏差的研究是加速度计能否将作用于卫星体的非保守力精确扣除的关键技术。GRACE卫星体和星载加速度计检验质量的质心偏差源主要来自于两个方面:第一,地面安装误差源:由于在地面安装时卫星体质心和加速度计检验质量质心存在偏移,导致了GRACE卫星加速度计的静电力和作用于卫星的非保守力存在固有偏差。第二,在轨飞行误差源:由于空间环境(温度、压力等)的复杂性导致在轨飞行的卫星发生形变以及对卫星进行实时轨道和姿态控制引起喷气燃料消耗(每2~3min喷气1次,每次喷气时间200~300ms),将会导致GRACE卫星体和星载加速度计检验质量的质心存在实时偏差。由于GRACE卫星体和星载加速度计检验质量的质心偏差和卫星姿态测量具有耦合效应,因此在地球重力场反演时会同时将卫星姿态测量误差引入卫星观测方程。GRACE星体和加速度计检验质量的质心偏差以及卫星姿态测量误差的引入必将会在加速度计的三轴测量中附加扰动误差,从而影响地球重力场反演的精度。因此,GRACE星体和星载加速度计检验质量质心偏差的***研究是提高地球重力场反演精度的重要保证。
国际卫星重力测量计划的成功实施对我国既存在机遇又不乏挑战。我国应尽快汲取国外长期积累的成功经验,积极推动我国卫星重力测量计划的实施,加快自主研制重力卫星的步伐,通过卫星重力测量计划的实现带动相关领域(地学、航天、电子、通信、材料等)的发展。目前国内大地测量学界紧跟国际卫星重力测量的动态,积极投身于重力卫星需求分析的研究当中。如果GRACE星体和星载加速度计检验质量的质心调整精度设计合理,在保证利用加速度计精确扣除非保守力的前提下,可以适当降低星载加速度计和质心调整***研制的难度以及避免不必要的人力、物力和财力的浪费。基于此目的,本发明开展了GRACE星体和星载加速度计检验质量的不同质心偏差影响GRACE地球重力场精度的研究论证,不仅为我国卫星重力测量计划中质心调整精度的优化选取提供了理论基础和技术保证,同时为将来研究工作的深入开展搭建了桥梁和纽带。
三、发明内容
本发明的目的是:基于新型扰动能量原理,通过建立重力卫星质心调整观测方程,有效降低重力卫星质心调整误差对地球引力位系数精度影响,进而高精度和高空间分辨率反演地球重力场。
本发明采用能量守恒原理反演地球重力场,并在此基础上引入参考数据,通过测量数据和参考数据的残差建立新型扰动能量观测方程,进一步提高了地球重力场的反演精度。
为达到上述目的,本发明采用了如下技术方案:
一种降低重力卫星质心调整误差对地球重力场精度影响的方法,包含下列步骤:
第一步,通过真实重力卫星采集测量数据和利用理想重力卫星获取参考数据,
所述测量数据包括,
通过K波段测距仪获取星间速度数据通过星载GPS接收机获取双星轨道位置数据(r1,r2)和双星轨道速度数据通过星载加速度计获取作用于双星的非保守力数据(f01,f02),通过星敏感器获取角速度数据Ω=(Ωx,Ωy,Ωz);
所述参考数据包括,
通过星载GPS接收机获取的双星轨道位置数据(r1,r2)和双星轨道速度数据的初始点,利用9阶Runge-Kutta线性单步法和12阶Adams-Cowell线性多步法数值模拟公式获取双星参考轨道位置数据和双星参考轨道速度数据
参考非保守力数据通过DTM2000阻力温度模型计算获得;
第二步,计算各个重力卫星由于单个卫星体和星载加速度计检验质量的质心偏移而在加速度计三轴附加的非保守力Δf
其中,L表示卫星体质心和星载加速度计检验质量质心偏移在星体坐标系OS-XSYSZS中的位置矢量,和分别表示L对时间的一阶导数和二阶导数;ω=RTΩR表示在星体坐标系中重力卫星绕星体坐标系(XS,YS,ZS)轴旋转的角速度矩阵,Ω表示在地心惯性系中重力卫星绕地心坐标系(XI,YI,ZI)轴旋转的角速度矩阵
R表示由星体坐标系到地心惯性系的转换矩阵
其中,xI,yI,zI分别表示重力卫星在地心惯性坐标系中位置矢量的3个分量,分别表示重力卫星在地心惯性坐标系中速度矢量的3个分量;K表示在星体坐标系中地球引力位的梯度
第三步,反演地球重力场
在地心惯性系中,双星扰动位差观测方程建立如下
Te12=Ek12-Ef12+Vω12-VT12-V012-E012 (2)
在地心惯性系中,双星参考扰动位差观测方程建立如下
其中,带上角标“o”的参数表示参考值;
由公式(2)-公式(3)得到双星相对扰动位差观测方程
其中, 表示双星各自的地心半径,x1(2),y1(2),z1(2)分别表示双星各自位置矢量r1(2)的三个分量,θ1和θ2分别表示双星各自的地心余纬度,λ1和λ2分别表示双星各自的地心经度;Re表示地球的平均半径,L0表示地球扰动位按球函数展开的最大阶数;表示规格化的Legendre函数,l表示阶数,m表示次数;表示待求的规格化地球引力位系数;利用参考数据同样计算参考地球扰动位
双星相对扰动位差观测方程(4)右边第一项表示双星的扰动动能差, 表示双星的动能差,表示双星的相对速度矢量,e12=r12/|r12|表示由第一颗卫星指向第二颗卫星的单位方向矢量,r12=r2-r1表示双星的相对位置矢量;利用参考数据同样计算双星的参考动能差
第三项表示双星的扰动旋转能差,表示双星的旋转能差,ωe表示地球的自转角速度,x12=x2-x1和y12=y2-y1分别表示双星的相对轨道位置,和分别表示双星的相对轨道速度;利用参考数据同样计算双星的参考旋转能差
基于最小二乘法解算双星相对扰动位差观测方程(4),进而获得地球引力位系数
本发明是基于扰动能量法有利于快速反演高精度和高空间分辨率地球重力场的特点而设计的,优点是:
1)卫星重力反演精度高;
2)质心调整误差对重力场精度影响小;
3)易于卫星重力***需求分析;
4)卫星观测方程物理含义明确;
5)计算机性能要求低。
四、附图说明
图1表示GRACE重力卫星。
图2表示GRACE星体和加速度计检验质量的质心偏移。
图3表示基于扰动能量法利用不同质心调整精度反演地球引力位系数精度。
图4表示基于未引入参考数据和引入参考数据反演地球重力场精度的对比。
五、具体实施方式
以下结合附图,对本发明作进一步的详细说明。
降低重力卫星质心调整误差对地球重力场精度影响的方法包含下列步骤:
步骤一:重力卫星测量数据采集
(1)通过K波段测距仪获取星间速度数据通过星载GPS接收机获取双星轨道位置数据(r1,r2)和双星轨道速度数据通过星载加速度计获取作用于双星的非保守力数据(f1,f2),通过星敏感器获取角速度数据Ω=(Ωx,Ωy,Ωz)。
(2)通过星载GPS接收机获取的双星轨道位置数据(r1,r2)和双星轨道速度数据的初始点,利用9阶Runge-Kutta线性单步法和12阶Adams-Cowell线性多步法数值模拟公式获取双星参考轨道位置数据和双星参考轨道速度数据
(5)分别将测量数据和参考数据代入国际公布模型DE-405、IERS96和CSR4.0联合计算获取作用于双星的三体摄动能差数据VT12和参考三体摄动能差数据
步骤二:重力卫星质心调整观测方程建立
如图2所示,在地心惯性系OI-XIYIZI中,GRACE卫星体质心和星载加速度计检验质量质心的位置矢量关系如下
rA=rS+RL (5)其中,rA和rS分别表示星载加速度计质心和卫星体质心在地心惯性系中的位置矢量,L表示卫星体质心和星载加速度计检验质量质心偏移在星体坐标系OS-XSYSZS中的位置矢量,和分别表示L对时间的一阶导数和二阶导数,R表示由星体坐标系到地心惯性系的转换矩阵
在公式(5)两边同时对时间t求导,可得速度运动方程
其中,Ω=(Ωx,Ωy,Ωz)表示在地心惯性系中GRACE绕(XI,YI,ZI)轴旋转的角速度矩阵。
在公式(6)两边同时对时间t求导,可得加速度运动方程
在公式(7)两边同时左乘RT,可得星体坐标系中的加速度运动方程
ω2=RTΩ2R表示在星体坐标系中作用于GRACE的离心角加速度矩阵
在星体坐标系中,将公式(8)由矩阵形式改写为矢量形式
在星体坐标系中,GRACE卫星的动力学方程为
在星体坐标系中,GRACE星载加速度计检验质量的动力学方程为
其中,和fA分别表示作用于GRACE星载加速度计检验质量质心的引力位梯度和静电力(等效为非保守力测量值)。
由公式(11)-公式(10)得
其中,Δf=fA-fS表示由于GRACE卫星体和星载加速度计检验质量的质心偏移引起的加速度计非保守力测量值和真值的偏差;在卫星质心OS处将星载加速度计检验质量质心的引力位梯度按泰勒展开(取零阶和一阶项)
由公式(9)和公式(12)联合可得,由于GRACE卫星体和星载加速度计检验质量的质心偏移而在加速度计三轴附加的非保守力
步骤三:地球重力场反演
本发明基于能量守恒原理反演地球重力场,同时引入参考数据,通过测量数据和参考数据的残差建立新型扰动能量观测方程,进一步提高了地球重力场的反演精度。
在地心惯性系中,双星扰动位差观测方程建立如下
Te12=Ek12-Ef12+Vω12-VT12-V012-E012 (16)
在地心惯性系中,双星参考扰动位差观测方程建立如下
其中,带上角标“o”的参数表示参考值。
由公式(16)-公式(17)可得
δTe12=δEk12-δEf12+δVω12-δVT12-δV012-δE012 (18)其中,表示双星相对扰动位差,Te12(r,θ,λ)表示地球扰动位
其中, 表示双星各自的地心半径,x1(2),y1(2),z1(2)分别表示双星各自位置矢量r1(2)的三个分量,θ1和θ2分别表示双星各自的地心余纬度,λ1和λ2分别表示双星各自的地心经度;Re表示地球的平均半径,L0表示地球扰动位按球函数展开的最大阶数;表示规格化的Legendre函数,l表示阶数,m表示次数;表示待求的规格化地球引力位系数;利用参考数据同样计算参考地球扰动位
双星相对扰动位差观测方程(18)右边第一项表示双星的扰动动能差, 表示双星的动能差,表示双星的相对速度矢量,e12=r12/|r12|表示由第一颗卫星指向第二颗卫星的单位方向矢量,r12=r2-r1表示双星的相对位置矢量;利用参考数据同样计算双星的参考动能差第二项表示双星的扰动耗散能差,表示双星的耗散能差;利用参考数据同样计算双星的参考耗散能差第三项表示双星的扰动旋转能差,表示双星的旋转能差,ωe表示地球的自转角速度,x12=x2-x1和y12=y2-y1分别表示双星的相对轨道位置,和分别表示双星的相对轨道速度;利用参考数据同样计算双星的参考旋转能差第四项表示双星的扰动三体摄动能差,VT12表示双星的三体摄动能差,表示双星的参考三体摄动能差;第五项表示双星的扰动中心引力位差,V012=GM/r2-GM/r1表示双星的中心引力位差;利用参考轨道数据同样计算双星的参考中心引力位差 最后一项 表示双星的扰动能量积分常数差,E012表示双星的能量积分常数差,通过卫星的初始位置和初始速度计算得到;表示双星的参考能量积分常数差,通过卫星的参考初始位置和参考初始速度计算得到。
GRACE星体和加速度计的质心调整精度影响地球重力场精度的具体原理和计算步骤如下:计算作用于卫星体的非保守力矢量和为和将位于星体系的和转化到地心系计算作用于双星的非保守力和将位于星体系的和转化到地心系计算作用于双星的参考非保守力和并代入公式(18)的双星耗散能差Ef12,然后基于最小二乘法解算公式(18),进而获得地球引力位系数在重力卫星飞行中,由于需要利用质心调节装置实时补偿星体和加速度计检验质量的质心偏移,因此会在加速度计三轴测量中引入新的误差源—质心偏移误差,进而影响地球重力场反演的精度。
如图3所示,实线表示德国波兹坦地学研究中心(GFZ)公布的120阶EIGEN-GRACE02S地球重力场模型的精度;圆圈线、三角线、虚线和十字线分别表示基于新型扰动能量法,当质心调整精度设计为0m、5×10-5m、1×10-4m和5×10-4m时,反演地球引力位系数精度。结果表明:第一,在120阶处,当质心调整精度设计为0m时,反演地球引力位系数精度为4.911×10-10;当质心调整精度分别设计为5×10-5m、1×10-4m和5×10-4m时,反演精度依次降低至6.551×10-10、1.171×10-9和4.760×10-9。第二,以EIGEN-GRACE02S模型的地球引力位系数精度为标准,当质心调整精度设计为(5~10)×10-5m时,其与K波段测距仪、GPS接收机、SuperSTAR加速度计等GRACE关键载荷的精度指标相匹配,对地球重力场反演精度的影响较小。
图4表示基于未引入参考数据(公式(16))和引入参考数据(公式(18))分别反演GRACE地球重力场精度对比。虚线表示基于未引入参考数据反演120阶GRACE地球重力场精度;在120阶处,累计大地水准面精度为1.762×10-1m。实线表示基于引入参考数据反演120阶GRACE地球重力场精度;在120阶处,累计大地水准面精度为6.113×10-2m。研究结果表明:基于引入参考数据反演120阶GRACE累计大地水准面精度较基于未引入参考数据的反演精度平均提高2~3倍。因此,引入参考数据可以大幅提高地球重力场反演精度,有利于建立下一代高精度和高空间分辨率地球重力场模型,优化选取卫星重力测量计划中质心调整精度,以及降低星载加速度计和质心调整***研制的难度。
Claims (1)
1.一种降低重力卫星质心调整误差对地球重力场精度影响的方法,包含下列步骤:
第一步,采集重力卫星测量数据和获取参考数据,
所述测量数据包括,
通过K波段测距仪获取星间速度数据通过星载GPS接收机获取双星轨道位置数据(r1,r2)和双星轨道速度数据通过星载加速度计获取作用于双星的非保守力数据(f01,f02),通过星敏感器获取角速度数据Ω=(Ωx,Ωy,Ωz);
所述参考数据包括,
通过星载GPS接收机获取的双星轨道位置数据(r1,r2)和双星轨道速度数据的初始点,利用9阶Runge-Kutta线性单步法和12阶Adams-Cowell线性多步法数值模拟公式获取双星参考轨道位置数据和双星参考轨道速度数据
第二步,计算各个重力卫星由于单个卫星体和星载加速度计检验质量的质心偏移而在加速度计三轴附加的非保守力Δf
R表示由星体坐标系到地心惯性系的转换矩阵
第三步,反演地球重力场
在地心惯性系中,双星扰动位差观测方程建立如下
Te12=Ek12-Ef12+Vω12-VT12-V012-E012 (2)
在地心惯性系中,双星参考扰动位差观测方程建立如下
由公式(2)-公式(3)得到双星相对扰动位差观测方程
其中, 表示双星各自的地心半径,x1(2),y1(2),z1(2)分别表示双星各自位置矢量r1(2)的三个分量,θ1和θ2分别表示双星各自的地心余纬度,λ1和λ2分别表示双星各自的地心经度;Re表示地球的平均半径,L0表示地球扰动位按球函数展开的最大阶数;表示规格化的Legendre函数,l表示阶数,m表示次数;表示待求的规格化地球引力位系数;利用参考数据同样计算参考地球扰动位
双星相对扰动位差观测方程(4)右边第一项表示双星的扰动动能差, 表示双星的动能差,表示双星的相对速度矢量,e12=r12/|r12|表示由第一颗卫星指向第二颗卫星的单位方向矢量,r12=r2-r1表示双星的相对位置矢量;利用参考数据同样计算双星的参考动能差
第三项表示双星的扰动旋转能差,表示双星的旋转能差,ωe表示地球的自转角速度,x12=x2-x1和y12=y2-y1分别表示双星的相对轨道位置,和分别表示双星的相对轨道速度;利用参考数据同样计算双星的参考旋转能差
第四项表示双星的扰动三体摄动能差,VT12表示双星的三体摄动能差,表示双星的参考三体摄动能差;
基于最小二乘法解算双星相对扰动位差观测方程(4),进而获得地球引力位系数
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310047888.XA CN103091723B (zh) | 2013-02-06 | 2013-02-06 | 降低重力卫星质心调整误差对地球重力场精度影响的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310047888.XA CN103091723B (zh) | 2013-02-06 | 2013-02-06 | 降低重力卫星质心调整误差对地球重力场精度影响的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103091723A true CN103091723A (zh) | 2013-05-08 |
CN103091723B CN103091723B (zh) | 2015-05-13 |
Family
ID=48204523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310047888.XA Expired - Fee Related CN103091723B (zh) | 2013-02-06 | 2013-02-06 | 降低重力卫星质心调整误差对地球重力场精度影响的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103091723B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901496A (zh) * | 2014-03-26 | 2014-07-02 | 哈尔滨工程大学 | 一种基于光纤陀螺sins与北斗的重力测量方法 |
CN106052694A (zh) * | 2016-07-11 | 2016-10-26 | 中南大学 | 基于重力矢量及其梯度张量对单个运动物体进行定位跟踪的方法 |
CN111483619A (zh) * | 2020-04-20 | 2020-08-04 | 中国科学院微小卫星创新研究院 | 航天器引力加速度计算方法及轨道控制方法 |
CN111551973A (zh) * | 2020-04-16 | 2020-08-18 | 北京踏歌智行科技有限公司 | 一种露天矿无人驾驶惯导***的故障检测和矫正方法 |
CN113805493A (zh) * | 2021-09-01 | 2021-12-17 | 哈尔滨工业大学 | 空间双星高精度跟踪指向演练装置与方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909381A (en) * | 1997-02-19 | 1999-06-01 | Itt Manufacturing Enterprises, Inc. | System of on board prediction of trajectories for autonomous navigation of GPS satellites |
US20050071118A1 (en) * | 2002-08-28 | 2005-03-31 | Hiroshi Usuda | Electronic apparatus, signal compensation device, and signal compensation method |
CN101793976A (zh) * | 2010-02-24 | 2010-08-04 | 中国测绘科学研究院 | 一种地球重力场数据的四维动态可视化分析方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102323450A (zh) * | 2011-05-19 | 2012-01-18 | 中国科学院测量与地球物理研究所 | 基于双星相邻能量差分原理的星载加速度计数据标校方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
-
2013
- 2013-02-06 CN CN201310047888.XA patent/CN103091723B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909381A (en) * | 1997-02-19 | 1999-06-01 | Itt Manufacturing Enterprises, Inc. | System of on board prediction of trajectories for autonomous navigation of GPS satellites |
US20050071118A1 (en) * | 2002-08-28 | 2005-03-31 | Hiroshi Usuda | Electronic apparatus, signal compensation device, and signal compensation method |
CN101793976A (zh) * | 2010-02-24 | 2010-08-04 | 中国测绘科学研究院 | 一种地球重力场数据的四维动态可视化分析方法 |
CN102323450A (zh) * | 2011-05-19 | 2012-01-18 | 中国科学院测量与地球物理研究所 | 基于双星相邻能量差分原理的星载加速度计数据标校方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
Non-Patent Citations (2)
Title |
---|
吴星等: "基于卫星重力梯度数据的线质量调和分析", 《大地测量与地球动力学》 * |
郑伟 等: "地球重力场模型研究进展和现状", 《大地测量与地球动力学》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901496A (zh) * | 2014-03-26 | 2014-07-02 | 哈尔滨工程大学 | 一种基于光纤陀螺sins与北斗的重力测量方法 |
CN106052694A (zh) * | 2016-07-11 | 2016-10-26 | 中南大学 | 基于重力矢量及其梯度张量对单个运动物体进行定位跟踪的方法 |
CN111551973A (zh) * | 2020-04-16 | 2020-08-18 | 北京踏歌智行科技有限公司 | 一种露天矿无人驾驶惯导***的故障检测和矫正方法 |
CN111551973B (zh) * | 2020-04-16 | 2022-04-05 | 北京踏歌智行科技有限公司 | 一种露天矿无人驾驶惯导***的故障检测和矫正方法 |
CN111483619A (zh) * | 2020-04-20 | 2020-08-04 | 中国科学院微小卫星创新研究院 | 航天器引力加速度计算方法及轨道控制方法 |
CN113805493A (zh) * | 2021-09-01 | 2021-12-17 | 哈尔滨工业大学 | 空间双星高精度跟踪指向演练装置与方法 |
CN113805493B (zh) * | 2021-09-01 | 2022-10-21 | 哈尔滨工业大学 | 空间双星高精度跟踪指向演练装置与方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103091723B (zh) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100585602C (zh) | 惯性测量***误差模型验证试验方法 | |
CN102621565B (zh) | 一种机载分布式pos的传递对准方法 | |
CN102879014B (zh) | 深空探测接近过程的光学成像自主导航半物理仿真试验*** | |
CN102305949B (zh) | 利用星间距离插值建立全球重力场模型的方法 | |
CN103018783B (zh) | 重力卫星编队轨道稳定性优化设计和精密反演地球重力场方法 | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN102393535B (zh) | 基于双星能量插值原理的卫星重力反演方法 | |
CN102736118B (zh) | 一种用于全球重力场测量的综合型卫星*** | |
CN103674030A (zh) | 基于天文姿态基准保持的垂线偏差动态测量装置和方法 | |
CN103076640A (zh) | 利用方差-协方差对角张量原理反演地球重力场的方法 | |
CN102313905B (zh) | 基于星间速度插值原理的地球重力反演方法 | |
CN102262248A (zh) | 基于双星空间三维插值原理的卫星重力反演方法 | |
CN103091723A (zh) | 降低重力卫星质心调整误差对地球重力场精度影响的方法 | |
CN102788598B (zh) | 基于三轴旋转的光纤捷联惯导***误差抑制方法 | |
Zheng et al. | Efficient and rapid estimation of the accuracy of future GRACE Follow‐On Earth's gravitational field using the analytic method | |
CN103017787A (zh) | 适用于摇摆晃动基座的初始对准方法 | |
CN104296908A (zh) | 三自由度气浮台干扰力矩组成测量装置 | |
CN103076639B (zh) | 基于残余星间速度原理反演地球重力场的方法 | |
CN102998713B (zh) | 基于功率谱半解析的卫星重力梯度反演方法 | |
CN103645489A (zh) | 一种航天器gnss单天线定姿方法 | |
CN108020866A (zh) | 一种星体重力场反演的方法和***、以及处理器 | |
CN104568652A (zh) | 一种高精度测量临近空间大气密度的方法及测量装置 | |
CN106997061A (zh) | 一种基于扰动星间相对速度提高重力场反演精度的方法 | |
CN102519455A (zh) | 基于紫外敏感器的自主导航半物理仿真试验*** | |
CN102323450B (zh) | 基于双星相邻能量差分原理的星载加速度计数据标校方法 |
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 |
Granted publication date: 20150513 Termination date: 20160206 |
|
CF01 | Termination of patent right due to non-payment of annual fee |