[发明专利]一种基于固定时间微分器预测的抗时延高精度自抗扰姿态控制方法有效
申请号: | 202010080580.5 | 申请日: | 2020-02-05 |
公开(公告)号: | CN111198570B | 公开(公告)日: | 2021-06-01 |
发明(设计)人: | 韦常柱;浦甲伦;李源;许河川;米长伟 | 申请(专利权)人: | 哈尔滨工业大学 |
主分类号: | G05D1/08 | 分类号: | G05D1/08 |
代理公司: | 哈尔滨龙科专利代理有限公司 23206 | 代理人: | 高媛 |
地址: | 150001 黑龙*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 固定 时间 微分 预测 抗时延 高精度 姿态 控制 方法 | ||
1.一种基于固定时间微分器预测的抗时延高精度自抗扰姿态控制方法,其特征在于,包括以下步骤:
步骤一:设计固定时间收敛微分器并获取姿态变化速率观测值;所述设计固定时间收敛微分器并获取姿态变化速率观测值的过程为:
第一步:利用速率陀螺获取延时kh的飞行器姿态信息其中h表示采样周期,k表示延时的周期个数,表示俯仰角,为偏航角,为滚转角;
第二步:分别针对俯仰、偏航和滚转三个通道构建如下所示的三阶固定时间收敛微分器:
其中下标分别表示对应参数为俯仰、偏航和滚转通道的参数;yi(t)表示不同通道的测量值,即t表示时间,z1i表示不同通道姿态角的观测值;z2i表示不同通道姿态变化速率的观测值;z3i表示不同通道姿态变化加速度观测值,分别表示z1i、z2i和z3i关于时间t的导数,α1i,α2i,α3i表示大于1的高幂次项幂指数;β1i,β2i,β3i表示小于等于1的低幂次项幂指数,k1i和κ1i分别对应于固定时间收敛微分器高幂次项和低幂次项的增益系数;
第三步:设定俯仰、偏航、滚转通道的微分器初值为z1i(0)=y1i(0),z2i=0,z3i=0,y1i(0)表示的是俯仰角、偏航角和滚转角的初始状态;
第四步:采用欧拉积分方法,利用俯仰、偏航和滚转通道的姿态初值以及微分器形式(1)进行解算,在每个解算周期内获取和并更新z1i、z2i和z3i,则z2i即为期望获取的三通道姿态角变化速率;
步骤二:基于预测器预测实时飞行状态;所述预测器预测实时飞行状态的过程为:
其中和为预测后的实时姿态角变化速率,和为预测后的实时姿态角;
步骤三:构建飞行器三通道姿态误差跟踪模型;所述三通道姿态误差跟踪模型的构建过程为:
第一步:基于绕质心转动动力学和姿态运动学构建飞行器姿态变化模型为:
其中ω=[ωx,ωy,ωz]为飞行器体轴的滚转、偏航、俯仰姿态角速率矢量,是俯仰、偏航和滚转角,δ表示飞行器的副翼、方向舵和升降舵,d表示飞行器的姿态通道扰动量,B1表示三通道的输入矩阵,J表示飞行器的转动惯量矩阵,表示Ω关于时间的导数,表示ω关于时间的导数,R表示转换矩阵,
其中Jx、Jy、Jz分别为体轴转动惯量,Jxy、Jxz、Jyz分别为不同惯性积;
第二步:基于输入输出线性化方法对公式(3)进行线性化,得到相对阶为2的精确线性化模型为:
其中,表示姿态角矢量Ω关于时间的二次导数,表示转换矩阵R关于时间的导数;
第三步:设定一阶状态量x1=Ω,二阶状态量为则得出飞行器姿态变化模型为如下形式:
其中表示x1关于时间的导数,表示x2关于时间的导数,M表示姿态变化模型的总扰动量;
第四步:设定三通道姿态角指令Ωc和姿态角指令变化率则飞行器姿态误差跟踪模型为:
其中X1=x1-Ωc,X1表示俯仰、偏航、滚转三通道姿态角与姿态角指令的偏差;X2和均表示X1关于时间的导数;表示X2关于时间的导数,D表示姿态误差跟踪模型的总扰动量;
步骤四:构建包含扰动状态观测器和非线性反馈律的自抗扰控制系统,利用飞行器实时飞行状态,通过自抗扰控制系统生成实时气动舵的摆动指令;飞行器自抗扰控制系统的构建过程为:
第一步:同样采用微分器(1),对(7)中X1和X2进行观测,得到估计值设定为和需要说明的是,代入X1中的x1项为(2)中经过预测更新后的俯仰、偏航、滚转姿态角和
第二步:以俯仰通道为例,采用如下形式的非线性反馈律:
其中x1、x2即为和中俯仰通道对应状态量,r、d为待设计量,h即为采样步长,为俯仰通道控制量,即矢量RJ-1B1δ中与俯仰通道的对应分量;偏航通道和滚转通道采用相同形式的反馈控制律,最后三通道解算出的控制量矢量形式设定为
第三步:以俯仰通道为例,采用如下形式的扰动状态观测器:
其中x2即为中俯仰通道对应状态量,z1为对x2的观测值,U1为上一控制解算周期俯仰通道的总控制量,z2即为对扰动D在俯仰通道的分量的估计值;偏航通道和滚转通道采用相同形式的扰动观测器,最后三通道扰动观测器输出的扰动量设定为分别为z1和z2关于时间的导数;β02,β03表示扰动状态观测器的增益系数;
第四步:综合公式(8)和公式(9),推导可得三通道的最终控制量为:
其中δ即为三通道的副翼、方向舵和升降舵。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工业大学,未经哈尔滨工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010080580.5/1.html,转载请声明来源钻瓜专利网。