CN106454750A - Multi-region indoor safety positioning method based on compressed sensing technology - Google Patents
Multi-region indoor safety positioning method based on compressed sensing technology Download PDFInfo
- Publication number
- CN106454750A CN106454750A CN201611047225.8A CN201611047225A CN106454750A CN 106454750 A CN106454750 A CN 106454750A CN 201611047225 A CN201611047225 A CN 201611047225A CN 106454750 A CN106454750 A CN 106454750A
- Authority
- CN
- China
- Prior art keywords
- node
- anchor
- rssi
- anchor node
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005516 engineering process Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000012544 monitoring process Methods 0.000 claims abstract description 10
- 238000011084 recovery Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 5
- 230000004807 localization Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract description 3
- 238000004891 communication Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000012795 verification Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 206010020400 Hostility Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/02—Systems for determining distance or velocity not using reflection or reradiation using radio waves
- G01S11/06—Systems for determining distance or velocity not using reflection or reradiation using radio waves using intensity measurements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W84/00—Network topologies
- H04W84/18—Self-organising networks, e.g. ad-hoc networks or sensor networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种基于压缩感知技术的多区域室内安全定位方法,该定位方法利用压缩感知的理论,包括离线、安全监测和在线三个阶段;1)在待定位节点所在设定区域内,按照等边三角形部署锚节点,并依据待定位节点距离锚节点之间的距离,构建待定位节点的RSSI安全密钥矩阵;利用待定位节点接收的实际RSSI值与RSSI安全密钥矩阵中的对应值进行比较,剔除恶意锚节点;待定位节点实时接收感知到的锚节点信息绘制重叠区域,利用信号恢复算法对正交化后的稀疏信号进行还原,进行精确的位置估算,获得待定位节点的定位信息。通过建立有自适应能力的随机数字典剔除掉恶意的定位信息,从而保证定位的安全。
The present invention provides a multi-region indoor safe positioning method based on compressed sensing technology. The positioning method utilizes the theory of compressed sensing and includes three stages: offline, safety monitoring and online; 1) within the set area where the node to be positioned is located, Deploy the anchor nodes according to an equilateral triangle, and construct the RSSI security key matrix of the node to be located according to the distance between the node to be located and the anchor node; use the actual RSSI value received by the node to be located and the correspondence in the RSSI security key matrix Values are compared to remove malicious anchor nodes; the node to be positioned receives the perceived anchor node information in real time to draw the overlapping area, and the signal recovery algorithm is used to restore the orthogonalized sparse signal to perform accurate position estimation and obtain the location of the node to be positioned positioning information. By establishing an adaptive random number dictionary, malicious positioning information is eliminated, so as to ensure the security of positioning.
Description
技术领域technical field
本发明属于无线传感器网络领域,特别涉及一种基于压缩感知技术的多区域室内安全定位方法。The invention belongs to the field of wireless sensor networks, in particular to a multi-region indoor safe positioning method based on compressed sensing technology.
背景技术Background technique
无线传感器网络室内定位近几年发展很迅速,其中包括智能感知技术,嵌入式计算技术和通信技术的应用提升了室内定位的精度拓宽了室内定位的应用范围。Wireless sensor network indoor positioning has developed rapidly in recent years, including the application of intelligent sensing technology, embedded computing technology and communication technology to improve the accuracy of indoor positioning and broaden the application range of indoor positioning.
室内定位在最近几年已经成功应用于各种领域,包括博物馆导航、仓库导航、停车场导航等。但是目前该领域的研究主要集中于如何提高定位精确度和能量有效性这方面,对无线传感器网络在室内定位过程中的安全性问题考虑不足。由于无线传感器节点是开放性部署的,导致它们很容易被攻击者俘获,从而定位所依据的信号的物理属性很容易被外界篡改。又因为信号的传输过程是在完全开放的环境中进行的,这导致了节点的信号在传输过程中容易被截获或者被篡改。无线传感器网络的室内定位过程遭受来自内部或者外部的攻击所产生的错误定位结果可能导致网络功能失效和监测结果出错,进而破坏网络应用的有效性。因此如何在有敌对可能的无线传感器网络应用中实现节点的安全定位,是一个必须解决的问题。Indoor positioning has been successfully applied in various fields in recent years, including museum navigation, warehouse navigation, parking lot navigation, etc. However, the current research in this field is mainly focused on how to improve the positioning accuracy and energy efficiency, and insufficient consideration is given to the security of wireless sensor networks in the process of indoor positioning. Since wireless sensor nodes are deployed openly, they are easily captured by attackers, and the physical properties of the signals based on positioning are easily tampered by the outside world. And because the signal transmission process is carried out in a completely open environment, this leads to the node's signal being easily intercepted or tampered with during transmission. The indoor positioning process of the wireless sensor network suffers from internal or external attacks, and the wrong positioning results may lead to the failure of network functions and errors in monitoring results, thereby destroying the effectiveness of network applications. Therefore, how to realize the safe positioning of nodes in the application of wireless sensor network with possible hostility is a problem that must be solved.
发明内容Contents of the invention
本发明所要解决的技术问题是,针对无线传感器网络室内定位的安全性问题,提供一种基于压缩感知技术的多区域室内安全定位的方法,保证定位精度的基础上,排除恶意的定位信息,从而保证定位环境的安全性。The technical problem to be solved by the present invention is to provide a multi-region indoor safe positioning method based on compressed sensing technology for the security of wireless sensor network indoor positioning, and to eliminate malicious positioning information on the basis of ensuring positioning accuracy, thereby Ensure the security of the positioning environment.
一种基于压缩感知技术的多区域室内安全定位方法,包括离线阶段、安全监测阶段和在线阶段,具体步骤如下:A multi-region indoor safe positioning method based on compressed sensing technology, including an offline stage, a safety monitoring stage and an online stage, the specific steps are as follows:
1)离线阶段1) Offline stage
在待定位节点所在设定区域内,按照等边三角形部署锚节点,并依据待定位节点距离锚节点之间的距离,构建待定位节点的RSSI安全密钥矩阵;In the set area where the node to be located is located, the anchor nodes are deployed according to an equilateral triangle, and the RSSI security key matrix of the node to be located is constructed according to the distance between the node to be located and the anchor node;
2)安全监测阶段2) Safety monitoring stage
利用待定位节点接收的实际RSSI值与RSSI安全密钥矩阵中的对应值进行比较,剔除恶意锚节点;Compare the actual RSSI value received by the node to be located with the corresponding value in the RSSI security key matrix to eliminate malicious anchor nodes;
3)在线阶段3) Online stage
待定位节点实时接收感知到的锚节点信息,将已感知的锚节点为圆心绘制圆形区域,得到所有锚节点对应的圆形区域形成的重叠区域;The node to be positioned receives the perceived anchor node information in real time, draws a circular area with the perceived anchor node as the center of the circle, and obtains the overlapping area formed by the circular areas corresponding to all anchor nodes;
依据重叠区域选取网格来构建在线测量矩阵,根据待定位节点感知到的锚节点信息构建在线观测矩阵;对在线测量矩阵和在线观测矩阵进行正交化预操作获得稀疏信号,利用信号恢复算法对正交化后的稀疏信号进行还原,进行精确的位置估算,获得待定位节点的定位信息。The grid is selected according to the overlapping area to construct the online measurement matrix, and the online observation matrix is constructed according to the anchor node information perceived by the node to be positioned; the orthogonal pre-operation is performed on the online measurement matrix and the online observation matrix to obtain sparse signals, and the signal recovery algorithm is used to The sparse signal after orthogonalization is restored, the precise position estimation is performed, and the positioning information of the node to be positioned is obtained.
在线测量矩阵Ψ中的元素依次为重叠区域中的网格节点的离线其维数是其中表示关联度为1的网格节点的数量;The elements in the online measurement matrix Ψ are in turn the offline Its dimension is in Indicates the number of grid nodes with an association degree of 1;
位于重叠区域中的网格节点的关联度即为1;The degree of association of grid nodes located in the overlapping area is 1;
表示第j个网格节点在第i个锚节点采样的离线RSS平均值; Indicates the offline RSS average value sampled by the j-th grid node at the i-th anchor node;
表示第j个网格节点接收到第i个锚节点的第τ个RSS采样值;如果未搜集到锚节点的信息,令 Indicates that the j-th grid node receives the τ-th RSS sample value of the i-th anchor node; if the information of the anchor node is not collected, let
在线观测矩阵为Φ, The online observation matrix is Φ,
Φ是M×L维的矩阵,Φ的每一行是一个1×L的向量,所有的元素满足φi,j∈{0,1},每行都只有一个元素为1,为1的元素所在列号表示选择的是L个锚节点中的第几个;如果φh,g=1,h=1,2…,M,g=1,2…,L,则表示选择的第h个锚节点是L个锚节点中的第g个,锚节点为选定锚节点,依据wi×RSSi,j,i=1,2,...,L,j=1,2,...,N的值,选取非零的锚节点;Φ is an M×L-dimensional matrix, each row of Φ is a 1×L vector, all elements satisfy φ i,j ∈{0,1}, and each row has only one element that is 1, where the element of 1 is The column number indicates which of the L anchor nodes is selected; if φ h, g = 1, h = 1, 2..., M, g = 1, 2..., L, it indicates the hth anchor selected The node is the gth among the L anchor nodes, and the anchor node is the selected anchor node, according to w i ×RSS i,j , i=1,2,...,L,j=1,2,... , the value of N, select a non-zero anchor node;
RSSi,j表示第j个网格节点在第i个锚节点采样的在线RSS平均值,wi表示第i个锚节点的权值,即每个锚节点在所有网格点采集RSS过程中的贡献值,Numk表示在第k个网格节点上是否感知到了该锚节点,即当在网格节点放置一测试传感节点时,锚节点是否在测试传感节点的通信范围之内,Numk∈{0,1}.RSS i, j represent the online RSS average value of the j grid node sampled at the i anchor node, w i represents the weight of the i anchor node, That is, the contribution value of each anchor node in the process of collecting RSS at all grid points, and Numk indicates whether the anchor node is perceived on the kth grid node, that is, when a test sensor node is placed on the grid node, the anchor node Whether the node is within the communication range of the test sensor node, Numk∈{0,1}.
进一步地,所述待定位节点的RSSI安全密钥矩阵为Num:Further, the RSSI security key matrix of the node to be located is Num:
其中,numi,j表示第i个待定位节点接收到的第j个锚节点发的随机数, 表示第i个待定位节点接收到第j个锚节点发送的RSSI值RSSi,j的阈值,若则的取值为k的取值范围为0-2;Among them, num i,j represents the random number sent by the jth anchor node received by the i-th node to be positioned, Indicates the threshold value of the RSSI value RSS i, j sent by the i-th node to be located received by the j-th anchor node, if but The value is The value range of k is 0-2;
表示距离锚节点为r1时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r 1 ;
表示距离锚节点为r2时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r2 ;
表示距离锚节点为r3时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r 3 ;
r1表示任意两个相邻锚节点之间的距离,r3为待定位节点能接收到RSSI信息的最大的范围半径,当待定位节点的与锚节点的距离超过了r3就将待定位节点接受到的RSSI值置零,r2为r1和r3的平均值;r 1 represents the distance between any two adjacent anchor nodes, and r 3 is the maximum range radius that the node to be positioned can receive RSSI information. When the distance between the node to be positioned and the anchor node exceeds r 3 , it will be positioned The RSSI value received by the node is set to zero, r 2 is the average value of r 1 and r 3 ;
i的取值范围为1-L,L表示待定位节点的个数;j表的取值范围为1-N,N表示锚节点的个数。The value range of i is 1-L, and L represents the number of nodes to be located; the value range of j table is 1-N, and N represents the number of anchor nodes.
进一步地,所述利用待定位节点接收的实际RSSI值与RSSI安全密钥矩阵中的对应值进行比较,剔除恶意锚节点的具体过程如下:Further, the actual RSSI value received by the node to be positioned is compared with the corresponding value in the RSSI security key matrix, and the specific process of removing malicious anchor nodes is as follows:
第i个待定位节点接收到第j个锚节点发送的随机数为numi,j,且感知到的RSSI值为RSSi,j,RSSi,j对应的阈值为 The i-th node to be positioned receives the random number num i,j sent by the j-th anchor node, and the perceived RSSI value is RSS i,j , and the threshold corresponding to RSS i,j is
判断随机数numi,j是否满足公式: Determine whether the random number num i, j satisfies the formula:
如果满足,则判定第i个待定位节点感知到的第j个锚节点的定位信息是正确的;If it is satisfied, it is determined that the location information of the jth anchor node perceived by the ith node to be positioned is correct;
如果不满足或者没有收到对应的随机数,则判断第i个待定位节点接受到的第j个锚节发送的定位信息是恶意的,在第i个待定位节点的定位过程中将不使用第j个锚节点的定位信息,删除该恶意锚节点发出的定位信息。If it is not satisfied or the corresponding random number is not received, it is judged that the positioning information sent by the j-th anchor node received by the i-th node to be positioned is malicious, and will not be used in the positioning process of the i-th node to be positioned. The positioning information of the jth anchor node, and delete the positioning information sent by the malicious anchor node.
进一步地,所述将已感知的锚节点为圆心绘制圆形区域时,圆形半径按照以下过程确定:Further, when drawing a circular area with the perceived anchor node as the center of the circle, the radius of the circle is determined according to the following process:
根据待定位节点接收到的RSSi,j取得的阈值得到对应的阈值下的半径等级ri,j,依据半径等级获取对应锚节点的圆形半径;The threshold obtained according to the RSS i,j received by the node to be positioned Obtain the radius level r i,j under the corresponding threshold, and obtain the circular radius of the corresponding anchor node according to the radius level;
若取值为则半径等级ri,j对应的半径为r1;like The value is Then the radius level r i,j corresponds to r 1 ;
若取值为则半径等级ri,j对应的半径为r2;like The value is Then the radius level r i,j corresponds to r 2 ;
若取值为则半径等级ri,j对应的半径为r3。like The value is Then the radius level r i,j corresponds to r 3 .
有益效果Beneficial effect
本发明的基于压缩感知技术的无线传感网络室内安全定位方法,该定位方法利用压缩感知的理论,包括离线、安全监测和在线三个阶段。通过建立有自适应能力的随机数字典剔除掉恶意的定位信息,从而保证定位的安全。与传统的定位方法相比,本专利的优点在于:The wireless sensor network indoor safe positioning method based on compressed sensing technology of the present invention utilizes the theory of compressed sensing and includes three stages of offline, safety monitoring and online. By establishing an adaptive random number dictionary, malicious positioning information is eliminated, so as to ensure the security of positioning. Compared with traditional positioning methods, the advantages of this patent are:
1)设计了一种多区域叠加的机制,提出了基于锚节点传输半径的有自适应能力的随机数密钥算法,来剔除掉室内定位环境中的恶意定位信息;1) A multi-area superposition mechanism is designed, and an adaptive random number key algorithm based on the transmission radius of the anchor node is proposed to eliminate malicious positioning information in the indoor positioning environment;
2)在在线线阶段进行了一步预处理的工作,对离线的RSS信息进行加权,来平衡每个锚节点的影响;2) A step of preprocessing is carried out in the online stage, and the offline RSS information is weighted to balance the influence of each anchor node;
3)引进了CS算法来提高定位的效率;3) The CS algorithm is introduced to improve the efficiency of positioning;
利用多区域覆盖和CS算法保证了定位的精确度,重点引入了实时校验机制,确保了定位环境的安全性。Using multi-area coverage and CS algorithm to ensure the accuracy of positioning, emphatically introducing a real-time verification mechanism to ensure the security of the positioning environment.
附图说明Description of drawings
图1为定位的场景图;Figure 1 is a scene diagram of positioning;
图2为本发明的定位方法流程图;Fig. 2 is the flow chart of positioning method of the present invention;
图3为不同算法在有攻击的环境中平均定位误差;Figure 3 shows the average positioning error of different algorithms in an attacked environment;
图4为不同的安全定位算法在遭受不同攻击时的定位误差。Figure 4 shows the positioning errors of different security positioning algorithms under different attacks.
具体实施方式detailed description
以下将结合附图和具体实施例对本发明做进一步详细说明:The present invention will be described in further detail below in conjunction with accompanying drawing and specific embodiment:
实施例1:Example 1:
一种基于压缩感知技术的多区域室内安全定位方法,包括离线阶段、安全监测阶段和在线阶段,具体步骤如下:A multi-region indoor safe positioning method based on compressed sensing technology, including an offline stage, a safety monitoring stage and an online stage, the specific steps are as follows:
一、离线阶段1. Offline stage
在待定位节点所在设定区域内,按照等边三角形部署锚节点,并依据待定位节点距离锚节点之间的距离,构建待定位节点的RSSI安全密钥矩阵;In the set area where the node to be located is located, the anchor nodes are deployed according to an equilateral triangle, and the RSSI security key matrix of the node to be located is constructed according to the distance between the node to be located and the anchor node;
1)建立离线测量矩阵:1) Establish offline measurement matrix:
定位区域内网格节点的总数为N,锚节点的个数为L;在每个网格上,进行q次采样,搜集网格点感知到的所有锚节点的RSS(接收信号强度)值,获取一个初始L×N维的离线测量矩阵 The total number of grid nodes in the positioning area is N, and the number of anchor nodes is L; on each grid, q times of sampling are performed to collect the RSS (received signal strength) values of all anchor nodes sensed by grid points, Obtain an initial L×N dimensional offline measurement matrix
其中表示第j个网格节点在第i个锚节点采样的离线RSS平均值;in Indicates the offline RSS average value sampled by the j-th grid node at the i-th anchor node;
表示网格节点j接收到锚节点i的第τ个RSS采样值;如果未搜集到锚节点的信息,令 Indicates that the grid node j has received the τth RSS sampling value of the anchor node i; if the information of the anchor node is not collected, let
2)加权锚节点:2) Weighted anchor nodes:
权值矩阵W=[w1,w2,…,wL]T,其中表示第i个锚节点的权值,即每个锚节点在所有网格节点采集RSS过程中的贡献值,Numk表示在第k个网格节点上是否感知到了该锚节点,即当在网格节点放置一测试传感节点时,锚节点是否在测试传感节点的通信范围之内,Numk∈{0,1};Weight matrix W=[w 1 ,w 2 ,…,w L ] T , where Indicates the weight of the i-th anchor node, that is, the contribution value of each anchor node in the RSS collection process of all grid nodes, and Numk indicates whether the anchor node is perceived on the k-th grid node, that is, when the grid node When the node places a test sensor node, whether the anchor node is within the communication range of the test sensor node, Numk∈{0,1};
3)建立离线指纹库:3) Establish offline fingerprint library:
参考节点j接收到测试点即锚节点i采样值的无偏估计表示为对于每个参考节点j,无偏估计矩阵表示为Δj=[Δ1,j,Δ2,j,…,ΔL,j]T;离线指纹库表示为其中(xj,yj)是参考节点j的坐标。The unbiased estimate of the sampled value of the anchor node i received by the reference node j is expressed as For each reference node j, the unbiased estimation matrix is expressed as Δ j = [Δ 1,j ,Δ 2,j ,…,Δ L,j ] T ; the offline fingerprint library is expressed as where (x j ,y j ) is the coordinate of reference node j.
二、安全监测阶段2. Safety monitoring stage
利用待定位节点接收的实际RSSI值与RSSI安全密钥矩阵中的对应值进行比较,剔除恶意锚节点;Compare the actual RSSI value received by the node to be located with the corresponding value in the RSSI security key matrix to eliminate malicious anchor nodes;
所述待定位节点的RSSI安全密钥矩阵为Num:The RSSI security key matrix of the node to be positioned is Num:
其中,numi,j表示第i个待定位节点接收到的第j个锚节点发的随机数, 表示第i个待定位节点接收到第j个锚节点发送的RSSI值RSSi,j的阈值,若则的取值为k的取值范围为0-2;Among them, num i, j represents the random number sent by the jth anchor node received by the i-th node to be positioned, Indicates the threshold value of the RSSI value RSS i,j sent by the i-th node to be located received by the j-th anchor node, if but The value is The value range of k is 0-2;
表示距离锚节点为r1时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r 1 ;
表示距离锚节点为r2时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r2 ;
表示距离锚节点为r3时,待定位节点接收到到的RSSI阈值; Indicates the RSSI threshold received by the node to be positioned when the distance from the anchor node is r 3 ;
r1表示任意两个相邻锚节点之间的距离,r3为待定位节点能接收到RSSI信息的最大的范围半径,当待定位节点的与锚节点的距离超过了r3就将待定位节点接受到的RSSI值置零,r2为r1和r3的平均值;r 1 represents the distance between any two adjacent anchor nodes, and r 3 is the maximum range radius that the node to be positioned can receive RSSI information. When the distance between the node to be positioned and the anchor node exceeds r 3 , it will be positioned The RSSI value received by the node is set to zero, and r 2 is the average value of r 1 and r 3 ;
i的取值范围为1-L,L表示待定位节点的个数;j表的取值范围为1-N,N表示锚节点的个数。The value range of i is 1-L, and L represents the number of nodes to be located; the value range of j table is 1-N, and N represents the number of anchor nodes.
所述利用待定位节点接收的实际RSSI值与RSSI安全密钥矩阵中的对应值进行比较,剔除恶意锚节点的具体过程如下:The actual RSSI value received by the node to be positioned is compared with the corresponding value in the RSSI security key matrix, and the specific process of removing malicious anchor nodes is as follows:
第i个待定位节点接收到第j个锚节点发送的随机数为numi,j,且感知到的RSSI值为RSSi,j,RSSi,j对应的阈值为 The i-th node to be positioned receives the random number num i,j sent by the j-th anchor node, and the perceived RSSI value is RSS i,j , and the threshold corresponding to RSS i,j is
判断随机数numi,j是否满足公式: Determine whether the random number num i, j satisfies the formula:
如果满足,则判定第i个待定位节点感知到的第j个锚节点的定位信息是正确的;If it is satisfied, it is determined that the location information of the jth anchor node perceived by the ith node to be positioned is correct;
如果不满足或者没有收到对应的随机数,则判断第i个待定位节点接受到的第j个锚节发送的定位信息是恶意的,在第i个待定位节点的定位过程中将不使用第j个锚节点的定位信息,删除该恶意锚节点发出的定位信息。If it is not satisfied or the corresponding random number is not received, it is judged that the positioning information sent by the j-th anchor node received by the i-th node to be positioned is malicious, and will not be used in the positioning process of the i-th node to be positioned. The positioning information of the jth anchor node, and delete the positioning information sent by the malicious anchor node.
三、在线阶段3. Online stage
待定位节点实时接收感知到的锚节点信息,将已感知的锚节点为圆心绘制圆形区域,得到所有锚节点对应的圆形区域形成的重叠区域;The node to be positioned receives the perceived anchor node information in real time, draws a circular area with the perceived anchor node as the center of the circle, and obtains the overlapping area formed by the circular areas corresponding to all anchor nodes;
依据重叠区域选取网格来构建在线测量矩阵,根据待定位节点感知到的锚节点信息构建在线观测矩阵;对在线测量矩阵和在线观测矩阵进行正交化预操作获得稀疏信号,利用信号恢复算法对正交化后的稀疏信号进行还原,进行精确的位置估算,获得待定位节点的定位信息。The grid is selected according to the overlapping area to construct the online measurement matrix, and the online observation matrix is constructed according to the anchor node information perceived by the node to be positioned; the orthogonal pre-operation is performed on the online measurement matrix and the online observation matrix to obtain sparse signals, and the signal recovery algorithm is used to The sparse signal after orthogonalization is restored, the precise position estimation is performed, and the positioning information of the node to be positioned is obtained.
在线测量矩阵Ψ中的元素依次为重叠区域中的网格节点的离线其维数是其中表示关联度为1的网格节点的数量;The elements in the online measurement matrix Ψ are in turn the offline Its dimension is in Indicates the number of grid nodes with an association degree of 1;
位于重叠区域中的网格节点的关联度即为1;The degree of association of grid nodes located in the overlapping area is 1;
表示第j个网格节点在第i个锚节点采样的离线RSS平均值; Indicates the offline RSS average value sampled by the j-th grid node at the i-th anchor node;
表示第j个网格节点接收到第i个锚节点的第τ个RSS采样值;如果未搜集到锚节点的信息,令 Indicates that the j-th grid node receives the τ-th RSS sample value of the i-th anchor node; if the information of the anchor node is not collected, let
在线观测矩阵为Φ, The online observation matrix is Φ,
Φ是M×L维的矩阵,Φ的每一行是一个1×L的向量,所有的元素满足φi,j∈{0,1},每行都只有一个元素为1,为1的元素所在列号表示选择的是L个锚节点中的第几个;如果φh,g=1,h=1,2…,M,g=1,2…,L,则表示选择的第h个锚节点是L个锚节点中的第g个,锚节点为选定锚节点,依据wi×RSSi,j,i=1,2,...,L,j=1,2,...,N的值,选取非零的锚节点;Φ is an M×L-dimensional matrix, each row of Φ is a 1×L vector, all elements satisfy φ i,j ∈{0,1}, and each row has only one element that is 1, where the element of 1 is The column number indicates which of the L anchor nodes is selected; if φ h , g =1, h=1,2...,M, g=1,2...,L, it indicates the hth anchor selected The node is the gth among the L anchor nodes, and the anchor node is the selected anchor node, according to w i ×RSS i,j , i=1,2,...,L,j=1,2,... , the value of N, select a non-zero anchor node;
RSSi,j表示第j个网格节点在第i个锚节点采样的在线RSS平均值,wi表示第i个锚节点的权值,即每个锚节点在所有网格点采集RSS过程中的贡献值,Numk表示在第k个网格节点上是否感知到了该锚节点,即当在网格节点放置一测试传感节点时,锚节点是否在测试传感节点的通信范围之内,Numk∈{0,1}.RSS i, j represent the online RSS average value of the j grid node sampled at the i anchor node, w i represents the weight of the i anchor node, That is, the contribution value of each anchor node in the process of collecting RSS at all grid points, and Numk indicates whether the anchor node is perceived on the kth grid node, that is, when a test sensor node is placed on the grid node, the anchor node Whether the node is within the communication range of the test sensor node, Numk∈{0,1}.
首先对Φ和Ψ进行正交化预操作,再利用信号恢复算法对稀疏信号进行还原;信号恢复采用凸优化算法中的l1-最小范式来处理,得到恢复后的稀疏信号 Firstly, the orthogonalization pre-operation is performed on Φ and Ψ, and then the signal recovery algorithm is used to perform the sparse signal Restore; the signal recovery is processed by the l 1 -minimum paradigm in the convex optimization algorithm, and the restored sparse signal is obtained
通过和节点的位置坐标进行加权,对待定位节点的位置进行更加精确的估算,表达式为: pass Weighted with the position coordinates of the node, the position of the node to be positioned is estimated more accurately, the expression is:
其中,px和py分别代表待定位节点最终的位置对应的横坐标和纵坐标;xk、yk和分别代表第k个网格的x坐标,y坐标和待定位节点出现在该网格的概率。Among them, px and py respectively represent the abscissa and ordinate corresponding to the final position of the node to be positioned; x k , y k and Represent the x-coordinate and y-coordinate of the kth grid and the probability that the node to be located appears in the grid, respectively.
在正交化预操作和信号恢复步骤中,假设T是测量值y的正交化预操作,对y进行正交化,得到正交化后的测量值y′=Ty;定义P=ΦΨ,Q=orth(PT)T,orth(A)表示对矩阵P的规范正交化操作,为P的伪逆矩阵,则定位问题描述为以下l1-最小范式模型:In the steps of orthogonalization pre-operation and signal recovery, assuming that T is the orthogonalization pre-operation of the measured value y, y is orthogonalized to obtain the measured value after orthogonalization y′=Ty; definition P=ΦΨ, Q=orth(P T ) T , orth(A) represents the normalized orthogonalization operation on the matrix P, is the pseudo-inverse matrix of P, then the localization problem is described as the following l 1 -minimum normal form model:
其中,θ是表示待定位节点可能出现的一组位置的向量,表示对于特定的误差值ε′的满足等式条件的特定的θ值,Z是计算时用到的目标矩阵函数,有Z=y,y为测量向量,是指在线阶段待定位节点采集的感知到的所有锚节点的RSS值,ε′是设定的误差,在区间[0,0.1]范围内取值;计算不等式z-Qθ≤ε'。Among them, θ is a vector representing a set of positions where the node to be located may appear, Indicates a specific value of θ that satisfies the condition of the equation for a specific error value ε′, Z is the target matrix function used in the calculation, Z=y, and y is the measurement vector, which refers to the perception collected by the node to be positioned in the online phase The RSS values of all the anchor nodes arrived, ε' is the set error, and the value is in the interval [0,0.1]; calculate the inequality z-Qθ≤ε'.
构建关联矩阵可以看做是一个聚类问题。多区域叠加算法被设计来解决该问题。多区域叠加算法的思想是利用叠加技术缩小定位区域。区域叠加算法在解决聚类问题上有很多优势。区域叠加是非迭代的。这样定位速度比一般的聚类算法更快;区域叠加算法很容易实现。在线阶段动态获取的RSS信息比离线数据库更能反映当前的环境状况,因此定位更加精确;区域叠加算法在线阶段对定位区域的缩小性能较好。因此,该算法比传统的聚类算法更有效。Building an incidence matrix can be seen as a clustering problem. A multi-region overlay algorithm is designed to solve this problem. The idea of the multi-area overlay algorithm is to use the overlay technology to narrow the positioning area. The region overlay algorithm has many advantages in solving clustering problems. Area overlays are non-iterative. In this way, the positioning speed is faster than the general clustering algorithm; the area superposition algorithm is easy to implement. The RSS information obtained dynamically in the online stage can better reflect the current environmental conditions than the offline database, so the positioning is more accurate; the area superposition algorithm has a better performance in shrinking the positioning area in the online stage. Therefore, this algorithm is more efficient than traditional clustering algorithms.
在聚类思想中,找到最小且最精确的定位区域是最关键的问题。区域叠加算法利用重叠机制来获取最优的区域。重叠的基本思想是将通信区域内的候选锚节点的通信半径区域全部叠加。候选节点通过最适函数来进行选取。在该算法中,候选节点有更高的在线RSS值和离线权值。In clustering thinking, finding the smallest and most accurate localization area is the most critical issue. The region overlay algorithm uses the overlap mechanism to obtain the optimal region. The basic idea of overlapping is to superimpose all the communication radius areas of the candidate anchor nodes in the communication area. Candidate nodes are selected through the optimal function. In this algorithm, candidate nodes have higher online RSS value and offline weight.
根据接收到的在线RSSI矩阵,得到非零的RSSI值的锚节点的标号,根据这些标号将相对应的传输半径矩阵进行融合得到待定位节点的位置估计值。最后根据得到的半径值进行区域叠加覆盖。接着正交化生成稀疏矩阵,再利用压缩感知和L1-最小范式回复稀疏信号,获取带定位节点的坐标。According to the received online RSSI matrix, the labels of the anchor nodes with non-zero RSSI values are obtained, and the corresponding transmission radius matrix is fused according to these labels to obtain the position estimation value of the node to be located. Finally, the area is superimposed and covered according to the obtained radius value. Then the sparse matrix is generated by orthogonalization, and then the sparse signal is recovered by using compressed sensing and L1-minimum paradigm, and the coordinates of the nodes with positioning are obtained.
图1为本发明的应用场景图。在场上定位区域中共部署了5个待定位节点,用红色的圆点表示。25个锚节点用蓝色三角形表示。4个恶意节点,用黄色的三角形表示。每个待定位节点感知周围的锚节点信息,并从数据库中下载锚节点的相关信息。各锚节点利用相关信息和定位方法确定自己的位置。恶意节点通过发送错误的RSS信息误导定位。FIG. 1 is a diagram of an application scenario of the present invention. A total of 5 nodes to be located are deployed in the positioning area on the field, indicated by red dots. The 25 anchor nodes are represented by blue triangles. 4 malicious nodes, represented by yellow triangles. Each node to be positioned perceives the surrounding anchor node information, and downloads the relevant information of the anchor node from the database. Each anchor node uses relevant information and positioning methods to determine its own position. Malicious nodes mislead positioning by sending wrong RSS information.
图2为本发明定位方法流程图。定位分为离线、安全监测和在线阶段。离线阶段构建初始数据库。再对离现阶段感知到的锚节点进行加权,使得的每个锚节点在定位过程中达到均衡。然后并建立安全密钥,设定锚节点的发送功率,根据发送功率不同给出不同范围的随机数。锚节点在发送RSS信息的同时发送随机数。每个待定位节点都能接收到不同的随机数,作为在线价段进行安全验证的密钥。在现阶段,根据待定位节点接受到的随机数集给出相对应的半径值。验证随机数和接收到的RSS信息是否对等。剔除掉恶意的定位信息。最后根据得到的半径值进行区域叠加覆盖。接着正交化生成稀疏矩阵,再利用压缩感知和L1-最小范式回复稀疏信号,获取带定位节点的坐标。Fig. 2 is a flowchart of the positioning method of the present invention. Positioning is divided into offline, safety monitoring and online phases. The initial database is built in the offline phase. Then weight the anchor nodes perceived from the current stage, so that each anchor node reaches a balance in the positioning process. Then establish a security key, set the sending power of the anchor node, and give random numbers in different ranges according to the sending power. Anchor nodes send random numbers while sending RSS information. Each node to be positioned can receive a different random number as a key for security verification in the online price segment. At this stage, the corresponding radius value is given according to the random number set received by the node to be positioned. Verify that the random number is equal to the received RSS information. Eliminate malicious positioning information. Finally, the area is superimposed and covered according to the obtained radius value. Then the sparse matrix is generated by orthogonalization, and then the sparse signal is recovered by using compressed sensing and L1-minimum paradigm, and the coordinates of the nodes with positioning are obtained.
图3为几种算法在遭受攻击时的定位误差,横坐标代表了恶意锚节点的个数,纵坐标代表定位的平均误差。CS_NSL算法是一种单一密钥校验算法。锚节点在发送定位信息的同时向所有定位区域内广播同一的密钥。当待定位节点接收到的密钥不是初始密钥时,则证明该定位信息是恶意的。CS_SL是一种基于压缩感知的多区域网格室内定位算法,它是利用网格将区域划分成为一个一个的小区域,提高定位的精确度。同时利用了压缩感知的方法,选取活性比较强锚节点进行定位,排除了部分恶意节点的干扰。其中当恶意锚节点数目逐渐增加时,CS_SL的定位误差急剧上升,而当恶意节点超过两个以后CS_NSL的定位误差也有明显的增加。而本专利提出的多区域的密钥验证算法的定位误差几乎没有变化。因为本发明的校验算法是实时性的动态验证的,几乎所有的恶意定位信息都可以被检验并排除。Figure 3 shows the positioning errors of several algorithms when they are under attack. The abscissa represents the number of malicious anchor nodes, and the ordinate represents the average positioning error. The CS_NSL algorithm is a single-key verification algorithm. The anchor node broadcasts the same key to all positioning areas while sending positioning information. When the key received by the node to be located is not the initial key, it proves that the positioning information is malicious. CS_SL is a multi-area grid indoor positioning algorithm based on compressed sensing. It uses grids to divide the area into small areas one by one to improve the accuracy of positioning. At the same time, the method of compressed sensing is used to select anchor nodes with relatively strong activity for positioning, and the interference of some malicious nodes is eliminated. Among them, when the number of malicious anchor nodes gradually increases, the positioning error of CS_SL increases sharply, and when the number of malicious anchor nodes exceeds two, the positioning error of CS_NSL also increases significantly. However, the positioning error of the multi-region key verification algorithm proposed in this patent hardly changes. Because the verification algorithm of the present invention is dynamically verified in real time, almost all malicious positioning information can be verified and eliminated.
图4为CS_NSL算法和本发明的算法在遭受外部攻击和内部攻击时的平均定位误差,其中选取了洪范攻击和伪造攻击作为内部和外部攻击的代表,对算法的安全性进行验证。图中,‘-*-’代表了本发明的定位方法(CSMR_SL),将它与其他两种算法CS_SL和CS_NSL进行对比。两种攻击结果都显示本发明所述的算法防御攻击的能力更好。外部攻击无法获取密钥,但是可以不断的发送错误的密钥去猜测安全密钥。内部攻击可以通过俘获锚节点,再发送错误的定位信息,从而来破坏定位过程。当密钥比较单一时,外部攻击可以通过不断试发送密钥,猜中算法的密钥,使得安全算法失效。而基于传输半径值的密钥实时变化,几乎无法猜中。而内部攻击中,根据密钥生成算法。如果接收到的随机数密钥不超过用于密钥生成的半径阈值,即定位信息正常。本发明所述的方法与传统的单一密钥算法相比,将以前的覆盖范围划分成几个子范围,在进行一一验证,很大程度上提高了算法的有效性。Fig. 4 is the average positioning error of the CS_NSL algorithm and the algorithm of the present invention when subjected to external attacks and internal attacks, wherein Hongfan attack and forgery attack are selected as representatives of internal and external attacks to verify the security of the algorithm. In the figure, '-*-' represents the positioning method (CSMR_SL) of the present invention, which is compared with other two algorithms CS_SL and CS_NSL. Both attack results show that the algorithm described in the present invention has a better ability to defend against attacks. External attacks cannot obtain the key, but can continuously send wrong keys to guess the security key. Internal attacks can destroy the positioning process by capturing anchor nodes and sending wrong positioning information. When the key is relatively single, an external attack can guess the key of the algorithm by continuously trying to send the key, making the security algorithm invalid. However, the key based on the transmission radius value changes in real time, making it almost impossible to guess. In the internal attack, according to the key generation algorithm. If the received random number key does not exceed the radius threshold used for key generation, the location information is normal. Compared with the traditional single-key algorithm, the method of the invention divides the previous coverage range into several sub-ranges and performs one-by-one verification, which greatly improves the effectiveness of the algorithm.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611047225.8A CN106454750B (en) | 2016-11-23 | 2016-11-23 | A kind of multizone indoor security localization method based on compressed sensing technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611047225.8A CN106454750B (en) | 2016-11-23 | 2016-11-23 | A kind of multizone indoor security localization method based on compressed sensing technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106454750A true CN106454750A (en) | 2017-02-22 |
CN106454750B CN106454750B (en) | 2019-05-17 |
Family
ID=58219857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611047225.8A Active CN106454750B (en) | 2016-11-23 | 2016-11-23 | A kind of multizone indoor security localization method based on compressed sensing technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106454750B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107046680A (en) * | 2017-05-22 | 2017-08-15 | 湖南大学 | An Indoor Positioning Method Based on Trusted Location Fingerprint |
CN108519577A (en) * | 2018-03-12 | 2018-09-11 | 中国矿业大学(北京) | A Distributed Localization Method Based on Compressed Sensing TOA Feature Signal Fingerprint Database |
CN108614237A (en) * | 2018-03-19 | 2018-10-02 | 西北大学 | Fingerprint base based on compressed sensing is established and RFID localization methods |
CN108693518A (en) * | 2018-05-15 | 2018-10-23 | 西北大学 | A kind of indoor orientation method |
CN108966369A (en) * | 2018-07-19 | 2018-12-07 | 广州泽祺信息科技有限公司 | A kind of home for destitute personnel positioning monitor system and its method |
CN112533134A (en) * | 2020-11-06 | 2021-03-19 | 浙江工业大学 | Wireless sensor network safety positioning method based on double detection |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103428850A (en) * | 2013-08-05 | 2013-12-04 | 湖南大学 | Compressed sensing based distributed multi-zone positioning method |
CN104898089A (en) * | 2015-04-03 | 2015-09-09 | 西北大学 | Device-free localization method based on space migration compressive sensing |
CN105005024A (en) * | 2015-07-21 | 2015-10-28 | 深圳市西博泰科电子有限公司 | Compression sensing and positioning method |
CN105682039A (en) * | 2016-03-15 | 2016-06-15 | 广州杰赛科技股份有限公司 | RF fingerprint positioning method and system |
-
2016
- 2016-11-23 CN CN201611047225.8A patent/CN106454750B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103428850A (en) * | 2013-08-05 | 2013-12-04 | 湖南大学 | Compressed sensing based distributed multi-zone positioning method |
CN104898089A (en) * | 2015-04-03 | 2015-09-09 | 西北大学 | Device-free localization method based on space migration compressive sensing |
CN105005024A (en) * | 2015-07-21 | 2015-10-28 | 深圳市西博泰科电子有限公司 | Compression sensing and positioning method |
CN105682039A (en) * | 2016-03-15 | 2016-06-15 | 广州杰赛科技股份有限公司 | RF fingerprint positioning method and system |
Non-Patent Citations (2)
Title |
---|
JUAN LUO,ZANYI HE: "Energy Confirmable Overlapping Target Tracking Based on Compressive Sensing in Wireless Sensor Networks", 《ADHOC & SENSOR WIRELESS NETWORKS》 * |
查俊莉: "基于压缩感知的RSS室内跟踪定位算法研究", 《万方数据知识服务平台》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107046680A (en) * | 2017-05-22 | 2017-08-15 | 湖南大学 | An Indoor Positioning Method Based on Trusted Location Fingerprint |
CN107046680B (en) * | 2017-05-22 | 2019-09-24 | 湖南大学 | A kind of indoor orientation method based on trustworthy location fingerprint |
CN108519577A (en) * | 2018-03-12 | 2018-09-11 | 中国矿业大学(北京) | A Distributed Localization Method Based on Compressed Sensing TOA Feature Signal Fingerprint Database |
CN108519577B (en) * | 2018-03-12 | 2023-09-15 | 中国矿业大学(北京) | Distributed positioning method based on compressed sensing TOA characteristic signal fingerprint library |
CN108614237A (en) * | 2018-03-19 | 2018-10-02 | 西北大学 | Fingerprint base based on compressed sensing is established and RFID localization methods |
CN108693518A (en) * | 2018-05-15 | 2018-10-23 | 西北大学 | A kind of indoor orientation method |
CN108693518B (en) * | 2018-05-15 | 2020-11-27 | 西北大学 | An indoor positioning method |
CN108966369A (en) * | 2018-07-19 | 2018-12-07 | 广州泽祺信息科技有限公司 | A kind of home for destitute personnel positioning monitor system and its method |
CN108966369B (en) * | 2018-07-19 | 2021-08-17 | 广州华创物联科技股份有限公司 | System and method for positioning and monitoring personnel in nursing home |
CN112533134A (en) * | 2020-11-06 | 2021-03-19 | 浙江工业大学 | Wireless sensor network safety positioning method based on double detection |
CN112533134B (en) * | 2020-11-06 | 2022-06-17 | 浙江工业大学 | A security location method for wireless sensor network based on double detection |
Also Published As
Publication number | Publication date |
---|---|
CN106454750B (en) | 2019-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106454750A (en) | Multi-region indoor safety positioning method based on compressed sensing technology | |
CN101676744B (en) | A high-precision tracking method for weak and small targets under complex background and low signal-to-noise ratio | |
CN103428850B (en) | Distributed multizone localization method based on compressed sensing | |
CN108601026A (en) | Perception data error attack detection method based on random sampling consistency | |
CN107046680B (en) | A kind of indoor orientation method based on trustworthy location fingerprint | |
CN109165671A (en) | Confrontation sample testing method based on sample to decision boundary distance | |
CN106411433A (en) | WLAN-based fine-grained indoor passive intrusion detection method | |
Meng et al. | Uav sensor spoofing detection algorithm based on gps and optical flow fusion | |
CN104581943A (en) | Node Localization Method for Distributed Wireless Sensor Networks | |
CN109688154B (en) | A method for establishing a network intrusion detection model and a network intrusion detection method | |
CN104732543A (en) | Infrared weak small target fast detecting method under desert and gobi background | |
Bi et al. | Detection and mitigation of position spoofing attacks on cooperative uav swarm formations | |
CN108834043B (en) | A priori knowledge-based compressed sensing multi-target passive localization method | |
Tian et al. | Indoor device-free passive localization for intrusion detection using multi-feature PNN | |
CN118566904A (en) | Airport periphery intrusion prevention method based on cooperative work of radar detection | |
Pala et al. | Advanced Deep Learning Solutions for Automated Diagnosis of Solar Panel Issues | |
CN107271991B (en) | A kind of optical electrical sensor target correlating method based on state estimation | |
Ning et al. | A Dual-Path Neural Network for High-Impedance Fault Detection | |
CN112287564B (en) | Electrode array optimization method based on goblet sea squirt group algorithm | |
Li et al. | Data-driven false data injection attacks on state estimation in smart grid | |
Li et al. | [Retracted] Abnormal Data Detection in Sensor Networks Based on DNN Algorithm and Cluster Analysis | |
CN109819397B (en) | An Approximate Triangle Interior Point Test Location Algorithm Against Sybil Attacks | |
Li et al. | Full-view coverage holes detection and healing solutions | |
Su et al. | Multi-sensor data fusion in wireless sensor network for target detection | |
Khanderiya et al. | A Novel Approach for Detection of Sybil Attack in Wireless Sensor Networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |