CN114513252A - 基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法 - Google Patents

基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法 Download PDF

Info

Publication number
CN114513252A
CN114513252A CN202210041295.1A CN202210041295A CN114513252A CN 114513252 A CN114513252 A CN 114513252A CN 202210041295 A CN202210041295 A CN 202210041295A CN 114513252 A CN114513252 A CN 114513252A
Authority
CN
China
Prior art keywords
point
density
points
density peak
clustering algorithm
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
Application number
CN202210041295.1A
Other languages
English (en)
Other versions
CN114513252B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202210041295.1A priority Critical patent/CN114513252B/zh
Publication of CN114513252A publication Critical patent/CN114513252A/zh
Application granted granted Critical
Publication of CN114513252B publication Critical patent/CN114513252B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/075Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
    • H04B10/079Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
    • H04B10/0795Performance monitoring; Measurement of transmission parameters
    • H04B10/07953Monitoring or measuring OSNR, BER or Q
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50Transmitters
    • H04B10/516Details of coding or modulation
    • H04B10/54Intensity modulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • H04B10/6161Compensation of chromatic dispersion
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • H04B10/6165Estimation of the phase of the received optical signal, phase error estimation or phase error correction
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Optical Communication System (AREA)

Abstract

本发明提供了一种用于概率整形相干光通信系统的非线性补偿方法,利用加权的密度峰值聚类算法,对接收到的概率整形信号进行处理,包括:计算接收数据中每个对象的密度和最小相对距离,并根据概率整形分布因子给每个对象赋予权值,将三项相乘后得到的参量按从大到小排列选出密度峰值点。然后根据密度峰值点给剩下的对象分配标签。本发明可以对不均匀的概率整形信号进行聚类,有效提升系统的性能,使系统误码率得到明显改善。

Description

