计算机测量与控制.201 6.24(3) · 84 · Computer Measurement&Control 控秘技术 中图分类号:TP273 文献标识码:A 文章编号:1671—4598(2016)03—0084—04 DOI:10.16526/j.cnki.11--4762/tp.2016.03.023 多旋翼无人机位置控制系统设计 赵立峰,张 凯,王 伟 (南京信息工程大学信息与控制学院,南京 210044) 摘要:为解决小型无人机位置控制系统成本高、控制精度低以及控制性能不稳定等缺点,以自主研发的多旋翼无人植保机作为研究 对象来设计位置控制器;设计过程中,飞行器的位置建模误差,模型参数变化、GPS定位精度以及外部扰动(如风速)成为有待解决的 难题;为解决上述问题,首先设计了卡尔曼滤波器对反馈所得速度和位置信息进行滤波,并采用参考模型滑模控制理论设计位置控制器, 提高控制系统的鲁棒性;为了验证位置控制器的性能,进行了户外飞行实验,实验结果表明,无论是在悬停还是目标点追踪情况下,所 设计控制系统都具有良好的跟踪性与鲁棒性。 关键词:位置控制;卡尔曼滤波;参考模型;滑模控制;鲁棒性 Position Control System Design for Muti--Rotor Type UAV Zhao Lifeng,Zhang Kai,Wang Wei (College of Information and Control,Nanjing University of Information Science and Technology, Nanjing 210044.China) Abstract:In order to achieve the good performance about position controller of small unmaned aerial vehicle(UAV),the self—research plants protect UAV is selected as the target.During the design processing,we need to solve the error of modeling,the change of parame— ter,the low accuracy of GPS and external interference.A Kalman filter is designed tO realize data filtering and estimate the unobserved state value.Then the reference model and slide mode controller is combined to minimize the tracking error between the state variables of the vchicle and the state variables of the reference mode1.The experimental result shows the good performance of designed controller under various ex— perimental conditions. Keywords:position control;kalman filter;reference model;slide mode controller;robustness O 引言 近年来,随着小型多旋翼无人机技术的不断发展和成熟, 其在治安监控、航拍测绘、农业植保等方面已经展现出了广泛 的应用前景。尤其是在农业植保方面,多旋翼无人机作业方式 流传感器并不能够用来确定无人机的位置,因而通常采用 GPS来获取无人机当前位置和速度信息,并结合PID或LQR 等控制方法设计位置控制器。随着植保机不断进行药物喷洒, 机体本身的重量时刻在变化,模型参数也随之改变,而采用 PID或LQG等方法无法解决因参数改变带来的模型误差,使 控制效果逐渐变差。 针对上述问题,以自主研发的植保机作为试验平台,研 究多旋翼无人机位置控制设计方法。本文以多旋翼无人机位 置控制为目标,提出采用GPS结合无人机位置模型,利用卡 相比于传统的人工作业方式,大大节省了人力成本,提高了作 业效率,同时也保证了作业人员不受农药的侵害;相比于传统 固定翼飞机大面积喷洒的方式,多旋翼飞行器作业可以显著降 低飞行成本,且同时具有更强的抗风抗干扰能力。在农药喷洒 的效果上,多旋翼无人机高速旋转的桨叶产生的下洗气流可以 吹开植物表面的叶片,使农药均匀地作用于农作物整体,有效 尔曼滤波器对速度和位置信息进行滤波。位置控制器设计方 面,采用鲁棒性较强的参考模型滑模控制方法来设计无人机 的提高了农药的施用效率,体现了精细农业的要求[1 ]。因 此,多旋翼无人机植保机正逐渐成为无人机领域研究的主流项 目之一。 一位置控制系统。飞行过程中,控制系统的动态性能由设计的 参考模型决定,而滑模控制可实现对外部干扰及模型参数摄 动的鲁棒性。 般来讲,无人机的位置控制系统可分为室内外两种。室 内位置控制通常是通过超声波和光流传感器来获取无人机的位 1 系统设计 多旋翼无人机本身具有6个自由度,分别为前后、左右、 偏航方向上的线运动和角运动,且各自由度之间具有耦合,因 此,该飞行系统实际上是一个多输入多输出的非线性系统。考 虑到实际飞行过程中,姿态角变化范围较小,可忽略各自由度 之间的耦合,将无人机控制系统简化为多个单输入单输出 系统。 置信息。得益于精准的反馈信息,室内位置控制器的设计较为 简单,且可以产生很好的控制效果。户外飞行时,超声波及光 收稿日期:2015—09—15;修回日期:2015—11—18。 作者简介:赵立峰(1991一),男,江苏盐城人,硕士研究生,主要从事 自动化控制方向的研究。 张凯(1965一),男,山东泰安人,教授,硕士研究生导师,主要从事 无人机的控制结构整体如图l所示,控制主要分为内环姿 智能信息检测与控制方向的研究。 王伟(1972一),男,黑龙江双鸭山人,教授,硕士研究生导师,主要 从事自动化控制方向的研究。 态控制和外环位置控制。其中,内环姿态控制是实现无人机飞 行的基础。根据运动学原理,可建立了从姿态角到位置的数学 模型,从而内外环之间的控制可通过姿态角联系为一体。整个 第3期 赵立峰,等:多旋翼无人机位置控制系统设计 · 85 · 控制流程为:外层位置控制器根据目标位置信息和反馈得到的 无人机当前位置的误差,计算出位置控制量,由于所建模型是 描述姿态角到位置之间的关系,故可将该位置控制量作为姿态 控制器的目标姿态角度。姿态控制器根据目标角度与当前角度 之间的误差,计算出角度控制量,驱使无人机完成相应位姿 (位置和姿态)运动。悬停情况下,目标点即为当前点,位置 误差为零,无人机将不产生动作。本文设计的位置控制器是建 立在内环姿态控制已完成的基础之上的,具体的姿态控制器设 计方法可参考 。 图1控制结构图 位置控制器具体设计框图如图2所示,其中,参考模型根 据目标输入产生各状态变量(速度,位置)的理想变化轨迹, 决定了控制器对目标输入的动态响应性能。参考模型的选取难 点在于既要考虑到模型的响应速度,又要使各状态量的变化轨 迹符合实际情况。考虑到GPS反馈所得信息可能存在噪声干 扰,设计卡尔曼滤波器对速度和位置数据进行滤波。控制器设 计方面,将目标输入,参考状态变量,及参考状态变量与无人 机实际状态之间的误差作为滑模控制器的输入,通过计算得到 控制量,使实际系统和参考模型的状态量之间的偏差逐渐趋近 于零。 图2位置控制器框图 1.1位置建模参考模型设计 多旋翼无人机机体采用对称结构,前后和左右方向上的位 置模型相同,因此建模过程中只考虑左右方向上的位置模型。 根据力学原理,可得出无人机飞行过程中左右方向上的运动方 程为 一嘲 一是∈跚 (1) 式(1)中,m为无人机总质量, 代表横向加速度, 为 横滚角,k∈表示风阻系数,s表示机身所受空气阻力面积, 为左右方向上速度。考虑到实际飞行中,多旋翼无人机飞行速 度一般较慢,且受风面积较小,故所受空气阻力可忽略,则式 (1)可简化为 a一 (2) 目标横滚角到实际横滚角之间的关系可表示为一阶惯性函 数,即 G 一 一志 ㈣ 其中,r为目标角度,k为增益系数, 表示延时系数。 结合(2)式和(3)式,可得目标角度到加速度的传递函数 口(s)一 r( ) (4) 对(4)式进行积分可得目标角度到速度的传递函数 (s)一 sr(s) (5) 将(5)式转化为状态空间方程可得 [ 广 ㈤ 为了方便拟合模型未知系数,将GPS获取的经度方向速 度 和纬度方向速度 转换为导航坐标系内的前后方向速度 和左右方向速度 ,转换方程 一"Oncos(yaw)一73esin(yaw) , Vy—Vn sin(yaw)+ cos(y ̄w) 其中:yaw为无人机的偏航角。将转换后的数据导入到 MATLAB系统辨识模块中,获取未知系数,得到速度模型。 图3为在同一姿态指令下,所得速度模型仿真结果与实际飞行 数据的对比,两条曲线基本重合,从而验证了该速度模型的可 靠性。对(6)式进行积分,可得如下位置模型,并选择速度 量作为输出结果。 1 0 0 kg 』 Td Ax+ 1 O O O O 1 O O [0 1O Q (8) 1.2位置建模参考模型设计 参考模型的设计要求模型可以稳定,准确,快速的追踪目 标输入,同时状态变化轨迹符合实际系统的的物理特性。根据 (8)式所得实际系统模型,可设参考模型为 z 一A X+B r ,n、 , 一( z 8 …… ! ~一 6 n ^ 喜2 i。 童_2 热 l_4 \f一 W W 一 三一6 ; V —8 0 l0 20 30 柏 50 60 70 80 90 t/s 图3速度模型辨识结果图 · 86 · 计算机测量与控制 一第24卷 Se (19) 参考模型中的状态变量直接对应系统实际模型的状态变 量。其中 —c,r为目标输入,为了使实际系统状态变量和 参考状态变量的变化趋势保持一致,引入状态误差e一-z—z , 结合(8)式与(9)式可得 e—A e+(A—A )e+ 一B r (1O) 对式(19)求导可得 一SA P 一SB (K1z+K2 r—H) (2O) 当状态变量之间的误差收敛于滑模面时,切换函数及其导 数均等于0,即 = 一0,结合式(19)可得等效控制输入 “ 一一(SB )~SA e +Kl 4-K2 r (11) (21) 设A --A=BK1,B =BKz,则式(1O)可写为 e:A e—B(K1z+ r—H) 考虑到外部扰动,建模误差及模型参数变化带来的影响 D亩,即 假设时间趋于无穷时,系统的状态变量将趋于稳定,则其 导数为0,即A 32+B r一0,则 z- 一一A B r (12) 将式(12)代人式(9)可得 Y 一~( A BK2 r (13) 系统稳定情况下,Y 可追踪目标输入r,即Y 一r,则由 式(13)可求出K2 K2一(--e B) (14) 将K2带人到B BK2,可得 B 一B(一C A B) (15) A 的值可通过MATLAB仿真实验调节得到。 即 广 1一 2一口3] A 一l l o1 0 1 0 o l J (16) 因为位置控制主要以稳定为主,不需要特别快的响应速 度,且廉价GPS所得的位置精度较低,反馈的速度数据也存 在噪声,控制太强可能会使无人机在目标点附近来回振动。当 选取a1—2.8, 一6.4, 一5.3时,系统的输出和状态变化 较为理想。图4为所设计参考模型的阶跃响应仿真结果图,追 踪时间为5秒左右,没有出现超调,且追踪波形较为平滑,满 足设计要求。 l- l①Reference Input l’ l②Reference ̄lodel Output 1· 一/ 一 ’ n 名 / 呈n 童n , n , , / 图4参考模型仿真结果图 1.3滑模控制器设计 设计滑模控制器的目的是使实际状态与参考模型状态之间 的偏差0,为加强控制器的追踪性能,使控制结果不出现稳态 误差,我们引入位置误差的积分 £一 — (17) 结合式(11),可得扩展系统为 [A m A e 其中“ 一一(Klz+K2 r ; 一A e +B “ +I (22) 其中:D为已知量,面为未知量,则等效控制输入为 一一(sB )SA (A e + ) (23) 将“ 带入式(18)可得 e =A e 一B (SB )一 S(A e + )+D一 (J—B (sB )S)A e +(J—B (SB ) S)腑 (24) 假设矩阵D满足:扩展矩阵[曰 ,D]和矩阵B 的秩相 同,则D可以表示为 D=D D (25) 则式(23)中受面影响的部分可表示成 (J—B (sB ) S) 一 (B 一B (SB )SB )D 亩一0 (26) 因此,滑模控制器的控制性能不受模型误差及外部干扰的 影响,具有较好的鲁棒性能。 为了使状态误差快速准确的收敛于滑模面,选择最优反馈 增益作为切换面s,即 S=B P (27) 其中矩阵P为下述黎卡提方程的解 PAs+AsTP一耶 B P+Q一0 (28) 为了抑制噪声及其它非线性因素影响,设计相应的非线性 控制部分得到非线性控制量,防止系统在切换面上下发生抖动 “ 一K ,( ) (29) 式中,K -为增益系数,f(d)为切换函数且一般为符号 函数。考虑到选取符号函数会使非线性控制量产生突变,影响 系统控制性能,现设计如下平滑函数取代 )一 (3o) 8为平滑函数权重,通常一个比较小的数字。结合之前的 等效控制,则整个控制器的产生的控制量为 “一“ +“ (31) 2实验结果及分析 为了验证设计的位置控制器性能,进行室外飞行实验。飞 行实验包括3个部分:1)卡尔曼滤波效果验证实验;2)悬停 和目标跟踪性能验证实验;3)鲁棒性实验。实验过程中,采 用无线通信模块Xbee—Pro900与地面站进行实时通信。 2.1卡尔曼滤波器验证实验 为了获取无噪声或噪声较小的反馈数据,需对所设计的卡 尔曼滤波器进行验证。图5给出GPS原始速度数据与滤波器 所得速度数据的对比,两条曲线基本吻合,且滤波器所得数据 更为平滑,验证了卡尔曼滤波器的正确性。 2.2悬停和跟踪性能实验 飞行实验过程中,通过操纵遥控手柄给出不同速度目标值 来测试所设计的控制器性能。图6中两条曲线分别表示所给出 第3期 赵立峰,等:多旋翼无人机位置控制系统设计 · 87 · = 苫 图5卡尔曼滤波结果 的速度目标值和无人机实际飞行速度,实验结果表明,在悬停 (目标速度为O)和有目标值情况下,实际速度都可以准确追 踪目标速度。输出的噪声是因为GPS的测量存在误差,因此 会导致无人机的实际速度有些波动,但波动较小,在可接受范 围内。图7给出了目标位置和实际位置信息的对比。由于目标 位置无法通过手柄实际给出,两条曲线分别是对目标速度积分 所得的目标位置和是通过对GPS经纬度换算所得的位置信息, 从而认为设计的位置控制器满足期望。 5 一0 5 t/0 图6跟踪性能速度验证结果图 l0 。 点 言一 0 昌20 10 。 暑 三一10 昌 一20 tls 图7跟踪性能位置验证位置结果图 2.3鲁棒性验证实验 为了验证所设计位置控制器的鲁棒性,往机身所载药箱内 注入5公斤水,且飞行过程中不断喷洒水来改变机身的质量。 实验结果如图7所示,由实验结果可知,改变整体重量引起模 型参数改变的情况下,控制器控制效果基本不受影响。 4 l : 。 “… ^: : ; z 。 营一2 — r :’; : 4 5 。 吾 5 0 20 40 60 80 100 120 140 160 180 200 t/s 图8鲁棒性验证速度结果图 3 总结 本文主要研究了多旋翼无人机位置控制器设计方法。根据 力学原理,对飞行过程中速度和位置建模,然后通过GPS获 取的速度和位置信息,对所得模型进行验证,并在此基础上设 计参考模型滑模控制器。最后通过室外飞行实验对该控制器性 能进行验证,实验结果表明,所设计的控制系统具有良好的追 踪性及鲁棒性。 参考文献: [1]Peng Y S,wu C F,Ma S H,Reseearch on key problems in as— signed recovery of UAV using parachute I-A-I.TENCON 2013— 2013 IEEE Region 10 Conference[C].2013. [2]Orsag M,Korpela C,Oh P,Modeling and control of MM—uAV: Mobile manipulating unmanned aerial vehicle[J].Journal of Intel— ligent&Robotic Systems,2013,69:227—240. I-3]Lee T,Leok M,McClamroch N,Nonlinear robust tracking control of a quadrotor uAv on SE[J].Asian Journal of Control,2013, 15:391—408. [4]Hemakumara P,Sukkarieh S,Learning uAV stability and control derivatives using Gaussian processes[J].Robotics。IEEE Trans— actions on,2013,29:813—824. [5]Yu H,Beard R,A vision—based collision avoidance technique for micro air vehicles using local——level frame mapping and path plan—- ning口].Autonomous Robots,2013,34:93—109. [6]Zhu S,Wang D,Low C,Ground target tracking using uav with in— put constraints[J].Journal of Intelligent&Robotic Systems, 2O13,69:417—429. [7]Wang W,Ma H,Xia M,Attitude and Altitude Controller Design for Quad—Rotor Type uAVs EJ].Mathematical Problems in Engi neering,2013.