1.一种基于非对称时不变障碍李雅普诺夫函数的四旋翼飞行器全状态受限反步控制方法，其特征在于，包括以下步骤：步骤1，建立四旋翼飞行器系统的动态模型，设定系统的初始值、采样时间以及相关控制参数，过程如下：

1.1确定从基于四旋翼飞行器系统的机体坐标系到基于地球的惯性坐标的转移矩阵T：其中φ，θ，ψ分别是四旋翼飞行器的翻滚角、俯仰角、偏航角，表示无人机依次绕惯性坐标系的各坐标轴旋转的角度；

1.2四旋翼飞行器平动过程中的动态模型如下：其中x，y，z分别表示四旋翼飞行器在惯性坐标系下的三个位置，Uf表示四旋翼飞行器的输入力矩，m为四旋翼飞行器的质量，g表示重力加速度，将式(1)代入式(2)得：

1.3四旋翼飞行器转动过程中的动态模型为：其中τx，τy，τz分别代表机体坐标系上各个轴的力矩分量，Ixx，Iyy，Izz分别表示机体坐标系下的各个轴的转动惯量的分量，×表示叉乘，ωp表示翻滚角速度，ωq表示俯仰角速度，ωr表示偏航角速度， 表示翻滚角加速度， 表示俯仰角加速度， 表示偏航角加速度；

考虑到无人机处于低速飞行或者悬停状态 ，姿态角变化较小，认为因此式(4)改写为：

联立式(3)和式(5)，得到四旋翼飞行器的动力学模型为：其中ux＝cosφsinθcosψ+sinφsinψ，uy＝cosφsinθsinψ-sinφcosψ；

1.4根据式(6)，定义φ，θ的期望值分别为：其中φd为φ的期望信号值，θd为θ期望信号值，arcsin为反正弦函数；

步骤2，在每一个采样时刻，计算位置跟踪误差及其一阶导数；计算姿态角跟踪误差及其一阶导数；设计位置和姿态角控制器，过程如下：

2.1定义z跟踪误差及其一阶导数：其中zd表示z的期望信号；

2.2定义q11：

2.3设计障碍李雅普诺夫函数V11：其中Ka1，Kb1为正常数：其中|e1|max为|e1|的最大值；

2.4求解式(10)一阶导数，得：其中 α1为虚拟控制量，其表达式为：其中k11为正常数；

将式(13)代入式(12)，得：

2.5定义q12：

2.6设计障碍李雅普诺夫函数V12：其中Kd1，Kc1为正常数：其中|s1|max为|s1|的最大值；

求解式(16)一阶导数，得：其中

将式(19)和式(6)代入式(18)，得：

2.7设计Uf：

其中k12为正常数；

2.8定义x，y跟踪误差分别为e2，e3，则有：其中xd，yd分别表示x，y的期望信号；

2.9定义q21，q31：

2.10设计障碍李雅普诺夫函数V21，V31：其中Ka2，Kb2，Ka3，Kb3为正常数：其中|e2|max为|e2|的最大值，|e3|max为|e3|的最大值；

2.11求解式(25)一阶导数，得：其中 α2，α3为虚拟控制量，其表达式为：其中k21，k31为正常数；

将式(28)代入式(27)，得：

2.12定义q22，q32：

2.13设计李雅普诺夫函数V22，V32：其中Kc2，Kd2，Kc3，Kd3为正常数：其中|s2|max为|s2|的最大值，|s3|max为|s3|的最大值；

求解式(32)一阶导数，得：其中

将式(35)和式(6)代入式(34)，得：

2.14设计ux，uy：其中k22，k32为正常数；

2.15定义姿态角跟踪误差及其一阶导数：其中j＝4，5，6，x4＝φ，x5＝θ，x6＝ψ，x4d表示φ的期望值，x5d表示θ的期望值，x6d表示ψ的期望值，e4表示φ的跟踪误差，e5表示θ的跟踪误差，e6表示ψ的跟踪误差；

2.16定义qj1：

2.17设计障碍李雅普诺夫函数Vj1：其中Kaj，Kbj为正常数：其中|ej|max为|ej|的最大值；

2.18求解式(40)一阶导数，得：其中 αj为虚拟控制量，其表达式为：其中kj1为正常数；

将式(43)代入式(42)，得：

2.19定义qj2：

2.20设计李雅普诺夫函数Vj2：其中Kdj，Kcj为正常数，满足其中|sj|max为|sj|的最大值；

求解式(46)一阶导数，得：其中

将式(48)和式(6)代入式(47)，得：

2.21通过式(49)设计τx，τy，τz：其中k42，k52，k62为正常数。

2.如权利要求1所述的一种基于非对称时不变障碍李雅普诺夫函数的四旋翼飞行器全状态受限反步控制方法，其特征在于：所述方法还包括以下步骤：步骤3，验证四旋翼飞行器系统的稳定性；

3.1将式(21)代入式(20)，得：

3.2将式(37)代入式(36)，得：

3.3将式(50)代入式(49)，得

3.4通过(51)，(52)，(53)可知四旋翼飞行器系统是稳定的。