-
LED灯将接收到的电信号转化为光信号进行传输,光电二极管(photo diode, PD)将接收到的光信号转化为电信号,使用空间角度和光强度的关系表示朗伯光源的光强分布[13-14]。参考文献[15]中通过计算得出光信号反射次数为3次~10次时对系统的影响可以忽略不计,由于室内情况复杂漫反射为主要反射源,大部分反射信号需3次及以上反射才能被PD接收端接收。因此本文中只考虑直射链路,对反射信号忽略不计。
在10m×10m×3m的室内空间中设置9盏LED灯作为信号发射器,将其呈等间隔对称分布安装在距离地面3m的天花板上。将PD作为接收器安装在移动端上。LED发射器和接收器定位系统模型如图 1所示。
接收光功率Pr与发射端功率Pt的关系为:
$ P_{\mathrm{r}}=H(0) P_{\mathrm{t}} $
(1) 式中,H(0)为信道增益,信道增益可以表示为[16]:
$ \begin{array}{l} H(0)= \\ \left\{\begin{array}{l} \frac{(l+1) A}{2 \pi d^{2}} \cos ^{l} \varphi \cos \psi f(\psi) g(\psi), \left(0 \leqslant \psi \leqslant \psi_{\max }\right) \\ 0, \left(\psi>\psi_{\max }\right) \end{array}\right. \end{array} $
(2) 式中,A为PD的面积; d为LED到PD的相对距离; φ为发散角; ψ为接收角; ψmax为PD最大视场角; f(ψ)为PD的光滤波器增益;g(ψ)为聚光器增益。将(1)式和(2)式联立即可推导出LED与PD间的距离d; l为朗伯散射系数,朗伯系数大小由LED灯的半功率角ψ1/2决定:
$ l=\frac{-\ln 2}{\ln \left(\cos \psi_{1 / 2}\right)} $
(3) $ I_{\theta}=I_{0} \cos ^{l} \psi $
(4) 式中, Iθ为发射角θ的发光强度;I0为垂直于发光面的发光强度。
RSSI定位是根据PD接收到来自多个发射器发出的光源信号强度得到欧氏距离,利用三边法求解出目标位置的对应坐标点。已知参考点与待测节点坐标,代入三边定位算法公式可计算出3个圆心到达待测点的距离d。由于室内充斥着电磁的干扰,三边测距法计算3个圆不会理想地交于一点, 因此通过黄金正弦算法寻找最优解。
-
黄金正弦算法(golden sine algorithm, Gold-SA)是一种新型元启发式算法[19-20]。黄金正弦算法是在正弦函数上沿着同一搜索路线,将单位圆上的每一个点访问一遍,当位置更新时加入黄金分割率进而去缩小解的范围,使得区域范围更小,进而可以提升算法搜索速度。Gold-SA算法在更新个体位置时,引入黄金分割率,分别设系数x1和x2, 将系数加入到位置更新公式中,缩小搜索空间使得寻优速度更快,黄金分割系数x1和x2为:
$ \left\{\begin{array}{l} x_{1}=a(1-c)+b c \\ x_{2}=a c+b(1-c) \end{array}\right. $
(5) 式中,a和b为黄金分割率的初始值,a=-π, b=π; c为分割率,c=($\sqrt 5 - 1$)/2。
黄金正弦算法中结合黄金分割率的系数,在搜索空间中寻找最优解,其位置更新公式为[21]:
$ V_{i}^{(t+1)}=V_{i}^{(t)}\left|\sin r_{1}\right|-r_{2} \sin r_{1}\left|x_{1} D_{i}^{(t)}-x_{2} V_{i}^{(t)}\right| $
(6) 式中,Vi(t)为第i个个体第t次迭代后的位置;Di(t)为第i个个体第t次迭代后的最优位置;r1为[0, 2π]内的随机数;r2为[0, π]内的随机数。
-
本文作者提出将Levy飞行变异机制结合自适移动因子改进黄金正弦算法,对位置更新公式进行变异,提高算法种群的多样性,能够快速跳出局部最优。Levy飞行机制的核心是短距离游走和长距离跳跃交替变化,当处于短距离游走时种群的多样性将会提高,处于长距离跳跃时种群搜索具有方向的多样性,搜索更为详细[22-23]。L为Levy飞行的移动步长参数,计算公式为:
$ L \approx \frac{\lambda \Gamma(\lambda) \sin (\pi \lambda / 2)}{\pi} \frac{1}{s^{1+\lambda}}, \left(s \gg s_{0}>0\right) $
(7) 式中,λ为缩放因子,参考文献[24]中建议取λ=1.5;Γ(λ)为标准的伽玛函数;s为移动步长;s0为起始步长,即s≫s0>0为s→∞。
在黄金正弦算法全局搜索中引入Levy飞行变异机制帮助算法跳出局部最优,但Levy飞行变异机制无法满足算法需要较高的收敛速度。本文中提出Levy飞行机制结合自适应移动因子对黄金正弦算法位置更新公式改进,进而提高算法的收敛速度。自适应移动因子α的公式为:
$ \alpha(t)=\sqrt{{\mathrm{e}}^{-q}+\frac{1}{j}} $
(8) 式中,j为当前的迭代次数;q=1+j/M,M为最大迭代次数。
使用Levy飞行变异机制结合自适应移动因子对Gold-SA位置更新公式进行变异,变异后的位置更新公式为:
$ \begin{array}{c} V_{i}^{(t+1)}=V_{i}^{(t)}\left|\sin r_{1}\right|-r_{2} \sin r_{1}\left|x_{1} D_{i}^{(t)}-x_{2} V_{i}^{(t)}\right|+ \\ \alpha {\boldsymbol{L}}(m, n) \times\left({\boldsymbol{x}}_{\text {best }}-{\boldsymbol{x}}_{i}\right) \end{array} $
(9) 式中,L(m, n)是关于多个未知量下的步长参数的矩阵,其中m表示控制常数,n为未知量个数,取m=1.5, n=3; xbest表示每次迭代最优解矩阵; xi表示当前迭代得出的解的矩阵。
基于Levy变异机制改进Gold-SA算法实现室内可见光定位的流程如下: (1)PD接收器接收LED光信号,计算出LED与PD间的距离d; (2)设计初始种群个体数为N, 最大迭代次数M, 空间维数D的初始种群; (3)计算每个个体的适应度的值Pbest,并将最优个体对应的适应度值存入xbest; (4)根据(9)式更新个体的位置; (5)更新个体,比较个体的适应度值取出最佳适应度值个体作为下一次迭代的解; (6)判断是否已经达到设置迭代次数,若达到则结束,若未达到则转向第(4)步; (7)得到全局最优解。
-
选取体积为10m×10m×3m的室内空间作为仿真测试空间,设置9盏LED灯按照3×3的网格排列方式排布。在室内随机放置一个待测点,分别选取最小二乘(least square mathod, LMS)的RSSI定位算法、PSO的RSSI定位算法、Gold-SA的RSSI定位算法和改进的Gold-SA的RSSI定位算法,比较各算法定位产生的定位误差。设置PD有效接收面积为1cm2。仿真参量表如表 1所示。
Table 1. Simulation parameters
parameter value room size 10m×10m×3m filter gain 1 effective receiving area of PD 1cm2 luminous intensity of LED lamp 0.73cd interval of LED lamp 3m field of view of receiving angle 60° height of receiving plane 0.73m 在3维空间中分别选取LSM定位算法、PSO定位算法、Gold-SA定位算法、改进后的Gold-SA定位算法进行仿真分析,仿真结果如图 2所示。
图 2a中采用LMS定位算法时,定位误差为2.5924m,由于室内环境复杂充满电磁干扰,导致精度低不能做到精确定位。由图 2b和图 2c看出, 基于PSO的定位算法误差为5.9993cm,基于Gold-SA的定位算法误差为5.9866cm,虽然上述两种算法相比于LMS算法精度提高很多,但仍存在收敛速度慢的缺点。改进Gold-SA定位算法如图 2d所示,定位算法误差为1.4703mm,相比PSO算法和未改进的Gold-SA算法,改进后Gold-SA算法的定位精度分别提升97.54%和95.52%。Levy变异机制增加算法搜索空间的多样性,自适应移动因子有效地避免群智能算法收敛早熟的缺点提高算法收敛速度,使算法误差明显减小,实现精确定位。在仿真参量不变的情况下对4个目标点进行同时定位,多目标仿真结果如图 3所示。
由图 3可以看出, 采用改进Gold-SA算法对4个目标点同时定位,定位误差分别为1.170cm, 1.243cm, 1.197cm, 1.470cm,平均定位误差为1.270cm。由图 3仿真分析可知, 改进后的Gold-SA定位算法可应用在多目标定位过程。
-
分别对PSO算法、Gold-SA算法和改进后的Gold-SA算法进行单目标迭代速度的分析,设置迭代数为200次,初始种群数量为20,采用3×3网格型LED灯源排布方案,在保证使用最少LED灯源且定位精度较高的前提下分析算法收敛速度,仿真结果如图 4所示。
由图 4可直观看出:采用Levy飞行变异机制结合自适应移动因子改进的Gold-SA算法收敛速度大幅度提升。改进后的Gold-SA算法经过40次迭代后就快速收敛到最低适应度值。基于PSO的RSSI定位算法110次迭代后收敛到最低适应度值,Gold-SA算法经过200次迭代后还未达到最低适应度值。因此,说明改进后的算法在快速定位方面更具有优势。
-
选取改进Gold-SA算法、未改进的Gold-SA算法及PSO算法,在定位误差为0cm~2cm区间内仿真分析3种定位算法单目标定位时间,结果如图 5所示。
由图 5中明显看出:3种算法在相同定位区间内,改进后的算法定位误差在0cm~2cm的范围内定位时间均比未改进的Gold-SA算法和PSO算法短,由此说明改进的Gold-SA算法可实现快速定位。当定位误差为0.2cm时,改进的Gold-SA定位时间为3.3s比传统Gold-SA算法快0.5s,比PSO算法快0.9s;定位误差为0.6cm时,改进后的算法定位时间为2.2s比传统Gold-SA快0.8s,比PSO算法快1.6s;定位误差为1.5cm时,改进后的算法定位时间为2.4s比传统Gold-SA算法快0.6s,比PSO算法快1s。由此可以说明改进后的Gold-SA算法相比于其它两种算法所需要的定位时间更短,改进后的Gold-SA算法整体优于未改进的Gold-SA算法和PSO算法。
-
为进一步验证改进后的黄金正弦算法优越性,分别对改进的Gold-SA算法、Gold-SA算法和PSO算法分别进行200次单目标定位误差计算。采用直方图更直观地观测定位算法的误差分布,如图 6所示。
由图 6可看出:与Gold-SA算法和PSO算法相比,改进后的Gold-SA算法误差大多集中在0cm~2cm范围内,当误差大于2cm时, 误差点个数呈大幅下降趋势,在6cm~8cm范围内没有误差点分布。Gold-SA算法误差集中在2cm~4cm范围内,PSO算法误差主要分布在2cm~4cm, 4cm~6cm范围内。改进后的Gold-SA算法相较与Gold-SA算法和PSO算法定位性能分别提高64.3%和68.9%。因此,改进后的算法相比于其它两种算法定位误差更小、定位效果更好。本文中对改进Gold-SA算法进行80次的误差计算,依然采用在10m×10m×3m的室内空间中3×3网格型LED灯源排布方案,误差折线图如图 7所示。
图 7中对改进的Gold-SA算法对单目标进行80次误差计算,经过多次仿真可以得出改进后的算法误差稳定在0cm~1cm区间内,最小误差可达1.4703mm。
由以上仿真综合分析可知,相比于PSO算法、Gold-SA算法和LMS算法,改进的Gold-SA算法具有收敛速度快、误差小趋于稳定等优点,改进后的Gold-SA算法相比PSO算法和Gold-SA算法定位性能更具优越性,满足室内高精度快速定位要求。
基于自适应Levy飞行的黄金正弦可见光定位研究
Research on golden sine visible light positioning based on adaptive Levy flight
-
摘要: 为了提高室内可见光定位的精度,采用了基于Levy飞行变异机制、结合自适应移动因子、改进黄金正弦算法的接收信号强度指示可见光定位方法, 将室内屋顶的发光二极管灯按3×3网格状排布,接收到光强信号通过朗伯模型得到未知节点与参考节点的距离,并采用Levy飞行变异机制提升算法搜索空间的多样性,结合自适应移动因子提高算法收敛速度,使得个体更新受局部极值约束力下降。结果表明, 改进算法平均定位误差为1cm,平均迭代次数40次~80次; 改进黄金正弦算法的定位速度和定位精度均得到提升。该研究对室内大型场所实时、快速精确定位有帮助。Abstract: To improve the accuracy of indoor visible light positioning, the received signal strength indication visible light positioning method based on Levy flight mutation mechanism combined with adaptive movement factor improved golden sine algorithm was adopted. The light-emithing diode lights on the indoor roof were arranged in a 3×3 grid, and the distance d between the unknown node and the reference node was obtained through the Lambertian model after receiving the light intensity signal. The Levy flight mutation mechanism was used to increase the diversity of the algorithm's search space, and the adaptive movement factor is used to increase the algorithm's convergence speed, so that individual updates were less restricted by local extreme values. The simulation shows that the average positioning error of the improved algorithm is 1cm, and the average number of iterations is 40~80. Therefore, the speed and precise positioning of the improved golden sine algorithm are improved, which is helpful for real-time and rapid positioning of large indoor places.
-
Table 1. Simulation parameters
parameter value room size 10m×10m×3m filter gain 1 effective receiving area of PD 1cm2 luminous intensity of LED lamp 0.73cd interval of LED lamp 3m field of view of receiving angle 60° height of receiving plane 0.73m -
[1] YUAN Sh. Study and implementation of mobile location in the WLAN[D]. Xi'an: Xidian University, 2014: 35-42 (in Chinese). [2] DENG Zh L, YU Y P, YUAN X, et al. Situation and development tendency of indoor positioning[J]. China Communications, 2013, 10(3): 42-55 (in Chinese). doi: 10.1109/CC.2013.6488829 [3] LIU R, YANG B, ZIO E, et al. Artificial intelligence for fault diagnosis of rotating machinery: A review[J]. Mechanical Systems & Signal Processing, 2018, 108(8): 33-47. [4] TANG Q Sh. Research on indoor localization algorithms based on visible light communication[D]. Xi'an: Xidian University, 2016: 98-102(in Chinese). [5] BAO J J, HUO Zh L, XU W, et al. A wireless location method with high precision for underground personnel tracking[J]. Industrial and Mining Automation, 2009, 35(10): 18-21 (in Chinese). [6] O'BRIEN D C, ZENG L B, LE-MINH H, et al. Visible light communications: challenges and possibilities[C]//2008 IEEE 19th International Symposium on Personal, Indoor and Mobile Radio Communications. New York, USA: IEEE, 2008: 1-5. [7] SHI W, QI X, LI J, et al. Simple solution to the optimal deployment of cooperative nodes in tow-dimensional TOA-based localization system[J]. Eurasip Joumal on Wireless Communication & Networking, 2017 (1): 1-16. [8] HOU Y, XUE Y, CAO C, et al. A RSS/AOA based indoor positioning system with a single LED lamp[C]//IEEE International Conference on Wireless Communication & Signal Processing (WCSP). New York, USA: IEEE, 2015: 45-52. [9] YANG G W, HUANG Zh B, FAN B, et al. Indoor positioning and orientating system based on visible light communication[J]. Journal of Communications, 2020, 41(12): 162-170 (in Chinese). [10] CHENG X, BAI Y H, GUO X, et al. Visible light indoor positioning technology based on improved particle swarm optimization algorithm[J]. Optical Communication, 2019, 43(4): 8-11(in Chinese). [11] LIU Zh G, ZHAO L, ZHU T, et al. Research of indoor positioning and illuminating systems based on visible light[J]. Laser Technology, 2018, 42(5): 646-650(in Chinese). [12] KARMY M, ELSAYED S, ZEKRY A. Performance enhancement of an indoor localization system based on visible light communication using RSSI/TDOA hybrid technique[J]. Journal of Communications, 2020, 15(5): 379-389. [13] CHENG Sh L, YAN K, LI Zh, et al. Visible-light indoor positioning technology based on parameters optimization of lambertain model[J]. Research on optical communication, 2018(5): 69-73 (in Chinese). [14] ZHOU P, ZHAO Y H, HU F Ch, et al. Research status of machine learning based signal processing in visible light communication[J]. Laser & Optoelectronics Process, 2020, 57(1): 010001(in Chinese). [15] DING J P. Visible light communications indoor channel modeling & performance optimization[D]. Beijing: Beijing University of Posts and Telecommunications, 2013: 45-49(in Chinese). [16] KAHN J M, BARRY J R. Wireless infrared communication[J]. Proceedings of the IEEE, 1997, 85(2): 265-298. doi: 10.1109/5.554222 [17] NADEEM U, HASSAN N U, PASHA M, et al. Indoor positioning system designs using visible LED lights: Performance comparison of TDM and FDM protocols[J]. Electronics Letters, 2015, 51(1): 72-74. doi: 10.1049/el.2014.1668 [18] NIU Y F, SUN D D. Application of visible light communication in indoor positioning of intelligent lighting[J]. Journal of Lighting Engineering, 2016, 27(6): 104-107 (in Chinese). [19] TANYILDIZI E, DEMIR G. Golden sine algorithm: A novel math-inspired algorithm[J]. Advances in Electrical and Computer Engineering, 2017, 17(2): 71-78. doi: 10.4316/AECE.2017.02010 [20] YU J F, LIU Sh. Golden sine simulated annealing algorithm for low carbon capacitated vehicle routing problem[J]. Science Technology and Engineering, 2020, 20(11): 4202-4209(in Chinese). [21] ZHOU Y R, LI N, ZHOU F H. Application of golden sine algorithm in hydrogeological parameter optimization[J]. People's Pearl River, 2020, 41(6): 117-120 (in Chinese). [22] GANDOMI A H, YANG X S, ALAVI A H. Cuckoo search algorithm: A metaheuristic approach to solve structural optimization problems[J]. Engineering with Computers, 2013, 29(1): 17-35. doi: 10.1007/s00366-011-0241-y [23] GHAEMI M, ZABIHINPOUR Z, ASGARI Y. Computer simulation study of the levy flight process[J]. Physica, 2010, A388(8): 1509-1514. [24] ZHAO Ch Y. Application research of flower pollination algorithm[D]. Nanning: Guangxi University for Nationalities, 2017: 23-26 (in Chinese).