CN112987567B - 非线性系统的固定时间自适应神经网络滑模控制方法 - Google Patents

非线性系统的固定时间自适应神经网络滑模控制方法 Download PDF

Info

Publication number
CN112987567B
CN112987567B CN202110175892.9A CN202110175892A CN112987567B CN 112987567 B CN112987567 B CN 112987567B CN 202110175892 A CN202110175892 A CN 202110175892A CN 112987567 B CN112987567 B CN 112987567B
Authority
CN
China
Prior art keywords
rotor
neural network
control
fixed time
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
CN202110175892.9A
Other languages
English (en)
Other versions
CN112987567A (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.)
Hebei University of Science and Technology
Original Assignee
Hebei University of Science and Technology
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 Hebei University of Science and Technology filed Critical Hebei University of Science and Technology
Priority to CN202110175892.9A priority Critical patent/CN112987567B/zh
Publication of CN112987567A publication Critical patent/CN112987567A/zh
Application granted granted Critical
Publication of CN112987567B publication Critical patent/CN112987567B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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

Abstract

本发明公开了一种非线性系统的固定时间自适应神经网络滑模控制方法,包括以下步骤:步骤S1,根据四轴飞行器的动态运动规律,建立四轴飞行器的动力学模型;步骤S2,基于非奇异快速终端滑模面提出双幂次固定时间控制律,以实现高鲁棒性和快速滑动率;步骤S3,仿真验证,以验证所提出的针对不确定非线性系统的固定时间自适应神经网络自适应控制的有效性;步骤S4,数值算例验证,使用Matlab仿真来验证固定时间自适应神经网络控制的有效性。基于非奇异快速终端滑模面设计了双幂次固定时间控制,以实现高鲁棒性和快速滑动率,并且在控制律中没有负指数项以有效避免奇异现象。本发明提出的算法可以大大提高系统的鲁棒性能,具有较高的实用价值。

Description

