CN105588568A - 火箭发射无动力飞行段弹道外推与滤值方法 - Google Patents
火箭发射无动力飞行段弹道外推与滤值方法 Download PDFInfo
- Publication number
- CN105588568A CN105588568A CN201510917868.2A CN201510917868A CN105588568A CN 105588568 A CN105588568 A CN 105588568A CN 201510917868 A CN201510917868 A CN 201510917868A CN 105588568 A CN105588568 A CN 105588568A
- Authority
- CN
- China
- Prior art keywords
- data
- trajectory
- extrapolation
- ballistic
- rocket
- 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
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/20—Instruments for performing navigational calculations
Abstract
本发明公开了一种火箭发射无动力飞行段弹道外推与滤值方法,根据一组可信的弹道测量数据或者两组位置/速度矢量数据,计算其瞬时入轨开普勒参数,然后根据轨道参数外推该时刻以后的弹道数据。根据多组外推数据按照可信度加权,计算出弹道数据结果,该结果可以用于弹道的显示以及实时测量弹道数据的异常点过滤。采用本发明的火箭无动力飞行段基于航天器定轨方法的弹道外推及滤值方法,显著提升了弹道显示的效果和精度。
Description
技术领域
本发明涉及火箭弹道处理与显示技术,特别是火箭发射无动力飞行段弹道外推与滤值方法,适用于火箭在无动力段飞行过程中依据已有测量数据(至少一组位置速度矢量或者两组时间位置或时间速度矢量),计算其瞬时入轨参数,然后进行弹道数据外推以及对后续实时测量数据的异常点过滤,采用高精度的定轨方法生成的弹道外推结果和滤值结果精度较传统方法有很大提升。
背景技术
在火箭飞行过程中,根据其飞行特点,经常会有发动机关机处于无动力飞行的滑行段,以保证火箭关键飞行动作能在测控范围内进行和节省火箭燃料。当火箭处于测控范围外时,无法获取其飞行弹道数据,为了保证数据的显示需求,必须对火箭弹道进行补弹道,传统的方法是根据处于测控范围内的速度作为其处于测控范围外的速度,理论弹道作为其处于测控范围外的弹道,显示效果及精度较差。在测控范围内由于测控设备的不稳定性,会出现一些数据异常点,如何对异常点进行过滤也是弹道处理的一个重要部分。在此背景下,我们提出一种火箭无动力飞行段基于航天器定轨方法的弹道外推及滤值方法,显著提升了弹道显示的效果和精度。
发明内容
本发明的目的是提供一种火箭在无动力飞行段的弹道数据外推及测量数据异常点过滤方法,方法概括:根据一组可信的弹道测量数据或者两组位置/速度矢量数据,计算其瞬时入轨开普勒参数,然后根据轨道参数外推该时刻以后的弹道数据。根据多组外推数据按照可信度加权,计算出弹道数据结果,该结果可以用于弹道的显示以及实时测量弹道数据的异常点过滤。
一种火箭发射无动力飞行段弹道外推与滤值方法,包括:
(1)从火箭飞行过程中的已有测量数据中选用一组可信的弹道测量数据(包含时间、位置矢量、速度矢量)或者两组位置/速度数据(包含时间、位置/速度矢量),计算火箭无动力段入轨的轨道根数;
(2)根据步骤(1)计算出的轨道根数外推以后某一时刻t的弹道数据(包含时间、位置矢量、速度矢量);
(3)另选若干组可信的弹道测量数据或位置/速度数据重复步骤(1),根据步骤(2)外推与步骤(2)相同时刻t的弹道数据;
(4)根据步骤(3)外推获得的时刻t的多组弹道数据,按照可信度加权方法(弹道数据外推时间差越小,可信度权值越高)计算出时刻t的最终的弹道数据点,作为显示和异常点过滤基础数据;
(5)根据弹道数据外推时间差确定外推数据精度标准差σ,按照3σ对实时测量数据进行异常点过滤,如果在范围内则认为该测量数据有效,将实时测量数据输出弹道,如果在误差范围外,则过滤实时测量该数据,将步骤(4)的计算结果作为测量数据值,输出弹道。
本发明方法的特点:此方法利用火箭发动机关机后的无动力飞行段受力相对简单,进行短时间弹道外推时仅需要考虑地球引力影响从而能够计算出入轨的参数,通过弹道外推和异常点过滤方法,避免了将飞行数据延后处理,从而实现火箭发射实时弹道数据处理。通过多组数据的可信度加权处理,能够降低测量设备所带来的系统误差和偶然误差。
通过本发明的方法处理后,可以显著提升弹道显示的效果和精度。
附图说明
图1为本发明的数据流程图;
图2为未经过过滤的传统补弹道数据的二维曲线图;
图3为未经过过滤的传统补弹道数据的三维曲线图;
图4为经过本发明计算得到的弹道数据的三维曲线图。
具体实施方式
参见图1,本发明方法具体包括:
(1)选用一组可信的弹道测量数据或者两组位置/速度数据计算无动力段入轨的轨道根数Kep:其中可信的弹道测量数据指在弹道曲线中,数据较为平滑,不属于异常点的数据,可以由外部其它数据处理部件提供。
通过一组弹道测量数据(r,v)计算轨道根数r表示位置矢量,v表示速度矢量。
a、计算面积速度矢量h以及它的模h:h=r×v;
b、计算过渡矢量W=h/h;
c、计算轨道倾角i;Wx、Wy、Wz分别为过渡矢量W的X、Y、Z方向的分量;
d、计算升交点赤经Ω;
e、计算过程值半通经p,其中GMe为万有引力常量与地球质量之积;
f、计算轨道半长轴a;r,v分别表示r和v的模,如果计算的结果是标量,那么r和v都是标量,也就是r和v的模,除了步骤a,b,其它都是采用标量计算;
g、计算过程值平运动角速度n;
h、计算轨道偏心率e;
i、计算偏近点角E;
j、计算平近点角M;M=E-esinE
k、计算纬度幅角u;x,y,z分别表示位置矢量r在x、y、z方向的分量;
l、计算真近点角v;
m、计算近地点幅角ω=u-v;
(2)根据计算出的轨道根数外推以后某时刻的弹道数据
通过轨道根数Kep(a,e,i,Ω,ω,M)以及与瞬时初轨时间差△t计算位置速度值(r,v):
a、计算中间值P、Q;
b、解开普勒方程E;
c、计算(t,r,v)值:
(3)选用其他可信测量数据重复(1)步骤,根据(2)外推同一时刻的弹道数据点(t,ri,vi);
(4)根据(2)和(3)步骤外推的多个同时刻弹道数据结果,按照可信度加权方法(时间差越小,可信度权值越高)计算出最终的弹道数据点,作为显示和异常点过滤基础数据(t,r,v);r=∑kiri,v=∑kivi
(5)根据弹道数据外推时间差确定外推数据精度标准差σ,按照3σ对实时测量数据进行异常点过滤,如果在范围内则认为该测量数据有效,如果在误差范围外,则过滤该数据点,并将(4)步骤计算结果作为测量值。
通过图2-图3可明显得出本发明的效果:由图2、图3可知,在测控设备跟踪即将结束时,获取的弹道数据质量比较差,而且在跟踪结束后无法获取新的弹道数据,采用传统的补弹道方法生成的弹道数据误差比较大,显示效果差;由图4可知,经过本发明处理的弹道数据比较平滑,数据精度高,能够满足火箭无动力飞行段测控盲区弹道显示需求。以上说明火箭发射无动力飞行段弹道外推与滤值方法是有效的。
Claims (2)
1.火箭发射无动力飞行段弹道外推与滤值方法,包括如下步骤:
(1)从火箭飞行过程中的已有测量数据中选用一组可信的弹道测量数据或者两组位置/速度数据,计算火箭无动力段入轨的轨道根数;
(2)根据步骤(1)计算出的轨道根数外推以后某一时刻t的弹道数据;
(3)另选若干组可信的弹道测量数据或位置/速度数据重复步骤(1),根据步骤(2)外推与步骤(2)相同时刻t的弹道数据;
(4)根据步骤(3)外推获得的时刻t的多个弹道数据,按照可信度加权方法计算出时刻t的最终的弹道数据点,作为显示和异常点过滤基础数据;
(5)根据弹道数据外推时间差确定外推数据精度标准差σ,按照3σ对实时测量数据进行异常点过滤,如果在范围内则认为该测量数据有效,将实时测量数据输出弹道,如果在误差范围外,则过滤实时测量该数据,将步骤(4)的计算结果作为测量数据值,输出弹道。
2.根据权利要求1所述火箭发射无动力飞行段弹道外推与滤值方法,其特征在于:所述可信度加权方法具体依据是:时间差越小,可信度权值越高。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917868.2A CN105588568B (zh) | 2015-12-11 | 2015-12-11 | 火箭发射无动力飞行段弹道外推与滤值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917868.2A CN105588568B (zh) | 2015-12-11 | 2015-12-11 | 火箭发射无动力飞行段弹道外推与滤值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105588568A true CN105588568A (zh) | 2016-05-18 |
CN105588568B CN105588568B (zh) | 2018-09-18 |
Family
ID=55928329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510917868.2A Expired - Fee Related CN105588568B (zh) | 2015-12-11 | 2015-12-11 | 火箭发射无动力飞行段弹道外推与滤值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105588568B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416093A (zh) * | 2018-01-30 | 2018-08-17 | 重庆大学 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
CN111881514A (zh) * | 2020-07-16 | 2020-11-03 | 北京航天自动控制研究所 | 一种发动机故障状态下的制导重构方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060169831A1 (en) * | 2004-01-06 | 2006-08-03 | Humphrey John M | Ballistic target defense system and methods |
CN105022035A (zh) * | 2015-07-31 | 2015-11-04 | 中国电子科技集团公司第三十八研究所 | 一种基于模型修正的弹道目标发射点估计装置及其方法 |
-
2015
- 2015-12-11 CN CN201510917868.2A patent/CN105588568B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060169831A1 (en) * | 2004-01-06 | 2006-08-03 | Humphrey John M | Ballistic target defense system and methods |
CN105022035A (zh) * | 2015-07-31 | 2015-11-04 | 中国电子科技集团公司第三十八研究所 | 一种基于模型修正的弹道目标发射点估计装置及其方法 |
Non-Patent Citations (2)
Title |
---|
俞建国等: "弹道目标航迹片段关联及优化", 《航空学报》 * |
杨海峰等: "《遥感卫星的轨道外推方法》", 《全球定位系统》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416093A (zh) * | 2018-01-30 | 2018-08-17 | 重庆大学 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
CN108416093B (zh) * | 2018-01-30 | 2021-10-12 | 重庆大学 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
CN111881514A (zh) * | 2020-07-16 | 2020-11-03 | 北京航天自动控制研究所 | 一种发动机故障状态下的制导重构方法 |
CN111881514B (zh) * | 2020-07-16 | 2024-04-05 | 北京航天自动控制研究所 | 一种发动机故障状态下的制导重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105588568B (zh) | 2018-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104457705B (zh) | 基于天基自主光学观测的深空目标天体初定轨方法 | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
Johnson et al. | Real-time terrain relative navigation test results from a relevant environment for Mars landing | |
CN105160125B (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN103235306B (zh) | 一种适用于高速机动飞行器载sar成像的运动补偿方法 | |
CN102620748B (zh) | 捷联惯导系统晃动基座条件下杆臂效应的估计和补偿方法 | |
CN104657559B (zh) | 基于圆柱型速度参数截面的地月自由返回轨道设计方法 | |
CN105021092A (zh) | 一种捷联寻的导引头的制导信息提取方法 | |
CN105974415A (zh) | 一种机载sar方位空变运动误差的高精度补偿方法 | |
CN109269511B (zh) | 未知环境下行星着陆的曲线匹配视觉导航方法 | |
CN106468554A (zh) | 一种非接触式的翻滚卫星的惯性参数的测算方法 | |
CN107255474A (zh) | 一种融合电子罗盘和陀螺仪的pdr航向角确定方法 | |
CN102944241A (zh) | 基于多胞型线性微分包含的航天器相对姿态确定方法 | |
CN105606846A (zh) | 一种基于姿态信息的加速度计校准方法 | |
CN103927289A (zh) | 一种依据天基卫星测角资料确定低轨目标卫星初始轨道的方法 | |
CN103017772A (zh) | 一种基于可观性分析的光学和脉冲星融合自主导航方法 | |
CN101929862A (zh) | 基于卡尔曼滤波的惯性导航系统初始姿态确定方法 | |
Changey et al. | Experimental validation | |
CN102944238A (zh) | 一种行星探测器接近目标过程中相对位置确定方法 | |
CN104864875B (zh) | 一种基于非线性h∞滤波的航天器自主定位方法 | |
Hill et al. | Ground-to-air flow visualization using Solar Calcium-K line Background-Oriented Schlieren | |
CN105588568A (zh) | 火箭发射无动力飞行段弹道外推与滤值方法 | |
CN106705967A (zh) | 一种基于行人航位推算的精度改善的室内定位和方法 | |
CN103364842B (zh) | 一种捷联式航空重力仪误差分离方法 | |
CN105333873B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180918 Termination date: 20191211 |
|
CF01 | Termination of patent right due to non-payment of annual fee |