CN109218961B - Multi-station cooperative interference positioning method and system based on virtual nodes - Google Patents

Multi-station cooperative interference positioning method and system based on virtual nodes Download PDF

Info

Publication number
CN109218961B
CN109218961B CN201810870985.1A CN201810870985A CN109218961B CN 109218961 B CN109218961 B CN 109218961B CN 201810870985 A CN201810870985 A CN 201810870985A CN 109218961 B CN109218961 B CN 109218961B
Authority
CN
China
Prior art keywords
node
positioning
coordinates
signal source
virtual
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.)
Active
Application number
CN201810870985.1A
Other languages
Chinese (zh)
Other versions
CN109218961A (en
Inventor
姚继明
郭经红
韦磊
卜宪德
吴鹏
刘世栋
张�浩
张刚
王玮
郭云飞
朱道华
陶静
吕立冬
喻强
刘川
孙云晓
郭雅娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
State Grid Corp of China SGCC
Original Assignee
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
State Grid Corp of China SGCC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Jiangsu Electric Power Co Ltd, Global Energy Interconnection Research Institute Co Ltd, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd, State Grid Corp of China SGCC filed Critical State Grid Jiangsu Electric Power Co Ltd
Priority to CN201810870985.1A priority Critical patent/CN109218961B/en
Publication of CN109218961A publication Critical patent/CN109218961A/en
Application granted granted Critical
Publication of CN109218961B publication Critical patent/CN109218961B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/023Services making use of location information using mutual or relative location information between multiple location based services [LBS] targets or of distance thresholds
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • H04W64/006Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses a multi-station cooperative interference positioning method and a system based on virtual nodes.A first three nodes of RSSI (received signal strength indicator) values of signals of a signal source to be positioned received by a monitoring node are used as positioning nodes according to the sequence of the RSSI values, the positioning node with the maximum RSSI value is used as a reference node, and the coordinates of a second/third virtual node are determined according to the coordinates of the second/third node, the transmission time difference between the second/third node and the reference node to reach the signal source to be positioned and the direction of the signals received by the second/third node; and constructing a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node, calculating initial positioning coordinates of each combination, and calculating the positioning coordinates of the signal source to be positioned according to the initial positioning coordinates and the corresponding weighting coefficients. By means of constructing the virtual nodes, node combinations of cooperative positioning are enriched, positioning errors caused by node distribution are reduced, and positioning accuracy is improved by means of weighting of positioning results of different cooperative combinations based on RSSI values.

Description

一种基于虚拟节点的多站协同干扰定位方法及系统A method and system for multi-station cooperative interference localization based on virtual nodes

技术领域technical field

本发明涉及无线通信领域,具体涉及一种基于虚拟节点的多站协同干 扰定位方法及系统。The present invention relates to the field of wireless communication, in particular to a method and system for multi-station cooperative interference positioning based on virtual nodes.

背景技术Background technique

无源定位采用的基本技术主要有时差(TDOA)定位技术、测角(AOA) 交叉定位技术、多普勒频差(FDOA)定位技术等。TDOA定位技术由于 可以获得较高的定位精度,同时对接收系统的要求较低,易于组网,因而 应用更加广泛。TDOA无源定位系统通常是由多个监测传感器组成的分 布式定位系统,多个监测传感器之间协同工作,在定位精度以及定位性能 方面相比于单站监测更优。TDOA定位计算过程主要分为两个部分,时 差计算和定位结果计算。时差计算主要考虑根据接收信号计算信源到两个 监测传感器的时间差;定位结果计算主要考虑根据时间差建立定位方程并 求解定位结果。The basic technologies used in passive positioning mainly include time difference (TDOA) positioning technology, angle measurement (AOA) cross positioning technology, and Doppler frequency difference (FDOA) positioning technology. TDOA positioning technology is more widely used because it can obtain higher positioning accuracy, at the same time, it has lower requirements on the receiving system and is easy to form a network. The TDOA passive positioning system is usually a distributed positioning system composed of multiple monitoring sensors, and the multiple monitoring sensors work together, which is better than single-station monitoring in terms of positioning accuracy and positioning performance. The TDOA positioning calculation process is mainly divided into two parts, time difference calculation and positioning result calculation. The time difference calculation mainly considers the time difference between the signal source and the two monitoring sensors according to the received signal; the positioning result calculation mainly considers establishing the positioning equation and solving the positioning result according to the time difference.

但是现有的定位算法主要存在以下问题:没有给出参考节点的选择依 据;定位结果的计算没有考虑协作节点接收到的RSSI(Received Signal Strength Indication,即接收的信号强度指示的)影响;当前算法的协作节 点选择都是基于实际部署的节点坐标,但在实际部署中,监测装置节点的 部署是不成规则的,使得定位结果精度不高。However, the existing positioning algorithms mainly have the following problems: the selection basis of the reference node is not given; the calculation of the positioning result does not consider the influence of the RSSI (Received Signal Strength Indication, that is, the received signal strength indication) received by the cooperative nodes; the current algorithm The selection of cooperative nodes is based on the node coordinates of the actual deployment, but in the actual deployment, the deployment of monitoring device nodes is irregular, so that the accuracy of the positioning results is not high.

发明内容SUMMARY OF THE INVENTION

因此,本发明提供一种基于虚拟节点的多站协同干扰定位方法及系统, 克服了现有技术中基于实际部署的节点坐标定位精度不高的问题。Therefore, the present invention provides a virtual node-based multi-station cooperative interference positioning method and system, which overcomes the problem of low accuracy of node coordinate positioning based on actual deployment in the prior art.

本发明实施例提供的一种基于虚拟节点的多站协同干扰定位方法,包括 如下步骤:获取监测节点中三个节点作为定位节点,分别为参考节点、第 二节点及第三节点;计算所述第二节点与所述参考节点到达待定位信号源 的传输时差,第三节点与所述参考节点到达待定位信号源的传输时差;根 据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源 的传输时差及所述第二节点接收信号的方向,确定第二虚拟节点的坐标; 根据所述第三节点的坐标、所述第三节点与所述参考节点到达待定位信号 源的传输时差及所述第三节点接收信号的方向,确定第三虚拟节点的坐标; 根据所述参考节点、所述第二节点与第二虚拟节点中之一、第三节点与第 三虚拟节点中之一构建定位组合;计算各所述定位组合对应的待定位信号 源的初始定位坐标;根据各定位节点接收待定位信号源信号的RSSI值分别 计算各所述定位组合对应初始定位坐标的加权系数;根据所述初始定位坐 标及与其对应的加权系数,计算所述待定位信号源的定位坐标。A method for locating multi-station cooperative interference based on virtual nodes provided by an embodiment of the present invention includes the following steps: acquiring three nodes in a monitoring node as positioning nodes, which are a reference node, a second node, and a third node; calculating the The transmission time difference between the second node and the reference node reaching the signal source to be located, the transmission time difference between the third node and the reference node reaching the signal source to be located; according to the coordinates of the second node, the second node and the Determine the coordinates of the second virtual node according to the transmission time difference of the reference node reaching the signal source to be located and the direction of the signal received by the second node; according to the coordinates of the third node, the arrival of the third node and the reference node The transmission time difference of the signal source to be located and the direction in which the third node receives the signal, determine the coordinates of the third virtual node; according to the reference node, one of the second node and the second virtual node, the third node and the One of the third virtual nodes constructs a positioning combination; calculates the initial positioning coordinates of the to-be-located signal source corresponding to each positioning combination; calculates the corresponding initial position of each positioning combination according to the RSSI value of the to-be-located signal source signal received by each positioning node Weighting coefficients of the positioning coordinates; calculating the positioning coordinates of the signal source to be positioned according to the initial positioning coordinates and the corresponding weighting coefficients.

优选地,所述获取监测节点中三个节点作为定位节点的步骤,包括: 获取所述监测节点中接收待定位信号源信号的RSSI值按高低排序的前三个 节点按高低排序的前三个节点作为所述定位节点;分别获取待定位信号源 信号到达各所述定位节点的RSSI值;将所述RSSI值最大的定位节点作为 所述参考节点;将另两个定位节点分别作为所述第二节点、第三节点。Preferably, the step of acquiring three nodes in the monitoring nodes as positioning nodes includes: acquiring the first three nodes in the monitoring node that receive RSSI values of the signal source signals to be located in order of high and low order The node is used as the positioning node; the RSSI value of the signal source signal to be located to each of the positioning nodes is obtained respectively; the positioning node with the largest RSSI value is used as the reference node; the other two positioning nodes are respectively used as the first positioning node. The second node and the third node.

优选地,根据所述第二节点的坐标、所述第二节点与所述参考节点到 达待定位信号源的传输时差及所述第二节点接收信号的方向,计算第二虚 拟节点的坐标的步骤,包括:根据所述第二节点与所述参考节点到达待定 位信号源的传输时差,计算所述第二节点和所述参考节点之间的传输距离 的差值;根据所述第二节点的坐标、所述第二节点和所述参考节点之间的 传输距离的差值及所述第二节点接收信号的方向,确定所述第二虚拟节点 的坐标。Preferably, the step of calculating the coordinates of the second virtual node is based on the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction in which the second node receives the signal , comprising: calculating the difference in the transmission distance between the second node and the reference node according to the transmission time difference between the second node and the reference node reaching the signal source to be located; The coordinates, the difference in the transmission distance between the second node and the reference node, and the direction in which the second node receives signals determine the coordinates of the second virtual node.

优选地,所述根据所述第二节点的坐标、所述第二节点与所述参考节 点到达待定位信号源的传输时差及所述第二节点接收信号的方向,确定第 二虚拟节点的坐标的步骤,具体包括:根据无线传播损耗计算公式,分别 计算第二节点在遮挡环境下的第一传输损耗及在无遮挡环境下的第二传输 损耗;根据所述第二节点与所述参考节点到达待定位信号源的传输时差计 算所述第二节点到达待定位信号源的第一距离;根据所述第二节点的RSSI 值、所述第一传输损耗、所述第一距离及所述参考节点与待定位信号源的 距离确定所述第二节点的RSSI值与所述第二虚拟节点的RSSI值的RSSI 差值;根据所述待定位信号源的功率、所述第二传输损耗及所述RSSI差值 计算在无遮挡环境下所述第二节点与所述待定位信号源的第二距离;根据 所述参考节点与待定位信号源的距离、第二距离及所述第二节点的坐标及 所述第二节点接收信号的方向确定所述第二虚拟节点的坐标。Preferably, the coordinates of the second virtual node are determined according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction in which the second node receives signals The step specifically includes: according to the wireless propagation loss calculation formula, respectively calculating the first transmission loss of the second node in the occlusion environment and the second transmission loss in the unobstructed environment; according to the second node and the reference node Calculate the first distance from the second node to the signal source to be located by the transmission time difference reaching the signal source to be located; according to the RSSI value of the second node, the first transmission loss, the first distance and the reference The distance between the node and the signal source to be located determines the RSSI difference between the RSSI value of the second node and the RSSI value of the second virtual node; according to the power of the signal source to be located, the second transmission loss and all The RSSI difference value calculates the second distance between the second node and the to-be-located signal source in an unobstructed environment; according to the distance between the reference node and the to-be-located signal source, the second distance and the The coordinates and the direction in which the second node receives the signal determine the coordinates of the second virtual node.

优选地,所述根据第三节点的坐标、所述第三节点与所述参考节点到 达待定位信号源的传输时差及所述第三节点接收信号的方向,计算第三虚 拟节点的坐标的步骤,具体包括:根据所述第三节点与所述参考节点到达 待定位信号源的传输时差,计算所述第三节点和所述参考节点之间的传输 距离的差值;根据所述第三节点的坐标、所述第三节点和所述参考节点之 间的传输距离的差值及所述第三节点接收信号的方向,确定所述第三虚拟节点的坐标。Preferably, the step of calculating the coordinates of the third virtual node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives signals , which specifically includes: calculating the difference between the transmission distances between the third node and the reference node according to the transmission time difference between the third node and the reference node reaching the signal source to be located; according to the third node The coordinates of the third virtual node are determined by the coordinates of the third node, the difference in the transmission distance between the third node and the reference node, and the direction in which the third node receives the signal.

优选地,所述根据第三节点的坐标、所述第三节点与所述参考节点到 达待定位信号源的传输时差及所述第三节点接收信号的方向,计算第三虚 拟节点的坐标的步骤,具体包括:根据无线传播损耗计算公式,分别计算 第三节点在遮挡环境下的第三传输损耗及在无遮挡环境下的第四传输损耗; 根据所述第三节点与所述参考节点到达待定位信号源的传输时差计算所述 第三节点到达待定位信号源的第三距离;根据所述第三节点的RSSI值、所 述第三传输损耗、所述第三距离及所述参考节点与待定位信号源的距离确 定所述第三节点的RSSI值与所述第三虚拟节点的RSSI值的RSSI差值;根 据所述待定位信号源的功率、所述第四传输损耗及所述RSSI差值计算在无 遮挡环境下所述第三节点与所述待定位信号源的第四距离;根据所述参考 节点与待定位信号源的距离、第四距离及所述第三节点的坐标及所述第三 节点接收信号的方向确定所述第三虚拟节点的坐标。Preferably, the step of calculating the coordinates of the third virtual node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives signals , which specifically includes: according to the wireless propagation loss calculation formula, respectively calculating the third transmission loss of the third node in the occlusion environment and the fourth transmission loss in the unobstructed environment; according to the third node and the reference node arrival to be determined The transmission time difference of the bit signal source calculates the third distance from the third node to the signal source to be located; according to the RSSI value of the third node, the third transmission loss, the third distance and the reference node and the The distance of the signal source to be located determines the RSSI difference between the RSSI value of the third node and the RSSI value of the third virtual node; according to the power of the signal source to be located, the fourth transmission loss and the RSSI The difference calculates the fourth distance between the third node and the to-be-located signal source in an unobstructed environment; according to the distance between the reference node and the to-be-located signal source, the fourth distance and the coordinates of the third node and The direction in which the third node receives the signal determines the coordinates of the third virtual node.

优选地,所述根据各定位节点接收待定位信号源信号的RSSI值分别计 算所述定位组合对应初始定位坐标的加权系数的步骤,具体包括:Preferably, the step of calculating the weighting coefficient of the corresponding initial positioning coordinates of the positioning combination according to the RSSI value that each positioning node receives the signal source signal to be positioned, specifically includes:

通过以下公式计算所述定位组合中各节点RSSI值的平均值:The average value of the RSSI values of each node in the positioning combination is calculated by the following formula:

Figure BDA0001752133240000041
Figure BDA0001752133240000041

其中,RSSIi表示第i个定位组合中各节点RSSI值的平均值,RSSIij表示第i个定位组合中第j个节点的RSSI值,n表示每组组合的节点的个 数;Wherein, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, RSSI ij represents the RSSI value of the jth node in the ith positioning combination, and n represents the number of nodes in each combination;

通过以下公式计算所述定位组合对应的加权系数,:The weighting coefficient corresponding to the positioning combination is calculated by the following formula:

Figure BDA0001752133240000051
Figure BDA0001752133240000051

其中,wi表示第i个定位组合的加权系数,RSSIi表示第i个定位组合 中各节点RSSI值的平均值,N表示定位组合的数量。Wherein, w i represents the weighting coefficient of the ith positioning combination, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, and N represents the number of positioning combinations.

优选地,所述根据所述初始定位坐标及与其对应的加权系数,计算所 述待定位信号源的定位坐标的步骤,通过以下公式计算所述待定位信号源 的定位坐标:Preferably, the step of calculating the positioning coordinates of the signal source to be located according to the initial positioning coordinates and the corresponding weighting coefficients, calculates the positioning coordinates of the signal source to be located by the following formula:

Figure BDA0001752133240000052
Figure BDA0001752133240000052

其中,wi表示第i个定位组合的加权系数,L表示待定位信号源信号的 定位坐标,Li表示第i个定位组合的初始定位坐标。Wherein, wi represents the weighting coefficient of the ith positioning combination, L represents the positioning coordinates of the signal source signal to be located, and Li represents the initial positioning coordinates of the ith positioning combination.

本发明实施例还提供一种基于虚拟节点的多站协同干扰定位系统,包括: 定位节点获取模块,用于获取监测节点中三个节点作为定位节点,分别为 参考节点、第二节点及第三节点;传输时差获取模块,用于计算所述第二 节点与所述参考节点到达待定位信号源的传输时差,第三节点与所述参考 节点到达待定位信号源的传输时差;第二虚拟节点坐标确定模块,用于根 据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源 的传输时差及所述第二节点接收信号的方向,确定第二虚拟节点的坐标; 第三虚拟节点坐标确定模块,用于根据所述第三节点的坐标、所述第三节 点与所述参考节点到达待定位信号源的传输时差及所述第三节点接收信号 的方向,确定第三虚拟节点的坐标;定位组合构建模块,用于根据所述参 考节点、所述第二节点与第二虚拟节点中之一、第三节点与第三虚拟节点 中之一构建定位组合;初始定位坐标计算模块,用于计算各所述定位组合对应的待定位信号源的初始定位坐标;加权系数获取模块,用于根据各定 位节点接收待定位信号源信号的RSSI值分别计算各所述定位组合对应初始 定位坐标的加权系数;定位坐标计算模块,用于根据所述初始定位坐标及 与其对应的加权系数,计算所述待定位信号源的定位坐标。An embodiment of the present invention further provides a virtual node-based multi-station cooperative interference positioning system, including: a positioning node acquisition module, configured to obtain three nodes in the monitoring nodes as positioning nodes, which are a reference node, a second node and a third node respectively. node; a transmission time difference acquisition module for calculating the transmission time difference between the second node and the reference node to the signal source to be located, and the transmission time difference between the third node and the reference node to the signal source to be located; the second virtual node The coordinate determination module is configured to determine the coordinate of the second virtual node according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction of the signal received by the second node. Coordinates; a third virtual node coordinate determination module, configured to, according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives the signal, determining the coordinates of the third virtual node; a positioning combination building module for building a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node; The initial positioning coordinate calculation module is used to calculate the initial positioning coordinates of the to-be-located signal source corresponding to each of the positioning combinations; the weighting coefficient acquisition module is used to calculate the RSSI values of the signals of the to-be-located signal sources received by each positioning node respectively. The positioning combination corresponds to the weighting coefficient of the initial positioning coordinates; the positioning coordinate calculation module is configured to calculate the positioning coordinates of the to-be-located signal source according to the initial positioning coordinates and the corresponding weighting coefficients.

本发明实施例还提供一种计算机设备,包括:至少一个处理器,以及 与所述至少一个处理器通信连接的存储器,其中,所述存储器存储有可被 所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行, 以使所述至少一个处理器执行上述的基于虚拟节点的多站协同干扰定位方 法。An embodiment of the present invention further provides a computer device, comprising: at least one processor, and a memory communicatively connected to the at least one processor, wherein the memory stores instructions executable by the at least one processor, The instructions are executed by the at least one processor, so that the at least one processor executes the above virtual node-based multi-station cooperative interference localization method.

本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储 介质存储有计算机指令,所述计算机指令用于使所述计算机执行上述的基 于虚拟节点的多站协同干扰定位方法。An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium stores computer instructions, and the computer instructions are used to cause the computer to execute the above-mentioned virtual node-based multi-station cooperative interference positioning method.

本发明技术方案,具有如下优点:The technical scheme of the present invention has the following advantages:

本发明实施例提供的基于虚拟节点的多站协同干扰定位方法及系统,将 监测节点接收待定位信号源信号的RSSI值最大的三个作为定位节点,将 RSSI值最大的定位节点作为参考节点,根据第二/三节点的坐标、第二/三 节点与参考节点到达待定位信号源的传输时差及第二/三节点接收信号的方 向,确定第二/三虚拟节点的坐标;根据参考节点、第二节点与第二虚拟节 点中之一、第三节点与第三虚拟节点中之一构建定位组合,计算各组合的 初始定位坐标,根据初始定位坐标及与对应的加权系数,计算待定位信号 源的定位坐标。通过构建虚拟节点的方式,丰富了协同定位的节点组合,降 低因节点分布所造成的定位误差,利用基于RSSI值的不同协同组合的定位结 果的加权方式,提高了定位精度。In the virtual node-based multi-station cooperative interference positioning method and system provided by the embodiment of the present invention, the monitoring nodes receive the three largest RSSI values of the signal source signal to be located as the positioning nodes, and the positioning nodes with the largest RSSI value are used as the reference nodes. Determine the coordinates of the second/third virtual node according to the coordinates of the second/third node, the transmission time difference between the second/third node and the reference node reaching the signal source to be located, and the direction of the signal received by the second/third node; according to the reference node, One of the second node and the second virtual node, and one of the third node and the third virtual node constructs a positioning combination, calculates the initial positioning coordinates of each combination, and calculates the signal to be positioned according to the initial positioning coordinates and the corresponding weighting coefficients. The location coordinates of the source. By constructing virtual nodes, the combination of nodes for co-location is enriched, the positioning error caused by node distribution is reduced, and the positioning accuracy is improved by using the weighting method of positioning results based on different cooperative combinations of RSSI values.

附图说明Description of drawings

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下 面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍, 显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普 通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获 得其他的附图。In order to illustrate the specific embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the specific embodiments or the prior art. Obviously, the accompanying drawings in the following description The drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

图1为本发明实施例中提供的基于虚拟节点的多站协同干扰定位方法 一个具体实施例的流程图;1 is a flowchart of a specific embodiment of a virtual node-based multi-station cooperative interference positioning method provided in an embodiment of the present invention;

图2为本发明实施例中确定定位节点及参考节点的一个具体示例的 流程图;Fig. 2 is a flow chart of a specific example of determining a positioning node and a reference node in an embodiment of the present invention;

图3为本发明实施例中计算第二虚拟节点坐标的一个具体示例的流 程图;3 is a flow chart of a specific example of calculating the coordinates of the second virtual node in the embodiment of the present invention;

图4为本发明实施例中构建虚拟节点的示意图;4 is a schematic diagram of constructing a virtual node in an embodiment of the present invention;

图5为本发明实施例中计算第二虚拟节点坐标的另一个具体示例的 流程图;5 is a flowchart of another specific example of calculating the coordinates of the second virtual node in the embodiment of the present invention;

图6为本发明实施例中计算第三虚拟节点坐标的一个具体示例的流 程图;6 is a flow chart of a specific example of calculating the coordinates of the third virtual node in the embodiment of the present invention;

图7为本发明实施例中计算第三虚拟节点坐标的另一个具体示例的 流程图;7 is a flowchart of another specific example of calculating the coordinates of the third virtual node in the embodiment of the present invention;

图8为本发明实施例中提供的基于虚拟节点的多站协同干扰定位系统 一个具体示例的流程图;8 is a flowchart of a specific example of a virtual node-based multi-station cooperative interference location system provided in an embodiment of the present invention;

图9为本发明实施例中提供的一种计算机设备的一个具体示例的流 程图。Fig. 9 is a flowchart of a specific example of a computer device provided in an embodiment of the present invention.

具体实施方式Detailed ways

下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然, 所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发 明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得 的所有其他实施例,都属于本发明保护的范围。The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present invention.

此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼 此之间未构成冲突就可以相互结合。In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

实施例1Example 1

本发明实施例提供一种基于虚拟节点的多站协同干扰定位方法,可用于 对干扰源进行定位,如图1所示,该方法包括如下步骤:An embodiment of the present invention provides a method for locating multi-station cooperative interference based on virtual nodes, which can be used to locate an interference source. As shown in Figure 1, the method includes the following steps:

步骤S1:获取监测节点中三个节点作为定位节点,分别为参考节点、 第二节点及第三节点。Step S1: Acquire three nodes in the monitoring nodes as positioning nodes, which are a reference node, a second node and a third node respectively.

在本发明实施例中,待定位信号源为一干扰源,在监测后台查询各个 监测节点接收到干扰源接收信号的强度,即各个监测节点RSSI值,来选取 定位节点和参考节点,因此上述的步骤S1,如图2所示,具体包括如下步 骤:In the embodiment of the present invention, the signal source to be located is an interference source, and in the monitoring background, the strength of the received signal received by each monitoring node from the interference source, that is, the RSSI value of each monitoring node, is used to select the positioning node and the reference node. Therefore, the above-mentioned Step S1, as shown in Figure 2, specifically includes the following steps:

步骤S11:获取监测节点中接收待定位信号源信号的RSSI值按高低排 序的前三个节点按高低排序的前三个节点作为定位节点。Step S11: obtain the first three nodes that receive the RSSI value of the signal source signal to be located in the monitoring node, and the first three nodes that are sorted by high and low are used as positioning nodes.

步骤S12:分别获取待定位信号源信号到达定位节点的RSSI值。Step S12: respectively acquiring the RSSI value of the signal source signal to be located arriving at the positioning node.

步骤S13:将RSSI值最大的定位节点作为参考节点,将另两个定位节 点分别作为第二节点、第三节点。Step S13: The positioning node with the largest RSSI value is used as the reference node, and the other two positioning nodes are respectively used as the second node and the third node.

在确定定位节点及参考节点后,根据监测后台获取各定位节点的信息 对干扰源进行定位,各定位节点的信息包括:各定位节点的坐标、干扰源 信号到达各定位节点的传输时间。在本发明实施例中,监测后台采用广义 互相关法来计算干扰源信号到达不同定位节点之间的传输时间。After the positioning node and reference node are determined, the interference source is located according to the information of each positioning node obtained from the monitoring background. The information of each positioning node includes: the coordinates of each positioning node and the transmission time of the interference source signal reaching each positioning node. In the embodiment of the present invention, the monitoring background adopts the generalized cross-correlation method to calculate the transmission time between the interference source signal reaching different positioning nodes.

步骤S2:计算第二节点与参考节点到达待定位信号源的传输时差,第 三节点与参考节点到达待定位信号源的传输时差。Step S2: Calculate the transmission time difference between the second node and the reference node to the signal source to be located, and the transmission time difference between the third node and the reference node to arrive at the signal source to be located.

在本发明实施例中,根据干扰源信号到达不同定位节点之间的传输时 间,来算第二节点与参考节点到达待定位信号源的传输时差,第三节点与 参考节点到达待定位信号源的传输时差。In the embodiment of the present invention, the transmission time difference between the second node and the reference node reaching the signal source to be located is calculated according to the transmission time between the interference source signal reaching different positioning nodes, and the time difference between the third node and the reference node reaching the signal source to be located is calculated. transmission time difference.

步骤S3:根据第二节点的坐标、第二节点与参考节点到达待定位信号 源的传输时差及第二节点接收信号的方向,确定第二虚拟节点的坐标。Step S3: Determine the coordinates of the second virtual node according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located and the direction in which the second node receives the signal.

本发明实施例中,参考节点用N1表示,第二节点用N2表示,第二虚 拟节点用N2’表示,考虑到接收节点有无受到遮挡等影响传输时延,可以分 4种场景来确定第二虚拟节点N2’的坐标:In the embodiment of the present invention, the reference node is represented by N1, the second node is represented by N2, and the second virtual node is represented by N2'. Considering whether the receiving node is blocked or not, the transmission delay is affected, and four scenarios can be used to determine the first The coordinates of the two virtual nodes N2':

(1)当参考节点N1无遮挡、第二节点N2无遮挡时,如图3所示,通过 以下步骤计算第二虚拟节点N2’的坐标:(1) When the reference node N1 is not blocked and the second node N2 is not blocked, as shown in Figure 3, the coordinates of the second virtual node N2' are calculated through the following steps:

步骤S311:根据第二节点与参考节点到达待定位信号源的传输时差, 计算第二节点和参考节点之间的传输距离的差值。Step S311: Calculate the difference between the transmission distances between the second node and the reference node according to the transmission time difference between the second node and the reference node when they arrive at the signal source to be located.

本发明实施例中,根据公式(1)来计算第二节点和参考节点之间的传 输距离的差值:In the embodiment of the present invention, the difference value of the transmission distance between the second node and the reference node is calculated according to formula (1):

Figure BDA0001752133240000101
Figure BDA0001752133240000101

其中,(xi,yi)表示第i个节点的坐标,(x,y)表示干扰源T的坐标,ri表示干扰源与第i个监测节点之间的距离,rj表示干扰源与第j个监测节点 之间的距离,dij表示第i个和第j个节点之间的距离的差值,TDOAij为第 i个和第j个节点到达待定位信号源的传输时差(例如,可采用现有技术中 的广义互相关法来计算信号到达不同监测节点之间的时间差),n表示监测节 点的数量,C表示电磁波传播速度。Among them, (x i , y i ) represents the coordinates of the ith node, (x, y) represents the coordinates of the interference source T, ri represents the distance between the interference source and the ith monitoring node, and r j represents the interference source The distance from the jth monitoring node, d ij represents the difference between the distances between the ith node and the jth node, and TDOA ij is the transmission time difference between the ith node and the jth node reaching the signal source to be located ( For example, the generalized cross-correlation method in the prior art can be used to calculate the time difference between signals arriving at different monitoring nodes), n represents the number of monitoring nodes, and C represents the propagation speed of electromagnetic waves.

因此,根据上述公式(1)可计算第二节点N2和参考节点N1的传输距 离的差值d12Therefore, the difference d 12 of the transmission distance between the second node N2 and the reference node N1 can be calculated according to the above formula (1).

步骤S312:根据第二节点的坐标、第二节点和参考节点之间的传输距 离的差值及第二节点接收信号的方向,确定第二虚拟节点的坐标。Step S312: Determine the coordinates of the second virtual node according to the coordinates of the second node, the difference in the transmission distance between the second node and the reference node, and the direction in which the second node receives the signal.

本发明实施例中,如图4所示,构建第二虚拟节点N2’,使得第二节点 N2与第二虚拟节点N2’的距离为d12,第二节点N2移动的方向依据其接收 到的干扰源T的信号方向确定(如果接收到左侧的信号,则确定移动方向 为向左侧移动,右侧同理),以第二节点N2的绝对坐标为参考点,根据第 二节点和参考节点之间的传输距离的差值d12,在横轴和纵轴方向分别移动 相同的距离,来确定第二虚拟节点N2’的坐标。In the embodiment of the present invention, as shown in FIG. 4 , the second virtual node N2' is constructed so that the distance between the second node N2 and the second virtual node N2' is d 12 , and the moving direction of the second node N2 is based on the received The signal direction of the interference source T is determined (if a signal from the left is received, the moving direction is determined to be moving to the left, and the same is true for the right), and the absolute coordinates of the second node N2 are used as the reference point. The difference d 12 of the transmission distance between the nodes is moved by the same distance in the direction of the horizontal axis and the vertical axis respectively to determine the coordinates of the second virtual node N2 ′.

(2)当参考节点N1无遮挡、第二节点N2有遮挡时,如图5所示,通过 以下步骤计算N2’的坐标:(2) When the reference node N1 is not blocked and the second node N2 is blocked, as shown in Figure 5, the coordinates of N2' are calculated through the following steps:

步骤S321:根据无线传播损耗计算公式,分别计算第二节点在遮挡环 境下的第一传输损耗及在无遮挡环境下的第二传输损耗。Step S321: According to the wireless propagation loss calculation formula, calculate the first transmission loss of the second node under the shielding environment and the second transmission loss under the unshielded environment, respectively.

本发明实施例中,在理想的无遮挡条件下,计算该第二传输损耗的无 线传播损耗公式如下:In the embodiment of the present invention, under ideal conditions without shielding, the wireless propagation loss formula for calculating the second transmission loss is as follows:

P=32.44+20lgd+20lgf (2),P loss =32.44+20lgd+20lgf (2),

有遮挡条件下,计算该第一传输损耗的无线传播损耗公式如下:Under the condition of occlusion, the wireless propagation loss formula for calculating the first transmission loss is as follows:

P (1)=32.44+20lgd+20lgf+Pw (3),P loss (1) = 32.44+20lgd+20lgf+P w (3),

其中,P表示无遮挡条件下的无线传输损耗、P (1)表示有遮挡条件下的 无线传输损耗,d为传输距离,f为干扰频率,Pw表示阻挡物、多径等有遮 挡条件下造成的损耗。Among them, P loss represents the wireless transmission loss under the condition of no occlusion, P loss (1) represents the wireless transmission loss under the occlusion condition, d is the transmission distance, f is the interference frequency, and P w represents the obstruction, multipath, etc. with occlusion loss due to conditions.

步骤S322:根据第二节点与参考节点到达待定位信号源的传输时差计 算第二节点到达待定位信号源的第一距离。Step S322: Calculate the first distance that the second node reaches the signal source to be located according to the transmission time difference between the second node and the reference node reaching the signal source to be located.

本发明实施例中,首先根据TDOA定位算法中的SC(shrinking-circle)定 位算法算出参考节点N1的最佳半径,然后根据参考节点N1的最佳半径及第 二节点N2与参考节点N1到达待定位信号源的传输时差计算得出d2In the embodiment of the present invention, the optimal radius of the reference node N1 is first calculated according to the SC (shrinking-circle) localization algorithm in the TDOA localization algorithm, and then according to the optimal radius of the reference node N1 and the arrival of the second node N2 and the reference node N1 to be determined The transit time difference of the bit signal source is calculated to obtain d 2 .

步骤S323:根据第二节点的RSSI值、第一传输损耗、第一距离及参考 节点与待定位信号源的距离确定第二节点的RSSI值与第二虚拟节点的 RSSI值的RSSI差值。Step S323: Determine the RSSI difference between the RSSI value of the second node and the RSSI value of the second virtual node according to the RSSI value of the second node, the first transmission loss, the first distance, and the distance between the reference node and the signal source to be located.

本发明实施例中,干扰源发射功率为Pf,第二节点N2接收的RSSI值 为P2,有遮挡的情况下,

Figure BDA0001752133240000121
根据公式(3)可得出:In the embodiment of the present invention, the transmit power of the interference source is P f , and the RSSI value received by the second node N2 is P 2 .
Figure BDA0001752133240000121
According to formula (3), it can be obtained:

Figure BDA0001752133240000122
Figure BDA0001752133240000122

其中,C为32.44+20lgf,因为监测的干扰频率f是已知的,所以C是 确定值。本发明实施例中,通过现有技术中的功率和信号强度的转换公式, 将干扰源发射功率Pf转换成信号强度的单位后进行差运算。Among them, C is 32.44+20lgf, because the monitored interference frequency f is known, so C is a certain value. In the embodiment of the present invention, a difference calculation is performed after converting the transmission power P f of the interference source into a unit of signal strength by using the conversion formula of power and signal strength in the prior art.

为使得构建的第二虚拟节点N2’与参考节点N1等节点组成的三角形的 几何中心为干扰源的位置,则第二虚拟节点N2’到达干扰源的距离应该与参 考节点N1到达干扰源的距离d1一致,设第二虚拟节点N2’接收RSSI的值为 P2 (1),则In order to make the geometric center of the triangle formed by the second virtual node N2' and the reference node N1 to be the position of the interference source, the distance between the second virtual node N2' and the interference source should be the same as the distance from the reference node N1 to the interference source. d 1 is the same, suppose the value of RSSI received by the second virtual node N2' is P 2 (1) , then

P2 (1)=Pf-20logd2'-C-Pw (5),P 2 (1) =P f -20logd 2 '-CP w (5),

其中d2’为第二虚拟节点N2’与干扰源之间的距离,和参考节点N1与干扰 源之间的距离d1相同。Wherein d 2 ′ is the distance between the second virtual node N2 ′ and the interference source, which is the same as the distance d 1 between the reference node N1 and the interference source.

因此,第二节点N2与第二虚拟节点N2’的RSSI差值P(Δ)为:Therefore, the RSSI difference P(Δ) between the second node N2 and the second virtual node N2' is:

P(Δ)=P2-P2 (1)=20lgd2'-20lgd2 (6),P(Δ)=P 2 -P 2 (1) =20lgd 2 '-20lgd 2 (6),

步骤S324:根据待定位信号源的功率、第二传输损耗及RSSI差值计算 在无遮挡环境下第二节点与待定位信号源的第二距离。Step S324: Calculate the second distance between the second node and the signal source to be located in an unobstructed environment according to the power of the signal source to be located, the second transmission loss and the RSSI difference.

在本发明实施例中,无遮挡的条件下的第二虚拟节点N2’的接收RSSI 值为P2 (2),由于第二虚拟节点N2’与干扰源的距离和参考节点N1与干扰源的 距离相等,且都在理想无遮挡条件下,所以In the embodiment of the present invention, the received RSSI value of the second virtual node N2' under the condition of no blocking is P 2 (2) , because the distance between the second virtual node N2' and the interference source and the distance between the reference node N1 and the interference source are The distances are equal, and they are all under ideal unobstructed conditions, so

P2 (2)=Pf-20lgd2'-C=P1 (7),P 2 (2) =P f -20lgd 2 '-C=P 1 (7),

其中,P1是参考节点N1接收到的RSSI值,可以通过监测后台直接监测 得到。Among them, P 1 is the RSSI value received by the reference node N1, which can be directly obtained by monitoring the background.

理想无遮挡条件下,第二节点N2接收的RSSI的值为P2′,则Under ideal unobstructed conditions, the RSSI value received by the second node N2 is P 2 ′, then

P2'=P2 (2)+P(Δ) (8),P 2 '=P 2 (2) +P(Δ) (8),

根据上述公式(4),可以计算出在理想场景下第二节点N2处距离干扰 源的距离d2 *,那么According to the above formula (4), the distance d 2 * from the second node N2 to the interference source can be calculated in the ideal scenario, then

Figure BDA0001752133240000131
Figure BDA0001752133240000131

上述的公式(9)结合公式(7),可以计算得出第二节点N2距离干扰源 的距离d2 *The above formula (9) and formula (7) can be used to calculate the distance d 2 * from the second node N2 to the interference source.

步骤S325:根据参考节点与待定位信号源的距离、第二距离及第二节 点的坐标及第二节点接收信号的方向确定第二虚拟节点的坐标。Step S325: Determine the coordinates of the second virtual node according to the distance between the reference node and the signal source to be located, the second distance and the coordinates of the second node and the direction in which the second node receives the signal.

本发明实施例中,第二节点N2到第二虚拟节点N2’需要移动的直线距 离为d2 *-d1,以第二节点N2的绝对坐标为参考点,依据第二节点N2接收 到的干扰信号方向确定(如果接收到左侧的信号,则确定移动方向为向左 侧移动,右侧同理),在横轴和纵轴方向分别移动相同的距离,来确定第二 虚拟节点N2’的坐标。In the embodiment of the present invention, the straight-line distance that needs to be moved from the second node N2 to the second virtual node N2' is d 2 * -d 1 , and the absolute coordinates of the second node N2 are used as the reference point. The direction of the interference signal is determined (if the signal from the left is received, the moving direction is determined to move to the left, and the same is true for the right), and the second virtual node N2' is determined by moving the same distance in the horizontal and vertical directions respectively. coordinate of.

(3)当参考节点N1有遮挡、第二节点N2无遮挡时,由于参考节点N1 有遮挡的情况下,接收的时延变长,导致参考节点N1预估到达干扰源的距 离要远,其他无遮挡的节点与其比较,计算的时延差相对变小,移动的距 离也就变短,但是考虑到N1是作为参考节点,其计算的最佳半径将为其他 虚拟节点的构建提供参考,所以构建的虚拟节点的半径需与参考节点N1一 致,而第二节点N2又是无遮挡环境下进行传输,虽然参考节点N1和构建的 第二虚拟节点N2’距离干扰源都变远了,但两者到达距离一致,对构建的三 角形的几何中心没有影响,因此这里确定第二虚拟节点N2’坐标的方法与场 景(1)的方法相同,在此不再赘述。(3) When the reference node N1 is blocked and the second node N2 is not blocked, the receiving delay becomes longer because the reference node N1 is blocked, so that the estimated distance from the reference node N1 to the interference source is longer, and other Compared with the unobstructed node, the calculated delay difference is relatively small, and the moving distance is also shortened. However, considering that N1 is used as a reference node, the calculated optimal radius will provide a reference for the construction of other virtual nodes, so The radius of the constructed virtual node must be consistent with the reference node N1, and the second node N2 transmits in an unobstructed environment. Although the reference node N1 and the constructed second virtual node N2' are farther away from the interference source, the two Therefore, the method for determining the coordinates of the second virtual node N2' is the same as the method for scene (1), and will not be repeated here.

(4)当参考节点N1有遮挡、第二节点N2有遮挡时,各个定位节点受 到的影响是一样的,相对时延差也是准确的,因此按照场景(1)的方法计算 第二虚拟节点N2’的坐标,在此不再赘述。(4) When the reference node N1 is blocked and the second node N2 is blocked, each positioning node is affected by the same, and the relative delay difference is also accurate. Therefore, the second virtual node N2 is calculated according to the method of scenario (1). ' coordinates, which will not be repeated here.

S4:根据第三节点的坐标、第三节点与参考节点到达待定位信号源的传 输时差及第三节点接收信号的方向,确定第三虚拟节点的坐标。S4: Determine the coordinates of the third virtual node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives the signal.

本发明实施例中,参考节点用N1表示,第三节点用N3表示,第三虚 拟节点用N3’表示,考虑到接收节点有无受到遮挡等影响传输时延,与确定 第二虚拟节点N2’的坐标的原理一样,分4种场景来确定第三虚拟节点用 N3’的坐标:In the embodiment of the present invention, the reference node is represented by N1, the third node is represented by N3, and the third virtual node is represented by N3'. The principle of the coordinates is the same, and there are 4 scenarios to determine the coordinates of the third virtual node with N3':

场景(一)当参考节点N1无遮挡、第三节点N3无遮挡时,(三)当 参考节点N1有遮挡、第三节点N3无遮挡时,(四)当参考节点N1有遮挡、 第三节点N3有遮挡时,如图6所示,包括如下步骤:Scenario (1) When the reference node N1 is not blocked and the third node N3 is not blocked, (3) When the reference node N1 is blocked and the third node N3 is not blocked, (4) When the reference node N1 is blocked and the third node is blocked When N3 is blocked, as shown in Figure 6, it includes the following steps:

步骤S411:根据第三节点与参考节点到达待定位信号源的传输时差, 计算第三节点和参考节点之间的传输距离的差值;Step S411: Calculate the difference in the transmission distance between the third node and the reference node according to the transmission time difference between the third node and the reference node reaching the signal source to be located;

步骤S412:根据第三节点的坐标、第三节点和参考节点之间的传输距 离的差值及第三节点接收信号的方向,确定第三虚拟节点的坐标。Step S412: Determine the coordinates of the third virtual node according to the coordinates of the third node, the difference in the transmission distance between the third node and the reference node, and the direction in which the third node receives the signal.

与上述的场景(1)、(3)、(4)确定第二虚拟节点N2’的坐标的原理一 样,不再赘述。The principle of determining the coordinates of the second virtual node N2' is the same as the above-mentioned scenarios (1), (3) and (4), and will not be repeated here.

场景(二)当参考节点N1无遮挡、第三节点N3有遮挡时,如图7所 示,通过以下步骤计算第三虚拟节点N3’的坐标:Scenario (two) when the reference node N1 is not blocked and the third node N3 is blocked, as shown in Figure 7, the coordinates of the third virtual node N3 ' are calculated by the following steps:

步骤S421:根据无线传播损耗计算公式,分别计算第三节点在遮挡环 境下的第三传输损耗及在无遮挡环境下的第四传输损耗。Step S421: According to the wireless propagation loss calculation formula, calculate the third transmission loss of the third node under the shielding environment and the fourth transmission loss under the unshielded environment, respectively.

步骤S422:根据第三节点与参考节点到达待定位信号源的传输时差计 算第三节点到达待定位信号源的第三距离。Step S422: Calculate the third distance that the third node reaches the signal source to be located according to the transmission time difference between the third node and the reference node reaching the signal source to be located.

步骤S323:根据第三节点的RSSI值、第一传输损耗、第三距离及参考 节点与待定位信号源的距离确定第三节点的RSSI值与第三虚拟节点的 RSSI值的RSSI差值。Step S323: Determine the RSSI difference between the RSSI value of the third node and the RSSI value of the third virtual node according to the RSSI value of the third node, the first transmission loss, the third distance and the distance between the reference node and the signal source to be located.

步骤S324:根据待定位信号源的功率、第四传输损耗及RSSI差值计算 在无遮挡环境下第三节点与待定位信号源的第四距离。Step S324: Calculate the fourth distance between the third node and the signal source to be located under the unobstructed environment according to the power of the signal source to be located, the fourth transmission loss and the RSSI difference.

步骤S325:根据参考节点与待定位信号源的距离、第四距离及第三节 点的坐标及第三节点接收信号的方向确定第三虚拟节点的坐标。Step S325: Determine the coordinates of the third virtual node according to the distance between the reference node and the signal source to be located, the fourth distance and the coordinates of the third node and the direction in which the third node receives the signal.

上述场景(二)当参考节点N1无遮挡、第三节点N3有遮挡时,计算 计算第三虚拟节点N3’的坐标与上述的场景(2)下确定第二虚拟节点N2’ 的坐标的原理一样,不再赘述。In the above scenario (2), when the reference node N1 is not blocked and the third node N3 is blocked, the calculation and calculation of the coordinates of the third virtual node N3' is the same as the principle of determining the coordinates of the second virtual node N2' in the above-mentioned scenario (2). ,No longer.

步骤S5:根据参考节点、第二节点与第二虚拟节点中之一、第三节点 与第三虚拟节点中之一构建定位组合。Step S5: construct a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node.

本发明实施例中,以(N1,N2,N3)、(N1,N2’,N3)、(N1,N2, N3’)、(N1,N2’,N3’)构建四个协同定位组合。In the embodiment of the present invention, four co-localization combinations are constructed by (N1, N2, N3), (N1, N2', N3), (N1, N2, N3'), and (N1, N2', N3').

步骤S6:计算各定位组合对应的待定位信号源的初始定位坐标。Step S6: Calculate the initial positioning coordinates of the signal source to be positioned corresponding to each positioning combination.

本发明实施例中,分别按照现有技术中的SC(shrinking-circle)定位算 法进行计算,得到上述四个协同定位组合对应的四个初始定位坐标L1,L2,L3,L4。In the embodiment of the present invention, the calculation is performed according to the SC (shrinking-circle) positioning algorithm in the prior art, respectively, to obtain four initial positioning coordinates L1, L2, L3, L4 corresponding to the above-mentioned four co-location combinations.

步骤S7:根据各定位节点接收待定位信号源信号的RSSI值分别计算各 定位组合对应初始定位坐标的加权系数。Step S7: according to the RSSI value that each positioning node receives the signal source signal to be positioned, calculate the weighting coefficient of each positioning combination corresponding to the initial positioning coordinates respectively.

本发明实施例中,第二虚拟节点N2’和第三虚拟节点N3’由于和参考 节点N1距离干扰源的距离相等,因此第二虚拟节点N2’和第三虚拟节点 N3’的RSSI值与参考节点的RSSI值相等。In this embodiment of the present invention, since the second virtual node N2' and the third virtual node N3' are at the same distance from the interference source as the reference node N1, the RSSI values of the second virtual node N2' and the third virtual node N3' are the same as the reference node N1. The RSSI values of the nodes are equal.

本发明实施例中,通过以下公式计算定位组合中各节点RSSI值的平 均值:In the embodiment of the present invention, the average value of the RSSI value of each node in the positioning combination is calculated by the following formula:

Figure BDA0001752133240000161
Figure BDA0001752133240000161

其中,RSSIi表示第i个定位组合中各节点RSSI值的平均值,RSSIij表 示第i个定位组合中第j个节点的RSSI值,n表示每个定位组合的节点的 个数。本发明实施例中每个定位组合中都包含三个节点,因此n为3。Among them, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, RSSI ij represents the RSSI value of the jth node in the ith positioning combination, and n represents the number of nodes in each positioning combination. In this embodiment of the present invention, each positioning combination includes three nodes, so n is 3.

通过以下公式计算定位组合对应的加权系数:The weighting factor corresponding to the positioning combination is calculated by the following formula:

Figure BDA0001752133240000162
Figure BDA0001752133240000162

其中,wi表示第i个定位组合的加权系数,RSSIi表示第i个定位组合 中各节点RSSI值的平均值,N表示定位组合的数量。本发明实施例中,有 四个定位组合,因此N为4。Wherein, w i represents the weighting coefficient of the ith positioning combination, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, and N represents the number of positioning combinations. In this embodiment of the present invention, there are four positioning combinations, so N is 4.

步骤S8:根据初始定位坐标及与其对应的加权系数,计算待定位信号 源的定位坐标。Step S8: according to the initial positioning coordinates and the corresponding weighting coefficients, calculate the positioning coordinates of the signal source to be positioned.

本发明实施例中,通过以下公式计算待定位信号源的定位坐标:In the embodiment of the present invention, the positioning coordinates of the signal source to be located are calculated by the following formula:

Figure BDA0001752133240000171
Figure BDA0001752133240000171

其中,wi表示第i个定位组合的加权系数,L表示待定位信号源信号的 定位坐标,Li表示第i个定位组合的初始定位坐标。Wherein, wi represents the weighting coefficient of the ith positioning combination, L represents the positioning coordinates of the signal source signal to be located, and Li represents the initial positioning coordinates of the ith positioning combination.

本发明实施例提供的基于虚拟节点的多站协同干扰定位方法,将监测节 点接收待定位信号源信号的RSSI值按高低排序的前三个节点作为定位节点, 将RSSI值最大的定位节点作为参考节点,根据第二/三节点的坐标、第二/ 三节点分别与参考节点到达待定位信号源的传输时差及第二/三节点接收信 号的方向,确定第二/三虚拟节点的坐标;根据参考节点、第二节点与第二 虚拟节点中之一、第三节点与第三虚拟节点中之一构建定位组合,计算各 组合的初始定位坐标,根据初始定位坐标及与对应的加权系数,计算待定 位信号源的定位坐标。通过构建虚拟节点的方式,丰富了协同定位的节点组 合,降低因节点分布所造成的定位误差,利用基于RSSI值的不同协同组合的 定位结果的加权方式,提高了定位精度。In the virtual node-based multi-station cooperative interference positioning method provided by the embodiment of the present invention, the first three nodes in which the monitoring node receives the RSSI value of the signal source signal to be located in the order of high and low are used as the positioning node, and the positioning node with the largest RSSI value is used as a reference. node, determine the coordinates of the second/third virtual node according to the coordinates of the second/third node, the transmission time difference between the second/third node and the reference node reaching the signal source to be located, and the direction in which the second/third node receives the signal; The reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node constructs a positioning combination, calculates the initial positioning coordinates of each combination, and calculates the initial positioning coordinates and the corresponding weighting coefficient according to the initial positioning coordinates. Positioning coordinates of the signal source to be positioned. By constructing virtual nodes, the combination of nodes for co-location is enriched, the positioning error caused by node distribution is reduced, and the positioning accuracy is improved by using the weighting method of positioning results based on different cooperative combinations of RSSI values.

实施例2Example 2

本发明实施例提供一种基于虚拟节点的多站协同干扰定位系统,如图8 所示,该基于虚拟节点的多站协同干扰定位系统包括:An embodiment of the present invention provides a virtual node-based multi-station cooperative interference positioning system. As shown in FIG. 8 , the virtual node-based multi-station cooperative interference positioning system includes:

定位节点获取模块1,用于获取监测节点中三个节点作为定位节点,分 别为参考节点、第二节点及第三节点。此模块具体执行实施例1中步骤S11~ S13的方法,在此不再赘述。The positioning node obtaining module 1 is used to obtain three nodes in the monitoring nodes as positioning nodes, which are a reference node, a second node and a third node respectively. This module specifically executes the methods of steps S11 to S13 in Embodiment 1, and details are not described herein again.

传输时差获取模块2,用于计算第二节点与参考节点到达待定位信号源 的传输时差,第三节点与参考节点到达待定位信号源的传输时差。此模块 具体执行实施例1中步骤S2的方法,在此不再赘述。The transmission time difference acquisition module 2 is used to calculate the transmission time difference between the second node and the reference node to the signal source to be located, and the transmission time difference between the third node and the reference node to the signal source to be located. This module specifically executes the method of step S2 in Embodiment 1, which is not repeated here.

第二虚拟节点坐标确定模块3,用于根据第二节点的坐标、第二节点与 参考节点到达待定位信号源的传输时差及所述第二节点接收信号的方向, 确定第二虚拟节点的坐标。此模块参见实施例1中记录的步骤S311~S312、 步骤S321~S325,在此不再赘述。The second virtual node coordinate determination module 3 is configured to determine the coordinates of the second virtual node according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction in which the second node receives the signal . For this module, refer to steps S311-S312 and steps S321-S325 recorded in Embodiment 1, and details are not repeated here.

第三虚拟节点坐标确定模块4,用于根据第三节点的坐标、第三节点与 参考节点到达待定位信号源的传输时差及所述第三节点接收信号的方向, 确定第三虚拟节点的坐标。此模块参见实施例1中记录的步骤S411~S412、 步骤S421~S425,在此不再赘述。The third virtual node coordinate determination module 4 is configured to determine the coordinates of the third virtual node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives the signal . For this module, refer to steps S411 to S412 and steps S421 to S425 recorded in Embodiment 1, which will not be repeated here.

定位组合构建模块5,用于根据参考节点、第二节点与第二虚拟节点中 之一、第三节点与第三虚拟节点中之一构建定位组合。此模块参见实施例1 中记录的步骤S5,在此不再赘述。The positioning combination building module 5 is configured to construct a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node. For this module, refer to step S5 recorded in Embodiment 1, and details are not repeated here.

初始定位坐标计算模块6,用于计算各定位组合对应的待定位信号源的 初始定位坐标。此模块参见实施例1中记录的步骤S6,在此不再赘述。The initial positioning coordinate calculation module 6 is used to calculate the initial positioning coordinates of the signal source to be positioned corresponding to each positioning combination. For this module, refer to step S6 recorded in Embodiment 1, and details are not repeated here.

加权系数获取模块7,用于根据各定位节点接收待定位信号源信号的 RSSI值分别计算各定位组合对应初始定位坐标的加权系数;此模块参见实 施例1中记录的步骤S7,在此不再赘述。The weighting coefficient acquisition module 7 is used to calculate the weighting coefficient of each positioning combination corresponding to the initial positioning coordinates according to the RSSI value of the signal source signal to be positioned received by each positioning node; this module refers to the step S7 recorded in the embodiment 1, and will not be repeated here. Repeat.

定位坐标计算模块8,用于根据初始定位坐标及与其对应的加权系数, 计算待定位信号源的定位坐标。此模块参见实施例1中记录的步骤S8,在 此不再赘述。The positioning coordinate calculation module 8 is configured to calculate the positioning coordinates of the signal source to be located according to the initial positioning coordinates and the corresponding weighting coefficients. For this module, refer to step S8 recorded in Embodiment 1, and details are not repeated here.

本发明实施例提供的基于虚拟节点的多站协同干扰定位系统,将监测节 点接收待定位信号源信号的RSSI值按高低排序的前三个节点作为定位节点, 将RSSI值最大的定位节点作为参考节点,根据第二/三节点的坐标、第二/ 三节点分别与参考节点到达待定位信号源的传输时差及第二/三节点接收信 号的方向,确定第二/三虚拟节点的坐标;根据参考节点、第二节点与第二 虚拟节点中之一、第三节点与第三虚拟节点中之一构建定位组合,计算各 组合的初始定位坐标,根据初始定位坐标及与对应的加权系数,计算待定 位信号源的定位坐标。通过构建虚拟节点的方式,丰富了协同定位的节点组 合,降低因节点分布所造成的定位误差,利用基于RSSI值的不同协同组合的 定位结果的加权方式,提高了定位精度。In the virtual node-based multi-station cooperative interference positioning system provided by the embodiment of the present invention, the first three nodes in which the RSSI values of the signal source signals to be located received by the monitoring nodes are sorted in order of high and low are used as the positioning nodes, and the positioning node with the largest RSSI value is used as a reference. node, determine the coordinates of the second/third virtual node according to the coordinates of the second/third node, the transmission time difference between the second/third node and the reference node reaching the signal source to be located, and the direction in which the second/third node receives the signal; The reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node constructs a positioning combination, calculates the initial positioning coordinates of each combination, and calculates the initial positioning coordinates and the corresponding weighting coefficient according to the initial positioning coordinates. Positioning coordinates of the signal source to be positioned. By constructing virtual nodes, the combination of nodes for co-location is enriched, the positioning error caused by node distribution is reduced, and the positioning accuracy is improved by using the weighting method of positioning results based on different cooperative combinations of RSSI values.

实施例3Example 3

本发明实施例提供一种计算机设备,如图9所示,包括:至少一个处 理器401,例如CPU(Central Processing Unit,中央处理器),至少一个通 信接口403,存储器404,至少一个通信总线402。其中,通信总线402用 于实现这些组件之间的连接通信。其中,通信接口403可以包括显示屏 (Display)、键盘(Keyboard),可选通信接口403还可以包括标准的有线 接口、无线接口。存储器404可以是高速RAM存储器(Ramdom Access Memory,易挥发性随机存取存储器),也可以是非不稳定的存储器 (non-volatile memory),例如至少一个磁盘存储器。存储器404可选的还 可以是至少一个位于远离前述处理器401的存储装置。其中处理器401可 以执行图1描述的基于虚拟节点的多站协同干扰定位方法,存储器404中存 储一组程序代码,且处理器401调用存储器404中存储的程序代码,以用 于执行实施例1中的基于虚拟节点的多站协同干扰定位方法。An embodiment of the present invention provides a computer device, as shown in FIG. 9 , including: at least one processor 401 , such as a CPU (Central Processing Unit, central processing unit), at least one communication interface 403 , memory 404 , and at least one communication bus 402 . Among them, the communication bus 402 is used to realize the connection communication between these components. Wherein, the communication interface 403 may include a display screen (Display) and a keyboard (Keyboard), and the optional communication interface 403 may also include a standard wired interface and a wireless interface. The memory 404 may be a high-speed RAM memory (Ramdom Access Memory, volatile random access memory), or may be a non-volatile memory (non-volatile memory), such as at least one disk memory. The memory 404 can optionally also be at least one storage device located remotely from the aforementioned processor 401. The processor 401 can execute the virtual node-based multi-station cooperative interference location method described in FIG. 1 , a set of program codes are stored in the memory 404 , and the processor 401 calls the program codes stored in the memory 404 for executing Embodiment 1 A virtual node-based multi-station cooperative interference localization method in .

其中,通信总线402可以是外设部件互连标准(peripheral componentinterconnect,简称PCI)总线或扩展工业标准结构(extended industry standardarchitecture,简称EISA)总线等。通信总线402可以分为地址总线、数据 总线、控制总线等。为便于表示,图9中仅用一条粗线表示,但并不表示 仅有一根总线或一种类型的总线。The communication bus 402 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus or the like. The communication bus 402 can be divided into an address bus, a data bus, a control bus, and the like. For ease of presentation, only one thick line is shown in Figure 9, but it does not mean that there is only one bus or one type of bus.

其中,存储器404可以包括易失性存储器(英文:volatile memory), 例如随机存取存储器(英文:random-access memory,缩写:RAM);存储 器也可以包括非易失性存储器(英文:non-volatile memory),例如快闪存 储器(英文:flash memory),硬盘(英文:harddisk drive,缩写:HDD) 或固态硬盘(英文:solid-state drive,缩写:SSD);存储器404还可以包括 上述种类的存储器的组合。The memory 404 may include volatile memory (English: volatile memory), such as random-access memory (English: random-access memory, abbreviation: RAM); the memory may also include non-volatile memory (English: non-volatile memory) memory), such as flash memory (English: flash memory), hard disk (English: harddisk drive, abbreviation: HDD) or solid-state drive (English: solid-state drive, abbreviation: SSD); the memory 404 may also include the above-mentioned types of memory The combination.