非线性系统的固定时间自适应神经网络滑模控制方法
技术领域
本发明涉及一种非线性系统的固定时间自适应神经网络滑模控制方法,属于无人机的智能控制技术领域。
背景技术
几十年来,由于其固有的特性,四旋翼直升机一直受到军事,民用和工程学学者的广泛关注。小型无人机(例如四旋翼飞机)的应用范围非常广泛,例如在恶劣环境下的军事侦察,民用物流,航空摄影以及农药喷洒。为了提高飞机在各种情况下的稳定性和可靠性,学者们在四旋翼无人机的智能控制研究中取得了很多研究成果。四旋翼已有众多控制策略,但由于四旋翼系统为复杂欠驱动系统,传统PID算法、滑模控制并不能满足实际控制需求。
PID的缺陷,概括起来就是信号处理太简单、未能充分发挥其优点,具体说来,有四个方面:
(1)产生误差的方式不太合理控制目标v在过程中可以“跳变”,但是被控对象输出Y的变化都有惯性,不可能跳变,要求让缓变的变量y来跟踪能够跳变的变量v,初始误差很大,易引起超调,很不合理。
(2)误差的微分信号的产生没有太好的办法由于微分器物理不可实现,只能近似实现,常用的近似微分器的形式为
(3)误差积分反馈的引入有很多负作用在PID控制中,误差积分反馈的作用是消除静差,提高系统响应的准确性,但同时误差积分反馈的引入,使闭环变得迟钝,容易产生振荡,易产生由积分饱和引起的控制量饱和。
(4)线性组合不一定是最好的组合方式PID控制器给出的控制量是误差的现在、过去、将来三者的线性组合。大量工程实践表明,线性组合不一定是最好的组合方式,能否在非线性领域找到更合适的组合方式是值得探索的。
滑模控制的缺点:当状态轨迹到达滑动模态面后,难以严格沿着滑动模态面向平衡点滑动,而是在其两侧来回穿越地趋近平衡点,从而产生抖振——滑模控制实际应用中的主要障碍。滑模控制的优点是能够克服系统的不确定性,对干扰和未建模动态具有很强的鲁棒性,尤其是对非线性系统的控制具有良好的控制效果。
发明内容
本发明所要解决的技术问题是提供一种非线性系统的固定时间自适应神经网络滑模控制方法,应用的算法能够有效解决非线性系统的不确定性,且具有优越的鲁棒性能;同时,通过设计固定时间神经网络自适应律,使系统可以达到固定时间收敛,具有更高的实用价值。
为解决上述问题,本发明所采取的技术方案是:
一种非线性系统的固定时间自适应神经网络滑模控制方法,包括以下步骤:
步骤S1,根据四轴飞行器的动态运动规律,建立四轴飞行器的动力学模型;
步骤S2,基于非奇异快速终端滑模面提出双幂次固定时间控制律,以实现高鲁棒性和快速滑动率;
步骤S3,仿真验证,以验证所提出的针对不确定非线性系统的固定时间自适应神经网络自适应控制的有效性;
步骤S4,数值算例验证,使用Matlab仿真来验证固定时间自适应神经网络控制的有效性。
作为本发明的进一步改进,步骤S1中所述动力学模型选择身体坐标系和地面坐标系,根据坐标系B和坐标系E的空间变换,利用牛顿-欧拉方程,计算四旋翼飞行器的数学模型。
作为本发明的进一步改进,步骤S2中,所述双幂次固定时间控制律的获取过程如下:考虑到以下动态系统:
Figure BDA0002940699430000021
可以改写为
Figure BDA0002940699430000022
其中y为统输出,dt为不确定干扰;
为使x能够跟踪xd定义跟踪误差
e1=x1-xd (54)
Figure BDA0002940699430000023
选取非奇异快速终端滑模面
Figure BDA0002940699430000031
其中α>0,β>0,1<p2<2,p1>p2,
Figure BDA0002940699430000032
g,h,p,q∈N均为奇数;
Figure BDA0002940699430000033
为使达到估计滑膜面,在不考虑干扰条件下,设计等效控制律
Figure BDA0002940699430000034
采用双幂指数切换控制率
Figure BDA0002940699430000035
其中k1>0,k2>0,且为系数。最终双幂次固定时间控制律为
u=ueq+usw (60)
Figure BDA0002940699430000036
然而,二阶非线性系统数学模型中存在未知的f=cos(x1)+dt,控制律(16)并不能完全的达到优越的控制效果。
作为本发明的进一步改进,为了逼近非线性系统中的不确定非线性部分,使用RBF设计神经网络控制。
作为本发明的进一步改进,步骤S3中,所述仿真验证过程如下:
参考轨迹为
xd=sin(t) (62)
控制器以及网络参数选择为
Figure BDA0002940699430000037
根据引理1可以得到固定时间Txmax=5.742;
系统的初始化条件以及假设存在的干扰dt
Figure BDA0002940699430000041
作为本发明的进一步改进,为了验证算法应用于实际系统模型的有效性,通过数值算例针对四旋翼飞行器跟踪控制进行有效性以及抗干扰性能的验证。
作为本发明的进一步改进,步骤S4中,所述数值算例验证过程如下:
考虑无人机动态系统(1),其参数为
m=2,l=0.2,g=9.8
ξx=ξx=ξx=1.2
ξφ=ξθ=ξx=1.2
Ix=1.25,Iy=1.25,Iz=2.5
参考轨迹的选择如下:
Figure BDA0002940699430000042
选择系统的初始条件为
x(0)=y(0)=z(0)=φ(0)=θ(0)=ψ(0)=0.5
然后可以给出控制器和神经网络自适应律为
Figure BDA0002940699430000043
Figure BDA0002940699430000044
其中ux是四旋翼位置虚拟控制器,
Figure BDA0002940699430000045
是固定时间自适应律,
Figure BDA0002940699430000046
是神经网络输出其中位置子系统控制器以及网络参数为
k1=3000,k2=0.1
α=50,β=3,
Figure BDA0002940699430000047
b=20,
Figure BDA0002940699430000048
σz=σx=1
根据引理1可以得到固定时间Txmax=Tymax=Tzmax=5.742;
本采用姿态解算,得到四旋翼飞行器的目标姿态角,假设期望
Figure BDA0002940699430000051
经过位置子系统的姿态计算,得到φd,θd,从而实现跟踪控制;
Qx=u1(sinψsinφ+cosψsinθcosφ)
Qy=u1(-cosψsinφ+sinψsinθcosφ)
Qz=u1(cosθcosφ)
Figure BDA0002940699430000052
Figure BDA0002940699430000053
其中
Figure BDA0002940699430000054
是四旋翼位置控制输入,Qx,Qy,Qz是虚拟控制输入。
然后可以给出控制器和神经网络自适应律为
Figure BDA0002940699430000055
Figure BDA0002940699430000056
Figure BDA0002940699430000057
Figure BDA0002940699430000058
Figure BDA0002940699430000059
Figure BDA00029406994300000510
其中uφ,uθ,uψ是四旋翼姿态控制器,
Figure BDA00029406994300000511
是固定时间自适应律,
Figure BDA00029406994300000512
是神经网络输出;姿态子系统控制器以及网络参数
k1=k3=3000,k2=k4=0.1
k5=k6=1,α3=1,β3=1,
Figure BDA0002940699430000061
α1=α2=50,β1=β2=3,
Figure BDA0002940699430000062
b=20,
Figure BDA0002940699430000063
σz=σx=1
根据引理1可以得到固定时间Tφmax=Tθmax=5.742,Tψmax=4.3784。
采用上述技术方案所产生的有益效果在于:
(1)设计基于非奇异快速终端滑模面的双幂次滑模控制律,能够有效解决传统滑模控制中存在的抖振和收敛速度慢的问题,使系统可以达到在有限时间内稳定,且解决了可能存在的奇异问题。
(2)设计固定时间神经网络自适应律,在不考虑理想权重和权重初始值的条件下,可以在有限时间内逼近不确定非线性系统。
(3)提供了精确的数学模型,算法能够保证在较大干扰条件下稳定,将所提算法在四旋翼飞行器算例验证,结果表明四旋翼飞行器具有非常强的鲁棒性能。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是现有技术中PID调节器的控制原理框图;
图2是四旋翼系统动态模型示意图;
图3是不确定非线性系统控制流程图;
图4是径向基函数神经网络结构图;
图5是仿真验证中位置轨迹与理想位置的误差轨迹示意图;
图6是仿真验证中速度轨迹与理想速度的误差轨迹示意图;
图7是神经网络近似轨迹图;
图8是四旋翼控制系统结构图;
图9是数值算例验证中位置轨迹x与理想位置xd的误差轨迹示意图;
图10是数值算例验证中位置轨迹y与理想位置yd的误差轨迹示意图;
图11是数值算例验证中位置轨迹z与理想位置zd的误差轨迹示意图;
图12是数值算例验证中实际位置与理想位置的3D效果示意图;
图13是数值算例验证中实际姿态φ与理想姿态φd的轨迹示意图;
图14是数值算例验证中实际姿态θ与理想姿态θd的轨迹示意图;
图15是数值算例验证中实际姿态ψ与理想姿态ψd的轨迹示意图;
图16是位置子系统网络近似轨迹示意图;
图17是姿态子系统网络近似轨迹示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的,决不作为对本申请及其应用或使用的任何限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示,PID调节器是一种线性调节器,这种调节器是将设定值r(t)与输出值c(t)进行比较构成控制偏差。
e(t)=r(t)-c(t)
将其按照比例、积分、微分运算后,并通过线性组合构成控制量,如图所示,所以简称P(比例)、I(积分)、D(微分)调节器。
应用PID算法对四旋翼飞行器的动态性能进行了算例验证,能够使达到稳定的控制,但PID控制策略在参数计算上具有很高的复杂性,应用的控制模型并不够精确,四旋翼飞行器难以在复杂环境下稳定运行。
滑模变结构控制的原理,是根据系统所期望的动态特性来设计系统的切换超面,通过滑动模态控制器使系统状态从超平面之外向切换超平面收束。系统一旦到达切换超平面,控制作用将保证系统沿切换超平面到达系统原点,这一沿切换超平面向原点滑动的过程称为滑模控制。
为了更好的解决非线性系统的控制问题,本发明所应用的算法能够有效解决非线性系统的不确定性,且具有优越的鲁棒性能。同时,通过设计固定时间神经网络自适应律,使系统可以达到固定时间收敛,具有更高的实用价值。经过算例验证能够有效应用于四旋翼跟踪控制。
根据四轴飞行器的动态运动规律,建立了四轴飞行器的动力学模型,并选择了身体坐标系和地面坐标系,如图2所示。考虑到如图2所示四旋翼无人机系统,选择身体坐标系以及地面坐标系,根据坐标系和系统动力学模型,利用牛顿欧拉方程,得到四旋翼系统的数学模型。姿态角表示为偏航角ψ,俯仰角θ,侧倾角φ。假设四旋翼为刚体,它的动态模型可以表示为:
总而言之,可以得到四轴飞行器位置和姿态动力学模型:
Figure BDA0002940699430000081
其中,字母上方带点表示所代表含义的导数,字母上方两个点表示所代表含义的二阶导数,m是四旋翼的质量,I=diag(Ix,Iy,Iz)是身体坐标系下三个坐标轴的惯性矩,g是所选取得重力加速度,wi,i=1,2,3,4是身体坐标系下的旋转角速度,ξ=diag(ξxyzφθψ)是空气阻力系数,Ir是转子惯性,
Figure BDA0002940699430000082
是总的剩余转子角度,da(·)=diag[dx dy dz]和dp(·)=diag[dφ dθ dψ]是位置和姿态系统中不确定干扰,Fat表示三个方向上的控制推力的合力,τ是转子产生的扭矩。
引理1假设V(·):Rn→R+∪{0}是一个连续的根本无界函数,并且满足以下两个条件:
Figure BDA0002940699430000091
Figure BDA0002940699430000092
其中a,b,p,q是Lyapunov函数V(x)的正实数系数且p∈(0,1),q∈(1,∞),0<c<∞。然后系统
Figure BDA0002940699430000093
的初始x=0则系统是几乎固定时间稳定。此外,具有以下不等式成立
V(x,t)≤ξ,t≥Tmax (68)
Figure BDA0002940699430000094
其中ξ是方程根,Tmax是最大收敛时间。
引理2为杨氏不等式,对于任意常数
Figure BDA0002940699430000095
可以得到不等式
Figure BDA0002940699430000096
其中pa>1,qa>1且
Figure BDA0002940699430000097
引理3。对于正常数
Figure BDA0002940699430000098
向量W*,
Figure BDA0002940699430000099
且满足
Figure BDA00029406994300000910
可以得到不等式
Figure BDA00029406994300000911
Figure BDA00029406994300000912
是理想权重,
Figure BDA00029406994300000913
是估计权重,
Figure BDA00029406994300000914
是近似误差,
Figure BDA00029406994300000915
控制器的不确定非线性系统控制流程如图3所示,
考虑到以下动态系统:
Figure BDA00029406994300000916
dt为不确定干扰。
为使x能够跟踪xd定义跟踪误差
e1=x1-xd (73)
Figure BDA0002940699430000101
选取非奇异快速终端滑模面
Figure BDA0002940699430000102
其中α>0,β>0,1<p2<2,p1>p2,
Figure BDA0002940699430000103
g,h,p,q∈N均为奇常数。
Figure BDA0002940699430000104
为使达到估计滑膜面,在不考虑干扰条件下,设计等效控制律
Figure BDA0002940699430000105
我们采用双幂指数切换控制率
Figure BDA0002940699430000106
其中k1>0,k2>0且为常数,最终控制律为
u=ueq+usw (79)
Figure BDA0002940699430000107
然而,二阶非线性系统数学模型中存在未知的f=cos(x1)+dt,所以控制律(16)并不能完全的达到优越的控制效果。因此,本实施例将采用RBF神经网络对非线性部分进行逼近,能够有效增强非线性系统的鲁棒性能。
如图4所示,RBF神经网络作为一种前馈型神经网络具有最佳的逼近效果,且不存在局部极小问题和学习收敛速度快等优势。因此,早已应用于大量学术研究之中,得到了有效的验证。
为了逼近非线性系统中的不确定非线性部分,使用RBFNNs设计神经网络控制。
FNN(x,W)=WΨ(x) (81)
Figure BDA0002940699430000108
其中
Figure BDA0002940699430000109
Ψ(x)=diag[ψ1(x) …ψl(x)]-节点向量,ψi(x)-高斯函数,μi=diag[μi1…μin]-基函数中心,ηi-高斯函数的标量宽度。
RNF神经网络可用于在一个紧急集合
Figure BDA0002940699430000111
上近似任何连续函数
F(x)=W*Ψ(x)+ε (83)
其中ε-神经网咯的近似误差,W*-理想权重
Figure BDA0002940699430000112
Figure BDA0002940699430000113
是估计权重和近似误差
Figure BDA0002940699430000114
选择神经网络控制器
Figure BDA0002940699430000115
其中
Figure BDA0002940699430000116
非线性系统滑模面可以写成
Figure BDA0002940699430000117
其中
f*=f+dt=W*Ψ(x)+ε (87)
计算系统模型的的近似误差
Figure BDA0002940699430000118
根据控制器(21)和滑模面(22)可以得到
Figure BDA0002940699430000119
其中k1>0,k2>0且为常数。
为验证滑模面稳定性,选择Lyapunov候选函数
Figure BDA00029406994300001110
Figure BDA00029406994300001111
未验证网络稳定,选择Lyapunov候选函数
Figure BDA0002940699430000121
Figure BDA0002940699430000122
Figure BDA0002940699430000123
其中ξ,σzx为大于0的常数。
为验证整体系统的稳定,选择Lyapunov函数
V=Vs+VNN (95)
Figure BDA0002940699430000124
Figure BDA0002940699430000125
根据引理2和引理3可有以下不等式成立
Figure BDA0002940699430000126
Figure BDA0002940699430000127
Figure BDA0002940699430000128
其中σ1>0,σ2>0,σ3>0,σ4>0且为常数,然后我们有
Figure BDA0002940699430000131
Figure BDA0002940699430000132
然后我们可以得到
Figure BDA0002940699430000133
其中
Figure BDA0002940699430000134
然后,根据引理1系统近似固定时间稳定,收敛时间为
Figure BDA0002940699430000135
其中σxz取决于ν,γ,
Figure BDA0002940699430000136
Γ<0。
仿真验证
在本节中,进行仿真以验证所提出的针对不确定非线性系统的固定时间自适应神经网络自适应控制的有效性。
参考轨迹为
xd=sin(t) (105)
控制器以及网络参数选择为
Figure BDA0002940699430000141
根据引理1可以得到固定时间Txmax=5.742。
系统的初始化条件以及假设存在的干扰dt
Figure BDA0002940699430000142
根据图5-7,能够验证固定时间自适应神经网络算法的有效性。为了验证算法应用于实际系统模型的有效性,下面将通过数值算例针对四旋翼飞行器跟踪控制进行有效性的验证。
四旋翼控制系统结构如图8所示,
四旋翼飞行器位置动力学模型:
Figure BDA0002940699430000143
其中
Figure BDA0002940699430000144
其中A是四旋翼质量对角矩阵,f1(·)是四旋翼位置函数矩阵,us(t)是四旋翼位置虚拟控制输入矩阵;
期望系统位置输出x能够跟踪xd,就有
z1=x-xd (109)
Figure BDA0002940699430000145
选取非奇异快速终端滑模面
Figure BDA0002940699430000146
其中α>0,β>0,1<p2<2,p1>p2,
Figure BDA0002940699430000151
g,h,p,q∈N均为奇数。
Figure BDA0002940699430000152
为使达到估计滑膜面,在不考虑干扰条件下,设计等效控制律
Figure BDA0002940699430000153
我们采用双幂指数切换控制率
Figure BDA0002940699430000154
其中k1>0,k2>0且为常数,最终控制律和固定时间自适应律为
us=ueq+usw (115)
Figure BDA0002940699430000155
Figure BDA0002940699430000156
其中ξ,σzx为大于0的常数。
根据(25)我们可以同理得到
Figure BDA0002940699430000157
为验证系统的稳定性,选取Lyapunov函数
Figure BDA0002940699430000158
Figure BDA0002940699430000159
Figure BDA00029406994300001510
根据引理2和引理3可有以下不等式成立
Figure BDA0002940699430000161
Figure BDA0002940699430000162
Figure BDA0002940699430000163
其中σ1>0,σ2>0,σ3>0,σ4>0且为常数然后我们有
Figure BDA0002940699430000164
Figure BDA0002940699430000165
然后我们可以得到
Figure BDA0002940699430000166
其中
Figure BDA0002940699430000167
然后,根据引理1系统近似固定时间稳定,最大收敛时间为
Figure BDA0002940699430000168
其中σxz取决于ν,γ,
Figure BDA0002940699430000171
Γ<0。
四旋翼飞行器姿态动力学模型:
Figure BDA0002940699430000172
其中
Figure BDA0002940699430000173
其中B是四旋翼姿态惯性矩矩阵,f2(·)是四旋翼姿态函数矩阵,ur(t)是四旋翼控制输入矩阵;
期望系统姿态p能够跟踪pd,从而有
z3=p-pd (130)
Figure BDA0002940699430000174
选取非奇异快速终端滑模面
Figure BDA0002940699430000175
其中α>0,β>0,1<p2<2,p1>p2,
Figure BDA0002940699430000176
g,h,p,q∈N均为奇数。
Figure BDA0002940699430000177
为使达到估计滑膜面,在不考虑干扰条件下,设计等效控制律
Figure BDA0002940699430000178
我们采用双幂指数切换控制率
Figure BDA0002940699430000179
其中k1>0,k2>0,最终控制律为
us=ueq+usw (136)
Figure BDA0002940699430000181
同理可证得姿态子系统固定时间稳定。
数值算例验证
在本节中,将使用Matlab仿真来验证固定时间自适应神经网络控制的有效性。
示例:考虑无人机动态系统(1),其参数为
m=2,l=0.2,g=9.8
ξx=ξx=ξx=1.2
ξφ=ξθ=ξx=1.2
Ix=1.25,Iy=1.25,Iz=2.5
参考轨迹的选择如下:
Figure BDA0002940699430000182
选择系统的初始条件为
x(0)=y(0)=z(0)=φ(0)=θ(0)=ψ(0)=0.5
然后可以给出控制器和神经网络自适应律为
Figure BDA0002940699430000183
Figure BDA0002940699430000184
其中ux是四旋翼位置虚拟控制器,
Figure BDA0002940699430000185
是固定时间自适应律,
Figure BDA0002940699430000186
是神经网络输出。
其中位置子系统控制器以及网络参数为
k1=3000,k2=0.1
α=50,β=3,
Figure BDA0002940699430000187
b=20,
Figure BDA0002940699430000188
σz=σx=1
根据引理1可以得到固定时间Txmax=Tymax=Tzmax=5.742。
本文采用姿态解算,得到四旋翼飞行器的目标姿态角,假设期望
Figure BDA0002940699430000191
经过位置子系统的姿态计算,可以得到φd,θd,从而实现跟踪控制。
Qx=u1(sinψsinφ+cosψsinθcosφ)
Qy=u1(-cosψsinφ+sinψsinθcosφ)
Qz=u1(cosθcosφ)
Figure BDA0002940699430000192
Figure BDA0002940699430000193
其中
Figure BDA0002940699430000194
是四旋翼位置控制输入,Qx,Qy,Qz是虚拟控制输入。
然后可以给出控制器和神经网络自适应律为
Figure BDA0002940699430000195
Figure BDA0002940699430000196
Figure BDA0002940699430000197
Figure BDA0002940699430000198
Figure BDA0002940699430000199
Figure BDA00029406994300001910
其中uφ,uθ,uψ是四旋翼姿态控制器,
Figure BDA00029406994300001911
是固定时间自适应律,
Figure BDA00029406994300001912
是神经网络输出。姿态子系统控制器以及网络参数
k1=k3=3000,k2=k4=0.1
k5=k6=1,α3=1,β3=1,
Figure BDA0002940699430000201
α1=α2=50,β1=β2=3,
Figure BDA0002940699430000202
b=20,
Figure BDA0002940699430000203
σz=σx=1
根据引理1可以得到固定时间Tφmax=Tθmax=5.742,Tψmax=4.3784。
仿真结果图9-15验证了算法在四旋翼跟踪控制中的有效性以及优越的收敛速度。图16-17显示了自适应神经网络近似效果。
本实施例提出了一种基于非奇异快速终端滑模面的自适应RBF神经网络双功率滑模固定时间控制。首先,基于非奇异快速终端滑模(NFTSM)面设计了双幂次固定时间控制,以实现高鲁棒性和快速滑动率,并且在控制律中没有负指数项以有效避免奇异现象。然后针对神经网络提出了一种固定时间自适应律,针对非线性系统提出了一种神经网络自适应固定时间控制方案,其收敛时间仅取决于控制参数,而不考虑初始条件。本文提出的算法可以大大提高系统的鲁棒性能,具有较高的实用价值。

Claims (3)

1.一种非线性系统的固定时间自适应神经网络滑模控制方法,其特征在于,包括以下步骤:
步骤S1,根据四轴飞行器的动态运动规律,建立四轴飞行器的动力学模型;
步骤S2,基于非奇异快速终端滑模面提出双幂次固定时间控制律,以实现高鲁棒性和快速滑动率;
步骤S3,仿真验证,以验证针对不确定非线性系统的固定时间自适应神经网络自适应控制的有效性;
步骤S4,数值算例验证,使用Matlab仿真来验证固定时间自适应神经网络控制的有效性;
步骤S1中所述动力学模型选择身体坐标系和地面坐标系,考虑四旋翼无人机系统,选择身体坐标系以及地面坐标系,根据坐标系和系统动力学模型,利用牛顿欧拉方程,得到四旋翼系统的数学模型,姿态角表示为偏航角ψ,俯仰角θ,侧倾角φ;假设四旋翼为刚体,它的动态模型可以表示为:
可以得到四轴飞行器位置和姿态动力学模型:
Figure FDA0003653748760000011
其中,字母上方带点表示所代表含义的导数,字母上方两个点表示所代表含义的二阶导数,m是四旋翼的质量,I=diag(Ix,Iy,Iz)是身体坐标系下三个坐标轴的惯性矩,g是所选取得重力加速度,wi,i=1,2,3,4是身体坐标系下的旋转角速度,ξ=diag(ξxyzφθψ)是空气阻力系数,Ir是转子惯性,
Figure FDA0003653748760000012
是总的剩余转子角度,da(·)=diag[dx dy dz]和dp(·)=diag[dφ dθ dψ]是位置和姿态系统中不确定干扰,Fat表示三个方向上的控制推力的合力,τ是转子产生的扭矩;
引理1假设V(·):Rn→R+∪{0}是一个连续的根本无界函数,并且满足以下两个条件:
Figure FDA0003653748760000013
Figure FDA0003653748760000021
其中a,b,p,q,c是所选定的Lyapunov函数V(x)的系数,且为正实数p∈(0,1),q∈(1,∞),c∈(0,∞);然后系统
Figure FDA0003653748760000022
x(0)=x0的初始x=0则系统是几乎固定时间稳定,此外,具有以下不等式成立
V(x,t)≤ξ,t≥Tmax (4)
其中ξ是方程根;
Figure FDA0003653748760000023
Tmax是固定时间收敛的最大时间;
引理2为杨氏不等式,对于任意常数
Figure FDA0003653748760000024
可以得到不等式:
Figure FDA0003653748760000025
其中pa>1,qa>1且
Figure FDA0003653748760000026
引理3对于正常数
Figure FDA0003653748760000027
向量W*,
Figure FDA0003653748760000028
且满足
Figure FDA0003653748760000029
可以得到不等式
Figure FDA00036537487600000210
Figure FDA00036537487600000211
是理想权重,
Figure FDA00036537487600000212
是估计权重,
Figure FDA00036537487600000213
是近似误差,
Figure FDA00036537487600000214
步骤S2中,所述双幂次固定时间控制律的获取过程如下:
考虑到以下动态系统:
Figure FDA00036537487600000215
dt为不确定干扰,u是系统控制输入;
为使x能够跟踪xd定义跟踪误差:
e1=x1-xd (9)
Figure FDA00036537487600000216
选取非奇异快速终端滑模面:
Figure FDA0003653748760000031
其中α>0,β>0,1<p2<2,p1>p2,
Figure FDA0003653748760000032
g,h,p,q∈N均为奇数且为常数;
Figure FDA0003653748760000033
为使达到估计滑模面,在不考虑干扰条件下,设计等效控制律:
Figure FDA0003653748760000034
采用双幂指数切换控制率:
Figure FDA0003653748760000035
其中k1>0,k2>0且为常数;
最终双幂次固定时间控制律为:
u=ueq+usw (15)
Figure FDA0003653748760000036
然而,二阶非线性系统数学模型中存在未知的f=cos(x1)+dt,控制律(16)并不能完全的达到优越的控制效果;
为了逼近非线性系统中的不确定非线性部分,使用RBF设计神经网络控制:
FNN(x,W)=Wh(x) (17)
Figure FDA0003653748760000037
其中x-网络输入,
Figure FDA0003653748760000038
-权重,Ψ(x)=diag[ψ1(x)…ψl(x)]-节点向量,hi(x)-高斯函数,μi=diag[μi1…μin]-基函数中心,ηi-高斯函数的标量宽度;
RBF神经网络可用于在一个紧急集合
Figure FDA0003653748760000039
上近似任何连续函数
F(x)=W*h(x)+ε (19)
其中ε-神经网咯的近似误差,W*-理想权重
Figure FDA00036537487600000310
Figure FDA00036537487600000311
是估计权重和近似误差
Figure FDA00036537487600000312
选择神经网络控制器
Figure FDA0003653748760000041
其中
Figure FDA0003653748760000042
-网络输出;
非线性系统滑模面可以写成
Figure FDA0003653748760000043
其中
f*=f+dt=W*Ψ(x)+ε (23)
计算系统模型的的近似误差
Figure FDA0003653748760000044
其中W*是理想权重,W*T是理想权重的转置,
Figure FDA0003653748760000045
是估计权重,
Figure FDA0003653748760000046
是估计权重转置,
Figure FDA0003653748760000047
是近似误差;
根据控制器(21)和滑模面(22)可以得到
Figure FDA0003653748760000048
其中k1>0,k2>0,且为常数;
为验证滑模面的稳定,选择Lyapunov候选函数
Figure FDA0003653748760000049
Figure FDA00036537487600000410
为验证网络算法的稳定性,选择Lyapunov候选函数
Figure FDA00036537487600000411
Figure FDA00036537487600000412
Figure FDA00036537487600000413
其中ξ,σzx为大于0的常数;
为验证整体系统的稳定,选择Lyapunov函数
V=Vs+VNN (31)
Figure FDA00036537487600000414
Figure FDA0003653748760000051
根据引理2和引理3可有以下不等式成立
Figure FDA0003653748760000052
Figure FDA0003653748760000053
Figure FDA0003653748760000054
其中σ1>0,σ2>0,σ3>0,σ4>0且为常数,有
Figure FDA0003653748760000055
Figure FDA0003653748760000056
然后可以得到
Figure FDA0003653748760000057
其中
Figure FDA0003653748760000061
然后,根据引理1系统近似固定时间稳定,收敛时间为
Figure FDA0003653748760000062
其中σxz取决于ν,γ,
Figure FDA0003653748760000063
Γ<0;
步骤S3中,所述仿真验证过程如下:
参考轨迹为
xd=sin(t) (41)
然后可以给出控制器和神经网络自适应律为
Figure FDA0003653748760000064
Figure FDA0003653748760000065
其中u是非线性系统控制器,
Figure FDA0003653748760000066
是固定时间自适应律,
Figure FDA0003653748760000067
是神经网络输出;
其中控制器以及网络参数选择为
Figure FDA0003653748760000068
根据引理1可以得到最大收敛时间Txmax=5.742;
系统的初始化条件以及假设存在的干扰dt
Figure FDA0003653748760000069
2.根据权利要求1所述的非线性系统的固定时间自适应神经网络滑模控制方法,其特征在于,为了验证算法应用于实际系统模型的有效性,通过数值算例针对四旋翼飞行器跟踪控制进行有效性以及抗干扰性能的验证,过程如下:
四旋翼飞行器位置动力学模型:
Figure FDA0003653748760000071
其中da(x,t)是四旋翼位置不确定干扰;
其中x代表四旋翼的三个位置状态x,y,z
Figure FDA0003653748760000072
其中A是四旋翼质量对角矩阵,f1(·)是四旋翼位置函数矩阵,us(t)是四旋翼位置虚拟控制输入;
根据四旋翼飞行器位置动力学模型,滑模函数(11),双幂次趋近率(14)以及神经网络控制律(30),可以得到以下虚拟控制器以及自适应律:
Figure FDA0003653748760000073
Figure FDA0003653748760000074
四旋翼飞行器姿态动力学模型:
Figure FDA0003653748760000075
其中db(p,t)是四旋翼姿态不确定干扰;
其中p代表四旋翼的三个姿态φ,θ,ψ
Figure FDA0003653748760000076
其中B是四旋翼姿态惯性矩矩阵,f2(·)是四旋翼姿态函数矩阵,ur(t)是四旋翼控制输入矩阵;根据四旋翼飞行器姿态动力学模型,利用位置子系统控制器设计同样的策略,可以得到姿态控制器为:
Figure FDA0003653748760000077
Figure FDA0003653748760000081
其中
Figure FDA0003653748760000082
是系统误差函数,sx,sp是四旋翼控制滑模面,
Figure FDA0003653748760000083
Figure FDA0003653748760000084
是神经网络输出,可证得姿态子系统固定时间稳定。
3.根据权利要求1所述的非线性系统的固定时间自适应神经网络滑模控制方法,其特征在于,步骤S4中,所述数值算例验证过程如下:
考虑无人机动态系统为公式(10)的四轴飞行器位置和姿态动力学模型,其参数为
m=2,l=0.2,g=9.8
ξx=ξx=ξx=1.2
ξφ=ξθ=ξx=1.2
Ix=1.25,Iy=1.25,Iz=2.5
参考轨迹的选择如下:
Figure FDA0003653748760000085
选择系统的初始条件为
x(0)=y(0)=z(0)=φ(0)=θ(0)=ψ(0)=0.5
然后可以给出虚拟控制器和神经网络自适应律为
Figure FDA0003653748760000086
Figure FDA0003653748760000087
其中ux是四旋翼位置虚拟控制器,
Figure FDA0003653748760000088
是固定时间自适应律,
Figure FDA0003653748760000089
是神经网络输出;
其中位置子系统控制器以及网络参数为
k1=3000,k2=0.1
α=50,β=3,
Figure FDA00036537487600000810
Figure FDA00036537487600000811
σz=σx=1
根据引理1可以得到最大收敛时间Txmax=Tymax=Tzmax=5.742;
采用姿态解算,得到四旋翼飞行器的目标姿态角,假设期望
Figure FDA0003653748760000091
经过位置子系统的姿态计算,得到φd,θd,从而实现跟踪控制;
Qx=u1(sinψsinφ+cosψsinθcosφ)
Qy=u1(-cosψsinφ+sinψsinθcosφ)
Qz=u1(cosθcosφ)
Figure FDA0003653748760000092
Figure FDA0003653748760000093
其中
Figure FDA0003653748760000094
是四旋翼位置控制输入,Qx,Qy,Qz是虚拟控制输入;
然后可以给出控制器和神经网络自适应律为
Figure FDA0003653748760000095
Figure FDA0003653748760000096
Figure FDA0003653748760000097
Figure FDA0003653748760000098
Figure FDA0003653748760000099
Figure FDA00036537487600000910
其中uφ,uθ,uψ是四旋翼姿态控制器,
Figure FDA00036537487600000911
是四旋翼固定时间自适应律,
Figure FDA00036537487600000912
是神经网络输出;
其中姿态子系统控制器以及网络参数
k1=k3=3000,k2=k4=0.1
k5=k6=1,α3=1,β3=1,
Figure FDA00036537487600000913
α1=α2=50,β1=β2=3,
Figure FDA00036537487600000914
b=20,
Figure FDA00036537487600000915
σz=σx=1
根据引理1可以得到最大收敛时间Tφmax=Tθmax=5.742,Tψmax=4.3784。
CN202110175892.9A 2021-02-09 2021-02-09 非线性系统的固定时间自适应神经网络滑模控制方法 Active CN112987567B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110175892.9A CN112987567B (zh) 2021-02-09 2021-02-09 非线性系统的固定时间自适应神经网络滑模控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110175892.9A CN112987567B (zh) 2021-02-09 2021-02-09 非线性系统的固定时间自适应神经网络滑模控制方法

