1.一种复杂风场下的桥梁检测特种无人机自主稳定方法，其特征在于，包括如下步骤：
步骤1：桥梁检测特种无人机采用多旋翼无人机，对多旋翼无人机进行动力学建模；
步骤1‑1：先进行如下假设：
(1)假设飞行区域的地面是一个平面，忽略地球自转，重力加速度设为常数；
(2)将多旋翼无人机机体结构和旋翼都视为刚体，忽略多旋翼无人机机体的弹性形变及振动；
(3)多旋翼无人机六个电机及旋翼均对称安装，且除正反极性外其他参数均相同；
(4)多旋翼无人机机体的质量分布均匀且质心与外形中心重合；
无人机在桥梁检测过程中受力包括：多旋翼无人机本体的重力、旋翼产生的升力、空气阻力；力矩包括：旋翼升力产生的空气动力矩、旋转产生的反扭力、空气阻力以及摩擦力产生的阻力矩；
步骤1‑2：建立机体坐标系和地面坐标系；
为确定无人机飞行中的位置，建立惯性坐标系E‑XYZ，地面坐标系原点是无人机起飞时的质心，X轴的正方向是起飞时无人机头的航向正方向，Y轴的正方向是起飞时无人机水平向左的方向，Z轴的正方向是起飞时无人机垂直向上的方向；
为确定无人机的姿态，建立机体坐标系B‑XbYbZb，机体坐标系原点是无人机的质心，Xb轴的正方向是机头的方向，Zb轴的正方向是垂直于机体平面向上，Yb轴的正方向由右手准则确定；
步骤1‑3：定义φ为无人机绕Xb轴转动的横滚角；θ为无人机绕Yb轴转动的俯仰角；ψ为旋T
翼无人机绕Zb轴转动的偏航角，则无人机姿态表示为Θ＝[φ θ ψ]；假设每个旋翼的升力和无人机平面垂直并指向上，此时每个旋翼的升力fi近似为：
其中，Ωi是旋翼转速；
此时六旋翼整体的推力表达为 方向沿着‑bZ的方向；总升力在惯性坐标系
3 T
中表示为‑fRe3∈R，其中e3＝[0，0，1]，R表示本体坐标系到惯性坐标系的旋转矩阵；
步骤1‑4：为了保证力矩平衡，在六旋翼中3个旋翼顺时针旋转，另外3个逆时针旋转；假设旋翼产生的扭矩与其升力成正比，此时得到每个旋翼的扭矩：
τi＝±cτfi    (2)
其中，cτ是常数，表示升力和扭矩的近似关系；
T
得到六旋翼总的升力f和总力矩M＝[M1 M2 M3]与各旋翼升力之间的关系：
其中d表示从螺旋桨中心到机体坐标系的距离；
步骤1‑5：最终得到多旋翼无人机的动力学模型表示为：
无人机位置动力学模型为：F1：
无人机姿态动力学模型：F2：
其中，x是无人机相对于惯性坐标系的位置，m为无人机质量，g为重力加速度，df为扰动矢量，J为转动惯量，ω为角速度，dτ为扰动力矩矢量；
步骤2：设计多旋翼无人机抗风扰控制算法；
步骤2‑1：采用基于几何控制理论的级联PID控制器作为无人机的控制器，由外控制环和内控制环构成，外控制环包括位置PID控制器和速度PID控制器，内控制环包括姿态PID控制器和角速度PID控制器，具体如下：
无人机的轨迹包括期望位置xd和期望航向ψd，位置跟踪误差定义为ex＝xd‑x；位置PID控制器设计为：
其中，vd表示期望速度，kPx、kDx和kIx分别表示位置PID控制器的比例、微分和积分参数；x表示实际位置；
给定期望速度后，定义速度跟踪误差ev＝vd‑v，速度PID控制器设计如下：
其中，Fd表示期望三轴力，kPv、kDv和kIv表示速度PID控制器的比例、微分和积分参数；v表示实际速度；
步骤2‑2：引入高增益加速度反馈算法，包括线加速度反馈和角加速度反馈；
在没有干扰的理想情况下无人机位置动力学模型中的力学方程表示为：
其中 是期望加速度，而考虑实际系统中存在扰动矢量df与设计的线加速度反馈项vf，则式(6)表示为：
式中，是无人机真实的线加速度；df和vf满足：
其中，Q(s)为滤波器；s为拉普拉斯算子；
结合以上式(7)和(8)得到设计的线加速度反馈项：
其中，a表示正常数；
此时，Fd与vf的和作为内控制环的输入，根据几何控制思想，Fd与vf被分解为模长f＝||Fd+vf||和方向 结合期望轨迹中的航向信息，计算出期望旋转矩阵Rd，根据几何控制，姿态跟踪误差表示为：
其中，映射V： 将SO(3)映射成为一个3维向量，表示向量的反对称运算的逆运算，SO(3)表示一种特殊正交群；则姿态PID控制器为：
其中，ωd表示期望加速度，kPR、kDR和kIR为姿态PID控制器的比例、微分和积分参数；
T
根据期望角速度ωd定义为角速度误差eω＝RRdωd‑ω，角加速度PID控制器为：
其中， 为期望角加速度，kPω、kDω和kIω为角加速度PID控制器的比例、微分和积分参数；ω表示实际角速度；
得到如下形式的角加速度反馈控制器：
其中，Md是期望力矩，是无人机机体角加速度矢量；vτ(t)表示角加速度反馈控制器；b为常数；
步骤3：设计GPS信号阶段性失锁条件下的状态估计算法；
步骤3‑1：采用GPS传感器和机载姿态确定装置获取无人机位姿信息，采用自适应卡尔曼滤波算法估计无人机的实时位置和姿态；
将无人机动力学模型统一转化为：
T
式中，f(x)＝[F1 F2]， w为系统高斯白噪声；
步骤3‑2：考虑到无人机作业时GPS信号阶段性失锁的情况，系统的观测方程分以下两种情况：
第一种情况：当GPS传感器输出信息可用时，此时无人机系统的测量方程为：
其中，xk，yk，zk， θk，ψk分别表示k时刻的各项值； 表示测量值；
考虑测量噪声，无人机系统的状态空间方程为：
式中，υ为测量传感器高斯白噪声，u表示控制输入，hk(.)表示厅时刻测量方程，f(.)表示系统状态函数，zk表示k时刻的测量函数；
第二种情况：当GPS信号失锁时：
由于桥梁检测无人机的飞行航向和速度阶段性保持不变，因此无人机的速度是先验已知的，即无人机飞行一定时间内的位移是先验已知的，因此，得到如下状态约束伪测量方程：
式中，vu为无人机飞行速度，tu为飞行时间， 为k+1时刻无人机位置状
态的期望值；xk‑1、yk‑1、zk‑1分别表示k‑1时刻的位置坐标；
步骤3‑2：令 将无人机的位置动力学模型变形为：
其中，
式(18)作为无人机另一个伪测量位置状态，即：
其中，wk表示k时刻的系统噪声； 表示k+1时刻无人机位置状态的期望值；
联立以上式(18)和(19)得到系统的观测方程为：
经离散后，系统的动力学模型为：
式中，ΔT为离散时间，wk‑1∈(0，Qk‑1)和vk∈(0，Rk)为离散后的高斯白噪声；Qk‑1表示系统噪声协方差矩阵，Rk表示测量噪声协方差矩阵；
步骤4：经过以上步骤，实现多旋翼无人机的自主稳定。
2.根据权利要求1所述的一种复杂风场下的桥梁检测特种无人机自主稳定方法，其特征在于，所述高增益加速度反馈具体如下：
首先利用欧拉‑拉格朗日方程建立通用的动力学方程：
其中 为广义位置向量，包括位置和角度； 为惯性矩阵，
为哥氏力项， 为重力项， 为摩擦阻力项， 为驱动力或力矩项，
为扰动力或力矩；
考虑无人机动力学模型中的任意一个自由度：
其中Jii是自由度的惯性系数；τi是控制输入；τui由各种耦合力以及扰动不确定项组成；
表示哥氏力项矩阵的第i行第j元素，Gii(q)表示重力项矩阵中第i个对角线元素，表示摩擦阻力项矩阵中第i个对角线元素；qj表示状态分量中第j个元素；
高增益加速度反馈设计如下：
其中ka是正常数，v是上一层控制器的输出；qi表示表示状态分量中第i个元素；
联立以上两式得到：
若ka足够大，即ka＞＞1且ka＞＞Jii，则上式近似为：
此时，得到了一个加速度跟踪器，直接跟踪给定加速度信号，同时消除了扰动项τui的影响；
引入前置滤波器，得到：
此时引入前置滤波器 得到：
无人机的任何一个自由度重新表示如下：
其中，A(s)、B(s)分别表示拉普拉斯变换后的多项式表达；
令ka→+∞得：
通过上式知不确定项τui经过一个滤波器Q(s)之后会对系统产生影响：
因为多旋翼无人机受到的扰动处于低频范围，因此选取Q(s)为高通滤波器的形式从而抑制低频扰动，选取：
其中a是正常数，系统表示为：
此时，扰动中的低频部分会被抑制，而高通滤波器的截止频率正好是a。
3.根据权利要求2所述的一种复杂风场下的桥梁检测特种无人机自主稳定方法，其特征在于，所述自适应扩展卡尔曼滤波算法步骤如下：
第一步：状态预测：
xk＝xk‑1+ΔTf(xk‑1)
其中，ΔT表示采样时间，f(xk‑1)表示k‑1时刻的状态方程；
第二步：协方差预测：
式中，In为单位矩阵；Pk‑1表示表示k‑1时刻的协方差矩阵，Qk‑1表示k‑1时刻的系统噪声方差矩阵， 表示k‑1时刻的预测状态；
第三步：计算卡尔曼滤波增益：
其中，h表示测量函数， 表示k时刻的预测状态；
第四步：更新状态变量，即计算状态变量的估计值：
第五步：更新状态估计误差协方差矩阵：
第六步：噪声协方差矩阵自适应调整策略：
新息定义为测量变量的实际值与预测值之间的误差，即：
新息的理论协方差矩阵：
T T
Ck＝E(vkvk)＝HkPkHk+Rk
新息的协方差矩阵计算方法如下：
式中，M为滑动窗口的长度；
当 时减小Rk；当 时，保持Rk不变；
定义的Rk调整因子为：
采用自适应协方差调整策略不更新Q和R，只在计算卡尔曼滤波增益和状态估计误差协方差矩阵给出Q和R修正量；
在更新滤波器增益时，自适应调整方法如下：
状态估计误差协方差矩阵的自适应调整方法为：
式中，ak表示波器增益自适应调整系数阵，κ表示协方差矩阵的自适应调整系数。