组合导航系统在四旋翼无人机上的实现
0 引言
近年来，随着智能化产业的迅速发展，无人机行业逐渐走进人们的视线。由于四旋翼无人机的结构简单可靠，具有垂直起降、定点悬停、稳定低速飞行和在小空间内飞行的能力［1］，所以受到广泛的关注并在多个领域得到了广泛的应用［2］。
在四旋翼无人机系统发展所面临的关键技术中，导航系统是其重要的组成部分之一，它承担着提供给飞行器位置、速度、姿态等参数状态数据的任务。目前常用的导航方法有：惯性导航、卫星导航、视觉导航以及它们的组合导航等。由于单一导航系统难以满足无人飞行器的发展要求［3］，于是产生了把两种或两种以上的不同导航系统以适当的方式组合在一起的组合导航，比单独使用任一导航系统时具有更高的系统性能［4］。本文是利用扩展卡尔曼滤波器将惯性导航系统和gps导航系统的组合来完成导航任务，建立了系统模型并将其在四旋翼无人飞行器上应用实现。
1 惯性导航系统
1.1 惯性传感器误差分析
影响导航系统工作精度的主要因素之一是惯性传感器误差的存在，为了减少传感器误差，需要先对传感器进行标定和误差补偿来保证导航系统的高精度和高性能。
1.2 惯性导航系统姿态矩阵计算
在导航系统的解算过程中，首先要处理的是坐标系的标定，通常有惯性坐标系（通常表示为i系） 、地理坐标系（通常表示为g系）、导航坐标系（通常表示为n系）、载体坐标系（通常表示为b系）等。本文在分析导航系统问题时是将导航坐标系设定为地理坐标系，地理坐标系采用x指东向、y指北向、z指天向的形式，惯导系统的原理方程可以用图1表示。
在导航系统中，最重要的计算就是姿态矩阵的求解，而常见的姿态矩阵算法主要有四元数法、欧拉角法、方向余弦法3种［5］。欧拉角法计算过程中方程式出现退化现象，方向余弦法通常计算量很大，因此经常采用的是四元数法求解姿态矩阵。但四元数法存在不可交换误差，为减小这种误差，本文采用的是bortz在1971年提出的等效旋转矢量算法［6］。
1.2.1 旋转矢量与姿态四元数的关系
设q（t+h）和q（t）分别为飞行器载体在t+h时刻和t时刻的姿态四元数，则：
1.2.3 旋转矢量的求解
通常，对于旋转矢量采用泰勒级数展开法来求解。若采用直线拟合角速度，等效旋转矢量的二子样算法为：
2 组合导航系统
虽然惯性导航系统能够连续工作并有效地提供姿态信息、位置信息和速度信息，但由于惯性传感器误差的积累，惯性导航系统的工作精度会随时间下降。而gps虽然可以提供长时间的误差为米级的高精度位置输出且用户设备成本较低，但由于gps信号会被遮挡或干扰，因此不能仅依赖gps提供连续导航参数［8］。
鉴于ins和gps系统的优缺点互补，将二者组合在一起，综合两系统的优势，能提供有效、长时、高精度、完整的导航参数。组合导航的一般结构如图2所示。
2.1 卡尔曼滤波器
kalman滤波是一种估计算法［9］，是导航系统中大多数状态估计算法的基础，如卫星导航结果的平滑、惯性导航系统的对准和标定、惯性导航系统与卫星或其他导航传感器间的组合导航等［10］，并已经成为从导航系统各种测量数据中获得最优估计结果的关键技术。
实际工程问题中，多数系统都是非线性的，因此采用扩展卡尔曼滤波器（ekf）。扩展卡尔曼滤波（ekf）为卡尔曼滤波的非线性形式［11］。
ekf的系统动态模型和观测模型分别为：
函数f（·）和h（·）不能直接应用在协方差中，取而代之的是可以计算其雅可比矩阵［12］，这个过程实质上将非线性的函数在当前估计值处线性化了。
2.2 系统模型与状态选择
本文将惯导系统和gps导航系统进行融合，在当地导航坐标系下建立系统模型。若卡尔曼滤波器估计的是相对于地球并投影到当地导航坐标系的姿态和速度误差，且估计的位置误差以纬度、经度和高度来表示，则状态向量变为：
式中，上标n表示投影到当地导航坐标系。
除了地球旋转角速率和陀螺测量值之外，姿态传播方程还引入了一个转移速率项，当地导航坐标系下的姿态误差为：
惯性导航系统的主要噪声源是由加速度计比力测量噪声导致的速度误差随机游走以及由陀螺角速度率测量噪声带来的姿态误差随机游走。如果单独估计加速度计和陀螺的动态零偏，则加速度计和陀螺零偏在运行时的变化可以近似为白噪声。
在ins/gps组合中，使用gps用户设备的测量输出与根据惯性导航参数预测的测量值之间的差来更新状态向量，采用哪些测量则依赖于组合结构。
3 仿真与实验结果
本次实验选用pixhawk飞控单元作为飞行器主控板，matlab作为仿真软件，分别对纯惯导系统与组合导航系统进行仿真，四旋翼无人飞行器搭建环境描述如下：
选用材料为尼龙加纤维的机架搭建飞行器，机架的对角轴距为35 cm；选取无刷电机型号为mt2312-960kv用于多旋翼飞行器中，提供动力输出；选取电池容量为5 000 ma，最大放电电流为30 a。该飞行器遥控器型号为乐迪at9，对应的接收器型号为2.4 g、9通道的r9d；地面站软件采用3dr推荐的专为px4/pixhawk设计的新的qgroundcontrol，在该环境下对飞行器进行校准和调试。
3.1 理想状态下的飞行轨迹
为了更好地理解导航系统的工作性能，在四旋翼无人飞行器飞行过程中，设定飞行顺序如图3所示。
具体过程可以描述为悬停时间为15 s；加速过程中，时间为10 s，加速度大小为0.5g（g为重力加速度）；爬升过程中，时间为25 s，仰角速度为2 °/s；俯冲过程中，时间为25 s，俯角速度为2 °/s；减速过程中，时间为5 s，加速度大小为1g；飞行周期为140 s，循环飞行3次，则无人飞行器理想状态下的飞行姿态参数如图4所示。
3.2 实验结果
（1）纯惯导系统导航系统姿态误差曲线如图5所示。
（2）组合导航系统的姿态误差曲线如图6所示。
由实验结果可知，理想飞行状态下，飞行器的横滚角和偏航角数值均为误差值（接近于0），只有俯仰角发生变化，如图4；若采用纯惯性导航，飞行器的姿态误差随着时间而增大，在飞行结束时，东向姿态误差（俯仰角误差）会达到1.1°左右，北向姿态误差（横滚角误差）达到约1.6°，天向姿态误差（偏航角误差）达到约1.5°，如图5；若采用ins/gps组合导航系统，飞行器东向（俯仰）和北向（横滚）的姿态误差均小于0.5°，而且在长时间内可以小于0.2°，而天向（偏航）姿态角误差虽然会在短时间内达到1°左右，但在其他时间段内都能保持在0.5°左右，如图6。因此，在相同的条件下，组合导航系统的性能优于纯惯导系统，误差较小，能够为四旋翼无人飞行器提供较为精确的导航信息。
4 结论
本文主要是分析了导航系统的原理和算法，在对传感器误差进行标定和补偿之后，采用了扩展卡尔曼滤波器将捷联惯性导航系统与gps导航系统融合，并成功应用于四旋翼无人飞行器上。仿真实验结果表明了系统方案在四旋翼无人飞行器上实现的可行性，且在长时间内可以提供较为准确的导航信息，误差较小，确保了无人飞行器的飞行效果。
参考文献
［1］ 马远超。四旋翼飞行器导航及控制技术研究［d］。哈尔滨：哈尔滨工程大学，2013.
［2］ milhim a b，zhang y m，rabbath c a.gain scheduling based pid controller for fault tolerant control of a quad-rotor uav［c］.aiaa infotech@aerospace，2010：1-13.
［3］ salychev o s.applied inertial navigation：problems and solutions［m］.moscow，russia：bmstu press，2004.
［4］ merwe r v d，wan e a.sigma-point kalman filters for integrated navigation［c］.proceedings of the 60th annual meeting of the institute of navigation（ion），2004：641-654.
［5］ 李晓峰。捷联惯导和组合导航的仿真研究［d］。西安：西安电子科技大学，2010.
［6］ bortz j e.a new mathematical formulation for strapdown inertial navigation［j］.aerospace and electronic systems，ieee transactions on，1971（1）：61-66.
［7］ 秦永元，张士邈。捷联惯导姿态更新的四子样旋转矢量优化算法研究［j］。中国惯性技术学报，2001，9（4）：1-7.
［8］ groves p d.principles of gnss，inertial， and multisensor integrated navigation systems［m］.artech house，2013.
［9］ mourikis a i，roumeliotis s i.a multi-state constraint kalman filter for vision-aided inertial navigation［c］.robotics and automation，2007 ieee international conference on.ieee，2007：3565-3572.
［10］ helal a，moore s e，ramachandran b.drishti：an integrated navigation system for visually impaired and disabled［c］.wearable computers，2001.proceedings.fifth international symposium on.ieee，2001：149-156.
［11］ haykin s s，haykin s s.kalman filtering and neural networks［m］.wiley interscience，2001.
［12］ julier s j，uhlmann j k，durrant-whyte h f.a new approach for filtering nonlinear systems［c］.american control conference，proceedings of the 1995.ieee，1995，3：1628-1632.
什么叫做电动机的过载保护?
功率放大器在无线电能传输系统的应用
Python与CANoe/CANoe4SW
工信部：加快存储器产业发展是补齐产业发展短板的必然要求
多协议微型天线设计案例分析
组合导航系统在四旋翼无人机上的实现
变得成熟稳重：雷蛇发布升级版Blade Stealth
一加6T外观渲染图和详细配置曝光
华为正式宣布：联合银联，正式启动手机POS产品首批应用试点合作
轻量级占用网络FlashOcc：主打实时性，高精度高效内存
SOLIDWORKS仿真计算是什么意思？
回顾梁家恩在第十届国际博士生论坛上的演讲：机遇与挑战并存
PCI-E接口的固态硬盘和SATA接口的区别是什么
在加入比特币潮流之前这几件事必须要知道
AI热并未给服务器市场带来预想中的振兴
7名亿嘉和智能“巡检卫士”在南方电网首座智能化运维的嵌入式附建变电站
OneNET多协议接入之HTTP协议的使用步骤
华为联合合作伙伴共助上海鲲鹏计算产业生态繁荣
智能制造设备的发展前景解析
与荣耀V9几乎同时推出的荣耀8青春版无人问津？
近年来，随着智能化产业的迅速发展，无人机行业逐渐走进人们的视线。由于四旋翼无人机的结构简单可靠，具有垂直起降、定点悬停、稳定低速飞行和在小空间内飞行的能力［1］，所以受到广泛的关注并在多个领域得到了广泛的应用［2］。
在四旋翼无人机系统发展所面临的关键技术中，导航系统是其重要的组成部分之一，它承担着提供给飞行器位置、速度、姿态等参数状态数据的任务。目前常用的导航方法有：惯性导航、卫星导航、视觉导航以及它们的组合导航等。由于单一导航系统难以满足无人飞行器的发展要求［3］，于是产生了把两种或两种以上的不同导航系统以适当的方式组合在一起的组合导航，比单独使用任一导航系统时具有更高的系统性能［4］。本文是利用扩展卡尔曼滤波器将惯性导航系统和gps导航系统的组合来完成导航任务，建立了系统模型并将其在四旋翼无人飞行器上应用实现。
1 惯性导航系统
1.1 惯性传感器误差分析
影响导航系统工作精度的主要因素之一是惯性传感器误差的存在，为了减少传感器误差，需要先对传感器进行标定和误差补偿来保证导航系统的高精度和高性能。
1.2 惯性导航系统姿态矩阵计算
在导航系统的解算过程中，首先要处理的是坐标系的标定，通常有惯性坐标系（通常表示为i系） 、地理坐标系（通常表示为g系）、导航坐标系（通常表示为n系）、载体坐标系（通常表示为b系）等。本文在分析导航系统问题时是将导航坐标系设定为地理坐标系，地理坐标系采用x指东向、y指北向、z指天向的形式，惯导系统的原理方程可以用图1表示。
在导航系统中，最重要的计算就是姿态矩阵的求解，而常见的姿态矩阵算法主要有四元数法、欧拉角法、方向余弦法3种［5］。欧拉角法计算过程中方程式出现退化现象，方向余弦法通常计算量很大，因此经常采用的是四元数法求解姿态矩阵。但四元数法存在不可交换误差，为减小这种误差，本文采用的是bortz在1971年提出的等效旋转矢量算法［6］。
1.2.1 旋转矢量与姿态四元数的关系
设q（t+h）和q（t）分别为飞行器载体在t+h时刻和t时刻的姿态四元数，则：
1.2.3 旋转矢量的求解
通常，对于旋转矢量采用泰勒级数展开法来求解。若采用直线拟合角速度，等效旋转矢量的二子样算法为：
2 组合导航系统
虽然惯性导航系统能够连续工作并有效地提供姿态信息、位置信息和速度信息，但由于惯性传感器误差的积累，惯性导航系统的工作精度会随时间下降。而gps虽然可以提供长时间的误差为米级的高精度位置输出且用户设备成本较低，但由于gps信号会被遮挡或干扰，因此不能仅依赖gps提供连续导航参数［8］。
鉴于ins和gps系统的优缺点互补，将二者组合在一起，综合两系统的优势，能提供有效、长时、高精度、完整的导航参数。组合导航的一般结构如图2所示。
2.1 卡尔曼滤波器
kalman滤波是一种估计算法［9］，是导航系统中大多数状态估计算法的基础，如卫星导航结果的平滑、惯性导航系统的对准和标定、惯性导航系统与卫星或其他导航传感器间的组合导航等［10］，并已经成为从导航系统各种测量数据中获得最优估计结果的关键技术。
实际工程问题中，多数系统都是非线性的，因此采用扩展卡尔曼滤波器（ekf）。扩展卡尔曼滤波（ekf）为卡尔曼滤波的非线性形式［11］。
ekf的系统动态模型和观测模型分别为：
函数f（·）和h（·）不能直接应用在协方差中，取而代之的是可以计算其雅可比矩阵［12］，这个过程实质上将非线性的函数在当前估计值处线性化了。
2.2 系统模型与状态选择
本文将惯导系统和gps导航系统进行融合，在当地导航坐标系下建立系统模型。若卡尔曼滤波器估计的是相对于地球并投影到当地导航坐标系的姿态和速度误差，且估计的位置误差以纬度、经度和高度来表示，则状态向量变为：
式中，上标n表示投影到当地导航坐标系。
除了地球旋转角速率和陀螺测量值之外，姿态传播方程还引入了一个转移速率项，当地导航坐标系下的姿态误差为：
惯性导航系统的主要噪声源是由加速度计比力测量噪声导致的速度误差随机游走以及由陀螺角速度率测量噪声带来的姿态误差随机游走。如果单独估计加速度计和陀螺的动态零偏，则加速度计和陀螺零偏在运行时的变化可以近似为白噪声。
在ins/gps组合中，使用gps用户设备的测量输出与根据惯性导航参数预测的测量值之间的差来更新状态向量，采用哪些测量则依赖于组合结构。
3 仿真与实验结果
本次实验选用pixhawk飞控单元作为飞行器主控板，matlab作为仿真软件，分别对纯惯导系统与组合导航系统进行仿真，四旋翼无人飞行器搭建环境描述如下：
选用材料为尼龙加纤维的机架搭建飞行器，机架的对角轴距为35 cm；选取无刷电机型号为mt2312-960kv用于多旋翼飞行器中，提供动力输出；选取电池容量为5 000 ma，最大放电电流为30 a。该飞行器遥控器型号为乐迪at9，对应的接收器型号为2.4 g、9通道的r9d；地面站软件采用3dr推荐的专为px4/pixhawk设计的新的qgroundcontrol，在该环境下对飞行器进行校准和调试。
3.1 理想状态下的飞行轨迹
为了更好地理解导航系统的工作性能，在四旋翼无人飞行器飞行过程中，设定飞行顺序如图3所示。
具体过程可以描述为悬停时间为15 s；加速过程中，时间为10 s，加速度大小为0.5g（g为重力加速度）；爬升过程中，时间为25 s，仰角速度为2 °/s；俯冲过程中，时间为25 s，俯角速度为2 °/s；减速过程中，时间为5 s，加速度大小为1g；飞行周期为140 s，循环飞行3次，则无人飞行器理想状态下的飞行姿态参数如图4所示。
3.2 实验结果
（1）纯惯导系统导航系统姿态误差曲线如图5所示。
（2）组合导航系统的姿态误差曲线如图6所示。
由实验结果可知，理想飞行状态下，飞行器的横滚角和偏航角数值均为误差值（接近于0），只有俯仰角发生变化，如图4；若采用纯惯性导航，飞行器的姿态误差随着时间而增大，在飞行结束时，东向姿态误差（俯仰角误差）会达到1.1°左右，北向姿态误差（横滚角误差）达到约1.6°，天向姿态误差（偏航角误差）达到约1.5°，如图5；若采用ins/gps组合导航系统，飞行器东向（俯仰）和北向（横滚）的姿态误差均小于0.5°，而且在长时间内可以小于0.2°，而天向（偏航）姿态角误差虽然会在短时间内达到1°左右，但在其他时间段内都能保持在0.5°左右，如图6。因此，在相同的条件下，组合导航系统的性能优于纯惯导系统，误差较小，能够为四旋翼无人飞行器提供较为精确的导航信息。
4 结论
本文主要是分析了导航系统的原理和算法，在对传感器误差进行标定和补偿之后，采用了扩展卡尔曼滤波器将捷联惯性导航系统与gps导航系统融合，并成功应用于四旋翼无人飞行器上。仿真实验结果表明了系统方案在四旋翼无人飞行器上实现的可行性，且在长时间内可以提供较为准确的导航信息，误差较小，确保了无人飞行器的飞行效果。
参考文献
［1］ 马远超。四旋翼飞行器导航及控制技术研究［d］。哈尔滨：哈尔滨工程大学，2013.
［2］ milhim a b，zhang y m，rabbath c a.gain scheduling based pid controller for fault tolerant control of a quad-rotor uav［c］.aiaa infotech@aerospace，2010：1-13.
［3］ salychev o s.applied inertial navigation：problems and solutions［m］.moscow，russia：bmstu press，2004.
［4］ merwe r v d，wan e a.sigma-point kalman filters for integrated navigation［c］.proceedings of the 60th annual meeting of the institute of navigation（ion），2004：641-654.
［5］ 李晓峰。捷联惯导和组合导航的仿真研究［d］。西安：西安电子科技大学，2010.
［6］ bortz j e.a new mathematical formulation for strapdown inertial navigation［j］.aerospace and electronic systems，ieee transactions on，1971（1）：61-66.
［7］ 秦永元，张士邈。捷联惯导姿态更新的四子样旋转矢量优化算法研究［j］。中国惯性技术学报，2001，9（4）：1-7.
［8］ groves p d.principles of gnss，inertial， and multisensor integrated navigation systems［m］.artech house，2013.
［9］ mourikis a i，roumeliotis s i.a multi-state constraint kalman filter for vision-aided inertial navigation［c］.robotics and automation，2007 ieee international conference on.ieee，2007：3565-3572.
［10］ helal a，moore s e，ramachandran b.drishti：an integrated navigation system for visually impaired and disabled［c］.wearable computers，2001.proceedings.fifth international symposium on.ieee，2001：149-156.
［11］ haykin s s，haykin s s.kalman filtering and neural networks［m］.wiley interscience，2001.
［12］ julier s j，uhlmann j k，durrant-whyte h f.a new approach for filtering nonlinear systems［c］.american control conference，proceedings of the 1995.ieee，1995，3：1628-1632.
什么叫做电动机的过载保护?
功率放大器在无线电能传输系统的应用
Python与CANoe/CANoe4SW
工信部：加快存储器产业发展是补齐产业发展短板的必然要求
多协议微型天线设计案例分析
组合导航系统在四旋翼无人机上的实现
变得成熟稳重：雷蛇发布升级版Blade Stealth
一加6T外观渲染图和详细配置曝光
华为正式宣布：联合银联，正式启动手机POS产品首批应用试点合作
轻量级占用网络FlashOcc：主打实时性，高精度高效内存
SOLIDWORKS仿真计算是什么意思？
回顾梁家恩在第十届国际博士生论坛上的演讲：机遇与挑战并存
PCI-E接口的固态硬盘和SATA接口的区别是什么
在加入比特币潮流之前这几件事必须要知道
AI热并未给服务器市场带来预想中的振兴
7名亿嘉和智能“巡检卫士”在南方电网首座智能化运维的嵌入式附建变电站
OneNET多协议接入之HTTP协议的使用步骤
华为联合合作伙伴共助上海鲲鹏计算产业生态繁荣
智能制造设备的发展前景解析
与荣耀V9几乎同时推出的荣耀8青春版无人问津？