CN106597368B - 室内干扰源的定位方法及系统 - Google Patents

室内干扰源的定位方法及系统 Download PDF

Info

Publication number
CN106597368B
CN106597368B CN201611146398.5A CN201611146398A CN106597368B CN 106597368 B CN106597368 B CN 106597368B CN 201611146398 A CN201611146398 A CN 201611146398A CN 106597368 B CN106597368 B CN 106597368B
Authority
CN
China
Prior art keywords
signal
interference
arrival
angle
interference signal
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
CN201611146398.5A
Other languages
English (en)
Other versions
CN106597368A (zh
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.)
Shenzhen University
Original Assignee
Shenzhen University
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 Shenzhen University filed Critical Shenzhen University
Priority to CN201611146398.5A priority Critical patent/CN106597368B/zh
Publication of CN106597368A publication Critical patent/CN106597368A/zh
Application granted granted Critical
Publication of CN106597368B publication Critical patent/CN106597368B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/04Position of source determined by a plurality of spaced direction-finders
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • H04W64/003Locating users or terminals or network equipment for network management purposes, e.g. mobility management locating network equipment

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)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Noise Elimination (AREA)

Abstract

本发明提供了一种室内干扰源的定位方法及系统,该定位方法包括:S1.根据无线接收端采集的无线信号,采用伪随机序列将工作信号与干扰信号分离;S2.利用功率延迟分布特性,将干扰信号中的直视路径分量分离;S3.确定干扰信号直视路径的到达角度;S4.接收端线性平移一段距离后,在新的位置上重复步骤S1、S2和S3,计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。本发明的有益效果是:本发明不需要预先部署多个接收端,只需要一台商用多天线接收端。通过移动可以消除多径的影响,提高了定位的准确性。在多个不同位置测量,等效多台设备协同工作进行定位,进一步提高了定位的准确性和效率。

Description