基于加权的密度峰值聚类算法的概率整形相干光通信系统非 线性补偿方法
技术领域
本发明涉及光纤通信领域,具体涉及一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法。
背景技术
在光通信系统中,非线性效应已经是目前通信系统最主要的受限因素,对于非线性损伤的补偿有数字反向传播(DBP),光学相位共轭等方法,但是这些方法要么复杂度过高,要么需要在链路中指定位置安装共轭装置导致成本过高,或者牺牲了频谱效率。近年来机器学习算法开始应用到光纤通信系统中,如K近邻、支持向量机(SVM)、基于密度的噪声应用空间聚类算法(DBSCAN)和K均值聚类算法(K-means)等。然而这些机器学习方法都是针对标准QAM信号的,在概率整形系统中,由于分布的不均匀带来了一些新的问题,不同的簇的数据量相差很大,在对其进行聚类的时候有可能会将同一个簇的数据错误的分类成不同的簇或者将不同的簇的数据分到同一簇中。针对这种情况,提出了一种加权的密度峰值聚类算法(Weighted-DPC),根据概率整形信号中不同幅度的点出现的概率不同这一特点设置一个权重,可以更准确的找到密度峰值,从而准确的将数据分到正确的簇中。
概率整形信号的星座图是不均匀的,不同类别的数据相差很大。当不同类别的数据数量相差太大时,DPC算法可能无法正确的找到密度峰值点。在概率整形信号的星座图中,靠近外部幅度更大的簇数据点更少,内部的簇数据点则更密集,在密度峰值点选取的时候,不像均匀分布的信号那样可以直接通过参量γ值排序后直接选出来。因此在计算γ时加入权值,更准确的对数据进行聚类。
发明内容
本发明提供了一种基于加权密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,根据概率整形信号分布的特点加入权值,对接收数据准确聚类从而有效弥补信号的非线性损伤,提高系统的误码率性能。
为了达到上述目的,总体思路是:提出一种加权的密度峰值聚类算法,由于恒成分分布匹配器方案生成的概率整形信号遵循Maxwell-Boltzmann分布,幅度越大的点出现的概率越小,在密度峰值聚类算法中,确定密度峰值点时加入权值,幅度更大的点赋予更高的权重,更准确的找到信号星座外围数据点分布稀疏的簇的密度峰值点。
具体的,本发明采用的是一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,包括以下步骤:
计算接收到的PS-MQAM信号数据集中每个点的密度ρi,其中PS-MQAM信号为概率整形正交幅度调制信号;
计算数据集每个点的相对最小距离δi,δi表示的是,在数据集中所有比点i密度大的点之中,距离点i最近的点与点i的距离;
引入加权因子αi,并将其与密度ρi以及相对最小距离δi相乘得到参量γi,然后将参量γi从大到小排序;
对于所述PS-MQAM信号数据集,取前M个γi对应的数据点作为密度峰值点,将选出的M个密度峰值点作为聚类中心,其中M为所述PS-MQAM信号数据集中簇的个数;
对剩余点按密度从大到小的顺序依次进行分配,分配时,将每个剩余点分配到与它最近邻且密度比它大的数据点所在的簇。
优选的,所述δi通过以下公式计算:
Figure BDA0003470356450000021
其中j为密度比i大的点,dij为点i与j之间的欧式距离。
优选的,对于全局密度最大的点,δi是数据集所有点中离该点的距离的最大值,此时所述δi通过以下公式计算:
Figure BDA0003470356450000022
优选的,
Figure BDA0003470356450000023
其中dc为截断距离。
优选的,所述加权因子αi由下式获取:
Figure BDA0003470356450000031
其中ν为恒成分分布匹配器概率整形方案中的Maxwell-Boltzmann分布整形因子,其中xi是接收的数据,xk是所述PS-MQAM信号数据集的M个标准星座点。
优选的,参量γi由下式获取:
γi=δi×ρi×αi
优选的,所述截断距离dc取值为0.05。
优选的,所述ρi使用高斯核方式计算。
相比于现有技术,本发明的有益效果在于:传统的聚类或分类算法如K-means、KNN等算法会受到数据分布不均匀的影响,本发明根据概率整形信号分布的特点,基于密度峰值聚类算法,在寻找密度峰值点时,对幅度更大出现概率更小的点赋予更大的权重,从而准确找到概率整形信号星座点外围相对稀疏的簇的密度峰值点,从而准确聚类,最终实现了信号的非线性补偿。
附图说明
图1是本发明实施例的装置示意图;
图2是密度峰值聚类算法的流程图;
图3是经典的密度峰值聚类算法和加权的密度峰值聚类算法在概率整形信号星座中聚类的结果对比;
图4是整形因子ν=0.2时PS-16QAM信号传输1600KM后的结果图;
图5是整形因子ν=0.1时PS-32QAM信号传输960KM后的结果图;
图6是整形因子ν=0.05时PS-64QAM信号传输320KM后的结果图;
具体实施方式
为了更详细的说明实施例和本发明的具体实现及有点,下面结合附图作进一步描述:
一种基于加权的密度峰值聚类的概率整形相干光通信系统非线性补偿方法,使用的装置如附图1所示,在发射端将比特流经马赫曾德调制器调制后进入光纤,经过长距离传输后由接收机接收,将光信号经光电转换后成为电信号,在电域进行色散补偿、偏振解复用和载波相位恢复后,使用本发明对信号的非线性损伤进行补偿。
传统的密度峰值聚类算法属于无监督算法,其准确聚类的关键在于找到每个簇的密度峰值点。作为密度峰值点,首先一定具有较大的密度值ρi,其次距离其他密度比自身大的点较远,因此有较大的最小相对距离δi。计算出综合考虑的参量γi,有M个簇就选出M个γi最大点作为密度峰值点。而在概率整形信号这种分布不均匀的数据中,数据分布稀疏的簇由于密度相对较低,计算出的参量γi较小,无法准确的找到密度峰值点。因此对分布稀疏的簇赋予更大的权重,使其计算出的参量γi增大,从而能准确的找到密度峰值点。
根据以上思路,本发明实施例提出一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,包括以下步骤:
计算接收到的PS-MQAM信号数据集中每个点的密度ρi,其中PS-MQAM信号为概率整形正交幅度调制信号;
计算数据集每个点的相对最小距离δi,δi表示的是,在数据集中所有比点i密度大的点之中,距离点i最近的点与点i的距离;
引入加权因子αi,并将其与密度ρi以及相对最小距离δi相乘得到参量γi,然后将参量γi从大到小排序,对于所述PS-MQAM信号数据集,取前M个γi对应的数据点作为密度峰值点,将选出的M个密度峰值点作为聚类中心,其中M为所述PS-MQAM信号数据集中簇的个数,然后对剩余点按密度从大到小的顺序依次进行分配,分配时,将每个剩余点分配到与它最近邻且密度比它大的数据点所在的簇。
具体的,所述δi通过以下公式计算:
Figure BDA0003470356450000041
具体的,对于全局密度最大的点,δi是数据集所有点中离该点的距离的最大值,此时所述δi通过以下公式计算:
Figure BDA0003470356450000042
具体的,
Figure BDA0003470356450000051
其中dij为两个数据点之间的欧式距离,dc为截断距离。
优选的,所述加权因子αi由下式获取:
Figure BDA0003470356450000052
其中ν为恒成分分布匹配器概率整形方案中的Maxwell-Boltzmann分布整形因子,其中xi是接收的数据,xk是所述PS-MQAM信号数据集的M个标准星座点。
具体的,参量γi由下式获取:
γi=δi×ρi×αi
具体的,所述截断距离dc取值为0.05。
具体的,所述ρi使用高斯核方式计算。
图3a是根据没有加权的参量γi画出的决策图及选出的峰值点,图3b是加权后的决策图和峰值点。图3c和图3e是未加权的聚类结果,可以看到由于分布不均匀导致密度峰值点选取错误,最终导致了错误的聚类结果。图3d和图3f是加权后的聚类结果,由于提高了星座图中外围幅度大但是分布稀疏的点的权重,找到了正确的密度峰值点。
图4是整形因子ν为0.2时PS-16QAM信号传输1600KM后的结果图,图5是ν=0.1时PS-32QAM信号传输960KM后的结果图,图6是ν=0.05时PS-64QAM信号传输320KM后的结果图。由图可以看到,入纤光功率较低时,系统性能主要受到信噪比和自发辐射噪声主导,误码率随功率的增大而降低,此时使用加权的密度峰值聚类算法后系统性能无明显改善;继续增大入纤功率,系统性能开始降低,此时信号损伤主要来自非线性效应,使用加权密度聚类算法后,系统性能得到明显的提高,说明此方法可以有效补偿概率整形信号的非线性损伤。