其中,处理器401可以是中央处理器(英文:central processing unit, 缩写:CPU),网络处理器(英文:network processor,缩写:NP)或者CPU 和NP的组合。The processor 401 may be a central processing unit (English: central processing unit, abbreviation: CPU), a network processor (English: network processor, abbreviation: NP), or a combination of CPU and NP.

其中,处理器401还可以进一步包括硬件芯片。上述硬件芯片可以是 专用集成电路(英文:application-specific integrated circuit,缩写:ASIC), 可编程逻辑器件(英文:programmable logic device,缩写:PLD)或其组合。 上述PLD可以是复杂可编程逻辑器件(英文:complex programmable logic device,缩写:CPLD),现场可编程逻辑门阵列(英文:field-programmable gate array,缩写:FPGA),通用阵列逻辑(英文:generic arraylogic,缩写: GAL)或其任意组合。The processor 401 may further include a hardware chip. The above hardware chip may be an application-specific integrated circuit (English: application-specific integrated circuit, abbreviation: ASIC), a programmable logic device (English: programmable logic device, abbreviation: PLD) or a combination thereof. The above-mentioned PLD may be a complex programmable logic device (English: complex programmable logic device, abbreviation: CPLD), a field programmable gate array (English: field-programmable gate array, abbreviation: FPGA), a general array logic (English: generic arraylogic , abbreviation: GAL) or any combination thereof.