Publications (2)

Publication Number Publication Date
CN112987567A CN112987567A (zh) 2021-06-18
CN112987567B true CN112987567B (zh) 2022-07-22

Family

ID=76392610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110175892.9A Active CN112987567B (zh) 2021-02-09 2021-02-09 非线性系统的固定时间自适应神经网络滑模控制方法

Country Status (1)

Country Link
CN (1) CN112987567B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113341732B (zh) * 2021-06-29 2023-11-21 南京工业大学 一种用于癫痫脑刺激的滑模控制方法
CN113867374B (zh) * 2021-07-23 2023-12-01 福州大学 基于滑模控制的四旋翼无人机参数预测和扰动的自适应轨迹跟踪控制器及其设计方法
CN114355762B (zh) * 2021-12-30 2023-09-26 上海电机学院 一种基于非奇异快速终端滑模的变桨控制方法
CN114371623B (zh) * 2022-01-10 2022-12-23 南京工业大学 一种结构热试验神经网络观测器分数阶全局滑模控制方法
CN114384808B (zh) * 2022-03-24 2022-07-15 西北工业大学 一种基于自适应神经网络的到达时间控制三维制导方法
CN115857342A (zh) * 2022-11-28 2023-03-28 北京理工大学 基于自适应神经网络的航天器观测载荷随动机构控制方法
CN117891281A (zh) * 2024-03-15 2024-04-16 天目山实验室 一种基于嵌套双模态的二阶非奇异终端滑模飞行控制方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104111618A (zh) * 2014-05-28 2014-10-22 河北科技大学 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法
CN107045557A (zh) * 2016-11-01 2017-08-15 长春工业大学 面向约束的可重构机械臂非奇异终端滑模力位置控制方法
CN108459611A (zh) * 2018-05-21 2018-08-28 济南大学 一种近空间飞行器的姿态跟踪控制方法
CN109240321A (zh) * 2018-09-28 2019-01-18 浙江工业大学 一种刚性飞行器的非奇异固定时间神经网络控制方法
CN109426146A (zh) * 2017-08-23 2019-03-05 中国空气动力研究与发展中心计算空气动力研究所 高超声速飞行器的高阶非奇异Terminal滑模控制方法
CN110109352A (zh) * 2019-04-11 2019-08-09 浙江工业大学 一种三自由度四旋翼飞行器的固定时间自适应姿态控制方法
CN110154028A (zh) * 2019-05-29 2019-08-23 天津大学 机械臂无模型自适应积分终端滑模控制方法
CN110221542A (zh) * 2019-06-04 2019-09-10 西北工业大学 一种二阶非线性多智能体系统固定时间协同跟踪控制方法
GB202019112D0 (en) * 2020-10-26 2021-01-20 Univ Guizhou Fractional-order MEMS gyroscope acceleration adaptive backstepping control method without accurate reference trajectory

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10184797B2 (en) * 2015-12-18 2019-01-22 Invensense, Inc. Apparatus and methods for ultrasonic sensor navigation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104111618A (zh) * 2014-05-28 2014-10-22 河北科技大学 一种活性污泥污水处理曝气环节的建模、滑模控制器及其控制方法
CN107045557A (zh) * 2016-11-01 2017-08-15 长春工业大学 面向约束的可重构机械臂非奇异终端滑模力位置控制方法
CN109426146A (zh) * 2017-08-23 2019-03-05 中国空气动力研究与发展中心计算空气动力研究所 高超声速飞行器的高阶非奇异Terminal滑模控制方法
CN108459611A (zh) * 2018-05-21 2018-08-28 济南大学 一种近空间飞行器的姿态跟踪控制方法
CN109240321A (zh) * 2018-09-28 2019-01-18 浙江工业大学 一种刚性飞行器的非奇异固定时间神经网络控制方法
CN110109352A (zh) * 2019-04-11 2019-08-09 浙江工业大学 一种三自由度四旋翼飞行器的固定时间自适应姿态控制方法
CN110154028A (zh) * 2019-05-29 2019-08-23 天津大学 机械臂无模型自适应积分终端滑模控制方法
CN110221542A (zh) * 2019-06-04 2019-09-10 西北工业大学 一种二阶非线性多智能体系统固定时间协同跟踪控制方法
GB202019112D0 (en) * 2020-10-26 2021-01-20 Univ Guizhou Fractional-order MEMS gyroscope acceleration adaptive backstepping control method without accurate reference trajectory

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
"Nonsingular Fast Terminal Sliding Mode Tracking Control for a Class of Uncertain Nonlinear Systems";Siyi Chen 等;《Journal of Control Science and Engineering》;20190530;第1-18页 *
"三自由度四旋翼飞行器抗干扰控制研究";陈荣;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20200831;第13-37页 *
《Adaptive Neural Network Sliding Mode Control for Quad Tilt Rotor Aircraft》;Yanchao Yin 等;《Hindawi》;20171011;第1-14页 *
《Adaptive_Neural_Network_Nonsingular_Fast_Terminal_Sliding_Mode_Control_for_Permanent_Magnet_Linear_Synchronous_Motor》;XIMEI ZHAO 等;《IEEE》;20191210;第180361-180372页 *
《RBFNN-based nonsingular fast terminal sliding mode control for robotic manipulators including actuator dynamics》;Ziyang Chen 等;《Neurocomputing》;20190719;第72-82页 *
《二阶非线性不确定系统分段非奇异快速终端滑模控制》;张贝贝 等;《高技术通讯》;20171231;第938-947页 *
《基于组合滑模的升力受限四旋翼轨迹跟踪》;陶少俊 等;《控制工程》;20201130;第1907-1914页 *
《基于高阶微分滑模面的不确定飞机控制系统自适应滑模控制》;武晓晶 等;《河北科技大学学报》;20181031;第430-437页 *

