CN109884895B - 基于饱和受限情况下的无人机自适应跟踪控制算法 - Google Patents

基于饱和受限情况下的无人机自适应跟踪控制算法 Download PDF

Info

Publication number
CN109884895B
CN109884895B CN201910179402.5A CN201910179402A CN109884895B CN 109884895 B CN109884895 B CN 109884895B CN 201910179402 A CN201910179402 A CN 201910179402A CN 109884895 B CN109884895 B CN 109884895B
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
control
subsystem
attitude
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
CN201910179402.5A
Other languages
English (en)
Other versions
CN109884895A (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.)
Nanjing University of Posts and Telecommunications
Original Assignee
Nanjing University of Posts and Telecommunications
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 Nanjing University of Posts and Telecommunications filed Critical Nanjing University of Posts and Telecommunications
Priority to CN201910179402.5A priority Critical patent/CN109884895B/zh
Publication of CN109884895A publication Critical patent/CN109884895A/zh
Application granted granted Critical
Publication of CN109884895B publication Critical patent/CN109884895B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于饱和受限情况下的无人机自适应跟踪控制算法,本发明方法包括以下步骤:首先,基于牛顿—欧拉模型建立非线性、欠驱动的四旋翼无人机模型,实现四旋翼无人机模型数学化处理;然后,基于输入饱和受限函数模型,建立系统框架和选取输入受限函数;其次,基于滑模变结构控制方法,实现四旋翼无人机位置子系统的控制处理;最后,基于自适应跟踪控制方法,引入自抗扰控制处理无人机不可控的扰动,实现四旋翼无人机姿态子系统的控制处理。本文所提方法能在无人机输入饱和受限、系统扰动不知且不可控等复杂情况下,采用双闭环的控制模型,实现无人机位置和姿态的自适应跟踪控制,最终使无人机飞行状态达到稳定。

Description

