CN104596520B - 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 - Google Patents
一种基于半长轴滤波值连续性的航天器轨道机动检测方法 Download PDFInfo
- Publication number
- CN104596520B CN104596520B CN201510027623.2A CN201510027623A CN104596520B CN 104596520 B CN104596520 B CN 104596520B CN 201510027623 A CN201510027623 A CN 201510027623A CN 104596520 B CN104596520 B CN 104596520B
- Authority
- CN
- China
- Prior art keywords
- semi
- spacecraft
- major axis
- detection
- maneuver
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 84
- 230000008859 change Effects 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims description 23
- 238000005070 sampling Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 11
- 238000013179 statistical model Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 13
- 238000009825 accumulation Methods 0.000 description 12
- 230000001133 acceleration Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000007123 defense Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开一种基于半长轴滤波值连续性的航天器轨道机动检测方法,包括:基于运动学模型的轨道半长轴估计;基于轨道半长轴连续性的脉冲机动检测;基于轨道半长轴变化率的连续推力机动检测。采用本发明的技术方案,可以实现对轨道机动的快速准确检测与辨识,提高航天器轨道估计精度和可靠性。
Description
技术领域
本发明涉及航空航天技术领域,尤其涉及一种基于半长轴滤波值连续性的航天器轨道机动检测方法。
技术背景
随着空间态势感知、攻防对抗和深空探测等领域任务的拓展,新型航天器的高机动、大范围运动特点使得对其进行实时高精度跟踪面临着更加严峻的挑战,尤其在空间攻防对抗等非合作目标跟踪中,未知的轨道机动极大的限制了航天器实时跟踪的精度,迫切需要有效的轨道机动检测和辨识技术。
非线性滤波是航天器实时轨道确定的主要方法,航天器轨道机动可以看作是其运动状态的突变,表现为一种***性的偏差。因此轨道机动检测就是在滤波过程中或滤波后,引入必要的方法来检测或辨识轨道是否产生了偏离。其检测过程可以参考一般性的运动目标机动检测方法,又有其特殊性。航天器常见机动模式有脉冲机动、连续推力机动和bangbang机动,其中bangbang机动可以看作脉冲机动与连续推力机动的随机组合,而且目前应用相对较少,因此研究的热点主要集中于前两种机动的检测。检测算法实现的关键是检测统计量的构造,通常有两类方法,一是从测量信息或其某种特征的变化规律判断,二是基于滤波信息或输入估计的假设检验问题。后者利用了目标运动先验信息,且可以实现测量数据的融合,应用较广,但算法的检测和跟踪精度均依赖于模型的选择。鉴于航天器的空间轨道运动特征,有动力学和运动学两类模型可以利用。前者的约束力较强,常用于高精度的跟踪滤波器设计,且基于该滤波器信息设计的检测算法对微小机动检测能力较强,但由于对检测参数的依赖性也较大,检测滞后和虚警之间存在矛盾;后者对测量信息变化敏感,由于轨道机动会引起测量数据某些特征的变化(测量数据是轨道在测量空间的投影),因此更利于航天器机动的检测,但其滤波结果的随机误差较大,不利于微小机动的检测与辨识。
常用机动检测和估计方法如等效噪声、输入检测与估计、模型转换等,但对于实时性和稳定性要求较高的实时轨道确定问题,上述方法对***改造较大,检测过程所需数据积累和计算量相对较多,且适用性因实际***状态噪声的不同而产生较大差异。此外,现有***多以航天器位置与速度为滤波状态,因此机动检测手段也经常直接针对位置与速度进行,受测量异常值,特别是成片野值的影响较大,可靠性难以保证。
发明内容:
本发明要解决的技术问题是,提供一种基于轨道半长轴连续性变化监测的轨道机动检测方法,实现对轨道机动的快速准确检测与辨识。
为解决上述问题,本发明采用如下的技术方案:
一种基于半长轴滤波值连续性的航天器轨道机动检测方法包括:
步骤S1、基于运动学模型的轨道半长轴估计
在每个采样时刻,航天器跟踪***利用对轨道运动学模型进行滤波,实时解算航天器位置与速度,进而获得该时刻的半长轴估计值;
步骤S2、基于轨道半长轴连续性的脉冲机动检测
利用窗口累积的轨道半长轴估值进行航天器脉冲机动检测,具体以下步骤:
步骤S2.1、利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴的统计特征;
步骤S2.2、采用逐点滑动的方式更新窗口数据,并通过检测新采样时刻的半长轴数据是否符合半长轴的统计特征征,进行脉冲机动检测;
步骤S2.3、若未发生机动,则利用新窗口数据计算半长轴的统计特征。
步骤S3、基于轨道半长轴变化率的连续推力机动检测
在脉冲机动检测基础上,利用窗口累积的轨道半长轴变化率估值进行航天器连续推力机动检测,具体包括以下步骤:
步骤S3.1、若未发生脉冲机动,则利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴变化率的统计特征;
步骤S3.2、采用逐点滑动的方式更新窗口数据,将窗口按时间分成前后两段,分别计算半长轴变化率,并通过检测两个变化率的差异是否符合半长轴变化率的统计特征,进行连续推力机动检测;
步骤S3.3、若未发生机动,则利用新窗口数据计算半长轴变化率的统计特征。
本发明提出利用轨道半长轴滤波值连续性的机动检测方法,具有实现简单,检测精度高、速度快和稳健性强特点,具体优点如下:
1、仅需在跟踪滤波器中引入简单的半长轴统计特征计算和连续性变化检测环节,即可分别实现脉冲机动和连续推力机动的检测,对现有跟踪***的改造较小,方法实现方便且适应性较强;
2、采用逐点滑动窗口方式,不存在累积数据自身引入的解算滞后,且仅使用了简单的多点差分算法,以及半长轴与半长轴变化率的统计特征计算,检测过程所需计算量很小,保证了机动检测的速度;
3、检测过程基于窗口内的轨道半长轴估值进行,充分利用了多个时刻的测量数据,可以降低测量异常的影响,提高检测精度的同时,保障机动检测的可靠性;
附图说明
图1是本发明的算法流程示意图;
图2是实施算例中航天器轨道的脉冲机动检测示意图;
图3是实施算例中航天器轨道的连续推力机动检测示意图。
具体实施方式
如图1所示,本发明公开一种基于半长轴滤波值连续性的航天器轨道机动检测方法,包括以下步骤:
步骤S1、基于运动学模型的航天器轨道半长轴估计
在每个采样时刻tk,航天器跟踪***利用对轨道运动学模型进行滤波,实时解算航天器位置与速度,进而获得该时刻的半长轴估计值ak。
从计算效率和精度出发,所述运动学模型为当前统计模型(Current StatisticModel,CS)。记tk时刻航天器惯性系位置为Xk=(xk,yk,zk)T,xk,yk,zk分别为在惯性系x,y,z轴的分量。当前统计模型认为机动加速度在各个方向是解耦的,并且在时间轴上是一阶时间相关的随机过程。记则状态转移方程
其中,分别为航天器tk时刻x方向的速度和加速度。分别为航天器tk时刻y方 向的速度和加速度,分别为航天器k时刻z方向的速度和加速度 dt为采样时间间隔,为机动加速度引入导致的状态噪声向量,其方差由机动加速度的概率密度函数确定。当前统计模型假设其满足修正的瑞利分布,三个方向相对独立,以x方向为例,tk时刻ωx的方差为
其中,表示tk-1时刻加速度估计值,与分别为加速度正向和负向最大值,一般根据目标运动规律先验知识获取。其它方向类似。获得轨道测量数据后,根据测量的非线性性,这里采用无迹滤波(the unscented filter,UKF)进行位置Xk=(xk,yk,zk)T与速度的实时估计,具体可过程可参考文献《A new extension of theKalman filter t0nonlinear systems》(Julier S J and Uhlmann J K.Proc AeroSense,Orlando,1997:54-65.)。进而由航天器位置Xk=(xk,yk,zk)T与速度的估值计算其轨道半长轴ak如下
其中,μ为地球引力常数,μ=398600.4418e9。
步骤S2、基于半长轴连续性的航天器轨道脉冲机动检测
利用窗口累积的轨道半长轴估值进行航天器脉冲机动检测,其检测过程包括如下步骤:
步骤S2.1、利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴的统计特征;
记ti时刻的半长轴估值为ai,若已连续获得了N个采样时刻的轨道估值{ai,ai+1,…,ai+N-1},N为积累窗口长度,采用序列差分检测来判断这些点是否满足初始化条件,计算公式如下:
其中,β为第一阈值,对于圆轨道可取值1.5~3,对于偏心率较大的轨道可取值3~6; 为半长轴的先验方差。
当窗口内任意相邻五点均满足公式(4),则初始积累完成,即认为这连续N个采样时刻对应的轨道比较平稳,无机动存在,可以作为后续检测的基础。
若初始积累完成,在计算窗口{ai,ai+1,…,ai+N-1}内,按照公式(5)和公式(6)计算窗口数据{ai,ai+1,…,ai+N-1}的中值和标准差中值和标准差为半长轴的统计特征,并转入Step 2.2;否则,窗口向前滑动一个采样点,即取{ai+1,ai+2,…,ai+N},重新初始积累。
其中,med(·)表示中值计算函数,std(·)表示标准差计算函数。
步骤S2.2、采用逐点滑动的方式更新窗口数据,并通过检测新采样时刻的半长轴数据是否符合半长轴的统计特征,进行脉冲机动检测;
当滤波获得窗口右侧采样时刻半长轴估值ai+N后,若下式成立,则认为该f+N时刻航天器潜在机动
其中,α为第二阈值,其反映了检测的灵敏度,通常取值为3~6。
为了减小测量异常的影响,保证检测可靠性,需要连续多个采样点超出阈值时才判断机动发生。这里选择四分之一的窗口数据,即当N/4个采样点满足公式(7)时,认为航天器在该时段内发生了脉冲机动,并将第一个潜在机动的时刻i+N做为脉冲机动起始时刻,窗口向右滑动一个采样点,转Step 2.1,重新开始窗口初始积累,并将初始积累结束的时刻作为脉冲机动的结束时刻;否则,窗口向右滑动一个采样点,转Step2.3。
步骤S2.3、若未发生机动,则利用新窗口数据计算半长轴的统计特征。
此时认为航天器并未处于脉冲机动段落,需要重新计算半长轴的统计特征,即中值和标准差,并转入Step 2.2,开始下一时刻的机动检测。为了更有效的抑制随机误差和异常值的影响,中值计算仍采用公式(5),即但标准差的计算引入如下的平滑方法:
其中,γ为渐消因子,取值越大则平滑效果越好,但会降低对轨道半长轴变化的敏感性,通常计算时取值0.8~0.95。
步骤S3、基于半长轴变化率连续性的航天器轨道连续推力机动检测
在脉冲机动检测基础上,利用窗口累积的轨道半长轴变化率估值进行航天器连续推力机动检测,其检测过程包括如下步骤:
步骤S3.1、在脉冲机动检测基础上,若未发生脉冲机动,则利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴变化率的统计特征;
设积累窗口的长度为N,采用二次多项式拟合方式计算轨道半长轴在窗口{ai,ai+1,…,ai+N-1}内的变化率。具体的,设采样时刻为T=[ti,ti+1,…,ti+N-1]T,则拟合系数为C=(DTD)-1DTL,变化率的计算方法为
其中,表示ai的变化率,
进一步,采用序列差分检测方法,判断是否满足初始化条件,差分检测公式为:
其中,β为第一阈值,对于圆轨道可取值1.5~3,对于偏心率较大的轨道可取值3~6; 为半长轴变化率的先验方差。
当窗口内任意相邻五点均满足公式(10),则初始积累完成,即认为这连续N个采样时刻对应的轨道比较平稳,无连续推力机动存在,可以作为后续检测的基础。
初始积累完成后,计算半窗口数据的标准差标准差为半长轴变化率的统计特征,其中,[·]表示取整运算,并转入Step 3.2。
标准差计算公式为:
需要指出的是,满足积累条件后只能说明目前航天运动特征相对平稳,其可能处于机动段落,也可能处于连续常值推力机动的段落。
步骤S3.2、采用逐点滑动的方式更新窗口数据,将窗口按时间分成前后两段,分别计算半长轴变化率,并通过检测两个变化率的差异是否符合半长轴变化率的统计特征,进行连续推力机动检测;
当滤波获得窗口右侧采样时刻半长轴估值ai+N后,将数据分为和两部分,若下式成立,则认为时刻f+[0.5N]航天器的机动特性发生变化(连续机动开始或结束)
其中,[·]表示取整运算,α为第二阈值,反映了检测的灵敏度,而为了兼顾测量异常的影响,通常取值为3~6。
公式(12)只能给出机动特性变化的判断,进一步,可以通过前后窗口比对判断时刻f+[0.5N]是机动起始还是终止。具体的如下
为了保证检测可靠性,仍采用步骤S2中Step 2.2的多点判读方式。当连续有四分之一的窗口数据(即N/4个采样点)满足公式(12)时,则认为航天器在该时段内的机动特性发生了改变,将第一个时刻f+[0.5N]作为变化起始时刻,窗口向右滑动一个采样点,转Step3.1,重新开始窗口初始积累;否则,窗口向右滑动一个采样点,转Step3.3。
步骤S3.3、若未发生机动,则利用新窗口数据计算半长轴变化率的统计特征。
当判断航天器未发生机动后,重新计算半长轴变化率计算方法见公式(9),并采用下(公式14)计算半长轴变化率标准差
统计特征更新后,转入Step 3.2,开始下一时刻的机动检测。
本发明采用的方法可以实现航天器脉冲和连续推力机动的检测与辨识,提高航天器跟踪的精度和可靠性。以某典型卫星轨道为例,给出本发明方法的一个验证实例。卫星的初始轨道根数和机动特性如表1所示,图2、3给出了航天器机动检测的示意图,检测结果如表2所示。可见在本算例中,方法的机动检测率达到了100%,其中脉冲机动开始时刻的检测并无明显滞后,机动阶数时刻相差在几秒以内;连续推力机动的开始和结束时刻滞后小于100s,远低于机动持续时间(1/10以下)。因此,算例验证了本发明方法的可行性和有效性。
表1:
表2:
说明书中描述的只是该发明的具体实施方式。虽然结合附图描述了本发明的实施方式,但是本领域内熟练的技术人员可以在所附权利要求的范围内做出各种变形或修改。
Claims (1)
1.一种基于半长轴滤波值连续性的航天器轨道机动检测方法,其特征在于,包括:
步骤S1、基于运动学模型的轨道半长轴估计;
在每个采样时刻,航天器跟踪***利用对轨道运动学模型进行滤波,实时解算航天器位置与速度,进而获得该时刻的半长轴估计值;
所述运动学模型为当前统计模型;设tk时刻航天器惯性系位置为Xk=(xk,yk,zk)T,其中,xk,yk,zk分别为在惯性系x,y,z轴的分量,设tk时刻航天器惯性系速度其中,为航天器tk时刻x方向的速度,为航天器tk时刻y方向的速度,为航天器k时刻z方向的速度;半长轴估计值ak计算公式如下:
其中,μ为地球引力常数,μ=398600.4418e9;
步骤S2、基于轨道半长轴连续性的脉冲机动检测;
利用窗口累积的轨道半长轴估值进行航天器脉冲机动检测,具体以下步骤:
步骤S2.1、利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴的统计特征;
步骤S2.2、采用逐点滑动的方式更新窗口数据,并通过检测新采样时刻的半长轴数据是否符合半长轴的统计特征,进行脉冲机动检测;
步骤S2.3、若未发生机动,则利用新窗口数据计算半长轴的统计特征;
步骤S3、基于轨道半长轴变化率的连续推力机动检测;
在脉冲机动检测基础上,利用窗口累积的轨道半长轴变化率估值进行航天器连续推力机动检测,具体包括以下步骤:
步骤S3.1、若未发生脉冲机动,则利用差分检测的方法实现窗口半长轴估值数据的累积,并计算半长轴变化率的统计特征;
步骤S3.2、采用逐点滑动的方式更新窗口数据,将窗口按时间分成前后两段,分别计算半长轴变化率,并通过检测两个变化率的差异是否符合半长轴变化率的统计特征,进行连续推力机动检测;
步骤S3.3、若未发生机动,则利用新窗口数据计算半长轴变化率的统计特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510027623.2A CN104596520B (zh) | 2015-01-20 | 2015-01-20 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510027623.2A CN104596520B (zh) | 2015-01-20 | 2015-01-20 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104596520A CN104596520A (zh) | 2015-05-06 |
CN104596520B true CN104596520B (zh) | 2017-07-28 |
Family
ID=53122456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510027623.2A Expired - Fee Related CN104596520B (zh) | 2015-01-20 | 2015-01-20 | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104596520B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501815B (zh) * | 2016-09-26 | 2021-08-27 | 中国人民解放军63920部队 | 一种仅天基测角跟踪的空间目标轨道机动融合检测方法 |
CN109189102B (zh) * | 2018-11-23 | 2021-10-01 | 上海航天控制技术研究所 | 一种星上高精度计算双星半长轴偏差的方法 |
CN113326620B (zh) * | 2021-06-03 | 2022-04-08 | 中国人民解放军32039部队 | 测距数据处理方法及装置 |
CN114690210A (zh) * | 2022-04-18 | 2022-07-01 | 山东大学 | 一种基于多普勒观测值的北斗卫星机动探测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101354251A (zh) * | 2008-09-12 | 2009-01-28 | 航天东方红卫星有限公司 | 一种深空探测器等效转移轨道确定方法 |
CN101893712A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 用于地球静止卫星精密定轨的选权拟合方法 |
CN102322862A (zh) * | 2011-06-29 | 2012-01-18 | 航天东方红卫星有限公司 | 一种编队飞行卫星绝对和相对轨道确定方法 |
CN102935898A (zh) * | 2012-11-30 | 2013-02-20 | 清华大学 | 空间纯引力轨道万有引力摄动在轨飞行验证方法 |
CN103676955A (zh) * | 2013-12-19 | 2014-03-26 | 北京航空航天大学 | 一种实现分布式编队飞行的卫星自主轨道控制*** |
CN103940431A (zh) * | 2014-04-11 | 2014-07-23 | 北京空间飞行器总体设计部 | 基于gnss精密定轨的圆轨道切向小推力在轨标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120026034A1 (en) * | 2010-08-02 | 2012-02-02 | Honda Motor Co., Ltd. | Position calculation method and apparatus with gps |
-
2015
- 2015-01-20 CN CN201510027623.2A patent/CN104596520B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101354251A (zh) * | 2008-09-12 | 2009-01-28 | 航天东方红卫星有限公司 | 一种深空探测器等效转移轨道确定方法 |
CN101893712A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 用于地球静止卫星精密定轨的选权拟合方法 |
CN102322862A (zh) * | 2011-06-29 | 2012-01-18 | 航天东方红卫星有限公司 | 一种编队飞行卫星绝对和相对轨道确定方法 |
CN102935898A (zh) * | 2012-11-30 | 2013-02-20 | 清华大学 | 空间纯引力轨道万有引力摄动在轨飞行验证方法 |
CN103676955A (zh) * | 2013-12-19 | 2014-03-26 | 北京航空航天大学 | 一种实现分布式编队飞行的卫星自主轨道控制*** |
CN103940431A (zh) * | 2014-04-11 | 2014-07-23 | 北京空间飞行器总体设计部 | 基于gnss精密定轨的圆轨道切向小推力在轨标定方法 |
Non-Patent Citations (6)
Title |
---|
A New Extension of the Kalman Filter to Nonlinear Systems;Simon J.Julier等;《Proceedings of SPIE - The International Society》;19991231;182-193 * |
不完备轨道信息下的LEO轨道面内机动检测方法;于大腾等;《宇航学报》;20130331;第34卷(第3期);314-319 * |
单天基平台测向的空间非合作目标跟踪;刘也等;《飞行器测控学报》;20140731;第33卷(第4期);328-333 * |
基于预报偏差的LEO航天器轨道异常检测;张涛涛等;《中国空间科学技术》;20121031(第5期);40-46.68 * |
编队卫星轨道参数的量测融合估计;刘也等;《电子与信息学报》;20071130;第29卷(第11期);2679-2682 * |
面向任务要求的编队轨道设计及基线性能分析;张锦绣等;《上海交通大学学报》;20070728(第07期);1195-1199 * |
Also Published As
Publication number | Publication date |
---|---|
CN104596520A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108333569B (zh) | 一种基于phd滤波的异步多传感器融合多目标跟踪方法 | |
Madhukar et al. | State estimation using extended kalman filter and unscented kalman filter | |
CN107145645B (zh) | 带不确定冲击的非平稳退化过程剩余寿命预测方法 | |
CN104596520B (zh) | 一种基于半长轴滤波值连续性的航天器轨道机动检测方法 | |
US9617018B2 (en) | Automated detection and characterization of earth-orbiting satellite maneuvers | |
CN110542885B (zh) | 一种复杂交通环境下的毫米波雷达目标跟踪方法 | |
CN110687555B (zh) | 一种导航卫星原子钟弱频率跳变在轨自主快速检测方法 | |
CN109886305B (zh) | 一种基于gm-phd滤波的多传感器非顺序量测异步融合方法 | |
JP2996956B1 (ja) | 追尾装置によるロケット軌道推定法、ロケット未来位置予測法、ロケット識別法、ロケット状況検知法 | |
CN103743401A (zh) | 基于多模型航迹质量的异步融合方法 | |
CN108344981B (zh) | 面向杂波的多传感器异步检测tsbf多目标跟踪方法 | |
Lee et al. | Interacting multiple model estimation for spacecraft maneuver detection and characterization | |
CN111127523A (zh) | 基于量测迭代更新的多传感器gmphd自适应融合方法 | |
CN115600051B (zh) | 基于短弧天基光学观测的轨道机动智能检测方法和装置 | |
Duraisamy et al. | Track level fusion algorithms for automotive safety applications | |
CN111623703A (zh) | 一种基于新型卡尔曼滤波的北斗变形监测实时处理方法 | |
CN108120452B (zh) | Mems陀螺仪动态数据的滤波方法 | |
CN103973263A (zh) | 一种新的逼近滤波方法 | |
CN110989341B (zh) | 一种约束辅助粒子滤波方法及目标跟踪方法 | |
Huang et al. | A bank of maximum a posteriori estimators for single-sensor range-only target tracking | |
CN111340853B (zh) | 基于ospa迭代的多传感器gmphd自适应融合方法 | |
CN116541668B (zh) | 一种游泳划水次数确定方法、装置、设备及存储介质 | |
Greaves et al. | Maneuver Detection for Cislunar Vehicles Using Optical Measurements | |
Li et al. | A new adaptive unscented Kalman filter based on covariance matching technique | |
Dubois et al. | Performance evaluation of a moving horizon estimator for multi-rate sensor fusion with time-delayed measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170728 |
|
CF01 | Termination of patent right due to non-payment of annual fee |