可选地,存储器404还用于存储程序指令。处理器401可以调用程序 指令,实现如本申请实施例1中提供的基于虚拟节点的多站协同干扰定位方 法。Optionally, memory 404 is also used to store program instructions. The processor 401 can invoke program instructions to implement the virtual node-based multi-station cooperative interference location method provided in Embodiment 1 of the present application.

本发明实施例还提供一种计算机可读存储介质,计算机可读存储介质 上存储有计算机可执行指令,该计算机可执行指令可执行上述任意方法实 施例中的基于虚拟节点的多站协同干扰定位方法。其中,所述存储介质可为 磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)、随机存储记忆 体(Random Access Memory,RAM)、快闪存储器(FlashMemory)、硬盘 (Hard Disk Drive,缩写:HDD)或固态硬盘(Solid-State Drive,SSD)等;所述存储介质还可以包括上述种类的存储器的组合。Embodiments of the present invention further provide a computer-readable storage medium, where computer-executable instructions are stored on the computer-readable storage medium, and the computer-executable instructions can execute the virtual node-based multi-station cooperative interference positioning in any of the foregoing method embodiments method. Wherein, the storage medium may be a magnetic disk, an optical disk, a read-only memory (Read-Only Memory, ROM), a random access memory (Random Access Memory, RAM), a flash memory (FlashMemory), a hard disk (Hard Disk) Drive, abbreviation: HDD) or solid-state drive (Solid-State Drive, SSD), etc.; the storage medium may also include a combination of the above-mentioned types of memories.

