四旋翼飞行器具有垂直起降、空中悬停、飞行姿态灵活变换等可操控性，其结构简单并且负载能力较大.其可以通过传感器收集各类信息，已被广泛应用于摄影摄像、军事勘察、天气监测以及救援活动等多个领域.四旋翼的飞行控制研究已发展近二十年，已取得丰硕的成果，如文献[1-3]中实现了四旋翼较高精度的姿态和位置控制，以及实现各种高难度避障、倒飞和复杂轨迹跟踪等动作.四旋翼飞行器是具有强耦合特性的典型非线性系统，一旦发生故障，飞行稳定性就会急剧下降，甚至导致无人机失控.如何保证四旋翼无人机在发生故障的情况下仍能得到有效控制，已成为四旋翼无人机领域的一个热点问题.

目前针对四旋翼的故障检测与隔离已经有诸多优秀成果，文献[4-5]中指出关于观测器的故障检测与隔离技术已经相当成熟，故本文的重点在研究成果较少的单翼损坏故障检测之后的飞行器控制上.对于四旋翼故障下的控制目前已经有很多研究成果，如文献[6]中的滑模控制算法，文献[7]中采用的反步法控制算法，文献[8]中采用的模型预测控制算法，大部分的研究偏向于传感器故障和执行器效率损失.对于执行器损坏故障研究甚少，文献[9]中提出的单翼损坏下的四旋翼飞行控制器采用鲁棒反馈线性化来设计，过程繁琐且计算任务繁重.

针对单翼损坏下的四旋翼飞行器，本文采用反步法来简化控制器设计.文献[10]中指出，四旋翼飞行器为欠驱动、强耦合、非线性、不稳定系统，在机体执行机构出现故障时，其欠驱动度增加，更加不利于飞行器的灵活控制.通过对四旋翼结构和运动机理的分析，滚转角和俯仰角的小变化可以引起位置的较大变化，会使机体失去平衡甚至引起机体撞击地面.而偏航能力的丧失则不会对飞行器的安全造成威胁，而且不会影响飞行器正常的位置控制.所以，采用反步法设计的单翼损坏下的四旋翼飞行控制器，以牺牲四旋翼的偏航控制为前提，解决因为单翼损坏而使得四旋翼欠驱动程度增加的难题，实现此故障下的四旋翼飞行控制.

1 四旋翼飞行器建模 1.1 无故障四旋翼飞行器动力学建模

四旋翼飞行器靠4个电机的转动带动螺旋桨旋转产生升力.文献[11]指出，通过改变不同电机的转速四旋翼可以实现俯仰、滚转、偏航等动作.如图 1所示，四旋翼的每个旋翼能产生升力和力矩，两对旋翼(1, 3)和(2, 4)按相反的旋转方向旋转，以抵消每对旋翼旋转带来的旋转力矩.

参照文献[12]，分别搭建如图 1所示两个右手坐标系，建立以四旋翼重心为原点的机体坐标系B={o_b, x_b, y_b, z_b}和标准地球坐标系E={o_e, x_e, y_e, z_e}.记四旋翼的姿态角定义在E系中，为Φ=(φ, θ, ψ)，其中φ、θ、ψ分别为四旋翼飞行器的滚转角、俯仰角和偏航角.则B系可以由E系按x-y-z的旋转顺序得到.

四旋翼在机体坐标系下的角速度ω=[ω₁, ω₂, ω₃]^T与在姿态角Φ=(φ, θ, ψ)的关系如下:

$ \omega = \left[ {\begin{array}{*{20}{c}} 1&0&{ - \sin \theta }\\ 0&{\cos \varphi }&{\cos \theta \sin \varphi }\\ 0&{ - \sin \varphi }&{\cos \varphi \cos \theta } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\dot \varphi }\\ {\dot \theta }\\ {\dot \psi } \end{array}} \right]. $

(1)

参照文献[13]，将4个旋翼的电机模型定义为

$ {F_i} = {K_f}\mathit{\Omega }_i^2,\left( {i = 1,2,3,4} \right). $

(2)

