基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法技术

一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法，包括以下步骤：步骤1，确定从基于四旋翼飞行器的机体坐标系到基于地球的惯性坐标系的转移矩阵；步骤2，根据牛顿欧拉公式分析四旋翼飞行器动力学模型；步骤3，计算跟踪误差，根据快速终端滑模面以及其一阶导数设计控制器。本发明专利技术针对四旋翼飞行器系统，结合基于对数增强型指数趋近律滑模控制以及快速终端滑模控制，不但能在远离滑模面时能增加趋近速度，并且能减小抖振，提高系统的快速性和鲁棒性，实现快速稳定控制，同时能实现跟踪误差的有限时间控制，解决了传统滑模面中只有当时间趋于无穷，跟踪误差才趋向0的问题。

【技术实现步骤摘要】
基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法
本专利技术涉及一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法。
技术介绍
四旋翼飞行器由于结构简单、机动性强、飞行方式独特的特点引起了国内外学者以及科研机构的广泛关注，并迅速成为目前国际上研究的热点之一。相比固定翼飞行器，旋翼飞行器可以垂直升降，对环境要求低，不需要跑道，降低了成本，有着巨大的商业价值。飞行器的发展使许多危险的高空作业变得轻松安全，在军事方面给其他国家造成威慑，在民用方面使工作效率大大增加。四旋翼飞行器具有较强的灵活性，可随时实现运动和悬停的快速过渡，并且能以较小的损坏风险胜任更具挑战性的飞行任务。在科学研究领域，由于四旋翼飞行器具有非线性、欠驱动、强耦合的动态特性，研究人员常将其作为理论研究、方法验证的实验载体。依托小型四旋翼飞行器，搭建飞行器飞行控制系统，进行飞行器高性能运动控制研究，是当前学术界的热点研究领域。趋近律滑模控制的特点是可以实现不连续控制，滑动模态是可设计的，而且与系统参数及扰动没有关联。趋近律滑模不仅可以合理设计到达滑模面的速度，减小趋近阶段的时间，提高系统的鲁棒性，而且能有效的减弱滑模控制中的抖振问题。目前，在四旋翼控制领域中，使用趋近律滑模控制比较少。增强型趋近律在传统趋近律的基础上进一步加快了系统到达滑模面的趋近速度同时使得抖振更小。
技术实现思路
为了克服传统滑模面无法实现有限时间控制以及进一步加快趋近律的趋近速度和减小抖振的问题，本专利技术采用了快速终端滑模控制以及基于对数增强型指数趋近律，通过切换控制的思想避免了奇异性问题，加快了系统到达滑模面的趋近速度，减小了抖振，实现了有限时间控制。为了解决上述技术问题提出的技术方案如下：一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法，包括以下步骤：步骤1，确定从基于四旋翼飞行器的机体坐标系到基于地球的惯性坐标系的转移矩阵；其中ψ、θ、φ分别是飞行器的偏航角、俯仰角、翻滚角，表示飞行器绕依次惯性坐标系各轴旋转的角度，Tψ表示ψ的转移矩阵，Tθ表示θ的转移矩阵，Tφ表示φ的转移矩阵；步骤2，根据牛顿欧拉公式分析四旋翼飞行器动力学模型，过程如下：2.1，平动过程中有：其中x、y、z分别表示四旋翼在惯性坐标系下的位置，m表示飞行器的质量，g表示重力加速度，mg表示四旋翼所受重力，四个旋翼产生的合力Ur；2.2，转动过程中有：其中τx、τy、τz分别代表机体坐标系上的各轴力矩分量，Ixx、Iyy、Izz分别代表机体坐标系上的各轴转动惯量分量，×表示叉乘，wp、wq、wr分别代表机体坐标系上的各轴姿态角速度分量，分别代表机体坐标系上的各轴姿态角加速度分量；考虑到飞行器处于低速飞行或者悬停状态下，姿态角变化较小，认为则转动过程中式(3)表示为式(4)2.3，联立式(1)，(2)，(4)，得飞行器的动力学模型如式(5)所示其中Ux、Uy、Uz分别为三个位置控制器的输入量；根据式(5)，对位置姿态关系进行解耦计算，结果如下：其中φd为φ的期望信号值，θd为θ的期望信号值，ψd为ψ的期望信号值，arcsin函数是反正弦函数，arctan函数是反正切函数；式(5)也可写成矩阵形式，如下：其中X1＝[x,y,z,φ,θ,ψ]T,B(X)＝diag(1,1,1,b1,b2,b3),U＝[Ux,Uy,Uz,τx,τy,τz]T步骤3，计算跟踪误差，根据快速终端滑模面以及其一阶导数设计控制器，过程如下：3.1，定义跟踪误差及其一阶微分和二阶微分：e＝X1-Xd(8)其中，Xd＝[xd,yd,zd,φd,θd,ψd]T，xd,yd,zd,φd,θd,ψd分别为x,y,z,φ,θ,ψ的可导期望信号；3.2，设计快速终端滑模面：其中，sigα(x)＝|x|α·sign(x)，α1＞α2＞1，λ1＞0，λ2＞0；对式(11)进行求导，得到：令式(12)简化为式(13)但由于α(e)中存在的负幂次项，当α(e)＝0且β(e)≠0会导致奇异性问题；考虑切换控制的方法：其中qi(e),αi(e),βi(e)分别为q(e),α(e),β(e)的第i个元素，i＝1,2,3,4,5,6；联立式(13)和式(14)，得到：联立式(7)、式(10)和式(15)，得到：3.3，设计增强型趋近律其中N-1(X)为N(X)的逆矩阵，k1＞0，k2＞0，0＜δ＜1，γ＞0，μ＞1，p为正整数；3.4，联立式(16)和式(17)，得到控制器其中B-1(X)为B(X)的逆矩阵。进一步，所述控制方法还包括以下步骤：步骤4，性质说明，过程如下：4.1，证明滑模的可达性：设计李雅普诺夫函数对此函数两边进行求导，得：因为恒大于0，所以式(18)恒小于0，满足滑模的可达性，系统能到达滑模面；4.2，增强型效果说明：系统远离滑模面时即|s|很大，N(s)趋近δ,系统的趋近速度加快；当系统接近滑模面时即|s|趋近0，N(s)趋近μ，系统的抖振减小。本专利技术的技术构思为：针对四旋翼飞行器系统，结合指数趋近律滑模控制以及快速终端滑模控制，设计了一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法。快速终端滑模面能实现跟踪误差的有限时间控制，解决了传统滑模面中时间趋于无穷，误差才趋向0的问题。基于对数增强型趋近律不但能在远离滑模面时能增加趋近速度，并且能减小抖振，提高系统的快速性和鲁棒性，实现快速稳定控制。本专利技术的有益效果为：增强了系统的鲁棒性，与传统的指数趋近律滑模控制相比，不但能在远离滑模面时能增加趋近速度，并且能减小抖振，缩短了滑动模态的到达时间，从而使系统更快地实现稳定收敛。除此之外，本专利技术利用快速终端滑模，解决了传统滑模面中时间趋于无穷，误差才趋向0的问题，实现了有限时间控制。附图说明图1为四旋翼飞行器的位置跟踪效果示意图，其中点线代表传统指数趋近律控制，虚线代表基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制。图2为四旋翼飞行器的位置跟踪误差示意图，其中点线代表传统指数趋近律控制，虚线代表基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制。图3为四旋翼飞行器的姿态角跟踪效果示意图，其中点线代表传统指数趋近律控制，虚线代表基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制。图4为四旋翼飞行器的姿态角跟踪误差示意图，其中点线代表传统指数趋近律控制，虚线代表基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制。图5为基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制下的位置控制器输入示意图。图6为四旋翼飞行器的传统指数趋近律控制下的位置控制器输入示意图。图7为基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制下的姿态角控制器输入示意图。图8为四旋翼飞行器的传统指数趋近律控制下的姿态角控制器输入示意图。图9为本专利技术的控制流程示意图。具体实施方式下面结合附图对本专利技术做进一步说明。参照图1-图9，一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法，包括以下步骤：步骤1，确定从基于四旋翼飞行器的机体坐标系到基于地球的惯性坐标系的转移矩阵本文档来自技高网...