Claims (8)

1.一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于,包括以下步骤:
计算接收到的PS-MQAM信号数据集中每个点的密度ρi,其中PS-MQAM信号为概率整形正交幅度调制信号;
计算数据集每个点的相对最小距离δi,δi表示的是,在数据集中所有比点i密度大的点之中,距离点i最近的点与点i的距离;
引入加权因子αi,并将其与密度ρi以及相对最小距离δi相乘得到参量γi,然后将参量γi从大到小排序;
对于所述PS-MQAM信号数据集,取前M个γi对应的数据点作为密度峰值点,将选出的M个密度峰值点作为聚类中心,其中M为所述PS-MQAM信号数据集中簇的个数;
对剩余点按密度从大到小的顺序依次进行分配,分配时,将每个剩余点分配到与它最近邻且密度比它大的数据点所在的簇。
2.根据权利要求1所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于,所述δi通过以下公式计算:
Figure FDA0003470356440000011
其中j为密度比i大的点,dij为点i与j之间的欧式距离。
3.根据权利要求1所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于:对于全局密度最大的点,δi是数据集所有点中离该点的距离的最大值,此时所述δi通过以下公式计算:
Figure FDA0003470356440000012
4.根据权利要求1所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于:
Figure FDA0003470356440000013
其中dc为截断距离。
5.根据权利要求1所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于:所述加权因子αi由下式获取:
Figure FDA0003470356440000021
其中ν为恒成分分布匹配器概率整形方案中的Maxwell-Boltzmann分布整形因子,其中xi是接收的数据,xk是所述PS-MQAM信号数据集的M个标准星座点。
6.根据权利要求1所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于,参量γi由下式获取:
γi=δi×ρi×αi
7.根据权利要求2所述的所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于,所述截断距离dc取值为0.05。
8.根据权利要求1所述的所述的一种基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法,其特征在于,所述ρi使用高斯核方式计算。
CN202210041295.1A 2022-01-14 2022-01-14 基于加权密度峰值聚类算法的光纤非线性补偿方法 Active CN114513252B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210041295.1A CN114513252B (zh) 2022-01-14 2022-01-14 基于加权密度峰值聚类算法的光纤非线性补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210041295.1A CN114513252B (zh) 2022-01-14 2022-01-14 基于加权密度峰值聚类算法的光纤非线性补偿方法

Publications (2)

Publication Number Publication Date
CN114513252A true CN114513252A (zh) 2022-05-17
CN114513252B CN114513252B (zh) 2023-06-20

Family

ID=81549115

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210041295.1A Active CN114513252B (zh) 2022-01-14 2022-01-14 基于加权密度峰值聚类算法的光纤非线性补偿方法

Country Status (1)

