CN110018637B - 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 - Google Patents
一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 Download PDFInfo
- Publication number
- CN110018637B CN110018637B CN201811611233.XA CN201811611233A CN110018637B CN 110018637 B CN110018637 B CN 110018637B CN 201811611233 A CN201811611233 A CN 201811611233A CN 110018637 B CN110018637 B CN 110018637B
- Authority
- CN
- China
- Prior art keywords
- attitude
- control
- tracking
- spacecraft
- alpha
- 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
-
- 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/042—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 parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种考虑完成时间约束的航天器姿态跟踪保性能控制方法,对刚性航天器的姿态跟踪控制问题,提出一种不需要预先知道系统信息(主要是惯量矩阵)的无模型保性能姿态控制技术,实现对姿态控制系统的鲁棒跟踪控制,并能保证技术指标的实现,还能兼顾任务的时间约束,保证系统状态能够在给定时间内完成对期望信号的跟踪。有益效果:能够保证系统在预设时间之前完成收敛,该控制方法具有无模型特点,即整个控制回路中均不需要模型具体参数(主要是姿态惯量矩阵)的信息,降低了控制器设计对系统的要求,姿态控制系统的技术指标(如稳态误差、瞬态收敛速度等)能够提前设计,并通过控制器中的边界函数进行保证,具有强鲁棒性的特点。
Description
技术领域
本发明属于航天器鲁棒控制领域,涉及一种考虑完成时间约束的航天器姿态跟踪保性能控制方法。
背景技术
目前已有很多基于模型信息的航天器姿态控制技术,如滑模控制技术,文献:CongB,Chen Z,Liu X.Improved adaptive sliding mode control for rigid spacecraftattitude tracking[J].Journal of Aerospace Engineering,2014, 27(4):04014004.;最优控制技术,文献:Pukdeboon C.Optimal output feedback controllers forspacecraft attitude tracking[J].Asian Journal of Control, 2012,15(5):1284-1294.;退步控制技术,文献:Ali I,Radice G,Kim J.Backstepping control design withactuator torque bound for spacecraft attitude maneuver[J]. Journal ofGuidance Control Dynamics,2012,33(1):254-259.等。这些技术在数 值仿真中取得了不错的效果,其中的一些甚至已经应用在合作航天器的在轨姿态控制 中。然而,随着空间技术的发展,新型航天任务向姿态控制系统提出了更高的要求。 大多数基于模型的控制方法需要航天器惯量矩阵的先验信息。但在一些新型任务,如 捕获非合作目标后的组合体姿态控制中,这些信息是难以提前获取的。
基于此,有学者开始考虑如何在未知模型信息的前提下,进行控制系统设计,即无模型控制技术。近期有学者提出了一种保性能控制技术,文献:Bechlioulis C P,Rovithakis G A.A low-complexity global approximation-free control scheme withprescribed performance for unknown pure feedback systems[J].Automatica,2014,50(4):1217-1226.和Theodorakopoulos A,Rovithakis G A.Low-Complexity PrescribedPerformance Control of Uncertain MIMO Feedback Linearizable Systems[J].IEEETransactions on Automatic Control,2016,61(7):1946-1952., 该技术能够在未知系统模型信息的前提下保证系统的收敛轨迹在预先给定的技术指标 以内,因此得到了极大的关注。在航天领域中,该技术也得到了相应的应用。文献Zhou Z G,Zhang Y A,Shi X N,etal.Robust attitude tracking for rigid spacecraft with prescribed transientperformance[J].International Journal of Control, 2017,90(11):2471-2479.中,作者将保性能控制技术引入了姿态稳定控制中,并严 格证明了预先设计的性能指标的实现。文献Luo J,Yin Z,Wei C,et al. Low-complexity prescribed performance controlfor spacecraft attitude stabilization and tracking[J].Aerospace Science andTechnology,2018,74: 173-183.中,作者考虑了系统参数(航天器的惯量矩阵)完全未知的情形,提出了一 种无模型的保性能控制技术。
然而,上述技术仍然存在技术缺陷。传统的保性能控制技术在理论上只能保证实际状态量在无穷时间处完成对姿态指令的跟踪,虽然在实际操作中能够在有限时间内 完成跟踪,但是该时间无法进行先验获取和设计。然而对于实际的航天控制任务,往 往要求姿态跟踪任务在有限时间完成,当超过时限未能完成任务时,可能会造成任务 的延误甚至失败。因此有必要提出一种考虑完成时间约束的航天器姿态跟踪控制方法。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种考虑完成时间约束的航天器姿态跟踪保性能控制方法,针对无法预先获取系统模型信息的姿态控制问题,考虑姿态控 制任务完成的时间约束,本发明方法,使得系统的收敛轨迹能够被约束在预先给定的 技术指标以内,实现保性能控制,且系统能够在给定时间约束内进入稳态,保证后续 任务的顺利完成。
技术方案
一种考虑完成时间约束的航天器姿态跟踪保性能控制方法,其特征在于步骤如下:
步骤1、姿态跟踪系统建模和组合状态量构造:
其中:σe=[σe1,σe2,σe3]T∈R3为姿态误差修正罗德里格斯参数(MRP),ωe=[ωe1,ωe2,ωe3]T∈R3为角速度误差, ωr=[ωr1,ωr2,ωr3]T∈R3为期望角速度,对于任意三维向量γ=[γ1 γ2 γ3]T,叉乘算子[γ×]定义 为:
航天器姿态跟踪控制系统的动力学方程为:
其中:ω=[ω1,ω2,ω3]T∈R3为航天器的真实角速度,J∈R3×3为航天器的惯量矩阵,在 本技术中设置为未知;τ∈R3和τd∈R3分别为控制力矩和未知上界的外部干扰
其中:λ>0为设计参数;
步骤2、技术指标边界函数设计:
其中:γ∈(0,1),α0和αT分别为边界函数的初值和终值,T为新型边界函数的收敛时间;
以边界函数,分别对姿态信息σei(t)和组合状态量信息xi(t)施加如下所示的边界函 数约束:
其中:δ∈(0,1)为抑制超调量的常数,αei(t)和αxi(t)的初值和终值分别定义为:αe0,i,αeT,i,αx0,i,αxT,i
步骤3、鲁棒无模型保性能姿态跟踪控制器设计:
其中:k>0为可调控制增益,Re和Rx分别定义为:Re=diag(Re1,Re2,Re3), Rx=diag(Rx1,Rx2,Rx3),且
所述映射跟踪误差se(t)=[se1(t),se2(t),se3(t)]T和sx(t)=[sx1(t),sx2(t),sx3(t)]T为:
-αxi(t)-4λαei(t)<ωei(t)<αxi(t)+4λαei(t)
在时间t≥T后,系统的姿态信息σei(t),组合状态量信息xi(t)和角速度信息ωei(t)将 分别收敛到所示的稳定域内:
Uωi:=(-αxT,i-4λαeT,i,αxT,i+4λαeT,i)。
有益效果
本发明提出的一种考虑完成时间约束的航天器姿态跟踪保性能控制方法,对刚性航天器的姿态跟踪控制问题,提出一种不需要预先知道系统信息(主要是惯量矩阵) 的无模型保性能姿态控制技术,实现对姿态控制系统的鲁棒跟踪控制,并能保证技术 指标的实现,还能兼顾任务的时间约束,保证系统状态能够在给定时间内完成对期望 信号的跟踪。
针对惯量矩阵未知的航天器的无模型姿态跟踪控制问题,本发明具有以下三点有益效果:
1、用户可以预先设计系统的收敛时间,该技术能够保证系统在预设时间之前完成收敛,这对于一些需要在规定时间内完成的任务有重大的实际意义;
2、该控制方法具有无模型特点,即整个控制回路中均不需要模型具体参数(主要是姿态惯量矩阵)的信息,降低了控制器设计对系统的要求,具有很强的工程实用价 值和理论意义;
3、姿态控制系统的技术指标(如稳态误差、瞬态收敛速度等)能够提前设计,并 通过控制器中的边界函数进行保证,具有强鲁棒性的特点。
附图说明
图1:不同收敛时间下边界函数的收敛情况
图2-4为鲁棒无模型保性能控制器作用下系统状态量的仿真结果图;
图2:姿态跟踪MRP误差示意图
图3:姿态跟踪角速度误差示意图
图4:姿态跟踪组合状态量示意图
图5:为鲁棒无模型保性能控制器计算得到的控制力矩的仿真结果图;
图6:为设置不同稳定时间时,姿态跟踪MRP误差响应的仿真结果图。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本实施例:以刚性航天器的姿态跟踪问题作为实施对象,其中航天器的惯量矩阵设置为J=diag(40.0,42.5,50.2)kg·m2,初始的MRP设置为σ(0)=[0.6,-0.5,0.3]T,初期的期望初始 MRP设置为σr(0)=[0.02,-0.03,0.05]T,初始角速度设计为:ω(0)=[-0.001,0.002,-0.0009]Trad/s, 时变期望角速度设置为:
ωr=2×10-4×[0.03sin(πt/200),0.03sin(πt/300),0.03sin(πt/250)]Trad/s
外部干扰设置为:
τd=2×10-4×[sin(0.8t),cos(0.5t),cos(0.3t)]TN·m
边界函数αei(t)和αxi(t)(i=1,2,3)的参数选取为:T=25,γ=0.6,δ=0.8,αei,0=1.5, αei,T=0.01,αxi,0=1.5和αxi,T=0.01。
控制器参数设计为:k=2,λ=0.3。
具体过程:
步骤一:姿态跟踪系统建模和组合状态量构造
本发明针对的航天器姿态跟踪控制系统的运动学方程为:
其中:σe=[σe1,σe2,σe3]T∈R3为姿态误差修正罗德里格斯参数(MRP),ωe=[ωe1,ωe2,ωe3]T∈R3为角速度误差, ωr=[ωr1,ωr2,ωr3]T∈R3为期望角速度,对于任意三维向量γ=[γ1 γ2 γ3]T,叉乘算子[γ×]定义 为:
本发明针对的航天器姿态跟踪控制系统的动力学方程为:
其中:ω=[ω1,ω2,ω3]T∈R3为航天器的真实角速度,J∈R3×3为航天器的惯量矩阵,在本技 术中设置为未知;τ∈R3和τd∈R3分别为控制力矩和未知上界的外部干扰。
构造如下所示的组合状态量:
其中:λ>0为设计参数。
步骤二:新型技术指标边界函数设计
区别于传统的指数稳定的技术指标边界函数,本发明设计了如下所示的新型边界函数:
其中:γ∈(0,1),α0和αT分别为边界函数的初值和终值,T为新型边界函数的收敛时间。
本发明设计的新型边界函数(4),可由用户任意给定其收敛时间。通过设置不同的收敛时间,可以对系统变量的性能起到很好的约束作用。该边界函数可通过图1进行 直观阐述。图中,TPF代表传统的边界函数收敛情况,可以看出其必须在时间趋于无 穷时才能完成收敛。ARPF代表本发明设计的新型边界函数。该函数可由用户任意设计 收敛时间,具有很强的实际意义和实用价值。
基于式(4)所示的新型边界函数,分别对姿态信息σei(t)和组合状态量信息xi(t)施加 如下所示的边界函数约束:
其中:δ∈(0,1)为抑制超调量的常数,αei(t)和αxi(t)的定义与式(3)相同,且其初值和终值 分别定义为:αe0,i,αeT,i,αx0,i,αxT,i。
步骤三:鲁棒无模型保性能姿态跟踪控制器设计
进行无约束化处理,定义映射跟踪误差se(t)=[se1(t),se2(t),se3(t)]T和sx(t)=[sx1(t),sx2(t),sx3(t)]T为:
基于此,式(5)和(6)中的有约束控制问题就转化为对映射跟踪误差se(t)和sx(t)的无 约束控制问题。
基于上述转化,本发明设计了如下所示的鲁棒无模型保性能姿态跟踪控制器:
其中:k>0为可调控制增益,Re和Rx分别定义为:Re=diag(Re1,Re2,Re3),Rx=diag(Rx1,Rx2,Rx3), 且
系统的稳定时间和稳定域分析
针对式(1)和(2)所示的姿态跟踪控制系统,在式(9)所示的鲁棒无模型保性能姿态 跟踪控制器的作用下,将满足如下所示的性能约束:
-αxi(t)-4λαei(t)<ωei(t)<αxi(t)+4λαei(t) (14)
结合式(4)所示的新型边界函数的定义,可以得出:在时间t≥T后,系统的姿态 信息σei(t),组合状态量信息xi(t)和角速度信息ωei(t)将分别收敛到如下所示的稳定域 内:
Uωi:=(-αxT,i-4λαeT,i,αxT,i+4λαeT,i) (17)
本发明针对刚性航天器的姿态跟踪控制问题,提出一种不需要预先知道系统信息(主要是惯量矩阵)的无模型保性能姿态控制技术,实现对姿态控制系统的鲁棒跟踪 控制,并能保证技术指标的实现,还能兼顾任务的时间约束,保证系统状态能够在给 定时间内完成对期望信号的跟踪。
Claims (2)
1.一种考虑完成时间约束的航天器姿态跟踪保性能控制方法,其特征在于步骤如下:
步骤1、姿态跟踪系统建模和组合状态量构造:
其中:σ e = [ σ e1 , σ e2 , σ e3 ]T ∈ R3 为姿态误差修正罗德里格斯参数(MRP),ωe=[ωe1,ωe2,ωe3]T∈R3为角速度误差,ωr=[ωr1,ωr2,ωr3]T∈R3为期望角速度,对于任意三维向量γ=[γ1 γ2 γ3]T,叉乘算子[γ×]定义为:
航天器姿态跟踪控制系统的动力学方程为:
其中:ω=[ω1,ω2,ω3]T∈R3为航天器的真实角速度,J∈R3×3为航天器的惯量矩阵,在本技术中设置为未知;τ∈R3和τd∈R3分别为控制力矩和未知上界的外部干扰
其中:λ>0为设计参数;
步骤2、技术指标边界函数设计:
其中:γ∈(0,1),α0和αT分别为边界函数的初值和终值,T为新型边界函数的收敛时间;
以边界函数,分别对姿态信息σei(t)和组合状态量信息xi(t)施加如下所示的边界函数约束:
其中:δ∈(0,1)为抑制超调量的常数,αei(t)和αxi(t)的初值和终值分别定义为:αe0,i,αeT,i,αx0,i,αxT,i
步骤3、鲁棒无模型保性能姿态跟踪控制器设计:
其中:k>0为可调控制增益,Re和Rx分别定义为:Re=diag(Re1,Re2,Re3),Rx=diag(Rx1,Rx2,Rx3),且
se(t)=[se1(t),se2(t),se3(t)]T和sx(t)=[sx1(t),sx2(t),sx3(t)]T为映射跟踪误差,定义为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811611233.XA CN110018637B (zh) | 2018-12-27 | 2018-12-27 | 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811611233.XA CN110018637B (zh) | 2018-12-27 | 2018-12-27 | 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110018637A CN110018637A (zh) | 2019-07-16 |
CN110018637B true CN110018637B (zh) | 2021-08-13 |
Family
ID=67188677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811611233.XA Active CN110018637B (zh) | 2018-12-27 | 2018-12-27 | 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110018637B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111442786B (zh) * | 2020-05-21 | 2021-08-20 | 合肥工业大学 | 一种飞行器陀螺仪的零点漂移偏差和姿态估计方法 |
CN113325861B (zh) * | 2021-06-02 | 2023-03-24 | 上海海事大学 | 一种非奇异预设时间四旋翼无人机姿态跟踪控制方法 |
CN113485395B (zh) * | 2021-06-29 | 2024-03-15 | 北京控制工程研究所 | 一种误差约束下固定时间相对姿轨跟踪控制方法 |
CN113485394B (zh) * | 2021-06-29 | 2023-06-30 | 北京控制工程研究所 | 一种高精度固定时间收敛的相对姿态容错跟踪控制方法 |
CN114019848B (zh) * | 2021-10-14 | 2024-04-30 | 南京航空航天大学 | 一种输入受限的机械系统保性能自调节控制方法 |
CN114329943B (zh) * | 2021-12-23 | 2023-01-24 | 哈尔滨工业大学(深圳) | 基于姿态旋转矩阵的控制性能边界设计方法、装置及介质 |
CN115629550B (zh) * | 2022-12-22 | 2023-04-18 | 西北工业大学 | 一种服务型航天器的自适应姿态跟踪控制及参数辨识方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7219014B2 (en) * | 2004-04-29 | 2007-05-15 | The Boeing Company | Method and apparatus for real-time star exclusion from a database |
CN102645892B (zh) * | 2012-03-16 | 2014-11-05 | 浙江工业大学 | 交通网络化控制系统的控制方法 |
CN104460678B (zh) * | 2014-07-25 | 2017-02-22 | 北京机械设备研究所 | 一种基于Sigmoid型非线性滑模函数的航天器姿态控制方法 |
CN104656447A (zh) * | 2015-01-16 | 2015-05-27 | 西北工业大学 | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 |
CN104527994B (zh) * | 2015-01-21 | 2016-08-24 | 哈尔滨工业大学 | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 |
CN106681140B (zh) * | 2016-12-23 | 2020-10-16 | 西北工业大学 | 一种不确定大系统低复杂度鲁棒保性能分散控制方法 |
CN108181913B (zh) * | 2017-12-06 | 2019-05-03 | 北京航空航天大学 | 一种具有指定跟踪性能的航天器自适应容错姿态跟踪控制方法 |
-
2018
- 2018-12-27 CN CN201811611233.XA patent/CN110018637B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110018637A (zh) | 2019-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110018637B (zh) | 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 | |
Wang et al. | Disturbance observer-based finite-time control design for a quadrotor UAV with external disturbance | |
Sun et al. | Fixed-time sliding mode disturbance observer-based nonsmooth backstepping control for hypersonic vehicles | |
Do et al. | Underactuated ship global tracking under relaxed conditions | |
Sun et al. | Fault-tolerant track control of hypersonic vehicle based on fast terminal sliding mode | |
CN111324142B (zh) | 一种导弹驾驶仪扰动补偿控制方法 | |
CN108594837A (zh) | 基于pd-smc和rise的无模型四旋翼无人机轨迹跟踪控制器及方法 | |
Wu et al. | Quantized fixed-time fault-tolerant attitude control for hypersonic reentry vehicles | |
Durmaz et al. | Sliding mode control for non-linear systems with adaptive sliding surfaces | |
Zhang et al. | Adaptive nonsingular terminal sliding mode control design for near space hypersonic vehicles | |
CN109634291B (zh) | 一种基于改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 | |
Jiao et al. | Adaptive mode switching of hypersonic morphing aircraft based on type-2 TSK fuzzy sliding mode control | |
CN109507890A (zh) | 一种基于eso的无人机动态逆广义预测控制器 | |
Van Oort et al. | Full-envelope modular adaptive control of a fighter aircraft using orthogonal least squares | |
Li et al. | Finite-time control for quadrotor based on composite barrier Lyapunov function with system state constraints and actuator faults | |
Qiao et al. | Robust H2 optimal depth control of an autonomous underwater vehicle with output disturbances and time delay | |
Ding et al. | Global smooth sliding mode controller for flexible air-breathing hypersonic vehicle with actuator faults | |
Huang et al. | Adaptive backstepping control for autonomous shipboard landing of a quadrotor with input saturation | |
Jiang et al. | Heading control of unmanned surface vehicle with variable output constraint model-free adaptive control algorithm | |
CN116873226A (zh) | 一种航天器姿态控制及rcs优化分配方法 | |
An et al. | Adaptive terminal sliding mode control for reentry vehicle based on nonlinear disturbance observer | |
Cao et al. | Robust multiple model predictive control for ascent trajectory tracking of aerospace vehicles | |
CN104460678B (zh) | 一种基于Sigmoid型非线性滑模函数的航天器姿态控制方法 | |
CN103809446B (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
Sun et al. | Tracking control of hypersonic vehicles with input saturation based on fast terminal sliding mode |
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 |