发明内容
本发明的目的是提供基于事件触发的绳系航天器复合体姿态稳定控制方法,实现了系统的在内部不确定性和外部干扰的情况下姿态的精确稳定控制,利用事件触发机制,达到节省通信资源的技术效果。
本发明采用以下技术方案:
本发明的实施例一提供了一种基于事件触发的绳系航天器复合体姿态稳定控制方法,包括以下步骤:
根据复合体系统的姿态误差计算复合体系统的第一状态变量;
根据复合体系统的角速度误差值以及第一状态变量计算复合体系统的第二状态变量;
根据第二状态变量计算复合体系统的自适应策略;
根据第二状态变量和自适应策略生成复合体系统的事件触发控制量;
将自适策略和事件触发控制量根据事件触发机制实时生成复合体系统的控制信号。
可选地,根据复合体系统的姿态误差计算复合体系统的第一状态变量包括:
根据实际姿态和期望姿态计算姿态误差;
将姿态误差值赋值给第一状态变量。
可选地,根据实际姿态和期望姿态计算姿态误差包括:
计算实际姿态与期望姿态的逆的卷积,得到姿态误差。
可选地,根据复合体系统的角速度误差值以及第一状态变量计算复合体系统的第二状态变量包括:
计算角速度误差值与角速度误差中间变量的值的差值,得到第二状态变量。可选地,角速度误差中间变量的计算方式为:
α1=-c1G-1z1,
其中,α1为角速度误差中间变量,z1为第一状态变量,c1为大于0的设计参数,G-1为系统向量的逆。
可选地,根据第二状态变量计算复合体系统的自适应策略包括:
通过第二状态变量计算得到复合体系统的自适应策略的计算方式为:
其中,z2为第二状态变量,为操作算子/>的转置,/>为转动惯量估计值,Γ为正定对角矩阵,λ为大于0的设计常数,/>为自适应策略。
可选地,根据第二状态变量和所述自适应策略生成所述复合体系统的事件触发控制量包括:
通过所述第二状态变量和所述自适应策略参与如下计算,得到所述复合体系统的所述事件触发控制量,该计算方式为:
其中,ue(t)为事件触发控制量,γ代表事件触发控制量触发前后的误差放缩系数,δ为大于0的设计常数,为所述自适应策略,ξT为向量ξ的转置。
可选地,根据所述事件触发控制量按照事件触发机制实时生成所述复合体系统的控制信号包括:
按照所述事件触发机制执行以下步骤:
计算当前时刻所述事件触发控制量的各个分量与对应触发后的控制量的各个分量之间的差值;
计算各个所述触发后的控制量与误差放缩系数的乘积;
进一步计算所述乘积与触发阈值的和值;
将断所述差值与所述和值进行比较,根据比较结果选择下一时刻的控制信号。
可选地,所述事件触发机制为:
其中,δ>0,和1>γ>0γ代表未触发的控制ue(t)和事件触发后控制u(t)之间的误差放缩系数;/>表示事件触发阈值。/>为未触发的事件控制量,其中uei(t),i=1,2,3为其三个方向的分量。/>触发后的事件控制量,其中ui(t),i=1,2,3为其三个方向的分量。
本发明的实施例二提供了一种基于事件触发的绳系航天器复合体姿态稳定控制系统,控制系统包括:
复合体姿态计算装置,该装置包括:
状态变量生成模块,根据复合体系统的姿态误差计算复合体系统的第一状态变量和根据复合体系统的角速度误差以及第一状态变量计算复合体系统的第二状态变量;
自适应策略生成模块,用于根据第二状态变量计算复合体系统的自适应策略;
传感器:用于接收状态变量生成模块发送的第一状态变量和第二状态变量;
控制器,用于接收传感器采集的第二状态变量,以及接收自适应策略生成模块生成的自适应策略,并生成复合体系统的事件触发控制量;
事件触发生成器,用于接收控制器发送的事件触发控制量,并根据事件触发机制实时生成复合体系统的控制信号。
可选地,控制系统还包括:执行器,用于接收事件触发器发送的控制信号,并基于控制信号实时作用复合体姿态计算装置。
本发明的有益效果是:本发明提供的基于事件触发的绳系航天器复合体姿态稳定控制方法具有以下几点有益效果:
1)、基于反步法,设计绳系复合体的高性能姿态稳定控制方案,实现多种不确定性和外界干扰存在下的快速姿态稳定控制;
2)、通过计算出的自适应策略,实现对系统多种不确定性信息的精确估计;
3)、设计事件触发机制,降低控制系统执行器的更新频次,节省通信资源。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
如图1所示,本申请提供的一种基于事件触发的绳系航天器复合体姿态稳定控制方法所采用的绳系复合体系统坐标定义图。
如图2所示,基于事件触发的绳系航天器复合体姿态稳定控制方法包括:
步骤S101、根据系统的姿态误差计算系统的第一状态变量;
步骤S102、根据复合体系统的角速度误差值以及第一状态变量计算复合体系统的第二状态变量;
步骤S103、根据第二状态变量计算系统的自适应策略;
步骤S104、根据第二状态变量生成系统的事件触发控制量;
步骤S105、根据事件触发控制量按照事件触发机制实时生成复合体系统的控制信号。
在一种实施例中,本实施例一的方法首先根据修正罗德里格参数法建立复合体系统的复合体姿态动力学方程
其中为复合体姿态角速度,ω1,ω2,ω3为角速度在本体坐标系下的分量,ω×为复合体姿态角速度角速度ω的斜对称矩阵,其可表示为ω=[0 -ω3ω2;ω3 0 -ω1;-ω2 ω1 0];对称正定矩阵/>为复合体系统的未知转动惯量,/>表示由太阳光压和重力场等因素造成的有界外部干扰;/>为复合体的实际姿态,其与系统姿态关系可以表示为/>其中,/>为欧拉转动矢量和/>为角度。非奇异矩阵/>其中/>为单位矩阵。
σ代表复合体的实际姿态,σd代表复合体的期望姿态,则姿态误差表示为同理,实际姿态角速度ω和期望角速度ωd的误差可表示为ωe=ω-R(σe)ωd=[ωe1,ωe2,ωe3]T。
基于上述复合体系统模型得出复合体的姿态误差动力学模型为:
其中
在一种实施例中,在上述姿态误差动力学模型的基础上,进一步基于反步法求解复合体的第一状态变量和第二装填变量。定义第一状态变量为z1,第二状态变量为z2:
步骤S101、根据系统的姿态误差计算系统的第一状态变量包括:根据实际姿态和期望姿态计算姿态误差;将姿态误差值赋值给第一状态变量。
其中,根据上述系统模型中的σ和σd,计算实际姿态与期望姿态的逆的卷积,得到姿态误差。
计算方式为:
其中,σ代表复合体的实际姿态和σd代表复合体的期望姿态,σe为姿态误差,将计算得到的姿态误差赋值给第一状态变量。
步骤S102、根据复合体系统的角速度误差值以及第一状态变量计算复合体系统的第二状态变量包括:
1)、计算角速度误差值与角速度误差中间变量的值的差值,得到第二状态变量。
角速度误差值计算方式为ωe=ω-R(σe)ωd=[ωe1,ωe2,ωe3]T,其中ω为实际角速度值和ωd为期望角速度值,ωe为角速度误差值;
可选地,角速度误差中间变量的计算方式为:
α1=-c1G-1z1,
其中,α1为角速度误差中间变量,z1为第一状态变量,c1为大于0的设计参数,G-1为系统向量的逆。
2)、计算角速度误差值与角速度误差中间变量的值的差值,得到第二状态变量。
需要说明的是,上述步骤1到步骤2的推导过程如下:
z1=σe,z2=ωe-α1, (3)
考虑李雅普诺夫函数其导函数可表示为:
为了对复合体姿态进行稳定控制,选择如下角速度误差中间变量:
α1=-c1G-1z1, (5)
其中c1>0为设计参数,将虚拟控制输入(5)代入可得:
可选地,根据系统的角速度误差值以及第一状态变量计算系统的第二状态变量包括:
可选地,角速度误差中间变量的计算方式为:
α1=-c1G-1z1,
在一种实施例中,同样基于反步法以及自适应求解技术,进行如下一系列计算操作:
为了驱使角速度误差ωe趋近α1,根据公式(2)得到复合体系统误差方程:
其中,为已知向量。由于转动惯量J具有不确定性,本发明引入操作算子Y将列向量s=[s1,s2,s3]T转化为3×6的如下矩阵形式
此外,令Jij代表转动惯量的分量并且有因此函数表示为:
Jω=Y(ω)Φ,Jφ=Y(φ)Φ, (9)
将(9)代入动态方程(7)中可得:
其中,且Φ为未知列向量。
对求取微分并基于Young不等式,可得:
其中,且/>为未知常数。构建如下形式的李雅普诺夫函数:
其中,为参数估计误差,Φ为系统未知参数向量,/>为未知参数向量的自适应估计值,/>为所设计的正定对角矩阵。将(6)和(11)代入/>中,可得:
其中,且ε2>0。
可选地,根据第二状态变量和所述自适应策略生成所述复合体系统的事件触发控制量包括:
通过所述第二状态变量和所述自适应策略参与如下计算,得到所述复合体系统的所述事件触发控制量,该计算方式为:
其中,ue(t)为事件触发控制量,γ代表事件触发控制量触发前后的误差放缩系数,δ为大于0的设计常数,为所述自适应策略,ξT为向量ξ的转置。
在一种实施例中,为了实现系统姿态的快速稳定控制,设计生成如下的事件触发控制量和自适应策略:
其中,ue(t)为所述事件触发控制量,γ代表事件触发控制量触发前后的误差放缩系数,δ为大于0的设计常数,为所述自适应策略,ξT为向量ξ的转置。
其中,z2为第二状态变量,为操作算子/>的转置,/>为转动惯量估计值,Γ为正定对角矩阵,λ为大于0的设计常数,/>为自适应策略。
步骤S103经过根据上述公式(6)-(13)的计算,最后,将各个系统参数以及求解出的第二状态变量计算参与公式(15)计算得到复合体系统的自适应策略,所述自适应策略是表示转动惯量的Jij的各个分量构成的矩阵,例如:
需要说明的是,基于反步法生成高精度姿态控制的自适应策略,并利用自适应技术实现对未知信息的精准估计。
步骤S104:根据第二状态变量和自适应策略生成复合体系统的事件触发控制量;
在一种实施例中,将上述求解出的自适应策略代入公式(16)中,求解出ξ,
再将ξ、ξ的转置和第二状态变量代入公式(14)进行计算,得到事件触发控制量。
步骤S105、根据事件触发控制量按照事件触发机制实时生成复合体系统的控制信号。
按照所述事件触发机制执行以下步骤:
计算当前时刻所述事件触发控制量的各个分量与对应触发后的控制量的各个分量之间的差值;
计算各个所述触发后的控制量与误差放缩系数的乘积;
进一步计算所述乘积与触发阈值的和值;
将断所述差值与所述和值进行比较,根据比较结果选择下一时刻的控制信号。
在一种实施例中,根据事件触发机制将当前时刻的事件触发控制量的各个分量基于阈值法进行选择,将满足阈值条件的触发事件量的分量作为下一时刻的复合体系统的控制信号。
具体事件触发机制设计如下:
其中,δ>0,和1>γ>0γ代表未触发的控制ue(t)和事件触发后控制u(t)之间的误差放缩系数;/>表示事件触发阈值。/>为未触发的事件控制量,其中uei(t),i=1,2,3为其三个方向的分量。/>触发后的事件控制量,其中ui(t),i=1,2,3为其三个方向的分量。
表示采样误差,/>代表控制器更新时刻。
在一种实施例中,当前时刻为tk,其中当前时刻的事件触发控制量的各个分量为ue1(t)、ue2(t)、ue3(t),t=tk,对应触发后的控制量的各个分量u1(t)、u2(t)、u3(t),t=tk,
计算当前时刻的事件触发控制量的各个分量与对应触发后的控制量的各个分量的差值,得到:ue1(t)-u1(t)、ue2(t)-u2(t)、ue3(t)-u3(t)。
计算各个触发后的控制量与误差放缩系数的乘积,例如:γ|u1(t)|、γ|u2(t)|、γ|u3(t)|,其中γ代表未触发的控制ue(t)和事件触发后控制u(t)之间的误差放缩系数;
进一步计算乘积与触发阈值的和值,例如: 表示事件触发阈值;
进一步计算乘积与触发阈值的和值;将差值与和值进行比较,根据比较结果选择下一时刻的控制信号。
在一种实施例中,将差值与和值进行比较,当满足时,将ue1(t)作为下一时刻的控制信号;当满足/>时,将ue2(t)作为下一时刻的控制信号;当满足/>时,将ue3(t)作为下一时刻的控制信号。
注意到,一旦上式被触发,时刻便被标记为tk+1并且更新后的控制量u(tk+1)便会传送到响应的执行器端,可以对传感器采集的姿态数据进行实时控制。除此之外,在时间段t∈[tk,tk+1)内,控制信号将保持常值u(t)=u(tk)。
本发明从步骤S104到步骤S105还针对事件触发控制量设计了对应的事件触发机制,其可以在控制精度得到保证的同时降低通信频次,减少通信资源的浪费。
基于本发明所述提供的基于事件触发的绳系航天器复合体姿态稳定控制方法,还提供过了一种对闭环系统进行李雅普诺夫稳定性证明方法。
根据所设计的事件触发机制,显然控制输入在tk+1时刻会被更新u(tk+1),这意味着在时间段t∈[tk,tk+1)内不等式成立。因此,在t∈[tk,tk+1)存在满足|μ1i(t)|≤1和|μ2i(t)|≤1的函数μ1(t),μ2(t)使得成立,即有/>其中/>和其矢量形式为:
上式中和/>
根据公式(18)和公式(15)代入的表达式中,可得
由引理1和Young不等式可得:
根据定义可知为有界函数,且满足/>其中δm为未知常数。将(21)和(22)代入(20)得:
将代入上式,可得
其中,λmin为矩阵Γ-1的最小特征值。
令从上式可得/>由此可得,z1,z2,/>均是有界的变量。同理,可以检查验证闭环系统信号均一致有界。且跟踪误差可以收敛到区域/>内。
现在,证明存在常数υ使得对于均有tk+1-tk≥υ。为了达到此目的,我们有/>且/>对任意t∈[tk,tk+1)均成立。此外,由上式可以看出u(t)为列向量/>的连续函数,因此/>以及/>均为有界函数。故存在常数/>使得下式成立。
其中ui和为向量u和/>的第i个分量。另一方面,由事件触发定义(18)可知,
定义并且结合(25)和(26)可得tk+1-tk≥υ成立。因此,Zeno现象可以有效避免,且系统稳定性可以证明。
基于上述方法的一种实验仿真结果及性能分析如下:
结合图3~图6。图3中可以看出姿态跟踪误差在10秒内快速收敛到原点附近;图4表示跟踪角速度误差也实现了快速跟踪收敛。图5为控制系统的控制输入量示意图;图6表示在事件触发机制下,控制输入保证高精度跟踪的同时减少了触发频率与通讯次数。
本发明的实施例二提供了一种基于事件触发的绳系航天器复合体姿态稳定控制系统,如图7所示,控制系统包括:
复合体姿态计算装置20,该装置包括:
状态变量生成模块201,根据复合体系统的姿态误差计算复合体系统的第一状态变量和根据复合体系统的角速度误差以及第一状态变量计算复合体系统的第二状态变量;
自适应策略生成模块202,用于根据第二状态变量计算复合体系统的自适应策略;
传感器203:用于接收状态变量生成模块发送的第一状态变量和第二状态变量;
控制器204,用于接收传感器采集的第二状态变量,以及接收自适应策略生成模块生成的自适应策略,并生成复合体系统的事件触发控制量;
事件触发生成器205,用于接收控制器发送的事件触发控制量,并根据事件触发机制实时生成复合体系统的控制信号。
可选地,控制系统还包括:执行器206,用于接收事件触发器发送的控制信号,并基于控制信号实时作用复合体姿态计算装置。
需要说明的是该系统所能实现的基于事件触发的绳系航天器复合体姿态稳定控制方法与实施例一的方法一致,此处不再赘述。