第8卷第3期 智能系统学报 Vol.8 No.3 2013年6月 CAAI Transactions on Intelligent Systems Jum.2013 D0I:10.3969/i.issn.1673-4785.201203012 网络出版地址:http://www.cnki.net/kcms/detail/23.1538.TP.20130310.2131.001.html 一 种基于卡尔曼滤波与模糊算法的 变电站机器人组合导航及控制系统设计 李大伟,贾鹏飞1,李卫国',张学龙1,渠展展 (1.华北电力大学电气与电子工程学院,北京102206:2.中国电力科学研究院电工与新材料技术研究所,北京100192) 摘要:针对不确定因素较多、外在干扰强的变电站工业现场环境,提出了一种新颖的基于卡尔曼能波与模糊控制 的变电站机器人组合导航与控制系统设计方法,实现机器人自主导航定位、自主行为控制,完成智能变电站无人巡 检任务.特别地,设计了一种基于高精度激光雷达的导航子系统,定位的同时兼有电子罗盘功能(即航角测量),有效 克服了变电站强电磁干扰环境下电子罗盘失效等问题现场试验证明,机器人行走轨迹与预设路径十分接近,可以实 现强电磁环境下的巡检要求,且具有设计成本低、精度高(亚米级)、实现方便等一系列优点,便于工业应用. 关键词:变电站:机器人;组合导航:控制系统:模糊控制:激光香达;卡尔曼滤波 中图分类号:TP242.2文献标志码:A文章编号:1673-4785(2013)03-0226-07 中文引用格式:李大伟,贾鹏飞,李卫国,等.一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计[J]智能系 统学报,2013,8(3):226-232. 英文引用格式:LI Dawei,JIA Pengfei,LI Weiguo,etal.A kind of integrated navigation and control system design for substation robot based on the Kalman filtering and fuzzy algorithm[J].CAAI Transactions on Intelligent Systems,2013,8(3):226-232. A kind of integrated navigation and control system design for substation robot based on the Kalman filtering and fuzzy algorithm LI Dawei',JIA Pengfei',LI Weiguo',ZHANG Xuelong',QU Zhanzhan2 (1.School of Electrical Electronic Engineering,North China Electric Power University,Beijing 102206,China;2.Institute of Elec- trical New Materials,China Electric Power Research Institute,Beijing 100192,China) Abstract:This paper aims at examining the substation's industrial environment which has a number of uncertainty factors and strong external disturbance.This paper proposes a novel design method for the integrated navigation and control system which is used in a substation intelligent robot based on the Kalman filtering and fuzzy control theory, for the purpose of fulfilling an autonomous navigation,positioning and behavior control,and finishing the unattend- ed inspection task.In particular,a navigation subsystem based on high-precision laser radar has the functions of po- sitioning,and that of electronic compass (for navigation angle measurement),which effectively overcomes the fail- ure of compass in substation with strong electromagnetic interference.The field test shows that the robot's actual track is close to the assigned track,achieving the independent inspection mission in a strong electromagnetic inter- ference environment.In addition,the design also has the advantages of lower cost,higher precision (sub-meter), and more easy implementation for industrial use. Keywords:transformer substation;robot;integrated navigation;control system;fuzzy control;laser radar;Kalman filte- ring 目前常见的机器人导航系统主要包括黑白线循 及激光雷达导航等.不同的导航系统各自有其优缺 迹导航、磁导航、GPS/DR组合导航、视觉图像导航 点,黑白线导航原理简单易行、稳定可靠但易受雨雪 收稿日期:2012-03-15.网络出版日期:2013-03-10. 天气影响:磁导航成本相对较高且易受强电磁环境 通信作者:李大伟.E-mail:myzone2009_vip@sina.com 干扰:GPS/DR组合导航)定位精度高、灵活方便
第 8 卷第 3 期 智 能 系 统 学 报 Vol.8 №.3 2013 年 6 月 CAAI Transactions on Intelligent Systems Jun. 2013 DOI:10.3969 / j.issn.1673⁃4785.201203012 网络出版地址:http: / / www.cnki.net / kcms/ detail / 23.1538.TP.20130310.2131.001.html 一种基于卡尔曼滤波与模糊算法的 变电站机器人组合导航及控制系统设计 李大伟1 ,贾鹏飞1 ,李卫国1 ,张学龙1 ,渠展展2 (1.华北电力大学 电气与电子工程学院,北京 102206; 2.中国电力科学研究院 电工与新材料技术研究所,北京 100192) 摘 要:针对不确定因素较多、外在干扰强的变电站工业现场环境,提出了一种新颖的基于卡尔曼滤波与模糊控制 的变电站机器人组合导航与控制系统设计方法,实现机器人自主导航定位、自主行为控制,完成智能变电站无人巡 检任务.特别地,设计了一种基于高精度激光雷达的导航子系统,定位的同时兼有电子罗盘功能(即航角测量),有效 克服了变电站强电磁干扰环境下电子罗盘失效等问题.现场试验证明,机器人行走轨迹与预设路径十分接近,可以实 现强电磁环境下的巡检要求,且具有设计成本低、精度高(亚米级)、实现方便等一系列优点,便于工业应用. 关键词:变电站;机器人;组合导航;控制系统;模糊控制;激光雷达;卡尔曼滤波 中图分类号:TP242.2 文献标志码:A 文章编号:1673⁃4785(2013)03⁃0226⁃07 中文引用格式:李大伟,贾鹏飞,李卫国,等.一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计[ J].智能系 统学报, 2013, 8(3): 226⁃232. 英文引用格式:LI Dawei, JIA Pengfei, LI Weiguo, et al. A kind of integrated navigation and control system design for substation robot based on the Kalman filtering and fuzzy algorithm[J]. CAAI Transactions on Intelligent Systems, 2013, 8(3): 226⁃232. A kind of integrated navigation and control system design for substation robot based on the Kalman filtering and fuzzy algorithm LI Dawei 1 , JIA Pengfei 1 , LI Weiguo 1 , ZHANG Xuelong 1 , QU Zhanzhan 2 ( 1. School of Electrical & Electronic Engineering, North China Electric Power University, Beijing 102206, China; 2. Institute of Elec⁃ trical & New Materials, China Electric Power Research Institute, Beijing 100192, China) Abstract:This paper aims at examining the substation’s industrial environment which has a number of uncertainty factors and strong external disturbance. This paper proposes a novel design method for the integrated navigation and control system which is used in a substation intelligent robot based on the Kalman filtering and fuzzy control theory, for the purpose of fulfilling an autonomous navigation, positioning and behavior control, and finishing the unattend⁃ ed inspection task. In particular, a navigation subsystem based on high⁃precision laser radar has the functions of po⁃ sitioning, and that of electronic compass (for navigation angle measurement), which effectively overcomes the fail⁃ ure of compass in substation with strong electromagnetic interference. The field test shows that the robot’ s actual track is close to the assigned track, achieving the independent inspection mission in a strong electromagnetic inter⁃ ference environment. In addition, the design also has the advantages of lower cost, higher precision (sub⁃meter), and more easy implementation for industrial use. Keywords:transformer substation; robot; integrated navigation; control system; fuzzy control; laser radar; Kalman filte⁃ ring 收稿日期:2012⁃03⁃15. 网络出版日期:2013⁃03⁃10. 通信作者:李大伟. E⁃mail: myzone2009_vip@ sina.com. 目前常见的机器人导航系统主要包括黑白线循 迹导航、磁导航、GPS / DR 组合导航、视觉图像导航 及激光雷达导航等.不同的导航系统各自有其优缺 点,黑白线导航原理简单易行、稳定可靠但易受雨雪 天气影响;磁导航成本相对较高且易受强电磁环境 干扰;GPS / DR 组合导航[1] 定位精度高、灵活方便
第3期 李大伟,等:一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计 ·227· 但成本非常高:视觉图像导航基于图像识别与神经 需要指出的是,由于机器人本身存在机械误差 网络技术,但其原理复杂,工业实现难度较大 以及里程计带来的测量误差,DR系统导航精度随 本文涉及的导航系统采用DR/LMS组合导航 时间具有发散性,定位精度有限,故必须辅助其他导 技术,DR导航采用常用的惯性导航方法,原理简单. 航子系统加以修正: 针对LMS激光导航系统的研究,文献[2]从最优路 1.2基于高精度激光雷达LMS的自主导航系统 径角度,给出了一种全局最优路径算法:文献[3]基 本系统采用了SICK公司生产的LMS激光雷达 于激光雷达进行环境建模、避障检测,提出了一种时 全向测距传感器,该测距系统采用双脉冲快速测距 变势场导航算法,但只定性而未定量描述了机器人 技术,以50Hz的扫描频率实现360°、20m范围内障 在视觉场内的坐标和位姿分布).文献[5]基于激 碍物的快速扫描,满足高精度(mm级)测距需求.测 光系统,提出了一种三边与三角测量法,虽然原理简 距系统与上位机Labview软件系统通过以太网接口 单、定位精度高,但仅适用于理想的实验室环境.文 实现相互通信,数据传输及处理速度能够满足实际 献[6]同样在实验室环境下给出了一种机器人速度 工业需求。 跟踪算法,但未能定量地描述位姿坐标问题。 1.2.1全局路径已知条件下的机器人直行段导航 综合目前导航系统的优缺点,本文提出的一种 定位 基于卡尔曼滤波[]的DR/LMS组合导航系统设计 直行段导航定位系统原理如图2所示,其定位 方法,是一种基于多传感器的定位精度优化估计方 与导航计算公式分别见式(2)和(3),式中各变量含 法,通过数据融合、滤波处理,实现了定位坐标定量 义及说明可以参见文献[11] 化管理,具备导航精度高的优点.特别地,本系统采 用的基于多级模糊算法[81的控制结构实现了机器 head 人运动系统的多级化闭环控制,优化机器人的行为 Left Right 轨迹的同时,增强了其对复杂工业现场环境的适应 (x1Y、 1△a 能力. 测量 1机器人导航系统原理及应用 (x,y)△ robot 本导航系统采用DR/LMS组合式设计,各传感 器相互融合,通过卡尔曼最优估计原理进行算法设 图2直行段导航定位原理 计,提高了导航的精确性与稳定性。 Fig.2 The principle diagram of navigation and positio- 1.1DR惯性导航系统 ning when going straight 图1是DR导航定位系统航位推算的原理图, 已知起始点A的定位坐标(x,y)、导航测量角 (x-x1)2+(y-y1)2=, 及光电位移传感器测量里程△S,可以推算出下一 (2) 时刻在B点的坐标估计值(x+1,y+1).式(1)给出了 (x-x)2+(y-为2)2= 机器人在X、Y轴方向的行驶位移,即 d2=√+ri-2rr2cos△a2; |xk+i=xg+△S4sin0, (1) +d品- yk+1=ye+△S4cos0. cos B= 2r2dn k1=(x1-x2)/(y1-y2): C(x,y) k2=0: (3) B(x,y) [T/2,y1=y2, Y= k1-k2 arctan ,y1≠y2; (xy) 1+kk2 X △am=T-a; 图1DR系统导航定位原理 △a,=T-B-Y; Fig.1 The principle of DR navigation and positioning sys- tem 0g=△am-△a
但成本非常高;视觉图像导航基于图像识别与神经 网络技术,但其原理复杂,工业实现难度较大. 本文涉及的导航系统采用 DR/ LMS 组合导航 技术,DR 导航采用常用的惯性导航方法,原理简单. 针对 LMS 激光导航系统的研究,文献[2]从最优路 径角度,给出了一种全局最优路径算法;文献[3]基 于激光雷达进行环境建模、避障检测,提出了一种时 变势场导航算法,但只定性而未定量描述了机器人 在视觉场内的坐标和位姿分布[4] .文献[5] 基于激 光系统,提出了一种三边与三角测量法,虽然原理简 单、定位精度高,但仅适用于理想的实验室环境.文 献[6]同样在实验室环境下给出了一种机器人速度 跟踪算法,但未能定量地描述位姿坐标问题. 综合目前导航系统的优缺点,本文提出的一种 基于卡尔曼滤波[7] 的 DR/ LMS 组合导航系统设计 方法,是一种基于多传感器的定位精度优化估计方 法,通过数据融合、滤波处理,实现了定位坐标定量 化管理,具备导航精度高的优点.特别地,本系统采 用的基于多级模糊算法[8⁃10]的控制结构实现了机器 人运动系统的多级化闭环控制,优化机器人的行为 轨迹的同时,增强了其对复杂工业现场环境的适应 能力. 1 机器人导航系统原理及应用 本导航系统采用 DR/ LMS 组合式设计,各传感 器相互融合,通过卡尔曼最优估计原理进行算法设 计,提高了导航的精确性与稳定性. 1.1 DR 惯性导航系统 图 1 是 DR 导航定位系统航位推算的原理图, 已知起始点 A 的定位坐标( xk,yk )、导航测量角 θk 及光电位移传感器测量里程 ΔSk,可以推算出下一 时刻在 B 点的坐标估计值(xk+1 ,yk+1 ).式(1)给出了 机器人在 X、Y 轴方向的行驶位移,即 xk+1 = xk + ΔSk sin θk, yk+1 = yk + ΔSk cos θk . { (1) 图 1 DR 系统导航定位原理 Fig.1 The principle of DR navigation and positioning sys⁃ tem 需要指出的是,由于机器人本身存在机械误差 以及里程计带来的测量误差,DR 系统导航精度随 时间具有发散性,定位精度有限,故必须辅助其他导 航子系统加以修正. 1.2 基于高精度激光雷达 LMS 的自主导航系统 本系统采用了 SICK 公司生产的 LMS 激光雷达 全向测距传感器,该测距系统采用双脉冲快速测距 技术,以50 Hz的扫描频率实现 360°、20 m范围内障 碍物的快速扫描,满足高精度(mm 级)测距需求.测 距系统与上位机 Labview 软件系统通过以太网接口 实现相互通信,数据传输及处理速度能够满足实际 工业需求. 1.2.1 全局路径已知条件下的机器人直行段导航 定位 直行段导航定位系统原理如图 2 所示,其定位 与导航计算公式分别见式(2)和(3),式中各变量含 义及说明可以参见文献[11]. 图 2 直行段导航定位原理 Fig.2 The principle diagram of navigation and positio⁃ ning when going straight (x - x1 ) 2 + (y - y1 ) 2 = r 2 1 , (x - x2 ) 2 + (y - y2 ) 2 = r 2 2 . { (2) d ~ 12 = r 2 1 + r 2 2 - 2r1 r2 cos Δα12 ; cos β = r 2 2 + d ~ 2 12 - r 2 1 2r2 d ~ 12 ; k1 = (x1 - x2 ) / (y1 - y2 ); k2 = 0; γ = π/ 2,y1 = y2 , arctan k1 - k2 1 + k1 k2 ,y1 ≠ y2 ; ì î í ï ï ïï Δαm = π - α; Δαt = π - β - γ; θk = Δαm - Δαt . ì î í ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï (3) 第 3 期 李大伟,等:一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计 ·227·
·228· 智能系统学报 第8卷 1.2.2全局路径已知条件下的机器人转弯段导航 理坐标分别为(x1,y1)和(x2,y2),1、2分别表示机 定位 器人距离对应参考对象的距离值,由此可以方便地 机器人定位坐标的求解同式(2).与直行段不同 解算出机器人当前的定位坐标(x,y).再根据定位坐 的是,转弯段预设行走路径不再是直线方向,而是围 标(x,y)及转弯条件下计算得到的圆心坐标(xo, 绕某一圆心的圆弧运动.如图3所示,机器人以0, yo,)可以计算出机器人距离圆心O1的距离值R 为圆心,半径为R(理论值)做圆周运动,其中圆心 及机器人朝向与半径方向的夹角值A.将这2个计算 0,的地理坐标(xo,yo,)与半径R的大小已事先 量分别与转弯段的理想半径值R和理想夹角值 由系统给定.假设由式(2)求得机器人的定位坐标为 (90°)作定量比较,得到机器人自主转弯的2个控 (x,y),则机器人自主导航的控制量△R与△8的求 制量:半径偏移量△R与航角偏移量△0,并将这2 解方程如式(4): 个控制量作为模糊控制器的2个控制输入参数,实 R=√x-o,)2+(y-yo,)尸, 现控制策略下的机器人导航定位闭环反馈调节. △R=Rnea-Ral, 本定位子系统的主要误差来源有:1)实测道路 两侧参考目标设备地理坐标带来的绝对误差:2)激 6=(T-w1)+arccos[(R2+r斤- 光雷达实测目标距离值时多点滤波等效带来的绝对 (x-xn,)2+(y1-yo,)2)/(2×R×r1)], 误差.鉴于上述两子系统的优缺点,为了提高系统定 △0=m/2-8. 位精度、保证系统运行的可靠性,引入基于卡尔曼滤 (4) 波的组合导航技术则显得尤为必要 式中:R表示机器人定位点与当前圆心O,的测量 L.3基于卡尔曼滤波技术的DR/LMS组合导航系 距离,π为机器人前进方向上激光雷达扫描点的对 统设计及应用 应角度值,w,为激光雷达在转弯内侧参考标记1 针对DR惯性导航子系统与LMS高精度激光 (x1,y:)处的扫描角度值,0表示机器人前进方向与 雷达视觉导航子系统各自的优缺点,为了进一步提 半径方向的夹角值,π/2则表示理想圆心运动条件 高系统的导航定位精度,改善系统的稳定性与可靠 下,机器人前进方向即切线方向与半径方向的垂直 性,可以以随机线性离散卡尔曼滤波为理论基础,将 夹角. DR与LMS两导航子系统优势互补,以增强系统的 (x2y) 鲁棒性与抗干扰性。 ◆6 4① ◇2 基于卡尔曼滤波模型,对DR/LMIS组合导航系 1rs 统的X、Y定位坐标进行卡尔曼滤波,可以得到如下 3④ Robot 方程 ◇5 (x.v) R 1)状态方程: O(x。y。) ①A +AStk-isin 01+I(-)W (1) y:=y1+ASL-icos -1+I(-W-Dy ⊙ 2)观测方程: 2,=g+V: 图3转弯段导航定位原理 Fig.3 The principle diagram of navigation and positio- (n =ye+Var ning when turning control 式中:xy为机器人在X轴、Y轴方向的状态位置 图3中序列1~6分别表示机器人视觉区域内 坐标,0-1为状态控制量(航向计算角),△S.-1为各 的参考对象,并且各个参考对象的地理坐标已知,机 采样阶段的转移控制系数,Wx、W-,为DR推 器人可以根据激光雷达,通过数据采样、过滤及滤波 算过程噪声序列,「k-)、Tk-w为DR推算过程 算法识别当前环境下的参考距离值,并建立参考对 噪声输入系数,之,之,为激光测距计算得到的机器人 象与参考距离值之间的一一映射关系,准确计算出 定位点的x、y坐标值,V、V为LMS测量系统等效 机器人当前定位坐标及航向角偏差.与直行段不同 噪声序列 的是,转弯段可以看作是机器人绕某一固定圆心点 系统定位坐标初始状态x。y。由系统给定,滤 的圆周运动,而不同的转弯段可以更新参考圆心坐 波误差初始值P。已知,系统过程噪声(估计噪声) 标来实现自主变轨运动.图3中参考对象1、2的地 方差Q.=E(E.W)和系统观测噪声(LMS测量噪
1.2.2 全局路径已知条件下的机器人转弯段导航 定位 机器人定位坐标的求解同式(2).与直行段不同 的是,转弯段预设行走路径不再是直线方向,而是围 绕某一圆心的圆弧运动.如图 3 所示,机器人以 O1 为圆心,半径为 Rmea(理论值)做圆周运动,其中圆心 O1 的地理坐标( xO1 ,yO1 )与半径 Rmea的大小已事先 由系统给定.假设由式(2)求得机器人的定位坐标为 (x,y),则机器人自主导航的控制量 ΔR 与 Δθ 的求 解方程如式(4): Rmea = (x - xO1 ) 2 + (y - yO1 ) 2 , ΔR = Rmea - Rcal, θ = (π - ω1 ) + arccos[(R 2 mea + r 2 1 - ((x1 - xO1 ) 2 + (y1 - yO1 ) 2 )) / (2 × Rcal × r1 )], Δθ = π/ 2 - θ. ì î í ï ï ï ï ï ï ï ï (4) 式中:Rcal表示机器人定位点与当前圆心 O1 的测量 距离,π 为机器人前进方向上激光雷达扫描点的对 应角度值,ω1为激光雷达在转弯内侧参考标记 1 (x1 ,y1 )处的扫描角度值,θ 表示机器人前进方向与 半径方向的夹角值,π/ 2 则表示理想圆心运动条件 下,机器人前进方向即切线方向与半径方向的垂直 夹角. 图 3 转弯段导航定位原理 Fig.3 The principle diagram of navigation and positio⁃ ning when turning control 图 3 中序列 1 ~ 6 分别表示机器人视觉区域内 的参考对象,并且各个参考对象的地理坐标已知,机 器人可以根据激光雷达,通过数据采样、过滤及滤波 算法识别当前环境下的参考距离值,并建立参考对 象与参考距离值之间的一一映射关系,准确计算出 机器人当前定位坐标及航向角偏差.与直行段不同 的是,转弯段可以看作是机器人绕某一固定圆心点 的圆周运动,而不同的转弯段可以更新参考圆心坐 标来实现自主变轨运动.图 3 中参考对象 1、2 的地 理坐标分别为(x1 ,y1 )和( x2 ,y2 ),r1 、r2 分别表示机 器人距离对应参考对象的距离值,由此可以方便地 解算出机器人当前的定位坐标(x,y).再根据定位坐 标(x,y) 及转弯条件下计算得到的圆心坐标( xO1 , yO1 )可以计算出机器人距离圆心 O1 的距离值 Rmea 及机器人朝向与半径方向的夹角值 θ.将这 2 个计算 量分别与转弯段的理想半径值 R 和理想夹角值 (90°)作定量比较,得到机器人自主转弯的 2 个控 制量:半径偏移量 ΔR 与航角偏移量 Δθ,并将这 2 个控制量作为模糊控制器的 2 个控制输入参数,实 现控制策略下的机器人导航定位闭环反馈调节. 本定位子系统的主要误差来源有:1)实测道路 两侧参考目标设备地理坐标带来的绝对误差;2)激 光雷达实测目标距离值时多点滤波等效带来的绝对 误差.鉴于上述两子系统的优缺点,为了提高系统定 位精度、保证系统运行的可靠性,引入基于卡尔曼滤 波的组合导航技术则显得尤为必要. 1.3 基于卡尔曼滤波技术的 DR/ LMS 组合导航系 统设计及应用 针对 DR 惯性导航子系统与 LMS 高精度激光 雷达视觉导航子系统各自的优缺点,为了进一步提 高系统的导航定位精度,改善系统的稳定性与可靠 性,可以以随机线性离散卡尔曼滤波为理论基础,将 DR 与 LMS 两导航子系统优势互补,以增强系统的 鲁棒性与抗干扰性. 基于卡尔曼滤波模型,对 DR/ LMS 组合导航系 统的 X、Y 定位坐标进行卡尔曼滤波,可以得到如下 方程. 1)状态方程: xk = xk-1 + ΔSk,k-1 sin θk-1 + Γ(k,k-1)xW(k-1)x, yk = yk-1 + ΔSk,k-1 cos θk-1 + Γ(k,k-1)yW(k-1)y . { 2)观测方程: zx k = xk + Vkx, zy k = yk + Vky . { 式中:xk、yk 为机器人在 X 轴、Y 轴方向的状态位置 坐标,θk-1为状态控制量(航向计算角),ΔSk,k-1为各 采样阶段的转移控制系数,W(k-1)x、W(k-1)y为 DR 推 算过程噪声序列,Γ(k,k-1) x、Γ(k,k-1)y为 DR 推算过程 噪声输入系数,zx k 、zy k 为激光测距计算得到的机器人 定位点的 x、y 坐标值,Vkx、Vky为 LMS 测量系统等效 噪声序列. 系统定位坐标初始状态 x0 、y0 由系统给定,滤 波误差初始值 P0 已知,系统过程噪声(估计噪声) 方差Qk = E(EkW T k ) 和系统观测噪声( LMS 测量噪 ·228· 智 能 系 统 学 报 第 8 卷
第3期 李大伟,等:一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计 ·229· 声)方差R=E(V)已确定,则得到如下滤波求解 而对于机器人自主行走控制器的另外一个控制 过程. 量一导航角,这里不对其作卡尔曼滤波处理,因为 1)状态一步预测: 其算法的优越性使得导航角的计算精度不受测距仪 ]- sin 0- 多点扫描带来的误差影响,导航精度高,仅需加以常 AS-1 用的限幅滤波、均值滤波等算法即可满足控制系统 2)状态估计: 的要求 1「(1-K. 2 基于模糊算法的机器人行走控制系 G. 0 统设计及应用 3)滤波增益: 通过提取上述机器人导航控制参数,将其作为 P(k,k-10 「K P(A-1)s +Ru 输入量,通过模糊算法设计,实现机器人的行走运动 控制,该算法极大地增强了系统对复杂工业现场环 P(k.k-1)3 境的适应能力. P(k.k-1)+Riy- 图4给出了机器人运动系统的模糊控制流程 4)步预测误差: 图新方法反映了机器人以导航反馈量,即距离偏移 P(k.t-1)x P(-D T-DQ- LP- , 量与航角偏移量为控制器输入参数,通过模糊推理、 -- 解模糊集,输出电机转速控制量,并将输出响应后的 5)滤波估计误差: 下一次导航定位数据反馈到原系统,通过给定值与 0 反馈值的比较,实现了多级模糊控制策略下的快速 响应特性与闭环控制, [P(6,k-1 「KRa P (k.k-1)y 数据库 规则库 给定值 距离偏差 左轮转速 模糊推理 右轮转速 化 糊 角度偏差 转向角度 机器人定位坐标以及朝向角 图4机器人模糊控制流程 Fig.4 The flow chart of robot fuzzy control 控制器为两输入三输出的模糊控制器,将真实 位u1和2的论域为[0,4],划分为5个等级;机器 的确定量输入转换为一个模糊矢量.控制器实时接 人转向角度u3的论域为[220,500],划分为9个等 收机器人的定位坐标和航向角,并将其与给定的目 级.以图5和图6为例,分别给出了基于Labview平 标值比较,计算当前距离偏差量和航角偏差量.2个 台的s和u3的隶属度函数图形. 偏差量作为控制器输入先进行模糊化处理,得到输 1.0m -NBI 入的模糊矢量后,根据预先设定的数据库和规则库 NBMI 0.8 NSI 进行模糊推理,推理后得到的结果为模糊矢量,进行 世0.6 NSMI ZI 解模糊后作为控制量输出,输出量分别为左电机转 0.4 PSMI 速、右电机转速以及转向电机调节角度.2个输入 0.2 PSI PBMI 中,机器人中心与道路中心线的距离偏差s的论域 -1.0-0.6 -0.20.20.6 1.0 为[-1.5,1.5],划分为9个等级:机器人朝向与道路 方向的角度偏差△p的论域为[-15,15],划分为5 图5距离偏移量s的隶属度函数 个等级.3个输出中,机器人左右2个电机的速度档 Fig.5 Membership functions of distance offset s
声)方差Rk =E(VkV T k )已确定,则得到如下滤波求解 过程.. 1)状态一步预测: x ^ k,k-1 y ^ k,k-1 é ë ê ê ù û ú ú = x ^ k-1 y ^ k-1 é ë ê ê ù û ú ú + ΔSk,k-1 sin θk-1 cos θk-1 é ë ê ê ù û ú ú . 2)状态估计: x ^ k y ^ k é ë ê ê ù û ú ú = (1 - Kkx) 0 0 (1 - Kky) é ë ê ê ù û ú ú x ^ k,k-1 y ^ k,k-1 é ë ê ê ù û ú ú + Kkx zx k Kky zy k é ë ê ê ù û ú ú . 3)滤波增益: Kkx Kky é ë ê ê ù û ú ú = P(k,k-1)x P(k,k-1)x + Rkx P(k,k-1)y P(k,k-1)y + Rky é ë ê ê ê ê ê ù û ú ú ú ú ú . 4)步预测误差: P(k,k-1)x P(k,k-1)y é ë ê ê ù û ú ú = P(k-1)x P(k-1)y é ë ê ê ù û ú ú + Γ 2 (k,k-1)xQ(k-1)x Γ 2 (k,k-1)yQ(k-1)y é ë ê ê ù û ú ú . 5)滤波估计误差: Pkx Pky é ë ê ê ù û ú ú = (1 - Kkx) 2 0 0 (1 - Kky) 2 é ë ê ê ù û ú ú · P(k,k-1)x P(k,k-1)y é ë ê ê ù û ú ú + K 2 kxRkx K 2 kyRky é ë ê ê ù û ú ú . 而对于机器人自主行走控制器的另外一个控制 量———导航角,这里不对其作卡尔曼滤波处理,因为 其算法的优越性使得导航角的计算精度不受测距仪 多点扫描带来的误差影响,导航精度高,仅需加以常 用的限幅滤波、均值滤波等算法即可满足控制系统 的要求. 2 基于模糊算法的机器人行走控制系 统设计及应用 通过提取上述机器人导航控制参数,将其作为 输入量,通过模糊算法设计,实现机器人的行走运动 控制,该算法极大地增强了系统对复杂工业现场环 境的适应能力. 图 4 给出了机器人运动系统的模糊控制流程 图.新方法反映了机器人以导航反馈量,即距离偏移 量与航角偏移量为控制器输入参数,通过模糊推理、 解模糊集,输出电机转速控制量,并将输出响应后的 下一次导航定位数据反馈到原系统,通过给定值与 反馈值的比较,实现了多级模糊控制策略下的快速 响应特性与闭环控制. 图 4 机器人模糊控制流程 Fig.4 The flow chart of robot fuzzy control 控制器为两输入三输出的模糊控制器,将真实 的确定量输入转换为一个模糊矢量.控制器实时接 收机器人的定位坐标和航向角,并将其与给定的目 标值比较,计算当前距离偏差量和航角偏差量.2 个 偏差量作为控制器输入先进行模糊化处理,得到输 入的模糊矢量后,根据预先设定的数据库和规则库 进行模糊推理,推理后得到的结果为模糊矢量,进行 解模糊后作为控制量输出,输出量分别为左电机转 速、右电机转速以及转向电机调节角度. 2 个输入 中,机器人中心与道路中心线的距离偏差 s 的论域 为[-1.5,1.5],划分为 9 个等级;机器人朝向与道路 方向的角度偏差 Δφ 的论域为[ -15,15],划分为 5 个等级.3 个输出中,机器人左右 2 个电机的速度档 位 u1 和 u2 的论域为[0,4],划分为 5 个等级;机器 人转向角度 u3 的论域为[220,500],划分为 9 个等 级.以图 5 和图 6 为例,分别给出了基于 Labview 平 台的 s 和 u3 的隶属度函数图形. 图 5 距离偏移量 s 的隶属度函数 Fig.5 Membership functions of distance offset s 第 3 期 李大伟,等:一种基于卡尔曼滤波与模糊算法的变电站机器人组合导航及控制系统设计 ·229·
·230· 智能系统学报 第8卷 为了得到合理的输出控制量,要求模糊计算方 图8定量地给出了不同滤波策略下本组合导航 法能够很好地表达输出隶属度函数的计算结果采 系统X坐标定位值的精度变化曲线,给定道路中心 用加权平均法,将论域中的每个元素x:(i=1,2,…, 线的直线方程约束条件为X=3.3m.从图8中可以 n)作为待判决输出模糊集合的隶属度u:的加权系 明显地看出,与X值理想定位变化趋势相比,经过 数,求得乘积山,再计算该乘积和三x叫,对于隶属 卡尔曼滤波后的定位波形具有更高的精度与可信 度经现场试验,引入卡尔曼滤波后,系统的定位精 度和的平均值x。,将其作为求得的判决结果最后, 度更加逼近于机器人的实际行为轨迹,其定位精度 将输出量化因子乘以xo得到所需要的控制量. NB5 可以达到10cm左右,并且经过多级参数整定后的模 NBM5 糊控制系统控制策略更加科学严谨.整个系统的导 0.8 NS5 NSM5 航控制频率在5Hz左右,其响应速度完全可以满足 0.6 Z5 现场应用需求. 0.4 PSM5 PS5 4.2m PBM5 10 4.0 限幅滤波X坐标曲线 2.2 26 3.03.43.84.24.65.0 理想X坐标曲线 43 3.8 图6转向角度输出山,的隶属度函数 E3.6 Fig.6 Membership function of steering angle output u 3.4 Kalman 3.2 滤波曲线 3.0 3 变电站机器人现场导航轨迹波形分 滤波前X坐标曲线 2.80 2040 6080100120140160 析及验证 定位个数n 本节对基于卡尔曼滤波与模糊控制的DR/LMS 图8定位X坐标滤波波形抽样分析 组合导航定位系统的机器人现场自主巡检过程进行 Fig.8 Sampling analysis on filtering wave of positio- 了轨迹跟踪与数据反馈.通过对全局路径下一个周 ning X coordinate 期路径定位数据的采集及波形分析,定量地验证了 图9定量地给出了机器人行为过程中本组合导 本导航及控制系统的动态稳定性和高精度性, 航系统定位Y坐标的精度变化曲线.由于复杂工业 图7给出了500kV变电站现场路径中的一段周 环境对采样数据带来的噪声干扰,定位Y坐标在某 期性路径,图中U形段为已知巡检路段,道路两侧 些定位点会出现幅度较大的突变;而经限幅滤波、均 的方框标识为视觉域内的参考点或自由空间点,其 值滤波等常用经典滤波方法有效地消除了外界噪声 中参考点的地理坐标已知,坐标系如图7所示理想 干扰.经过卡尔曼滤波后的Y坐标变化曲线的线性 的机器人行走轨迹应该是围绕中心路径(图7中虚 度与经典滤波波形曲线类似,区别仅仅在于一定程 线)的U形运动. 度上提高了系统定位精度与可靠性.而实际导航系 ◇ 0 统的控制量主要取决于X坐标的偏移量和导航角, 为了简化控制策略及程序的复杂性、提高响应速度, 对于Y坐标采用常用经典滤波的定位算法完全可 以满足现场要求 ◇ 30 经典滤波后的Y坐标 ◇ 中 号2 ◇ ◇ 由 18 滤波前的Y坐标测量值 (其余部分与经典滤波曲线重合) 0 2040 6080100120140160 测量个数n 图7自主巡检机器人现场理想轨迹模拟 图9定位Y坐标滤波波形抽样分析 Fig.7 Field track simulation of the independent inspec- Fig.9 Sampling analysis on filtering wave of positio- tion robot ning Y coordinate
为了得到合理的输出控制量,要求模糊计算方 法能够很好地表达输出隶属度函数的计算结果.采 用加权平均法,将论域中的每个元素 xi( i = 1,2,…, n)作为待判决输出模糊集合的隶属度 ui 的加权系 数,求得乘积 xiui,再计算该乘积和∑ n i = 1 xiui 对于隶属 度和的平均值 x0 ,将其作为求得的判决结果.最后, 将输出量化因子乘以 x0 得到所需要的控制量. 图 6 转向角度输出 u3 的隶属度函数 Fig.6 Membership function of steering angle output u3 3 变电站机器人现场导航轨迹波形分 析及验证 本节对基于卡尔曼滤波与模糊控制的 DR/ LMS 组合导航定位系统的机器人现场自主巡检过程进行 了轨迹跟踪与数据反馈.通过对全局路径下一个周 期路径定位数据的采集及波形分析,定量地验证了 本导航及控制系统的动态稳定性和高精度性. 图 7 给出了500 kV变电站现场路径中的一段周 期性路径,图中 U 形段为已知巡检路段,道路两侧 的方框标识为视觉域内的参考点或自由空间点,其 中参考点的地理坐标已知,坐标系如图 7 所示.理想 的机器人行走轨迹应该是围绕中心路径(图 7 中虚 线)的 U 形运动. 图 7 自主巡检机器人现场理想轨迹模拟 Fig.7 Field track simulation of the independent inspec⁃ tion robot 图 8 定量地给出了不同滤波策略下本组合导航 系统 X 坐标定位值的精度变化曲线,给定道路中心 线的直线方程约束条件为 X = 3.3 m.从图 8 中可以 明显地看出,与 X 值理想定位变化趋势相比,经过 卡尔曼滤波后的定位波形具有更高的精度与可信 度.经现场试验,引入卡尔曼滤波后,系统的定位精 度更加逼近于机器人的实际行为轨迹,其定位精度 可以达到10 cm左右,并且经过多级参数整定后的模 糊控制系统控制策略更加科学严谨.整个系统的导 航控制频率在5 Hz左右,其响应速度完全可以满足 现场应用需求. 图 8 定位 X 坐标滤波波形抽样分析 Fig.8 Sampling analysis on filtering wave of positio⁃ ning X coordinate 图 9 定量地给出了机器人行为过程中本组合导 航系统定位 Y 坐标的精度变化曲线.由于复杂工业 环境对采样数据带来的噪声干扰,定位 Y 坐标在某 些定位点会出现幅度较大的突变;而经限幅滤波、均 值滤波等常用经典滤波方法有效地消除了外界噪声 干扰.经过卡尔曼滤波后的 Y 坐标变化曲线的线性 度与经典滤波波形曲线类似,区别仅仅在于一定程 度上提高了系统定位精度与可靠性.而实际导航系 统的控制量主要取决于 X 坐标的偏移量和导航角, 为了简化控制策略及程序的复杂性、提高响应速度, 对于 Y 坐标采用常用经典滤波的定位算法完全可 以满足现场要求. 图 9 定位 Y 坐标滤波波形抽样分析 Fig.9 Sampling analysis on filtering wave of positio⁃ ning Y coordinate ·230· 智 能 系 统 学 报 第 8 卷