室内干扰源的定位方法及系统
技术领域
本发明涉及无线通信技术领域,尤其涉及室内干扰源的定位方法及系统。
背景技术
随着无线通信技术的不断发展,室内环境下的终端设备越来越多,因此无线接入点的数量也随着增长。如何为每个用户提供高质量的无线通信服务,是网络提供商关注的重点。由于无线频谱资源有限,现有的大量Wi-Fi网络设备产生的相同频率的信号之间会产生严重的干扰,降低了通信的质量。在室内环境下,同种类型的无线设备之间相同频率的干扰也难以消除,使有用信号和干扰信号混叠,降低网络的性能。为了消除同频设备之间的干扰,往往需要确定对设备造成干扰的干扰源的位置,然后进行干扰的消除。为了实现对室内设备的定位,人们提出了基于接收信号强度指示(Received-Signal-Strength-Indicator,RSSI)、到达角度(Angle-of-Arrival,AoA)等方法,获取预先部署的多个接收端接收到的信号相关属性值,合并这些信息确定设备的位置。然而,在室内环境下,多径的影响使采用这些方法的系统都存在种种不足。基于RSSI的系统利用距离与信号强度的关系,算出每个接收端与干扰源的距离,最后用三边测量的方法确定设备的位置,室内多径会使RSSI变得不稳定,影响定位的精度。而基于AoA的方法综合多个接收端的AoA信息和RSSI信息来定位目标,但多径同样会影响AoA的测量准确性,引入定位的误差。现有定位方法还有人工定位的方式,测试人员利用扫频仪和方向性天线在受到干扰的设备附近扫描,确定干扰源的方向。再根据干扰信号的强度估计干扰源的距离,从而确定干扰源的位置。该方法具有较高的准确性,但测试人员需要携带专门的检测设备,会造成不便,同时人工定位方式效率低。
上述的定位方法都有准确性和效率较低的问题,而且需要部署较多的设备共同完成定位。
发明内容
本发明提供了一种室内干扰源的定位方法,包括如下步骤:
S1.根据无线接收端采集的无线信号,采用伪随机序列将工作信号与干扰信号分离;
S2.利用功率延迟分布特性,将干扰信号中的直视路径分量分离;
S3.确定干扰信号直视路径的到达角度;
S4.接收端线性平移一段距离后,在新的位置上重复步骤S1、S2和S3,计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
作为本发明的进一步改进,所述步骤S1分离干扰信号包括:
S11.接收端采集收到的无线信号,根据设备的工作频率,使用带通滤波器过滤噪声及带外干扰信号,得到待处理信号;
S12.利用已知伪随机序列与所述待处理信号模二加,得到工作信号,再从接收端接收的混合信号中减去工作信号,得到干扰信号。
作为本发明的进一步改进,所述步骤S2分离干扰信号直视路径分量包括:
S21.利用多径信号在时域上的色散,功率延迟分布特性被用来分离通过不同路径达到接收天线的多径成分,设定一个功率阈值,将上述干扰信号各路径信号大于该阈值的认为是主路径信息;
S22.将所述主路径信息通过傅里叶变换从时域转到频域。
作为本发明的进一步改进,在所述步骤S3中,运用MUSIC算法确定干扰信号直视路径的到达角度。
作为本发明的进一步改进,在所述步骤S3中,根据信号达到天线阵列的相位偏移,计算信号到达不同天线的时间差,利用信号在天线阵列到达时间的不同,运用MUSIC算法将干扰信号直视路径的到达角度θ求出;在所述步骤S4中,直线平移距离d到达新的位置,重复步骤S1,S2和S3,得到新位置的干扰信号的到达角,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
本发明还提供了一种室内干扰源的定位系统,包括信号获取模块、信号分析模块和定位模块,所述信号分析模块包括分离模块和计算模块,
信号获取模块:根据无线接收端采集的无线信号,采用伪随机序列将工作信号与干扰信号分离;
分离模块:利用功率延迟分布特性,将干扰信号中的直视路径分量分离;
计算模块:确定干扰信号直视路径的到达角度;
定位模块:接收端线性平移一段距离后,在新的位置上重复信号获取模块、分离模块和计算模块,计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
作为本发明的进一步改进,所述信号获取模块包括:
采集单元:接收端采集收到的无线信号,根据设备的工作频率,使用带通滤波器过滤噪声及带外干扰信号,得到待处理信号;
数据处理单元:利用已知伪随机序列与所述待处理信号模二加,得到工作信号,再从接收端接收的混合信号中减去工作信号,得到干扰信号。
作为本发明的进一步改进,所述分离模块包括:
提取主路径信息单元:利用多径信号在时域上的色散,功率延迟分布特性被用来分离通过不同路径达到接收天线的多径成分,设定一个功率阈值,将上述干扰信号各路径信号大于该阈值的认为是主路径信息;
转换单元:将主路径信息通过傅里叶变换从时域转到频域。
作为本发明的进一步改进,在所述计算模块中,运用MUSIC算法确定干扰信号直视路径的到达角度。
作为本发明的进一步改进,在所述计算模块中,根据信号达到天线阵列的相位偏移,计算信号到达不同天线的时间差,利用信号在天线阵列到达时间的不同,运用MUSIC算法将干扰信号直视路径的到达角度θ求出;
在所述定位模块中,直线平移距离d到达新的位置,重复信号获取模块、分离模块和计算模块,得到新位置的干扰信号的到达角,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
本发明的有益效果是:本发明不需要预先部署多个接收端,只需要一台商用多天线接收端。通过移动可以消除多径的影响,提高了定位的准确性。在多个不同位置测量,等效多台设备协同工作进行定位,进一步提高了定位的准确性和效率。
附图说明
图1是本发明的定位方法一实施例的流程图;
图2是本发明的干扰源与接收器之间的位置示意图。
具体实施方式
本发明公开了一种室内干扰源的定位方法,包括如下步骤:
S1.根据无线接收端采集的无线信号,采用伪随机序列(Pseudo-RandomSequence,PN序列)将工作信号与干扰信号分离;
S2.利用功率延迟分布(Power-Delay-Profile,PDP)特性,将干扰信号中的直视路径分量分离;
S3.运用MUSIC算法确定干扰信号直视路径的到达角度(Angle-of-Arrival,AOA);
S4.接收端线性平移一段距离后,在新的位置上重复步骤S1,S2和S3,计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
本发明的接收端数目为1,接收端的天线数为2个或2个以上。所采用的设备为现有无线商用设备,不需要安装专用设备。
在实际应用中,我们使用多天线无线网卡来接收无线信号,干扰源为无线路由器AP。该方法是基于AP的移动,在不同的位置上测量干扰信号直视路径的到达角度,由到达角度确定的路径轨迹的交点,作为干扰源的位置。工作于同一信道的AP之间会产生较大的干扰,带通滤波器并不能滤除同频干扰。为了定位该干扰AP,需要将干扰信号直视路径分量分离,利用MUSIC算法得到波束的到达角度。通过移动,多次测量上述干扰的到达角度。三角测量法被用来确定干扰源的位置。
所述步骤S1分离干扰信号包括:
S11.接收端采集收到的无线信号,根据设备的工作频率,使用带通滤波器过滤噪声及带外干扰信号,得到待处理信号;
S12.利用已知伪随机序列与所述待处理信号模二加,得到工作信号,再从接收端接收的混合信号中减去工作信号,得到干扰信号。
当本系统开始工作时,无线接收端采集附近区域的接收到的无线信号,以802.11nWiFi的第1信道为例,该信道的中心频率为2.412GHz,信道宽度为20MHz。可选取中心频率为2.412GHz,通频带为20MHz的带通滤波器将带外的信号和干扰去除,得到待处理的混合信号。当物理层使用直接序列扩频时,利用已知的PN序列可以将有用信号解调出来,在混合信号中将有用信号去除,得到干扰信号。
所述步骤S2分离干扰信号直视路径分量包括:
S21.利用多径信号在时域上的色散,功率延迟分布特性被用来分离通过不同路径达到接收天线的多径成分,设定一个功率阈值,将上述干扰信号各路径信号大于该阈值的认为是主路径信息;
S22.将所述主路径信息通过傅里叶变换(Fast-Fourier-Transform,FFT)从时域转到频域。
室内环境下,因为墙壁、地板等反射物的存在,信号从发射端到接收端的过程会经历许多的路径,不同的路径在达到时间、到达角度上都会有所不同。信号在不同的路径经历的衰减也是不同的,他们在时域上会发生色散。较短路径上的信号在到达天线阵列的时间较早,而较长的路径的信号到达天线阵列的时间较晚。所以功率延迟分布可以被用来分离直视路径成分和非直视路径成分。通过预先设定的一个功率阈值,将大于该阈值的信号认为是可能的直视路径成分。将选取延时较短部分的信号作为直视路径分量,通过FFT将信号从时域转换到频率,作为下一步骤的输入。
在所述步骤S3中,根据信号达到天线阵列的相位偏移,计算信号到达不同天线的时间差,利用信号在天线阵列到达时间的不同,运用MUSIC算法将干扰信号直视路径的到达角度θ求出;
在阵列信号模型中,假设接收端一共有M根天线,天线间距为d,涉及空间存在I个信号。多径的数量为ki条。假设信号源是远场源,所以可将接收信号看作平面波。在第m根天线接收到的信号可以表示为
其中,代表第i个信号的第k个成分,θ1k表示对应的到达角。通过调用标准MUSIC算法,可以得到MUSIC谱,该谱中的尖峰可得到相应的到达角度θ。
在所述步骤S4中,直线平移距离d到达新的位置,重复步骤S1,S2和S3,得到新位置的干扰信号的到达角,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
如图2所示,干扰源的位置由接收端在不同位置的到达角度确定。
具体地,如图1所示,实现室内干扰源定位的流程包括:
1)无线发射端(AP)发送无线信号,无线接收端(带无线网卡的电
脑)接收无线信号,收集数据;
2)带通滤波器过滤噪声及带外干扰信号;
3)利用已知伪随机序列与上述信号模二加,得到有用信号;
4)从接收端接收的混合信号中减去有用信号,得到干扰信号;
5)通过功率延迟分布获得主路径信息;
6)通过MUSIC算法得到干扰信号到达角;
7)多次移动接收端得到多个到达角,用三角测量法确定干扰源的位置。
本发明不需要预先部署多个接收端,只需要一台商用多天线接收端。通过移动可以消除多径的影响,提高了定位的准确性。在多个不同位置测量,等效多台设备协同工作进行定位,进一步提高了定位的准确性和效率。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (10)