Country Link
CN (1) CN114513252B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108667523A (zh) * 2018-03-06 2018-10-16 苏州大学 基于无数据辅助的knn算法的光纤非线性均衡方法
CN108965178A (zh) * 2018-06-29 2018-12-07 苏州大学 基于机器学习的智能自适应均衡器及均衡解调方法
CN109104248A (zh) * 2018-07-06 2018-12-28 西南交通大学 一种基于SOFM神经网络聚类的DRoF前端信号量化方法
CN110190906A (zh) * 2019-04-15 2019-08-30 浙江工业大学 一种针对高阶qam相干光系统的非线性均衡方法
US20190356527A1 (en) * 2018-05-15 2019-11-21 Cable Television Laboratories, Inc. Carrier-phase recovery system and method
CN111371503A (zh) * 2020-03-17 2020-07-03 华南师范大学 一种概率整形星座调制相干光通信系统中盲偏振解复用的方法
CN113346957A (zh) * 2021-06-01 2021-09-03 北京理工大学 一种oam-qpsk传输的聚类非线性补偿方法
CN113541817A (zh) * 2021-07-02 2021-10-22 华南师范大学 星座概率整形相干光通信系统的偏振相关损耗补偿方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108667523A (zh) * 2018-03-06 2018-10-16 苏州大学 基于无数据辅助的knn算法的光纤非线性均衡方法
US20190356527A1 (en) * 2018-05-15 2019-11-21 Cable Television Laboratories, Inc. Carrier-phase recovery system and method
CN108965178A (zh) * 2018-06-29 2018-12-07 苏州大学 基于机器学习的智能自适应均衡器及均衡解调方法
CN109104248A (zh) * 2018-07-06 2018-12-28 西南交通大学 一种基于SOFM神经网络聚类的DRoF前端信号量化方法
CN110190906A (zh) * 2019-04-15 2019-08-30 浙江工业大学 一种针对高阶qam相干光系统的非线性均衡方法
CN111371503A (zh) * 2020-03-17 2020-07-03 华南师范大学 一种概率整形星座调制相干光通信系统中盲偏振解复用的方法
CN113346957A (zh) * 2021-06-01 2021-09-03 北京理工大学 一种oam-qpsk传输的聚类非线性补偿方法
CN113541817A (zh) * 2021-07-02 2021-10-22 华南师范大学 星座概率整形相干光通信系统的偏振相关损耗补偿方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ELIAS GIACOUMIDIS等: "A blind nonlinearity compensator using DBSCAN clustering for coherent optical transmission systems" *
JUNFENG ZHANG等: "Blind and Noise-Tolerant Modulation Format Identification" *
董颖娣;彭进业;: "量子密钥分配系统中非线性偏振耦合的研究", no. 06 *
谢连妮;曹蕾;张璐;余玉揆;: "光纤通信中的克尔非线性补偿技术", no. 06 *

Also Published As

Publication number Publication date
CN114513252B (zh) 2023-06-20

Similar Documents

Publication Publication Date Title
CN108880692B (zh) 一种相干光通信系统调制格式识别和光信噪比监控方法
CN107707494B (zh) 用于64-qam相干光通信系统的光纤非线性均衡方法
Li et al. Nonparameter nonlinear phase noise mitigation by using M-ary support vector machine for coherent optical systems
CN108667523A (zh) 基于无数据辅助的knn算法的光纤非线性均衡方法
US11456910B1 (en) Carrier-phase recovery system and method
CN101438517A (zh) 光电场接收器、光多值信号接收器以及光传送系统
CN110233653B (zh) 基于加权集成聚类算法的mimo系统的盲多径识别方法及系统
CN108965178A (zh) 基于机器学习的智能自适应均衡器及均衡解调方法
CN108173599A (zh) 基于DENCLUE聚类的Stokes空间相干光调制格式识别方法
CN107359940A (zh) 相位噪声估计的方法和装置
Yu et al. A modified PSO assisted blind modulation format identification scheme for elastic optical networks
CN114513394B (zh) 一种基于注意力机制图神经网络的信号调制格式识别方法、系统、装置及存储介质
CN109547376A (zh) 一种基于调制特性的高斯混合模型的调制识别方法
CN113346957B (zh) 一种oam-qpsk传输的聚类非线性补偿方法
CN103580740A (zh) 一种信号处理方法以及相关设备
US20230291633A1 (en) Probabilistic shaping qam dynamic equalization and digital signal processing method
CN114553315A (zh) 基于CNN-biRNN的光纤非线性均衡方法及系统
CN114513252A (zh) 基于加权的密度峰值聚类算法的概率整形相干光通信系统非线性补偿方法
Zhang et al. Transformer-based detector for OFDM with index modulation
WO2013051244A1 (ja) 信号処理装置及び信号処理方法
CN101496303A (zh) 在无信道状态信息反馈的mimo系统中的功率分配
CN116389207A (zh) 一种基于信号幅度直方图的调制格式识别方法
CN109067683B (zh) 无线通信中盲检测及调制星座优化方法、存储介质
CN105763265A (zh) 一种相干光接收装置中的载波相位估计方法
CN109412670A (zh) Gsm-mbm系统中基于松弛迭代的低复杂度检测方法

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