式中：F_i为各旋翼所受到的升力; K_f为升力系数; Ω_i为旋翼旋转速度.即在忽略空气扰动情况下，视旋翼升力与其旋转速度成正比.

由此可以得到，飞行器在机体坐标系下的总升力为

$ {F_b} = \left( {0,0,\sum\limits_{i = 1}^4 {{F_i}} } \right). $

(3)

根据电机模型，考虑四旋翼飞行器的运动原理，可以建立输入向量为

$ \left( {\begin{array}{*{20}{c}} {{u_1}}\\ {{u_2}}\\ {{u_3}}\\ {{u_4}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{u_f}}\\ {{u_p}}\\ {{u_q}}\\ {{u_r}} \end{array}} \right) = \left[ {\begin{array}{*{20}{c}} 1&1&1&1\\ 0&{ - l}&0&l\\ { - l}&0&l&0\\ d&{ - d}&d&{ - d} \end{array}} \right]\left( {\begin{array}{*{20}{c}} {{F_1}}\\ {{F_2}}\\ {{F_3}}\\ {{F_4}} \end{array}} \right). $

(4)

在假设飞行器为理想完全对称的刚体的前提下，可近似得到机体的转动惯量矩阵为

$ \mathit{\boldsymbol{I}} = {\rm{diag}}\left[ {{I_x},{I_y},{I_z}} \right]. $

(5)

结合式(1)~(5)，由牛顿运动定律，可以得到四旋翼飞行器在地球坐标系下的动力方程如下:

$ \left\{ \begin{array}{l} \ddot x = \frac{{{u_1}}}{m}\left( {\cos \psi \sin \theta \cos \varphi + \sin \psi sin\varphi } \right) - \frac{{{k_1}\dot x}}{m},\\ \ddot y = \frac{{{u_1}}}{m}\left( {\sin \psi \sin \theta \cos \varphi - \sin \varphi \cos \psi } \right) - \frac{{{k_2}\dot y}}{m},\\ \ddot z = \frac{{{u_1}}}{m}\left( {\cos \theta \cos \varphi } \right) - g - \frac{{{k_3}\dot z}}{m},\\ \dot p = qr\frac{{{I_y} - {I_z}}}{{{I_x}}} + \frac{{{I_r}}}{{{I_x}}}q{\mathit{\Omega }_r} + \frac{l}{{{I_x}}}\left( { - {F_2} + {F_4}} \right) - \frac{{l{k_p}}}{{{I_x}}}p,\\ \dot q = pr\frac{{{I_z} - {I_x}}}{{{I_y}}} - \frac{{{I_r}}}{{{I_y}}}p{\mathit{\Omega }_r} + \frac{l}{{{I_x}}}\left( { - {F_1} + {F_3}} \right) - \frac{{l{k_q}}}{{{I_y}}}q,\\ \dot r = pq\frac{{{I_x} - {I_y}}}{{{I_z}}} + \frac{c}{{{I_z}}}\left( {{F_1} - {F_2} + {F_3} - {F_4}} \right) - \frac{{{k_r}}}{{{I_z}}}r,\\ \ddot \varphi = \dot \theta \dot \psi \frac{{{I_y} - {I_z}}}{{{I_x}}} + \frac{l}{{{I_x}}}{u_2} - \frac{{l{k_p}}}{{{I_x}}}\dot \varphi ,\\ \ddot \theta = \dot \varphi \dot \psi \frac{{{I_z} - {I_x}}}{{{I_y}}} + \frac{l}{{{I_y}}}{u_3} - \frac{{l{k_q}}}{{{I_y}}}\dot \theta ,\\ \ddot \psi = \frac{c}{{{I_z}}}{u_4} - \frac{{{k_r}}}{{{I_z}}}\dot \psi . \end{array} \right. $

1.2 单翼损坏情况下的故障建模

当飞行器的执行器发生一个旋翼处发动机损坏故障时，基于四旋翼飞行器的结构对称性，不失一般性地假定编号为2的翼受损，无法提供正常的升力.如图 2所示，由于2翼的损坏导致其不能提供升力，相当于F₂=0.由式(4)可知此时有

$ {u_p} = l{F_4}, $

可重新建立输入向量为

$ \mathit{\boldsymbol{u = }}{\left( {{u_1},{u_2},{u_3}} \right)^{\rm{T}}} = {\left( {{u_f},{u_q},{u_r}} \right)^{\rm{T}}} = \left[ {\begin{array}{*{20}{c}} 1&1&1\\ { - l}&l&0\\ d&d&{ - d} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{F_1}}\\ {{F_3}}\\ {{F_4}} \end{array}} \right]. $

