CN108363298A - 一种基于四元数描述的四旋翼无人机快速收敛控制方法 - Google Patents
一种基于四元数描述的四旋翼无人机快速收敛控制方法 Download PDFInfo
- Publication number
- CN108363298A CN108363298A CN201810043419.3A CN201810043419A CN108363298A CN 108363298 A CN108363298 A CN 108363298A CN 201810043419 A CN201810043419 A CN 201810043419A CN 108363298 A CN108363298 A CN 108363298A
- Authority
- CN
- China
- Prior art keywords
- quadrotor drone
- formula
- indicate
- quadrotor
- drone
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明属于飞行控制技术领域,特别涉及一种基于四元数描述的四旋翼无人机快速收敛控制方法。本发明首先建立四旋翼无人机系统模型,四旋翼无人机系统分为位置子系统和姿态子系统,设计四旋翼无人机有限时间位置控制器,设计四旋翼无人机有限时间姿态控制器,最终得到四旋翼无人机的电机转速控制器,根据四旋翼无人机的电机转速控制器控制飞行器飞行。通过采用有限时间控制技术,提高了闭环系统的收敛速度,使得系统能够在有限时间内更快的收敛到期望值;本四旋翼无人机快速收敛控制方法相比于传统的PID控制方法,具备更好的抗干扰性能和鲁棒性,有效地减小了四旋翼无人机在实际飞行过程中各种外界干扰的影响。
Description
技术领域
本发明属于飞行控制技术领域,特别涉及一种基于四元数描述的四旋翼无人机快速收敛控制方法。
背景技术
四旋翼无人机是一种拥有六个自由度能够垂直升降,室内飞行,定点悬停的无人无人机,由于其结构简单、机动性能好、维护成本低等优势,已经被广泛地应用于航空拍摄、环境监测、地理探索、搜索及救援等领域,具有重要的研究意义和应用前景。在控制领域界,四旋翼无人机是一个典型的非线性模型,具有多变量、强耦合和欠驱动的特点,而且易受各种外界环境的干扰,如:风力干扰,因此其控制问题受到了国内外学者的广泛关注。
目前,针对四旋翼无人无人机的控制问题,国内外的许多大学和研究机构利用各种非线性控制理论设计出了不同的控制器。其控制算法一般只保证了闭环系统的渐近稳定性,但是无人机收敛到期望状态的调节时间是无穷大,从收敛时间角度分析,如果无人机能够在有限时间内收敛到期望值显然具有重要意义,因此亟需提出一种拥有更快的收敛速度的四旋翼无人机快速收敛控制方法。
发明内容
本发明为了克服上述现有技术的不足,提供了一种基于四元数描述的四旋翼无人机快速收敛控制方法,本发明能够在有限时间内收敛到期望值,本方法使四旋翼无人机具备更好的抗干扰性能和鲁棒性,可以有效地抑制外界干扰。
为实现上述目的,本发明采用了以下技术措施:
一种基于四元数描述的四旋翼无人机快速收敛控制方法,包括以下步骤:
S1、建立四旋翼无人机系统模型,四旋翼无人机系统分为位置子系统和姿态子系统;
S2、设计四旋翼无人机有限时间位置控制器;
S3、设计四旋翼无人机有限时间姿态控制器;
S4、设计四旋翼无人机的电机转速控制器;
S5、根据四旋翼无人机的电机转速控制器控制飞行器飞行。
优选的,步骤S1的具体操作步骤包括:
四旋翼无人机系统的位置子系统模型描述如公式(1)所示:
χ=[x,y,z]T∈R3表示四旋翼无人机的实际位置,Θ=[φ,θ,ψ]T∈R3表示四旋翼无人机的实际欧拉角姿态,R3表示三维的向量空间,其中,
x表示四旋翼无人机的x坐标,表示x的一阶导数,表示x的二阶导数;
y表示四旋翼无人机的y坐标,表示y的一阶导数,表示y的二阶导数;
z表示四旋翼无人机的z坐标,表示z的一阶导数,表示z的二阶导数;
φ表示四旋翼无人机的横滚角,
θ表示四旋翼无人机的俯仰角,
ψ表示四旋翼无人机的偏航角,
m表示四旋翼无人机的质量,g表示重力加速度,F表示由四旋翼无人机产生的总的升力,K1,K2,K3表示空气动力阻尼系数;
四旋翼无人机系统的姿态子系统模型描述如公式(2)所示:
其中,表示基于四元数描述的四旋翼无人机姿态,即实际的四元数姿态;[]T为转置矩阵,qv=[q1,q2,q3]T为向量,q0、q1、q2、q3分别为四元数的四个参数,J∈R3×3表示四旋翼无人机对称的惯性矩阵,Ω=[Ω1,Ω2,Ω3]∈R3表示x、y、z三个方向上的角速度,R3×3为3×3的向量空间,为Ω的一阶导数,Ω×为斜对称矩阵,τ∈R3表示由电机转动产生的力矩,矩阵E(q)如公式(3)所示:
其中,I3表示3×3的单位矩阵,当Ω=(Ω1,Ω2,Ω3)T时,()T为转置矩阵,斜对称矩阵Ω×为公式(4):
优选的,步骤S2的具体操作步骤包括:
对于四旋翼无人机的位置子系统设计有限时间位置控制器,使四旋翼无人机的实际位置在有限时间跟踪上期望位置,并通过李雅普诺夫函数对位置误差系统进行稳定性分析,其设计的具体步骤包括:
对于公式(1)中的四旋翼无人机的位置子系统模型,定义三通道的虚拟控制信号如公式(5)所示,
其中,μx为x通道的虚拟控制信号,μy为y通道的虚拟控制信号,μz为z通道的虚拟控制信号;
设定四旋翼无人机的期望位置:
χd=[xd,yd,zd]T (6)
其中,xd表示四旋翼无人机的x轴期望位置,yd表示四旋翼无人机的y轴期望位置,zd表示四旋翼无人机的z轴期望位置;
位置误差χe为四旋翼无人机的期望位置减去四旋翼无人机的实际位置,如公式(7)所示:
χe=χd-χ=[xd-x,yd-y,zd-z]T=[ex,ey,ez]T (7)
其中,ex表示x轴方向的位置误差,ey表示y轴方向的位置误差,ez表示z轴方向的位置误差;
设计四旋翼无人机的有限时间位置控制器,定义新的饱和函数sat(s,α),如公式(8)所示:
其中,sign()为符号函数,||为绝对值符号,s为函数输入量,α为实数,对于公式(1)设计四旋翼无人机的有限时间位置控制器,如公式(9)所示:
其中,为xd的一阶导数,为xd的二阶导数,为yd的一阶导数,为yd的二阶导数,为zd的一阶导数,为zd的二阶导数,α1、α2、Kpx、Kpy、Kpz、Kdx、Kdy、Kdz均为控制器参数,且0<α1<1,α2=2α1/(1+α1),Kpx>0,Kpy>0,Kpz>0,Kdx>0,Kdy>0,Kdz>0;
四旋翼无人机沿x轴方向的位置误差ex在有限时间内收敛到零;
取李雅普诺夫函数V1,如公式(10)所示:
其中,V1为所选取的李雅普诺夫函数,ρ为积分变量;
对李雅普诺夫函数进行求导,如公式(11)所示:
四旋翼无人机沿x轴方向的位置误差ex及其导数在有限时间内进入区域|ex|<1,根据公式(8)可知,四旋翼无人机沿x轴方向的位置误差系统变为公式(12);
其中,为ex的二阶导数,公式(12)的高阶项部分为对于任意的ε>0和所有ex≠0,满足公式(13):
其中,ε为任意大于零的实数,r2为公式(12)的加权扩张量,n为公式(12)的齐次度;
结合公式(11)、公式(12)、公式(13),四旋翼无人机沿x轴方向的位置误差ex在有限时间内收敛到零,四旋翼无人机沿y轴方向的位置误差ey以及四旋翼无人机沿z轴方向的位置误差ez均在有限时间内收敛到零,根据公式(8)可知,位置误差χe在有限时间内收敛到零。
优选的,步骤S3的具体操作步骤包括:
对于四旋翼无人机的姿态子系统模型设计有限时间控制器,使四旋翼无人机的实际姿态在有限时间内跟踪上期望姿态,并通过李雅普诺夫函数对姿态误差系统进行稳定性分析,其设计的具体步骤包括:
对于公式(5)中的虚拟控制信号,反解出四旋翼无人机的期望欧拉角姿态为公式(14);
其中,φd为期望的横滚角,θd为期望的俯仰角,ψd为期望的偏航角;
设定期望的偏航角姿态为零,即ψd≡0;
根据欧拉角至四元数转换方程,获得期望的四元数姿态如公式(15)所示;
其中,qd0、qd1、qd2、qd3均为期望的四元数的参数,R4为四维的向量空间;
由公式(2)的四旋翼无人机系统的姿态子系统模型可知,期望的角速度Ωd和期望的四元数转换关系为公式(16);
其中,为期望的四元数的一阶导数;E(qd)为矩阵;
根据公式(16)得出期望的角速度Ωd;
定义期望的四元数姿态qd和实际的四元数姿态q之间的差值为四元数姿态误差e,
获得基于四元数姿态误差的四旋翼无人机的姿态误差方程为公式(18);
其中,角速度误差Ωe=Ω-CΩd,参数矩阵E(e)满足公式(3);
对于公式(18)设计有限时间姿态控制器,如公式(19);
其中,β1,β2为正增益,r2=1+r,r3=1+2r,r为偶数和奇数之比;
定义
构造李雅普诺夫函数为公式(21);
其中,e0,ev满足公式(17),Ωe,j为参数Ωe的第j个分量,j={1,2,3};
其中虚拟角速度误差为公式(22),
其中,为虚拟角速度误差,介于实际角速度误差Ωe和虚拟角速度误差之间定义一个误差量ξ,
则,
其中c1,c2均为大于零的常数,ξj为误差量ξ的第j个分量,ej参数e的第j个分量,uj为参数u的第j个分量;
将公式(20)中的控制信号u,写成如公式(25)形式:
其中β2≥1+c1+c2,
将公式(25)带入公式(24),得到公式(26);
根据公式(21)的李雅普诺夫函数的定义可知:
其中,由此可知:
四元数姿态误差和角速度误差在有限时间内收敛到零。
进一步的,步骤S4的具体操作步骤包括:
获得四旋翼无人机的四个电机的转速如公式(29)所示:
其中,d为四旋翼无人机的旋翼中心到飞行器重心的距离,κ为旋翼的扭矩系数,b为旋翼的扭矩升力系数,κ和b均为正实数。
进一步的,将[φ,θ,ψ]T∈R3用于位置子系统中,将q=[q0,q1,q2,q3]T∈R4用于姿态子系统中,[φ,θ,ψ]T∈R3和q=[q0,q1,q2,q3]T∈R4的转化关系分别为公式(30)和公式(31);
本发明的有益效果在于:
1)、本发明通过采用有限时间控制技术,提高了闭环系统的收敛速度,使得系统能够在有限时间内更快的收敛到期望值;本四旋翼无人机快速收敛控制方法相比于传统的PID控制方法,具备更好的抗干扰性能和鲁棒性,有效地减小了四旋翼无人机在实际飞行过程中各种外界干扰的影响。
2)、本发明四旋翼无人机的姿态控制器设计通过采用四元数描述法,可以有效地避免姿态奇异现象发生,能够为四旋翼无人机提供一种全局姿态跟踪控制算法。
附图说明
图1为本发明采用的双环控制策略示意图;
图2为本发明在无扰动情况下,四旋翼无人机系统在有限时间控制器(FC)和PD控制器作用下的位置响应曲线;
图3为本发明在无扰动情况下,四旋翼无人机系统在有限时间控制器(FC)和PD控制器作用下的姿态响应曲线;
图4为本发明在存在扰动情况下,四旋翼无人机系统在有限时间控制器(FC)和PD控制器作用下的位置响应曲线;
图5为本发明在存在扰动情况下,四旋翼无人机系统在有限时间控制器(FC)和PD控制器作用下的姿态响应曲线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的仿真和实验都通过仿真软件MATLAB R2010中的Simulink工具箱搭建四旋翼无人机的模型来实现,图1为本发明采用的双环控制策略示意图。
基于四元数描述的四旋翼无人机快速收敛控制方法的具体步骤如下:
S1、建立四旋翼无人机系统模型,四旋翼无人机系统分为位置子系统和姿态子系统;
如图1所示,四旋翼无人机闭环控制系统采用内外双环控制结构,输入量为期望的空间位置,输出量为四旋翼无人机的实际位置及实际姿态,内环为姿态跟踪控制环,实现期望角度的跟踪,外环为位置控制环,实现四旋翼无人机的悬停控制。
四旋翼无人机系统模型描述如下:
其中,参数选取m=0.468kg,g=9.8m/s2,空气阻尼系数K1=0.01Ns/m,K2=0.01Ns/m,K3=0.01Ns/m,
位置初始值为(x,y,z)T=[0,0,0]T,姿态初始值为(q0,q1,q2,q3)T=[0.8832,0.3,-0.2,-0.3]T;
S2、设计四旋翼无人机有限时间位置控制器;
对于四旋翼无人机的位置子系统设计有限时间控制器,使得四旋翼无人机的实际位置在有限时间跟踪上期望位置,其具体步骤包括:
首先,对于公式(1)中位置模型,定义三通道的虚拟控制信号如公式(2)所示;
设定四旋翼无人机的期望位置:
位置误差χe为四旋翼无人机的期望位置减去四旋翼无人机的实际位置,如公式(4)所示:
χe=χd-χ=[xd-x,yd-y,zd-z]T (4)
因此,结合公式(2),对于公式(1)中四旋翼无人机的位置子系统,设计有限时间位置控制器如下:
其中,控制器参数选取为Kpx=4.2,Kpy=Kpz=4.5,Kdx=Kdy=Kdz=3.8,α1=3/4,α2=6/7,则四旋翼无人机可以在有限时间内收敛到期望的位置。
S3、设计四旋翼无人机有限时间姿态控制器;
对于四旋翼无人机的姿态子系统模型设计有限时间控制器,使得四旋翼无人机的实际在有限时间跟踪上期望姿态,其具体步骤包括:
结合公式(2),获得四旋翼无人机的期望欧拉角姿态为公式(6);
根据欧拉角至四元数转换方程,获得期望的四元数姿态如公式(7)所示;
期望的角速度Ωd可以根据公式(8)得出;
获得时间姿态和期望姿态支架的四元数姿态误差e:
对于公式(1)中四旋翼无人机的姿态子系统,设计有限时间姿态控制器,如公式(10);
其中,角速度误差Ωe=Ω-CΩd,参数β1=6,β2=5,r=-2/7,即r2=5/7,r3=3/7,则四旋翼无人机可以在有限时间内跟踪上期望的姿态;
本发明分别考虑不存在外部干扰和存在外部干扰两种情况下系统的性能,给出四旋翼无人机在有限时间FC控制器和经典PD控制器作用下的位置、姿态对比曲线。
情况一:不存在外部干扰情况下的系统响应;
不存在外部干扰时,四旋翼无人机在有限时间控制器FC和经典PD控制器作用下的位置响应曲线见图2,姿态响应曲线见图3,比较结果见表1:
表1:
情况二:存在外部干扰情况下的系统响应;
存在外部干扰时,四旋翼无人机在有限时间控制器FC和经典PD控制器作用下的位置响应曲线见图4,姿态响应曲线见图5。
本发明提供的一种基于四元数描述的四旋翼无人机快速收敛控制方法,实现四旋翼无人机的定点悬停控制,上述具体实施例,为进行对比分别采用有限时间控制以及经典PD控制。图2和图3是系统无扰动时,在有限时间控制器FC和PD控制器作用下的位置、姿态响应曲线。对比两者,有限时间控制器可以提供一个更小超调量和更快收敛性能的控制效果,图4和图5是系统加入了扰动时,在有限时间控制器FC和PD控制器作用下的位置、姿态响应曲线,结果表明,相对于PD控制,有限时间控制算法FC具有较快的收敛性和较强的抗扰动能力。
综上所述,本发明可以拥有更快的收敛速度,更好的抗干扰能力和鲁棒性。
Claims (6)
1.一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于,包括以下步骤:
S1、建立四旋翼无人机系统模型,四旋翼无人机系统分为位置子系统和姿态子系统;
S2、设计四旋翼无人机有限时间位置控制器;
S3、设计四旋翼无人机有限时间姿态控制器;
S4、设计四旋翼无人机的电机转速控制器;
S5、根据四旋翼无人机的电机转速控制器控制飞行器飞行。
2.如权利要求1所述的一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于,步骤S1的具体操作步骤包括:
四旋翼无人机系统的位置子系统模型描述如公式(1)所示:
χ=[x,y,z]T∈R3表示四旋翼无人机的实际位置,Θ=[φ,θ,ψ]T∈R3表示四旋翼无人机的实际欧拉角姿态,R3表示三维的向量空间,其中,
x表示四旋翼无人机的x坐标,表示x的一阶导数,表示x的二阶导数;
y表示四旋翼无人机的y坐标,表示y的一阶导数,表示y的二阶导数;
z表示四旋翼无人机的z坐标,表示z的一阶导数,表示z的二阶导数;
φ表示四旋翼无人机的横滚角,
θ表示四旋翼无人机的俯仰角,
ψ表示四旋翼无人机的偏航角,
m表示四旋翼无人机的质量,g表示重力加速度,F表示由四旋翼无人机产生的总的升力,K1,K2,K3表示空气动力阻尼系数;
四旋翼无人机系统的姿态子系统模型描述如公式(2)所示:
其中,表示基于四元数描述的四旋翼无人机姿态,即实际的四元数姿态;[]T为转置矩阵,qv=[q1,q2,q3]T为向量,q0、q1、q2、q3分别为四元数的四个参数,J∈R3×3表示四旋翼无人机对称的惯性矩阵,Ω=[Ω1,Ω2,Ω3]∈R3表示x、y、z三个方向上的角速度,R3×3为3×3的向量空间,为Ω的一阶导数,Ω×为斜对称矩阵,τ∈R3表示由电机转动产生的力矩,矩阵E(q)如公式(3)所示:
其中,I3表示3×3的单位矩阵,当Ω=(Ω1,Ω2,Ω3)T时,()T为转置矩阵,斜对称矩阵Ω×为公式(4):
3.如权利要求2所述的一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于,步骤S2的具体操作步骤包括:
对于四旋翼无人机的位置子系统设计有限时间位置控制器,使四旋翼无人机的实际位置在有限时间跟踪上期望位置,并通过李雅普诺夫函数对位置误差系统进行稳定性分析,其设计的具体步骤包括:
对于公式(1)中的四旋翼无人机的位置子系统模型,定义三通道的虚拟控制信号如公式(5)所示,
其中,μx为x通道的虚拟控制信号,μy为y通道的虚拟控制信号,μz为z通道的虚拟控制信号;
设定四旋翼无人机的期望位置:
χd=[xd,yd,zd]T (6)
其中,xd表示四旋翼无人机的x轴期望位置,yd表示四旋翼无人机的y轴期望位置,zd表示四旋翼无人机的z轴期望位置;
位置误差χe为四旋翼无人机的期望位置减去四旋翼无人机的实际位置,如公式(7)所示:
χe=χd-χ=[xd-x,yd-y,zd-z]T=[ex,ey,ez]T (7)
其中,ex表示x轴方向的位置误差,ey表示y轴方向的位置误差,ez表示z轴方向的位置误差;
设计四旋翼无人机的有限时间位置控制器,定义新的饱和函数sat(s,α),如公式(8)所示:
其中,sign()为符号函数,| |为绝对值符号,s为函数输入量,α为实数,对于公式(1)设计四旋翼无人机的有限时间位置控制器,如公式(9)所示:
其中,为xd的一阶导数,为xd的二阶导数,为yd的一阶导数,为yd的二阶导数,为zd的一阶导数,为zd的二阶导数,α1、α2、Kpx、Kpy、Kpz、Kdx、Kdy、Kdz均为控制器参数,且0<α1<1,α2=2α1/(1+α1),Kpx>0,Kpy>0,Kpz>0,Kdx>0,Kdy>0,Kdz>0;
四旋翼无人机沿x轴方向的位置误差ex在有限时间内收敛到零;
取李雅普诺夫函数V1,如公式(10)所示:
其中,V1为所选取的李雅普诺夫函数,ρ为积分变量;
对李雅普诺夫函数进行求导,如公式(11)所示:
四旋翼无人机沿x轴方向的位置误差ex及其导数在有限时间内进入区域|ex|<1,根据公式(8)可知,四旋翼无人机沿x轴方向的位置误差系统变为公式(12);
其中,为ex的二阶导数,公式(12)的高阶项部分为对于任意的ε>0和所有ex≠0,满足公式(13):
其中,ε为任意大于零的实数,r2为公式(12)的加权扩张量,n为公式(12)的齐次度;
结合公式(11)、公式(12)、公式(13),四旋翼无人机沿x轴方向的位置误差ex在有限时间内收敛到零,四旋翼无人机沿y轴方向的位置误差ey以及四旋翼无人机沿z轴方向的位置误差ez均在有限时间内收敛到零,根据公式(8)可知,位置误差χe在有限时间内收敛到零。
4.如权利要求3所述的一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于,步骤S3的具体操作步骤包括:
对于四旋翼无人机的姿态子系统模型设计有限时间控制器,使四旋翼无人机的实际姿态在有限时间内跟踪上期望姿态,并通过李雅普诺夫函数对姿态误差系统进行稳定性分析,其设计的具体步骤包括:
对于公式(5)中的虚拟控制信号,反解出四旋翼无人机的期望欧拉角姿态为公式(14);
其中,φd为期望的横滚角,θd为期望的俯仰角,ψd为期望的偏航角;
设定期望的偏航角姿态为零,即ψd≡0;
根据欧拉角至四元数转换方程,获得期望的四元数姿态如公式(15)所示;
其中,qd0、qd1、qd2、qd3均为期望的四元数的参数,R4为四维的向量空间;
由公式(2)的四旋翼无人机系统的姿态子系统模型可知,期望的角速度Ωd和期望的四元数转换关系为公式(16);
其中,为期望的四元数的一阶导数;E(qd)为矩阵;
根据公式(16)得出期望的角速度Ωd;
定义期望的四元数姿态qd和实际的四元数姿态q之间的差值为四元数姿态误差e,
获得基于四元数姿态误差的四旋翼无人机的姿态误差方程为公式(18);
其中,角速度误差Ωe=Ω-CΩd,参数矩阵E(e)满足公式(3);
对于公式(18)设计有限时间姿态控制器,如公式(19);
其中,β1,β2为正增益,r2=1+r,r3=1+2r,r为偶数和奇数之比;
定义
构造李雅普诺夫函数为公式(21);
其中,e0,ev满足公式(17),Ωe,j为参数Ωe的第j个分量,j={1,2,3};
其中虚拟角速度误差为公式(22),
其中,为虚拟角速度误差,介于实际角速度误差Ωe和虚拟角速度误差之间定义一个误差量ξ,
则,
其中c1,c2均为大于零的常数,ξj为误差量ξ的第j个分量,ej参数e的第j个分量,uj为参数u的第j个分量;
将公式(20)中的控制信号u,写成如公式(25)形式:
其中β2≥1+c1+c2,
将公式(25)带入公式(24),得到公式(26);
根据公式(21)的李雅普诺夫函数的定义可知:
其中,由此可知:
四元数姿态误差和角速度误差在有限时间内收敛到零。
5.如权利要求4所述的一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于,步骤S4的具体操作步骤包括:
获得四旋翼无人机的四个电机的转速如公式(29)所示:
其中,d为四旋翼无人机的旋翼中心到飞行器重心的距离,κ为旋翼的扭矩系数,b为旋翼的扭矩升力系数,κ和b均为正实数。
6.如权利要求5所述的一种基于四元数描述的四旋翼无人机快速收敛控制方法,其特征在于:将[φ,θ,ψ]T∈R3用于位置子系统中,将q=[q0,q1,q2,q3]T∈R4用于姿态子系统中,[φ,θ,ψ]T∈R3和q=[q0,q1,q2,q3]T∈R4的转化关系分别为公式(30)和公式(31);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810043419.3A CN108363298A (zh) | 2018-01-17 | 2018-01-17 | 一种基于四元数描述的四旋翼无人机快速收敛控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810043419.3A CN108363298A (zh) | 2018-01-17 | 2018-01-17 | 一种基于四元数描述的四旋翼无人机快速收敛控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108363298A true CN108363298A (zh) | 2018-08-03 |
Family
ID=63006296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810043419.3A Pending CN108363298A (zh) | 2018-01-17 | 2018-01-17 | 一种基于四元数描述的四旋翼无人机快速收敛控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108363298A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885074A (zh) * | 2019-02-28 | 2019-06-14 | 天津大学 | 四旋翼无人机有限时间收敛姿态控制方法 |
CN110032205A (zh) * | 2019-04-29 | 2019-07-19 | 河海大学常州校区 | 一种具有抗干扰能力的无人机姿态控制方法 |
CN110275542A (zh) * | 2019-06-14 | 2019-09-24 | 合肥工业大学 | 一种基于自适应有限时间控制的四旋翼飞行器控制方法 |
CN110673623A (zh) * | 2019-11-04 | 2020-01-10 | 苏州经贸职业技术学院 | 一种基于双环pd控制算法控制的四旋翼无人机着陆方法 |
CN111580537A (zh) * | 2020-05-28 | 2020-08-25 | 西北工业大学 | 一种无人机特技飞行控制系统及方法 |
WO2020177267A1 (zh) * | 2019-03-07 | 2020-09-10 | 中国科学院深圳先进技术研究院 | 四旋翼无人机的控制方法、装置、设备及可读介质 |
CN111984027A (zh) * | 2019-05-23 | 2020-11-24 | 南京航空航天大学 | 异构多无人机一致性控制方法 |
CN113759706A (zh) * | 2021-11-05 | 2021-12-07 | 普宙科技(深圳)有限公司 | 基于轴角矢量的多旋翼姿态自抗扰控制方法及系统 |
CN114995483A (zh) * | 2022-06-28 | 2022-09-02 | 四川大学 | 基于对偶四元数的倾转旋翼多旋翼无人机跟踪控制方法 |
CN115309172A (zh) * | 2022-07-01 | 2022-11-08 | 中国人民解放军国防科技大学 | 一种用于运输血液的无人机的姿态控制方法 |
CN117055602A (zh) * | 2023-09-05 | 2023-11-14 | 广东工业大学 | 一种四旋翼无人机的全姿态控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2324624C1 (ru) * | 2006-11-17 | 2008-05-20 | ОАО "ОКБ им. А.С. Яковлева" | Самолет с системой формирования резервной информации о пространственном положении летательного аппарата |
JP5061264B1 (ja) * | 2012-03-23 | 2012-10-31 | 国立大学法人 千葉大学 | 小型姿勢センサ |
CN106527462A (zh) * | 2016-12-02 | 2017-03-22 | 广西师范大学 | 无人机控制装置 |
CN106774373A (zh) * | 2017-01-12 | 2017-05-31 | 哈尔滨工业大学 | 一种四旋翼无人机有限时间姿态跟踪控制方法 |
CN107479371A (zh) * | 2017-07-03 | 2017-12-15 | 浙江工业大学 | 一种基于快速非奇异终端滑模的四旋翼无人机有限时间自适应控制方法 |
CN107490966A (zh) * | 2017-08-22 | 2017-12-19 | 浙江工业大学 | 一种基于改进幂次趋近律的飞行器有限时间自适应姿态控制方法 |
-
2018
- 2018-01-17 CN CN201810043419.3A patent/CN108363298A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2324624C1 (ru) * | 2006-11-17 | 2008-05-20 | ОАО "ОКБ им. А.С. Яковлева" | Самолет с системой формирования резервной информации о пространственном положении летательного аппарата |
JP5061264B1 (ja) * | 2012-03-23 | 2012-10-31 | 国立大学法人 千葉大学 | 小型姿勢センサ |
CN106527462A (zh) * | 2016-12-02 | 2017-03-22 | 广西师范大学 | 无人机控制装置 |
CN106774373A (zh) * | 2017-01-12 | 2017-05-31 | 哈尔滨工业大学 | 一种四旋翼无人机有限时间姿态跟踪控制方法 |
CN107479371A (zh) * | 2017-07-03 | 2017-12-15 | 浙江工业大学 | 一种基于快速非奇异终端滑模的四旋翼无人机有限时间自适应控制方法 |
CN107490966A (zh) * | 2017-08-22 | 2017-12-19 | 浙江工业大学 | 一种基于改进幂次趋近律的飞行器有限时间自适应姿态控制方法 |
Non-Patent Citations (3)
Title |
---|
HAIBODU,ET AL.: "Finite-time formation control for a group of quadrotor aircraft", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
WENWU ZHU,ET AL.: "Hovering control for quadrotor aircraft based on finite-time control algorithm", 《NONLINEARDYNAMICS》 * |
朱文武等: "基于有限时间姿态控制算法的四旋翼飞行器悬停控制", 《系统科学与数学》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885074A (zh) * | 2019-02-28 | 2019-06-14 | 天津大学 | 四旋翼无人机有限时间收敛姿态控制方法 |
CN109885074B (zh) * | 2019-02-28 | 2022-02-15 | 天津大学 | 四旋翼无人机有限时间收敛姿态控制方法 |
WO2020177267A1 (zh) * | 2019-03-07 | 2020-09-10 | 中国科学院深圳先进技术研究院 | 四旋翼无人机的控制方法、装置、设备及可读介质 |
CN110032205A (zh) * | 2019-04-29 | 2019-07-19 | 河海大学常州校区 | 一种具有抗干扰能力的无人机姿态控制方法 |
CN110032205B (zh) * | 2019-04-29 | 2021-09-28 | 河海大学常州校区 | 一种具有抗干扰能力的无人机姿态控制方法 |
CN111984027A (zh) * | 2019-05-23 | 2020-11-24 | 南京航空航天大学 | 异构多无人机一致性控制方法 |
CN111984027B (zh) * | 2019-05-23 | 2022-07-15 | 南京航空航天大学 | 异构多无人机一致性控制方法 |
CN110275542B (zh) * | 2019-06-14 | 2022-04-08 | 合肥工业大学 | 一种基于自适应有限时间控制的四旋翼飞行器控制方法 |
CN110275542A (zh) * | 2019-06-14 | 2019-09-24 | 合肥工业大学 | 一种基于自适应有限时间控制的四旋翼飞行器控制方法 |
CN110673623A (zh) * | 2019-11-04 | 2020-01-10 | 苏州经贸职业技术学院 | 一种基于双环pd控制算法控制的四旋翼无人机着陆方法 |
CN110673623B (zh) * | 2019-11-04 | 2022-11-22 | 苏州经贸职业技术学院 | 一种基于双环pd控制算法控制的四旋翼无人机着陆方法 |
CN111580537A (zh) * | 2020-05-28 | 2020-08-25 | 西北工业大学 | 一种无人机特技飞行控制系统及方法 |
CN113759706B (zh) * | 2021-11-05 | 2022-02-08 | 普宙科技(深圳)有限公司 | 基于轴角矢量的多旋翼姿态自抗扰控制方法及系统 |
CN113759706A (zh) * | 2021-11-05 | 2021-12-07 | 普宙科技(深圳)有限公司 | 基于轴角矢量的多旋翼姿态自抗扰控制方法及系统 |
CN114995483A (zh) * | 2022-06-28 | 2022-09-02 | 四川大学 | 基于对偶四元数的倾转旋翼多旋翼无人机跟踪控制方法 |
CN115309172A (zh) * | 2022-07-01 | 2022-11-08 | 中国人民解放军国防科技大学 | 一种用于运输血液的无人机的姿态控制方法 |
CN117055602A (zh) * | 2023-09-05 | 2023-11-14 | 广东工业大学 | 一种四旋翼无人机的全姿态控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108363298A (zh) | 一种基于四元数描述的四旋翼无人机快速收敛控制方法 | |
CN107479567B (zh) | 动态特性未知的四旋翼无人机姿态控制器及方法 | |
CN106094855B (zh) | 一种四旋翼无人机的终端协同控制方法 | |
He et al. | A Simple Attitude Control of Quadrotor Helicopter Based on Ziegler‐Nichols Rules for Tuning PD Parameters | |
CN109884895A (zh) | 基于饱和受限情况下的无人机自适应跟踪控制算法 | |
Kuantama et al. | Feedback linearization LQR control for quadcopter position tracking | |
Li et al. | Robust tracking control strategy for a quadrotor using RPD-SMC and RISE | |
CN107368088B (zh) | 一种基于误差指数型函数的四旋翼飞行器非线性滑模位姿控制方法 | |
CN107561931B (zh) | 一种基于单指数型函数的四旋翼飞行器非线性滑模位姿控制方法 | |
CN111258216B (zh) | 一种适用于四旋翼飞行器的滑模重复控制器 | |
CN107368089B (zh) | 一种基于双指数型函数的四旋翼飞行器非线性滑模位姿控制方法 | |
CN106707754A (zh) | 一种基于切换系统的货物搬运旋翼无人机建模及自适应控制方法 | |
Jiao et al. | Analysis and design the controller for quadrotors based on PID control method | |
CN107992082A (zh) | 基于分数阶幂次切换律的四旋翼无人机飞行控制方法 | |
CN107263483B (zh) | 二自由度关节机器人轨迹的协调控制方法 | |
CN109885074A (zh) | 四旋翼无人机有限时间收敛姿态控制方法 | |
CN109283932A (zh) | 一种基于积分反步滑模的四旋翼飞行器姿态控制方法 | |
Zhu et al. | Back‐stepping sliding mode control method for quadrotor UAV with actuator failure | |
CN115639830B (zh) | 一种空地智能体协同编队控制系统及其编队控制方法 | |
CN110333733A (zh) | 一种四旋翼飞行器的串级变论域模糊pid姿态控制系统及方法 | |
CN110275542A (zh) | 一种基于自适应有限时间控制的四旋翼飞行器控制方法 | |
CN108438256A (zh) | 一种基于永磁动量交换球的对地凝视卫星姿态控制方法 | |
CN109976364B (zh) | 一种六旋翼飞行器姿态解耦控制方法 | |
CN107977011A (zh) | 基于分数阶控制算法的四旋翼无人机飞行控制方法 | |
El Houm et al. | Optimal new sliding mode controller combined with modified supertwisting algorithm for a perturbed quadrotor UAV |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180803 |
|
RJ01 | Rejection of invention patent application after publication |