基于饱和受限情况下的无人机自适应跟踪控制算法
技术领域
本发明涉及一种无人机控制算法,具体的说是一种四旋翼无人机跟踪控制算法,属于工业自动控制技术领域。
背景技术
四旋翼无人机是一种具有六个自由度的垂直起降飞机,可以突破跑道限制并自由悬停。与其他无人机平台相比,它具有结构简单,稳定,超强的灵活性。四旋翼无人机广泛应用于战术侦察,航空摄影,电磁干扰,短途货运等领域。对四旋翼无人机的研究具有重要的现实意义,已成为国内外研究的热点之一。
目前,大量学者将四旋翼无人机作为复杂的非线性,多变量,欠驱动控制系统来研究。并且,已有许多关于四旋翼无人机跟踪控制的论文发表出来。近年来,针对四旋翼无人机系统的位置和姿态控制的研究,不同学者分别提出了自抗扰(ADRC),滑模控制,预测控制,反步法控制,自适应控制以及鲁棒控制等控制方法。然而,四旋翼无人机具有非线性,高阶,强耦合和欠驱动的特点,导致其稳定性分析和控制器设计难以实现。
在四旋翼无人机控制系统中,滑模控制和自适应跟踪控制被广泛采用。本发明提出了一种基于滑模控制和自适应跟踪控制的双闭环控制器,用于控制四旋翼无人机的位置和姿态。在本发明所提出的双闭环结构中,滑模控制用于控制外环的位置子系统,自适应跟踪控制用于控制内环的姿态子系统。
然而,在实际飞行中,无人机必须进行大量的机动性飞行,并且所需的控制量通常超过执行器可提供的最大输出(执行器的饱和特性)。这样的限制往往降低实际的飞行控制性能,在严重的情况下,它可能会导致整个系统的不稳定。因此,在飞行控制中需要考虑致动器饱和的影响。本发明提出了一类具有输入约束的不确定MIMO非线性系统的自适应跟踪控制。引入辅助设计系统来分析输入约束的影响,并将其状态用于自适应跟踪控制设计。
本发明基于全局稳定轨迹跟踪控制算法,提出了一种欠驱动四旋翼无人机控制模型。将四旋翼无人机系统将分为欠驱动位置子系统和姿态子系统。对于这样的双闭环控制系统,位置子系统控制器采用滑模控制方法设计。考虑到输入饱和约束的影响,姿态子系统采用自抗扰和自适应跟踪控制方法设计。通过这样的控制器设计,可以顺利跟踪四旋翼无人机系统的位置和姿态角。并且,双闭环控制系统也将具有全局稳定性。
发明内容
本发明为解决四旋翼无人机系统在输入饱和受限的控制情况下,提出了一种基于饱和受限情况下的无人机自适应跟踪控制算法,采用双闭环的控制框架,对外环采用滑模控制方法控制位置子系统,对内外采用自适应控制方法控制姿态子系统,从而达到了四旋翼无人机系统在输入饱和受限的情况下所期望要求的效果。
本发明提供一种基于饱和受限情况下的无人机自适应跟踪控制算法,包括以下步骤:
步骤A、在考虑扰动存在的情况下,建立四旋翼无人机的二阶系统数学模型;
步骤B、基于四旋翼无人机的二阶系统模型,对于外环位置子系统,提出滑模控制器对无人机系统进行控制,以实现对无人机位置状态的控制;
步骤C、基于四旋翼无人机的二阶系统模型,对于内环姿态子系统,提出自适应控制策略对无人机系统进行控制,以实现对无人机姿态状态的控制。
作为本发明的进一步限定,步骤A具体为:所述二阶系统数学模型为
Figure BDA0001990660910000031
其中,状态x,y,z表示四旋翼的位置,状态θ,φ,ψ由三个姿态角组成,g是重力加速度,Ii(i=1,2,3)代表惯性矩,Ki(i=1,2…6)代表阻力系数,di(i=1,2…6)是系统扰动,滚动角φ和俯仰角θ的范围(-2/π<φ、θ<2/π),偏航角ψ的范围也限于(-π<ψ<π),ui(i=1,2,3,4)表示四旋翼的四个控制输入量。
将公式(1)所表示的四旋翼无人机数学模型拆分成位置和姿态两部分,以便于分布式处理:
Figure BDA0001990660910000032
Figure BDA0001990660910000033
将四旋翼无人机模型写成矩阵形式,因此,公式(2)可以用以下形式书写:
Figure BDA0001990660910000034
其中,ξ=[x y z]T是四旋翼无人机在惯性坐标系中的质心位置,
Figure BDA0001990660910000035
是垂直单位矢量,g是重力加速度,dF=[md1 md2 md3]T,m是四旋翼无人机的质量;
Figure BDA0001990660910000036
是旋转矩阵,其中c·=cos(·)、s·=sin(·)。
作为本发明的进一步限定,步骤B具体为:所述滑模控制器为:
定义一个滑模控制函数
Figure BDA0001990660910000041
其中,λ1是一个正数并满足Hurwitz条件,并且s1(t)表示平动的滑模面。
选择指数收敛的控制律,滑模面的导数应如下
Figure BDA0001990660910000042
其中ε1≥DF/m和k1是正常数,符号函数定义如下
Figure BDA0001990660910000043
考虑公式(5),公式(6)的导数是
Figure BDA0001990660910000044
因此,对比公式(5)和公式(8),位置子系统的滑模控制律可以设计如下
Figure BDA0001990660910000045
作为本发明的进一步限定,步骤C具体为:所述自适应控制器为:
假设预期跟踪的参考姿态是δd=[θd φd ψd]T并且虚拟控制律是αi=[α1 α2 α3]T;定义错误变量e1=δ-δd和e2=δ1i,虚拟控制定律αi可以定义如下
Figure BDA0001990660910000046
考虑公式(3)并求导,我们得到可用位置子系统的误差方程
Figure BDA0001990660910000047
其中δ=[θ φ ψ]T为无人机姿态角,U=[u2 u3 u4]T,
Figure BDA0001990660910000048
I=[I1 I2 I3]T代表惯性矩,dF=[d4 d5 d6]T代表阻力系数。
考虑Lyapunov函数
Figure BDA0001990660910000051
Vi求导可得
Figure BDA0001990660910000052
引入αi,可得
Figure BDA0001990660910000053
为方便输入饱和受限影响的分析,给出了辅助设计系统
Figure BDA0001990660910000054
其中
Figure BDA0001990660910000055
η∈R3是辅助设计系统的状态,设计参数是正常数,应根据跟踪性能的要求选择适当的值;
考虑输入饱和的影响,请选择以下控制律
Figure BDA0001990660910000056
作为本发明的进一步限定,其对于四旋翼无人机在获得虚拟输入之后,还得计算实际升程控制输入并将其传递到姿态子系统中间命令信号,具体包括以下步骤:
将虚拟控制输入写为向量Uξ=[ux,uy,uz]T,将中间命令信号设置为Θd=[φddd]T
考虑
Figure BDA0001990660910000057
Uξ=[ux,uy,uz]T
Figure BDA0001990660910000058
可以得到
Figure BDA0001990660910000061
进一步
Figure BDA0001990660910000062
变为矩阵形式
Figure BDA0001990660910000063
使用Θd代替姿态子系统的中间命令信号Θ,获得
Figure BDA0001990660910000064
因为
Figure BDA0001990660910000065
可得
Figure BDA0001990660910000066
滚转角指令信号是
Figure BDA0001990660910000067
俯仰角指令信号是
Figure BDA0001990660910000068
偏航角ψd可用于跟踪给定指令信号的任意偏航角。
作为本发明的进一步限定,其对于四旋翼无人机的稳定性采用如下方式验证:是利用Lyapunov稳定性理论,设计Lyapunov函数证明在四旋翼无人机位置子系统和姿态子系统的稳定性;具体包括以下步骤:
步骤a:对于采用滑模控制策略的位置子系统,考虑Lyapunov函数
Figure BDA0001990660910000071
V1的导数如下:
Figure BDA0001990660910000072
由于ε1≥DF/m,可得
Figure BDA0001990660910000073
从上面的分析我们可以看出s1是指数收敛,即eξ指数收敛。类似地,位置误差子系统是指数稳定的。
步骤b:对于采用自适应控制策略的姿态子系统,考虑Lyapunov函数
Figure BDA0001990660910000074
根据公式(27),V的导数
Figure BDA0001990660910000075
将公式(15)带入公式(28),可得
Figure BDA0001990660910000081
为了确保系统的稳定性,选择适当的参数k22来确保k22-I3×3≥0。
作为本发明的进一步限定,在外界扰动和内在误差无法估计的时候,提出非线性ES0控制策略,对内部不可建模和外部不可估计的扰动进行估计和补偿:ESO可用于姿态子系统,以估计和补偿外部环境的总干扰,对于这种情况,系统扩张状态方程是
Figure BDA0001990660910000082
状态观察器,即ESO是为了扩张状态而建立的
Figure BDA0001990660910000083
选择适当的参数βi(i=1,2,3),d4可以被估计和补偿;
类似地,φ和ψ也可以使用扩展状态观察器来估计误差,这样也可以通过ESO进行补偿;
通过将自适应跟踪控制方法与ESO相结合,可以得到新的姿态角控制律,如下所示
Figure BDA0001990660910000084
其中di是对系统总扰动的估计。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)本发明考虑了四旋翼无人机的输入饱和受限的情况,因为无人机的转速不可能无限大,转速的限定决定了无人机的输入是受限的,所以考虑四旋翼无人机的输入饱和受限更具现实意义,符合实际的飞行模型;
(2)本发明选用了双闭环的控制框架,在内环子系统采用滑模控制策略实现四旋翼无人机的位置控制,在外环子系统采用自适应控制策略实现四旋翼无人机的姿态控制。双闭环的控制框架,能够对四旋翼无人机的位置和姿态分别控制,实现了不同控制策略的具体应用,对比传统的控制方法,更有创新性和应用前景;
(3)本发明针对四旋翼无人机的姿态控制提出了自适应和非线性ESO控制方法,对不可建模和无法估计的扰动进行控制,并通过ESO进行补偿,从而使无人机的姿态控制更加稳定,无位置控制形成双闭环,最终达到无人机的稳定性控制要求;
(4)本发明对无人机的输入饱和受限,位置和姿态一共六个自由度提出了相对应的滑模控制和自适应控制,并在姿态子系统中选取了虚拟控制器和辅助方程,成功的实现了对四旋翼无人机的控制,最终通过仿真验证本发明的正确性,具有一定的创新性。
附图说明
图1为本发明中控制算法示意图一。
图2为本发明中控制算法示意图二。
图3为本发明中位置跟踪结果图。
图4为本发明中带有输入饱和控制器的姿态角稳定性结果图。
图5为本发明中没有输入饱和控制器对姿态角稳定性结果图。
图6为本发明中控制输入结果图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
下面结合附图对本发明的技术方案进行详细说明:
本发明的目的是提供一种基于饱和受限情况下的无人机自适应跟踪控制算法,其实现思路为:首先,基于牛顿—欧拉模型建立非线性、欠驱动的四旋翼无人机模型,实现四旋翼无人机数学化处理;然后,基于输入饱和受限函数模型,建立系统框架和选取输入受限函数;其次,基于滑模变结构控制方法,实现四旋翼无人机位置子系统的控制处理;最后,基于自适应跟踪控制方法,引入自抗扰控制处理无人机不可控的扰动,实现四旋翼无人机姿态子系统的控制处理。本文所提方法能在无人机输入饱和受限、系统扰动不知且不可控等复杂情况下,采用双闭环的控制模型,实现无人机位置和姿态的自适应跟踪控制,最终使无人机飞行状态达到稳定。
如图1-2所示,一种基于饱和受限情况下的无人机自适应跟踪控制算法,其特征在于,包括以下步骤:
步骤A、在考虑扰动存在的情况下,建立四旋翼无人机的二阶系统数学模型:
Figure BDA0001990660910000101
其中,状态x,y,z表示四旋翼的位置,状态θ,φ,ψ由三个姿态角组成,g是重力加速度,Ii(i=1,2,3)代表惯性矩,Ki(i=1,2…6)代表阻力系数,di(i=1,2…6)是系统扰动,滚动角φ和俯仰角θ限于(-2/π<φ、θ<2/π),偏航角ψ也限于(-π<ψ<π),ui(i=1,2,3,4)表示四旋翼的四个控制输入量。
将公式(1)所表示的四旋翼无人机数学模型拆分成位置和姿态两部分,以便于分布式处理:
Figure BDA0001990660910000111
Figure BDA0001990660910000112
将四旋翼无人机模型编写成矩阵形式,因此,公式(2)可以用以下形式书写:
Figure BDA0001990660910000113
其中,ξ=[x y z]T是四旋翼无人机在惯性坐标系中的质心位置,
Figure BDA0001990660910000114
是垂直单位矢量,g是重力加速度,dF=[md1 md2 md3]T,m是四旋翼无人机的质量;
Figure BDA0001990660910000115
是旋转矩阵,其中c·=cos(·)、s·=sin(·);
步骤B、基于四旋翼无人机的二阶系统模型,对于外环位置子系统,提出滑模控制器对无人机系统进行控制,以实现对无人机位置状态的控制,所述滑模控制器为:
定义一个滑模控制函数
Figure BDA0001990660910000116
其中,λ1是一个正数并满足Hurwitz条件,并且s1(t)表示平动的滑模面。
选择指数收敛的控制律,滑模面的导数应如下
Figure BDA0001990660910000117
其中ε1≥DF/m和k1是正常数,符号函数定义如下
Figure BDA0001990660910000121
考虑公式(5),公式(6)的导数是
Figure BDA0001990660910000122
因此,对比公式(5)和公式(8),位置子系统的SMC定律可以设计如下
Figure BDA0001990660910000123
步骤C、基于四旋翼无人机的二阶系统模型,对于内环姿态子系统,提出自适应控制策略对无人机系统进行控制,以实现对无人机姿态状态的控制,所述自适应控制器为:
假设预期跟踪的参考姿态是δd=[θd φd ψd]T并且虚拟控制律是αi=[α1 α2 α3]T;定义错误变量e1=δ-δd和e2=δ1i,虚拟控制定律αi可以定义如下
Figure BDA0001990660910000124
考虑公式(3)并求导,我们得到可用位置子系统的误差方程
Figure BDA0001990660910000125
其中δ=[θ φ ψ]T为无人机姿态角,U=[u2 u3 u4]T,
Figure BDA0001990660910000126
I=[I1I2I3]T代表惯性矩,dF=[d4d5d6]T代表阻力系数。
考虑Lyapunov函数
Figure BDA0001990660910000127
Vi求导可得
Figure BDA0001990660910000128
引入αi,可得
Figure BDA0001990660910000131
为方便输入饱和受限影响的分析,给出了辅助设计系统
Figure BDA0001990660910000132
其中
Figure BDA0001990660910000133
η∈R3是辅助设计系统的状态,设计参数是正常数,应根据跟踪性能的要求选择适当的值;
考虑输入饱和的影响,请选择以下控制律
Figure BDA0001990660910000134
其对于四旋翼无人机在获得虚拟输入之后,还得计算实际升程控制输入并将其传递到姿态子系统中间命令信号,具体包括以下步骤:
将虚拟控制输入写为向量Uξ=[ux,uy,uz]T,将中间命令信号设置为Θd=[φddd]T
考虑
Figure BDA0001990660910000135
Uξ=[ux,uy,uz]T
Figure BDA0001990660910000136
可以得到
Figure BDA0001990660910000137
进一步
Figure BDA0001990660910000138
变为矩阵形式
Figure BDA0001990660910000141
使用Θd代替姿态子系统的中间命令信号Θ,获得
Figure BDA0001990660910000142
因为
Figure BDA0001990660910000143
可得
Figure BDA0001990660910000144
滚转角指令信号是
Figure BDA0001990660910000145
俯仰角指令信号是
Figure BDA0001990660910000146
偏航角ψd可用于跟踪给定指令信号的任意偏航角。
其对于四旋翼无人机的稳定性采用如下方式验证:是利用Lyapunov稳定性理论,设计Lyapunov函数证明在四旋翼无人机位置子系统和姿态子系统的稳定性;具体包括以下步骤:
步骤a:对于采用滑模控制策略的位置子系统,考虑Lyapunov函数
Figure BDA0001990660910000147
V1的导数如下:
Figure BDA0001990660910000151
由于ε1≥DF/m,可得
Figure BDA0001990660910000152
从上面的分析我们可以看出s1是指数收敛,即eξ指数收敛。类似地,位置误差子系统是指数稳定的。
步骤b:对于采用自适应控制策略的姿态子系统,考虑Lyapunov函数
Figure BDA0001990660910000153
根据公式(27),V的导数
Figure BDA0001990660910000154
将公式(15)带入公式(28),可得
Figure BDA0001990660910000161
为了确保系统的稳定性,选择适当的参数k22来确保k22-I3×3≥0。
在外界扰动和内在误差无法估计的时候,提出非线性ESO控制策略,对内部不可建模和外部不可估计的扰动进行估计和补偿:
ESO可用于姿态子系统,以估计和补偿外部环境的总干扰,对于这种情况,系统扩张状态方程是
Figure BDA0001990660910000162
状态观察器,即ESO是为了扩张状态而建立的
Figure BDA0001990660910000163
选择适当的参数βi(i=1,2,3),d4可以被估计和补偿;
类似地,φ和ψ也可以使用扩展状态观察器来估计误差,这样也可以通过ESO进行补偿;
通过将自适应跟踪控制方法与ESO相结合,可以得到新的姿态角控制律,如下所示
Figure BDA0001990660910000164
其中di是对系统总扰动的估计,δ=[θ φ ψ]T为无人机姿态角,U=[u2 u3 u4]T,
Figure BDA0001990660910000171
I=[I1 I2 I3]T代表惯性矩,dF=[d4 d5 d6]T代表阻力系数。
为了验证本方案的效果以及可行性,使用Matlab/Simulink平台构建系统仿真框图,通过仿真结果分析并验证本发明所提出算法的正确性及有效性,在仿真环境中,设置四旋翼无人机控制系统的初始位置[2 0 1 0 0 0],并且系统的初始角度为[0 0 0 0 0 0]。为了保证闭环系统的稳定性,内环收敛速度必须大于外环收敛速度。因此,取内环控制器的增益远远大于外环控制器的增益,选择
Figure BDA0001990660910000172
选择预期值ξd=[xd yd zd]T=[1 2 3]Td=π/3。
如图3所示,本发明用滑模控制方法来控制四旋翼无人机的位置子系统,通过图3可知,位置坐标x,y,z快速收敛到它们对应的期望值。
如图4所示,θd和φd是上面提到的虚拟姿态角,这里选择ψd=π/3。本发明提出自适应跟踪控制方法来控制四旋翼无人机的姿态子系统,由图4可知,姿态角δ=[θ φ ψ]T快速收敛到它们对应的期望值。在这种情况下,输入的饱和极限是umax=100,如果选择较小的值,则系统将不稳定。
为了反映本发明对输入饱和度控制的效果,设置不同的参数进行比较。在输入饱和受限依然存在的情况下,控制器不采用本发明所提出的饱和控制方法,选择umax=60,姿态角的跟踪如图5所示。比较图4和图5,图4的控制效果明显优于图5。从该结果可以分析出当输入饱和度受限时,本发明所提出饱和控制的方法可以有效地控制四旋翼无人机的姿态。
系统发输入控制规律如图6所示,控制输入u1控制位置子系统,其余三个输入控制姿态子系统。
本发明公开了具有输入约束的四旋翼直升机系统的滑模控制和自适应跟踪控制,考虑到输入约束和系统干扰,提出了具有Lyapunov稳定性分析的自适应跟踪控制和ESO方法;在输入约束研究的发展过程中,引入辅助设计系统来分析输入约束的影响,并将辅助设计系统的状态引入自适应跟踪控制方法中;并且通过Lyapunov稳定性分析,证明了所提出的滑模控制方法及具有输入约束的自适应跟踪控制在双闭环控制系统中,所有信号都具有渐近稳定性。最后,通过仿真,验证了所提出的滑模控制方法和具有输入约束的自适应控制方法的有效性。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (4)

