CN110262241B - 基于高斯过程预测控制的航天器轨道控制方法 - Google Patents
基于高斯过程预测控制的航天器轨道控制方法 Download PDFInfo
- Publication number
- CN110262241B CN110262241B CN201910578209.9A CN201910578209A CN110262241B CN 110262241 B CN110262241 B CN 110262241B CN 201910578209 A CN201910578209 A CN 201910578209A CN 110262241 B CN110262241 B CN 110262241B
- Authority
- CN
- China
- Prior art keywords
- satellite
- orbit
- control
- relative
- gaussian process
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/041—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a variable is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
Description
技术领域
本发明涉及一种卫星自主轨道控制技术领域,特别是一种基于高斯过程预测控制的航天器轨道控制方法。
背景技术
近年来,针对静止轨道卫星星座的研究呈现迅猛的态势。星座中的每颗卫星由于初始入轨点存在误差或者运行过程中受轨道摄动因素影响会存在偏差,运行一段时间之后,卫星便会偏离所设计的运行轨道,可能导致星座失效,甚至卫星间的碰撞。为使静止轨道卫星维持在一定范围内运行,必须对卫星进行控制,将此称为位置保持。
当前对于静止轨道卫星的位置保持控制方法主要有脉冲位置保持和小推力位置保持两种方式。脉冲位置保持方法是由地面站发出指令,约15天进行一次东西位置保持机动,2个月左右进行一次南北位置保持机动。但是这种控制方式精度不是很高,而且脉冲位置保持方法适用于配备化学推进系统的卫星,对于电推进系统的卫星,脉冲假设不成立,容易产生误差。针对小推力的系统的卫星,常见的是基于优化算法的控制策略,在算法设计过程对于卫星的所受的摄动多采用理想的数值模型,但实际中卫星所受到的摄动影响十分复杂,难以精确建模,在实际操作中会产生各种各样的误差,这就导致很多算法多停留在理论验证阶段。
发明内容
本发明的目的在于克服现有技术的上述不足,提供基于高斯过程预测控制的航天器轨道控制方法,提高了静止轨道卫星位保方案的控制精度,且燃料消耗并无显著增加。
本发明的上述目的是通过如下技术方案予以实现的:
基于高斯过程预测控制的航天器轨道控制方法,包括如下步骤
步骤一、获取卫星所在轨道当前的轨道参数;根据轨道参数计算t时刻卫星与卫星位置保持误差范围中心点之间的相对轨道根数x(t);
步骤二、将相对轨道根数x(t)转换成高斯变分方程;并将高斯变分方程简化,得到卫星相对轨道根数状态方程对卫星相对轨道根数状态方程进行离散,得到设计控制器所需要的卫星相对轨道根数离散迭代模型x(k),其中k为离散后的当前时刻;
步骤三、采用高斯过程机器学习的方法对摄动项进行学习和预测;记录历史控制量以及卫星的历史相对轨道根数,设当前时刻为k,取前Nt个时刻的摄
步骤四、设定目标函数为J;设Xlb为卫星相对轨道根数的下界;设定Xub为卫星相对轨道根数的上界;即卫星位置保持误差范围的边界;设定Ulb为控制量的下界;设定Uub为控制量上界;即卫星加速度的约束条件;根据Xlb、Xub、Ulb和Uub,求解一组控制量序列使得目标函数J最小;则该组控制量序列即为最优的控制量序列;
在上述的基于高斯过程预测控制的航天器轨道控制方,所述的步骤一中,轨道参数包括轨道的半长轴a、偏心率e、倾角i,升交点赤经Ω和近地点幅角ω。
在上述的基于高斯过程预测控制的航天器轨道控制方,所述的步骤一中,相对轨道根数x(t)的计算方法为:
式中,E1为卫星与卫星位置保持误差盒中心点之间的相对角速率;
E2、E3分别为偏心率向量的两个分量;
E4、E5分别为倾角向量的两个分量;
E6为卫星与卫星位置保持误差盒中心点之间的平经度误差;
n为卫星轨道的平均角速率;
ne为地球自转角速率;
ε为平经度;
t0为初始时刻;
α为卫星位置保持误差盒中心点的赤经。
在上述的基于高斯过程预测控制的航天器轨道控制方,所述的步骤二中,将相对轨道根数x(t)转换成高斯变分方程的方法为:
式中,υ为卫星轨道的真近点角;
p为卫星轨道的半通径;
λ为卫星轨道的真经度;
r为卫星轨道半径;
E为卫星的偏近点角;
ur为卫星在RTN卫星本体坐标系oxyz下x轴的加速度;
ut为卫星在RTN卫星本体坐标系oxyz下y轴的加速度;
un为卫星在RTN卫星本体坐标系oxyz下z轴的加速度。
在上述的基于高斯过程预测控制的航天器轨道控制方,所述的RTN卫星本体坐标系oxyz的建立方法为:
卫星质心为原点o;x轴方向为指向地心方向;y轴方向为指向卫星速度方向;z轴方向由右手定则确定。
式中,A为时不变系统矩阵;
α(t)为t时刻卫星位置保持误差中心点的赤经;
B(*)为时变的控制矩阵;
uc(t)为t时刻卫星的控制输入;
ud(t)为t时刻卫星的摄动输入。
在上述的基于高斯过程预测控制的航天器轨道控制方,步骤二中,卫星相对轨道根数离散迭代模型x(k):
uc(k-1)为上一时刻的控制输入;
在上述的基于高斯过程预测控制的航天器轨道控制方,步骤四中,目标函数J为:
式中,x(i)为i时刻卫星的相对轨道根数;
uc(i)为i时刻卫星的应的控制量;
Q为x(i)权值矩阵;
R为uc(i)权值矩阵;
T为转置。
本发明与现有技术相比具有如下优点:
(1)本发明采用了相对轨道根数的卫星轨道表示方法,简化了控制器结构,降低了控制器设计难度;
(2)本发明采用了模型预测控制器进行静止轨道卫星的位置保持控制,能够有效的处理多种约束,提高卫星的控制精度;
(3)本发明利用高斯过程进行摄动项的学习与预测,提高卫星轨道预报精度。
附图说明
图1为本发明航天器轨道控制流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
本发明提供一种基于高斯过程预测控制的航天器轨道控制方法,在自主导航结果的基础上,获得轨道当前的轨道平根数。我们将误差盒的中心点看作是一颗虚拟卫星,通过计算可以得到实际卫星同虚拟卫星间的相对轨道根数模型。将相对轨道根数模型写为高斯变分方程的形式,通过化简可以得到一个具有时不变状态矩阵A的非线性状态方程,其中输入部分包括控制输入和摄动输入,分别为控制输入量与摄动输入量在卫星RTN方向的控制/摄动加速度。接下来,通过对状态方程的离散化可以得到系统的迭代模型。
利用高斯过程对摄动项以及不确定性进行学习。卫星运行过程中会受到各种摄动力的影响,由于其数学形式复杂,且具有不确定性,难以精确建模。本发明采用了机器学习的方法,利用高斯过程,对系统的摄动因素以及不确定性进行学习,同时估计和预测出未来一段时间的摄动值,用于补偿预测控制器中的模型误差。
本发明基于模型预测控制设计静止轨道卫星位置保持的控制器,利用离散化的迭代模型,对系统未来N步内的状态进行预测,并且在考虑系统约束的情况下,求解优化问题。通过求解最优控制序列从而得到当前时刻的控制量。
如图1所示为航天器轨道控制流程图,由图可知,基于高斯过程预测控制的航天器轨道控制方法,包括如下步骤
步骤一、获取卫星所在轨道当前的轨道参数;轨道参数包括轨道的半长轴 a、偏心率e、倾角i,升交点赤经Ω和近地点幅角ω。根据轨道参数计算t时刻卫星与卫星位置保持误差范围中心点之间的相对轨道根数x(t);静止轨道卫星位置保持目标是控制卫星在以定点位置为中心的一定的空间范围内,我们称这个空间范围为卫星位置保持的误差范围。利用获得的卫星的轨道参数计算当前卫星与卫星位置保持误差范围中心点之间的相对轨道根数。
相对轨道根数x(t)的计算方法为:
式中,E1为卫星与卫星位置保持误差盒中心点之间的相对角速率;
E2、E3分别为偏心率向量的两个分量;
E4、E5分别为倾角向量的两个分量;
E6为卫星与卫星位置保持误差盒中心点之间的平经度误差;
n为卫星轨道的平均角速率;
ne为地球自转角速率;
ε为平经度;
t0为初始时刻;
α为卫星位置保持误差盒中心点的赤经。
步骤二、将相对轨道根数x(t)转换成高斯变分方程;并将高斯变分方程简化,得到卫星相对轨道根数状态方程对卫星相对轨道根数状态方程进行离散,得到设计控制器所需要的卫星相对轨道根数离散迭代模型x(k),其中k为离散后的当前时刻;
将相对轨道根数x(t)转换成高斯变分方程的方法为:
式中,υ为卫星轨道的真近点角;
p为卫星轨道的半通径;
λ为卫星轨道的真经度;
r为卫星轨道半径;
E为卫星的偏近点角;
ur为卫星在RTN卫星本体坐标系oxyz下x轴的加速度;
ut为卫星在RTN卫星本体坐标系oxyz下y轴的加速度;
un为卫星在RTN卫星本体坐标系oxyz下z轴的加速度。
RTN卫星本体坐标系oxyz的建立方法为:
卫星质心为原点o;x轴方向为指向地心方向;y轴方向为指向卫星速度方向;z轴方向由右手定则确定。
式中,A为时不变系统矩阵;
α(t)为t时刻卫星位置保持误差中心点的赤经;
B(*)为时变的控制矩阵;
uc(t)为t时刻卫星的控制输入;
ud(t)为t时刻卫星的摄动输入。
卫星相对轨道根数离散迭代模型x(k)为:
uc(k-1)为上一时刻的控制输入;
步骤四、设定目标函数为J;设Xlb为卫星相对轨道根数的下界;设定Xub为卫星相对轨道根数的上界;即卫星位置保持误差范围的边界;设定Ulb为控制量的下界;设定Uub为控制量上界;即卫星加速度的约束条件;那么如此以来,静止轨道卫星的位置保持控制问题可以转化为以下的一个优化问题,即根据Xlb、 Xub、Ulb和Uub,求解一组控制量序列使得目标函数J最小;则该组控制量序列即为最优的控制量序列;上述优化问题可以将其转化为标准的二次规划形式进行求解。
目标函数J为:
式中,x(i)为i时刻卫星的相对轨道根数;
uc(i)为i时刻卫星的应的控制量;
Q为x(i)权值矩阵;
R为uc(i)权值矩阵;
T为转置。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (7)
1.基于高斯过程预测控制的航天器轨道控制方法,其特征在于:包括如下步骤
步骤一、获取卫星所在轨道当前的轨道参数;根据轨道参数计算t时刻卫星与卫星位置保持误差范围中心点之间的相对轨道根数x(t);
步骤二、将相对轨道根数x(t)转换成高斯变分方程;并将高斯变分方程简化,得到卫星相对轨道根数的状态方程;对卫星相对轨道根数状态方程进行离散,得到设计控制器所需要的卫星相对轨道根数离散迭代模型x(k),其中k为离散后的当前时刻;
uc(k-1)为上一时刻的控制输入;
步骤四、设定目标函数为J;设Xlb为卫星相对轨道根数的下界;设定Xub为卫星相对轨道根数的上界;Xlb、Xub即卫星位置保持误差范围的边界;设定Ulb为控制量的下界;设定Uub为控制量上界;Ulb、Uub即卫星加速度的约束条件;根据Xlb、Xub、Ulb和Uub,求解一组控制量序列使得目标函数J最小;则该组控制量序列即为最优的控制量序列;其中,目标函数J为:
式中,x(i)为i时刻卫星的相对轨道根数;
uc(i)为i时刻卫星的应的控制量;
Q为x(i)权值矩阵;
R为uc(i)权值矩阵;
T为转置;
2.根据权利要求1所述的基于高斯过程预测控制的航天器轨道控制方法 ,其特征在于:所述的步骤一中,轨道参数包括轨道的半长轴a、偏心率e、倾角i,升交点赤经Ω和近地点幅角ω。
5.根据权利要求4所述的基于高斯过程预测控制的航天器轨道控制方法 ,其特征在于:所述的RTN卫星本体坐标系oxyz的建立方法为:
卫星质心为原点o;x轴方向为指向地心方向;y轴方向为指向卫星速度方向;z轴方向由右手定则确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910578209.9A CN110262241B (zh) | 2019-06-28 | 2019-06-28 | 基于高斯过程预测控制的航天器轨道控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910578209.9A CN110262241B (zh) | 2019-06-28 | 2019-06-28 | 基于高斯过程预测控制的航天器轨道控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110262241A CN110262241A (zh) | 2019-09-20 |
CN110262241B true CN110262241B (zh) | 2022-03-04 |
Family
ID=67923094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910578209.9A Active CN110262241B (zh) | 2019-06-28 | 2019-06-28 | 基于高斯过程预测控制的航天器轨道控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110262241B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111547274B (zh) * | 2020-03-19 | 2023-08-29 | 上海航天控制技术研究所 | 一种航天器高精度自主目标预报方法 |
CN113641949B (zh) * | 2021-08-05 | 2023-03-28 | 中国西安卫星测控中心 | 一种地球同步转移段轨道根数高精度拟合方法 |
CN114355976B (zh) * | 2021-12-31 | 2023-12-05 | 中山大学 | 一种基于学习的控制无人机在风扰下完成轨迹跟踪的方法 |
CN114967453B (zh) * | 2022-05-25 | 2024-07-16 | 北京理工大学 | 一种基于神经网络的卫星东西位保协态初值估计方法 |
CN115320891B (zh) * | 2022-10-12 | 2023-01-24 | 北京航天驭星科技有限公司 | 一种基于虚拟卫星的近圆标称轨道控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103123488A (zh) * | 2013-01-18 | 2013-05-29 | 西北工业大学 | 空间绳系机器人系统逼近目标协调控制方法 |
CN104019818A (zh) * | 2014-06-19 | 2014-09-03 | 北京理工大学 | 一种基于预测轨迹的行星导航轨道器布局优化方法 |
CN107798187A (zh) * | 2017-10-24 | 2018-03-13 | 北京理工大学 | 一种高效卫星星座多目标优化方法 |
CN108614420A (zh) * | 2018-05-04 | 2018-10-02 | 天津大学 | 基于非线性规划的星簇级卫星容错控制方法 |
-
2019
- 2019-06-28 CN CN201910578209.9A patent/CN110262241B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103123488A (zh) * | 2013-01-18 | 2013-05-29 | 西北工业大学 | 空间绳系机器人系统逼近目标协调控制方法 |
CN104019818A (zh) * | 2014-06-19 | 2014-09-03 | 北京理工大学 | 一种基于预测轨迹的行星导航轨道器布局优化方法 |
CN107798187A (zh) * | 2017-10-24 | 2018-03-13 | 北京理工大学 | 一种高效卫星星座多目标优化方法 |
CN108614420A (zh) * | 2018-05-04 | 2018-10-02 | 天津大学 | 基于非线性规划的星簇级卫星容错控制方法 |
Non-Patent Citations (2)
Title |
---|
《连续小推力作用下航天器机动轨道设计》;孙冲;《中国博士学位论文全文数据库(电子期刊)工程科技II辑》;20190115;全文 * |
微小卫星姿轨自主确定技术研究;郁丰;《中国博士学位论文全文数据库(电子期刊)工程科技II辑》;20090515;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110262241A (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110262241B (zh) | 基于高斯过程预测控制的航天器轨道控制方法 | |
EP3429926B1 (en) | Method and system for controlling operation of spacecraft, spacecraft | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN106697333A (zh) | 一种航天器轨道控制策略的鲁棒性分析方法 | |
EP3765369B1 (en) | Spacecraft and vehicle | |
Lee et al. | Robust position and attitude control for spacecraft formation flying | |
CN106297422B (zh) | 一种基于最小二乘的非线性相对轨迹预报方法 | |
CN104252548A (zh) | 一种燃料最优的火星探测器入轨目标点设计方法 | |
Kasiri et al. | Coupled position and attitude control of a servicer spacecraft in rendezvous with an orbiting target | |
CN110209190B (zh) | 一种卫星标称轨道无偏飞行控制的方法 | |
Lippe et al. | Minimization of delta-v for satellite swarm maintenance using a virtual chief | |
Biggs et al. | In-situ tracking of a solar sail’s characteristic acceleration using a robust active disturbance estimator | |
WO2023234420A1 (en) | Eigenmotion control for near rectilinear halo orbits | |
Seelbinder | On-board trajectory computation for mars atmospheric entry based on parametric sensitivity analysis of optimal control problems | |
Leomanni | Attitude and orbit control techniques for spacecraft with electric propulsion | |
De Florio et al. | Precise autonomous orbit control in low earth orbit | |
Mehrjardi et al. | Integrated attitude-orbit dynamics and control of spacecraft systems: State of the art and future trends | |
Tsujita | A feasible study on the model predictive control for docking approach of small spacecraft using thrusters and a control moment gyro | |
Sternberg et al. | Mapping a chaser satellite's feasibility space for soft docking | |
Zhang et al. | Space target surveillance based on non-linear model predictive control | |
Ramezani et al. | Autonomous Orbital Correction for Nano Satellites Using J2 Perturbation and LSTM Networks | |
Xie et al. | Autonomous guidance, navigation, and control of spacecraft | |
Morantea et al. | Low-Thrust Trajectory Optimization and Autonomy Analysis for a Medium-Earth-Orbit Constellation Deployment | |
Laguardia | Design of an optimization algorithm for in-orbit inspection relative trajectories | |
Sternberg | Multiple model trajectory generation for uncertain target spin direction |
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 | ||
GR01 | Patent grant |