CN113608540B - 一种比例积分型事件触发航天器姿态控制方法 - Google Patents

一种比例积分型事件触发航天器姿态控制方法 Download PDF

Info

Publication number
CN113608540B
CN113608540B CN202110726476.3A CN202110726476A CN113608540B CN 113608540 B CN113608540 B CN 113608540B CN 202110726476 A CN202110726476 A CN 202110726476A CN 113608540 B CN113608540 B CN 113608540B
Authority
CN
China
Prior art keywords
attitude
spacecraft
signal
feedback controller
proportional integral
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
Application number
CN202110726476.3A
Other languages
English (en)
Other versions
CN113608540A (zh
Inventor
史大威
宋继良
李双汐
李春晖
王军政
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202110726476.3A priority Critical patent/CN113608540B/zh
Publication of CN113608540A publication Critical patent/CN113608540A/zh
Application granted granted Critical
Publication of CN113608540B publication Critical patent/CN113608540B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0833Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using limited authority control

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Security & Cryptography (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (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

一种比例积分型事件触发航天器姿态控制方法
技术领域
本发明属于事件触发控制与航天器姿态控制的技术领域,具体涉及一种比例积分型事件触发航天器姿态控制方法。
背景技术
航天器姿态控制是航天器在轨运行顺利执行各类任务的关键前提。在轨运行的航天器,为了完成各类任务,需要运行在一定的参考姿态下,保证航天器上的载荷能够工作。同时,由于各类空间环境摄动以及航天器自身不确定性的影响,航天器姿态容易因干扰而发生改变,需要定期进行控制以保持姿态稳定。随着航天器编队技术的发展,编队航天器之间稳定通信的需要为航天器姿态控制提出了更高的要求。现有航天器大多星上资源有限,且刚性航天器的运动学和动力学方程具有高度耦合的非线性特性,控制系统可能会受到外部干扰及内部不确定性的影响,这使得航天器的高精度姿态控制极具挑战。
针对考虑事件触发的航天器姿态控制方法,据调查与了解,目前公开的现有技术如下:申请公布号CN109471448的发明专利公开了一种基于事件驱动的柔性航天器姿态控制方法,但该方法所设计事件触发条件相对复杂,且并未考虑存在外部干扰的情况;申请公布号CN107589671的发明专利公开了一种基于事件驱动的卫星姿态控制方法,采用了混合事件驱动条件,但是同样未考虑外部干扰存在的情况;申请公布号CN112363522的发明专利公布了一种基于混杂强制型感测器的网络化航天器姿态控制方法,其利用静态事件触发机制,避免了观测器中大量数据的传输,并克服了非线性不确定项的影响,提高了航天器姿态跟踪系统的鲁棒性。
发明内容
有鉴于此,本发明提供了一种比例积分型事件触发航天器姿态控制方法,能够对抗外部环境干扰与内部不确定性对航天器姿态的干扰,同时减少航天器姿态控制中不必要的信号传递,在节省有限的星上通信资源的同时保证姿态控制性能。
实现本发明的技术方案如下:
一种比例积分型事件触发航天器姿态控制方法,包括以下步骤:
步骤一、基于四元数的卫星姿态动力学与运动学模型,构建卫星姿态跟踪误差系统模型;
步骤二、根据卫星姿态跟踪误差系统模型,确定姿态跟踪误差系统扩张状态,设计扩张状态观测器;
步骤三、利用扩张状态观测器得到系统观测状态,设计基于系统观测状态的反馈控制器;
步骤四、基于反馈控制器输出信号的采样误差,设计比例积分型事件触发策略用于决定是否更新航天器姿态的控制信号。
进一步地,若比例积分型事件触发条件满足时,采样器对反馈控制器的输出信号进行采样,若比例积分型事件触发条件未满足,则航天器姿态的控制信号通过零阶保持器对上次触发传递而来的信号进行保持得到。
进一步地,所述比例积分型事件触发条件通过对比反馈控制器输出信号采样误差和输出信号采样误差关于时间的积分值的大小,判断何时进行触发,并更新控制信号,利用经零阶保持器保持后的控制信号实现对卫星姿态的控制,使得姿态跟踪误差收敛。
进一步地,所述比例积分型事件触发条件如下:
Figure BDA0003138866750000031
tk+1,i=min{tk,ik,i,tk,i+Ti},i=1,...,n
其中,inf表示下确界,n为控制器输出信号的维度,θi为第i维信号两次触发之间满足条件的可能时间间隔,θk,i表示第k时刻到下次触发满足条件的最小时间间隔,比例系数δ1,i、积分系数δ2,i、高增益参数εi和阈值常数Mi均为设定参数,且都为正值,Ti为预设的两次触发之间的最大允许间隔,tk,i和tk+1,i分别表示第k次的触发时刻和第k+1次触发时刻,σi(·)为反馈控制器输出的第i维信号的采样误差,定义为:
Figure BDA0003138866750000032
其中
Figure BDA0003138866750000033
表示反馈控制器输出的第i维信号,
Figure BDA0003138866750000034
表示tk,i时刻的反馈控制器输出采样值经过零阶保持器保持后的信号。
有益效果:
本发明提供一种比例积分事件触发的航天器姿态控制方法。首先,通过设计扩张状态观测器与反馈控制律,有效克服外部干扰与内部不确定性的影响,提高姿态控制的精度。其次,通过对比控制器输出信号采样误差与其关于时间的积分值,决定姿态控制信号的更新时刻,实现在保证控制性能的前提下,降低控制信号的传输频率,节省通信资源的效果。最后,本发明所提供的方法相比普通比例型事件触发机制,具有更低的通信率与更高的采样效率。
附图说明
图1为本发明方法流程图。
图2为本发明原理框图。
图3卫星姿态角速度跟踪结果图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种比例积分型事件触发的姿态控制方法,包括姿态跟踪误差系统模型建立,扩张状态观测器设计,反馈控制信号设计,比例积分型事件触发机制设计。如附图1所示,首先建立四元数描述的姿态跟踪误差系统模型;其次定义系统扩张状态,设计扩张状态观测器;然后根据扩张状态观测器输出的观测值,设计反馈控制器,计算反馈控制律;最后设计基于控制信号采样误差的比例积分事件触发机制,通过事件触发条件决定是否更新航天器姿态控制律。
本方法的具体步骤如下:
(1)构建姿态跟踪误差系统模型。
考虑四元数形式的刚体航天器动力学与运动学模型:
Figure BDA0003138866750000041
Figure BDA0003138866750000042
其中ω=[ω1 ω2 ω3]∈R3为航天器角速度,u∈R3为控制力矩,d∈R3为外部干扰,J∈R3×3为航天器的惯量矩阵,I3∈R3×3为单位对角矩阵,单位四元数
Figure BDA0003138866750000043
表示航天器本体坐标系相对于惯性坐标系的姿态方向,其由矢量部分qv=[q1 q2 q3]T∈R3和标量部分q4∈R组成,并且满足
Figure BDA0003138866750000044
的约束条件。此外,
Figure BDA0003138866750000045
且对于其他向量具有类似的计算法则。
定义航天器期望的姿态四元数与角速度分别为
Figure BDA0003138866750000051
和ωd,那么进一步地,可以定义航天器姿态四元数跟踪误差与角速度跟踪误差分别为
Figure BDA0003138866750000052
Figure BDA0003138866750000053
和ωe,且
Figure BDA0003138866750000054
ωe=ω-Cωd,
其中
Figure BDA0003138866750000055
且‖C‖=1以及
Figure BDA0003138866750000056
进一步地,可以得到姿态跟踪误差系统:
Figure BDA0003138866750000057
Figure BDA0003138866750000058
由于燃料消耗和有效载荷变化等因素,航天器姿态控制惯量矩阵存在不确定性,将带有参数不确定性的惯量矩阵描述为J=J0+ΔJ,其中J0表示非奇异惯量矩阵的标称值,ΔJ表示与J相关的不确定性。由矩阵逆的性质,(J0+ΔJ)-1可以表示为:
Figure BDA0003138866750000059
其中,
Figure BDA00031388667500000510
进一步地,定义辅助状态变量z=ωe+kqev,则可以将姿态跟踪误差系统表示为:
Figure BDA00031388667500000511
其中K∈R3×3为正定矩阵,F∈R3为姿态跟踪误差系统中的非线性部分,G∈R3为姿态跟踪误差系统中的惯量矩阵不确定部分,
Figure BDA00031388667500000512
表示标称惯量矩阵的逆矩阵,
Figure BDA00031388667500000513
表示外部干扰,具体的,
Figure BDA00031388667500000514
Figure BDA0003138866750000061
Figure BDA0003138866750000062
(2)构建姿态跟踪误差系统扩张状态,设计扩张状态观测器。由步骤(1)构建的姿态误差跟踪系统模型,定义x1=z,定义扩张状态
Figure BDA0003138866750000063
则姿态跟踪误差系统可以表示为
Figure BDA0003138866750000064
Figure BDA0003138866750000065
其中
Figure BDA0003138866750000066
则对于上述包含扩张状态的系统,设计扩张状态观测器形式为:
Figure BDA0003138866750000067
Figure BDA0003138866750000068
其中
Figure BDA0003138866750000069
Figure BDA00031388667500000610
分别表示对状态x1和x2的观测值,且
Figure BDA00031388667500000611
Figure BDA00031388667500000612
b=[b1 b2 b3]T和c=[c1 c2 c3]T,其中b1~b3和c1~c3为正常数,ε=[ε1 ε2ε3]T为观测器的高增益参数,η1=[η1 η2 η3]T,且
Figure BDA00031388667500000613
Figure BDA00031388667500000614
定义为:
Figure BDA00031388667500000615
(3)根据步骤(2)设计的扩张状态观测器观测得到的系统状态,设计反馈控制律。控制律的形式如下所示:
Figure BDA00031388667500000616
其中
Figure BDA00031388667500000617
为正常数,
Figure BDA00031388667500000618
Figure BDA00031388667500000619
为饱和函数,定义为:
Figure BDA0003138866750000071
需要说明的是,控制器计算得到的控制信号
Figure BDA0003138866750000072
有三个维度,分别对应卫星姿态三个方向角速度的控制力矩。
(4)步骤(3)所设计的控制器计算得到的控制量为连续的,但是在实际的姿态控制场景中,在航天器本体姿态或参考姿态变化并不大,控制量没有明显改变时,并不需要实时更新控制信号。出于降低控制信号更新频率,节省通信资源的目的,如附图2所示,本发明采用事件触发的机制决定控制信号是否更新。设计比例积分形式的事件触发条件如下:
Figure BDA0003138866750000073
tk+1,i=min{tk,ik,i,tk,i+Ti},i=1,2,3
其中inf表示下确界,θi为第i维信号两次触发之间满足条件的可能时间间隔,θk,i表示第k时刻到下次触发满足条件的最小时间间隔,比例系数δ1,i、积分系数δ2,i、高增益参数εi和阈值常数Mi均为设计参数,且都为正值,Ti为预设的两次触发之间的最大允许间隔,tk,i和tk+1,i分别表示第k次的触发时刻和第k+1次触发时刻,σi(·)为反馈控制器输出的第i维信号的采样误差,定义为:
Figure BDA0003138866750000074
其中
Figure BDA0003138866750000075
表示tk,i时刻的采样值经零阶保持器保持后的信号。
需要说明的是,上述事件触发条件是通过连续监测反馈控制器输出信号,计算控制信号的采样误差,当控制信号的采样误差与其积分值的差满足预先设定的阈值时,才更新控制信号。如附图2所示,在航天器姿态控制执行机构端,通过零阶保持器将离散的采样的信号还原为连续信号,以满足姿态控制要求。此外,本发明所设计的比例积分事件触发条件不存在芝诺现象,且两次采样间隔不会超过在触发条件中预先设定的Ti
如图3所示,为本发明实施例的仿真结果。图3展示了角速度ω2的跟踪控制效果,可以看到在本发明所提出的控制方法下,卫星角速度跟踪误差收敛,并且控制信号只在触发时刻进行更新,有效降低了信号传输速率。
本发明针对刚性航天器,在存在外部干扰与内部惯量矩阵不确定性的情况下,提出了比例积分型的事件触发姿态控制方法,在保证姿态控制性能基础上,降低了姿态跟踪控制系统的通信消耗。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种比例积分型事件触发航天器姿态控制方法,其特征在于,包括以下步骤:
步骤一、基于四元数的卫星姿态动力学与运动学模型,构建卫星姿态跟踪误差系统模型;
步骤二、根据卫星姿态跟踪误差系统模型,确定姿态跟踪误差系统扩张状态,设计扩张状态观测器;
步骤三、利用扩张状态观测器得到系统观测状态,设计基于系统观测状态的反馈控制器;
步骤四、基于反馈控制器输出信号的采样误差,设计比例积分型事件触发策略用于决定是否更新航天器姿态的控制信号;
比例积分型事件触发条件如下:
Figure FDA0003630080380000011
tk+1,i=min{tk,ik,i,tk,i+Ti},i=1,…,n
其中,inf表示下确界,n为控制器输出信号的维度,θi为第i维信号两次触发之间满足条件的可能时间间隔,θk,i表示第k时刻到下次触发满足条件的最小时间间隔,比例系数δ1,i、积分系数δ2,i、高增益参数εi和阈值常数Mi均为设定参数,且都为正值,Ti为预设的两次触发之间的最大允许间隔,tk,i和tk+1,i分别表示第k次的触发时刻和第k+1次触发时刻,σi(·)为反馈控制器输出的第i维信号的采样误差,定义为:
Figure FDA0003630080380000012
其中
Figure FDA0003630080380000013
表示反馈控制器输出的第i维信号,
Figure FDA0003630080380000014
表示tk,i时刻的反馈控制器输出采样值经过零阶保持器保持后的信号。
2.如权利要求1所述的一种比例积分型事件触发航天器姿态控制方法,其特征在于,若比例积分型事件触发条件满足时,采样器对反馈控制器的输出信号进行采样,若比例积分型事件触发条件未满足,则航天器姿态的控制信号通过零阶保持器对上次触发传递而来的信号进行保持得到。
3.如权利要求1所述的一种比例积分型事件触发航天器姿态控制方法,其特征在于,所述比例积分型事件触发条件通过对比反馈控制器输出信号采样误差和输出信号采样误差关于时间的积分值的大小,判断何时进行触发,并更新控制信号,利用经零阶保持器保持后的控制信号实现对卫星姿态的控制,使得姿态跟踪误差收敛。
CN202110726476.3A 2021-06-29 2021-06-29 一种比例积分型事件触发航天器姿态控制方法 Active CN113608540B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110726476.3A CN113608540B (zh) 2021-06-29 2021-06-29 一种比例积分型事件触发航天器姿态控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110726476.3A CN113608540B (zh) 2021-06-29 2021-06-29 一种比例积分型事件触发航天器姿态控制方法

Publications (2)

Publication Number Publication Date
CN113608540A CN113608540A (zh) 2021-11-05
CN113608540B true CN113608540B (zh) 2022-07-05

Family

ID=78336921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110726476.3A Active CN113608540B (zh) 2021-06-29 2021-06-29 一种比例积分型事件触发航天器姿态控制方法

Country Status (1)

Country Link
CN (1) CN113608540B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114326393B (zh) * 2021-12-16 2024-03-12 安徽大学 非完整小车轨迹跟踪控制方法及装置
CN114942592B (zh) * 2022-06-01 2024-07-02 曲阜师范大学 一种柔性航天器的双端事件触发自适应抗干扰控制方法
CN115509245B (zh) * 2022-09-22 2024-08-27 哈尔滨工业大学 一种基于比例变换的连续力矩航天器姿态跟踪控制方法
CN115857329B (zh) * 2022-11-07 2024-03-08 西北工业大学深圳研究院 基于事件触发的绳系航天器复合体姿态稳定控制方法
CN117075635B (zh) * 2023-09-21 2024-07-05 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) 基于触发策略的多航天器控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108445748A (zh) * 2018-02-02 2018-08-24 北京航空航天大学 一种基于事件触发机制的自适应航天器姿态跟踪控制方法
CN109002058A (zh) * 2018-09-10 2018-12-14 北京航空航天大学 基于事件触发的航天器编队飞行相对位置协同控制方法
CN109828597A (zh) * 2019-03-02 2019-05-31 哈尔滨理工大学 一种基于事件触发的航天器姿态控制系统
CN112180734A (zh) * 2020-10-15 2021-01-05 杭州电子科技大学 一种基于分布式自适应事件触发的多智能体一致性方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7117045B2 (en) * 2001-09-08 2006-10-03 Colorado State University Research Foundation Combined proportional plus integral (PI) and neural network (nN) controller

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108445748A (zh) * 2018-02-02 2018-08-24 北京航空航天大学 一种基于事件触发机制的自适应航天器姿态跟踪控制方法
CN109002058A (zh) * 2018-09-10 2018-12-14 北京航空航天大学 基于事件触发的航天器编队飞行相对位置协同控制方法
CN109828597A (zh) * 2019-03-02 2019-05-31 哈尔滨理工大学 一种基于事件触发的航天器姿态控制系统
CN112180734A (zh) * 2020-10-15 2021-01-05 杭州电子科技大学 一种基于分布式自适应事件触发的多智能体一致性方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于事件触发的航天器姿态自适应容错控制;张嘉芮等;《飞控与探测》;20200331(第02期);第17-25页 *

Also Published As

Publication number Publication date
CN113608540A (zh) 2021-11-05

Similar Documents

Publication Publication Date Title
CN113608540B (zh) 一种比例积分型事件触发航天器姿态控制方法
Zou et al. Fixed-time attitude tracking control for rigid spacecraft without angular velocity measurements
Liu et al. Adaptive neural network control of a flexible spacecraft subject to input nonlinearity and asymmetric output constraint
Roy et al. Artificial-delay adaptive control for underactuated Euler–Lagrange robotics
Zhang et al. Learning observer based and event-triggered control to spacecraft against actuator faults
Zhang et al. Neural-networks and event-based fault-tolerant control for spacecraft attitude stabilization
CN110347170B (zh) 可重复使用运载器再入段鲁棒容错制导控制系统及工作方法
Shin Adaptive dynamic surface control for a hypersonic aircraft using neural networks
CN105468007A (zh) 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法
Shahbazi et al. Robust constrained attitude control of spacecraft formation flying in the presence of disturbances
Zhang et al. Fault-tolerant attitude stabilization for spacecraft with low-frequency actuator updates: An integral-type event-triggered approach
Lian et al. Adaptive backstepping control based autopilot design for reentry vehicle
Xu et al. Fuzzy-logic-based robust attitude control of networked spacecraft via event-triggered mechanism
An et al. Differentiator based full-envelope adaptive control of air-breathing hypersonic vehicles
CN108646564A (zh) 一种基于事件触发的不确定再入飞行器模型的设计方法
CN111439392A (zh) 一种航天器编队位置协同控制方法
Xie et al. Adaptive neural network model-based event-triggered attitude tracking control for spacecraft
CN113126491A (zh) 基于t-s模糊干扰建模的抗干扰跟踪控制设计方法
Wang et al. Disturbance attenuation predictive optimal control for quad-rotor transporting unknown varying payload
CN116500891A (zh) 基于强化学习的航天器容错控制方法
Zhang et al. Semiglobal finite-time stabilization of saturated spacecraft rendezvous system by dynamic event-triggered and self-triggered control
CN117311377B (zh) 一种基于复合干扰分离估计的航天器精细姿态控制方法
Zhou et al. Event-triggered integral sliding mode controller for rigid spacecraft attitude tracking with angular velocity constraint
CN116923730B (zh) 一种具有自调节预设性能约束的航天器姿态主动容错控制方法
You et al. Extended state observer based robust position tracking control using nonlinear damping gain for quadrotors with external disturbance

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