CN111650954B - 基于深度学习的四旋翼无人机地面效应补偿降落控制方法 - Google Patents

基于深度学习的四旋翼无人机地面效应补偿降落控制方法 Download PDF

Info

Publication number
CN111650954B
CN111650954B CN202010546066.6A CN202010546066A CN111650954B CN 111650954 B CN111650954 B CN 111650954B CN 202010546066 A CN202010546066 A CN 202010546066A CN 111650954 B CN111650954 B CN 111650954B
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
ground effect
controller
neural network
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
CN202010546066.6A
Other languages
English (en)
Other versions
CN111650954A (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN202010546066.6A priority Critical patent/CN111650954B/zh
Publication of CN111650954A publication Critical patent/CN111650954A/zh
Application granted granted Critical
Publication of CN111650954B publication Critical patent/CN111650954B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及四旋翼无人机降落控制技术,为实现四旋翼无人机降落过程中能够较好地抑制地面效应对无人机的干扰,同时对深度学习网络进行谱归一化,从而保证本发明提出的非线性控制器的稳定性。为此,本发明采用的技术方案是,基于深度学习的四旋翼无人机地面效应补偿降落控制方法,包括如下步骤:建立四旋翼无人机降落过程的非线性动力学模型,利用谱归一化的深度神经网络拟合地面效应,进而设计非线性控制器实现四旋翼无人机的降落控制。本发明主要应用于无人机自动控制降落场合。

Description

基于深度学习的四旋翼无人机地面效应补偿降落控制方法
技术领域
本发明涉及四旋翼无人机降落控制技术,具体讲,涉及针对四旋翼无人机在地面效应影响下的降落控制方法。
背景技术
四旋翼无人机具有能垂直起降、结构简单与灵活度高的特点,可以实现垂直着陆回收。因此,四旋翼无人机的精准自主降落回收受到了很多研究人员的关注。四旋翼无人机在降落过程中,旋翼逐渐靠近地面,气流经地面反弹会对无人机产生较强的扰动。地面效应对无人机产生的干扰会降低无人机降落的精度。
目前国内外研究人员的相关研究工作主要集中在无人机相对于降落位置的定位问题,对于四旋翼无人机的地面效应抑制降落控制问题研究较少。
对于四旋翼无人机地面效应的抑制控制问题,目前主要使用两种策略:第一种为使用基于Cheeseman-Bennett模型改进的地面效应模型,结合不同的控制算法进行地面效应抑制;第二种为使用鲁棒控制、滑模控制等方法将地面效应作为扰动进行抑制。
对于以上两种控制策略,国内外的研究人员均进行了研究。Li Danjun(会议:201534th Chinese Control Conference(CCC);著者:Li D.J.,Zhou Y.,Shi Z.Y.等;出版年月:2015年;文章题目:Autonomous landing of quadrotor based on ground effectmodelling;页码:5647-5652)等人对Cheeseman-Bennett模型进行了改进使其更适合于四旋翼无人机,并设计了高度方向的鲁棒控制器。Keshavarzian H等人(期刊:AerospaceScience and Technology;著者:Keshavarzian.H,Daneshjou.K;出版年月:2019年;文章题目:Modified under-actuated quadrotor model for forwarding flight in thepresence of ground effect;页码:242-252)在改进Cheeseman_Bennett模型的基础上基于反步法设计了控制器,该研究使用预先采集的飞行数据计算改进模型的参数。Lee D等人(会议:2012IEEE international conference on robotics and automation;著者:LeeD.,Ryan T.,Kim H.J.;出版年月:2012年;文章题目:Autonomous landing of a VTOL UAVon a moving platform using image-based visual servoing;页码:971-976)设计了自适应滑模控制器以抑制地面效应。
但是,目前的控制策略均有其自身的局限性。被广泛使用的Cheeseman-Bennett模型主要针对旋翼以恒定速率旋转的大型直升机。而四旋翼无人机由于旋翼直径小且转速变化较大,所以难以使用该模型进行描述。此外,由于四旋翼无人机各旋翼之间的相互作用,更增加了四旋翼无人机地面效应的复杂程度(会议:2018IEEE International Conferenceon Robotics and Automation(ICRA);著者:Conyers S.A.,Rutherford M.J.,ValavanisK.P.;出版年月:2018年;文章题目:An empirical evaluation of ground effect forsmall-scale rotorcraft;页码:1244-1250)。虽然很多研究都对Cheeseman-Bennett模型进行了改进使其适用于四旋翼无人机,但实际效果一般,且其依然需要基于气流场平稳的假设。直接设计滑模控制器抑制地面效应的方法响应较慢,且容易引起抖震。此外,基于飞行数据拟合地面效应模型的方法容易对控制器的稳定性产生影响。
发明内容
为克服现有技术的不足,本发明旨在提出一种基于深度学习的非线性控制器,实现四旋翼无人机降落过程中能够较好地抑制地面效应对无人机的干扰,同时对深度学习网络进行谱归一化,从而保证本发明提出的非线性控制器的稳定性。为此,本发明采用的技术方案是,基于深度学习的四旋翼无人机地面效应补偿降落控制方法,包括如下步骤:建立四旋翼无人机降落过程的非线性动力学模型,利用谱归一化的深度神经网络拟合地面效应,进而设计非线性控制器实现四旋翼无人机的降落控制。
具体步骤如下:
首先对四旋翼无人机降落过程进行受力分析,进而获得四旋翼无人机降落过程的非线性动力学模型:
Figure BDA0002540797940000021
在式(1)中各变量定义如下:P,v为无人机在世界坐标系下的坐标和速度,g=[0,0,-g]T为重力加速度向量,ω为无人机的角速度,R为机体坐标系到世界坐标系的变换矩阵,且R∈SO(3),S(·)为一个斜对称矩阵,
Figure BDA0002540797940000022
分别表示P,v,R,ω的微分,fu为升力向量,fu=[0,0,T]T,且T为机体坐标系下的升力,fg为无人机受到的地面效应,τu为转矩向量,τu=[τx,τy,τz]T,τx,τy,τz分别为无人机三个姿态角方向的力矩,采用η=[T,τx,τy,τz]T表示无人机自身的动力,选取四个电机的转速的平方向量
Figure BDA0002540797940000023
Figure BDA0002540797940000024
作为控制输入,则控制输入u与无人机动力η的关系为η=Bu,输入系数矩阵B的定义如下:
Figure BDA0002540797940000025
其中,cT,cQ为升力系数与转矩系数,larm为无人机臂长;
其次,使用ReLU深度神经网络对地面效应进行拟合,ReLU深度神经网络用下列等式表示:
f(x,θ)=WL+1φ(WLφ(WL-1(…φ(W1x)…))) (3)
式(3)中,θ=[W1,W2,…,WL+1]为神经网络的权值,φ(x)=max(0,x)为ReLU激励函数;
在神经网络训练过程中需对神经网络进行谱归一化,即对网络权值做如下式的处理:
Figure BDA0002540797940000026
其中γ为一个正常数,对于式(3)定义的ReLU深度神经网络f(x,θ),利用式(4)进行谱归一化,则神经网络的利普希茨常数满足下列不等式:
Figure BDA0002540797940000031
式(5)中,
Figure BDA0002540797940000032
为网络谱归一化后的参数,
Figure BDA0002540797940000033
γ是一个正常数;
使用ReLU深度神经网络对四旋翼无人机地面效应的拟合,采用随机梯度下降法优化网络参数,并在训练过程中对网络进行谱归一化的过程表示为:
Figure BDA0002540797940000034
其中,yt为第t个地面效应数据,xt为无人机状态量与控制输入;
然后进行控制器设计:定义位置跟踪误差e=[ex,ey,ez]T如下:
e=Pd-P (7)
这里Pd=[pdx,pdy,pdz]T代表无人机降落过程中的期望轨迹pdx,pdy,pdz为Pd在坐标系x,y,z轴上的投影。定义滑模面s如下:
Figure BDA0002540797940000035
其中,
Figure BDA0002540797940000036
为的e微分
Figure BDA0002540797940000037
为对称增益矩阵,角标x,y,z分表表示三维坐标系下三个方向上的投影。定义无人机总推力fd=Rfu,通过设计fd进行控制器设计。
无人机的位置控制器设计为:
Figure BDA0002540797940000038
式(9)中,
Figure BDA0002540797940000039
Figure BDA00025407979400000310
式(9)中
Figure BDA00025407979400000322
分别表示Pd的二阶微分与一阶微分,K1=diag{k1x,k1y,k1z}为增益矩阵,
Figure BDA00025407979400000313
为深度神经网络对地面效应的估计值,ξg为无人机状态量,由无人机的高度,姿态等与地面效应相关的状态组成,u为控制量,式(11)中
Figure BDA00025407979400000314
表示Y的微分,K2=diag{k1x,k2y,k2z}为增益矩阵,角标x,y,z均表示在三维坐标系下三个坐标方向的投影。定义辅助函数向量
Figure BDA00025407979400000315
如下式所示:
Figure BDA00025407979400000316
并且有
Figure BDA00025407979400000317
定义
Figure BDA00025407979400000318
为机体坐标系z轴(朝上)上的单位向量在惯性坐标系下的表示,则控制器输出的总推力期望为
Figure BDA00025407979400000319
进一步地,根据fd的方向,结合偏航角的目标值解算姿态角的目标值Rd,无人机内环控制器将根据姿态角的目标值Rd计算无人机姿态角的控制转矩τd,进而解得控制量
Figure BDA00025407979400000320
即:
Figure BDA00025407979400000321
式(13)所示系统为一个非仿射系统,u无法直接解出,由于实际的控制系统均为离散系统,采用下式求解控制量:
Figure BDA0002540797940000041
其中,uk-1为上一时刻的控制量;
若控制器中的控制增益k1,k2满足如下不等式时,误差e在有限时间内收敛到零;
Figure BDA0002540797940000042
其中μ=-k1p12
Figure BDA0002540797940000043
Figure BDA0002540797940000044
p11=1,
Figure BDA0002540797940000048
a,Γ均为正常数且β>1。
还包括验证步骤:根据谱归一化后的神经网络利普西斯常数小于1,结合Lyapunov分析方法可以证明误差e在有限时间内收敛到零。
还包括进行地面效应数据收集以及实验验证步骤,地面效应数据收集方法如下:
地面效应的真值通过机载传感器获得,结合飞行控制器上的惯性测量单元IMU、飞行控制器的输出信号、以及动作捕捉系统的信息,解算无人机在整个飞行期间受到的地面效应,地面效应的数值由如下所示的公式进行解算:
fg=ma+mg-Rfu (16)
其中,m为无人机质量,a为无人机的加速度,fu为总升力,且有
Figure BDA0002540797940000045
CT为升力系数,
Figure BDA0002540797940000046
为四个桨的转速平方之和,由于实际飞控为PWM输出,PWM波占空比与转速平方成正比,即
Figure BDA0002540797940000047
因此fu=[0,0,CTρpwm·PWMtotal]T,由于无人机悬停时产生的升力等于自身重力,基于此对CTρpwm进行标定,式(16)中a可由无人机机载加速度计得到,R可由机载陀螺仪数据得到,fu可由PWM输出结合标定结果得出,因此可以由式(16)计算出无人机整个飞行期间的地面效应;
在获得地面效应真值后,使用神经网络进行拟合,结合设计的非线性控制器在无人机平台上进行实验,验证控制器的合理性和可行性。
本发明的特点及有益效果是:
本发明设计基于深度学习的四旋翼无人机降落控制器,对于抑制四旋翼无人机在降落过程中的地面效应干扰具有较好的效果。同时,对深度神经网络进行了谱归一化,保证了控制器的稳定性,使四旋翼无人机在有地面干扰的情况下高度跟踪误差在有限时间内收敛到零。
附图说明:
图1为本发明的无人机实验平台。
图2为实验过程中采集的地面效应数据,用于神经网络的训练。
图3为无补偿的无人机高度变化曲线。
图4为使用本发明控制器的无人机高度变化曲线。
图5为本发明的工作流程。
具体实施方式
为克服现有技术的不足,本发明旨在提出一种基于深度学习的非线性控制器,实现四旋翼无人机降落过程中能够较好地抑制地面效应对无人机的干扰,同时对深度学习网络进行谱归一化,从而保证本发明提出的非线性控制器的稳定性。本发明采用的技术方案是,建立非线性动力学模型,利用谱归一化的深度神经网络拟合地面效应,进而设计非线性控制器实现四旋翼无人机的降落控制。
进一步具体步骤是,首先对四旋翼无人机降落过程进行受力分析,进而获得四旋翼无人机降落过程的非线性动力学模型:
Figure BDA0002540797940000051
在式(1)中各变量定义如下:P,v为无人机在世界坐标系下的坐标和速度,g=[0,0,-g]T为重力加速度向量,ω为无人机的角速度,R为机体坐标系到世界坐标系的变换矩阵,且R∈SO(3),S(·)为一个斜对称矩阵,
Figure BDA0002540797940000052
分别表示P,v,R,ω的微分,fu为升力向量,fu=[0,0,T]T,且T为机体坐标系下的升力,fg为无人机受到的地面效应,τu为转矩向量,τu=[τx,τy,τz]T,τx,τy,τz分别为无人机三个姿态角方向的力矩,采用η=[T,τx,τy,τz]T表示无人机自身的动力,选取四个电机的转速的平方向量
Figure BDA0002540797940000053
Figure BDA0002540797940000054
作为控制输入,则控制输入u与无人机动力η的关系为η=Bu,输入系数矩阵B的定义如下:
Figure BDA0002540797940000055
其中,cT,cQ为升力系数与转矩系数,larm为无人机臂长;
其次,使用ReLU深度神经网络对地面效应进行拟合,ReLU深度神经网络用下列等式表示:
f(x,θ)=WL+1φ(WLφ(WL-1(…φ(W1x)…))) (3)
式(3)中,θ=[W1,W2,…,WL+1]为神经网络的权值,φ(x)=max(0,x)为ReLU激励函数;
在神经网络训练过程中需对神经网络进行谱归一化,即对网络权值做如下式的处理:
Figure BDA0002540797940000056
其中γ为一个正常数,对于式(3)定义的ReLU深度神经网络f(x,θ),利用式(4)进行谱归一化,则神经网络的利普希茨常数满足下列不等式:
Figure BDA0002540797940000057
式(5)中,
Figure BDA0002540797940000058
为网络谱归一化后的参数,
Figure BDA0002540797940000059
γ是一个正常数;
使用ReLU深度神经网络对四旋翼无人机地面效应的拟合,采用随机梯度下降法优化网络参数,并在训练过程中对网络进行谱归一化的过程表示为:
Figure BDA0002540797940000061
其中,yt为第t个地面效应数据,xt为无人机状态量与控制输入;
然后进行控制器设计:定义位置跟踪误差e=[ex,ey,ez]T如下:
e=Pd-P (7)
这里Pd=[pdx,pdy,pdz]T代表无人机降落过程中的期望轨迹pdx,pdy,pdz为Pd在坐标系x,y,z轴上的投影。定义滑模面s如下:
Figure BDA0002540797940000062
其中,
Figure BDA0002540797940000063
为的e微分
Figure BDA0002540797940000064
为对称增益矩阵,角标x,y,z分表表示三维坐标系下三个方向上的投影。定义无人机总推力fd=Rfu,通过设计fd进行控制器设计。
无人机的位置控制器设计为:
Figure BDA0002540797940000065
式(9)中,
Figure BDA0002540797940000066
Figure BDA0002540797940000067
式(9)中
Figure BDA00025407979400000620
分别表示Pd的二阶微分与一阶微分,K1=diag{k1x,k1y,k1z}为增益矩阵,
Figure BDA00025407979400000610
为深度神经网络对地面效应的估计值,ξg为无人机状态量,由无人机的高度,姿态等与地面效应相关的状态组成,u为控制量,式(11)中
Figure BDA00025407979400000619
表示Y的微分,K2=diag{k1x,k2y,k2z}为增益矩阵,角标x,y,z均表示在三维坐标系下三个坐标方向的投影。定义辅助函数向量
Figure BDA00025407979400000611
如下式所示:
Figure BDA00025407979400000612
并且有
Figure BDA00025407979400000613
定义
Figure BDA00025407979400000614
为机体坐标系z轴(朝上)上的单位向量在惯性坐标系下的表示,则控制器输出的总推力期望为
Figure BDA00025407979400000615
进一步地,根据fd的方向,结合偏航角的目标值解算姿态角的目标值Rd,无人机内环控制器将根据姿态角的目标值Rd计算无人机姿态角的控制转矩τd,进而解得控制量
Figure BDA00025407979400000616
即:
Figure BDA00025407979400000617
式(13)所示系统为一个非仿射系统,u无法直接解出,由于实际的控制系统均为离散系统,采用下式求解控制量:
Figure BDA00025407979400000618
其中,uk-1为上一时刻的控制量;
若控制器中的控制增益k1,k2满足如下不等式时,误差e在有限时间内收敛到零;
Figure BDA0002540797940000071
其中μ=-k1p12
Figure BDA0002540797940000072
Figure BDA0002540797940000073
p11=1,
Figure BDA0002540797940000074
a,Γ均为正常数且β>1。
验证的具体步骤为,根据谱归一化后的神经网络利普西斯常数小于1,结合Lyapunov分析方法可以证明误差e在有限时间内收敛到零。
还包括进行地面效应数据收集以及实验验证步骤,地面效应数据收集方法如下:
地面效应的真值通过机载传感器获得,结合飞行控制器上的惯性测量单元IMU、飞行控制器的输出信号、以及动作捕捉系统的信息,解算无人机在整个飞行期间受到的地面效应,地面效应的数值由如下所示的公式进行解算:
fg=ma+mg-Rfu (16)
其中,为无人机质量,a为无人机的加速度,fu为总升力,且有
Figure BDA0002540797940000075
CT为升力系数,
Figure BDA0002540797940000076
为四个桨的转速平方之和,由于实际飞控为PWM输出,PWM波占空比与转速平方成正比,即
Figure BDA0002540797940000077
因此fu=[0,0,CTρpwm·PWMtotal]T,由于无人机悬停时产生的升力等于自身重力,基于此对CTρpwm进行标定,式(16)中a可由无人机机载加速度计得到,R可由机载陀螺仪数据得到,fu可由PWM输出结合标定结果得出,因此可以由式(16)计算出无人机整个飞行期间的地面效应;
在获得地面效应真值后,使用神经网络进行拟合,结合设计的非线性控制器在无人机平台上进行实验,验证控制器的合理性和可行性。
本发明所要解决的技术问题是:四旋翼无人机在降落过程中会受到地面效应的干扰,增加了降落控制难度,降低了降落的精度。
本发明采取的技术方案是:建立非线性动力学模型,利用谱归一化的深度神经网络拟合地面效应,进而设计非线性控制器实现四旋翼无人机的降落控制。包括如下步骤:
首先对四旋翼无人机降落过程进行受力分析,进而获得四旋翼无人机降落过程的非线性动力学模型:
Figure BDA0002540797940000078
在式(1)中各变量定义如下:P,v为无人机在世界坐标系下的坐标和速度,g=[0,0,-g]T为重力加速度向量,ω为无人机的角速度,R为机体坐标系到世界坐标系的变换矩阵,且R∈SO(3),S(·)为一个斜对称矩阵,
Figure BDA0002540797940000079
分别表示P,v,R,ω的微分,fu为升力向量,fu=[0,0,T]T,且T为机体坐标系下的升力,fg为无人机受到的地面效应,τu为转矩向量,τu=[τx,τy,τz]T,τx,τy,τz分别为无人机三个姿态角方向的力矩,采用η=[T,τx,τy,τz]T表示无人机自身的动力,选取四个电机的转速的平方向量
Figure BDA0002540797940000081
Figure BDA0002540797940000082
作为控制输入,则控制输入u与无人机动力η的关系为η=Bu,输入系数矩阵B的定义如下:
Figure BDA0002540797940000083
其中,cT,cQ为升力系数与转矩系数,1arm为无人机臂长;
其次,使用ReLU深度神经网络对地面效应进行拟合,ReLU深度神经网络用下列等式表示:
f(x,θ)=WL+1φ(WLφ(WL-1(…φ(W1x)…))) (3)
式(3)中,θ=[W1,W2,…,WL+1]为神经网络的权值,φ(x)=max(0,x)为ReLU激励函数;
在神经网络训练过程中需对神经网络进行谱归一化,即对网络权值做如下式的处理:
Figure BDA0002540797940000084
其中γ为一个正常数,对于式(3)定义的ReLU深度神经网络f(x,θ),利用式(4)进行谱归一化,则神经网络的利普希茨常数满足下列不等式:
Figure BDA0002540797940000085
式(5)中,
Figure BDA0002540797940000086
为网络谱归一化后的参数,
Figure BDA0002540797940000087
γ是一个正常数;
使用ReLU深度神经网络对四旋翼无人机地面效应的拟合,采用随机梯度下降法优化网络参数,并在训练过程中对网络进行谱归一化的过程表示为:
Figure BDA0002540797940000088
其中,yt为第t个地面效应数据,xt为无人机状态量与控制输入;
然后进行控制器设计:定义位置跟踪误差e=[ex,ey,ez]T如下:
e=Pd-P (7)
这里Pd=[pdx,pdy,pdz]T代表无人机降落过程中的期望轨迹pdx,pdy,pdz为Pd在坐标系x,y,z轴上的投影。定义滑模面s如下:
Figure BDA0002540797940000089
其中,
Figure BDA00025407979400000810
为的e微分
Figure BDA00025407979400000811
为对称增益矩阵,角标x,y,z分表表示三维坐标系下三个方向上的投影。定义无人机总推力fd=Rfu,通过设计fd进行控制器设计。
无人机的位置控制器设计为:
Figure BDA00025407979400000812
式(9)中,
Figure BDA0002540797940000091
Figure BDA0002540797940000092
式(9)中
Figure BDA0002540797940000093
分别表示Pd的二阶微分与一阶微分,K1=diag{k1x,k1y,k1z}为增益矩阵,
Figure BDA0002540797940000094
为深度神经网络对地面效应的估计值,ξg为无人机状态量,由无人机的高度,姿态等与地面效应相关的状态组成,u为控制量,式(11)中
Figure BDA0002540797940000095
表示Y的微分,K2=diag{k1x,k2y,k2z}为增益矩阵,角标x,y,z均表示在三维坐标系下三个坐标方向的投影。定义辅助函数向量
Figure BDA0002540797940000096
如下式所示:
Figure BDA0002540797940000097
并且有
Figure BDA0002540797940000098
定义
Figure BDA0002540797940000099
为机体坐标系z轴(朝上)上的单位向量在惯性坐标系下的表示,则控制器输出的总推力期望为
Figure BDA00025407979400000910
进一步地,根据fd的方向,结合偏航角的目标值解算姿态角的目标值Rd,无人机内环控制器将根据姿态角的目标值Rd计算无人机姿态角的控制转矩τd,进而解得控制量
Figure BDA00025407979400000911
即:
Figure BDA00025407979400000912
式(13)所示系统为一个非仿射系统,u无法直接解出,由于实际的控制系统均为离散系统,采用下式求解控制量:
Figure BDA00025407979400000913
其中,uk-1为上一时刻的控制量;
若控制器中的控制增益k1,k2满足如下不等式时,误差e在有限时间内收敛到零;
Figure BDA00025407979400000914
其中μ=-k1p12
Figure BDA00025407979400000915
Figure BDA00025407979400000916
p11=1,
Figure BDA00025407979400000917
a,Γ均为正常数且β>1。
验证的具体步骤为,根据谱归一化后的神经网络利普西斯常数小于1,结合Lyapunov分析方法可以证明误差e在有限时间内收敛到零。
还包括进行地面效应数据收集以及实验验证步骤,地面效应数据收集方法如下:
地面效应的真值通过机载传感器获得,结合飞行控制器上的惯性测量单元IMU、飞行控制器的输出信号、以及动作捕捉系统的信息,解算无人机在整个飞行期间受到的地面效应,地面效应的数值由如下所示的公式进行解算:
fg=ma+mg-Rfu (16)
其中,m为无人机质量,a为无人机的加速度,fu为总升力,且有
Figure BDA00025407979400000918
CT为升力系数,
Figure BDA00025407979400000919
为四个桨的转速平方之和,由于实际飞控为PWM输出,PWM波占空比与转速平方成正比,即
Figure BDA0002540797940000101
因此fu=[0,0,CTρpwm·PWMtotal]T,由于无人机悬停时产生的升力等于自身重力,基于此对CTρpwm进行标定,式(16)中a可由无人机机载加速度计得到,R可由机载陀螺仪数据得到,fu可由PWM输出结合标定结果得出,因此可以由式(16)计算出无人机整个飞行期间的地面效应;
为验证本文提出本发明的基于深度学习的四旋翼无人机地面效应补偿降落控制器的有效性,搭建实验平台进行实验验证。
一、实验平台简介
实验平台由两台计算机、一架四旋翼无人机以及动作捕捉系统组成。其中,捕捉系统实时捕获无人机的位置信息,使用一台计算机作为数据处理计算机,将结果通过WIFI无线网络发送到无人机。另一台计算计算机用于地面效应的拟合训练,并根据无人机的位置及控制量将地面效应的估计值实时地发送给无人机。四旋翼无人机由机架,飞行控制器,电机与螺旋桨构成。无人机质量为0.91kg,机架直径为35cm,螺旋桨直径为24cm,螺距为12.7cm。如图1所示。
二、降落控制实验
首先进行地面效应数据的获取,手动控制无人机飞行获取飞行数据。利用飞行数据计算地面效应真值,如图2所示。
无人机飞行控制器参数选取为Λ=diag{1.25,1.25,2.43},K1=diag{0.83,0.83,0.993},K2=diag{0.03,0.03,0.073},进行从100cm高度降落实验。分别进行无补偿的降落实验以及使用本发明的控制器的降落实验,实验结果分别如图3、图4所示。
图3为无补偿控制中无人机的高度变化曲线。无人机在6.6秒起飞,目标值为100cm;在48.82秒降落,目标值为0cm。无人机在6.6秒起飞后,在控制器的作用下在11.86秒进入稳态范围。无人机在48.82秒开始降落,并且在52.08秒开始进入到无法降落的状态。结果表明由于受到强烈的地面效应作用,无人机降落时无法像起飞时一样快速收敛到目标值。图4为使用本发明控制器的无人机的高度变化曲线。无人机在3.2秒起飞,目标值为100cm;在8.3秒进入稳态区间;在40.2秒开始降落,在45.7秒降落成功。可以看出,无人机从100cm向下降落过程的收敛时间(5.3秒)与起飞至100cm收敛时间(5.1秒)一致,地面效应对无人机降落过程影响较小。
经过上述分析与实验,证明了本发明所提算法的有效性与可行性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于深度学习的四旋翼无人机地面效应补偿降落控制方法,其特征是,包括如下步骤:建立四旋翼无人机降落过程的非线性动力学模型,利用谱归一化的深度神经网络拟合地面效应,进而设计非线性控制器实现四旋翼无人机的降落控制;具体步骤如下:
首先对四旋翼无人机降落过程进行受力分析,进而获得四旋翼无人机降落过程的非线性动力学模型:
Figure FDA0003641256040000011
在式(1)中各变量定义如下:P,v为无人机在世界坐标系下的坐标和速度,g=[0,0,-g]T为重力加速度向量,ω为无人机的角速度,R为机体坐标系到世界坐标系的变换矩阵,且R∈SO(3),S(·)为一个斜对称矩阵,
Figure FDA0003641256040000012
分别表示P,v,R,ω的微分,fu为升力向量,fu=[0,0,T]T,且T为机体坐标系下的升力,fg为无人机受到的地面效应,τu为转矩向量,τu=[τx,τy,τz]T,τx,τy,τz分别为无人机三个姿态角方向的力矩,采用η=[T,τx,τy,τz]T表示无人机自身的动力,选取四个电机的转速的平方向量
Figure FDA0003641256040000013
作为控制输入,则控制输入u与无人机动力η的关系为η=Bu,输入系数矩阵B的定义如下:
Figure FDA0003641256040000014
其中,cT,cQ为升力系数与转矩系数,larm为无人机臂长;
其次,使用ReLU深度神经网络对地面效应进行拟合,ReLU深度神经网络用下列等式表示:
f(x,θ)=WL+1φ(WLφ(WL=1(…φ(W1x)…))) (3)
式(3)中,θ=[W1,W+,…,WL+1]为神经网络的权值,φ(x)=max(0,x)为ReLU激励函数;
在神经网络训练过程中需对神经网络进行谱归一化,即对网络权值做如下式的处理:
Figure FDA0003641256040000015
其中γ为一个正常数,对于式(3)定义的ReLU深度神经网络f(x,θ),利用式(4)进行谱归一化,则神经网络的利普希茨常数满足下列不等式:
Figure FDA0003641256040000016
式(5)中,
Figure FDA0003641256040000017
为网络谱归一化后的参数,
Figure FDA0003641256040000018
γ是一个正常数;
使用ReLU深度神经网络对四旋翼无人机地面效应的拟合,采用随机梯度下降法优化网络参数,并在训练过程中对网络进行谱归一化的过程表示为:
Figure FDA0003641256040000021
其中,yt为第t个地面效应数据,xt为无人机状态量与控制输入;
然后进行控制器设计:定义位置跟踪误差e=[ex,ey,ez]T如下:
e=Pd-P (7)
这里Pd=[pdx,pdy,pdz]T代表无人机降落过程中的期望轨迹pdx,pdy,pdz为Pd在坐标系x,y,z轴上的投影;定义滑模面s如下:
Figure FDA0003641256040000022
其中,
Figure FDA0003641256040000023
为e的微分,
Figure FDA0003641256040000024
为对称增益矩阵,角标x,y,z分表表示三维坐标系下三个方向上的投影;定义无人机总推力fd=Rfu,通过设计fd进行控制器设计;无人机的位置控制器设计为:
Figure FDA0003641256040000025
式(9)中,
Figure FDA0003641256040000026
Figure FDA0003641256040000027
式(9)中
Figure FDA0003641256040000028
分别表示Pd的二阶微分与一阶微分,K1=didg{k1x,k1y,k1z}为增益矩阵,
Figure FDA0003641256040000029
为深度神经网络对地面效应的估计值,ξg为无人机状态量,由无人机的高度,姿态与地面效应相关的状态组成,u为控制量,式(11)中
Figure FDA00036412560400000210
表示Y的微分,K2=didg{k1x,k2y,k2z}为增益矩阵,角标x,y,z均表示在三维坐标系下三个坐标方向的投影;定义辅助函数向量
Figure FDA00036412560400000211
如下式所示:
Figure FDA00036412560400000212
并且有
Figure FDA00036412560400000213
定义
Figure FDA00036412560400000214
为机体坐标系z轴上的单位向量在惯性坐标系下的表示,则控制器输出的总推力期望为
Figure FDA00036412560400000215
根据fd的方向,结合偏航角的目标值解算姿态角的目标值Rd,无人机内环控制器将根据姿态角的目标值Rd计算无人机姿态角的控制转矩τd,进而解得控制量
Figure FDA00036412560400000216
即:
Figure FDA00036412560400000217
式(13)所示系统为一个非仿射系统,u无法直接解出,由于实际的控制系统均为离散系统,采用下式求解控制量:
Figure FDA00036412560400000218
其中,uk=1为上一时刻的控制量;
若控制器中的控制增益k1,k2满足如下不等式时,误差e在有限时间内收敛到零;
Figure FDA0003641256040000031
其中μ=-k1p12
Figure FDA0003641256040000032
Figure FDA0003641256040000033
p11=1,
Figure FDA0003641256040000034
a,Γ均为正常数且β>1。
2.如权利要求1所述的基于深度学习的四旋翼无人机地面效应补偿降落控制方法,其特征是,还包括验证步骤:根据谱归一化后的神经网络利普西斯常数小于1,结合Lyapunov分析方法证明误差e在有限时间内收敛到零。
3.如权利要求1所述的基于深度学习的四旋翼无人机地面效应补偿降落控制方法,其特征是,还包括进行地面效应数据收集以及实验验证步骤,地面效应数据收集方法如下:
地面效应的真值通过机载传感器获得,结合飞行控制器上的惯性测量单元IMU、飞行控制器的输出信号、以及动作捕捉系统的信息,解算无人机在整个飞行期间受到的地面效应,地面效应的数值由如下所示的公式进行解算:
fg=md-mg-Rfu (16)
其中,m为无人机质量,d为无人机的加速度,fu为总升力,且有
Figure FDA0003641256040000035
CT为升力系数,
Figure FDA0003641256040000036
为四个桨的转速平方之和,由于实际飞控为PWM输出,PWM波占空比与转速平方成正比,即
Figure FDA0003641256040000037
因此fu=[0,0,CTρpwm·PWMtotal]T,由于无人机悬停时产生的升力等于自身重力,基于此对CTρpw4进行标定,式(16)中d由无人机机载加速度计得到,R由机载陀螺仪数据得到,fu由PWM输出结合标定结果得出,因此由式(16)计算出无人机整个飞行期间的地面效应;
在获得地面效应真值后,使用神经网络进行拟合,结合设计的非线性控制器在无人机平台上进行实验,验证控制器的合理性和可行性。
CN202010546066.6A 2020-06-16 2020-06-16 基于深度学习的四旋翼无人机地面效应补偿降落控制方法 Active CN111650954B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010546066.6A CN111650954B (zh) 2020-06-16 2020-06-16 基于深度学习的四旋翼无人机地面效应补偿降落控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010546066.6A CN111650954B (zh) 2020-06-16 2020-06-16 基于深度学习的四旋翼无人机地面效应补偿降落控制方法

Publications (2)

Publication Number Publication Date
CN111650954A CN111650954A (zh) 2020-09-11
CN111650954B true CN111650954B (zh) 2022-07-12

Family

ID=72351481

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010546066.6A Active CN111650954B (zh) 2020-06-16 2020-06-16 基于深度学习的四旋翼无人机地面效应补偿降落控制方法

Country Status (1)

Country Link
CN (1) CN111650954B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113391555B (zh) * 2021-07-13 2022-08-19 天津大学 基于非线性扰动观测器的四旋翼无人机降落控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107479371A (zh) * 2017-07-03 2017-12-15 浙江工业大学 一种基于快速非奇异终端滑模的四旋翼无人机有限时间自适应控制方法
CN107491083A (zh) * 2017-09-11 2017-12-19 北京航空航天大学 一种基于饱和自适应滑模控制的四旋翼自主着船方法
CN109885074A (zh) * 2019-02-28 2019-06-14 天津大学 四旋翼无人机有限时间收敛姿态控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107479371A (zh) * 2017-07-03 2017-12-15 浙江工业大学 一种基于快速非奇异终端滑模的四旋翼无人机有限时间自适应控制方法
CN107491083A (zh) * 2017-09-11 2017-12-19 北京航空航天大学 一种基于饱和自适应滑模控制的四旋翼自主着船方法
CN109885074A (zh) * 2019-02-28 2019-06-14 天津大学 四旋翼无人机有限时间收敛姿态控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Adaptive Super-Twisting Control for Automatic Carrier Landing of Aircraft;ZIYANG ZHEN 等;《IEEE Transactions on Aerospace and Electronic Systems》;20190711;第56卷(第2期);984-997 *
Guanya Shi 等.Neural Lander:Stable Drone Landing Control Using Learned Dynamics.《2019 International Conference on Robotics and Automation (ICRA)》.2019,9784-9790. *
Neural Lander:Stable Drone Landing Control Using Learned Dynamics;Guanya Shi 等;《2019 International Conference on Robotics and Automation (ICRA)》;20190812;9784-9790 *
Vision-based Fast-terminal Sliding Mode Super Twisting Controller for Autonomous Landing of a Quadrotor on a Static Platform;Archit Krishna Kamath 等;《2019 28th IEEE International Conference on Robot and Human Interactive Communication (RO-MAN)》;20200113;1-6 *

Also Published As

Publication number Publication date
CN111650954A (zh) 2020-09-11

Similar Documents

Publication Publication Date Title
CN106249745B (zh) 四轴无人机的控制方法
CN109116860B (zh) 三旋翼无人机的非线性鲁棒控制方法
Castillo et al. Real-time stabilization and tracking of a four-rotor mini rotorcraft
CN112346470A (zh) 一种基于改进自抗扰控制的四旋翼姿态控制方法
CN106844887B (zh) 旋翼无人机的动力学建模方法及装置
CN111399531B (zh) 高超声速飞行器滑翔段制导与姿态控制一体化设计方法
CN107491083B (zh) 一种基于饱和自适应滑模控制的四旋翼自主着船方法
Zhu et al. Adaptive backstepping control for a miniature autonomous helicopter
CN108563897B (zh) 一种基于极值搜索算法的四旋翼无人机参数辨识方法
CN106043695B (zh) 一种油动多旋翼无人机定桨距变转速系统及控制技术
Amelin et al. Recursive identification of motion model parameters for ultralight UAV
CN112965512B (zh) 一种基于螺旋桨模型的无人机抗风控制方法
Ali et al. Attitude and altitude control of trirotor UAV by using adaptive hybrid controller
CN110723309A (zh) 一种四旋翼无人机转动惯量测量方法
Subudhi et al. Modeling and trajectory tracking with cascaded PD controller for quadrotor
CN114489094A (zh) 一种基于加速度反馈增强的旋翼无人机抗风扰控制方法
CN111650954B (zh) 基于深度学习的四旋翼无人机地面效应补偿降落控制方法
CN107678442B (zh) 一种基于双模型下的四旋翼自主着船控制方法
Santiaguillo-Salinas et al. Observer-based time-varying backstepping control for parrot’s ar. drone 2.0
Hegde et al. Transition flight modeling and robust control of a VTOL unmanned quad tilt-rotor aerial vehicle
CN111061282A (zh) 基于能量法的四旋翼无人机吊挂飞行系统控制方法
CN108279562A (zh) 一种基于滑模pid控制的飞行机械臂
Katigbak et al. Autonomous trajectory tracking of a quadrotor UAV using PID controller
Zhan et al. Control system design and experiments of a quadrotor
Alcan et al. Robust hovering control of a quadrotor using acceleration feedback

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