CN102968124B - 基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 - Google Patents
基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 Download PDFInfo
- Publication number
- CN102968124B CN102968124B CN201210499072.6A CN201210499072A CN102968124B CN 102968124 B CN102968124 B CN 102968124B CN 201210499072 A CN201210499072 A CN 201210499072A CN 102968124 B CN102968124 B CN 102968124B
- Authority
- CN
- China
- Prior art keywords
- detector
- control
- model
- error
- landing
- 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
Abstract
本发明涉及一种基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法,尤其涉及一种在行星着陆过程中探测器受到多源扰动条件下的轨迹跟踪控制方法,属于深空探测制导与控制技术领域。本方法考虑目标天体引力模型与外界干扰等不确定因素对轨迹跟踪的影响,引入模型不确定界的设计理念,通过设计鲁棒反馈控制律,并兼顾控制机构的非线性特性,确保实际状态轨迹与前馈制导设计轨迹的偏离量不超过期望的界,从而优化多源扰动下的行星探测器着陆轨迹跟踪性能;能有效地克服多源扰动和系统不确定性对探测器着陆飞行轨迹的影响,提高了探测器系统的控制性能,使行星探测器着陆飞行控制精度得到了保障。
Description
技术领域
本发明涉及一种基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法,尤其涉及一种在行星着陆过程中探测器受到多源扰动条件下的轨迹跟踪控制方法,属于深空探测制导与控制技术领域。
背景技术
在行星着陆过程中,存在着天体引力场不确知、外界干扰等相关不确定性因素,且这些不确知因素带来的摄动将随着探测器高度的下降不断增大,这就需要探测器的跟踪控制算法对未知因素具有一定的处理能力。通常对于月球探测任务,由于对目标天体的长时间观测,事先获得了较为完备的天体特性信息,因此基于地面站的控制策略已成为对其进行接近操作过程中的常规模式。
然而,对于深空探测任务,目前并不具备对其进行长期观测的条件,在下降段探测任务中动力学环境复杂,可利用观测手段有限。尤其对于具有不规则引力场特性的行星着陆控制而言,一个很小的控制偏差会导致最终着陆器与目标着陆点的大范围偏离,因此整个阶段的控制精度要达到百米以内的可靠着陆性能就变得极具挑战性,同时控制机构的非线性特性进一步加剧了问题解决的难度。这就需要具有较强自主性和鲁棒性的跟踪控制方法,以实现行星表面高精度的着陆任务。
在先技术[1]中(参见J.Zhang,L.Y.Yang,G.Z.Shen.Modeling andAttitude Control of Aircraft with Center of Gravity Variations.IEEEAerospace Conference Proceedings.2009:3108-3118.),着重针对系统模型中的不确知参数特性,在控制律的设计中引入对模型参数的估计,并给出适当的参数估计调整律,使得控制系统得到稳定。该方法在一定程度上缓解了天体引力场模型参数不确定性对着陆轨迹跟踪的不利影响,然而该方法所做假设与简化条件过多,且并未考虑控制机构的开启关闭时间约束,很难得到广泛应用。
发明内容
本发明的目的在于针对由于行星表面环境模型不精确、参数时变,以及形状不规则、质量分布不均匀等复杂多变特性,给行星着陆控制系统带来的影响,提出一种基于模型不确定界的行星探测器着陆轨迹跟踪鲁棒控制方法。
本方法考虑目标天体引力模型与外界干扰等不确定因素对轨迹跟踪的影响,引入模型不确定界的设计理念,通过设计鲁棒反馈控制律,并兼顾控制机构的非线性特性,确保实际状态轨迹与前馈制导设计轨迹的偏离量不超过期望的界,从而优化多源扰动下的行星探测器着陆轨迹跟踪性能。
本发明的技术方案具体包括如下步骤:
步骤1,在目标天体质心固连坐标系下,建立探测器误差动力学模型:
其中,δx=χ-x为探测器状态误差,χ为真实的探测器状态,x为探测器标称状态,为探测器的反馈控制输入,φ为真实的引力加速度,g为对真实引力加速度φ的建模,d为外界扰动加速度,ψ=φ(Crχ)-φ(Crx)与为天体引力模型误差,A为系统矩阵,B为输入矩阵,Cr为状态系数矩阵。
所述目标天体质心固连坐标系的原点在小行星的质量中心,z轴沿小行星最大惯量轴方向,x轴沿历元时刻小行星最小惯量轴所指方向,y轴与x轴、z轴之间满足右手法则。
步骤2,考虑所需的推力器关闭时间,反馈策略的设计基于跟踪误差的离散动力学模型,设定反馈控制输入在给定时间间隔内与前馈控制同时作用。对误差动力学模型离散化,得到
其中,k为离散时间点,Ek为离散化后的扰动输入矩阵。
步骤3,设定误差约束集为一个满足的多胞体:
其中,δxmax为允许的最大跟踪误差,p表示多胞体的顶点,q为总的顶点数,ap为根据任务要求设定的状态约束系数矩阵。
步骤4,将天体引力模型误差ψ、与外界扰动加速度d描述为如下形式:
其中,Ω={Θ|ΘTΘ≤γ2I,γ>0},Θ是集合Ω中的变量;γ、σ与δ为设计任务给定的关于允许模型误差界与外界扰动界的指标参数。
步骤5,设计反馈控制律δuk,使探测器飞行状态控制误差δxk在步骤3得到的误差约束集内,且满足控制机构非线性特性||δuk||≤δUmax(δUmax为由设计任务给定反馈控制输入的约束界)。δuk的具体形式为:
δuk=LQ-1δxk
其中,矩阵Q=QT>0,矩阵L和Q满足:
α、ζ为标量,且1>α>0,ζ>0。
步骤6,将步骤5得到的控制律δuk=LQ-1δxk作为反馈控制输入,输入步骤1建立的探测器误差动力学模型,从而实现行星探测器着陆轨迹跟踪鲁棒控制。
有益效果
本发明的控制方法有效地克服了多源扰动和系统不确定性对探测器着陆飞行轨迹的影响,提高了探测器系统的控制性能,使行星探测器着陆飞行控制精度得到了保障。
附图说明
图1为本发明方法的流程图;
图2为具体实施方式中的探测器飞行状态时间历程;其中(a)为探测器位置三轴分量,(b)探测器速度三轴分量;
图3为具体实施方式中的探测器三轴反馈控制输入曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实施例对本发明内容作进一步说明。
①在目标天体质心固连坐标系下建立探测器着陆动力学标称模型:
其中,为探测器在目标天体质心固连坐标系下的矢径,为目标天体自旋角速率,u为探测器前馈控制输入。
设真实系统动力学描述为:
其中, Cr=[0 I],为向量叉乘ω×(·)的矩阵表示,uc为实际的探测器控制输入。
将式(8)与式(9)相减可得误差动力学模型为
②考虑所需的推力器关闭时间Δt,对误差动力学模型离散化,得到
其中,
③设计反馈控制律δuk=Kδxk(K为状态反馈控制增益矩阵),保证真实轨迹在标称轨迹周围的不变椭球εP={δx|δxTPδx≤1,P>0}内,且 为误差约束集。
具体过程为:设定天体引力与外界扰动的模型误差满足:
则在反馈控制律的作用下,误差状态可达域为:
存在Lyapunov函数满足
则椭球
根据式(14)得到如下不等式
应用S方法对上式进行处理,则存在τi≥0(i=1,…,4)使
对上式左右两边同时取负,有
令r1=1-α,τ2=β, 将式(17)展开为
从而得到
由Schur补引理,将式(19)等价表示为线性矩阵不等式形式:
令K=LP,Q=P-1,ζ=1β,并用矩阵diag([Q ζI I I I I])分别左右乘式(20),得到式(5)。将状态约束与控制约束||δuk||≤δUmax分别描述为线性不等式约束,如式(6)与式(7)所示。则可得到,若存在矩阵Q=QT>0,L和标量1>α>0,ζ>0,使由式(5)-式(7)构成的线性矩阵不等式组成立,则有控制律δuk=LQ-1δxk使探测器飞行状态控制误差且满足控制机构非线性特性||δuk||≤δUmax。
以小行星Eros着陆接近段任务为例,对上述控制方法进行验证分析。目标天体自旋角速度为1639.4°/day,名义半径为16km,引力常数为4.4621×105m3/s2,在仿真中采用其四阶引力势函数模型。探测器初始状态位置为[8950 20 -50]m,初始速度为[1.5 2 0]m/s,期望末端位置为[8450 0 0]m,期望末端速度为[0 0 0]ms,期望着陆时间为300s。
当标称系统模型与真实系统存在1%~5%的不确定性时,为确保探测器真实轨迹与标称轨迹的各轴位置偏差不超过3m,各轴速度偏差不超过0.2m/s,设定反馈控制误差状态约束集为设定反馈控制输入幅值上限为δUmax=20N。
考虑模型不确定性的存在,在行星着陆过程中加入反馈控制。图2中(a)和(b)分别为探测器三轴位置与速度分量,由图可见,探测器真实状态与标称轨迹基本相符,其位置偏差在3m范围内,末速度偏差在0.05m/s左右;进一步观察控制推力情况,由图3可见,探测器进行了9次反馈控制,每次反馈控制作用时间为10s,最大控制输入幅值不超过12N,可见其满足推力器在有限开启时间内的推力幅值上限约束,保证了跟踪控制的整体设计要求。
Claims (2)
1.基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法,其特征在于:包括以下步骤:
步骤1,在目标天体质心固连坐标系下,建立探测器误差动力学模型:
其中,δx=χ-x为探测器状态误差,χ为真实的探测器状态,x为探测器标称状态,为探测器的反馈控制输入,φ为真实的引力加速度,g为对真实引力加速度φ的建模,d为外界扰动加速度,ψ=φ(Crχ)-φ(Crx)与为天体引力模型误差,A为系统矩阵,B为输入矩阵,Cr为状态系数矩阵;
步骤2,设定反馈控制输入在给定时间间隔内与前馈控制同时作用;对误差动力学模型离散化,得到
其中,k为离散时间点,Ek为离散化后的扰动输入矩阵;
步骤3,设定误差约束集为一个满足的多胞体:
其中,δxmax为允许的最大跟踪误差,p表示多胞体的顶点,q为总的顶点数,ap为状态约束系数矩阵;
步骤4,将天体引力模型误差ψ、与外界扰动加速度d描述为如下形式:
其中,Ω={Θ|ΘTΘ≤γ2I,γ>0},γ、σ与δ为设计任务给定的关于允许模型误差界与外界扰动界的指标参数;
步骤5,设计反馈控制律δuk,使探测器飞行状态控制误差δxk在步骤3得到的误差约束集内,且满足控制机构非线性特性||δuk||≤δUmax,δUmax为由设计任务给定反馈控制输入的约束界;δuk的具体形式为:
δuk=LQ-1δxk
其中,矩阵Q=QT>0,矩阵L和Q满足:
α、ζ为标量,且1>α>0,ζ>0;
步骤6,将步骤5得到的控制律δuk=LQ-1δxk作为反馈控制输入,输入步骤1建立的探测器误差动力学模型,实现行星探测器着陆轨迹跟踪鲁棒控制。
2.根据权利要求1所述的基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法,其特征在于:所述目标天体质心固连坐标系的原点在小行星的质量中心,z轴沿小行星最大惯量轴方向,x轴沿历元时刻小行星最小惯量轴所指方向,y轴与x轴、z轴之间满足右手法则。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210499072.6A CN102968124B (zh) | 2012-11-29 | 2012-11-29 | 基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210499072.6A CN102968124B (zh) | 2012-11-29 | 2012-11-29 | 基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102968124A CN102968124A (zh) | 2013-03-13 |
CN102968124B true CN102968124B (zh) | 2015-04-15 |
Family
ID=47798313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210499072.6A Expired - Fee Related CN102968124B (zh) | 2012-11-29 | 2012-11-29 | 基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102968124B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103453907B (zh) * | 2013-09-06 | 2016-01-13 | 北京理工大学 | 基于分层大气模型的行星进入段导航滤波方法 |
CN104494845B (zh) * | 2014-12-03 | 2016-04-13 | 长春工业大学 | 无模型干扰在线估计的小天体探测器着陆控制方法 |
CN105159315B (zh) * | 2015-09-01 | 2017-10-31 | 西北工业大学 | 面向空间绳系机器人的目标逼近过程鲁棒规划控制方法 |
CN105447228B (zh) * | 2015-11-11 | 2018-10-19 | 中国人民解放军国防科学技术大学 | 一种地球扰动引力赋值模型适用性能评估方法 |
CN107168056B (zh) * | 2017-05-23 | 2019-10-11 | 浙江大学 | 一种自适应的月球软着陆轨道快速优化控制器 |
CN107203133B (zh) * | 2017-05-23 | 2019-10-11 | 浙江大学 | 一种智能的月球软着陆轨道控制器 |
CN107202584B (zh) * | 2017-07-06 | 2020-02-14 | 北京理工大学 | 一种行星精确着陆抗扰制导方法 |
CN108562293B (zh) * | 2018-03-13 | 2021-05-28 | 北京理工大学 | 基于干扰观测器的行星着陆有限时间控制方法 |
CN109211246B (zh) * | 2018-11-07 | 2020-08-21 | 北京理工大学 | 不确定环境下行星着陆轨迹规划方法 |
CN109976156B (zh) * | 2019-03-13 | 2021-08-06 | 南京航空航天大学 | 固定翼无人机栖落机动轨迹的建模与预测控制方法 |
CN117631547B (zh) * | 2024-01-26 | 2024-04-26 | 哈尔滨工业大学 | 一种小天体不规则弱引力场下的四足机器人着陆控制方法 |
-
2012
- 2012-11-29 CN CN201210499072.6A patent/CN102968124B/zh not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
Hazard Detection and Avoidance for lanetary Landing Based on Lyapunov Control Method;Shengying ZHU等;《Proceedings of the 10th World Congress on Intelligent Control and Automation》;20120708;第2822-2826页 * |
Reaserch on Optical/Inertial Integrated Navigation for Soft Landing on Asteroids Using Mixed H2/H∞ Filter;Ai Gao等;《International Conference on Ontelligent Control and Information Processing》;20100815;第20-23页 * |
实用飞行器软着陆控制律设计;马克茂等;《导弹与航天运载技术》;20011231(第2期);第39-43页 * |
着陆小行星的滑模变结构控制;李爽等;《宇航学报》;20051130;第26卷(第6期);第808-812页 * |
附着小行星的视线制导规律;崔祜涛等;《空间科学学报》;20020731;第22卷(第3期);第256-260页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102968124A (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102968124B (zh) | 基于模型不确定界的行星着陆轨迹跟踪鲁棒控制方法 | |
Hanover et al. | Performance, precision, and payloads: Adaptive nonlinear mpc for quadrotors | |
CN108287476B (zh) | 基于高阶滑模控制和扰动观测器的空间翻滚非合作目标自主交会制导方法 | |
CN104457705B (zh) | 基于天基自主光学观测的深空目标天体初定轨方法 | |
Gao et al. | Robust fault tolerant tracking control design for a linearized hypersonic vehicle with sensor fault | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN106778012A (zh) | 一种小天体附着探测下降轨迹优化方法 | |
CN110376882A (zh) | 基于有限时间扩张状态观测器的预定性能控制方法 | |
CN102981507A (zh) | 一种软着陆自主障碍规避常推力器控制方法 | |
CN105629732A (zh) | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 | |
CN104656447A (zh) | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 | |
CN104215244A (zh) | 基于发射惯性坐标系的空天飞行器组合导航鲁棒滤波方法 | |
Allibert et al. | Velocity aided attitude estimation for aerial robotic vehicles using latent rotation scaling | |
CN102582850B (zh) | 提高卫星磁控精度的方法 | |
CN105138808A (zh) | 基于摄动理论的滑翔弹道误差传播分析方法 | |
CN104192322A (zh) | 一种行星动力下降段轨迹在线生成的抗干扰制导控制方法 | |
CN103123487A (zh) | 一种航天器姿态确定方法 | |
CN109613827A (zh) | 一种相对速度未知的平动点轨道交会控制方法 | |
Kang et al. | Finite-memory-structured online training algorithm for system identification of unmanned aerial vehicles with neural networks | |
CN105527974A (zh) | 一种缺失径向控制的欠驱动航天器悬停渐近控制方法 | |
CN103412573A (zh) | 基于级联方程的椭圆轨道航天器相对位置退步控制方法 | |
CN107678284A (zh) | 高速飞行器的鲁棒补偿控制方法及高速飞行器 | |
CN102998975B (zh) | 一种欠驱动航天器角速度稳定的鲁棒控制方法 | |
CN102636817B (zh) | 一种纯引力轨道万有引力干扰抑制的质量补偿设计方法 | |
Zhang et al. | Robust trajectory tracking guidance for low L/D lunar return vehicles using command filtered backstepping approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150415 Termination date: 20171129 |
|
CF01 | Termination of patent right due to non-payment of annual fee |