本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、 或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施 例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个 或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不 限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的 形式。As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序 产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流 程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中 的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专 用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个 机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产 生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方 框中指定的功能的装置。The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each process and/or block in the flowchart and/or block diagrams, and combinations of processes and/or blocks in the flowchart and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing device to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing device produce Means for implementing the functions specified in a flow or flow of a flowchart and/or a block or blocks of a block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理 设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存 储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个 流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture comprising instruction means, the instructions The apparatus implements the functions specified in the flow or flow of the flowcharts and/or the block or blocks of the block diagrams.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上, 使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现 的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流 程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能 的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing device to cause a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process such that The instructions provide steps for implementing the functions specified in the flow or blocks of the flowcharts and/or the block or blocks of the block diagrams.

显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方 式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可 以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予 以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保 护范围之中。Obviously, the above-mentioned embodiments are only examples for clear description, rather than limiting the implementation manner. For those of ordinary skill in the art, on the basis of the above description, other different forms of changes or modifications can also be made. There is no need and cannot be exhaustive of all implementations here. And the obvious changes or changes derived from this are still within the protection scope of the present invention.

Claims (10)

1.一种基于虚拟节点的多站协同干扰定位方法,其特征在于,包括如下步骤:1. a kind of multi-station cooperative interference positioning method based on virtual node, is characterized in that, comprises the steps: 获取监测节点中三个节点作为定位节点,分别为参考节点、第二节点及第三节点,包括:Acquire three nodes in the monitoring nodes as positioning nodes, namely the reference node, the second node and the third node, including: 获取所述监测节点中接收待定位信号源信号的RSSI值按高低排序的前三个节点作为所述定位节点;Acquiring the first three nodes in the monitoring node that receive the RSSI values of the signal source signals to be located in order of high and low as the positioning nodes; 分别获取待定位信号源信号到达各所述定位节点的RSSI值;Respectively obtain the RSSI value of the signal source signal to be located to each of the positioning nodes; 将所述RSSI值最大的定位节点作为所述参考节点;Using the positioning node with the largest RSSI value as the reference node; 将另两个定位节点分别作为所述第二节点、第三节点;Using the other two positioning nodes as the second node and the third node respectively; 计算所述第二节点与所述参考节点到达待定位信号源的传输时差,第三节点与所述参考节点到达待定位信号源的传输时差;Calculate the transmission time difference between the second node and the reference node when they arrive at the signal source to be located, and the transmission time difference between the third node and the reference node when they arrive at the signal source to be located; 根据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源的传输时差及所述第二节点接收信号的方向,确定第二虚拟节点的坐标;Determine the coordinates of the second virtual node according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction in which the second node receives the signal; 根据所述第三节点的坐标、所述第三节点与所述参考节点到达待定位信号源的传输时差及所述第三节点接收信号的方向,确定第三虚拟节点的坐标;Determine the coordinates of the third virtual node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction in which the third node receives the signal; 根据所述参考节点、所述第二节点与第二虚拟节点中之一、第三节点与第三虚拟节点中之一构建定位组合;constructing a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node; 计算各所述定位组合对应的待定位信号源的初始定位坐标;Calculate the initial positioning coordinates of the to-be-located signal source corresponding to each of the positioning combinations; 根据各定位节点接收待定位信号源信号的RSSI值分别计算各所述定位组合对应初始定位坐标的加权系数;According to the RSSI value of the signal source signal to be located received by each positioning node, the weighting coefficient of each of the positioning combinations corresponding to the initial positioning coordinates is calculated respectively; 根据所述初始定位坐标及与其对应的加权系数,计算所述待定位信号源的定位坐标。According to the initial positioning coordinates and the corresponding weighting coefficients, the positioning coordinates of the to-be-located signal source are calculated. 2.根据权利要求1所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,根据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源的传输时差及所述第二节点接收信号的方向,计算第二虚拟节点的坐标的步骤,包括:2. The virtual node-based multi-station cooperative interference positioning method according to claim 1, wherein, according to the coordinates of the second node, the transmission of the second node and the reference node to the signal source to be located The time difference and the direction in which the second node receives the signal, and the step of calculating the coordinates of the second virtual node, including: 根据所述第二节点与所述参考节点到达待定位信号源的传输时差,计算所述第二节点和所述参考节点之间的传输距离的差值;Calculate the difference in transmission distance between the second node and the reference node according to the transmission time difference between the second node and the reference node reaching the signal source to be located; 根据所述第二节点的坐标、所述第二节点和所述参考节点之间的传输距离的差值及所述第二节点接收信号的方向,确定所述第二虚拟节点的坐标。The coordinates of the second virtual node are determined according to the coordinates of the second node, the difference between the transmission distances between the second node and the reference node, and the direction in which the second node receives signals. 3.根据权利要求1所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,所述根据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源的传输时差及所述第二节点接收信号的方向,确定第二虚拟节点的坐标的步骤,具体包括:3 . The method for locating multi-station cooperative interference based on virtual nodes according to claim 1 , wherein the said second node and the reference node reach the signal source to be located according to the coordinates of the second node, the second node and the reference node. 4 . The step of determining the coordinates of the second virtual node specifically includes: 根据无线传播损耗计算公式,分别计算第二节点在遮挡环境下的第一传输损耗及在无遮挡环境下的第二传输损耗;According to the wireless propagation loss calculation formula, respectively calculate the first transmission loss of the second node in the occlusion environment and the second transmission loss in the unobstructed environment; 根据所述第二节点与所述参考节点到达待定位信号源的传输时差计算所述第二节点到达待定位信号源的第一距离;Calculate the first distance from the second node to the signal source to be located according to the transmission time difference between the second node and the reference node reaching the signal source to be located; 根据所述第二节点的RSSI值、所述第一传输损耗、所述第一距离及所述参考节点与待定位信号源的距离确定所述第二节点的RSSI值与所述第二虚拟节点的RSSI值的RSSI差值;Determine the RSSI value of the second node and the second virtual node according to the RSSI value of the second node, the first transmission loss, the first distance, and the distance between the reference node and the signal source to be located The RSSI difference of the RSSI value; 根据所述待定位信号源的功率、所述第二传输损耗及所述RSSI差值计算在无遮挡环境下所述第二节点与所述待定位信号源的第二距离;Calculate the second distance between the second node and the to-be-located signal source in an unobstructed environment according to the power of the to-be-located signal source, the second transmission loss and the RSSI difference; 根据所述参考节点与待定位信号源的距离、第二距离及所述第二节点的坐标及所述第二节点接收信号的方向确定所述第二虚拟节点的坐标。The coordinates of the second virtual node are determined according to the distance between the reference node and the signal source to be located, the second distance, the coordinates of the second node, and the direction in which the second node receives signals. 4.根据权利要求1所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,所述根据第三节点的坐标、所述第三节点与所述参考节点到达待定位信号源的传输时差所述第三节点接收信号的方向,计算第三虚拟节点的坐标的步骤,具体包括:4. The method for locating multi-station cooperative interference based on virtual nodes according to claim 1, characterized in that, according to the coordinates of the third node, the third node and the reference node arrive at the transmission of the signal source to be located. The step of calculating the coordinates of the third virtual node according to the time difference of the direction in which the third node receives the signal, specifically includes: 根据所述第三节点与所述参考节点到达待定位信号源的传输时差,计算所述第三节点和所述参考节点之间的传输距离的差值;Calculate the difference in the transmission distance between the third node and the reference node according to the transmission time difference between the third node and the reference node when they reach the signal source to be located; 根据所述第三节点的坐标、所述第三节点和所述参考节点之间的传输距离的差值及所述第三节点接收信号的方向,确定所述第三虚拟节点的坐标。The coordinates of the third virtual node are determined according to the coordinates of the third node, the difference between the transmission distances between the third node and the reference node, and the direction in which the third node receives signals. 5.根据权利要求1所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,所述根据第三节点的坐标、所述第三节点与所述参考节点到达待定位信号源的传输时差及所述第三节点接收信号的方向,计算第三虚拟节点的坐标的步骤,具体包括:5. The method for locating multi-station cooperative interference based on virtual nodes according to claim 1, characterized in that, according to the coordinates of the third node, the third node and the reference node arrive at the transmission of the signal source to be located. The time difference and the direction in which the third node receives the signal, and the step of calculating the coordinates of the third virtual node, specifically includes: 根据无线传播损耗计算公式,分别计算第三节点在遮挡环境下的第三传输损耗及在无遮挡环境下的第四传输损耗;According to the wireless propagation loss calculation formula, respectively calculate the third transmission loss of the third node in the occlusion environment and the fourth transmission loss in the unobstructed environment; 根据所述第三节点与所述参考节点到达待定位信号源的传输时差计算所述第三节点到达待定位信号源的第三距离;Calculate the third distance from the third node to the signal source to be located according to the transmission time difference between the third node and the reference node reaching the signal source to be located; 根据所述第三节点的RSSI值、所述第三传输损耗、所述第三距离及所述参考节点与待定位信号源的距离确定所述第三节点的RSSI值与所述第三虚拟节点的RSSI值的RSSI差值;Determine the RSSI value of the third node and the third virtual node according to the RSSI value of the third node, the third transmission loss, the third distance, and the distance between the reference node and the signal source to be located The RSSI difference of the RSSI value; 根据所述待定位信号源的功率、所述第四传输损耗及所述RSSI差值计算在无遮挡环境下所述第三节点与所述待定位信号源的第四距离;Calculate the fourth distance between the third node and the to-be-located signal source in an unobstructed environment according to the power of the to-be-located signal source, the fourth transmission loss and the RSSI difference; 根据所述参考节点与待定位信号源的距离、第四距离及所述第三节点的坐标及所述第三节点接收信号的方向确定所述第三虚拟节点的坐标。The coordinates of the third virtual node are determined according to the distance between the reference node and the signal source to be located, the fourth distance, the coordinates of the third node, and the direction in which the third node receives signals. 6.根据权利要求1-5任一所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,所述根据各定位节点接收待定位信号源信号的RSSI值分别计算所述定位组合对应初始定位坐标的加权系数的步骤,具体包括:6. The virtual node-based multi-station cooperative interference positioning method according to any one of claims 1 to 5, wherein, according to the RSSI value of the signal source signal to be located received by each positioning node, the corresponding positioning combination is calculated respectively. The steps of initially positioning the weighting coefficient of the coordinates include: 通过以下公式计算所述定位组合中各节点RSSI值的平均值:The average value of the RSSI values of each node in the positioning combination is calculated by the following formula:
Figure FDA0002453868150000041
Figure FDA0002453868150000041
其中,RSSIi表示第i个定位组合中各节点RSSI值的平均值,RSSIij表示第i个定位组合中第j个节点的RSSI值,n表示每组组合的节点的个数;Wherein, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, RSSI ij represents the RSSI value of the jth node in the ith positioning combination, and n represents the number of nodes in each combination; 通过以下公式计算所述定位组合对应的加权系数,:The weighting coefficient corresponding to the positioning combination is calculated by the following formula:
Figure FDA0002453868150000042
Figure FDA0002453868150000042
其中,wi表示第i个定位组合的加权系数,RSSIi表示第i个定位组合中各节点RSSI值的平均值,N表示定位组合的数量。Wherein, w i represents the weighting coefficient of the ith positioning combination, RSSI i represents the average value of the RSSI values of each node in the ith positioning combination, and N represents the number of positioning combinations.
7.根据权利要求6所述的基于虚拟节点的多站协同干扰定位方法,其特征在于,所述根据所述初始定位坐标及与其对应的加权系数,计算所述待定位信号源的定位坐标的步骤,通过以下公式计算所述待定位信号源的定位坐标:7. The virtual node-based multi-station cooperative interference positioning method according to claim 6, wherein, according to the initial positioning coordinates and the corresponding weighting coefficients, the calculation of the positioning coordinates of the to-be-located signal source is performed. Step, calculate the positioning coordinates of the signal source to be located by the following formula:
Figure FDA0002453868150000051
Figure FDA0002453868150000051
其中,wi表示第i个定位组合的加权系数,L表示待定位信号源信号的定位坐标,Li表示第i个定位组合的初始定位坐标。Wherein, wi represents the weighting coefficient of the ith positioning combination, L represents the positioning coordinates of the signal source signal to be located, and Li represents the initial positioning coordinates of the ith positioning combination.
8.一种基于虚拟节点的多站协同干扰定位系统,其特征在于,包括:8. A virtual node-based multi-station cooperative interference positioning system, characterized in that, comprising: 定位节点获取模块,用于获取监测节点中三个节点作为定位节点,分别为参考节点、第二节点及第三节点,包括:The positioning node obtaining module is used to obtain three nodes in the monitoring nodes as positioning nodes, which are the reference node, the second node and the third node respectively, including: 获取所述监测节点中接收待定位信号源信号的RSSI值按高低排序的前三个节点作为所述定位节点;Acquiring the first three nodes in the monitoring node that receive the RSSI values of the signal source signals to be located in order of high and low as the positioning nodes; 分别获取待定位信号源信号到达各所述定位节点的RSSI值;Respectively obtain the RSSI value of the signal source signal to be located to each of the positioning nodes; 将所述RSSI值最大的定位节点作为所述参考节点;Using the positioning node with the largest RSSI value as the reference node; 将另两个定位节点分别作为所述第二节点、第三节点;Using the other two positioning nodes as the second node and the third node respectively; 传输时差获取模块,用于计算所述第二节点与所述参考节点到达待定位信号源的传输时差,第三节点与所述参考节点到达待定位信号源的传输时差;a transmission time difference acquisition module, configured to calculate the transmission time difference between the second node and the reference node when they reach the signal source to be located, and the transmission time difference between the third node and the reference node when they arrive at the signal source to be located; 第二虚拟节点坐标确定模块,用于根据所述第二节点的坐标、所述第二节点与所述参考节点到达待定位信号源的传输时差及所述第二节点接收信号的方向,确定第二虚拟节点的坐标;The second virtual node coordinate determination module is configured to determine the second node according to the coordinates of the second node, the transmission time difference between the second node and the reference node reaching the signal source to be located, and the direction of the signal received by the second node. The coordinates of the two virtual nodes; 第三虚拟节点坐标确定模块,用于根据所述第三节点的坐标、所述第三节点与所述参考节点到达待定位信号源的传输时差及所述第三节点接收信号的方向,确定第三虚拟节点的坐标;The third virtual node coordinate determination module is configured to determine the third node according to the coordinates of the third node, the transmission time difference between the third node and the reference node reaching the signal source to be located, and the direction of the signal received by the third node. The coordinates of the three virtual nodes; 定位组合构建模块,用于根据所述参考节点、所述第二节点与第二虚拟节点中之一、第三节点与第三虚拟节点中之一构建定位组合;a positioning combination building module, configured to construct a positioning combination according to the reference node, one of the second node and the second virtual node, and one of the third node and the third virtual node; 初始定位坐标计算模块,用于计算各所述定位组合对应的待定位信号源的初始定位坐标;an initial positioning coordinate calculation module for calculating the initial positioning coordinates of the signal source to be positioned corresponding to each of the positioning combinations; 加权系数获取模块,用于根据各定位节点接收待定位信号源信号的RSSI值分别计算各所述定位组合对应初始定位坐标的加权系数;a weighting coefficient obtaining module, configured to calculate the weighting coefficients corresponding to the initial positioning coordinates of each of the positioning combinations according to the RSSI value of the signal source signal to be positioned received by each positioning node; 定位坐标计算模块,用于根据所述初始定位坐标及与其对应的加权系数,计算所述待定位信号源的定位坐标。The positioning coordinate calculation module is configured to calculate the positioning coordinates of the signal source to be located according to the initial positioning coordinates and the corresponding weighting coefficients. 9.一种计算机设备,其特征在于,包括:至少一个处理器,以及与所述至少一个处理器通信连接的存储器,其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行上述权利要求1-7中任一所述的基于虚拟节点的多站协同干扰定位方法。9. A computer device, comprising: at least one processor, and a memory communicatively connected to the at least one processor, wherein the memory stores instructions executable by the at least one processor, The instructions are executed by the at least one processor, so that the at least one processor executes the method for locating multi-station cooperative interference based on a virtual node according to any one of the preceding claims 1-7. 10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于使所述计算机执行上述权利要求1-7中任一所述的基于虚拟节点的多站协同干扰定位方法。10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores computer instructions, the computer instructions are used to cause the computer to Multi-station cooperative interference localization method for virtual nodes.
CN201810870985.1A 2018-08-02 2018-08-02 Multi-station cooperative interference positioning method and system based on virtual nodes Active CN109218961B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810870985.1A CN109218961B (en) 2018-08-02 2018-08-02 Multi-station cooperative interference positioning method and system based on virtual nodes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810870985.1A CN109218961B (en) 2018-08-02 2018-08-02 Multi-station cooperative interference positioning method and system based on virtual nodes

Publications (2)

Publication Number Publication Date
CN109218961A CN109218961A (en) 2019-01-15
CN109218961B true CN109218961B (en) 2020-07-28

Family

ID=64988124

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810870985.1A Active CN109218961B (en) 2018-08-02 2018-08-02 Multi-station cooperative interference positioning method and system based on virtual nodes

Country Status (1)

Country Link
CN (1) CN109218961B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111999696B (en) * 2020-08-02 2023-07-04 中国人民解放军海军工程大学 A multi-platform direction-finding cross-localization optimization method
CN112698267B (en) * 2020-12-05 2023-03-31 中海电信有限公司 Interference source testing and positioning method
CN114630416B (en) * 2020-12-14 2023-03-21 腾讯科技(深圳)有限公司 Method, device and equipment for positioning mobile terminal

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1862286A (en) * 2006-06-15 2006-11-15 北京邮电大学 Method for precisely positioning sensor node
CN101576615A (en) * 2008-05-05 2009-11-11 北京银易通网络科技有限公司 System and method model for hybrid positioning of WSN
CN102158956A (en) * 2011-03-08 2011-08-17 哈尔滨工业大学 Improved weighting trilateral positioning method based on RSSI (received signal strength indicator) in wireless sensor network
CN104039011A (en) * 2014-07-02 2014-09-10 保定市天河电子技术有限公司 Positioning method and device
CN104159292A (en) * 2013-05-14 2014-11-19 山东交通学院 Three-dimensional wireless sensing network positioning method based on virtual center node
JP2017067569A (en) * 2015-09-29 2017-04-06 株式会社Nttドコモ Terminal device and positioning program
CN107271957A (en) * 2017-06-20 2017-10-20 南京理工大学 Indoor 3-D positioning method based on TDOA and TOA
CN107656243A (en) * 2017-08-25 2018-02-02 天津大学 Combine DOA/TOA oceans multi-path environment localization method in inhomogeneous medium
CN108051779A (en) * 2017-10-13 2018-05-18 西安电子科技大学 A kind of positioning node preferred method towards TDOA

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7899006B2 (en) * 2006-12-05 2011-03-01 Zebra Enterprise Solutions Corp. Location system for wireless local area network (WLAN) using RSSI and time difference of arrival (TDOA) processing

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1862286A (en) * 2006-06-15 2006-11-15 北京邮电大学 Method for precisely positioning sensor node
CN101576615A (en) * 2008-05-05 2009-11-11 北京银易通网络科技有限公司 System and method model for hybrid positioning of WSN
CN102158956A (en) * 2011-03-08 2011-08-17 哈尔滨工业大学 Improved weighting trilateral positioning method based on RSSI (received signal strength indicator) in wireless sensor network
CN104159292A (en) * 2013-05-14 2014-11-19 山东交通学院 Three-dimensional wireless sensing network positioning method based on virtual center node
CN104039011A (en) * 2014-07-02 2014-09-10 保定市天河电子技术有限公司 Positioning method and device
JP2017067569A (en) * 2015-09-29 2017-04-06 株式会社Nttドコモ Terminal device and positioning program
CN107271957A (en) * 2017-06-20 2017-10-20 南京理工大学 Indoor 3-D positioning method based on TDOA and TOA
CN107656243A (en) * 2017-08-25 2018-02-02 天津大学 Combine DOA/TOA oceans multi-path environment localization method in inhomogeneous medium
CN108051779A (en) * 2017-10-13 2018-05-18 西安电子科技大学 A kind of positioning node preferred method towards TDOA

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A virtual TDOA localization scheme of Chinese DTMB signal in radio monitoring networks;Liu Shang等;《China Communications 》;20151225;全文 *
An energy efficient localization algorithm for wireless sensor networks using a mobile anchor node;Baoli Zhang;《 2008 International Conference on Information and Automation》;20080826;全文 *
一种新型三维传感器网络质心定位算法;刘清等;《微电子学与计算机》;20140531;全文 *
基于虚拟节点的BP无线传感器网络定位算法;孙凯等;《传感器与微系统》;20100930;全文 *