1.一种基于饱和受限情况下的无人机自适应跟踪控制算法,其特征在于,包括以下步骤:
步骤A、在考虑扰动存在的情况下,建立四旋翼无人机的二阶系统数学模型,步骤A具体为:所述二阶系统数学模型为
Figure FDA0003466360540000011
其中,状态x,y,z表示四旋翼的位置,状态θ,φ,ψ由三个姿态角组成,g是重力加速度,Ii(i=1,2,3)代表惯性矩,Ki(i=1,2…6)代表阻力系数,di(i=1,2…6)是系统扰动,滚动角φ和俯仰角θ的范围(-2/π<φ、θ<2/π),偏航角ψ的范围也限于(-π<ψ<π),ui(i=1,2,3,4)表示四旋翼的四个控制输入量;
将公式(1)所表示的四旋翼无人机数学模型拆分成位置和姿态两部分,以便于分布式处理:
Figure FDA0003466360540000012
Figure FDA0003466360540000013
将四旋翼无人机模型写成矩阵形式,因此,公式(2)可以用以下形式书写:
Figure FDA0003466360540000014
其中,ξ=[x y z]T是四旋翼无人机在惯性坐标系中的质心位置,
Figure FDA0003466360540000021
是垂直单位矢量,g是重力加速度,令dF=[md1 md2 md3]T,m是四旋翼无人机的质量;
Figure FDA0003466360540000022
Figure FDA0003466360540000023
是旋转矩阵,其中c·=cos(·)、s·=sin(·);
步骤B、基于四旋翼无人机的二阶系统模型,对于外环位置子系统,提出滑模控制器对无人机系统进行控制,以实现对无人机位置状态的控制,步骤B具体为:所述滑模控制器为:
定义一个滑模控制函数
Figure FDA0003466360540000024
其中,λ1是一个正数并满足Hurwitz条件,并且s1(t)表示平动的滑模面;
选择指数收敛的控制律,滑模面的导数应如下
Figure FDA0003466360540000025
其中ε1≥DF/m和k1是正常数,符号函数定义如下
Figure FDA0003466360540000026
考虑公式(5),公式(6)的导数是
Figure FDA0003466360540000027
因此,对比公式(5)和公式(8),位置子系统的滑模控制律可以设计如下
Figure FDA0003466360540000028
步骤C、基于四旋翼无人机的二阶系统模型,对于内环姿态子系统,提出自适应控制策略对无人机系统进行控制,以实现对无人机姿态状态的控制,步骤C具体为:所述自适应控制器为:
假设预期跟踪的参考姿态是δd=[θd φd ψd]T并且虚拟控制律是αi=[α1 α2 α3]T;定义错误变量e1=δ-δd和e2=δ1i,虚拟控制定律αi可以定义如下
Figure FDA0003466360540000031
其中I3×3表示单位矩阵;
考虑公式(3)并求导,我们得到可用位置子系统的误差方程
Figure FDA0003466360540000032
其中δ=[θ φ ψ]T为无人机姿态角,U=[u2 u3 u4]T,
Figure FDA0003466360540000033
I=[I1 I2 I3]T代表惯性矩;
考虑Lyapunov函数
Figure FDA0003466360540000034
Vi求导可得
Figure FDA0003466360540000035
引入αi,可得
Figure FDA0003466360540000036
为方便输入饱和受限影响的分析,给出了辅助设计系统
Figure FDA0003466360540000037
其中
Figure FDA0003466360540000038
η∈R3是辅助设计系统的状态,设计参数是正常数,应根据跟踪性能的要求选择适当的值;
考虑输入饱和的影响,请选择以下控制律
Figure FDA0003466360540000041
2.根据权利要求1所述的基于饱和受限情况下的无人机自适应跟踪控制算法,其特征在于,其对于四旋翼无人机在获得虚拟输入之后,还得计算实际升程控制输入并将其传递到姿态子系统中间命令信号,具体包括以下步骤:
将虚拟控制输入写为向量Uξ=[ux,uy,uz]T,将中间命令信号设置为Θd=[φddd]T表示无人机的虚拟姿态角;
考虑
Figure FDA0003466360540000042
Uξ=[ux,uy,uz]T
Figure FDA0003466360540000043
可以得到
Figure FDA0003466360540000044
进一步
Figure FDA0003466360540000045
变为矩阵形式
Figure FDA0003466360540000046
使用Θd代替姿态子系统的中间命令信号Θ,获得
Figure FDA0003466360540000047
因为
Figure FDA0003466360540000048
可得
Figure FDA0003466360540000051
滚转角指令信号是
Figure FDA0003466360540000052
俯仰角指令信号是
Figure FDA0003466360540000053
偏航角ψd可用于跟踪给定指令信号的任意偏航角。
3.根据权利要求2所述基于饱和受限情况下的无人机自适应跟踪控制算法,其特征在于,其对于四旋翼无人机的稳定性采用如下方式验证:是利用Lyapunov稳定性理论,设计Lyapunov函数证明在四旋翼无人机位置子系统和姿态子系统的稳定性;具体包括以下步骤:
步骤a:对于采用滑模控制策略的位置子系统,考虑Lyapunov函数
Figure FDA0003466360540000054
V1的导数如下:
Figure FDA0003466360540000055
由于ε1≥DF/m,可得
Figure FDA0003466360540000056
从上面的分析我们可以看出s1是指数收敛,即eξ指数收敛;类似地,位置误差子系统是指数稳定的;
步骤b:对于采用自适应控制策略的姿态子系统,考虑Lyapunov函数
Figure FDA0003466360540000061
根据公式(27),V的导数
Figure FDA0003466360540000062
将公式(15)带入公式(28),可得
Figure FDA0003466360540000063
为了确保系统的稳定性,选择适当的参数k22来确保k22-I3×3≥0。
4.根据权利要求3所述基于饱和受限情况下的无人机自适应跟踪控制算法,其特征在于,在外界扰动和内在误差无法估计的时候,提出非线性ESO控制策略,对内部不可建模和外部不可估计的扰动进行估计和补偿:
ESO可用于姿态子系统,以估计和补偿外部环境的总干扰,对于这种情况,系统扩张状态方程是
Figure FDA0003466360540000071
状态观察器,即ESO是为了扩张状态而建立的
Figure FDA0003466360540000072
选择适当的参数βi(i=1,2,3),d4可以被估计和补偿;
类似地,φ和ψ也可以使用扩展状态观察器来估计误差,这样也可以通过ESO进行补偿;
通过将自适应跟踪控制方法与ESO相结合,可以得到新的姿态角控制律,如下所示
Figure FDA0003466360540000073
其中di是对系统总扰动的估计,δ=[θ φ ψ]T为无人机姿态角,U=[u2 u3 u4]T,
Figure FDA0003466360540000074
I=[I1 I2 I3]T代表惯性矩,dF=[d4 d5 d6]T代表阻力系数。
CN201910179402.5A 2019-03-11 2019-03-11 基于饱和受限情况下的无人机自适应跟踪控制算法 Active CN109884895B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910179402.5A CN109884895B (zh) 2019-03-11 2019-03-11 基于饱和受限情况下的无人机自适应跟踪控制算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910179402.5A CN109884895B (zh) 2019-03-11 2019-03-11 基于饱和受限情况下的无人机自适应跟踪控制算法