(6)

当飞行器处于悬停状态下时，为了保证飞行器处于水平状态，由于四旋翼的对称结构和运动机理，需同时取F₄=0.由式(6)可知，此时的u_r=dF₁+dF₄=u_f*d>0，此式显示不可变号，即悬停状态下不可任意控制偏航角速度.由于四旋翼飞行器突发单翼损坏时，其首要工作是使飞行器尽快达到稳定状态而不至于撞机，故只要能使飞行器的滚转角和俯仰角以及位置得到控制即可保证，所以此时的偏航控制并不是很重要.又由于四旋翼本身为欠驱动系统，4个旋翼转速的输入要给出空间6个自由度的输出，故在一个旋翼损坏时，其欠驱动程度增加，所以此时舍弃偏航控制无非是最好的选择.

对于故障下的飞行器状态空间描述，选择状态量为

$ \begin{array}{*{20}{l}} {x = {{\left( {{x_1},{x_2},{x_3},{x_4},{x_5},{x_6},{x_7},{x_8},{x_9},{x_{10}},{x_{11}},{x_{12}}} \right)}^{\rm{T}}} = }\\ {\;\;\;\;\;\;{{\left( {x,y,z,\dot x,\dot y,\dot z,\varphi ,\theta ,\psi ,p,q,r} \right)}^{\rm{T}}},} \end{array} $

则在单翼损坏故障下的飞行器状态空间表述为：

$ {{\dot x}_1} = {x_4}, $

$ {{\dot x}_2} = {x_5}, $

$ {{\dot x}_3} = {x_6}, $

$ {{\dot x}_4} = \frac{{\cos {x_7}\sin {x_8} + \sin {x_7}\sin {x_9}}}{m}{u_1} - \frac{{{k_1}}}{m}{x_4}, $

$ {{\dot x}_5} = \frac{{\cos {x_7}\sin {x_8}\sin {x_9} - \sin {x_7}\cos {x_9}}}{m}{u_1} - \frac{{{k_2}}}{m}{x_5}, $

$ {{\dot x}_6} = \frac{{\cos {x_7}\cos {x_8}}}{m}{u_1} - \frac{{{k_3}}}{m}{x_5} - g, $

$ {{\dot x}_7} = {x_{10}} + {x_{11}}\sin {x_7}\tan {x_8} + {x_{12}}\cos {x_7}\tan {x_8}, $

$ {{\dot x}_8} = {x_{11}}\cos {x_7} - {x_{12}}\sin {x_7}, $

$ {{\dot x}_9} = \frac{{{x_{11}}\sin {x_7} + {x_{12}}\cos {x_7}}}{{\cos {x_8}}}, $

$ {{\dot x}_{10}} = \frac{1}{{{I_x}}}\left( { - {k_1}{x_{10}} - {x_{11}}{x_{12}}\left( {{I_z} - {I_y}} \right) + \frac{l}{2}\left( {{u_1} - \frac{{{u_3}}}{d}} \right)} \right), $

$ {{\dot x}_{11}} = \frac{1}{{{I_y}}}\left( { - {k_2}{x_{11}} - {x_{10}}{x_{12}}\left( {{I_x} - {I_z}} \right) + {u_2}} \right), $

$ {{\dot x}_{12}} = \frac{1}{{{I_z}}}\left( { - {k_3}{x_{12}} + {u_3}} \right), $

(7)

考虑四旋翼的对称性，选取I_x=I_y.

2 控制器设计

本文的控制器设计建立在单翼损坏的情况视为已经发生的基础上，不做与故障的检测和隔离相关的工作和分析，重点放在故障发生后的系统控制上.

经分析，在2翼损坏的情况下，四旋翼飞行器的欠驱动程度增加，为了使四旋翼在不可避免的单翼损坏情况下能尽量减小损失和保证飞行器不增加外部损坏，本文设计思路为牺牲偏航控制，依赖飞行器无故障的剩余3个旋翼来达到近悬停状态和实现轨迹跟踪.为了解决欠驱动和实现解耦，在四旋翼飞行器的控制过程中，文献[14-15]中，均采用划分子系统来设计控制器，本文也将使用这种方法，将系统划分为内、外环分层控制.设计的内环控制器主要控制飞行器姿态，使单翼损坏的四旋翼飞行器得到稳定的近悬停状态.这里所提到的近悬停状态是指飞行器保证恒定的高度和较稳定的姿态角，同时以一定速度绕垂直轴转动的状态.设计的外环控制器是飞行器位置控制的主体，建立在内环控制器效果良好的基础之上，外环控制器通过期望位置指令信号解算出内环的跟踪姿态信息.

由于飞行器的水平位置向量(x, y)的加速度为机体所受总升力u_f的水平分量提供，而此水平分力的大小受滚转角φ和俯仰角θ的影响.单翼损坏下的四旋翼飞行器的飞行控制器结构如图 3所示，水平位置的控制也依赖已经形成的内环控制，所以内环的必要工作是获取滚转角和俯仰角的期望值，以此使飞行器达到期望的水平位置.机体的垂直位置分量z的控制，是通过调节机体总体升力来实现的，需要关注的是，此时还涉及到偏航角旋转速度r.对应于z、u_f以及r三者之间的关系以及水平位置控制中涉及到的姿态控制，故在新的位置指令下，内环将从原有的平衡状态达到另一个平衡状态.

2.1 内环控制器设计

内环主要用来控制飞行器的滚转角φ、俯仰角θ和偏航角转动速度r，为了解决姿态的耦合问题，选用反步法来设计控制器，使其实现李雅普诺夫意义下的稳定.

设定控制输入为u=[F₁, F₃, F₄]^T，选取状态量为

$ {x_1} = \left[ {\begin{array}{*{20}{c}} \varphi \\ \theta \\ {\int {r{\rm{d}}t} } \end{array}} \right],{x_2} = \left[ {\begin{array}{*{20}{c}} p\\ q\\ r \end{array}} \right]. $

四旋翼飞行器的姿态动力学，即式(7)，用上述输入和状态可写成状态空间形式为

$ \left\{ \begin{array}{l} {{\dot x}_1} = {f_1}{x_2}\\ {{\dot x}_2} = {f_2}u + {g_c} \end{array} \right., $

式中各矩阵如下：

$ {f_1} = \left[ {\begin{array}{*{20}{c}} 1&{\sin \varphi \tan \theta }&{\cos \varphi \tan \theta }\\ 0&{\cos \varphi }&{ - \sin \varphi }\\ 0&0&1 \end{array}} \right], $

$ {f_2} = \left[ {\begin{array}{*{20}{c}} 0&0&{\frac{l}{{{I_x}}}}\\ { - \frac{l}{{{I_x}}}}&{\frac{l}{{{I_x}}}}&0\\ {\frac{c}{{{I_z}}}}&{\frac{c}{{{I_z}}}}&{ - \frac{c}{{{I_z}}}} \end{array}} \right],{g_c} = \left[ {\begin{array}{*{20}{c}} {\frac{{{I_y} - {I_z}}}{{{I_x}}}qr - \frac{{l{k_4}}}{{{I_x}}}p}\\ {\frac{{{I_z} - {I_x}}}{{{I_y}}}pr - \frac{{l{k_5}}}{{{I_y}}}q}\\ { - \frac{{{k_6}}}{{{I_z}}}r} \end{array}} \right], $

且注意到f₂可逆，当$ \left( {-\frac{\pi }{2} < \varphi < \frac{\pi }{2}} \right) $和$ \left( {-\frac{\pi }{2} < \theta < \frac{\pi }{2}} \right)$时f₁可逆，记两者的逆矩阵分别为f₁^－1和f₂^－1.

为了使姿态达到期望值，可以通过设计一系列的虚拟控制器来稳定每个状态，用反步法设计的过程如下.

Step 1  考虑虚拟系统

$ {{\dot x}_1} = {f_1}{v_1}, $

确定该虚拟系统的跟踪误差为

$ {e_1} = {x_{1d}} - {x_1}, $

选取Lyapunov函数为

$ {V_1} = \frac{1}{2}e_1^{\rm{T}}{e_1}, $

(8)

对式(8)求导有

$ {{\dot V}_1} = e_1^{\rm{T}}{{\dot e}_1} = e_1^{\rm{T}}\left( {{{\dot x}_{1d}} - {{\dot x}_1}} \right) = e_1^{\rm{T}}\left( {{{\dot x}_{1d}} - {f_1}{v_1}} \right). $

为了使此虚拟系统稳定，虚拟控制器选为

$ {v_1} = f_1^{ - 1}\left( {{A_1}{e_1} + {{\dot x}_{1d}}} \right), $

其中A₁为正定阵，此时可得

$ {{\dot V}_1} = - e_1^{\rm{T}}{A_1}{e_1} < 0. $

(9)

Step 2  考虑虚拟系统

$ {{\dot x}_2} = {f_2}u + {g_c}, $

确定该虚拟系统的误差式为

$ \begin{array}{l} {e_2} = {v_1} - {x_2} = f_1^{ - 1}\left( {{A_1}{e_1} + {{\dot x}_{1d}} - {f_1}{{\dot x}_1}} \right) = \\ \;\;\;\;\;\;\;f_1^{ - 1}\left( {{A_1}{e_1} + {{\dot e}_1}} \right). \end{array} $

因此可得

$ {{\dot e}_1} = - {A_1}{e_1} + {f_1}{e_2}. $

选取Lyapunov函数为

$ {V_2} = \frac{1}{2}\sum\limits_{i = 1}^2 {e_i^{\rm{T}}{e_i}} , $

求导有

$ \begin{array}{l} {{\dot V}_2} = e_1^{\rm{T}}{{\dot e}_1} + e_2^{\rm{T}}{{\dot e}_2} = e_1^{\rm{T}}\left( { - {A_1}{e_1} + {f_1}{e_2}} \right) + e_2^{\rm{T}}\left( {{{\dot v}_1} - {{\dot x}_2}} \right) = \\ \;\;\;\;\;\;\; - e_1^{\rm{T}}{A_1}{e_1} + e_2^{\rm{T}}\left( {f_1^{\rm{T}}{e_1} + {{\dot v}_1} - {f_2}u - {g_c}} \right). \end{array} $

为了使整个内环系统稳定，内环控制器选为

$ u = f_2^{ - 1}\left( {{A_2}{e_2} + f_1^{\rm{T}}{e_1} + {{\dot v}_1} - {g_c}} \right), $

其中A₂正定时有

$ {{\dot V}_2} = - \sum\limits_{i = 1}^2 {e_i^{\rm{T}}{A_i}{e_i} < 0} . $

(10)

综上所述，内环控制器为

$ \left\{ \begin{array}{l} {v_1} = f_1^{ - 1}\left( {{A_1}{e_1} + {{\dot x}_{1d}}} \right),\\ u = f_2^{ - 1}\left( {{A_2}{e_2} + f_1^{\rm{T}}{e_1} + {{\dot v}_1} - {g_c}} \right). \end{array} \right. $

2.2 外环控制器设计

假定内环是Lyapunov稳定的，且工作在平衡状态，取其工作点为φ→φ_d, θ→θ_d, r→r_d，且φ_d和θ_d是小角度，达到近悬停状态.

为了保证近悬停状态，必取F₄=0，此时有$ {u_f} = \frac{{{u_r}}}{d} $，且升力用来平衡机体重力，即u_f=mg，故此时可得在平衡工作点时：

$ {u_r} = mgd, $

$ {r_d} = \frac{{mgd}}{{{k_6}}}, $

由此反推，当r→r_d时，u_r=k₆r_d，即$ {u_f} = \frac{{{k_6}{r_d}}}{d} $.

对于水平控制，摘取其位置动力学方程如下

$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {\ddot x}\\ {\ddot y}\\ {\ddot z} \end{array}} \right] = \frac{{{u_f}}}{m}\left[ {\begin{array}{*{20}{c}} {\cos \psi }&{\sin \psi }&0\\ {\sin \psi }&{ - \cos \psi }&0\\ 0&0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\sin {\theta _d}\cos {\varphi _d}}\\ {\sin {\varphi _d}}\\ 1 \end{array}} \right] + \\ \;\;\;\;\;\;\;\;\left[ {\begin{array}{*{20}{c}} { - \frac{{{k_1}\dot x}}{m}}\\ { - \frac{{{k_1}\dot y}}{m}}\\ { - \frac{{{k_1}\dot z}}{m} - g} \end{array}} \right]. \end{array} $

