CN113296525A - 基于数据驱动的倾转三旋翼无人机非线性控制方法 - Google Patents
基于数据驱动的倾转三旋翼无人机非线性控制方法 Download PDFInfo
- Publication number
- CN113296525A CN113296525A CN202011334242.6A CN202011334242A CN113296525A CN 113296525 A CN113296525 A CN 113296525A CN 202011334242 A CN202011334242 A CN 202011334242A CN 113296525 A CN113296525 A CN 113296525A
- Authority
- CN
- China
- Prior art keywords
- unmanned aerial
- aerial vehicle
- rotor unmanned
- tilting
- controller
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000003044 adaptive effect Effects 0.000 claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000013461 design Methods 0.000 claims description 25
- 239000013598 vector Substances 0.000 claims description 20
- 238000011156 evaluation Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000012938 design process Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 230000006641 stabilisation Effects 0.000 abstract description 3
- 238000011105 stabilization Methods 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 description 19
- 238000004088 simulation Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 235000019606 astringent taste Nutrition 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000002459 sustained effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 206010034719 Personality change Diseases 0.000 description 1
- 241001334139 Rugopharynx chi Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Abstract
本发明涉及倾转三旋翼无人机非线性控制,为实现倾转三旋翼无人机姿态和高度方向的镇定控制,并对外部扰动具有一定的鲁棒性。为此,本发明采取的技术方案是,基于数据驱动的倾转三旋翼无人机非线性控制方法,在倾转三旋翼无人机动力学模型基础上,1)在动力学模型存在参数未知的前提下,使用基于数据驱动的无模型自适应控制方法保证系统状态跟踪收敛;2)融入离散自适应滑模控制增强倾转三旋翼无人机对系统外部扰动的鲁棒性,同时补偿1)中的模型逼近误差,从而实现倾转三旋翼无人机非线性控制。本发明主要应用于倾转三旋翼无人机非线性控制场合。
Description
技术领域
本发明涉及倾转三旋翼无人机非线性控制,具体涉及基于数据驱动的倾转三旋翼无人机 非线性控制方法。
背景技术
近年来,随着微机电系统和自动控制的快速发展,多旋翼无人机由于其具有良好的机动 性和稳定性,受到了越来越多的关注。目前已经在诸如高空摄影,桥梁勘测,农业灌溉等多 个领域得到了广泛的应用。而相比于其他的多旋翼无人机,倾转三旋翼无人机保留了垂直起 降的功能,其机械结构更加紧凑,这大大提高了系统的续航能力,更加适合在工业应用(会 议:In Proceedings of the 2005 IEEE International Conference onRobotics and Automation;著者: S.Salazar-Cruz and R.Lozano.;出版年月:2005年4月;文章题目:Stabilization and nonlinear control for a novel trirotor mini-aircraft;页码:2612–2617)。
倾转三旋翼无人机包含三个呈Y型分布的无刷直流电机,以及一个用以控制航向的伺服 舵机。相比于其他对称结构的多旋翼无人机,该种机型动力学模型更为复杂。目前针对倾转 三旋翼飞行控制,学者们提出了多种非线性控制方法,取得了一定了成果。如一些学者采用 反步滑模控制方法保证系统稳定性,并采用模糊控制补偿系统的模型中的不确定项,最终通 过数值仿真验证了算法的有效性(期刊:Microsystem Technologies;著者:S.Wang,J.Zhang, Q.Zhang,and C.Pei;出版年月:2017年12月;文章题目:Aninnovative fuzzy backstepping sliding mode controller for a tri-rotorunmanned aerial vehicle;页码:5621–5630)。也有一些学 者采用反馈线性化和H∞的控制方法处理系统的非线性和耦合性,同样采用数值仿真的方式验 证了算法的有效性(会议:In Proceedings of 2012 UKACC International Conference on Control; 著者:M.K.Mohamed and A.Lanzon.;出版年月:2012年9月;文章题目:Design and control ofnovel trirotor UAV;页码:304–309)。
另一方面,由系统的输入输出产生的数据可以实时反应系统的实际动态。可以采用基于 数据驱动的控制方式实时逼近系统的真实动态模型(期刊:IEEE Transactions onIndustrial Electronics;著者:Z.Hou,R.Chi,and H.Gao;出版年月:2017年5月;文章题目:An overview of dynamic-linearization based data-driven control andapplications;页码:4076–4090)。目前, 基于数据驱动的无模型自适应控制方法已经成功应用于无线通信系统,可植入心脏泵系统, 单旋翼直升机飞控系统等。考虑到倾转三旋翼无人机系统模型复杂且较难获得精确的动力学 模型,采用基于数据驱动的无模型自适应控制方法将会弥补上述不足。
综上,关于倾转三旋翼无人机控制的研究,目前研究人员已经取得了一定成就,但也存 在一些局限:1)较大程度依赖倾转三旋翼无人机精确的动力学模型,但是当考虑到实际飞行 过程时,一方面较难获得系统精确的动力学模型,另一方面控制器未考虑系统中的高频未建 模动态,这都将导致上述控制方法很难保证较好的飞行性能。2)多数非线性控制算法仅仅停 留在数值仿真层面,并未考虑到实际飞行过程中系统外部扰动。
发明内容
为克服现有技术的不足,本发明旨在实现倾转三旋翼无人机姿态和高度方向的镇定控制, 并对外部扰动具有一定的鲁棒性。为此,本发明采取的技术方案是,基于数据驱动的倾转三 旋翼无人机非线性控制方法,在倾转三旋翼无人机动力学模型基础上,1)在动力学模型存在 参数未知的前提下,使用基于数据驱动的无模型自适应控制方法保证系统状态跟踪收敛;2) 融入离散自适应滑模控制增强倾转三旋翼无人机对系统外部扰动的鲁棒性,同时补偿1)中 的模型逼近误差,从而实现倾转三旋翼无人机非线性控制。
具体步骤如下:
1)建立倾转三旋翼无人机动力学模型
为了更好的描述倾转三旋翼无人机的动力学和运动学模型,{I}={OI,XI,YI,ZI}和 {B}={OB,XB,YB,ZB}分别表示原点在地面和机体中心的惯性坐标系和机体坐标系, {Bd}={OBd,XBd,YBd,ZBd}表示期望位置的机体坐标系。定义为相 对于{I}表示在{B}中的单位四元数,为相对于{I}表示在{Bd}中的单位四元数,为相对于{B}表示在{Bd}中的单位四元数满足 下式:
其中q0,q0d,e0和qv,qvd,ev分别表示单位四元数的标量部分和矢量部分,定义R(q),Rd(q)和分别表示从{I}到{B},{I}到{Bd}和{Bd}到{B}的旋转矩阵,其中R(q),Rd(q) 和用四元数形式表示如下:
定义fi为第i个电机旋转所产生的升力,li为第i个电机到倾转三旋翼无人机中心的距离。 由于倾转三旋翼无人机存在机械对称,i=1,2,3,有l1=l2=l;存在一条虚线穿过2号电机和 无人机的中心,一条虚线穿过1号电机和2号电机,则δ表示两条虚线的夹角,同时用α(t)表 示舵机偏离XBOBZB平面的夹角,并定义顺时针为正方向,显然,α(t)的存在当导致倾转三 旋翼无人机在偏航方向进行旋转;
定义(·)B表示定义在{B}中的元素,(·)I表示定义在{I}中的元素。
则倾转三旋翼无人机相对于{B}中的四元素形式的运动学模型可以表示如下:
倾转三旋翼无人机的动力学模型表示为:
其中表示系统惯性矩阵,表示 外部未知有界扰动矢量,且Da为正常数,表示由三 个电机和舵机产生的控制输入矢量,z表示倾转三旋翼无人机在{I}中的高度值,φ,θ,ψ表示 滚转角,俯仰角和偏航角,m表示倾转三旋翼无人机的质量,g表示当地的重力加速度值,表示在z方向的总推力,表示在z方向的外部未知有界扰动,且||dτ||<Dz, Dz为正常数;
τB(t),Tm(t)和fi(t)的关系用下式表示:
其中ki,表示力矩和升力之间的常系数。联立(4)式和(5)式,倾转三旋翼无人机的动 力学模型改写成:
2)姿态控制器设计
2.1姿态系统的动态线性化
将r(t)对时间求导,并将(7)式和(10)式代入结果得:
将(12)式进行离散化处理得:
由(15)式得:
其中I表示单位矩阵,Ξ(k)定义如下:
然后,定义r(k+1)关于r(k)和U(k)的偏导数为:
则(16)式改写成:
Ξ(k)=[χ1(k),χ2(k)]△L(k). (20)
基于式(19),(20),倾转三旋翼无人机的姿态动力学模型写为:
△r(k+1)=Ω1(k)△r(k)+Ω2(k)△U(k). (21)
其中
△r(k+1)=Ω(k)△L(k). (23)
2.2无模型自适应控制器设计
控制性能依赖于实变的更新矩阵Ω(k),所以需要针对Ω(k)设计最优更新率,考虑到更新 速率与更新精度,最优评价函数设计如下:
同理,无模型自适应控制器强调最优控制器输出和系统输出,故最优评价函数设计如下:
J(U(k))=||rd(k+1)-r(k+1)||2+σ||U(k)-U(k+1)||2 (25)
其中rd(k+1)表示系统的等价期望输出信号,σ>0表示更新权重因子。
式(24)(25)分别对Ω(k)和U(k)求偏导,并另结果等于零,求其极小值得到:
2.3离散滑模控制器设计
为了增加系统对外部扰动的鲁棒性,同时补偿系统动态模型逼近误差,在控制器中需要 加入离散滑模控制器,其表达式如下:
其中T表示采用周期,q>0为可调参数,满足1-qT>0。设计离散滑模面如下:
s(k)=e(k)=rd(k)-r(k), (29)
结合式(21),得:
s(k+1)=rd(k+1)-Ω1(k)△r(k)-Ω2(k)△U(k)-r(k). (30)
将式(28)代入式(30)得
令离散滑模控制器为:
UASM(k)=△U(k). (32)
联立式(27)(32),式(7)中总的姿态控制器设计如下:
τB(k)=UMFAC(k)+ΓUASM(k) (33)
3)高度控制器设计
定义倾转三旋翼无人机高度通道跟踪误差ez(t)和其滤波误差rz(t)为:
其中ξ为正的增益系数,rz(t)对时间求导,并将式(7)代入结果得:
其中辅助函数Nz(t)定义如下:
与倾转三旋翼无人机姿态控制器设计流程类似,式(35)化为如下离散形式:
因此,高度控制器设计如下:
Tm(k)=UMZ(k)+ΓZUAZ(k) (40)
△Tm(k)=Tm(k)-Tm(k-1)。
本发明的特点及有益效果是:
本发明针对倾转三旋翼无人机的姿态通道和高度通道,建立了含有模型不确定性和外部 扰动的动力学模型。考虑到精确模型较难获取,以及外部未知扰动的存在,设计了一种基于 数据驱动的MFASC控制方法。该方法不依赖倾转三旋翼无人机的精确数学模型,对外部扰 动具有较强的鲁棒性。飞行实验证明该方法具有较好的控制精度,对外部扰动具有较强的鲁 棒性。
附图说明:
图1是本发明采用的倾转三旋翼无人机坐标系统示意图;
图2是本发明采用的倾转三旋翼无人机机体参数定义示意图;
图3是本发明所使用的倾转三旋翼无人机硬件在环仿真平台;
图4是采用控制方案后倾转三旋翼无人机受到外部扰动前后的三个姿态角的曲线图;
图5是采用控制方案后倾转三旋翼无人机受到外部扰动前后的四元数曲线图;
图6是采用控制方案后倾转三旋翼无人机受到外部扰动前后的控制器输出曲线图;
图7是采用控制方案后倾转三旋翼无人机受到外部扰动前后的高度状态输出曲线图;
图8是采用LQR控制器后倾转三旋翼无人机受到外部扰动前后的三个姿态角的曲线图;
图9是采用LQR控制器后倾转三旋翼无人机受到外部扰动前后的控制器输出曲线图;
图10是采用LQR控制器后倾转三旋翼无人机受到外部扰动前后的高度状态输出曲线图;
图11是采用控制方案后倾转三旋翼无人机受到外部持续阵风扰动前后的三个姿态角的 曲线图;
图12是采用控制方案后倾转三旋翼无人机受到外部持续阵风扰动前后的四元数曲线图;
图13是采用控制方案后倾转三旋翼无人机受到外部持续阵风扰动前后的控制器输出曲 线图;
图14是采用控制方案后倾转三旋翼无人机受到外部持续阵风扰动前后的高度状态输出 曲线图;
图15是采用LQR控制器后倾转三旋翼无人机受到外部持续阵风扰动前后的三个姿态角 的曲线图;
图16是采用LQR控制器后倾转三旋翼无人机受到外部持续阵风扰动前后的控制器输出 曲线图。
具体实施方式
本发明涉及倾转三旋翼无人机非线性控制。针对倾转三旋翼无人机系统精确模型较难获 得以及存在外界扰动的情况,提出一种无模型自适应滑模控制方法(MFASC)。该方法采用基 于数据驱动的无模型自适应控制方法实时估计三旋翼无人机的姿态和高度数学模型,并采用 离散滑模控制方法补偿估计误差和外部扰动,实现了倾转三旋翼无人机姿态和高度方向的镇 定控制,并对外部扰动具有一定的鲁棒性。具体涉及基于数据驱动的倾转三旋翼无人机的姿 态和高度控制方法。
本发明的目的在于克服现有技术的上述不足,减小倾转三旋翼无人机对于系统模型的依 赖,提高倾转三旋翼无人机在实际飞行过程中的稳定性和鲁棒性。具体来讲:1)在系统模型 存在不确定和参数未知的前提下,使用基于数据驱动的无模型自适应控制方法保证系统状态 跟踪收敛。2)融入离散自适应滑模控制增强倾转三旋翼无人机对系统外部扰动的鲁棒性,同 时补偿1)中的模型逼近误差。3)将提出的控制算法在实验平台上进行实际飞行验证,实验 结果显示提出的控制算法具有较高的控制精度以及对外部扰动的鲁棒性。本发明采取的技术 方案实施步骤如下:
1)建立倾转三旋翼无人机动力学模型
为了更好的描述倾转三旋翼无人机的动力学和运动学模型,首先定义两个坐标系统如图 1所示。图1中{I}={OI,XI,YI,ZI}和{B}={OB,XB,YB,ZB}分别表示原点在地面和机体中心的 惯性坐标系和机体坐标系。{Bd}={OBd,XBd,YBd,ZBd}表示期望位置的机体坐标系。定义为相对于{I}表示在{B}中的单位四元数,为相对于{I}表示在{Bd}中的单位四元数,为相对于{B}表示在{Bd}中的单位四元数满足下式:
其中q0,q0d,e0和qv,qvd,ev分别表示单位四元数的标量部分和矢量部分。定义R(q),Rd(q)和分别表示从{I}到{B},{I}到{Bd}和{Bd}到{B}的旋转矩阵,其中R(q),Rd(q) 和可以用四元数形式表示如下:
图2描述了倾转三旋翼无人机的动力学系数含义。其中定义fi(i=1,2,3)为第i个电机旋 转所产生的升力,li(i=1,2,3)为第i个电机到倾转三旋翼无人机中心的距离。由于倾转三旋翼 无人机存在机械对称,有l1=l2=l。假设存在一条虚线穿过2号电机和无人机的中心,一条 虚线穿过1号电机和2号电机,则δ表示两条虚线的夹角。同时用α(t)表示舵机偏离XBOBZB平面的夹角,并定义顺时针为正方向。显然,α(t)的存在当导致倾转三旋翼无人机在偏航方 向进行旋转。
定义(·)B表示定义在{B}中的元素,(·)I表示定义在{I}中的元素。
则倾转三旋翼无人机相对于{B}中的四元素形式的运动学模型可以表示如下(期刊:IEEE Transactions on Industrial Informatics;著者:B.Xian and W.Hao;出版年月: 2019年4月;文章题目:Nonlinear robust fault-tolerant control of the tilttrirotor uav under rear servo’s stuck fault;页码:2158–2166):
倾转三旋翼无人机的动力学模型可以表示为:
其中表示系统惯性矩阵,表示 外部未知有界扰动矢量,且Da为正常数。表示由三 个电机和舵机产生的控制输入矢量,z表示倾转三旋翼无人机在{I}中的高度值,φ,θ,ψ表示 滚转角,俯仰角和偏航角,m表示倾转三旋翼无人机的质量,g表示当地的重力加速度值, 表示在z方向的总推力,表示在z方向的外部未知有界扰动,且||dτ||<Dz,Dz为正常数。
τB(t),Tm(t)和fi(t)的关系可以用下式表示:
其中ki,(i=1,2,3)表示力矩和升力之间的常系数。联立(4)式和(5)式,倾转三旋翼无 人机的动力学模型可以改写成:
2)姿态控制器设计
2.1姿态系统的动态线性化
将r(t)对时间求导,并将(7)式和(10)式代入结果得:
将(12)式进行离散化处理得:
由(15)式可得:
其中I表示单位矩阵,Ξ(k)定义如下:
然后,定义r(k+1)关于r(k)和U(k)的偏导数为:
则(16)式可以改写成:
Ξ(k)=[χ1(k),χ2(k)]△L(k). (20)
基于式(19),(20),倾转三旋翼无人机的姿态动力学模型可以写为:
△r(k+1)=Ω1(k)△r(k)+Ω2(k)△U(k). (21)
其中
△r(k+1)=Ω(k)△L(k). (23)
从式(23)可以看出,倾转三旋翼无人机的姿态动力学模型满足广义的Lipschitz条件, 这也符合一般机械运动系统的系统特性。
2.2无模型自适应控制器设计
从式(23)可以看出,控制性能依赖于实变的更新矩阵Ω(k),所以需要针对Ω(k)设计最 优更新率,考虑到更新速率与更新精度,最优评价函数设计如下:
同理,无模型自适应控制器强调最优控制器输出和系统输出,故最优评价函数设计如下:
J(U(k))=||rd(k+1)-r(k+1)||2+σ||U(k)-U(k+1)||2 (25)
其中rd(k+1)表示系统的等价期望输出信号,σ>0表示更新权重因子。
式(24)(25)分别对Ω(k)和U(k)求偏导,并另结果等于零,求其极小值得到:
2.3离散滑模控制器设计
为了增加系统对外部扰动的鲁棒性,同时补偿系统动态模型逼近误差,在控制器中需要 加入离散滑模控制器,其表达式如下:
其中T表示采用周期,q>0为可调参数,满足1-qT>0。设计离散滑模面如下:
s(k)=e(k)=rd(k)-r(k), (29)
结合式(21),可得
s(k+1)=rd(k+1)-Ω1(k)△r(k)-Ω2(k)△U(k)-r(k). (30)
将式(28)代入式(30)得
令离散滑模控制器为:
UASM(k)=△U(k). (32)
联立式(27)(32),式(7)中总的姿态控制器可以设计如下:
τB(k)=UMFAC(k)+ΓUASM(k) (33)
3)高度控制器设计
定义倾转三旋翼无人机高度通道跟踪误差ez(t)和其滤波误差rz(t)为:
其中ξ为正的增益系数。将rz(t)对时间求导,并将式(7)代入结果可得:
其中辅助函数Nz(t)定义如下:
与倾转三旋翼无人机姿态控制器设计流程类似,式(35)可以化为如下离散形式:
因此,高度控制器设计如下:
Tm(k)=UMZ(k)+ΓZUAZ(k) (40)
△Tm(k)=Tm(k)-Tm(k-1)。
本发明采用的技术方案是:针对倾转三旋翼无人机的姿态通道和高度通道,建立了含有 模型不确定性和外部扰动的动力学模型。并考虑到精确模型较难获取,以及外部未知扰动的 存在,设计了一种基于数据驱动的MFASC控制方法。包括如下步骤:
1)建立倾转三旋翼无人机动力学模型
为了更好的描述倾转三旋翼无人机的动力学和运动学模型,首先定义两个坐标系统如图 1所示。图1中{I}={OI,XI,YI,ZI}和{B}={OB,XB,YB,ZB}分别表示原点在地面和机体中心的 惯性坐标系和机体坐标系。{Bd}={OBd,XBd,YBd,ZBd}表示期望位置的机体坐标系。定义 为相对于{I}表示在{B}中的单位四元数,为相对于{I}表示在{Bd}中的单位四元数,为相对于{B}表示在{Bd}中的单位四元数满足下式:
其中q0,q0d,e0和qv,qvd,ev分别表示单位四元数的标量部分和矢量部分。定义R(q),Rd(q)和分别表示从{I}到{B},{I}到{Bd}和{Bd}到{B}的旋转矩阵,其中R(q),Rd(q) 和可以用四元数形式表示如下:
图2描述了倾转三旋翼无人机的动力学系数含义。其中定义fi(i=1,2,3)为第i个电机旋 转所产生的升力,li(i=1,2,3)为第i个电机到倾转三旋翼无人机中心的距离。由于倾转三旋翼 无人机存在机械对称,有l1=l2=l。假设存在一条虚线穿过2号电机和无人机的中心,一条 虚线穿过1号电机和2号电机,则δ表示两条虚线的夹角。同时用α(t)表示舵机偏离XBOBZB平面的夹角,并定义顺时针为正方向。显然,α(t)的存在当导致倾转三旋翼无人机在偏航方 向进行旋转。
定义(·)B表示定义在{B}中的元素,(·)I表示定义在{I}中的元素。
则倾转三旋翼无人机相对于{B}中的四元素形式的运动学模型可以表示如下(期刊:IEEE Transactions on Industrial Informatics;著者:B.Xian and W.Hao;出版年月: 2019年4月;文章题目:Nonlinear robust fault-tolerant control of the tilttrirotor uav under rear servo’s stuck fault;页码:2158–2166):
倾转三旋翼无人机的动力学模型可以表示为:
其中表示系统惯性矩阵,表示 外部未知有界扰动矢量,且Da为正常数。表示由三 个电机和舵机产生的控制输入矢量,z表示倾转三旋翼无人机在{I}中的高度值,φ,θ,ψ表示 滚转角,俯仰角和偏航角,m表示倾转三旋翼无人机的质量,g表示当地的重力加速度值, 表示在z方向的总推力,表示在z方向的外部未知有界扰动,且||dτ||<Dz, Dz为正常数。
τB(t),Tm(t)和fi(t)的关系可以用下式表示:
其中ki,(i=1,2,3)表示力矩和升力之间的常系数。联立(4)式和(5)式,倾转三旋翼无 人机的动力学模型可以改写成:
2)姿态控制器设计
2.1姿态系统的动态线性化
将r(t)对时间求导,并将(7)式和(10)式代入结果得:
将(12)式进行离散化处理得:
由(15)式可得:
其中I表示单位矩阵,Ξ(k)定义如下:
然后,定义r(k+1)关于r(k)和U(k)的偏导数为:
则(16)式可以改写成:
Ξ(k)=[χ1(k),χ2(k)]△L(k). (20)
基于式(19),(20),倾转三旋翼无人机的姿态动力学模型可以写为:
△r(k+1)=Ω1(k)△r(k)+Ω2(k)△U(k). (21)
其中
△r(k+1)=Ω(k)△L(k). (23)
从式(23)可以看出,倾转三旋翼无人机的姿态动力学模型满足广义的Lipschitz条件, 这也符合一般机械运动系统的系统特性。
2.2无模型自适应控制器设计
从式(23)可以看出,控制性能依赖于实变的更新矩阵Ω(k),所以需要针对Ω(k)设计最 优更新率,考虑到更新速率与更新精度,最优评价函数设计如下:
同理,无模型自适应控制器强调最优控制器输出和系统输出,故最优评价函数设计如下:
J(U(k))=||rd(k+1)-r(k+1)||2+σ||U(k)-U(k+1)||2 (25)
其中rd(k+1)表示系统的等价期望输出信号,σ>0表示更新权重因子。
式(24)(25)分别对Ω(k)和U(k)求偏导,并另结果等于零,求其极小值得到:
2.3离散滑模控制器设计
为了增加系统对外部扰动的鲁棒性,同时补偿系统动态模型逼近误差,在控制器中需要 加入离散滑模控制器,其表达式如下:
其中T表示采用周期,q>0为可调参数,满足1-qT>0。设计离散滑模面如下:
s(k)=e(k)=rd(k)-r(k), (29)
结合式(21),可得
s(k+1)=rd(k+1)-Ω1(k)△r(k)-Ω2(k)△U(k)-r(k). (30)
将式(28)代入式(30)得
令离散滑模控制器为:
UASM(k)=△U(k). (32)
联立式(27)(32),式(7)中总的姿态控制器可以设计如下:
τB(k)=UMFAC(k)+ΓUASM(k) (33)
3)高度控制器设计
定义倾转三旋翼无人机高度通道跟踪误差ez(t)和其滤波误差rz(t)为:
其中ξ为正的增益系数。将rz(t)对时间求导,并将式(7)代入结果可得:
其中辅助函数Nz(t)定义如下:
与倾转三旋翼无人机姿态控制器设计流程类似,式(35)可以化为如下离散形式:
因此,高度控制器设计如下:
Tm(k)=UMZ(k)+ΓZUAZ(k) (40)
△Tm(k)=Tm(k)-Tm(k-1)。
下面给出具体实施实例:
一、实验平台介绍
本发明利用图3所示的实验平台验证所设计的无模型自适应控制滑模控制器的效果。本 实验平台为倾转三旋翼无人机硬件在环仿真平台。该平台采用真实的倾转三旋翼无人机为被 控对象,并在倾转三旋翼无人机上加载了自主设计的姿态传感器采集电路板,由此可以得到 真实且直观的倾转三旋翼无人机姿态控制效果,也使结果更贴近实际飞行情况。而倾转三旋 翼无人机的高度通道则采用仿真模拟的形式完成,在高度通道加入白噪声信号模拟真实环境 中的传感器信号。该平台通过网络接口与上位机Matlab/Simulink模块进行通信,采用X86 结构下位机完成数据的实时采集,显示,同时完成控制算法的实现。具有方便调试,结果直 观,真实倾转三旋翼无人机模拟的特点。
二、飞行实验结果
为验证本发明提出的基于数据驱动的倾转三旋翼无人机非线性控制方法的有效性和可实 施性,在上述实验平台上进行半实物飞行验证实验。实验共分为两部分,分别为外部人为扰 动实验和外部抗风实验,同时采用传统LQR控制器作为对比。
实验过程中,倾转三旋翼无人机的系统参数为m=0.5kg,g=9.8m/s2。本发明所提出的 控制器参数选择为:
λ=diag{400,314,1500},ρ=diag{0.42,0.49,0.07},ξ=2.5,ρz=0.8,σ=σz=1,
ε=εz=0.01,q=qz=0.01,Γ=0.01,Γz=0.5.期望的单位四元数为qd=[1,0,0,0]T,期望的角速度为ωd=[0,0,0]T。针对虚拟的高度通道,给定的参考轨迹zd(t)为 从地面10秒钟上升到离地面10m高度。
为得到对应的LQR控制器参数,将式(5)中倾转三旋翼无人机动力学模型通过Matlab/Simulink工具箱进行线性化处理。LQR控制器形式为u(t)=-Kx(t),其中可调参数K由Matlab工具箱自动计算,并通过实际飞行实验微调,最终为:
实验1:抗扰对比实验
本实验中,首先倾转三旋翼无人机处于悬停状态,然后在滚转,俯仰和偏航方向施加外 部人为扰动。实验结果如图4-7所示。同样的LQR实验结果如图8-10所示。
从图4可以看出,外部扰动加入为25秒到65秒之间,扰动使得最大的姿态角误差为: 滚转和俯仰通道10度,偏航通道20度。然后,倾转三旋翼无人机能够迅速返回到平衡点, 并保持姿态角度误差1度以内。从图8可以看出,在相同情况下,滚转和俯仰方向的最大姿 态角误差只能达到4度,偏航方向最大为6度。当返回到平衡位置后,稳态误差接近1度。
从图7和图10可以看出,在高度通道两种控制器均达到了较好的轨迹跟踪效果。其中本 发明提出的控制算法具有较小的超调量和更高的跟踪精度。
实验2:外部风扰对比实验
本实验中,倾转三旋翼无人机在外部阵风中飞行,阵风速度为6m/s-7m/s,采用本发 明提出的控制算法得到的实验结果如图11-14所示,相应的采用LQR控制器得到的实验结果 如图15-17所示。
从图11中可以看出,外部阵风扰动在第23秒时加入,引起在滚转,俯仰和偏航方向的 姿态突变,然后在10秒之内,倾转三旋翼无人机姿态回到1度以内。而从图15可以看出,阵风加入后,姿态角度无法回复到1度以内,这将导致倾转三旋翼无人机较难在空间保持到一点,其控制效果较差。
从图14可以看出,本发明提出的控制方法在高度通道的稳态误差为0.05m,而从17可 以看出,采用LQR控制方法得到的稳态误差为0.1m。
综上,本发明提出的基于数据驱动的倾转三旋翼无人机非线性控制方法具有较高的控制 精度和较好的可实施性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之 内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于数据驱动的倾转三旋翼无人机非线性控制方法,其特征是,在倾转三旋翼无人机动力学模型基础上,1)在动力学模型存在参数未知的前提下,使用基于数据驱动的无模型自适应控制方法保证系统状态跟踪收敛;2)融入离散自适应滑模控制增强倾转三旋翼无人机对系统外部扰动的鲁棒性,同时补偿1)中的模型逼近误差,从而实现倾转三旋翼无人机非线性控制。
2.如权利要求1所述的基于数据驱动的倾转三旋翼无人机非线性控制方法,其特征是,具体步骤如下:
1)建立倾转三旋翼无人机动力学模型
为了更好的描述倾转三旋翼无人机的动力学和运动学模型,{I}={OI,XI,YI,ZI}和{B}={OB,XB,YB,ZB}分别表示原点在地面和机体中心的惯性坐标系和机体坐标系,{Bd}={OBd,XBd,YBd,ZBd}表示期望位置的机体坐标系。定义为相对于{I}表示在{B}中的单位四元数,为相对于{I}表示在{Bd}中的单位四元数,为相对于{B}表示在{Bd}中的单位四元数满足下式:
其中q0,q0d,e0和qv,qvd,ev分别表示单位四元数的标量部分和矢量部分,定义R(q),Rd(q)和分别表示从{I}到{B},{I}到{Bd}和{Bd}到{B}的旋转矩阵,其中R(q),Rd(q)和用四元数形式表示如下:
定义fi为第i个电机旋转所产生的升力,li为第i个电机到倾转三旋翼无人机中心的距离。由于倾转三旋翼无人机存在机械对称,i=1,2,3,有l1=l2=l;存在一条虚线穿过2号电机和无人机的中心,一条虚线穿过1号电机和2号电机,则δ表示两条虚线的夹角,同时用α(t)表示舵机偏离XBOBZB平面的夹角,并定义顺时针为正方向,显然,α(t)的存在当导致倾转三旋翼无人机在偏航方向进行旋转;
定义(·)B表示定义在{B}中的元素,(·)I表示定义在{I}中的元素。
则倾转三旋翼无人机相对于{B}中的四元素形式的运动学模型可以表示如下:
倾转三旋翼无人机的动力学模型表示为:
其中表示系统惯性矩阵,表示外部未知有界扰动矢量,且Da为正常数,表示由三个电机和舵机产生的控制输入矢量,z表示倾转三旋翼无人机在{I}中的高度值,φ,θ,ψ表示滚转角,俯仰角和偏航角,m表示倾转三旋翼无人机的质量,g表示当地的重力加速度值,表示在z方向的总推力,表示在z方向的外部未知有界扰动,且||dτ||<Dz,Dz为正常数;
τB(t),Tm(t)和fi(t)的关系用下式表示:
其中ki,表示力矩和升力之间的常系数。联立(4)式和(5)式,倾转三旋翼无人机的动力学模型改写成:
2)姿态控制器设计
2.1姿态系统的动态线性化
将r(t)对时间求导,并将(7)式和(10)式代入结果得:
将(12)式进行离散化处理得:
由(15)式得:
其中I表示单位矩阵,Ξ(k)定义如下:
然后,定义r(k+1)关于r(k)和U(k)的偏导数为:
则(16)式改写成:
Ξ(k)=[χ1(k),χ2(k)]ΔL(k). (20)
基于式(19),(20),倾转三旋翼无人机的姿态动力学模型写为:
Δr(k+1)=Ω1(k)Δr(k)+Ω2(k)ΔU(k). (21)
其中
Δr(k+1)=Ω(k)ΔL(k). (23)
2.2无模型自适应控制器设计
控制性能依赖于实变的更新矩阵Ω(k),所以需要针对Ω(k)设计最优更新率,考虑到更新速率与更新精度,最优评价函数设计如下:
同理,无模型自适应控制器强调最优控制器输出和系统输出,故最优评价函数设计如下:
J(U(k))=||rd(k+1)-r(k+1)||2+σ||U(k)-U(k+1)||2 (25)
其中rd(k+1)表示系统的等价期望输出信号,σ>0表示更新权重因子。
式(24)(25)分别对Ω(k)和U(k)求偏导,并另结果等于零,求其极小值得到:
2.3离散滑模控制器设计
为了增加系统对外部扰动的鲁棒性,同时补偿系统动态模型逼近误差,在控制器中需要加入离散滑模控制器,其表达式如下:
其中T表示采用周期,q>0为可调参数,满足1-qT>0。设计离散滑模面如下:
s(k)=e(k)=rd(k)-r(k), (29)
结合式(21),得:
s(k+1)=rd(k+1)-Ω1(k)Δr(k)-Ω2(k)ΔU(k)-r(k). (30)
将式(28)代入式(30)得
令离散滑模控制器为:
UASM(k)=ΔU(k). (32)
联立式(27)(32),式(7)中总的姿态控制器设计如下:
τB(k)=UMFAC(k)+ΓUASM(k) (33)
3)高度控制器设计
定义倾转三旋翼无人机高度通道跟踪误差ez(t)和其滤波误差rz(t)为:
其中ξ为正的增益系数,rz(t)对时间求导,并将式(7)代入结果得:
其中辅助函数Nz(t)定义如下:
与倾转三旋翼无人机姿态控制器设计流程类似,式(35)化为如下离散形式:
因此,高度控制器设计如下:
Tm(k)=UMZ(k)+ΓZUAZ(k) (40)
ΔTm(k)=Tm(k)-Tm(k-1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011334242.6A CN113296525B (zh) | 2020-11-24 | 2020-11-24 | 基于数据驱动的倾转三旋翼无人机非线性控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011334242.6A CN113296525B (zh) | 2020-11-24 | 2020-11-24 | 基于数据驱动的倾转三旋翼无人机非线性控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113296525A true CN113296525A (zh) | 2021-08-24 |
CN113296525B CN113296525B (zh) | 2024-01-23 |
Family
ID=77318461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011334242.6A Active CN113296525B (zh) | 2020-11-24 | 2020-11-24 | 基于数据驱动的倾转三旋翼无人机非线性控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113296525B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114003053A (zh) * | 2021-11-02 | 2022-02-01 | 东南大学 | 一种基于ArduPilot的固定翼无人机自动驾驶自适应控制系统 |
CN114167720A (zh) * | 2021-10-18 | 2022-03-11 | 天津大学 | 基于观测器的倾转式三旋翼无人机轨迹跟踪控制方法 |
CN114488819A (zh) * | 2022-02-16 | 2022-05-13 | 哈尔滨逐宇航天科技有限责任公司 | 一种弱气动模型依赖的飞行器数据驱动控制方法 |
CN114995483A (zh) * | 2022-06-28 | 2022-09-02 | 四川大学 | 基于对偶四元数的倾转旋翼多旋翼无人机跟踪控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107357166A (zh) * | 2017-04-29 | 2017-11-17 | 天津大学 | 小型无人直升机的无模型自适应鲁棒控制方法 |
CN109991992A (zh) * | 2019-03-27 | 2019-07-09 | 天津大学 | 基于浸入不变集的小型无人直升机无模型自适应控制方法 |
-
2020
- 2020-11-24 CN CN202011334242.6A patent/CN113296525B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107357166A (zh) * | 2017-04-29 | 2017-11-17 | 天津大学 | 小型无人直升机的无模型自适应鲁棒控制方法 |
CN109991992A (zh) * | 2019-03-27 | 2019-07-09 | 天津大学 | 基于浸入不变集的小型无人直升机无模型自适应控制方法 |
Non-Patent Citations (6)
Title |
---|
BIN XIAN: "Nonlinear Robust Fault-Tolerant Control of the Tilt Trirotor UAV Under Rear Servo\'s Stuck Fault: Theory and Experiments", IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS, pages 2158 * |
GU XUN: "model free adaptive control design for a tilt trirotor unmanned aerial vehicle with quaternion feedback: Theory and implementation", INTERNATIONAL JOURNAL OF ADAPTIVE CONTROL AND SIGNAL PROCESSING, pages 122 - 137 * |
LI ZHI: "Chattering free sliding adaptive attitude contorl for quadrotor", 2016 IEEE CHINESE GUIDANCE, NAVIGATION AND CONTORL CONFERENCE, pages 707 - 712 * |
夏青元;徐锦法;张梁;: "倾转旋翼飞行器无模型自适应姿态控制", 系统工程与电子技术, no. 01 * |
徐晓旭: "基于浸入—不变集的三旋翼无人机高度系统 自适应控制器设计", 计算机测量与控制 * |
潘小龙: "小型无人直升机的无模型自适应鲁棒控制设计", 控制理论与应用, pages 1171 - 1178 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114167720A (zh) * | 2021-10-18 | 2022-03-11 | 天津大学 | 基于观测器的倾转式三旋翼无人机轨迹跟踪控制方法 |
CN114003053A (zh) * | 2021-11-02 | 2022-02-01 | 东南大学 | 一种基于ArduPilot的固定翼无人机自动驾驶自适应控制系统 |
CN114003053B (zh) * | 2021-11-02 | 2023-12-01 | 东南大学 | 一种基于ArduPilot的固定翼无人机自动驾驶自适应控制系统 |
CN114488819A (zh) * | 2022-02-16 | 2022-05-13 | 哈尔滨逐宇航天科技有限责任公司 | 一种弱气动模型依赖的飞行器数据驱动控制方法 |
CN114488819B (zh) * | 2022-02-16 | 2022-09-16 | 哈尔滨逐宇航天科技有限责任公司 | 一种弱气动模型依赖的飞行器数据驱动控制方法 |
CN114995483A (zh) * | 2022-06-28 | 2022-09-02 | 四川大学 | 基于对偶四元数的倾转旋翼多旋翼无人机跟踪控制方法 |
WO2024001583A1 (zh) * | 2022-06-28 | 2024-01-04 | 四川大学 | 基于对偶四元数的倾转旋翼多旋翼无人机跟踪控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113296525B (zh) | 2024-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113296525A (zh) | 基于数据驱动的倾转三旋翼无人机非线性控制方法 | |
Bisheban et al. | Geometric adaptive control with neural networks for a quadrotor in wind fields | |
Derafa et al. | Super twisting control algorithm for the attitude tracking of a four rotors UAV | |
Madani et al. | Adaptive control via backstepping technique and neural networks of a quadrotor helicopter | |
Cui et al. | Adaptive super-twisting trajectory tracking control for an unmanned aerial vehicle under gust winds | |
Gonzalez-Vazquez et al. | A new nonlinear PI/PID controller for quadrotor posture regulation | |
Koo et al. | Nonlinear control of a helicopter based unmanned aerial vehicle model | |
CN112947518B (zh) | 一种基于干扰观测器的四旋翼鲁棒姿态控制方法 | |
Liu et al. | Robust command filtered adaptive backstepping control for a quadrotor aircraft | |
CN112578805A (zh) | 一种旋翼飞行器的姿态控制方法 | |
Sun et al. | Nonlinear robust compensation method for trajectory tracking control of quadrotors | |
Wang et al. | Event driven model free control of quadrotor | |
CN115533915A (zh) | 一种不确定环境下空中作业机器人主动接触检测控制方法 | |
Kumar et al. | Quaternion feedback based autonomous control of a quadcopter uav with thrust vectoring rotors | |
Zhang et al. | Discrete-time adaptive neural tracking control and its experiments for quadrotor unmanned aerial vehicle systems | |
Vallejo-Alarcon | Robust backstepping control for highly demanding quadrotor flight | |
CN113961010A (zh) | 基于抗饱和有限时间自适应神经网络容错技术的四旋翼植保无人机跟踪控制方法 | |
Zhan et al. | Control system design and experiments of a quadrotor | |
Song et al. | Research on attitude control of quadrotor uav based on active disturbance rejection control | |
Meradi et al. | A predictive sliding mode control for quadrotor’s tracking trajectory subject to wind gusts and uncertainties | |
Guillén-Bonilla et al. | Vision-based nonlinear control of quadrotors using the photogrammetric technique | |
Bangura et al. | An open-source implementation of a unit quaternion based attitude and trajectory tracking for quadrotors | |
Abrougui et al. | Flight Controller Design Based on Sliding Mode Control for Quadcopter Waypoints Tracking | |
CN112034872B (zh) | 一种四旋翼无人机积分补偿确定性策略梯度控制方法 | |
CN111338369B (zh) | 一种基于非线性逆补偿的多旋翼飞行控制方法 |
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 |