Publications (2)

Publication Number Publication Date
CN109884895A CN109884895A (zh) 2019-06-14
CN109884895B true CN109884895B (zh) 2022-03-11

Family

ID=66931555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910179402.5A Active CN109884895B (zh) 2019-03-11 2019-03-11 基于饱和受限情况下的无人机自适应跟踪控制算法

Country Status (1)

Country Link
CN (1) CN109884895B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111736550B (zh) * 2019-06-15 2024-04-05 沈阳工业大学 一种植保机单钟摆喷杆悬架系统的非线性控制方法
CN110376898A (zh) * 2019-08-14 2019-10-25 西南石油大学 四旋翼无人机的快速终端滑模自适应控制系统及其方法
CN110502027B (zh) * 2019-09-16 2022-08-12 南京邮电大学 一种基于自适应终端滑模的四旋翼无人机姿态容错控制方法
CN110673479A (zh) * 2019-09-30 2020-01-10 中国科学院长春光学精密机械与物理研究所 轨迹跟踪优化控制方法、装置、多旋翼无人机及存储装置
CN110673619B (zh) * 2019-10-21 2022-06-17 深圳市道通智能航空技术股份有限公司 一种飞行姿态的控制方法、装置、无人机和存储介质
CN111413994B (zh) * 2020-03-13 2023-08-08 浙江树人学院(浙江树人大学) 一种四旋翼无人机直接自适应模糊控制方法
CN111443724B (zh) * 2020-04-20 2022-03-11 中南大学 一种基于自抗扰混合切换控制的四旋翼无人机的控制方法
CN111694278B (zh) * 2020-06-28 2021-07-23 南京工业大学 一种四旋翼无人机的鲁棒跟踪控制方法及系统
CN111880553B (zh) * 2020-08-11 2022-10-11 北京航空航天大学 一种考虑惯量不确定性的四旋翼无人机姿态控制方法
CN112462796B (zh) * 2020-11-28 2023-01-17 中国人民解放军海军航空大学青岛校区 一种刚体飞行器姿态角稳定的自适应反演控制系统及方法
CN112835372B (zh) * 2020-11-30 2024-05-07 天津大学 四旋翼无人机的固定时间控制方法
CN112947513B (zh) * 2021-01-27 2022-10-21 西北工业大学 一种基于容错和抗饱和机制的四旋翼无人机姿态控制方法
CN113359824A (zh) * 2021-05-31 2021-09-07 杭州电子科技大学 基于模糊模型的无人机集群控制方法
CN113867374B (zh) * 2021-07-23 2023-12-01 福州大学 基于滑模控制的四旋翼无人机参数预测和扰动的自适应轨迹跟踪控制器及其设计方法
CN114237270B (zh) * 2021-12-08 2023-11-07 南京航空航天大学 一种考虑输入饱和的无人直升机跟踪控制方法
CN114967718B (zh) * 2022-05-06 2023-03-28 四川大学 一种控制受限下的四旋翼轨迹跟踪最优控制方法
CN117891177A (zh) * 2024-03-15 2024-04-16 国网浙江省电力有限公司宁波供电公司 无人机控制器模型构建方法、装置、设备及存储介质
CN117908576B (zh) * 2024-03-19 2024-05-28 南京信息工程大学 一种抗饱和的四旋翼无人机有限时间容错控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106774373A (zh) * 2017-01-12 2017-05-31 哈尔滨工业大学 一种四旋翼无人机有限时间姿态跟踪控制方法
CN106842953A (zh) * 2017-03-13 2017-06-13 贾杰 一种无人直升机自适应低阶控制器
CN107491083A (zh) * 2017-09-11 2017-12-19 北京航空航天大学 一种基于饱和自适应滑模控制的四旋翼自主着船方法
CN109283932A (zh) * 2018-09-18 2019-01-29 浙江工业大学 一种基于积分反步滑模的四旋翼飞行器姿态控制方法
CN109343369A (zh) * 2018-11-19 2019-02-15 南京邮电大学 一种基于非线性观测器的四旋翼容错控制器设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106774373A (zh) * 2017-01-12 2017-05-31 哈尔滨工业大学 一种四旋翼无人机有限时间姿态跟踪控制方法
CN106842953A (zh) * 2017-03-13 2017-06-13 贾杰 一种无人直升机自适应低阶控制器
CN107491083A (zh) * 2017-09-11 2017-12-19 北京航空航天大学 一种基于饱和自适应滑模控制的四旋翼自主着船方法
CN109283932A (zh) * 2018-09-18 2019-01-29 浙江工业大学 一种基于积分反步滑模的四旋翼飞行器姿态控制方法
CN109343369A (zh) * 2018-11-19 2019-02-15 南京邮电大学 一种基于非线性观测器的四旋翼容错控制器设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Trajectory Tracking Control for Quadrotor Robot Subject to Payload Variation and Wind Gust Disturbance;chen wang 等;《Journal of Intelligent & Robotic Systems》;20160118;第319-325页 *
四旋翼无人机鲁棒抗扰控制技术研究;杨翼;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180315(第03期);第34-36、46页 *