1.一种室内干扰源的定位方法,其特征在于,包括如下步骤:
S1.根据无线接收端采集的无线信号,采用伪随机序列将工作信号与干扰信号分离;
S2.利用功率延迟分布特性,将干扰信号中的直视路径分量分离;
S3.确定干扰信号直视路径的到达角度;
S4.接收端线性平移一段距离后,在新的位置上重复步骤S1、S2和S3, 计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
2.根据权利要求1所述的定位方法,其特征在于,所述步骤S1分离干扰信号包括:
S11.接收端采集收到的无线信号,根据设备的工作频率,使用带通滤波器过滤噪声及带外干扰信号,得到待处理信号;
S12.利用已知伪随机序列与所述待处理信号模二加,得到工作信号,再从接收端接收的混合信号中减去工作信号,得到干扰信号。
3.根据权利要求1所述的定位方法,其特征在于,所述步骤S2分离干扰信号直视路径分量包括:
S21.利用多径信号在时域上的色散,功率延迟分布特性被用来分离通过不同路径达到接收天线的多径成分,设定一个功率阈值,将上述干扰信号各路径信号大于该阈值的认为是主路径信息;
S22.将所述主路径信息通过傅里叶变换从时域转到频域。
4.根据权利要求1至3任一项所述的定位方法,其特征在于,在所述步骤S3中,运用MUSIC算法确定干扰信号直视路径的到达角度。
5.根据权利要求4所述的定位方法,其特征在于,
在所述步骤S3中,根据信号达到天线阵列的相位偏移,计算信号到达不同天线的时间差,利用信号在天线阵列到达时间的不同,运用MUSIC算法将干扰信号直视路径的到达角度θ求出。
6.一种室内干扰源的定位系统,其特征在于,包括信号获取模块、信号分析模块和定位模块,所述信号分析模块包括分离模块和计算模块,
信号获取模块:根据无线接收端采集的无线信号,采用伪随机序列将工作信号与干扰信号分离;
分离模块:利用功率延迟分布特性,将干扰信号中的直视路径分量分离;
计算模块:确定干扰信号直视路径的到达角度;
定位模块:接收端线性平移一段距离后,在新的位置上重复信号获取模块、分离模块和计算模块, 计算上述干扰信号的到达角度,综合多个位置的到达角度,用三角测量法确定干扰源的位置。
7.根据权利要求6所述的定位系统,其特征在于,所述信号获取模块包括:
采集单元:接收端采集收到的无线信号,根据设备的工作频率,使用带通滤波器过滤噪声及带外干扰信号,得到待处理信号;
数据处理单元:利用已知伪随机序列与所述待处理信号模二加,得到工作信号,再从接收端接收的混合信号中减去工作信号,得到干扰信号。
8.根据权利要求6所述的定位系统,其特征在于,所述分离模块包括:
提取主路径信息单元:利用多径信号在时域上的色散,功率延迟分布特性被用来分离通过不同路径达到接收天线的多径成分,设定一个功率阈值,将上述干扰信号各路径信号大于该阈值的认为是主路径信息;
转换单元:将主路径信息通过傅里叶变换从时域转到频域。
9.根据权利要求6至8任一项所述的定位系统,其特征在于,在所述计算模块中,运用MUSIC算法确定干扰信号直视路径的到达角度。
10.根据权利要求9所述的定位系统,其特征在于,在所述计算模块中,根据信号达到天线阵列的相位偏移,计算信号到达不同天线的时间差,利用信号在天线阵列到达时间的不同,运用MUSIC算法将干扰信号直视路径的到达角度θ求出。
CN201611146398.5A 2016-12-13 2016-12-13 室内干扰源的定位方法及系统 Active CN106597368B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611146398.5A CN106597368B (zh) 2016-12-13 2016-12-13 室内干扰源的定位方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611146398.5A CN106597368B (zh) 2016-12-13 2016-12-13 室内干扰源的定位方法及系统