Also Published As

Publication number Publication date
CN109218961A (en) 2019-01-15

Similar Documents

Publication Publication Date Title
CN107659893B (en) Error compensation method and device, electronic equipment and readable storage medium
WO2021203871A1 (en) Cooperative positioning method and apparatus, device, and storage medium
EP2769233B1 (en) Time of arrival based wireless positioning system
US20080080441A1 (en) Localization apparatus for recognizing location of node in sensor network and method thereof
US9660740B2 (en) Signal strength distribution establishing method and wireless positioning system
US20130010617A1 (en) Relative position determination of wireless network devices
US8824325B2 (en) Positioning technique for wireless communication system
CN106772229B (en) Indoor positioning method and related equipment
EP2759176A1 (en) Time of arrival based positioning system
CN109218961B (en) Multi-station cooperative interference positioning method and system based on virtual nodes
WO2013043664A1 (en) Hybrid positioning system based on time difference of arrival (tdoa) and time of arrival (toa)
EP3391073A1 (en) Positioning of mobile device
WO2017071271A1 (en) Positioning method and device
WO2013043685A1 (en) Positioning system based on time - difference - of -arrival (tdoa)
WO2014168635A1 (en) Determining an angle of direct path of a signal
CN113923772A (en) Base station calibration method and device, computer equipment and storage medium
WO2014164101A1 (en) Method to obtain accurate vertical component estimates in 3d positioning
CN106885567B (en) A kind of inertial navigation cooperative positioning method and positioning device
CN112540347A (en) Method and device for judging distance of sound source, terminal equipment and storage medium
CN109996205B (en) Sensor data fusion method and device, electronic equipment and storage medium
CN108834053B (en) Positioning method, device and equipment
US20170038454A1 (en) Method and system for estimating signal generation position based on signal strength
CN114521017A (en) Method and device for positioning, electronic equipment and storage medium
Khalaf-Allah An extended closed-form least-squares solution for three-dimensional hyperbolic geolocation
KR20190043023A (en) The method for calculating ambiguity probability of the direction of arrival

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