CN108445748B - 一种基于事件触发机制的自适应航天器姿态跟踪控制方法 - Google Patents
一种基于事件触发机制的自适应航天器姿态跟踪控制方法 Download PDFInfo
- Publication number
- CN108445748B CN108445748B CN201810107715.5A CN201810107715A CN108445748B CN 108445748 B CN108445748 B CN 108445748B CN 201810107715 A CN201810107715 A CN 201810107715A CN 108445748 B CN108445748 B CN 108445748B
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- control
- event trigger
- matrix
- trigger mechanism
- 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
技术领域
本发明涉及一种基于事件触发机制的自适应航天器姿态跟踪控制方法,主要应用于模块化航天器的姿态跟踪控制,属于航天器控制技术领域。
背景技术
近年来,模块化航天器受到重视,模块化航天器采用模块化开放式网络架构,利用即插即用系统的基本思想,将航天器的各个组件模块化而进行组装,可以降低航天器成本,加快组装速度,利于有缺陷或过时的部件轻松拆卸和更换。模块化航天器的核心技术是无线数据通信和无线电力传输,它与笨重、庞大、不灵活的传统电缆互连不同,在此技术下,航天器的所有功能部件都相互独立,并通过低成本无线网络连接。但是,负责执行器模块和控制模块之间的数据传输的无线网络的带宽和计算能力是有限的。因此,如何设计航天器的控制方案,以减少通信量而不影响稳定性和控制性能是非常重要的。
传统控制方法都是在连续时间框架内开发的,控制算法在一个固定的、周期的和足够小的采样周期内执行,有很好的理论和分析的方便性。但是,在这种方式中,无论是否必要,在每一个采样时间内,执行器和控制器之间都要进行通信以及更新控制算法,这将造成网络资源的浪费和执行机构的磨损。显然,这种方法在网络资源有限的条件下是不合适的。而事件触发机制可以依靠设定控制更新条件,来减小控制更新频率,大大减少通信量。
发明内容
为了解决上述已有技术存在的不足,本发明提供一种基于事件触发机制的自适应航天器姿态跟踪控制方法,包括以下步骤:
(1)建立含有外部干扰的航天器动力学模型与相对姿态跟踪误差模型;
(2)定义控制过程中应用的事件触发机制,基于事件触发机制构建传输控制力矩信号与控制器设计的控制力矩信号之间的关系式;
(3)基于步骤(1)建立的航天器动力学模型与相对姿态跟踪误差模型、步骤(2)的基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式,利用反步法对航天器姿态进行自适应控制;
具体的:
步骤(1)的航天器动力学模型与相对姿态跟踪误差模型为:
其中,J为航天器的转动惯量矩阵,且是3×3的正定对称矩阵;w=[w1,w2,w3]T为航天器在本体坐标系下相对于惯性坐标系的姿态角速度,w1,w2,w3分别为姿态角速度在本体坐标系的x轴、y轴和z轴上的分量;we=[we1,we2,we3]T为航天器在本体坐标系上相对于期望坐标系的姿态角速度,为we相对时间的导数;wd=[wd1,wd2,wd3]T为航天器在期望坐标系相对于惯性坐标系的姿态角速度,为wd相对时间的导数;S(w)、S(we)是斜对称矩阵,其形式分别为t为时间;τ(t)=[τ1(t),τ2(t),τ3(t)]T为航天器的执行机构实际产生的控制力矩信号,因为不考虑所述航天器的执行机构受到的干扰,所以等于基于所述事件触发机制构建的传输控制力矩信号;d为航天器受到的空间环境中的干扰力矩信号,其值未知但是有界;q0e、qe=[qe1,qe2,qe3]T分别为航天器由期望坐标系相对本体坐标系的相对姿态跟踪误差四元数的标量和向量,分别是q0e、qe相对时间的导数,是的转置矩阵,且满足C为期望姿态到航天器实际姿态的转换矩阵,是3×3的矩阵,且||C||为C的二范数,且||C||=1,为C的导数,且S(qe)是斜对称矩阵,其形式为I是3×3的单位矩阵;
步骤(2)的基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式为:
其中,τi(t)为基于事件触发机制构建的传输控制力矩信号τ(t)的分量;ki(t)和为的基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式的系数,均为时变标量;ui(t)是的控制器设计的控制力矩信号u(t)=[u1(t),u2(t),u3(t)]T的分量。
步骤(3)为基于步骤(1)建立的航天器动力学模型与相对姿态跟踪误差模型、步骤(2)基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式,利用反步法对航天器进行自适应控制设计:
基于反步法,首先引入两个新变量z1=qe,z2=we-α1,设计虚拟控制量为:
α1=-c1z1
其中,c1>0为虚拟控制器的增益;
基于虚拟控制量,设计自适应控制器和自适应更新律:
其中,为书写简便记 为ud的转置矩阵;为书写简便记||z2||为新变量z2的二范数,为z2的转置矩阵;为书写简便记且f是3×1维矩阵,表示为f=[f1,f2,f3]T;L为线性算子,将f从3×1维空间映射到3×6维空间,表示为L(f)T为L(f)的转置矩阵;L(w)表示为L(w)T为L(w)的转置矩阵;S(w)T是S(w)的转置矩阵;Γ1是控制器增益矩阵,为6×6的正对称矩阵;α1为虚拟控制量,为α1相对时间的导数;为惯性矩阵J的线性回归量,为的估计值,为估计值的自适应律;与步骤(2)的基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式相结合,构建时变控制增益矩阵K(t)=diag{k1(t),k2(t),k3(t)},则基于所述事件触发机制构建的传输控制力矩信号即航天器受到执行机构实际的控制力矩信号与控制器设计的控制力矩信号之间的关系式可表示为为方便分析,构建航天器受到的外界干扰和为对时变控制增益矩阵K(t)进行界估计,则存在常数km使在t≥0的任何时间内都有K(t)的最小特征值λmin(K(t))≥km>0;表示对时变控制增益矩阵K(t)的上界估计,是l的估计值,为估计值的自适应律;表示外界干扰和Td上界估计的范数的与时变系数K(t)的上界估计的乘积,是p的估计值,为估计值的自适应律;ε(t)是有界的光滑辅助信号,ε(t)>0,满足取δ1和δ2是辅助信号的常数,满足δ1>0,δ2>0;γ1,γ2是控制器增益参数,是正的常数。
优选地,姿态角速度w的初始值为w(0)=[-0.2,0.15,-0.3]T,航天器受到的空间环境中的干扰力矩信号d的初始值为d=0.15×[sin(0.1πt),sin(0.2πt),sin(0.1πt)]T,期望的单位四元数的标量q0d的初始值为qd(0)=[0.2,-0.2,0.3]T,航天器姿态的单位四元数的标量和向量q和q0的初始值为航天器在期望坐标系上相对于惯性坐标系的姿态角速度矢量wd的初始值为wd=0.1×[sin(t/40),-sin(t/50),2sin(t/30)]T。
优选地,选取的有界的光滑辅助信号为ε(t)=0.2×e-0.01t。
本发明的优点在于:
(1)与传统控制方法相比,将事件触发机制与自适应控制相结合,减小航天器执行器和控制器之间的通信以及更新控制算法次数,大大减少通信量和控制执行频率,并能够保证控制的精准性。
(2)引入了界估计和辅助信号处理航天器受未知惯性参数和外部干扰的影响,使系统渐进稳定。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明的一种基于事件触发机制的自适应航天器姿态跟踪控制方法的流程图;
图2为本发明的一种基于事件触发机制的自适应航天器姿态跟踪控制方法的系统框图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
下面将结合附图和实施例对本发明作进一步的详细说明。
如图1所示,本发明的一种基于事件触发机制的自适应航天器姿态跟踪控制方法步骤为:
(1)建立含有外部干扰的航天器动力学模型与相对姿态跟踪误差模型;
(2)定义控制过程中应用的事件触发机制,基于所述事件触发机制构建传输控制力矩信号与控制器设计的控制力矩信号之间的关系式;
(3)基于步骤(1)建立的航天器动力学模型与相对姿态跟踪误差模型、步骤(2)所述的基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式,利用反步法对航天器姿态进行自适应控制。
图2为本发明的一种基于事件触发机制的自适应航天器姿态跟踪控制方法的系统框图,如图所示,自适应控制器设计出的控制力矩信号经过事件触发机制处理,作为传输控制力矩信号传输给通信网络,通信网络再将传输控制力矩信号传输给执行机构,执行机构输出控制力矩信号作用到航天器上。
具体的实施步骤为:
步骤(1):用四元数法建立航天器的动力学模型与姿态模型:
其中,J为航天器的转动惯量矩阵,且是3×3的正定对称矩阵;w=[w1,w2,w3]T是航天器在本体坐标系下相对于惯性坐标系的姿态角速度,w1,w2,w3分别为姿态角速度在本体坐标系的x轴、y轴和z轴上的分量,姿态角速度的初始值可任意选取,优选地,姿态角速度的初始值为w(0)=[-0.2,0.15,-0.3]T,为w相对时间的导数;t为时间;τ(t)=[τ1(t),τ2(t),τ3(t)]T为航天器的执行机构实际产生的控制力矩信号,因为不考虑执行机构受到的干扰,所以等于基于所述事件触发机制构建的传输控制力矩信号;d为航天器受到的空间环境中的干扰力矩信号,其值未知但是有界,初始值可任意选取,优选地,航天器受到的空间环境中的干扰力矩信号初始值为d=0.15×[sin(0.1πt),sin(0.2πt),sin(0.1πt)]T;S(w)为斜对称矩阵,其形式为q0和q为航天器姿态的单位四元数的标量和向量,表示卫星在本体坐标系下相对于惯性坐标系的旋转姿态,q0与绕欧拉轴旋转的角度有关,α表示绕着欧拉轴转过的一个角度,和分别为q0和q相对时间的导数,q=[q1,q2,q3]T与欧拉轴方向有关,ex,ey,ez代表欧拉轴三个方向上的旋转轴,满足在满足前述条件的情况下q和q0的初始值可任意选取,优选地,q和q0初始值选为S(q)为斜对称矩阵,其形式为I为3×3的单位矩阵。
为研究航天器姿态跟踪问题,引入从期望坐标系到本体坐标系的航天器的相对姿态误差q0e和qe=[qe1,qe2,qe3]T分别是航天器由期望坐标系相对本体坐标系的相对姿态跟踪误差四元数的标量和向量,q0e=q0q0d-qT qd,qe=-q0qd+qq0d+S(q)qd,和qT分别为qe和q的转置矩阵;q0d和qd为期望的单位四元数的标量与向量,且满足q2 0d+q2 d=1,在满足前述条件的情况下q0d的初始值可任意选取,优选地,q0d初始值为qd(0)=[0.2,-0.2,0.3]T,wd=[wd1,wd2,wd3]T为航天器在期望坐标系上相对于惯性坐标系的姿态角速度矢量,初始值可任意选取,优选地,航天器在期望坐标系上相对于惯性坐标系的姿态角速度矢量的初始值为wd=0.1×[sin(t/40),-sin(t/50),2sin(t/30)]T。
引入航天器在本体坐标系上相对于期望坐标系的姿态角速度:
we=w-Cwd
其中,we=[we1,we2,we3]T为航天器在本体坐标系上相对于期望坐标系的姿态角速度;wd=[wd1,wd2,wd3]T为航天器在本体坐标系相对于惯性坐标系的期望角速度;C是期望姿态到航天器实际姿态的转换矩阵,是3×3的矩阵,||C||为C的二范数,且||C||=1,为C的导数,且S(qe)是斜对称矩阵,其形式为推导出航天器动力学模型与相对姿态跟踪误差模型如下:
步骤(2):定义控制中应用的事件触发机制,基于事件触发机制构建传输控制力矩信号与控制器设计的控制力矩信号之间的关系式:
首先,定义触发事件:
其中,τi(t)是基于事件触发机制构建的传输控制力矩信号τ的分量,|τi(t)|是τi(t)的绝对值;ui(t)是控制器设计的控制力矩信号u(t)=[u1(t),u2(t),u3(t)]T的分量;ei(t)=ui(t)-τi(t)是控制误差,|ei(t)|是ei(t)的绝对值;m1,m2,δ是定义事件触发阈值的设计参数,是正数,且δ∈(0,1),D是基于事件触发机制构建的传输控制力矩信号阈值,是正数,根据控制效果分别取值为m1=0.01,m2=1.5,D=4和δ=0.3;tk(k=1,2,3,…)是基于事件触发机制的传输控制力矩信号的更新时间,tk+1为下一更新时间。
事件触发机制的触发原理为:当满足|τi(t)|≤D,取满足|ei(t)|≥δ|τi(t)|+m1的最小时间t为下一更新时间tk+1;当满足|τi(t)|>D,取满足|ei(t)|≥m2的最小时间t为下一更新时间tk+1。当满足触发条件,时间将被记为下一更新时间tk+1,同时在时刻tk+1控制器设计的控制力矩信号ui(tk+1)将会传递给航天器的执行器模块。在时间间隔[tk,tk+1)中,基于事件触发机制构建的传输控制力矩信号保持为上一更新时间控制器设计的控制力矩信号ui(tk)。
进一步地,在定义的触发事件中包括以下两种控制更新策略:
1)相对阈值策略。当|τi|≤D时,事件触发阈值与控制器设计的控制力矩信号成比例。
2)固定阈值策略。当|τi|>D时,事件触发阈值为常数。
将以上两种阈值策略相结合得到切换阈值策略,具有以下优点:可以根据控制需求设定基于事件触发机制构建的传输控制力矩信号的阈值,当|τi|≤D,相对阈值策略可以保证事件触发阈值跟随控制器设计的控制力矩信号,能够进行更精确的控制,达到更好的控制性能。当|τi|>D,固定阈值策略能够防止基于事件触发机制构建的传输控制力矩信号突然以较大的脉冲跳动而导致的控制性能的恶化。
在所定义的触发事件中所有与事件相关的参数都应该基于保证闭环系统解的存在性和唯一性的指导原则,由于系统未知,假设这些参数是未知的。同时可知,减小δ,m1和m2,事件触发阈值也相应减小,得到更精确的控制,但过小会增加网络通信量。在时间间隔t∈[tk,tk+1)内,可得:
其中,ρi(t),ηi(t)和μi(t)是控制器设计的控制力矩信号与基于事件触发机制构建的传输控制力矩信号之间的关系的系数,都是连续时变标量,|ρi(t)|、|ηi(t)|、|μi(t)|分别为ρi(t)、ηi(t)、μi(t)的绝对值,并且当i=1,2,3时,满足|ρi(t)|≤1,|ηi(t)|≤1和|μi(t)|≤1。
进一步推得基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式:
步骤(3):利用反步法设计自适应控制器:
首先,引入两个新的变量z1=qe,z2=we-α1,分别求得z1、z2相对时间的一阶导数如下:
其中,引入的新变量z1为3×1维矩阵,z1=[z11,z12,z13]T,是z1相对时间的导数,是z1的转置矩阵;S(z1)为斜对称矩阵,表达形式为 是引入的新变量z2相对时间的导数;与步骤(2)的基于事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式相结合,构建时变控制增益矩阵K(t)=diag{k1(t),k2(t),k3(t)},则基于所述事件触发机制构建的传输控制力矩信号即航天器受到执行机构实际的控制力矩信号与控制器设计的控制力矩信号之间的关系式可表示为为方便分析,构建航天器受到的外界干扰和为
基于反步法,设计李亚普诺夫方程:
对李亚普诺夫方程求导数:
设计虚拟控制量:
α1=-c1z1
其中,c1>0为虚拟控制器的增益。优选地,c1取值为c1=1。
利用非线性回归矩阵转化模型得到:
其中,为书写简便记且f是3×1维矩阵,表示为f=[f1,f2,f3]T;L为线性算子,将f从3×1维空间映射到3×6维空间,表示为L(f)T为L(f)的转置矩阵;L(w)表示为是惯性矩阵J的线性回归量。
因为d虽然未知但是有界,可知Td有界,结合触发事件可知存在一个常数km使在t≥0的任何时间内都有K(t)的最小特征值λmin(K(t))≥km>0。为了处理事件触发的时变性与d的未知性,进行界估计,定义:
其中,l表示对控制增益矩阵K(t)的上界估计;p表示对外界干扰和Td上界估计的范数与时变控制增益矩阵K(t)的上界估计的乘积。并令 是l的估计值,是l与其估计值的误差,是p的估计值,是p与其估计值的误差。
对未知变量Td进行处理:
其中,||z2||是z2的二范数。
其中,为书写简便记ε(t)是有界的光滑辅助信号,ε(t)>0,满足取δ1和δ2是辅助信号的常数,满足δ1>0,δ2>0,辅助信号将用于保证系统的渐进收敛,只要满足条件可以任意选取,优选地,取ε(t)=0.2×e-0.01t。令 是的估计值,是与其估计值的误差。从而构建李雅普诺夫二次函数V2:
推导出V2的导数范围:
其中,c2,γ1,γ2是控制器增益参数,是正的常数。优选地,c2,γ1,γ2的取值分别为c2=1,γ1=12,γ2=12;Γ1是控制器增益矩阵,为6×6的正对称矩阵,优选地,Γ1的取值为Γ=20;为书写简便记 和分别为估计值和的自适应律。
根据李雅普诺夫二次函数V2得到更新律:
控制器设计的控制力矩信号u(t):
控制信号可以表示为:
综上,可以得到V2的导数的范围:
得到闭环李亚普诺夫函数:
V=V1+V2
得到V的导数范围如下:
结合辅助信号ε(t)的定义,可得系统渐进稳定。
通过Matlab仿真,得到基于事件触发机制的自适应航天器姿态跟踪控制方法,可以实现在惯性矩阵不确定、外界干扰、以及引入的事件触发机制具有时变性的影响下,卫星相对姿态跟踪误差渐进收敛;同时引入的事件触发使实际传输的控制力矩信号的更新次数大幅度减小,有效减小通信量,其效果可达到在100赫兹的采样频率下,基于事件触发机制构建的传输的控制力矩信号更新次数可以减小约97%。基于事件触发机制的自适应航天器姿态跟踪控制方法,通过定义一个事件触发机制,从而大大减少通信量和控制执行频率。同时,在控制中引入了界估计解决时变参数问题,引入辅助信号并利用界估计,解决了未知惯性参数和外部干扰的影响。
以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于事件触发机制的自适应航天器姿态跟踪控制方法,其特征在于包括以下步骤:
(1)建立含有外部干扰的航天器动力学模型与相对姿态跟踪误差模型;
(2)定义控制过程中应用的事件触发机制,基于所述事件触发机制构建传输控制力矩信号与控制器设计的控制力矩信号之间的关系式;
(3)基于步骤(1)建立的航天器动力学模型与相对姿态跟踪误差模型、步骤(2)所述的基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式,利用反步法对航天器姿态进行自适应控制;
所述步骤(1)的航天器动力学模型与相对姿态跟踪误差模型为:
其中,J为所述航天器的转动惯量矩阵,且是3×3的正定对称矩阵;w=[w1,w2,w3]T为所述航天器在本体坐标系下相对于惯性坐标系的姿态角速度,w1,w2,w3分别为姿态角速度在本体坐标系的x轴、y轴和z轴上的分量;we=[we1,we2,we3]T为所述航天器在本体坐标系上相对于期望坐标系的姿态角速度,为we相对时间的导数;wd=[wd1,wd2,wd3]T为所述航天器在期望坐标系相对于惯性坐标系的姿态角速度,为wd相对时间的导数;S(w)、S(we)是斜对称矩阵,其形式分别为t为时间;τ(t)=[τ1(t),τ2(t),τ3(t)]T为所述航天器的执行机构实际产生的控制力矩信号,因为不考虑所述航天器的执行机构受到的干扰,所以等于基于所述事件触发机制构建的传输控制力矩信号;d为所述航天器受到的空间环境中的干扰力矩信号,其值未知但是有界;q0e、qe=[qe1,qe2,qe3]T分别为所述航天器由期望坐标系相对本体坐标系的相对姿态跟踪误差四元数的标量和向量,分别是q0e、qe相对时间的导数,是的转置矩阵,且满足C为期望姿态到所述航天器实际姿态的转换矩阵,是3×3的矩阵,且||C||为C的二范数,且||C||=1,为C的导数,且S(qe)是斜对称矩阵,其形式为I是3×3的单位矩阵;
所述步骤(2)的基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式为:
2.根据权利要求1所述的一种基于事件触发机制的自适应航天器姿态跟踪控制方法,其特征在于:基于步骤(1)建立的航天器动力学模型与相对姿态跟踪误差模型、步骤(2)基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式,利用反步法对所述航天器进行自适应控制设计:
基于反步法,首先引入两个新变量z1=qe,z2=we-α1,设计虚拟控制量为:
α1=-c1z1
其中,c1>0为虚拟控制器的增益;
基于所述虚拟控制量,设计自适应控制器和自适应更新律:
其中,为书写简便记 为ud的转置矩阵;为书写简便记||z2||为新变量z2的二范数,为z2的转置矩阵;为书写简便记且f是3×1维矩阵,表示为f=[f1,f2,f3]T;L为线性算子,将f从3×1维空间映射到3×6维空间,表示为L(f)T为L(f)的转置矩阵;L(w)表示为L(w)T为L(w)的转置矩阵;S(w)T是S(w)的转置矩阵;Γ1是控制器增益矩阵,为6×6的正对称矩阵;α1为虚拟控制量,为α1相对时间的导数;θ=[J11,J22,J33,J23,J13,J12]T为惯性矩阵J的线性回归量,为θ的估计值,为估计值的自适应律;与所述步骤(2)的基于所述事件触发机制构建的传输控制力矩信号与控制器设计的控制力矩信号之间的关系式相结合,构建时变控制增益矩阵K(t)=diag{k1(t),k2(t),k3(t)},则基于所述事件触发机制构建的传输控制力矩信号即航天器受到执行机构实际的控制力矩信号与控制器设计的控制力矩信号之间的关系式可表示为为方便分析,构建航天器受到的外界干扰和为对时变控制增益矩阵K(t)进行界估计,则存在常数km使在t≥0的任何时间内都有K(t)的最小特征值λmin(K(t))≥km>0;表示对时变控制增益矩阵K(t)的上界估计,是l的估计值,为估计值的自适应律;表示外界干扰和Td上界估计的范数与时变控制增益矩阵K(t)的上界估计的乘积,是p的估计值,为估计值的自适应律;ε(t)是有界的光滑辅助信号,ε(t)>0,满足取δ1和δ2是辅助信号的常数,满足δ1>0,δ2>0;γ1,γ2是控制器增益参数,是正的常数。
4.根据权利要求2所述的一种基于事件触发机制的自适应航天器姿态跟踪控制方法,其特征在于:选取的有界的光滑辅助信号为ε(t)=0.2×e-0.01t。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810107715.5A CN108445748B (zh) | 2018-02-02 | 2018-02-02 | 一种基于事件触发机制的自适应航天器姿态跟踪控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810107715.5A CN108445748B (zh) | 2018-02-02 | 2018-02-02 | 一种基于事件触发机制的自适应航天器姿态跟踪控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108445748A CN108445748A (zh) | 2018-08-24 |
CN108445748B true CN108445748B (zh) | 2020-07-03 |
Family
ID=63191528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810107715.5A Active CN108445748B (zh) | 2018-02-02 | 2018-02-02 | 一种基于事件触发机制的自适应航天器姿态跟踪控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108445748B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109358497B (zh) * | 2018-09-14 | 2020-04-21 | 北京航空航天大学 | 一种基于b样条函数的卫星路径规划和预测控制的跟踪方法 |
CN109471448A (zh) * | 2018-12-27 | 2019-03-15 | 西北工业大学 | 一种基于事件驱动的柔性航天器姿态控制方法 |
CN109933078B (zh) * | 2019-03-29 | 2022-02-08 | 西北工业大学 | 鲁棒微分博弈控制的多微卫星姿态接管失效航天器方法 |
CN110109351B (zh) * | 2019-04-08 | 2021-05-18 | 广东工业大学 | 一种基于指定性能的多智能体一致性控制方法 |
CN111284732B (zh) * | 2020-03-13 | 2021-08-03 | 北京航空航天大学 | 一种基于事件触发通信的航天器抗干扰姿态协同控制方法 |
CN111596678B (zh) * | 2020-05-07 | 2021-07-06 | 北京航空航天大学 | 一种航天器逼近空间翻滚目标的相对位置自适应控制方法 |
CN112558471B (zh) * | 2020-11-24 | 2022-05-24 | 西北工业大学 | 一种基于动态事件触发的航天器编队离散分布式非合作博弈方法 |
CN113050493B (zh) * | 2021-03-19 | 2022-03-04 | 大连理工大学 | 一种网络化环境下的小车倒立摆系统输出反馈控制方法 |
CN113219826B (zh) * | 2021-04-09 | 2022-03-15 | 杭州电子科技大学 | 基于事件触发机制的非线性振荡器的自适应跟踪控制方法 |
CN113608540B (zh) * | 2021-06-29 | 2022-07-05 | 北京理工大学 | 一种比例积分型事件触发航天器姿态控制方法 |
CN115686048B (zh) * | 2022-10-31 | 2023-06-02 | 哈尔滨工业大学 | 执行器受限航天器交会系统的动态触发有限时间控制方法 |
CN115857329B (zh) * | 2022-11-07 | 2024-03-08 | 西北工业大学深圳研究院 | 基于事件触发的绳系航天器复合体姿态稳定控制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105159307B (zh) * | 2015-08-27 | 2018-01-09 | 北京天航华创科技股份有限公司 | 一种带有执行器饱和问题的运动体姿态事件触发控制方法 |
CN106557088A (zh) * | 2017-02-15 | 2017-04-05 | 大连海事大学 | 一种基于事件触发机制的船舶航向控制器及其控制方法 |
CN107577145B (zh) * | 2017-08-25 | 2020-06-09 | 湘潭大学 | 编队飞行航天器反步滑模控制方法 |
-
2018
- 2018-02-02 CN CN201810107715.5A patent/CN108445748B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108445748A (zh) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108445748B (zh) | 一种基于事件触发机制的自适应航天器姿态跟踪控制方法 | |
Shao et al. | Robust dynamic surface trajectory tracking control for a quadrotor UAV via extended state observer | |
Song et al. | 1 bit encoding–decoding-based event-triggered fixed-time adaptive control for unmanned surface vehicle with guaranteed tracking performance | |
CN107577145B (zh) | 编队飞行航天器反步滑模控制方法 | |
CN113589689B (zh) | 一种基于多参数自适应神经网络的滑模控制器设计方法 | |
Qi et al. | Stable indirect adaptive control based on discrete-time T–S fuzzy model | |
Lu et al. | Adaptive fuzzy output feedback control for nonlinear systems based on event-triggered mechanism | |
CN110543184B (zh) | 一种刚性飞行器的固定时间神经网络控制方法 | |
CN110116409B (zh) | 一种基于扰动观测器的四通道遥操作双边控制方法 | |
CN109062240B (zh) | 一种基于神经网络估计的刚性飞行器固定时间自适应姿态跟踪控制方法 | |
CN111198570B (zh) | 一种基于固定时间微分器预测的抗时延高精度自抗扰姿态控制方法 | |
You et al. | Reliable control for flexible spacecraft systems with aperiodic sampling and stochastic actuator failures | |
CN107203138B (zh) | 一种输入输出饱和的飞行器鲁棒控制方法 | |
Zhang et al. | A novel event-triggered robust neural formation control for USVs with the optimized leader–follower structure | |
An et al. | Robust fixed-time tracking control for underactuated AUVs based on fixed-time disturbance observer | |
Zhang et al. | Anti-disturbance control for dynamic positioning system of ships with disturbances | |
Dian et al. | Gain scheduled dynamic surface control for a class of underactuated mechanical systems using neural network disturbance observer | |
Wang et al. | Globally stable adaptive dynamic surface control for cooperative path following of multiple underactuated autonomous underwater vehicles | |
Xie et al. | Finite-time command-filtered approximation-free attitude tracking control of rigid spacecraft | |
CN110687796B (zh) | 基于神经网络的自适应有限时间命令滤波反步控制方法 | |
CN111439392B (zh) | 一种航天器编队位置协同控制方法 | |
Zeng et al. | Finite-time coordinated formation control of discrete-time multi-AUV with input saturation under alterable weighted topology and time-varying delay | |
CN111427264B (zh) | 一种复杂遥操作技术的神经自适应固定时间控制方法 | |
CN113220007A (zh) | 执行机构故障的挠性航天器有限时间姿态协同控制方法 | |
Tan et al. | Rotation matrix-based finite-time attitude coordinated control for spacecraft |
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 |