Publications (2)

Publication Number Publication Date
CN106597368A CN106597368A (zh) 2017-04-26
CN106597368B true CN106597368B (zh) 2018-07-24

Family

ID=58801091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611146398.5A Active CN106597368B (zh) 2016-12-13 2016-12-13 室内干扰源的定位方法及系统

Country Status (1)

Country Link
CN (1) CN106597368B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102018215675B4 (de) * 2018-09-14 2022-10-06 Vega Grieshaber Kg Fremdstrahlungserkennung mit Gamma-Modulator
CN111278040A (zh) * 2018-12-05 2020-06-12 中国移动通信集团四川有限公司 干扰源定位方法、装置、设备及计算机存储介质
CN111601332A (zh) * 2019-02-21 2020-08-28 大唐移动通信设备有限公司 一种干扰源位置的确定方法及装置
CN110646807A (zh) * 2019-09-19 2020-01-03 上海兰宝传感科技股份有限公司 一种基于激光三角测量系统抗环境光干扰的处理方法
CN110855385A (zh) * 2019-11-27 2020-02-28 深圳格林帕科技有限公司 针对物联网移动设备的动态无线环境的硬仿真设备
CN113009414B (zh) * 2019-12-20 2024-03-19 中移(成都)信息通信科技有限公司 信号源位置确定方法、装置、电子设备及计算机存储介质
CN112601249B (zh) * 2020-12-08 2023-02-28 重庆紫光华山智安科技有限公司 干扰源识别方法、装置、网络设备及存储介质
CN114817862B (zh) * 2022-05-26 2023-09-26 江苏省地震局 定位电磁干扰源的方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008151533A (ja) * 2006-12-14 2008-07-03 Samsung Yokohama Research Institute Co Ltd 側位システムおよび無線通信装置
CN104181503A (zh) * 2014-08-04 2014-12-03 武汉邮电科学研究院 一种基于rss的led定位系统及方法
CN105247934A (zh) * 2013-04-12 2016-01-13 惠普发展公司,有限责任合伙企业 移动设备的位置确定
CN106125046A (zh) * 2016-08-20 2016-11-16 华南理工大学 一种基于cdma调制的rss/aoa三维定位方法
CN106130602A (zh) * 2016-08-20 2016-11-16 华南理工大学 一种基于可见光室内定位的cdma编码方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140329536A1 (en) * 2013-05-01 2014-11-06 Qualcomm Incorporated Synthetic wideband ranging design

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008151533A (ja) * 2006-12-14 2008-07-03 Samsung Yokohama Research Institute Co Ltd 側位システムおよび無線通信装置
CN105247934A (zh) * 2013-04-12 2016-01-13 惠普发展公司,有限责任合伙企业 移动设备的位置确定
CN104181503A (zh) * 2014-08-04 2014-12-03 武汉邮电科学研究院 一种基于rss的led定位系统及方法
CN106125046A (zh) * 2016-08-20 2016-11-16 华南理工大学 一种基于cdma调制的rss/aoa三维定位方法
CN106130602A (zh) * 2016-08-20 2016-11-16 华南理工大学 一种基于可见光室内定位的cdma编码方法

