![Rank: 8](images/default/star_level3.gif) ![Rank: 8](images/default/star_level3.gif)
- UID
- 1023166
- 性别
- 男
- 来自
- 燕山大学
|
![](http://images.eccn.com/silabs/silicon_chip_980x60_202203.jpg)
关键字:无人机 组合导航 卡尔曼滤波 扩展卡尔曼滤波
四轴飞行器是一种特殊结构的无人机(Unmanned AerialVehicle,UAV),可以只通过调节四个旋翼的转速,实现对飞行姿态的控制,它可为海上、废墟等不适合人员进入的环境,提供侦察、救援、绘图等服务。UAV 多数使用GPS 或惯性导航系统,然而,在GPS 信号不足或室内环境下,惯性导航系统误差又随时间积累,所以本身不能长时间地提供高精度的导航,需要辅助导航 。视觉导航方法,即由机载的视觉传感器获得周围环境的图像信息,然后通过图像算法及相机的位置标定解算出载体的位置,有更大的发展空间,它成本低、质量轻,且易于实现。
视觉导航方法主要集中分为两类:一类是以视觉伺服为主的导航方法,它基于视觉传感器,直接利用图像特征控制飞行器到指定位置,通常有两个控制环,外环控制实时特征与目标特征间的误差,估计飞行器的理想速度,然后,通过内环控制电机转速到理想速度,最后,使特征的图像匹配到目标位置。另一类是视觉辅助导航方法,即融合视觉传感器和惯性传感器给出载体位置,融合的方法多见Kalman 滤波方法。
本文基于四轴飞行器的数学模型,研究扩展卡尔曼滤波( Extended Kalman Filter, EKF) 和不敏卡尔曼滤波(Unscented Kalman Filter,UKF)两种非线性滤波算法在四轴飞行器视觉辅助导航方法中的应用情况,并对比Kalman 方法,UKF 和EKF 算法融合四轴飞行器多传感器信息的效果。
1 四轴飞行器数学模型
四轴飞行器的动力学模型满足如下假设:四轴飞行器是均匀对称的刚体结构;惯性系的原点与飞行器的几何中心一致;飞行器所受阻力及重力与飞行高度无关,保持不变;四轴各向拉力与电机的转速成正比。四轴飞行器的结构模型如图1 所示。
![](http://image.ednchina.com/2015/UAVEKFUKF01.JPG)
图1 四轴飞行器的数学模型示意图
图1 中涉及到惯性坐标系- E 系和载体坐标系- B 系,φ,θ,ψ 分别表示惯性坐标系下载体的欧拉角,定义Θ =[φ θ ψ]T 。则根据坐标系变换原则,载体系到惯性系的变换矩阵R 如式1 所示,其中c 表示余弦函数, s 表示正弦函数。
![](http://image.ednchina.com/2015/UAVEKFUKF02.JPG)
四轴飞行器电机的推力F 与它的转速ω 的平方成正比:
![](http://data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAcFBQYFBAcGBQYIBwcIChELCgkJChUPEAwRGBUaGRgVGBcbHichGx0lHRcYIi4iJSgpKywrGiAvMy8qMicqKyr/2wBDAQcICAoJChQLCxQqHBgcKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKir/wAARCAAdAawDASIAAhEBAxEB/8QAGwABAQEAAwEBAAAAAAAAAAAAAAUGAQMEBwj/xAAtEAABBAIABQIFBAMAAAAAAAABAAIDBAURBhITITFBYQcUIlFxMkJSoRdTcv/EABUBAQEAAAAAAAAAAAAAAAAAAAAB/8QAHhEBAQACAgIDAAAAAAAAAAAAAAERQSHwMVFhcdH/2gAMAwEAAhEDEQA/AP0iiL5px/m2VfiFhIJBdkq4+jZyN6GnzkytAEcbXBvpzOJ79vpCbk77NWvpaKbh8UzG/NSsdKPnJGymF8rnthPI1pa3Z8bBPbXcnsqSUnIihW+KoaluWu7E5mUxuLS+HHSPY73DgNEKxVsC3UisNjliEjQ4MmYWPb7Fp7g+yayO1FieJ7sH+RuH68gLjSrz3nMj11JXHUUcYHrsvcdeNs2fG1Rk47xEE0jbZfTZVrR2b8lrUYpCT9DJO/6yf2jfj3G05me94N4WsjfGOqGY17Fk70Iq7OZ57bPkgeAfJ/vS7KN2DJY+veqPL4LMTZY3EEba4bB0e47FZ+7lbHEM2WwOCtQ0chRdD1pLdd00ckMjeYFoY9h7gOHkEaPsVoKMMtfHwQzuidJHGGuMMfTZsD9rdnlH2Gz+Unjk+neiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgISANk6ARZ7M2zb4rxWAB1DNDLdsj/ZHGWNDPwXSNJ+4aR4KTm4Gh8jYQkAbJA/KicR4CrnIOXLXbUGPhjc50da1JW27+bnxuaSAPTevU77a8HAlO/J8P61XiKSW4ZOq1jrXeR9cvd0ucnuXdPl3vv9+6ap6apFD4UuzWKNulakdNPjLklN0rjsyNbpzHE+p5HN2fvtXESUWbocLzV+OMvxBetw2W368VaKDoEGGOMuOtlx3suJPYei0iJvK6wlZDJNh4gxeO3p1nqy7EvKdMaBrWvq7vHbY8b9FVXGgSDobHgrlBm73A9DIX5rcuUz8T5ncxZXzlqKNv8Ayxsga0ewCu0qrKNGGrHJNIyFgY188rpXuA9XPcSXH3J2u9EnEwfLH5vgefK5i5lIckyC1NJS6LnQc3RjrydQs/UCeZxOyNdtDuvHb+G8tixlwMsx9bK2a9yaOxV6jjNCGcoLg4AxkxtJboHzoja3iKy2eO+PyF5Y44xvCNHJZW1Y+cyeZtwttWWnoAbLY2hg+rlaxpJAJJ89+62K4IBGiAfyuVNAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgKHmaUkWcx2drROldUbJXsRsG3GCTlJLR6lrmMOvOubWzoK4ibyM7xXh89mmVWYLMU8dFG/nmZZomwJv4g6kboA99d9kBezFx5PF4aR/EGQblLYcXOfWq9JpHgNbHzOP8AZ2SqyJoSOHMXLjcfM+3r5u7Ykt2A07DXvPZoPrytDW79eXaroiD/2Q==)
其中: i = 1,2,3,4 分别代表4 个电机, k 为推力系数。则四轴飞行器的总推力为:
![](http://image.ednchina.com/2015/xUAVEKFUKF04.JPG.pagespeed.ic.IsC4RixQAC.jpg)
定义P = [X Y Z] T表示载体在惯性系下的位置矢量,根据牛顿第二定律知,载体的运动加速度与所受力成正比:
![](http://data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAcFBQYFBAcGBQYIBwcIChELCgkJChUPEAwRGBUaGRgVGBcbHichGx0lHRcYIi4iJSgpKywrGiAvMy8qMicqKyr/2wBDAQcICAoJChQLCxQqHBgcKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKir/wAARCAAqAgIDASIAAhEBAxEB/8QAHAABAAIDAQEBAAAAAAAAAAAAAAUGAwQHAgEI/8QAOxAAAgEEAQIEAgYGCwEAAAAAAAECAwQFEQYSIRMiMUEHFDJRYWJxgRUWIzORoQgXJDRCU2SSlLHh4v/EABYBAQEBAAAAAAAAAAAAAAAAAAABAv/EACgRAQABAQYFBAMAAAAAAAAAAAABEQIhMUGh8DJRYbHRA3GB4RKRwf/aAAwDAQACEQMRAD8A/SIAAAGplVkXirhYSVtG/cf2Du1J0lL73T31+AG2DRtY5VZS6d7O0di4U/lo0oyVVS0+vrbemt61r8zeAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFG+IauL7L8Vw2PvbqzuL3JeJUqW1aUJKhSi6lTaT002oR7p/SGcQZTK8lQ+J/MbjgXB7jP2tK3rzt5wj4FdteL1SUdJp9n337+hRL+te1ON8ovLHJ5FQu85DF4Wmryq/BqdUKMqqfV1PzqctNteVvXdmH433MOVZfhHCrSvKpTyt+q9ecXpulDy7/Pc2n90Um1+MRnMa0ntK1izM1yrp9rFT+KOTxvI+I4vkeKt4vlFuqkHazl1Ws2lqMk/pLzJbTXv2OnFZsuA4a15JRz907nJZS3peDb3N7V6/Ah37Qikox9X3S39pYbr+6Vv3i8kv3X0/T/AA/b9RbUxT9/TNmJrf0+2UHMPE/1HxF/4/8A8HRcZ3xVq93L/Yx73S1VfZfT+99f2ily5o3mecrca4bk83bqg5WFvKv0V99M+lb6ez7N+iZzep8aMxi+D8c5Xm8LafJ5u68D5W3qyValFuWpraaltR3rt6ruZv6R+Vq0fh7a4Gzf9qzt9StoR3puKfU/y2or8yA4zhqGG+PdvxHlTlk7SyxlOtx5XU24W7jFdXTDtFvcZ92m/Kh6d838+0VnQt3WbuXeaRq6l8Q+V3fEOL08njaFG6uatzRtqNtW2vFlUkopJp9n337+huu65DQzOOoVIWF1aV6k4XVSjTnCVvqm5J95NPbSXfXqVP4g0P1j+JfDOMqc40qNWrlrmVOWpQVKPTTf++X8i58Yw1bBYX5S6u6l5cSr1a1S4qz6pTc5trb0vSLivT2Fnhr1nx3qlrGI6b0olwa2Qd8rCo8Srd3el4auXJU/Xvvp7+myCpy514kfFo8e6Nrq6atfevfXlJC5LMA9ae/Q4hR/SF7xyle2eYyVO75DyfwsW1d1JK2toVO8opyaknToyfm2vN6d3tF9qm8YjfsThXeEz/Hbwc+47l6Fhyzl1xdZGtTxVjKnbwVzWnUip0qXiV6u23r95FPXbyoisbjuQcmtsTSuby/pyr3FevmruncSpxlQfV4dClp9k04acdNKL20292IrT41Ju3v2dWB5pwVOnGnHeopJbbb7fa/U9EAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACGz+Wq2dawxti18/kqzpUm1tUoRXVUqNfdiu33nFe5MlTzVN0Pidxq9rPVCpa3lnFt9lVl4c4r840p/wERWYje5wMpWbqpWNk5VqrVKhDc6tWe9JLu23/2VDiHLclyDnPJrC8oRt7LHwtZWdNx1NxqQlLqn7pvS8vt+OyV5bjs3k7S1pYKpYKMK6qXFK+6+irFLcY+Tvrq02vfWvTaKRwb9Zf66OX/pL9E66LL535fxP8mXR4e/57/Itm+ZrvAnh3zXq/yFfDcksVWnKpjspP5fzd/Ar9LcNfdkk1r2ko69WTpU+fQnd0MFj7dv5i4zNrOHT6qNKfizf4dMH/Etgjh+fHkni+PPgABAME7K1q31K8qW1Kd1RhKFKtKCc4RlrqSfqk9Lf4IzgCOnx7D1MbHHzxdo7OFXxo0PBj0KfU5dSWtb6m3v62zWueG8avMnHI3eAxta9hrpualrB1I67LUtbWiaAAAAAABEZXinHs7dU7rNYTH5CvSWqdW6toVJQW99m12PV3xfA5C8tbu+w1jcXNnFRt61W3jKdFL0UW1tL8CVAwEdLj2GlmVl5YuzeSS6VeOjHxdfV1a3okQAAAA+NKUXGSTTWmn7mnSwuLoQsoUcda044/fyijRilb7XS+jt5eza7ezN0AVvkfGIXWEyVPBWlnQvMhKLupOmo/Mw6l4kJS12coJx37b2ROH+H2Hx/IrHK4bjdvx6drKTqSoTUZ14uEo+G4wbi47kpbb3uK7e5egIuwJvigAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGnlMXQy9i7a56o6lGpTqQep0pxe4zi/rT/wDexuADzSVSNGEa01OoopSko9Kk/d69j7pJtpLb9X9Z9AGlHGQeYlkq8vFrRp+FQTWlRg9OSX2tpbf2Je3fdAAAAD//2Q==)
则载体在惯性系下的线运动方程为:
![](http://image.ednchina.com/2015/UAVEKFUKF06.JPG)
定义p,q,r 表示四轴飞行器在载体系中的运动角速度,则欧拉角与运动角速度的关系为:
![](http://image.ednchina.com/2015/UAVEKFUKF07.JPG)
角运动方程为:
![](http://image.ednchina.com/2015/xUAVEKFUKF08.JPG.pagespeed.ic.JYVEe8EC6F.jpg)
其中: I 为惯性矩阵, MR为旋转力矩, M 为飞行器所受力矩,由于四轴飞行器是对称的刚体,且其重心与几何中心重合,所以惯性矩阵为对角阵,对角线元素为Ix、Iy、Iz , 考虑沿水平方向的旋转效应,旋转力矩为:
![](http://image.ednchina.com/2015/xUAVEKFUKF09.JPG.pagespeed.ic.Kg5eoE7wZa.jpg)
飞行器所受力矩为:
![](http://image.ednchina.com/2015/UAVEKFUKF10.JPG)
则将式(6)展开为:
![](http://image.ednchina.com/2015/UAVEKFUKF11.JPG)
定义控制量如下:
![](http://image.ednchina.com/2015/UAVEKFUKF12.JPG)
2 非线性滤波算法
2.1 扩展卡尔曼滤波
EKF 算法利用泰勒展开的一次项来对非线性方程作线性化处理, 再结合经典的卡尔曼滤波进行滤波估计。EKF算法简单,易于实现。
设系统的状态方程为:
![](http://image.ednchina.com/2015/UAVEKFUKF13.JPG)
其中: x 为系统状态矢量, z 为系统观测量; w(k),v(k) 为不相关零均值的高斯噪声; Q(k),R(k) 分别表示其方差矩阵。φ(x),r(x) 表示泰勒展开算子,即线性化过程。
测量更新:
![](http://image.ednchina.com/2015/UAVEKFUKF14.JPG)
其中J 为雅可比算子。
2.2 不敏卡尔曼滤波(UKF)
不敏卡尔曼滤波是基于采样方法的近似法,用带有权值的样本集来近似目标的状态后验概率密度,采用的是确定的样本点,因而避免由线性化而导致的误差。设系统的状态方程为:
![](http://image.ednchina.com/2015/xUAVEKFUKF15.JPG.pagespeed.ic.zxYnJ-H8PH.jpg)
其中: xk ∈Rn 为系统状态矢量, zk ∈Rm 为测量量, qk-1 ∈Rn和rk-1 ∈ Rm是高斯噪声,且分别qk-1 ~ N(0,Qk-1 ),rk-1 ~N(0,Rk-1 ) 分布。初始状态x0 的均值和协方差分别记为m0和P0 。UKF 算法的根本是无迹变换(Unscented Transform,UT),其过程如下:
先根据系统状态的均值x 和方差矩阵Pxx 选择一组样本点x0 ,…,xi 和相应权值w0 ,…,wi :
![](http://image.ednchina.com/2015/UAVEKFUKF16.JPG)
其中参数k 是标量的常数。
由样本点x0 ,…,xi 经过非线性函数得到一组新的样本点z0 ,…,zi ,zi = h(xi )(i = 0,…,2n)。
用这些带有权值的新样本点计算观测状态的均值和方差
![](http://image.ednchina.com/2015/xUAVEKFUKF17.JPG.pagespeed.ic.FwVYewA_AL.jpg)
UKF 滤波算法同样包含时间更新与测量更新两个过程:
时间更新:由k - 1 预测值经过上述的无迹变换过程计算k 时刻预测状态 的均值 和方差矩阵![](http://data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAcFBQYFBAcGBQYIBwcIChELCgkJChUPEAwRGBUaGRgVGBcbHichGx0lHRcYIi4iJSgpKywrGiAvMy8qMicqKyr/2wBDAQcICAoJChQLCxQqHBgcKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKir/wAARCAAZACgDASIAAhEBAxEB/8QAGAABAAMBAAAAAAAAAAAAAAAAAAMFBgf/xAAnEAABAwQCAgIBBQAAAAAAAAABAgMEAAUGERIxEyEiQQcVMlFhYv/EABQBAQAAAAAAAAAAAAAAAAAAAAD/xAAeEQABAgcBAAAAAAAAAAAAAAAAARECIUFRkaHRYf/aAAwDAQACEQMRAD8A7NnOZjFo0ePEMP8AU5y/HFTPeLTO/wDSgCSfoJSCSSOhs1bY3f42R2RmdGcaLmuEhpt1K/A6P3Nkg9g1zTOVTMk/PmL2O2tsvixRnLq+2+4UNlRPFGyEqOwQD191sMCwh3Enb3Onzky7he5ypkjxIKWmyd6QkE/W+z3SGcLrV9K3cCKSsnm58ybClKiMlgSxFL7YkKQXAzzHMpBAKtd62QN/3QEtKUoCghYTZIGWycmjMyBdpaPG8+qW6oLT9J4lXHQ0NDXrVSZMu8JiMCwpeLoeQpzxBv23yAWPn63xKiPXYGzr0bulLJYVczVuk3qJJukm8IluQ2Wy7HQQztSeI2n4++W0q1sgaUP59ReO7y81i3JqO81bVxUNgkIStO9rVy5EkA7QkgDe09gd6qlBRhSlKA//2Q==)
![](http://image.ednchina.com/2015/UAVEKFUKF18.JPG)
测量更新:计算k 时刻测量量zk 及其均值zk 、方差矩阵Pzz,k 及协方差矩阵Pxz,k 。
![](http://image.ednchina.com/2015/UAVEKFUKF19.JPG)
3 仿真模型建立
四轴飞行器状态矢量![](http://image.ednchina.com/2015/xuuUAVEKFUKF04.JPG.pagespeed.ic.oWkeUxDqkW.jpg) , 四轴飞行器的速度和欧拉角由惯性导航单元(InertialMeasurement Unit, IMU)测量及解算。根据四轴飞行器的运动学方程,建立如式(25)所示的状态方程:
![](http://image.ednchina.com/2015/UAVEKFUKF20.JPG)
4 仿真结果
根据四轴飞行器的数学模型及各滤波算法的递推公式,使用Matlab 对系统进行动态仿真。四轴飞行器的空间状态初始值设为:三个方向的位置初值x(0) = 20 m, y(0) =30 m, z(0) =2 m, 三个方向线速度、角速度和角加速初值为零。图2(a) ~ (c)分别为X,Y,Z 三个方向的载体位置曲线。图中的真实值为视觉导航系统给出的位置估计。
![](http://image.ednchina.com/2015/UAVEKFUKF21.JPG)
![](http://image.ednchina.com/2015/UAVEKFUKF22.JPG)
![](http://image.ednchina.com/2015/UAVEKFUKF23.JPG)
图2 X、Y、Z 方向位置曲线
Kalman 滤波在对非线性特征的四轴飞行器的位置滤波上存在明显的滤波偏差,相比之下,EKF 和UKF 有较好的滤波效果,三个方向的估计值误差曲线如图5、图6、图7 所示。从X,Y,Z 三个方向的估计误差曲线图上,更能直观地看出三种滤波算法的滤波结果,滤波估计误差的均方差如表1 所示。
![](http://image.ednchina.com/2015/UAVEKFUKF24.JPG)
![](http://image.ednchina.com/2015/UAVEKFUKF25.JPG)
![](http://image.ednchina.com/2015/UAVEKFUKF26.JPG)
![](http://image.ednchina.com/2015/UAVEKFUKF27.JPG)
图3 X、Y、Z 方向估计误差曲线
Kalman 滤波非线性环境下滤波性能较差, 甚至会出现滤波发散点,而EKF 和UKF 能很好地处理四轴飞行器的非线性特性,从误差均方差上看,UKF 优于EKF,这也是由于EKF是将非线性化的系统线性化造成的偏差,而UKF 是基于样本点的滤波方式,对噪声的适应性更强。
5 结语
本文研究EKF 和UKF 算法在四轴飞行器视觉组合导航中的应用,针对四轴飞行器的非线性特点,传统的Kalman 滤波方法存在一定的的误差,而EKF 和UKF 在融合视觉位置信息时具有明显的优势,且UKF 比EKF 更适应于四轴飞行器的组合导航系统。
|
|