Also Published As

Publication number Publication date
CN109884895A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN109884895B (zh) 基于饱和受限情况下的无人机自适应跟踪控制算法
Labbadi et al. Robust adaptive backstepping fast terminal sliding mode controller for uncertain quadrotor UAV
Ramirez-Rodriguez et al. Robust backstepping control based on integral sliding modes for tracking of quadrotors
CN109765918B (zh) 一种无人直升机鲁棒自适应补偿控制方法
Mokhtari et al. A new robust control for minirotorcraft unmanned aerial vehicles
CN107765553B (zh) 针对旋翼无人机吊挂运输系统的非线性控制方法
Liu et al. Robust fault-tolerant formation control for tail-sitters in aggressive flight mode transitions
CN109597426A (zh) 基于l1自适应控制的四旋翼直升机轨迹跟踪控制方法
CN112684805B (zh) 一种考虑姿态约束的高机动微型无人机控制方法
Jiang et al. Composite adaptive finite-time control for quadrotors via prescribed performance
CN112327926B (zh) 一种无人机编队的自适应滑模控制方法
Lv et al. Extended state observer based MPC for a quadrotor helicopter subject to wind disturbances
Sun et al. Nonlinear robust compensation method for trajectory tracking control of quadrotors
Xiong et al. Self-tuning sliding mode control for an uncertain coaxial octorotor UAV
CN108459611B (zh) 一种近空间飞行器的姿态跟踪控制方法
Aruneshwaran et al. Neural adaptive flight controller for ducted-fan UAV performing nonlinear maneuver
CN111435253A (zh) 四旋翼飞行器有界输出控制
CN109634110B (zh) 一种基于非线性干扰观测器的尾座式无人机过渡段鲁棒控制方法
Schwung et al. Quadrotor tracking control: Design and experiments
Chen et al. Decoupling attitude control of a hypersonic glide vehicle based on a nonlinear extended state observer
El Ferik et al. Robust FOSMC of quadrotor in the presence of slung load
CN115079715B (zh) 角速度及输入约束下分布式航天器编队姿态同步控制方法
Gao et al. Event-triggered tracking control scheme for quadrotors with external disturbances: theory and validations
Housny et al. Robust sliding mode control for quadrotor UAV
CN112835372B (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