CN112415559A - 一种基于多项式展开技术的高阶容错卫星轨道确定方法 - Google Patents
一种基于多项式展开技术的高阶容错卫星轨道确定方法 Download PDFInfo
- Publication number
- CN112415559A CN112415559A CN202011164514.2A CN202011164514A CN112415559A CN 112415559 A CN112415559 A CN 112415559A CN 202011164514 A CN202011164514 A CN 202011164514A CN 112415559 A CN112415559 A CN 112415559A
- Authority
- CN
- China
- Prior art keywords
- state
- measurement
- satellite orbit
- value
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005516 engineering process Methods 0.000 title claims abstract description 12
- 238000005259 measurement Methods 0.000 claims abstract description 73
- 239000011159 matrix material Substances 0.000 claims abstract description 55
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 18
- 239000000126 substance Substances 0.000 claims description 3
- 238000012938 design process Methods 0.000 claims description 2
- 239000011541 reaction mixture Substances 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 3
- 238000013459 approach Methods 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000546 chi-square test Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
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
- G01S19/50—Determining position whereby the position solution is constrained to lie upon a particular curve or surface, e.g. for locomotives on railway tracks
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于多项式展开技术的高阶容错卫星轨道确定方法,涉及卫星的轨道确定领域。通过动态地调整估计过程中的测量误差协方差矩阵,确保其与实际的测量误差相匹配,从而提高非线性高阶容错扩展卡尔曼滤波器的滤波一致性和对错误测量的鲁棒性,提高卫星轨道确定的精度。本发明方法综合考虑了动力学方程和测量方程的非线性特性,在滤波器预测阶段,采用高阶泰勒多项式逼近高精度预测卫星状态,在保证估计精度的要求下,有助于减少滤波器对测量新息的更新需求和降低滤波器对测量频率的要求;容错策略的设计进一步保证了滤波器的鲁棒性。能够普遍用于卫星、空间碎片等的轨道确定,可缓解测量设备的测量负担并保持对错误测量的鲁棒性。
Description
技术领域
本发明属于卫星轨道确定领域,涉及一种基于多项式展开技术的高阶容错卫星轨道确定方法。
背景技术
近年来,随着航天活动的精细化,要求在轨卫星轨道确定的精度不断增加,致使传统轨道确定方法需要计算大量的测量数据,加剧了与地面有限的观测系统之间的矛盾。相比于卡尔曼滤波器、扩展卡尔曼滤波器、无味卡尔曼滤波器和粒子滤波器等传统的滤波器,基于泰勒多项式展开技术的非线性高阶扩展卡尔曼滤波器很好地平衡了滤波器的精度和计算效率。同时,由于非线性高阶扩展卡尔曼滤波器采用高阶泰勒多项式逼近非线性的卫星轨道动力学方程,卫星状态的预测精度得到了明显的提升,因此可以降低对测量频率的要求。这个优点很好的缓和了大量在轨卫星的测量需求和有限测控站之间的矛盾。
一般而言,航天器轨道的观测数据由于受到观测过程误差甚至是偶然错误的影响,往往会引起滤波器的估计错误。类似于传统的扩展卡尔曼滤波器、无味卡尔曼滤波器和粒子滤波器等,非线性高阶扩展卡尔曼滤波器也面临着相同的问题。测量过程中,随机测量错误或者测量设备的固有偏差,会导致滤波器无法对卫星在某一时刻的状态进行准确估计。如何克服测量错误对非线性高阶扩展卡尔曼滤波器滤波过程的影响,具有非常重要的实际意义。
发明内容
本发明的目的在于克服上述现有技术中,非线性高阶扩展卡尔曼滤波器因存在随机测量错误而无法准确估计卫星在某一时刻的状态的缺点,提供一种基于多项式展开技术的高阶容错卫星轨道确定方法。
为了达到上述目的,本发明采用以下技术方案予以实现:
一种基于多项式展开技术的高阶容错卫星轨道确定方法,包括如下步骤:
步骤一,建立卫星轨道动力学方程;
步骤二,建立航天器的观测方程;
步骤三,基于泰勒展开多项式分别对步骤一的卫星轨道动力学方程和步骤二的观测方程进行非线性预测,得到卫星轨道的状态预测值和航天器的观测预测值,并结合真实观测值计算预测误差,得到协方差矩阵;
步骤四,根据步骤三得到的测量误差,设计缩放矩阵因子,自适应调整测量协方差矩阵,保证滤波一致性;
步骤五,设置观测时间,对步骤三的卫星轨道的状态预测值进行更新,得到任意时刻航天器的轨道估计状态。
优选地,步骤一所述的卫星轨道动力学方程为:
式(1)中,x表示n维状态矢量,t表示时间。
优选地,步骤二所述的观测方程为:
z=h(x,t)+u (2)
其中,z表示t时刻的观测量,x表示t时刻的状态值,u表示t时刻的观测噪声。
优选地,步骤三的具体操作为:
得到卫星轨道的预测状态为:
优选地,步骤三所述的航天器的预测观测值为:
zk+1=h(xk+1,tk+1)+uk+1 (6)
其中,zk+1表示tk+1时刻的观测量,xk+1表示tk+1时刻的状态预测值,uk+1表示 tk+1时刻的观测噪声;
式(8)中,E[]表示期望值。
优选地,步骤四所述的更新具体为:
其中,为卫星轨道预测状态的协方差矩阵;γ=[γ1,…,γn]T和l=[l1,…,ln]T,表示变量偏差分量的指数,p和q表示测量矢量的索引值,i表示状态的索引值,表示观测噪声的协方差矩阵,系数与表示泰勒展开项的系数。
优选地,步骤四的缩放矩阵的设计过程具体为:
定义测量新息矢量为
与现有技术相比,本发明具有以下有益效果:
本发明公开了一种基于多项式展开技术的高阶容错卫星轨道确定方法,通过动态地监测测量新息,动态调整协方差矩阵用以匹配实际的测量误差特性,提出一种自适应容错策略。在滤波器使用高阶非线性泰勒展开式对状态预测过程结束后,通过计算滤波过程中的测量新息并检验滤波一致性特性,动态地设计一个缩放矩阵用以调节测量误差协方差矩阵,保证其匹配真实的测量误差和测量错误。提出的滤波器充分利用了高阶泰勒展开多项式能够精确逼近非线性过程的特点,减少滤波器在预测过程和更新过程中的精度丢失,进而提高精度。同时,高阶泰勒展开是一种半解析方法,相比于传统的非线性滤波方法可以加快运算速度。本发明提出的容错设计方法改善了传统的非线性高阶扩展卡尔曼滤波器对错误测量的鲁棒性,不仅保持原有非线性高阶扩展卡尔曼滤波特性,即高精度、高效率,适用于强非线性系统和低测量频率情况,而且可以规避由于测量错误导致的滤波发散问题。综合考虑动力学方程和测量方程的非线性特性,在滤波器预测阶段,采用高阶泰勒多项式逼近高精度预测卫星状态,在保证估计精度的要求下,有助于减少滤波器对测量新息的更新需求和降低滤波器对测量频率的要求;容错策略的设计进一步保证了滤波器的鲁棒性。本发明提出的基于多项式展开技术的高阶容错卫星轨道确定方法普遍用于卫星、空间碎片等的轨道确定,可缓解测量设备的测量负担并保持对错误测量的鲁棒性。
进一步地,缩放矩阵Gk+1本身是一个对角矩阵,且由于测量元器件在任务执行过程中一般不可以进一步提升性能,各元素均大于1,因此,最终的缩放矩阵设计为I表示m维单位矩阵。缩放矩阵可以进一步用于计算卡尔曼滤波器的增益。
附图说明
图1为标准的高阶非线性扩展卡尔曼滤波器和本发明提出的高阶非线性容错扩展卡尔曼滤波器得到的位置误差示意图;
图2为标准的高阶非线性扩展卡尔曼滤波器和本发明提出的高阶非线性容错扩展卡尔曼滤波器得到的速度误差示意图;
图3为标准的高阶非线性扩展卡尔曼滤波器和本发明提出的高阶非线性容错扩展卡尔曼滤波器夫人新息平方归一化值示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述:
为了更好地说明本发明的目的和优点,下面结合附图和实例对本发明内容做进一步说明。
实施例1:
本实例公开一种基于多项式展开技术的高阶容错卫星轨道确定方法,具体实现步骤如下:
步骤一:建立无量纲的卫星轨道状态方程为
其中分别表示无量纲的时间,卫星的无量纲位置矢量和角速度值,t,r,ω表示相应的有量纲值,表示卫星无量纲位置矢量的范数;另外,卫星的标称角速度表示卫星在半长轴为a=8788km的轨道上运行的平均的轨道角速度,μe表示地球的引力常数。
步骤二:建立非线性测量方程,包括径向距离,赤经和赤纬,
其中u=[u1,u2,u3]T表示测量噪声矢量,径向距离的测量噪声满足零均值标准差为1米的高斯分布,赤经和赤纬的测量噪声满足零均值标准差为1.745×10-6 rad的高斯分布。同时,我们假设每个轨道可以获得12个等间距的测量值。
步骤三:基于泰勒多项式展开的非线性状态预测:假设tk时刻航天器的标称状态为初始偏差为δxk,状态多项式变量为微分方程在tk+1时刻的解为对该解在标称状态处做泰勒展开,得到高阶展开式其中表示状态t1状态邻域[xk+1]对初始邻域大小δxk的非线性依赖关系。
航天器的预测观测值为:
zk+1=h(xk+1,tk+1)+uk+1 (6)
其中,zk+1表示tk+1时刻的观测量,xk+1表示tk+1时刻的状态预测值,uk+1表示 tk+1时刻的观测噪声;
式(8)中,E[]表示期望值。
步骤四:非线性状态更新
其中,为卫星轨道预测状态的协方差矩阵;γ=[γ1,…,γn]T和l=[l1,…,ln]T,表示变量偏差分量的指数,p和q表示测量矢量的索引值,i表示状态的索引值,表示观测噪声的协方差矩阵,系数与表示泰勒展开项的系数。
进一步的,定义测量新息矢量为
通过设计一个缩放矩阵Gk+1,来动态的调整估计得到的测量协方差矩阵进而匹配真实的测量误差。具体的过程如下,首先定义
由于缩放矩阵Gk+1的作用是作用于测量误差协方差矩阵的对角元素,因此,缩放矩阵Gk+1本身是一个对角矩阵,且由于测量元器件在任务执行过程中一般不可以进一步提升性能,各元素均大于1,因此,最终的缩放矩阵 I表示m维单位矩阵。缩放矩阵进一步用于计算卡尔曼滤波器的增益,
步骤五:反复执行上述算法,得到任意时刻航天器的轨道状态。
本发明方法的实例验证:
运行在一个低地球轨道上的卫星的初始状态为
初始状态误差为真实状态的0.1%,同时假设在时间区间[136.6,182.2]min上,赤经赤纬的测量值上存在一个常值的测量偏差0.00001rad,其他时间测量只存在测量误差而没有测量偏差。
图1和图2展示了分别由标准的高阶非线性扩展卡尔曼滤波器和本发明提出的高阶非线性容错扩展卡尔曼滤波器得到的位置和速度误差。结果表明高阶非线性容错扩展卡尔曼滤波器具有良好的估计精度并对错误测量具有强健的鲁棒性,相反,标准的高阶非线性扩展卡尔曼滤波器无法处理错误测量并表现出一定程度的精度损失。
图3展示了新息平方归一化值(NIS)并与卡方检验临界值进行了比较,结果表明,标准的高阶非线性扩展卡尔曼滤波器在[136.6,182.2]min上的NIS值超出了临界值,表明在这一段时间内,由于错误测量的影响,滤波器的一致性遭到了破坏,导致了标准的高阶非线性扩展卡尔曼滤波器精度损失。
综上所述,本发明提出的基于多项式展开技术的高阶容错卫星轨道确定方法,通过动态地调整估计过程中的测量误差协方差矩阵,确保其与实际的测量误差相匹配,从而提高非线性高阶容错扩展卡尔曼滤波器的滤波一致性和对错误测量的鲁棒性,提高卫星轨道确定的精度。本发明方法不仅保持原有非线性高阶扩展卡尔曼滤波特性,即高精度、高效率,适用于强非线性系统和低测量频率情况,而且可以规避由于测量错误导致的滤波发散问题。该方法不仅适用于可操作的自主人造航天器,同样适用于空间碎片的轨道确定。本发明所提出的滤波方法不仅适用于卫星定轨,同样适用于强非线性系统的状态估计问题。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (8)
1.一种基于多项式展开技术的高阶容错卫星轨道确定方法,其特征在于,包括如下步骤:
步骤一,建立卫星轨道动力学方程;
步骤二,建立航天器的观测方程;
步骤三,基于泰勒展开多项式分别对步骤一的卫星轨道动力学方程和步骤二的观测方程进行非线性预测,得到卫星轨道的状态预测值和航天器的观测预测值,并结合真实观测值计算预测误差,得到协方差矩阵;
步骤四,根据步骤三得到的测量误差,设计缩放矩阵因子,自适应调整测量协方差矩阵,保证滤波一致性;
步骤五,设置观测时间,对步骤三的卫星轨道的状态预测值进行更新,得到任意时刻航天器的轨道估计状态。
3.根据权利要求1所述的高阶容错卫星轨道确定方法,其特征在于,步骤二所述的观测方程为:
z=h(x,t)+u (2)
其中,z表示t时刻的观测量,x表示t时刻的状态值,u表示t时刻的观测噪声。
4.根据权利要求1所述的高阶容错卫星轨道确定方法,其特征在于,步骤三的具体操作为:
得到卫星轨道的预测状态为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011164514.2A CN112415559B (zh) | 2020-10-27 | 2020-10-27 | 一种基于多项式展开技术的高阶容错卫星轨道确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011164514.2A CN112415559B (zh) | 2020-10-27 | 2020-10-27 | 一种基于多项式展开技术的高阶容错卫星轨道确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112415559A true CN112415559A (zh) | 2021-02-26 |
CN112415559B CN112415559B (zh) | 2024-06-25 |
Family
ID=74841165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011164514.2A Active CN112415559B (zh) | 2020-10-27 | 2020-10-27 | 一种基于多项式展开技术的高阶容错卫星轨道确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112415559B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113269309A (zh) * | 2021-05-31 | 2021-08-17 | 长春市榣顺科技有限公司 | 适用于微弱信号提取的rbf神经网络改进ukf方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987355A (zh) * | 2006-12-22 | 2007-06-27 | 北京航空航天大学 | 一种基于自适应扩展卡尔曼滤波的地球卫星自主天文导航方法 |
US20110291881A1 (en) * | 2010-05-27 | 2011-12-01 | Seiko Epson Corporation | Position calculating method and position calculating device |
KR20120075961A (ko) * | 2010-12-29 | 2012-07-09 | 세종대학교산학협력단 | 위성 데이터 스무딩을 이용한 정지궤도위성의 정밀 궤도 결정장치 및 결정 방법 |
JP2012154769A (ja) * | 2011-01-26 | 2012-08-16 | Seiko Epson Corp | 加速度検出方法、位置算出方法及び加速度検出装置 |
EP2813864A2 (en) * | 2013-05-24 | 2014-12-17 | O2 Micro, Inc. | Receivers and methods for multi-mode navigation |
CN109032176A (zh) * | 2018-07-25 | 2018-12-18 | 西北工业大学 | 一种基于微分代数的地球同步轨道确定和参数确定方法 |
CN109255096A (zh) * | 2018-07-25 | 2019-01-22 | 西北工业大学 | 一种基于微分代数的地球同步卫星轨道不确定演化方法 |
US20190251215A1 (en) * | 2018-02-15 | 2019-08-15 | Regents Of The University Of Minnesota | Accurate estimation of upper atmospheric density using satellite observations |
CN110779532A (zh) * | 2019-11-18 | 2020-02-11 | 河南工业大学 | 一种应用于近地轨道卫星的地磁导航系统及方法 |
WO2020133711A1 (zh) * | 2018-12-28 | 2020-07-02 | 上海海积信息科技股份有限公司 | 一种卫星定轨方法、装置及电子设备 |
-
2020
- 2020-10-27 CN CN202011164514.2A patent/CN112415559B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1987355A (zh) * | 2006-12-22 | 2007-06-27 | 北京航空航天大学 | 一种基于自适应扩展卡尔曼滤波的地球卫星自主天文导航方法 |
US20110291881A1 (en) * | 2010-05-27 | 2011-12-01 | Seiko Epson Corporation | Position calculating method and position calculating device |
KR20120075961A (ko) * | 2010-12-29 | 2012-07-09 | 세종대학교산학협력단 | 위성 데이터 스무딩을 이용한 정지궤도위성의 정밀 궤도 결정장치 및 결정 방법 |
JP2012154769A (ja) * | 2011-01-26 | 2012-08-16 | Seiko Epson Corp | 加速度検出方法、位置算出方法及び加速度検出装置 |
EP2813864A2 (en) * | 2013-05-24 | 2014-12-17 | O2 Micro, Inc. | Receivers and methods for multi-mode navigation |
US20190251215A1 (en) * | 2018-02-15 | 2019-08-15 | Regents Of The University Of Minnesota | Accurate estimation of upper atmospheric density using satellite observations |
CN109032176A (zh) * | 2018-07-25 | 2018-12-18 | 西北工业大学 | 一种基于微分代数的地球同步轨道确定和参数确定方法 |
CN109255096A (zh) * | 2018-07-25 | 2019-01-22 | 西北工业大学 | 一种基于微分代数的地球同步卫星轨道不确定演化方法 |
WO2020133711A1 (zh) * | 2018-12-28 | 2020-07-02 | 上海海积信息科技股份有限公司 | 一种卫星定轨方法、装置及电子设备 |
CN110779532A (zh) * | 2019-11-18 | 2020-02-11 | 河南工业大学 | 一种应用于近地轨道卫星的地磁导航系统及方法 |
Non-Patent Citations (1)
Title |
---|
汪栋 等: "HY-2卫星雷达高度计时标偏差估算", 《海洋学报》, vol. 35, no. 5, 30 September 2013 (2013-09-30) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113269309A (zh) * | 2021-05-31 | 2021-08-17 | 长春市榣顺科技有限公司 | 适用于微弱信号提取的rbf神经网络改进ukf方法 |
CN113269309B (zh) * | 2021-05-31 | 2022-09-13 | 长春市榣顺科技有限公司 | 适用于微弱信号提取的rbf神经网络改进ukf方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112415559B (zh) | 2024-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109032176B (zh) | 一种基于微分代数的地球同步轨道确定和参数确定方法 | |
CN109255096B (zh) | 一种基于微分代数的地球同步卫星轨道不确定演化方法 | |
EP0856784B1 (de) | Verfahren und Vorrichtung zur bordautonomen Bestimmung der Position eines Satelliten | |
CN106873611A (zh) | 一种多通道线性自抗扰控制器的设计方法 | |
CN104950898A (zh) | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 | |
Lee et al. | Sigma point filtering for sequential orbit estimation and prediction | |
He et al. | Iterative learning disturbance observer based attitude stabilization of flexible spacecraft subject to complex disturbances and measurement noises | |
CN111238489B (zh) | 一种低轨卫星大气阻力摄动建模和计算方法 | |
US6356815B1 (en) | Stellar attitude-control systems and methods with weighted measurement-noise covariance matrices | |
CN112415559A (zh) | 一种基于多项式展开技术的高阶容错卫星轨道确定方法 | |
CN110706265A (zh) | 一种改进srckf强跟踪滤波的机动目标跟踪方法 | |
Petrova et al. | Medium-term predictions of F10. 7 and F30 cm solar radio flux with the adaptive Kalman filter | |
Cui et al. | A unified model for transfer alignment at random misalignment angles based on second-order EKF | |
CN113919194A (zh) | 一种基于滤波误差法的非线性飞行动力学系统辨识方法 | |
CN112319859B (zh) | 一种基于自主滤波阶切换的非线性卫星轨道确定方法 | |
US6317662B1 (en) | Stable and verifiable state estimation methods and systems with spacecraft applications | |
Hough | Improved performance of recursive tracking filters using batch initialization and process noise adaptation | |
Su et al. | Variational Bayesian adaptive high‐degree cubature Huber‐based filter for vision‐aided inertial navigation on asteroid missions | |
CN113670315B (zh) | 一种基于变分迭代卡尔曼滤波的李群重尾干扰噪声动态飞行器姿态估计方法 | |
Zhao et al. | Unscented Kalman filtering for SINS attitude estimation | |
Ye et al. | Iterative noise estimation-based cubature Kalman filtering for distributed pos in aerial earth observation imaging | |
Seren et al. | Adaptive extended Kalman filtering for virtual sensing of longitudinal flight parameters | |
CN114200491A (zh) | 一种基于导航数据的应急航天器星历确定方法及系统 | |
Hajiyev et al. | Attitude filtering with uncertain process and measurement noise covariance using SVD‐aided adaptive UKF | |
Xiao et al. | A meticulous covariance adaptive Kalman filter for satellite attitude estimation |
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 |