Also Published As

Publication number Publication date
CN106597368A (zh) 2017-04-26

Similar Documents

Publication Publication Date Title
CN106597368B (zh) 室内干扰源的定位方法及系统
CN104812061B (zh) 一种基于mimo-ofdm信道状态信息的室内测距及定位方法
CN106909779B (zh) 基于分布式处理的mimo雷达克拉美罗界计算方法
Amar et al. Direct position determination (DPD) of multiple known and unknown radio-frequency signals
CN103713277B (zh) 一种基于位置信息场的辐射源定位方法
CN106792808A (zh) 一种基于信道状态信息的室内环境下视距路径识别方法
CN104076349B (zh) 一种基于多普勒频移的被动式移动目标定位方法
CN102364885B (zh) 一种基于信号频谱包络的频谱感知方法
CN107132505A (zh) 直达与非直达混合场景中的多目标直接定位方法
CN105954712A (zh) 联合无线电信号复包络和载波相位信息的多目标直接定位方法
CN107843910B (zh) 一种适于复杂环境下的虚拟多站tdoa定位方法与装置
CN114415210B (zh) 一种多维信息联合的卫星导航欺骗式干扰源直接定位方法
Banin et al. Next generation indoor positioning system based on WiFi time of flight
CN105182322A (zh) 基于反射信号相位差的被动式定位方法
CN103945526A (zh) 基于诱发探测技术的无线设备定位方法及系统
CN109613480A (zh) 一种可消除多径无噪声信道下符号间干扰的超宽带定位系统
KR101835975B1 (ko) 무선 통신 시스템에서 단말의 위치를 측정하기 위한 방법, 로케이션 서버 및 단말
CN111381226B (zh) 一种基于多频段接收信号强度的测距增强方法
CN205880195U (zh) 无线电测向系统
CN111770527A (zh) 基于二维特征的可视与非可视信道识别方法及装置
Astafiev et al. Algorithm for preliminary processing channel state information of the WIFI communication channel for building indoor positioning systems
CN113203985B (zh) 一种短波同频信号直接定位方法
Pereira et al. Multi-technology RF fingerprinting with leaky-feeder in underground tunnels
Ghoraishi et al. Azimuth and delay dispersion of mobile radio wave propagation through vegetation
CN114401055A (zh) 一种智能频谱检测系统

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