CN108152789A - 利用rss信息的无源多站多目标数据关联与定位方法 - Google Patents
利用rss信息的无源多站多目标数据关联与定位方法 Download PDFInfo
- Publication number
- CN108152789A CN108152789A CN201810003502.8A CN201810003502A CN108152789A CN 108152789 A CN108152789 A CN 108152789A CN 201810003502 A CN201810003502 A CN 201810003502A CN 108152789 A CN108152789 A CN 108152789A
- Authority
- CN
- China
- Prior art keywords
- data
- target
- point
- receiving station
- corelation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于电子对抗技术领域,具体的说是一种利用RSS信息的无源多站多目标数据关联与定位方法。本发明在是在利用多个接收站接收的来自不同辐射源发射信号的信号强度,对多个站的多个接收数据进行关联时,针对出现虚假目标和真实目标交叉点,所提出的解决方法:首先采用多目标数据关联算法,将多个接收站的多个信号接受强度数据进行关联,根据多站多目标数据的关联,对于所得到的代价矩阵,采用自适应k‑medoids聚类算法估计出目标坐标。
Description
技术领域
本发明属于电子对抗技术领域,具体涉及一种无源定位中基于接收信号强度(Received Signal Strength,RSS)的多传感器多目标数据融合方法。
背景技术
电子侦察过程,准确估计目标辐射源位置有助于获取辐射源信息,是做好高层次上的态势估计和威胁估计的关键和主要依据,也是对目标实现精准打击的重要保证。在无源定位系统中,其中基于RSS信息的定位方法由于对传输带宽要求低,无需额外硬件的优点得到广泛应用。基于RSS目标定位方法利用了接收信号强度和路径衰减的关系,可以实现对单目标的定位。但针对多站多目标定位时该方法会面临多数据关联问题和如何去除多站数据关联产生的鬼影点即虚假目标问题。同时,鬼影点的个数会随着目标个数的增加而呈现出几何级数的增长如何快速高效的剔除虚假交叉点。诸如此类的问题,使得获得的目标数目估计和位置的估计之间存在较大的偏差,因此需要有新的方法,针对性的解决这类问题。
数据聚类是在无标记样本的条件下将数据分组,从而发现数据的天然结构,在数据分析中扮演了重要的角色。基于划分的聚类方法是一种被广泛研究和应用的数据聚类方法,其中k-means算法被选为数据挖掘领域的十大算法之一。但是在实际情况下,k-means由于需要诸如类别个数等先验信息、对离群点敏感等问题而有所缺点,于是国内外学者研究出类k-means算法,其中k-medoids算法是一种比k-means算法能更好处理离群点的算法,同时,于k-means算法不同是,k-medoids算法只选取真实的数据点作为中心点,而k-means算法中的中心点被设定为当前簇中所有数据点的平均值,是一个虚拟点,相比簇均值,真实的数据点受到离群点的影响会更小一点。以k-medoids算法为代表的聚类算法在无源定位中的研究目前还很少,是值得研究的一个方向。
发明内容
本发明的目的,就是针对上述问题,提出了一种结合多辐射源数据关联和自适应k-medoids两种算法相结合来实现基于RSS信息的无源多站多目标数据融合,从而完成目标位置估计。。
本发明所采用的技术方案为:
基于RSS信息的无源多站多目标数据融合,所述的是在利用多个接收站接收的来自不同辐射源发射信号的信号强度,对多个站的多个接收数据进行关联时,针对出现虚假目标和真实目标交叉点,所提出的解决方法。基于RSS信息的无源多站多目标定位示意图如图1所示,为简单起见,图1只给出了两个接收站在对两个目标进行交叉定位时的鬼影位置和真实目标所在的位置。
基于RSS信息实现定位的原理是:以目标1为例,设目标1的坐标为(xm,yn),距接收站1、2的距离分别为接收站1、2接收目标1发射信号的接收功率分别为 则有经典的强度-距离损耗公式:
Pr=P0-10βlg(d/d0)
计算可得到其中,Pr表示接收强度,P0表示在参考传输距离为d0处的接收强度,d为信号传输距离。β为路径衰减指数,一般介于2至5之间,自由空间通常取为2。再由公式:
计算可得到目标1的坐标(xm,yn)。以上是针对单个目标定位的解决方法,但是当涉及多个接收站对对多个目标定位时,还需要考虑多个目标的数据的关联问题,即不同接收站的哪些数据是来自同一个目标,多数据关联即会出现多个鬼影点,即虚假目标。
鬼影点是由两接收站分别对两个不同的真实目标信号接收强度的错误关联组合得到的,所以鬼影点的个数会随着目标个数的增加而呈现出几何级数的增长,当真实目标个数为n时,鬼影点的个数则为n2-n。要剔除鬼影点,减少鬼影点对真实目标定位的影响,首先要做的是多站多目标数据关联,然后在进行多目标定位,所以发明如下算法:
首先采用多目标数据关联算法,将多个接收站的多个信号接受强度数据进行关联,算法路程图如图2所示,具体做法如下:
S1、假设目标和接收站都位于XOY平面上,设接收站的位置坐标为(xk,yk),其中k=1,2,...K,K是接收站总数,每个接收站所接收的信号强度为:
其中j=1,2...J,J是目标总数,表示测量噪声。
S2、将目标平面划分为M×N范围的网格,每个格点代表目标平面中的一个位置坐标((xm,yn),其中m=1,2,...M,n=1,2,...N,遍历网格平面中的每个网格点,假设每个点上都有一个目标发出信号,计算每个接收站接收到的每个点上的接收强度:
S3、计算每个搜索点(xm,yn)相对于接收站的RSS值接收站所接收到的真实RSS值之间的误差
其中,k=1,2,...K,K是接收站总数,j=1,2...J,J是目标总数,m=1,2,...M,n=1,2,...N。
S4、计算由每次搜索点与真实值之间的总误差,组成代价矩阵C(m,n):
其中,k=1,2,...K,K是接收站总数,m=1,2,...M,n=1,2,...N。
S5、S1~S4即完成了多站多目标数据的关联,对于所得到的代价矩阵C(m,n),采用自适应k-medoids聚类算法估计出目标坐标。
进一步的,所述执行步骤S5中自适应k-medoids聚类算法流程图如图3所示,具体如下:
S6、数据预处理:将代价矩阵中大于阈值的数据剔除,把剩下所有数据在矩阵中所对应的索引作为元素组成新的包含n个对象的数据库X;
S7、自适应聚类个数为J(目标总数),随机选择J个代表对象作为初始中心点;
S8、将数据库中其余每个非代表对象指派到距离其最近的代表对象所属的簇,即对于数据集中的其他每一个点,计算它与每个聚类中心(指已选择的初始中心)的距离D(xj),j=1,2...J,将其分配到最小的D(xj)值对应的代表对象所属的簇。
S9、剔除步骤S8中空簇集合;
S10、在每个簇集合中随机选取一个非代表对象Or,计算用Or替代其所属簇中的代表对象Oj所产生的代价S,即计算各簇集合中,距簇内各样本点距离的总的绝度误差构成总代价S;
S11、如果Or对应的S最小,则由Or替代Oj构成新的簇集合,作为新的中心点;
S12、重复S7~S11,直到形成的簇集合的中心点不再发生变化。
本发明的有益效果为,本发明可以准确对多站数据完成关联,并最终准确估计出目标的位置,方法简单,效果良好。
附图说明
图1为基于RSS信息多站多目标定位示意图;
图2为多目标数据关联算法流程图;
图3为自适应k-medoids聚类算法流程图;
图4为定位模型图;
图5为多目标数据关联后的目标分布图;
图6为理想情况下定位正确率随门限大小变化曲线图;
图7为定位正确率随误差变化曲线图;
图8为定位性能RMSE值随误差变化曲线图;
图9为定位正确率随目标间距变化曲线图;
图10为定位性能RMSE值随目标间距变化曲线图。
具体实施方式
下面结合附图对上述基于RSS信息的无源多站多目标定位算法方案进行验证,为简化起见,对算法模型作如下假设:
1.接收站和目标都在XOY平面内;
2.所有的工程误差都叠加到RSS值测量误差中;
3.假设目标静止或运动速度极低;
4.P0=1W,β=2,d0=1m;
5.搜索间隔为0.5km;
如图4所示,假设目标区域为50km×50km的方形区域,利用4个固定接收站对4个目标进行数据关联并且定位,接收站采用最优布局,即在目标区域的4个方向,坐标分别为(0,0),(0,50),(50,0),(50,50),单位为km,目标平行与X轴分布,目标位置为[20+(i-1)*R,30],i=1,2,3,4,单位为km,其中R为目标间距。所有接收站的测量误差服从均值为零的高斯分布,且误差之间相互独立。1000次Monte Carlo仿真后得出多站数据关联效果和定位效果。
多站多目标数据关联效果:
如图5所示,对多个站所观测到的关于各个目标的RSS信息进行关联后,从图中可以看出有4个明显的谱峰,可以认为这4个明显的谱峰对应的位置正是4个目标位置,由此可以知算法能够准确的将各个接收站接收的信号进行正确关联,从而能够减少鬼影点对后续定位的影响。
自适应k-medoids聚类算法定位效果:
仿真情况1:目标间距为5km,门限为0:0.5:2.5,单位km,测量误差为0.1dbm;
研究理想情况时,该算法在不同门限值下的定位正确率。由于在观测过程中始终存在测量误差,因此当目标估计位置与目标真实位置之间的距离小于某一门限时,可以认为已经获得正确的定位结果。假设目标呈分散分布,在RSS值测量误差非常小几乎可以忽略不计,由图6可以看出,在理想情况下,即目标间距5km,测量误差为0.1dbm时,该算法在估计位置和真实位置之间的距离差为0.5km时就可以达到很高的正确率。
仿真情况2:目标间距R为5km,门限为2km,测量误差为0~20,单位dbm;
在此情况下,假设真实目标位置与估计目标位置之间的距离差小于2km时即表示目标位置估计正确,依次增加各个站测量误差,得到该算法下测量误差对定位正确率的影响。由图7可以看出,当误差逐渐增大时,定位正确率也逐渐下降,当测量误差在0~10dbm之间时,正确率曲线下降速度快,10dbm以后曲线变化平缓。同时由图8可以看出,随着误差的增加,其对应的平均距离相对误差逐渐增加,误差在0~10dbm曲线呈陡峭上升趋势,10dbm以后曲线变化缓慢。图7和图8说明在误差允许的范围内,该算法对实现多站多目标测量数据关联以及定位可行。
仿真情况3:目标间距R为1~5km,门限为0.5km,测量误差为0.1dbm;
在此情况下,假设真实目标位置与估计目标位置之间的距离小于0.5km时即表示目标位置估计正确,依次增加目标之间的距离,得到该算法下目标间距对定位正确率的影响。由图9可以看出,当目标间距逐渐增大时,定位正确率也逐渐增加,当目标间距为4km时定位正确率极高,同时图10可以看出平均距离相对误差随着间距的增加逐渐减小。说明在误差允许的范围内,该算法对实现多站多目标测量数据关联以及定位可行。
Claims (1)
1.利用RSS信息的无源多站多目标数据关联与定位方法,其特征在于,包括:
采用多目标数据关联算法,将多个接收站的多个信号接受强度数据进行关联,具体为:
S1、假设目标和接收站都位于XOY平面上,设接收站的位置坐标为(xk,yk),其中k=1,2,...K,K是接收站总数,每个接收站所接收的信号强度为:
其中j=1,2...J,J是目标总数,表示测量噪声;
S2、将目标平面划分为M×N范围的网格,每个格点代表目标平面中的一个位置坐标((xm,yn),其中m=1,2,...M,n=1,2,...N,遍历网格平面中的每个网格点,假设每个点上都有一个目标发出信号,计算每个接收站接收到的每个点上的接收强度:
S3、计算每个搜索点(xm,yn)相对于接收站的RSS值和接收站所接收到的真实RSS值之间的误差
S4、计算由每次搜索点与真实值之间的总误差,组成代价矩阵C(m,n):
S5、S1~S4即完成了多站多目标数据的关联,对于所得到的代价矩阵C(m,n),采用自适应k-medoids聚类算法估计出目标坐标,具体包括:
S51、数据预处理:将代价矩阵中大于阈值的数据剔除,把剩下所有数据在矩阵中所对应的索引作为元素组成新的包含n个对象的数据库X;
S52、自适应聚类个数为J,即目标总数,随机选择J个代表对象作为初始中心点;
S53、将数据库中其余每个非代表对象指派到距离其最近的代表对象所属的簇,即对于数据集中的其他每一个点,计算它与每个聚类中心(指已选择的初始中心)的距离D(xj),j=1,2...J,将其分配到最小的D(xj)值对应的代表对象所属的簇。
S54、剔除步骤S53中空簇集合;
S55、在每个簇集合中随机选取一个非代表对象Or,计算用Or替代其所属簇中的代表对象Oj所产生的代价S,即计算各簇集合中,距簇内各样本点距离的总的绝度误差构成总代价S;
S56、如果Or对应的S最小,则由Or替代Oj构成新的簇集合,作为新的中心点;
S57、重复S51~S56,直到形成的簇集合的中心点不再发生变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810003502.8A CN108152789B (zh) | 2018-01-03 | 2018-01-03 | 利用rss信息的无源多站多目标数据关联与定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810003502.8A CN108152789B (zh) | 2018-01-03 | 2018-01-03 | 利用rss信息的无源多站多目标数据关联与定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108152789A true CN108152789A (zh) | 2018-06-12 |
CN108152789B CN108152789B (zh) | 2019-11-01 |
Family
ID=62460522
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810003502.8A Active CN108152789B (zh) | 2018-01-03 | 2018-01-03 | 利用rss信息的无源多站多目标数据关联与定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108152789B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110412504A (zh) * | 2019-08-12 | 2019-11-05 | 电子科技大学 | 基于角度与时差信息的无源多站多目标关联与定位方法 |
CN111079859A (zh) * | 2019-12-31 | 2020-04-28 | 哈尔滨工程大学 | 一种无源多站多目标测向交叉定位与虚假点去除方法 |
CN111273228A (zh) * | 2019-05-09 | 2020-06-12 | 哈尔滨工程大学 | 一种基于遍历搜索策略的多相干源定位方法 |
CN111352087A (zh) * | 2020-03-25 | 2020-06-30 | 电子科技大学 | 基于dbscan的被动mimo雷达多目标定位方法 |
CN111398896A (zh) * | 2020-04-02 | 2020-07-10 | 哈尔滨工程大学 | 基于基准线成对交叉定位点距离的虚假点去除预处理方法 |
CN113390406A (zh) * | 2021-06-16 | 2021-09-14 | 电子科技大学 | 基于被动多传感器系统的多目标数据关联与定位方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102075686A (zh) * | 2011-02-10 | 2011-05-25 | 北京航空航天大学 | 一种鲁棒的实时在线摄像机跟踪方法 |
CN102761913A (zh) * | 2011-04-26 | 2012-10-31 | 航天信息股份有限公司 | 基于区域划分的无线信号传输参数确定的定位方法 |
CN103854014A (zh) * | 2014-02-25 | 2014-06-11 | 中国科学院自动化研究所 | 一种基于上下文稀疏表示的恐怖视频识别方法及装置 |
CN103985120A (zh) * | 2014-05-13 | 2014-08-13 | 西北工业大学 | 一种遥感图像多目标关联的方法 |
CN104363653A (zh) * | 2014-10-13 | 2015-02-18 | 西北大学 | 一种消除环境噪声的被动式定位方法 |
CN105184314A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于像素聚类的wrapper式高光谱波段选择方法 |
CN107436427A (zh) * | 2017-01-12 | 2017-12-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 空间目标运动轨迹与辐射信号关联方法 |
-
2018
- 2018-01-03 CN CN201810003502.8A patent/CN108152789B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102075686A (zh) * | 2011-02-10 | 2011-05-25 | 北京航空航天大学 | 一种鲁棒的实时在线摄像机跟踪方法 |
CN102761913A (zh) * | 2011-04-26 | 2012-10-31 | 航天信息股份有限公司 | 基于区域划分的无线信号传输参数确定的定位方法 |
CN103854014A (zh) * | 2014-02-25 | 2014-06-11 | 中国科学院自动化研究所 | 一种基于上下文稀疏表示的恐怖视频识别方法及装置 |
CN103985120A (zh) * | 2014-05-13 | 2014-08-13 | 西北工业大学 | 一种遥感图像多目标关联的方法 |
CN104363653A (zh) * | 2014-10-13 | 2015-02-18 | 西北大学 | 一种消除环境噪声的被动式定位方法 |
CN105184314A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于像素聚类的wrapper式高光谱波段选择方法 |
CN107436427A (zh) * | 2017-01-12 | 2017-12-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 空间目标运动轨迹与辐射信号关联方法 |
Non-Patent Citations (2)
Title |
---|
SREE DIVYA CHITTE ET AL.: "Distance Estimation From Received Signal Strength Under Log-Normal Shadowing: Bias and Variance", 《IEEE SIGNAL PROCESSING LETTERS》 * |
李万春: "加权多维标量的接收信号强度定位方法", 《信号处理》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111273228A (zh) * | 2019-05-09 | 2020-06-12 | 哈尔滨工程大学 | 一种基于遍历搜索策略的多相干源定位方法 |
CN110412504A (zh) * | 2019-08-12 | 2019-11-05 | 电子科技大学 | 基于角度与时差信息的无源多站多目标关联与定位方法 |
CN111079859A (zh) * | 2019-12-31 | 2020-04-28 | 哈尔滨工程大学 | 一种无源多站多目标测向交叉定位与虚假点去除方法 |
CN111079859B (zh) * | 2019-12-31 | 2020-12-04 | 哈尔滨工程大学 | 一种无源多站多目标测向交叉定位与虚假点去除方法 |
CN111352087A (zh) * | 2020-03-25 | 2020-06-30 | 电子科技大学 | 基于dbscan的被动mimo雷达多目标定位方法 |
CN111352087B (zh) * | 2020-03-25 | 2022-05-24 | 电子科技大学 | 基于dbscan的被动mimo雷达多目标定位方法 |
CN111398896A (zh) * | 2020-04-02 | 2020-07-10 | 哈尔滨工程大学 | 基于基准线成对交叉定位点距离的虚假点去除预处理方法 |
CN111398896B (zh) * | 2020-04-02 | 2020-12-22 | 哈尔滨工程大学 | 基于基准线成对交叉定位点距离的虚假点去除预处理方法 |
CN113390406A (zh) * | 2021-06-16 | 2021-09-14 | 电子科技大学 | 基于被动多传感器系统的多目标数据关联与定位方法 |
CN113390406B (zh) * | 2021-06-16 | 2022-05-24 | 电子科技大学 | 基于被动多传感器系统的多目标数据关联与定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108152789B (zh) | 2019-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108152789B (zh) | 利用rss信息的无源多站多目标数据关联与定位方法 | |
CN108089148B (zh) | 一种基于时差信息的无源多站多目标测向交叉定位方法 | |
CN108061877B (zh) | 一种基于角度信息的无源多站多目标测向交叉定位方法 | |
CN109901153A (zh) | 基于信息熵权和最近邻域数据关联的目标航迹优化方法 | |
CN103648164B (zh) | 一种基于到达时间差和Gossip算法的无线传感器网络分布式定位方法 | |
CN111479231A (zh) | 一种用于毫米波大规模mimo系统的室内指纹定位方法 | |
CN102325370A (zh) | 一种高精度的无线传感器网络节点的三维定位器 | |
CN111460597B (zh) | 一种基于改进的多目标粒子群优化算法的雷达布站方法 | |
CN108802674A (zh) | 一种针对直接定位的联合搜索方法及装置 | |
Will et al. | The membership degree min-max localization algorithm | |
Fan et al. | A structured bidirectional LSTM deep learning method for 3D terahertz indoor localization | |
CN104185272A (zh) | 基于优化平均每跳距离和优选锚节点的wsn定位方法 | |
Wu et al. | TDOA based indoor positioning with NLOS identification by machine learning | |
CN107063259A (zh) | 一种航迹关联方法及电子设备 | |
CN110049549A (zh) | 基于WiFi指纹的多融合室内定位方法及其系统 | |
CN106934324A (zh) | 基于简化多假设算法的雷达数据关联方法 | |
CN102288938B (zh) | 一种有效的无线传感器网络节点的三维定位方法 | |
Huang et al. | Analysis of TOA localization with heteroscedastic noises | |
WO2019148801A1 (zh) | 一种wsn节点的定位方法及装置 | |
CN110673090A (zh) | 基于dbscan的无源多站多目标定位方法 | |
CN108737952A (zh) | 基于rssi测距改进的多边形加权质心定位方法 | |
CN110412504A (zh) | 基于角度与时差信息的无源多站多目标关联与定位方法 | |
CN101308206B (zh) | 一种白噪声背景下的圆周轨迹机动目标跟踪方法 | |
CN115952691A (zh) | 多站无源时差交叉联合定位系统的优化布站方法及装置 | |
CN104469939B (zh) | 基于rss统计分布分块区域误差限的wlan定位网络优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |