[0004] 本发明针对现有技术在只能获取航天器交会系统相对位置信息时,提出一种基于事件触发的航天器交会系统有限时间输出反馈控制方法,实现航天器交会系统有效控制。
[0005] 本发明基于低增益反馈、事件触发控制,考虑到实际系统中执行器饱和以及系统状态不可直接量测的情况,同时,为了实现系统在有限时间内稳定,设计了基于事件触发的航天器交会系统有限时间输出反馈控制器。
[0006] 本发明的具体步骤是:
[0007] 步骤1、建立航天器交会系统状态空间模型
[0008] 考虑航天器交会系统的C‑W方程:
[0009]
[0010] 假设目标航天器运行在半径为R的圆形轨道上,建立目标航天器轨道坐标系o‑xyz,坐标轴x是圆轨道半径的方向,坐标轴y是追踪航天器运行的方向,坐标轴z垂直于目标航天器相对地球质心运动的平面,并且方向与坐标轴x和y构成右手坐标系,原点o是目标航天器的质心。x,y,z分别表示追踪航天器与目标航天器在x、y、z轴方向上的相对距离。ax,ay,az分别为坐标轴x、y、z方向上的加速度分量。ωx,ωy,ωz分别为这三个加速度分量的最大值。sat(·)表示单位饱和函数。引力常数μ=GM,M是被环绕的星球质量,G为万有引力常数。可以计算得到目标航天器的轨道角速度
[0011] 选择如下状态向量,
[0012]
[0013] 建立如下航天器交会系统状态方程
[0014]
[0015] 其中 M=diag{ωx,ωy,ωz}, 矩阵A和矩阵B如下所示:
[0016]
[0017]
[0018] 考虑到只能获取航天器交会系统的相对位置信息,因此
[0019] Y=CX,
[0020] Y=[x y z]T表示航天器交会系统的控制输出。其中
[0021]
[0022] 进一步得到航天器交会系统的状态空间模型
[0023]
[0024] 步骤2、事件触发条件设计
[0025] 事件触发条件设计为
[0026]
[0027] 其中, 表示观测器对系统状态X的观测值,
[0028] eY(t)=Y(t)‑Y(tk),
[0029] 表示系统当前控制输出Y(t)与事件触发时刻的控制输出Y(tk)的差值。
[0030] t∈[tk,tk+1),k∈N,N表示自然数集合,tk是事件触发时刻。 是事件触发的参数,且
[0031] 步骤3、时变参数设计
[0032] 时变参数设计为
[0033]
[0034] 有限时间T定义为
[0035]
[0036] 其中
[0037]
[0038] ξ0=ξ(t0),t0表示系统初始时刻。θc=θc(ξ0)≥1是一个常数,且标量θc(ξ0)可以通过下式得到
[0039] θc(ξ0)=6ξ0λmax(U(ξ0)W(ξ0)‑1),
[0040] W(ξ0)和U(ξ0)分别是下列Lyapunov方程的唯一对称正定解
[0041]
[0042]
[0043] θo=θo(ξ0)≥1是一个常数,且标量θo(ξ0)可以通过下式得到
[0044] θc(ξ0)=6ξ0λmax(U(ξ0)Wo(ξ0)‑1),
[0045] Wo(ξ0)和Uo(ξ0)分别是下列Lyapunov方程的唯一对称正定解
[0046]
[0047]
[0048] 步骤4、设计状态观测器
[0049] 当(A,C)是可观测的,设计如下状态观测器
[0050]
[0051] 其中 是系统状态X的观测值,L是观测器的增益矩阵。
[0052] L=Q(ξ(t))CT,
[0053] Q(ξ(t))是下列参量Lyapunov方程的唯一对称正定解
[0054] ATQ(ξ(t))+Q(ξ(t))A‑Q(ξ(t))CTCQ(ξ(t))=‑ξ(t)Q(ξ(t))。
[0055] 令 根据航天器交会系统模型以及状态观测器,进一步得到
[0056]
[0057] 步骤5、控制器设计
[0058] 设计如下有限时间输出反馈控制器,
[0059] U=‑BTP(ξ(t))(X‑e),
[0060] P(ξ(t))∈R6×6是下列参量Lyapunov方程的解
[0061] ATP(ξ(t))+P(ξ(t))A‑P(ξ(t))BBTP(ξ(t))=‑ξ(t)P(ξ(t)),
[0062] 步骤6、设计椭球集合
[0063] 首先,定义如下两个集合,
[0064]
[0065]
[0066] ‖‖表示矩阵或向量的2范数, 是一个椭球集,当X(t)包含于集合 中时,执行器不发生饱和。
[0067] 通过计算可以得知,
[0068]
[0069] 也就是说,对于任意的 执行器不会发生饱和。
[0070] 即,
[0071] sat(BTP(ξ(t))(X‑e))=BTP(ξ(t))(X‑e)
[0072] 步骤7、建立闭环系统状态空间模型
[0073] 将所设计的有限时间输出反馈控制器代入航天器交会系统的状态空间模型中,得到如下闭环系统状态空间模型
[0074]
[0075] 考虑到对于任意的 执行器不会发生饱和。进一步化简得到如下闭环系统状态空间模型
[0076]
[0077] 步骤8、闭环系统的稳定性分析
[0078] 根据Lyapunov稳定性理论,选择如下多Lyapunov函数
[0079] V(X,e)=πX(ξ(t))XTP(ξ(t))X+πe(ξ(t))eTQ‑1(ξ(t))e,
[0080] 其中 πe(ξ(t))=78ξ(t)tr(P(ξ(t)))tr(Q(ξ(t)))。
[0081] 令
[0082] VX(X,t)=πX(ξ(t))XTP(ξ(t))X,
[0083] VX(X,t)对时间t∈[tk,tk+1)求导可以得到
[0084]
[0085] 令
[0086] Ve(e,t)=πe(ξ(t))eTQ‑1(ξ(t))e,
[0087] Ve(e,t)对时间t∈[tk,tk+1)求导可以得到
[0088]
[0089] 最终可以得到
[0090]
[0091] 其中,
[0092]
[0093]
[0094] 将设计的时变参量ξ(t)带入 得到
[0095]
[0096] 这也就说明了闭环系统有限时间T内稳定。
[0097] 本发明考虑到航天器交会系统只能获得航天器相对位置信息的情况,基于低增益反馈,事件触发控制和全阶状态观测器设计了一种基于事件触发的有限时间输出反馈控制器。所设计的控制器避免了执行器饱和并且节约了系统计算资源,使得航天器交会系统在有限时间内稳定并且所设计的状态观测器的状态估计值趋于系统状态真实值。利用本发明的方法,可以在航天器交会系统状态不可直接测量的情况下,实现两个航天器在有限时间内完成交会任务。