CN104599300B - 一种适应多尺度显示的航天器轨道实时绘制方法 - Google Patents
一种适应多尺度显示的航天器轨道实时绘制方法 Download PDFInfo
- Publication number
- CN104599300B CN104599300B CN201410818029.0A CN201410818029A CN104599300B CN 104599300 B CN104599300 B CN 104599300B CN 201410818029 A CN201410818029 A CN 201410818029A CN 104599300 B CN104599300 B CN 104599300B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- mrow
- msup
- display
- orbital
- 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
Links
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明公开了属于航天器轨道设计及计算机图形学技术领域的一种适应多尺度显示的航天器轨道实时绘制方法。包括轨道数据的动态坐标基表示方法和轨控参数优先显示方法,在航天仿真及飞行测控的轨道监视过程中,利用变坐标基表示、轨道线快速采样以及参数突出表现等技术方法,通过计算机图形学方式展现航天器状态,为航天测控提供信息保障。
Description
技术领域
本发明属于航天器轨道设计及计算机图形学技术领域,具体涉及一种适应多尺度显示的航天器轨道实时绘制方法。
背景技术
航天器飞行控制过程中,决策人员需要根据下传的飞行器位置、姿态等数据,监视航天器的飞行轨道、姿态等状态参数,是航天器飞行控制的关键环节之一。监视的实现途径是利用计算机图形学进行轨道曲线绘制,其中常用的方法是用地球坐标系表示轨道线,在此坐标系下绘制轨道线的每一个点,由于航天器尺度为米级,因此近距离观察航天器近景工况和状态时,地球坐标系下的轨道线会出现在航天器部位抖动的现象,造成了近景状态下不能同步显示轨道线,不利于其姿态、轨道参数的监视。对于轨道的微量控制、精确控制等具体的任务,常规的绘制显示方法,卫星的轨道参数变化量在轨道显示上变化不明显,不能有效地发挥轨控过程的监视作用。
发明内容
本发明的目的在于提供一种适应多尺度显示的航天器轨道实时绘制方法。采用了天体坐标系与航天器本体坐标系相结合的轨道参数表达方法,在航天器近景监视条件下采用航天器体坐标系进行轨道数据表示和计算机绘制,在绕飞天体尺度显示条件下采用天体坐标系进行数据处理,同时在轨控过程展示方面,提出了轨控变量优先的显示方法,保证了轨控监视的有效性。
一种适应多尺度显示的航天器轨道实时绘制方法,包括轨道数据的动态坐标基表示方法和轨控参数优先显示方法;
所述轨道数据的动态坐标基表示方法步骤为:
航天器近景显示条件下,轨道线所有数据点转移到以航天器质心为原点,轨道面所在的平面为坐标平面,航天器运动方向为X轴方向,轨道面内垂直于X轴方向远离中心天体的方向为Z方向的坐标系;
天体宏观显示条件下轨道显示,该状态下航天器采用放大显示或航天器不可见状态,轨道线所有数据点转移到以天体质心为原点的坐标系;
依据轨道参数对轨道线上的点进行采样,绘制并显示;
所述轨控参数优先显示方法的步骤为:
轨道线所有数据点转移到以航天器质心为原点,轨道面所在的平面为坐标平面,航天器运动方向为X轴方向,轨道面内垂直于X轴方向远离中心天体的方向为Z方向的坐标系;
分别以根据变轨时间和轨道参数放大倍数m、n,以及轨道变化参数Δa、Δb,确定变轨过程中时刻t的轨道形状;
对时刻t的轨道进行轨道点采样,进行轨道线绘制。
本发明的有益效果:本发明的方法实现了天体尺度轨道展示与航天器尺度轨道、姿态展示的兼容,实现了航天器轨控过程的同步展示与轨控参数的突出展示、提高了轨控监视的直观性。
附图说明
图1为本发明中航天器体坐标系下椭圆轨道的示意图;
图中,o-xy表示航天器坐标系的轨道平面,E所示实体表示航天器所饶飞的天体,a,b分别是椭圆轨道曲线的长半轴和短半轴。
图2为轨控参数优先显示方法示意图;
图中,o-xy表示航天器坐标系的轨道平面,E所示实体表示航天器所饶飞的天体,a1,b1分别是变轨前椭圆轨道曲线的长半轴和短半轴,a2,b2分别是实际变轨后椭圆轨道曲线的长半轴和短半轴,分别是优化显示的变轨后椭圆轨道曲线的长半轴和短半轴。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
实施例1
一种适应多尺度显示的航天器轨道实时绘制方法,包括轨道数据的动态坐标基表示方法和轨控参数优先显示方法;
所述轨道数据动态坐标基表示方法为:
(1)从天体坐标系到航天器体坐标系转换
利用天体坐标系下航天器运行轨道根数(a,e,i,v,φ,Ω,),计算航天器轨道椭圆曲线短半轴,其中椭圆长半轴即为轨道根数a,轨道偏心率e,则短半轴b计算公式如下:
以航天器体坐标系为基础,建立轨道曲线的描述方程,在体坐标系下的轨道线描述如下:
(2)沿着轨道曲线逐点采样,依次获取轨道线上所有点(xi,yi),其具体方法参考由另一个专利给出。
(3)将所有采样点(xi,yi),送到计算机渲染通道,将依次将所有采样点用设定的曲线连接,即可完成曲线绘制。
所述轨控参数优先显示方法为:
(1)轨道参数确定
实际变轨过程为,在时间T内,轨道参数的半长轴、半短轴(a1,b1)变化为参数(a2,b2),设定优化显示过程中半长轴和半短轴最终变化为其关系为:
Δa=(a2-a1)
Δb=(b2-b1)
其中m,n为根据优化显示的需求所确定的大于1的数值。
在变轨过程中的时刻t,轨道参数的计算方法为:
(2)沿着轨道曲线逐点采样,依次获取轨道线上所有点(xi,yi)。
(3)将所有采样点(xi,yi),送到计算机渲染通道,将依次将所有采样点用设定的曲线连接,即可完成曲线绘制。
本发明的方法在航天器运行、测量与控制、飞行轨道设计等三维可视化任务中,解决了天体尺度下航天器近景显示的轨道线抖动问题,满足了任务多种显示需求下的轨道显示问题;该实时绘制方法提高了轨道绘制的实时性、大量节省了计算资源;有效提供了变轨过程监视,为航天控制监控决策提供了重要的信息保障。该方法可应用于空间探测、轨道设计、飞行仿真等领域的三维可视化。
Claims (1)
1.一种适应多尺度显示的航天器轨道实时绘制方法,其特征在于,包括轨道数据的动态坐标基表示方法和轨控参数优先显示方法;
所述轨道数据的动态坐标基表示方法步骤为:
航天器近景显示条件下,轨道线所有数据点转移到以航天器质心为原点,轨道面所在的平面为坐标平面,航天器运动方向为X轴方向,轨道面内垂直于X轴方向远离中心天体的方向为Z方向的坐标系;
天体宏观显示条件下轨道显示,航天器采用放大显示或航天器不可见状态,轨道线所有数据点转移到以天体质心为原点的坐标系;依据轨道参数对轨道线上的点进行采样,绘制并显示;
所述轨控参数优先显示方法为:
(1)轨道参数确定
实际变轨过程为,在时间T内,轨道参数的半长轴、半短轴(a1,b1)变化为参数(a2,b2),设定优化显示过程中半长轴和半短轴最终变化为(a′2,b′2),其关系为:
Δa=(a2-a1)
Δb=(b2-b1)
a′2=a1+mΔa
b′2=b1+nΔb
其中m,n为根据优化显示的需求所确定的大于1的数值;
在变轨过程中的时刻t,轨道参数的计算方法为:
<mrow>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>a</mi>
<mi>t</mi>
<mo>/</mo>
<mi>T</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>a</mi>
<mi>t</mi>
<mo>/</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>b</mi>
<mi>t</mi>
<mo>/</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>=</mo>
<mn>1</mn>
</mrow>
(2)沿着轨道曲线逐点采样,依次获取轨道线上所有点(xi,yi);
(3)将所有采样点(xi,yi),送到计算机渲染通道,将依次将所有采样点用设定的曲线连接,即可完成曲线绘制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410818029.0A CN104599300B (zh) | 2014-12-26 | 2014-12-26 | 一种适应多尺度显示的航天器轨道实时绘制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410818029.0A CN104599300B (zh) | 2014-12-26 | 2014-12-26 | 一种适应多尺度显示的航天器轨道实时绘制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104599300A CN104599300A (zh) | 2015-05-06 |
CN104599300B true CN104599300B (zh) | 2018-03-02 |
Family
ID=53125050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410818029.0A Active CN104599300B (zh) | 2014-12-26 | 2014-12-26 | 一种适应多尺度显示的航天器轨道实时绘制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104599300B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109615134A (zh) * | 2018-12-10 | 2019-04-12 | 中国民用航空飞行学院 | 一种离场航线数据处理系统及方法、信息数据处理终端 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968401A (zh) * | 2012-11-08 | 2013-03-13 | 北京空间飞行器总体设计部 | 一种航天器精测数据处理系统 |
CN103593181A (zh) * | 2013-10-23 | 2014-02-19 | 中国运载火箭技术研究院 | 一种可配置分布式三维视景仿真系统 |
CN103605869A (zh) * | 2013-11-27 | 2014-02-26 | 北京卫星环境工程研究所 | 面向航天器总装的三维实做数字样机系统及其构建方法 |
CN103914385A (zh) * | 2014-04-11 | 2014-07-09 | 北京航空航天大学 | 面向安全苛刻系统并行测试的测试细则可视化动态构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8743064B2 (en) * | 2010-09-17 | 2014-06-03 | A.I. Solutions, Inc. | Gesture orbit design |
-
2014
- 2014-12-26 CN CN201410818029.0A patent/CN104599300B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968401A (zh) * | 2012-11-08 | 2013-03-13 | 北京空间飞行器总体设计部 | 一种航天器精测数据处理系统 |
CN103593181A (zh) * | 2013-10-23 | 2014-02-19 | 中国运载火箭技术研究院 | 一种可配置分布式三维视景仿真系统 |
CN103605869A (zh) * | 2013-11-27 | 2014-02-26 | 北京卫星环境工程研究所 | 面向航天器总装的三维实做数字样机系统及其构建方法 |
CN103914385A (zh) * | 2014-04-11 | 2014-07-09 | 北京航空航天大学 | 面向安全苛刻系统并行测试的测试细则可视化动态构建方法 |
Non-Patent Citations (3)
Title |
---|
《航天器运行原理》课程的可视化教学探索;闻新 等;《沈阳航空航天大学学报》;20140504;第27-31页 * |
空间态势中航天器轨道的可视化方法;汪荣峰 等;《计算机应用》;20130308;第32卷(第2期);第255-257页 * |
航天器飞行及测控过程可视化系统设计与实现;程委 等;《计算机仿真》;20101103;第27卷(第8期);第22-27页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104599300A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103112600B (zh) | 一种星际转移轨道设计方法 | |
WO2017113567A1 (zh) | 火星探测器自主导航方法 | |
CN104570742B (zh) | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 | |
CN104176275B (zh) | 一种使用动量轮与磁力矩器联合的速率阻尼方法 | |
CN110018691A (zh) | 小型多旋翼无人机飞行状态估计系统和方法 | |
CN105652880B (zh) | 用于飞行器大空域飞行的非线性抗饱和高度指令生成方法 | |
CN103940442A (zh) | 一种采用加速收敛算法的定位方法及装置 | |
CN108663936B (zh) | 模型不确定航天器无退绕姿态跟踪有限时间控制方法 | |
WO2021078206A1 (zh) | 流量校正方法及装置、无人飞行器 | |
CN105865455A (zh) | 一种利用gps与加速度计计算飞行器姿态角的方法 | |
CN110116409A (zh) | 一种基于扰动观测器的四通道遥操作双边控制方法 | |
CN106248082B (zh) | 一种飞行器自主导航系统及导航方法 | |
CN108021138A (zh) | 一种地磁场模型简化设计方法 | |
CN104599300B (zh) | 一种适应多尺度显示的航天器轨道实时绘制方法 | |
CN102508492B (zh) | 一种飞行器在等高航路点间的定高度大圆飞行实现方法 | |
CN103941592A (zh) | 飞行机器人动力学模型的在线建模方法 | |
CN103438892A (zh) | 一种改进的基于ekf的天文自主定轨算法 | |
CN106094848B (zh) | 一种无角速度测量对地定向空间站的角动量管理方法 | |
Ting et al. | Sliding mode control design for oblique wing aircraft in wing skewing process | |
CN112129272B (zh) | 视觉里程计的实现方法和实现装置 | |
CN104537698B (zh) | 一种曲率自适应的轨道绘制方法 | |
CN103979116A (zh) | 一种飞机平显最佳攻角框动态定位方法 | |
CN110703286B (zh) | 一种用于航空gps信号延迟的估算系统及其方法 | |
CN105955287B (zh) | 一种对地定向空间站的姿态控制和角动量管理方法 | |
CN118133724B (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 |