高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于卡尔曼滤波的接收信号强度指示差值定位算法

耿友林 解成博 尹川 郭兰图 王先义

胡淼, 张子石, 章越峰, 许蒙蒙, 周雪芳, 沈成竹. 用于雷电测向的三磁环天线设计和研究[J]. 电子与信息学报, 2022, 44(8): 2878-2883. doi: 10.11999/JEIT210476
引用本文: 耿友林, 解成博, 尹川, 郭兰图, 王先义. 基于卡尔曼滤波的接收信号强度指示差值定位算法[J]. 电子与信息学报, 2019, 41(2): 455-461. doi: 10.11999/JEIT180268
HU Miao, ZHANG Zishi, ZHANG Yuefeng, XU Mengmeng, ZHOU Xuefang, SHEN Chengzhu. Research of Three-magnetic Loop Antenna for Lightning Direction Finding[J]. Journal of Electronics & Information Technology, 2022, 44(8): 2878-2883. doi: 10.11999/JEIT210476
Citation: Youlin GENG, Chengbo XIE, Chuan YIN, Lantu GUO, Xianyi WANG. Received Signal Strength Indication Difference Location Algorithm Based on Kalman Filter[J]. Journal of Electronics & Information Technology, 2019, 41(2): 455-461. doi: 10.11999/JEIT180268

基于卡尔曼滤波的接收信号强度指示差值定位算法

doi: 10.11999/JEIT180268
详细信息
    作者简介:

    耿友林:男,1966年生,教授,研究方向为复杂介质目标电磁散射的解析和数值方法、复杂介质中的电波传播特性以及小天线设计和辐射特性等

    解成博:男,1994年生,硕士生,研究方向为信号源定位

    尹川:男,1987年生,讲师,研究方向为电磁成像

    通讯作者:

    解成博 chengbo_x@163.com

  • 中图分类号: TN911.7