【技术保护点】
1.一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法，其特征在于，所述控制方法包括以下步骤：步骤1，确定从基于四旋翼飞行器的机体坐标系到基于地球的惯性坐标系的转移矩阵；

【技术特征摘要】
1.一种基于对数增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法，其特征在于，所述控制方法包括以下步骤：步骤1，确定从基于四旋翼飞行器的机体坐标系到基于地球的惯性坐标系的转移矩阵；其中ψ、θ、φ分别是飞行器的偏航角、俯仰角、翻滚角，表示飞行器绕依次惯性坐标系各轴旋转的角度，Tψ表示ψ的转移矩阵，Tθ表示θ的转移矩阵，Tφ表示φ的转移矩阵；步骤2，根据牛顿欧拉公式分析四旋翼飞行器动力学模型，过程如下：2.1，平动过程中有：其中x、y、z分别表示四旋翼在惯性坐标系下的位置，m表示飞行器的质量，g表示重力加速度，mg表示四旋翼所受重力，四个旋翼产生的合力Ur；2.2，转动过程中有：其中τx、τy、τz分别代表机体坐标系上的各轴力矩分量，Ixx、Iyy、Izz分别代表机体坐标系上的各轴转动惯量分量，×表示叉乘，wp、wq、wr分别代表机体坐标系上的各轴姿态角速度分量，分别代表机体坐标系上的各轴姿态角加速度分量；考虑到飞行器处于低速飞行或者悬停状态下，认为则转动过程中式(3)表示为式(4)2.3，联立式(1)，(2)，(4)，得飞行器的动力学模型如式(5)所示其中Ux、Uy、Uz分别为三个位置控制器的输入量；根据式(5)，对位置姿态关系进行解耦计算，结果如下：其中φd为φ的期望信号值，θd为θ的期望信号值，ψd为ψ的期望信号值，arcsin函数是反正弦函数，arctan函数是反正切函数；式(5)也可写成矩阵形式，如下：其中X1＝[x,y,z,φ,θ,ψ]T,B(X)＝diag(1,1,1,b1,b2,b3),U＝[Ux,Uy,Uz,τx,τy,τz]T步骤3，...

【专利技术属性】
技术研发人员：陈强，陈凯杰，陶玫玲，胡轶，吴春，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江,33