设水平控制系统的输入为

$ \mathop u\limits^ \gg = \left[ {\begin{array}{*{20}{c}} {\frac{{\sin {\theta _d}\cos {\varphi _d}}}{m}{u_f}}\\ {\frac{{\sin {\varphi _d}}}{m}{u_f}}\\ {\frac{1}{m}{u_f}} \end{array}} \right]. $

取状态量为

$ \mathop {{x_1}}\limits^ \gg = \left[ {\begin{array}{*{20}{c}} x\\ y\\ z \end{array}} \right],\mathop {{x_2}}\limits^ \gg = \left[ {\begin{array}{*{20}{c}} {\dot x}\\ {\dot y}\\ {\dot z} \end{array}} \right]. $

由式(7)，可得系统状态空间表达为

$ \left\{ \begin{array}{l} \mathop {{x_1}}\limits^{\dot \gg } = \mathop {{x_2}}\limits^ \gg ,\\ \mathop {{x_2}}\limits^{\dot \gg } = {f_\psi }\mathop u\limits^ \gg + {g_k}. \end{array} \right. $

其中：

$ {f_\psi } = \left[ {\begin{array}{*{20}{c}} {\cos \psi }&{\sin \psi }&0\\ {\sin \psi }&{ - \cos \psi }&0\\ 0&0&1 \end{array}} \right],{g_k} = \left[ {\begin{array}{*{20}{c}} { - \frac{{{k_t}\dot x}}{m}}\\ { - \frac{{{k_t}\dot y}}{m}}\\ { - \frac{{{k_t}\dot z}}{m} - g} \end{array}} \right], $

且注意到f_ψ恒可逆.

定义外环虚拟系统为$ \mathop {{v_1}}\limits^ \gg $，则两个跟踪误差为：

$ \mathop {{e_1}}\limits^ \gg = \mathop {{x_{1d}}}\limits^ \gg - \mathop {{x_1}}\limits^ \gg , $

$ \mathop {{e_2}}\limits^ \gg = \mathop {{v_1}}\limits^ \gg - \mathop {{x_2}}\limits^ \gg . $

用反步法设计外环控制率，推导过程类似内环设计，故在此不做累述，其中，Lyapunov函数设计为：

$ \mathop {{V_1}}\limits^ \gg = \frac{1}{2}\mathop {e_1^{\rm{T}}}\limits^ \gg \mathop {{e_1}}\limits^ \gg , $

$ \mathop {{V_2}}\limits^ \gg = \frac{1}{2}\mathop {e_1^{\rm{T}}}\limits^ \gg \mathop {{e_1}}\limits^ \gg + \frac{1}{2}\mathop {e_2^{\rm{T}}}\limits^ \gg \mathop {{e_2}}\limits^ \gg . $

当设计外环控制率为

$ \left\{ \begin{array}{l} \mathop {{v_1}}\limits^ \gg = {B_1}\mathop {{e_1}}\limits^ \gg + \mathop {{x_{1d}}}\limits^{\dot \gg } ,\\ \mathop u\limits^ \gg = f_\psi ^{ - 1}\left( {\mathop {{e_1}}\limits^ \gg + {B_2}\mathop {{e_2}}\limits^ \gg + \mathop {{v_1}}\limits^{\dot \gg } - {g_k}} \right), \end{array} \right. $

可得：

$ \mathop {{V_1}}\limits^{\dot \gg } = \mathop {e_1^{\rm{T}}}\limits^ \gg \mathop {{e_1}}\limits^{\dot \gg } = - \mathop {e_1^{\rm{T}}}\limits^ \gg {B_1}\mathop {{e_1}}\limits^ \gg < 0, $

(11)

$ \mathop {{V_2}}\limits^{\dot \gg } = \mathop {e_1^{\rm{T}}}\limits^ \gg \mathop {{e_1}}\limits^{\dot \gg } + \mathop {e_2^{\rm{T}}}\limits^ \gg \mathop {{e_2}}\limits^{\dot \gg } = - \mathop {e_1^{\rm{T}}}\limits^ \gg {B_1}\mathop {{e_1}}\limits^ \gg - \mathop {e_2^{\rm{T}}}\limits^ \gg {B_2}\mathop {{e_2}}\limits^ \gg < 0, $

(12)

式中B₁、B₂为正定阵，使得外环为Lyapunov稳定.

由于系统的位置控制过程是：外环控制器解算出姿态信息作为内环的输入，通过位置与姿态的耦合关系来实现位置跟踪.在内环渐进稳定的基础上如果外环也是Lyapunov稳定的，则整个系统的输出是稳定的.又根据式(9)~(10)、式(11)~(12)可得，内、外环分别是Lyapunov稳定的.因此，给出有界的系统期望值时，在内外环作用下，四旋翼飞行系统将保持渐进稳定，并且能追踪期望位置和轨迹.本文的仿真结果也说明了这一点.

3 结果和分析

在MATLAB/Simulink中搭建控制系统，对本文提出的反步控制算法进行仿真.实验采用的是微型四旋翼飞行器，模型结构参数见表 1.

设定飞行器的起始位置坐标(x₀, y₀, z₀)=(0, 0, 0)，初始姿态角(φ₀, θ₀, ψ₀)=(0, 0, 0)，给定期望位置为(x_d, y_d, z_d)=(2, 3, 1).

位置和姿态的仿真结果如图 4所示.可以看出俯仰角和滚转角在6 s内稳定于稳定值，滚转角和俯仰角的超调量分别为7°和2°，基本满足飞行器的控制性能要求.偏航角速度在波动后趋于恒值.位置控制环无超调量，响应时间约为5 s，在完成位置跟踪之后，飞行器可以在期望位置最终到达近悬停飞行状态.

上述仿真针对单个目标点的跟踪能力，为了考查控制器对于连续多个目标点的简单路径的跟踪能力，设定多点直线轨迹，验证其跟踪和降落能力.设计4个点为A(0.5, 0.5, 0.6)、B(1.0, 0.5, 0.6)、C(1.0, 1.0, 0.6)和D(0.5, 1.0, 0.6).期望路径设定为飞行器以路径A→B→C→D→A飞行，然后降落至坐标原点.仿真效果如图 5所示，飞行器的飞行轨迹与期望大致重合，简单路径跟踪效果较好.

4 结论

1) 本文针对单翼损坏下四旋翼飞行器，用反步法设计了姿态控制器和位置控制器.该控制器牺牲飞行器的偏航控制，来有效处理失去一个翼的升力情况下，四旋翼的欠驱动程度增加导致姿态和位置不可控的问题.

2) 仿真结果说明，偏航控制能力的失去不会对四旋翼的安全造成威胁，也不会对飞行器的正常轨迹跟踪造成较大影响，且保证飞行器能在以恒定速度绕垂直轴转动的情况下进行稳定飞行，同时还能以较快速度跟踪简单的期望轨迹.

3) 本文采用的算法对解决在单翼损坏下的四旋翼位置控制和轨迹跟踪具有较强的理论价值，比已有的此类控制器的设计过程更加简单.通过进一步实践探讨，此算法将可用于实际环境中四旋翼突发单翼损坏故障下飞行器紧急降落问题.