Received Signal Strength Indication Difference Location Algorithm Based on Kalman Filter

  • 摘要:

    针对频谱监测系统中被监测信号无法控制并且没有任何先验知识,只能通过对信号被动监测,即接收与处理信号来估计信号源位置的要求,该文提出一种基于接收信号强度指示差值(RSSID)的定位算法,并利用卡尔曼滤波提高其定位精度。该文将两监测站之间的RSSID转换成信号源到两监测站的距离之比,根据距离之比构造定位方程矩阵,进而利用最小二乘法求取信号源位置。仿真结果表明:所提算法比经典RSSI定位算法性能更优,降低了环境因素对定位精度的影响,并且能更好地满足参数较少的定位服务需求,可以有效地应用于频谱监测系统中。同时,卡尔曼滤波可以有效改善系统的定位精度,达到预期的定位效果。

  • 雷电是一种发生在大气中的天气现象,发生时会产生瞬态大电流、高电压和强电磁脉冲辐射等[1]。根据卫星监测和气象部门的数据统计,全球每秒大约有2000 多个雷电发生[2]。近年来,由于电子、通信设备发展迅速,导致雷电灾害由以往的森林火灾和人员伤亡形式逐渐转向对电子设备以及大型通讯网络设备的干扰与破坏,对军事、铁路、航天航空等领域造成严重影响。因此开展雷电探测研究,提供雷电预警,对防雷减灾的工作具有十分重要的指导意义[3-6]。直至今日,人类对雷电探测技术的研究已经发展了200多年,从20世纪20年代起,美国开始研究利用磁环天线来对雷电位置进行精确探测,成为雷电探测技术发展史上的一个里程碑,代表了雷电探测技术逐渐走向成熟[7-11]

    目前雷电定位系统(Lightning Location System, LLS)主要采用正交磁环天线(Orthogonal Magnetic Loop Antenna, OMLA)作为核心探测器件[12]。为提高测向精度,以往对OMLA的研究多注重于减小其测向角度误差,如对由雷电回击通道和地面不垂直造成的极化误差研究[13]、由雷电测向设备附近地形地势造成的场地误差研究[14]、由天线磁环非正交性造成的误差研究[15]以及OMLA自身结构误差的研究等[16],而对采用新型磁环天线结构设计提高天线测向精度的研究则相对较少。随着LLS探测精度要求的不断提高,采用正交磁环天线探测雷电方向的测量精度达到瓶颈。具体而言,在排除外界条件造成天线的测角误差后,虽然可以通过引入修正矩阵来降低磁环的非正交性误差[15],或引入补偿系数来矫正正交磁环天线的一致性误差[16],在一定程度上提高测向精度,但经矫正后的正交磁环天线仍存在一定的结构误差角,导致测角误差随雷电方向角呈波浪形变化。即当雷电信号方向平行于其中一个磁环时,测角误差达到最大,最大值约等于天线结构误差角的值。因此研究如何设计新型磁环天线结构以减轻结构误差角对测角的影响,对于提升雷电测向精度具有一定的意义。

    本文提出一种由3个磁环两两呈60°组成的新型三环天线,由于其特殊结构,可以有效降低结构误差角对测角的影响,提高天线测向精度。首先根据三磁环天线特殊结构推导了低频信号情况下三环天线的测角公式,建立了结构误差角与测角误差的仿真模型;然后通过实验比对三环天线与同尺寸正交磁环天线测角精度。本文最后部分给出分析结论。

    雷电与地面的回击通道可以被看作一个垂直的电偶极子,当地面为理想导体时,雷电电磁脉冲仅激发横磁波[17]。三磁环天线结构如图1所示,假设H为垂直地面的闪电回击通道,3个磁环两两呈60°组成天线来进行雷电源方位角的探测,中轴线O与回击通道H平行。

    图 1  三磁环天线结构示意图

    设HO平面(即雷电来向)与A环、B环、C环的夹角分别为θ, 60°–θ, 120°–θ,3个磁环受雷电电磁脉冲激发产生的自感电动势峰值分别为EA, EB, EC。再假设磁环中电动势、电流沿逆时针为正方向,且磁环自感磁链的正方向和电流的正方向之间符合右手螺旋关系,根据法拉第电磁感应定律,A环感应电动势计算公式为

    EA=2πλBNSμrsinθ
    (1)

    式中,λ为信号波长;B为信号磁感应强度;N为天线匝数;S为天线面积;μr为磁环的相对磁导率;θ为雷电信号与A环的夹角,同理将B环、C环与雷电来向夹角代入式(1)可得其各自感应电动势EB, EC

    考虑3个磁环间还存在互感电动势:3个磁环受雷电电磁脉冲激发产生自感电动势后,每个磁环内变化的感应电流会在其他两个磁环上激发互感电动势。设VA, VB, VC分别为磁环A, B, C两端实际电动势;eAB为A环电流在B环激发的感应电动势,同理eBC, eAC分别为B环电流在C环激发的感应电动势以及A环电流在C环激发的感应电动势,则有

    VA=EA+eBA+eCAVB=EB+eAB+eCBVC=EC+eAC+eBC}
    (2)

    考虑雷电信号一部分处于200 kHz以下的低频频段,在此情况下,磁环间互感电动势很小,因此3个磁环间互感可忽略不计。根据欧姆定律将3个磁环的感应电动势两两组合进行计算,分别可得3组雷电测向角度θ1,θ2,θ3

    θ1 = arctan(3EAEA2EB)θ2 = arctan(3EAEA+2EC)θ3 = arctan(3(ECEB)EC+EB)}
    (3)

    对求得的3个角度求算数平均值作为最终测向结果。此外由于三角函数的周期性特征,由式(3)所计算得到的方向角θ存在180°的二义性,可使用垂直极化的电场天线判断雷电信号极性,进一步确定方向角符号。

    为验证三磁环天线结构误差角对测角误差的影响,建立测角误差关于结构误差角的仿真模型。结构误差角由天线磁环夹角不理想造成,理论上天线的3个磁环保持两两呈60°,但实际情况下由于加工、磁环扭曲等问题,使得磁环间存在等效结构误差角σ。首先推导只有单个磁环存在结构误差角的特殊情况时测角误差计算式,进而拓展到一般情况。如图2所示,设A环准确东西摆放,角度自西向东逆时针由–180°~180°递进。C环准确朝向120°,B环与理想方向60°之间的夹角σ即为等效结构误差角。

    图 2  单磁环结构误差角示意图

    定义雷电来向与正东方向实际夹角为θ,通过式(3)得到的计算角度为θ。根据几何关系,理想情况下B环感应电动势EB与实际情况下感应电动势EB的比值为

    EBEB=sin(60θ)sin(60(θσ))=K(σ,θ)
    (4)

    根据A环、B环电压关系可进一步求得计算角度θ1

    θ1=arctan(3EAEA2EB/2EBK(σ,θ)K(σ,θ))
    (5)

    代入EA消去EB即可化简得到计算角度θ1关于结构误差角σ和实际角度θ的函数关系

    θ1=arctan(3tanθ(1+3sinσcosσ)tanθ+3cosσ+sinσ)
    (6)

    同理可得θ3

    θ3=arctan(3(3sinσcosσ1)tanθ + 3cosσ + sinσ3(3sinσcosσ + 1)tanθ + 3cosσ + sinσ + 3)
    (7)

    由于C环不存在结构误差角,所以θ2=θ。定义计算角度θ与实际角度θ的差为测角误差δ,则有

    δ=θ1+θ32θ3
    (8)

    为更直观地体现结构误差角σ对测角误差δ的影响,仿真测角误差δ关于实际角度θ的变化情况,如图3所示。

    图 3  三磁环天线测角误差趋势图

    图3可知,在B环结构误差角度σ一定的情况下,测角误差δ随实际角度θ呈波浪形变化,当θ在60°和–120°附近时测角误差最大,最大值约等于0.7倍的结构误差角σ

    进一步讨论一般情况,当B环、C环均存在结构误差角时,如图4所示,B环、C环与60°,120°的夹角分别为σ1,σ2,规定当磁环沿逆时针偏离理想方向时,其结构误差角为正,反之为负。定义σ1,σ2中绝对值较大的为三磁环天线结构误差角σ。为便于分析,先设定σ1,σ2大小相等。

    图 4  两磁环结构误差角示意图

    类似单磁环存在结构误差角时求解方法,可分别求得两种情况时测角误差δ与三磁环天线结构误差角σ的函数关系,进一步仿真测角误差δ关于实际角度θ的变化情况,如图5所示。

    图 5  三磁环天线测角误差趋势图

    图5可知,当三磁环天线结构误差角σ一定时,测角误差δ总体呈波浪形变化趋势。分析图5(a)可知,当σ1, σ2同号时,测角误差最大值出现在90°和–90°附近,最大值约为0.75σ。当σ1, σ2大小不相等时,保持σ1值不变,改变σ2的大小在0~σ1之间变化,通过仿真可得,测角误差最大值在0.7σ1~0.75σ1之间变化,当σ2为0时,测角误差最大值为0.7σ1。因此在σ1, σ2同号的情况下,测角误差最大值为0.75σ

    分析图5(b)可知,当σ1, σ2异号时,测角误差最大值出现在±60°和±120°附近,最大值的绝对值约等于0.6σ。当σ1, σ2大小不相等时,保持σ1值不变,改变σ2的大小在0~σ1变化,通过仿真可得,测角误差最大值在0.6σ1~0.7σ1之间变化,当σ2为0时,测角误差最大值为0.7σ1。因此这种情况下三磁环天线测角误差最大值为0.7σ

    上述推导需满足在同等电磁场条件下,3个磁环的感应电动势峰值相同。然而在实际情况中由于制作工艺等原因会造成磁环的面积、阻抗不同,进而导致各磁环感应电动势峰值不同,由此引发的误差可称为一致性误差。一致性误差可通过引入补偿系数进行矫正,设三磁环感应电动势峰值分别为VA, VB, VC,选定A环感应电压值为标准值,对B环、C环电压值进行补偿,则B环补偿系数ηB可表示为

    ηB=VBVA
    (9)

    同理可得C环补偿系数ηC。在多个电磁场强度下,分别测量三磁环在相同角度时感应电动势,对磁环A,B和磁环A,C两组电压数据分别作拟合直线,两直线斜率即为B环和C环的补偿系数。通过上述方法能够基本矫正磁环的一致性误差。

    考虑正交磁环天线磁环间夹角90°与三磁环天线夹角60°比值为1.5,因此在同等结构误差角条件下,正交磁环天线结构误差角约为三磁环天线的1.5倍,即正交磁环天线结构误差角约为1.5σ。已知正交磁环天线测角误差最大值约等于其结构误差角[16],仿真结果如图6所示。综上所述,不考虑一致性误差的影响,理论上采用三磁环天线进行雷电测向,测角误差能够降低约50%。

    图 6  正交磁环天线测角误差趋势图

    以下通过对照实验对比三磁环天线与同尺寸正交磁环天线测角精度。实验选取半径为7.5 cm的磁环,根据关于雷电测向的正交磁环天线的接特性研究[18],可确定该磁环天线在最佳接收效益时的线圈匝数为60匝。设计参数如表1所示的三磁环天线,实物如图7所示。

    表 1  磁环天线参数
    磁环半径r(cm)线圈宽度d(cm)导线直径
    Φ(mm)
    磁环匝数
    (N)
    7.51.50.360
    下载: 导出CSV 
    | 显示表格
    图 7  三磁环天线实物图

    利用信号发生器产生200 kHz的正弦信号,示波器连接信号发射线圈,发射线圈匝数为30匝,半径为40 cm,面积远大于天线线圈。将天线固定在位于发射线圈正前方10 cm的转台上,通过转动天线来改变信号方向角。将天线3个磁环分别接入终端阻抗为1 MΩ、耦合方式为直流的MDO-3052型数字示波器的3个通道,记录天线各磁环的响应电压。

    首先对三磁环天线的一致性误差进行检验和矫正,利用发射线圈从0°~90°间隔10°分别对3个磁环进行响应测试并记录3组电压值数据,采用第3节所述方法可求得B环、C环的补偿系数ηB, ηC分别为1.03, 1.02。

    接着利用发射线圈从–180°~180°间隔10°对天线进行响应测试,将3个磁环的响应电压值通过补偿系数修正后代入式(3)求得方向角计算值θ,与实际方向角θ相减得到测角误差δ。用相同参数的磁环制作正交磁环天线,先对其一致性误差进行矫正,然后通过对照响应试验求得不同角度时正交磁环天线的测角误差δ。实验结果表明,相比于正交磁环天线,三磁环天线测角误差明显降低,误差最大值由3.3°下降为1.19°,降低65%,测角误差平均值从1.67°下降为0.49°,降低70.6%,标准差由1.01°下降为0.32°,降低68.3%。为更直观体现优化程度,将两种天线测角误差取模并作多项式拟合,得出曲线图如图8所示。

    图 8  两种天线测角误差对比图

    图8可知,相比于正交磁环天线,采用三磁环天线对信号源进行方位角探测的测角误差显著降低,整体误差降低约50%。进一步分析,正交磁环天线测角误差峰值分布在 ±90°附近,误差峰值约为3°;三磁环天线的误差角度峰值分布在 ±120°和±60°附近,误差峰值约为1°,结合第3节仿真可知三磁环天线两结构误差角σ1,σ2为异号,测角误差变化趋势与仿真结果基本相符。

    本文分析了正交磁环天线由于结构误差角造成的测角误差较大的问题,提出一种改进的三磁环天线,该天线通过两两呈60°的磁环进行测向。通过建立测角误差仿真模型分析了三磁环天线测角误差优化效果,对比实验表明,在同等制作工艺下,相较于正交磁环天线,三磁环天线测角误差整体降低约50%,提高了测向精度。

  • 图  1  RSSID定位原理图

    图  2  随信噪比变化的定位性能

    图  3  随监测站间距变化的定位性能

    图  4  卡尔曼滤波效果图

    图  5  单目标定位误差对比

    图  6  多目标定位平均误差对比

    表  1  单目标定位10000次误差统计分析(km)

    是否预处理定位方法最大误差最小误差平均误差
    RSSI定位5.08020.02391.3993
    RSSID定位4.62240.00760.8527
    RSSI定位1.45370.22730.6249
    RSSID定位0.88010.00680.2683
    下载: 导出CSV

    表  2  多目标定位10000次平均误差统计分析(km)

    是否预处理定位方法最大误差最小误差平均误差
    RSSI定位1.86021.35991.5911
    RSSID定位1.10150.76200.9170
    RSSI定位1.53120.94701.1530
    RSSID定位0.82840.19480.2930
    下载: 导出CSV
  • 王先义, 陈丹俊, 刘斌, 等. 复杂电磁环境战场频谱管理[J]. 中国电子科学研究院学报, 2008, 3(4): 338–344. doi: 10.3969/j.issn.1673-5692.2008.04.002

    WANG Xianyi, CHEN Danjun, LIU Bin, et al. Battlefield spectrum management in complex electromagnetic environment[J]. Journal of China Academy of Electronics and Information Technology, 2008, 3(4): 338–344. doi: 10.3969/j.issn.1673-5692.2008.04.002
    史小红. 基于TDOA 的无线定位方法及其性能分析[J]. 东南大学学报(自然科学版), 2013, 43(2): 252–257. doi: 10.3969/j.issn.1001-0505.2013.02.005

    SHI Xiaohong. TDOA-based wireless positioning method and its performance[J]. Journal of Southeast University (Natural Science Edition), 2013, 43(2): 252–257. doi: 10.3969/j.issn.1001-0505.2013.02.005
    WARNER G A and DOSSO S E. Bowhead whale localization using time-difference-of-arrival data from asynchronous recorders[J]. Journal of the Acoustical Society of America, 2017, 141(3): 1921. doi: 10.1121/1.4978438
    WEI Junming and XU Changbin. Improvement of software defined radio based TDOA source localization[C]. Annual Conference of the IEEE Industrial Electronics Society, Dallas, USA, 2015: 5307–5313.
    LIU Congfeng, YANG Jie, and WANG Fengshuai. Joint TDOA and AOA location algorithm[J]. Journal of Systems Engineering & Electronics, 2013, 24(2): 183–188. doi: 10.1109/JSEE.2013.00023
    NOROOZI A. Algebraic solution of source location estimation using TDOA and AOA measurements[C]. Iranian Conference on Electrical Engineering (ICEE), Tehran, Iran, 2017: 1609–1614.
    曲付勇, 孟祥伟. 基于约束总体最小二乘方法的到达时差到达频差无源定位算法[J]. 电子与信息学报, 2014, 36(5): 1075–1081. doi: 10.3724/SP.J.1146.2013.01019

    QU Fuyong and MENG Xiangwei. Source localization using TDOA and FDOA measurements based on constrained total least squares algorithm[J]. Journal of Electronics &Information Technology, 2014, 36(5): 1075–1081. doi: 10.3724/SP.J.1146.2013.01019
    LU Zhiyu, BA Bin, WANG Jianhui, et al. A direct position determination method with combined TDOA and FDOA based on particle filter[J]. Chinese Journal of Aeronautics, 2018, 31(1): 161–168. doi: 10.1016/j.cja.2017.11.007
    陈淑敏, 乔晓田, 毛佳, 等. 基于接收信号强度(RSSI) 的室内二次定位方法[J]. 传感技术学报, 2015, 28(4): 572–577. doi: 10.3969/j.issn.1004-1699.2015.04.020

    CHEN Shumin, QIAO Xiaotian, MAO Jia, et al. RSSI-based twice location method for indoor positioning[J]. Chinese Journal of Sensors and Actuators, 2015, 28(4): 572–577. doi: 10.3969/j.issn.1004-1699.2015.04.020
    ZHENG Jungang, WU Chengdong, HAO Chu, et al. Localization algorithm based on RSSI and distance geometry constrain for wireless sensor network[C]. IEEE International Conference on Electrical and Control Engineering, Wuhan, China, 2010: 2836–2839.
    TAZAWA R, HONMA N, MIURA A, et al. RSSI-based localization using wireless beacon with three-element array[J]. IEICE Transactions on Communications, 2018, 101(2). doi: 10.1587/transcom.2017ISP0027
    MARC A K J and KAZUNORI O. LRD: A distributed and accurate localization technique for wireless sensors networks[C]. TENCON 2010-2010 IEEE Region 10 Conference, Fukuoka, Japan, 2010: 234–239.
    谢树果, 叶知秋, 王磊, 等. 一种基于场强测量的辐射源定位方法[J]. 电波科学学报, 2013, 28(6): 1133–1138. doi: 10.13443/j.cjors.2013.06.024

    XIE Shuguo, YE Zhiqiu, WANG Lei, et al. Radio source locating method based on signal strength measurement[J]. Chinese Journal of Radio Science, 2013, 28(6): 1133–1138. doi: 10.13443/j.cjors.2013.06.024
    卞国龙, 黄海松, 王安忆, 等. 无线传感器网络定位技术的研究[J]. 实验室研究与探索, 2017, 36(6): 122–127. doi: 10.3969/j.issn.1006-7167.2017.06.029

    BIAN Guolong, HUANG Haisong, WANG Anyi, et al. Research on localization technology of wireless sensor networks[J]. Research and Exploration in Laboratory, 2017, 36(6): 122–127. doi: 10.3969/j.issn.1006-7167.2017.06.029
    JUNG J Y, KANG D O, CHOI J H, et al. D2D distance measurement using Kalman filter algorithm for distance-based service in an office environment[C]. International Conference on Advanced Communication Technology (ICACT), Seoul, South Korea, 2015: 221–224.
    ZHANG Kaibi, ZHANG Yangchuan, and WAN Subo. Research of RSSI indoor ranging algorithm based on Gaussian - Kalman linear filtering[C]. IEEE Advanced Information Management, Communicates, Electronic and Automation Control Conference (IMCEC), Xi'an, China, 2016: 1628–1632.
    YU Zongzuo and GUO Gaizhi. Improvement of positioning technology based on RSSI in ZigBee networks[J]. Wireless Personal Communications, 2016, 95(3): 1–20. doi: 10.1007/s11277-016-3860-1
    ZHOU Cheng, YUAN Jiazheng, LIU Hongzhe, et al. Bluetooth indoor positioning based on RSSI and Kalman filter[J]. Wireless Personal Communications, 2017, 96(3): 1–16. doi: 10.1007/s11277-017-4371-4
    方震, 赵湛, 郭鹏, 等. 基于RSSI测距分析[J]. 传感技术学报, 2007, 20(11): 2526–2530. doi: 10.3969/j.issn.1004-1699.2007.11.037

    FANG Zhen, ZHAO Zhan, GUO Peng, et al. Analysis of distance measurement based on RSSI[J]. Chinese Journal of Sensors and Actuators, 2007, 20(11): 2526–2530. doi: 10.3969/j.issn.1004-1699.2007.11.037
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  1986
  • HTML全文浏览量:  955
  • PDF下载量:  101
  • 被引次数: 1
出版历程
  • 收稿日期:  2018-03-21
  • 修回日期:  2018-08-13
  • 网络出版日期:  2018-08-31
  • 刊出日期:  2019-02-01

目录

/

返回文章
返回