Also Published As

Publication number Publication date
CN112987567A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112987567B (zh) 非线性系统的固定时间自适应神经网络滑模控制方法
Lian et al. Adaptive attitude control of a quadrotor using fast nonsingular terminal sliding mode
Miranda-Colorado et al. Reduction of power consumption on quadrotor vehicles via trajectory design and a controller-gains tuning stage
Shin Adaptive dynamic surface control for a hypersonic aircraft using neural networks
CN110347170B (zh) 可重复使用运载器再入段鲁棒容错制导控制系统及工作方法
CN112965371B (zh) 基于固定时间观测器的水面无人艇轨迹快速跟踪控制方法
Sierra et al. Wind and payload disturbance rejection control based on adaptive neural estimators: application on quadrotors
CN113031446A (zh) 不确定时滞非线性系统非奇异性神经自适应跟踪控制方法
Zeghlache et al. Fault tolerant control based on neural network interval type-2 fuzzy sliding mode controller for octorotor UAV
CN115576341A (zh) 基于函数微分和自适应变增益的无人机轨迹跟踪控制方法
Sun et al. An integrated backstepping and sliding mode tracking control algorithm for unmanned underwater vehicles
Suresh et al. An on-line learning neural controller for helicopters performing highly nonlinear maneuvers
Wang et al. Neural network adaptive inverse model control method for quadrotor UAV
CN113377121A (zh) 一种基于深度强化学习的飞行器智能抗扰动控制方法
Qin et al. Dual-loop robust attitude control for an aerodynamic system with unknown dynamic model: Algorithm and experimental validation
Akbar et al. Adaptive modified super-twisting control for a quadrotor helicopter with a nonlinear sliding surface
CN109108964B (zh) 一种基于自适应动态规划Nash博弈的空间机械臂协调控制方法
Aruneshwaran et al. Neural adaptive flight controller for ducted-fan UAV performing nonlinear maneuver
Zhang et al. A safety planning and control architecture applied to a quadrotor autopilot
Boudjedir et al. Dual neural network for adaptive sliding mode control of quadrotor helicopter stabilization
Song et al. Enhanced fireworks algorithm-auto disturbance rejection control algorithm for robot fish path tracking
Han et al. Fuzzy modeling and control for a nonlinear quadrotor under network environment
Choi et al. A simple fuzzy-approximation-based adaptive control of uncertain unmanned helicopters
An et al. Prediction of Angular Acceleration for Practical Realization of Incremental Nonlinear Trajectory-tracking Control for Aircrafts
Ma et al. Finite-time trajectory tracking control of quadrotor UAV via adaptive RBF neural network with lumped uncertainties

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