基于聚类和K近邻算法的井下人员定位算法.pdf
扫码移动阅读 第45卷 第4期 2 0 1 9年4月 工矿自动化 In d u s t r y a n d M i n e A u t o m a t i o n Vo l . 4 5N o . 4 A p r . 2 0 1 9 文章编号1 6 7 1-25 1 X(2 0 1 9)0 4-00 4 3-07 D O I1 0. 1 3 2 7 2/ j . i s s n. 1 6 7 1-25 1 x. 2 0 1 8 1 1 0 0 7 2 基于聚类和 K近邻算法的井下人员定位算法 莫树培1, 2, 唐琎 2, 汪郁 1, 赖普坚 2, 金礼模 1 ( 1. 贵州工业职业技术学院 图书与信息中心,贵州 贵阳 5 5 1 4 0 0; 2. 中南大学 信息科学与工程学院,湖南 长沙 4 1 0 0 8 3) 摘要 针对现有基于指纹模的井下定位算法存在的计算量大、 实时性低、 定位精度较低的问题, 提出了基 于聚类和 K近邻算法的井下人员定位算法。用二分k-me a n s聚类算法对采集的 R S S I数据进行分类, 建立 离线指纹数据库; 无线移动终端和动态修正器实时采集 R S S I值, 分别存储到在线定位数据库和动态修正数 据库; 根据待测点和动态修正器的离线数据和实时数据, 采用软硬件动态修正加权 K 近邻算法计算权重值, 结合离线指纹数据库中待测点的物理位置信息估算其实时位置。实验分析结果表明, 所提定位算法的最小 标准误差为0. 4 6m, 最大标准误差为3. 2 6m, 平均误差为1. 6 2m。对比分析结果表明, 与未进行聚类分析 的算法相比, 本文算法的精度更高, 实时性更好; 与未动态修正权重值的算法相比, 本文算法的运算时间略有 增加, 但定位精度提高了3 7. 2 1%。 关键词 井下人员定位;指纹定位;二分k-me a n s聚类算法;软硬件动态修正加权K近邻算法;动态 修正 中图分类号T D 6 5 5. 3 文献标志码A 收稿日期2 0 1 8-11-30; 修回日期2 0 1 9-03-10; 责任编辑 胡娴。 基金项目 贵州省科技厅项目( 黔科合 LH 字〔 2 0 1 6〕7 0 6 9) ; 贵州工业职业技术学院校级科研课题(2 0 1 8 0 0 9) 。 作者简介 莫树培(1 9 8 2-) , 男, 布依族, 贵州都匀人, 副教授, 硕士, 主要研究方向为无线网络传感器、 模式识别与人工智能, E-ma i l4 1 8 5 1 3 4 8 @q q . c o m。 引用格式 莫树培, 唐琎, 汪郁, 等.基于聚类和K近邻算法的井下人员定位算法[J]. 工矿自动化,2 0 1 9,4 5(4) 4 3-48. MO S h u p e i,TANG J i n,WANG Y u,e t a l . U n d e r g r o u n d p e r s o n n e l p o s i t i o n i n g a l g o r i t h m b a s e d o n c l u s t e r i n g a n d K-ne a r e s t n e i g h b o r a l g o r i t h m[J]. I n d u s t r y a n d M i n e A u t o m a t i o n,2 0 1 9,4 5(4) 4 3-48. Un d e r g r o u n d p e r s o n n e l p o s i t i o n i n g a l g o r i t h m b a s e d o n c l u s t e r i n g a n d K-n e a r e s t n e i g h b o r a l g o r i t h m MO S h u p e i 1,2, TANG J i n 2, WANG Y u 1, LA I P u j i a n 2, JI N L i m o 1 ( 1 . B o o k a n d I n f o r m a t i o n C e n t e r,G u i z h o u I n d u s t r y P o l y t e c h n i c C o l l e g e,G u i y a n g 5 5 1 4 0 0,C h i n a; 2. S c h o o l o f I n f o r m a t i o n S c i e n c e a n d E n g i n e e r i n g,C e n t r a l S o u t h U n i v e r s i t y,C h a n g s h a 4 1 0 0 8 3,C h i n a) A b s t r a c tI n v i e w o f p r o b l e m s o f l a r g e a m o u n t o f c a l c u l a t i o n,l o w r e a l -ti m e p e r f o r m a n c e a n d l o w p o s i t i o n i n g a c c u r a c y o f e x i s t i n g f i n g e r p r i n t -ba s e d u n d e r g r o u n d p o s i t i o n i n g a l g o r i t h m,u n d e r g r o u n d p e r s o n n e l p o s i t i o n i n g a l g o r i t h m b a s e d o n c l u s t e r i n g a n d K-ne a r e s t n e i g h b o r a l g o r i t h m w a s p r o p o s e d . B i s e c t i n g k-me a n s c l u s t e r i n g a l g o r i t h m i s u s e d t o c l a s s i f y c o l l e c t e d R S S I d a t a t o e s t a b l i s h a n o f f l i n e f i n g e r p r i n t d a t a b a s e .R e a l t i m e R S S I v a l u e s a r e c o l l e c t e d b y w i r e l e s s m o b i l e t e r m i n a l a n d d y n a m i c c o r r e c t o r a n d s t o r e d i n o n l i n e p o s i t i o n i n g d a t a b a s e a n d d y n a m i c c o r r e c t i o n d a t a b a s e r e s p e c t i v e l y . A c c o r d i n g t o o f f l i n e d a t a a n d r e a l-ti m e d a t a,w e i g h t v a l u e i s c a l c u l a t e d u s i n g s o f t w a r e a n d h a r d w a r e d y n a m i c c o r r e c t i o n w e i g h t e d K-ne a r e s t n e i g h b o r a l g o r i t h m,a n d r e a l-ti m e p o s i t i o n i s e s t i m a t e d b y c o m b i n i n g t h e p h y s i c a l l o c a t i o n i n f o r m a t i o n o f t h e p o i n t t o b e m e a s u r e d i n t h e o f f l i n e f i n g e r p r i n t d a t a b a s e . T h e e x a m p l e a n a l y s i s r e s u l t s s h o w t h a t t h e m i n i m u m s t a n d a r d e r r o r o f t h e p r o p o s e d p o s i t i o n i n g a l g o r i t h m i s 0. 4 6m, t h e m a x i m u m s t a n d a r d e r r o r i s 3. 2 6 m,a n d t h e a v e r a g e e r r o r i s 1. 6 2 m.T h e r e s u l t s o f c o m p a r a t i v e a n a l y s i s s h o w t h a t t h e p r o p o s e d a l g o r i t h m h a s h i g h e r p r e c i s i o n a n d b e t t e r r e a l-ti m e p e r f o r m a n c e t h a n t h e a l g o r i t h m w i t h o u t c l u s t e r i n g a n a l y s i s .C o m p a r e d w i t h t h e a l g o r i t h m w i t h o u t d y n a m i c c o r r e c t i o n o f we i g h t s,t h e c o m p u t a t i o n t i m e o f t h e p r o p o s e d a l g o r i t h m i s s l i g h t l y i n c r e a s e d,b u t t h e p o s i t i o n i n g a c c u r a c y i s i n c r e a s e d b y 3 7. 2 1%. K e y w o r d sU n d e r g r o u n d p e r s o n n e l p o s i t i o n i n g;f i n g e r p r i n t p o s i t i o n i n g;b i s e c t i n g k-me a n s c l u s t e r i n g a l g o r i t h m;s o f t w a r e a n d h a r d w a r e d y n a m i c c o r r e c t i o n w e i g h t e d K-ne a r e s t n e i g h b o r a l g o r i t h m; d y n a m i c c o r r e c t i o n 0 引言 在井下无线传感器网络中, 可通过网络 MA C 地址、 接收信号强度信息(R S S I) 和传输时间等进行 人员定位[ 1]。在井下无线定位方法中, 基于 AOA ( 到达角度) 的方法[ 2]计算复杂度高, 基于 T OA( 到 达时间) [3]、 TD OA( 到达时间差) [4]的方法对设备有 严格的时间同步要求, 而基于 R S S I的方法获取定 位信号强度方便, 功耗低、 成本低, 是一种较好的井 下定位方式[ 5-7]。 基于 R S S I的井下定位方法有三边定位、 质心 定位和指纹定位等。尚俊龙等[ 8]提出了井下加权三 边定位算法, 用加权三边测算出待测点位置, 定位精 度为3m; 赵彤等[ 9]提出了井下节点合作加权质心 定位算法, 利用改进的小区域加权质心定位算法对 待测点进行定位, 定位精度为1. 5m, 但其实验场地 为楼道, 与井下环境差距较大; 王永星等[ 10]提出了 基于信息传输模型的井下目标定位算法, 利用最小 二乘法建立指纹模型, 定位精度为2. 7 3m。董建平 等[ 11]提出了基于指纹模的井下定位算法, 利用采样 点信息形成指纹数据库, 并利用K近邻算法进行在 线匹配定位, 但该方法存在以下问题 指纹数据库未 进行数据分类, 定位时必须计算出全部参考点的欧 氏距离, 才能得到K个最近的参考点, 计算量较大, 实时性不好; 在K近邻算法中采用K个点的平均 值计算定位位置, 未考虑距离较远和较近的参考点 的权重值问题, 从而影响到定位精度。 针对上述问题, 本文提出了基于聚类和K近邻 算法的井下人员定位算法, 用二分k-me a n s聚类算 法[ 12-13]对井下巷道中采集的数据进行分类, 建立离 线指纹数据库; 用动态修正器反馈井下巷道中的电磁 变化量, 实时修正权重值; 再利用软硬件动态修正加 权 K 近 邻 算 法 (S o f t w a r e a n d H a r d w a r e D y n a m i c C o r r e c t i o n WK N N,S AH D C-WK N N) 进行在线实时位 置估算。实验结果表明, 该无线定位算法提高了定位 系统的稳定性和实时性, 定位精度达到 1 . 6 2m。 1 煤矿井下人员定位系统结构及原理 1. 1 系统结构 煤矿井下人员定位系统由地面控制中心和井下 的光纤交换机、A P( 接入点) 、 动态修正器和井下作 业人员携带的无线移动终端等组成, 如图1所示。 图1 井下人员定位系统结构 Fi g . 1 S t r u c t u r e o f u n d e r g r o u n d p e r s o n n e l p o s i t i o n i n g s y s t e m 1. 2 系统原理 井下人员定位过程包括离线建库和实时定位 2个阶段, 如图2所示。 图2 井下人员定位系统原理 Fi g . 2 P r i n c i p l e o f u n d e r g r o u n d p e r s o n n e l p o s i t i o n i n g s y s t e m 离线建库阶段包含采样和聚类生成离线指纹数 据库。 采样 在井下巷道中选取采样点, 测量该点的物 理位置; 利用自制无线移动终端多次采集该点的 44 工矿自动化 2 0 1 9年第4 5卷 RS S I值, 通过无线网络传输采样数据并存储到原始 信息数据库。按照这种方法采集所有采样点的数 据, 由服务器对所有采样点R S S I求平均值, 将R S S I 值和采样点物理位置信息进行组合并存入原始信息 数据库。在采样的同时, 动态修正器也采集 R S S I 数据, 通过无线网络传输并存储到动态修正数据库。 聚类生成离线指纹数据库 利用二分 k-me a n s 聚类算法对原始信息数据库中所有采样点的数据进 行聚类分析, 得到离线指纹数据库。 实时定位阶段 一方面, 待测点通过无线移动终 端采集 实时 R S S I值 并 传 送 到 在 线 定 位 数 据 库; 另一方面, 动态修正器也实时采集 R S S I值, 传送到 动态修正数据库。利用待测点和动态修正器的离线 数据和实时数据, 用 S AHD C-WKNN 算法计算权 重值, 结合离线指纹数据库中待测点的物理位置信 息估算其实时位置。 2 二分k-me a n s聚类算法 2. 1 原始信息数据库的建立 设m个采 样 点 的 位 置 坐 标 为 (x1,y1), ( x2,y2), , (xm,ym), 在每个采样点能够获得n个 AP节点的R S S I值, 则第i个采样点接收的R S S I 值的组合可表示为 Ri=(Ri1,Ri2,,Rin) ( 1) 式中Ri n为第i个采样点接收的第n 个 A P节点的 R S S I值,i=1,2, ,m。 将采样点坐标及其接收的R S S I值组合起来, 表示为 Mi=(Ri1,Ri2,,Rin,xi,yi) ( 2) 将m个采样点的数据都以式(2) 所示的数据结 构存储到原始信息数据库中, 以备聚类运算时调用。 2. 2 标准k-me a n s聚类算法原理 设聚类数目为k(k≤m), 经过标准k-me a n s聚 类算法生成的簇集合为C={C1,C2,,Ck}, 每个簇 的中心点Ok为 Ok= Mi∈C k Mi |Ck| ( 3) 由式(3) 可得标准k-me a n s聚类算法的目标函 数为 S(C)= k j=1 Mi∈C k ‖Mi-Oj‖ 2 ( 4) 标准k-me a n s聚类算法的目标是找到最小化误 差平方和的聚类结果。具体来说, 首先设置聚类个 数k , 从 数据集中选择k个初始中心点; 然后将每一 个数据对象指派到距离最近的中心点所属的类中, 从而构成k个簇; 计算每个簇中包含的数据对象的 平均值, 将其作为簇中心点的新值; 重复进行指派和 更新这2个步骤, 对簇的分布进行迭代更新, 直到簇 中样本不再发生变化, 整个聚类过程完成。 2. 3 二分k-me a n s聚类算法流程 二分k-m e a n s聚类算法将所有点看作一个簇, 将这个簇分裂为2个簇, 之后迭代这个过程, 每次都 选取一个簇进行分裂, 并将其分为2个新簇, 最终得 到k个簇。 用二分k-me a n s聚类算法生成井下人员定位离 线指纹数据库时, 输入为m个采样点的组合数据 Mi和聚类个数k , 输 出为最终聚类结果。具体算法 流程如下 ( 1)初始化簇表, 将m个采样点组合成一个簇。 ( 2)从簇表中取出一个簇, 设聚类数目为2, 使 用标准k-m e a n s聚类算法对选定的簇进行聚类。 ( 3)从聚类结果中选取误差平方和最小的一组 簇, 将其添加到簇表中。 ( 4)判断簇数量是否达到k个, 若达到则聚类 结束, 否则跳转到步骤( 2) 。 与标准k-me a n s聚类算法相比, 二分k-me a n s 聚类算法的聚类结果更加稳定, 但计算时间也更长。 因为建立的是离线指纹数据库, 在实时定位阶段不 需要重复计算, 离线建库阶段增加的计算时间可以 不用考虑, 所以, 二分k-me a n s聚类算法更适用于建 立无线定位离线指纹数据库[ 14-15]。 3 SA H D C-WK N N算法 S AHD C-WKNN算法中的权重值wi由2个部 分组成[ 16-17] 一 部分是与待测点欧氏距离最小的K (K≤P,P为巷道中安装的动态修正器的数量) 个 采样点的权重值ws i ; 另 一部分是与待测点欧氏距离 最小的K个动态修正器的权重值wh i。 ( 1)设实时定位阶段待测点D接收到的R S S I 信号的平均值集合为( RD1,RD2,,RDn) , 则 D点和 第i(i=1,2, ,K) 个采样点的欧氏距离为 di= n j=1 ( RDj-R i j 槡 )( 5) 与待测点欧氏距离最小的K个采样点的权重 值ws i为 wsi= w′si k i=1w ′si ( 6) w′si= 1 d 2 i ( 7) (2)与待测点欧氏距离最小的K个动态修正 54 2 0 1 9年第4期莫树培等 基于聚类和 K近邻算法的井下人员定位算法 器接收的数据包含离线数据和实时数据。离线数据 集合GOL为 GOL= R′11R′12 R′1n R′21R′22 R′2n R′K1R′K2 R′ 熿 燀 燄 燅 Kn ( 8) 式中R′ Kn为离线建库阶段第K个动态修正器接收 的第n 个 A P节点的 R S S I值。 对GOL中的数据求平均值 珚 R′i= 1 n n j=1R ′ i j ( 9) 实时数据集合GRT为 GRT= R″11R″12 R″1n R″21R″22 R″2n R″K1R″K2 R″ 熿 燀 燄 燅 Kn ( 10) 式中R″Kn为实时定位阶段第K个动态修正器接收 的第n 个 A P节点的R S S I值。 对GRT中的数据求珚 R′i的标准差si si= 1 n n j=1 ( R″ i j -珚R′i) 槡 2 ( 11) 与待测点欧氏距离最小的K个动态修正器的 权重值wh i为 whi= w′hi K j=1w ′hj ( 12) w′hi= si 珚 R′i ( 13) 由式(6) 和式(1 2) 得到权重值wi为 wi= 1 2( wsi+whi) ( 1 4) (3)结合离线指纹数据库中待测点的物理位置 信息, 求待测点的实时位置(X, Y) (X, Y)= K i=1w i(xi,yi) ( 1 5) 4 实验分析 4. 1 设备选用及测点布置 考虑到煤矿井下各类机电设备电磁干扰严重、 杂物堆放、 人员设备移动等因素, 选取贵州某煤矿井 下回风巷、 切眼和运输巷共计8 9 6m的巷道进行实 验, 如图3所示。巷道中共安装了9个动态修正器, 回风巷中3个, 切眼中2个, 运输巷中4个。 根据巷道的作业环境和走向, 将A P之间的距 离控制在8 0m内, 实现井下巷道无线网络全覆盖 且无盲点。选用本质安全型无线 A P,2. 4GH z频 段, 同时限制W i F i发射功率不超过1 0 0mW。自制 图3 井下实验巷道 F i g . 3 U n d e r g r o u n d e x p e r i m e n t r o a d w a y 动态 修 正 器 和 无 线 移 动 终 端 主 要 由 主 控 芯 片 S T C 8 A 8 K 6 4 S 4 A 1 2、W i F i模块 E S P 8 2 6 6 E X和电源 模块等组成, 经过防爆测试, 符合煤矿井下设备本质 安全性要求[ 18-22]。 根据巷道截面宽度, 笔者设计并定制了一种由 尼龙棒组成的田字格形状的可拆卸塑料架子( 可拆 解成日字形或口字形) , 长宽均为3 m, 一次可放 9个无线移动终端, 能同时采集9个采样点, 大大提 高了采集效率。在实际采样过程中, 采样点与采样 点的间距基本控制在1. 5~2m, 每个采样点采集 1 0次R S S I数据, 共1 7 2 4个采样点。同时, 测试人 员另外采集了1 0 2次测试数据, 组成测试数据集。 4. 2 定位结果 将采样点数据全部存放到原始信息数据库, 进 行相应处理后, 通过二分k-me a n s聚类算法建立离 线指纹数据库。设置k值为19, 用WKNN算法 和S AHD C-WKNN算法进行实时定位, 定位成功 率( 误差在2m以下) 对比见表1。 表1 W KNN算法和S AHD C-WKNN算法定位成功率对比 T a b l e 1 C o m p a r i s o n o f p o s i t i o n i n g s u c c e s s r a t e o f WKNN a l g o r i t h m a n d S AHD C-WKNN a l g o r i t h m k 定位成功率/% WKNN算法S AHD C-WKNN算法 1 32. 0 4 4 5. 2 6 2 4 2. 1 5 6 7. 5 8 3 6 3. 7 6 8 1. 0 5 4 5 0. 0 1 7 6. 5 6 5 4 7. 2 4 7 2. 9 1 6 4 0. 2 6 6 4. 3 8 7 3 9. 2 6 5 8. 3 2 8 3 6. 5 9 5 4. 6 4 9 33. 7 1 4 8. 5 9 从表1可看出, 当k值取3时能达到较好的定 位效果, 同时计算量适中, 因此, S AHD C-WKNN 算 法中k值定为3。 64 工矿自动化 2 0 1 9年第4 5卷 实时定位阶段, 作业人员携带无线移动终端进 入测试区域后, 将实时接收到的 R S S I值发送回服 务器, 系统将该 R S S I值存放到在线定位数据库中, 再调用权值修正程序、 动态修正 WKNN 程序进行 井下作业人员位置估算。设( xs,ys) 为井下作业人 员实际的物理位置坐标( 实测坐标) , ( xe,ye ) 为 定位 系统估算出的井下作业人员物理位置坐标( 估算坐 标) , 用式( 1 6) 、 式(1 7) 计算系统定位标准误差E(i) 和平均误差珚 E , 结果见表2 。 E(i)= [ xs-xe] 2 +[ys-ye]槡 2 ( 16) 珚 E= n i=1E ( i) n ( 17) 表2 SAHD C-WKNN算法定位结果 T a b l e 2 p o s i t i o n i n g r e s u l t s o f S AHD C-WKNN a l g o r i t h m m 序号 实测坐标 (xs, ys) 估算坐标 (xe, ye) 标准误差 1 (3 4. 0 0,4. 0 0)(3 4. 6 6,3. 8 2) 0. 6 8 2 (1 2 5. 0 0,3. 5 0)(1 2 2. 1 2,3. 3 1) 2. 8 9 3 (2 4 8. 0 0,2. 0 0)(2 4 9. 9 2,2. 2 8) 1. 9 4 4 (3 2 1. 0 0,0. 5 0)(3 1 9. 1 1,0. 5 5) 1. 8 9 5 (4 5 2. 0 0,2. 5 0)(4 5 3. 7 3,2. 7 2) 1. 7 4 6 (5 7 9. 0 0,1. 5 0)(5 8 0. 5 4,1. 6 5) 1. 5 5 7 (6 3 5. 0 0,2. 0 0)(6 3 4. 2 1,1. 9 8) 0. 7 9 8 (7 2 4. 0 0,2. 5 0)(7 2 3. 0 1,2. 3 7) 1. 0 0 9 (8 0 5. 0 0,2. 0 0)(8 0 1. 7 6,2. 3 4) 3. 2 6 1 0 (8 7 8. 0 0,3. 0 0)(8 7 8. 4 2,2. 8 1) 0. 4 6 从表2可看出, 最小标准误差为0. 4 6m, 最大 标准误差为3. 2 6m, 系统的平均误差为1. 6 2m。 4. 3 对比分析 通过在线定位数据库中所有待测点的数据, 采 用NN( 最近邻) ,KNN(K近邻) ,WKNN( 加权K近 邻) ,S AHD C-WKNN等4种近邻算法进行定位, 用 每种算法进行5 0 0次独立定位并计算误差, 得到定 位误差对比, 如图4所示。 从图4可看出,S AHD C-WKNN算法的误差比 最低, 定位精度最高。 实时性也是井下人员定位系统的一个很重要的 性能指标, 将“ 未聚类分析+K N N” [11]、 “ 聚类分析+ WKNN”与 本 文 提 出 的 “ 聚 类 分 析 +S AHD C- WKNN” 算法的精度和实时性进行对比, 结果 见 表3。 从表3可看出, 与“ 未聚类分析+KNN” 相比, 本文算法的精度更高, 实时性更好; 与“ 聚类分析+ WKNN” 相比, 本文算法的运算时间略有增加, 但定 图4 定位误差对比 F i g . 4 C o m p a r i s o n o f p o s i t i o n i n g e r r o r 表3 算法定位精度与实时性比较 T a b l e 3 C o m p a r i s o n o f p o s i t i o n i n g a c c u r a c y a n d r e a l-ti m e p e r f o r m a n c e o f t h e a l g o r i t h m s 算法平均误差/m运算时间/s 未聚类分析+KNN 3. 1 5 1. 1 6 3 聚类分析+WKNN 2. 5 8 0. 4 0 9 聚类分析+S AHD C-WKNN 1. 6 2 0. 4 9 5 位精度提高了3 7. 2 1%。因此, 本文提出的“ 聚类分 析+S AHD C-WKNN” 算法的定位效果最好, 更能 满足煤矿井下人员定位要求。 5 结论 ( 1)针对煤矿井下特定的地理和通信环境, 提 出了基于聚类和K近邻算法的井下人员定位算法。 通过二分k-me a n s聚类算法对离线数据进行分类, 提高了定位算法实时性; 通过在巷道中增加动态修 正器来实时反馈井下电磁环境情况, 动态修正权重 值, 提高了定位算法精度。 ( 2)实验分析结果表明, 所提定位算法的最小 标准误差为0. 4 6m, 最大标准误差为3. 2 6m, 平均 误差为1. 6 2m。对比分析结果表明, 与未进行聚类 分析的算法相比, 本文算法的精度更高, 实时性更 好; 与未动态修正权重值的算法相比, 本文算法的运 算时间略有增加, 但定位精度提高了3 7. 2 1%。 ( 3)在接下来的工作中, 应继续优化无线定位 算法, 提高定位精度, 在定位性能不变的前提下减少 采样点数量。 参考文献(R e f e r e n c e s) [1] 孙继平.煤矿信息化自动化新技术与发展[J].煤炭科 学技术,2 0 1 6, 4 4(1) 1 9-23. S UN J i p i n g .N e w t e c h n o l o g y a n d d e v e l o p m e n t o f m i n e i n f o r m a t i z a t i o n a n d a u t o m a t i o n[J]. C o a l S c i e n c e 74 2 0 1 9年第4期莫树培等 基于聚类和 K近邻算法的井下人员定位算法 an d T e c h n o l o g y,2 0 1 6,4 4(1) 1 9-23. [2] 闫 雷 兵, 陆 音, 张 业 荣.基 于 改 进 最 小 二 乘 算 法 的 T D OA/AOA定 位 方 法 [J].电 波 科 学 学 报,2 0 1 6, 3 1(2) 3 9 4-40 0. YAN L e i b i n g,L U Y i n,Z HANG Y e r o n g . I m p r o v e d l e a s t-s q ua r e s a l g o r i t h m f o r T D OA/AOA-ba s e d l o c a l i z a t i o n[J]. C h i n e s e J o u r n a l o f R a d i o S c i e n c e, 2 0 1 6,3 1(2) 3 9 4-40 0. [3] 孙 继 平, 蒋 恩 松.基 于W i F i信 号 二 次 扩 频 的 矿 井 T OA测距方法[J].工矿自动化,2 0 1 7,4 3(1 0) 5 3-58. S UN J i p i n g,J I ANG E n s o n g . M i n e T OA r a n g i n g m e t h o d b a s e d o n r e-s p re a d s p e c t r u m t o W i F i s i g n a l [J]. I n d u s t r y a n d M i n e A u t o m a t i o n,2 0 1 7,4 3(1 0) 5 3-58. [4] 孙哲星, 孙继平.异步测时矿井人员精确定位方法 [J].煤炭学报, 2 0 1 8,4 3(5) 1 4 6 4-14 7 0. S UN Z h e x i n g,S UN J i p i n g . U n d e r g r o u n d c o a l m i n e a c c u r a t e p e r s o n n e l p o s i t i o n i n g m e t h o d b a s e d o n a s y n c h r o n o u s t i m e-me a s u r i n g[J]. J o u r n a l o f C h i n a C o a l S o c i e t y,2 0 1 8,4 3(5) 1 4 6 4-14 7 0. [5] F ANG X u m i n g,J I ANG Z o n g h u a,NAN L e i,e t a l . O p t i m a l w e i g h t e d K-ne a r e s t n e i g h b o r a l g o r i t h m f o r w i r e l e s s s e n s o r n e t w o r k f i n g e r p r i n t l o c a t i o n i n n o i s y e n v i r o n m e n t[J]. I E T C o mm u n i c a t i o n s,2 0 1 8,1 2(1 0) 1 1 7 1-11 7 7. [6] 孙哲星.基于时间测距的矿井人员定位方法研究[J]. 工矿自动化,2 0 1 8, 4 4(4) 3 0-33. S UN Z h e x i n g . R e s e a r c h o n m i n e p e r s o n n e l p o s i t i o n i n g m e t h o d b a s