CN108599949B - 一种安全密钥速率计算方法及系统 - Google Patents
一种安全密钥速率计算方法及系统 Download PDFInfo
- Publication number
- CN108599949B CN108599949B CN201810883479.6A CN201810883479A CN108599949B CN 108599949 B CN108599949 B CN 108599949B CN 201810883479 A CN201810883479 A CN 201810883479A CN 108599949 B CN108599949 B CN 108599949B
- Authority
- CN
- China
- Prior art keywords
- transmission efficiency
- channel transmission
- confidence interval
- estimation
- system noise
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000005540 biological transmission Effects 0.000 claims abstract description 163
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000003321 amplification Effects 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L9/00—Cryptographic mechanisms or cryptographic arrangements for secret or secure communications; Network security protocols
- H04L9/08—Key distribution or management, e.g. generation, sharing or updating, of cryptographic keys or passwords
- H04L9/0816—Key establishment, i.e. cryptographic processes or cryptographic protocols whereby a shared secret becomes available to two or more parties, for subsequent use
- H04L9/0852—Quantum cryptography
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Theoretical Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
- Storage Device Security (AREA)
Abstract
本发明公开一种安全密钥速率计算方法及系统,所述方法包括:对系统噪声和通道传输效率进行区间估计,得到系统噪声置信区间和通道传输效率置信区间;选取系统噪声置信区间的上边界为系统噪声估计值;若所述通道传输效率置信区间的最大值小于峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值,若所述通道传输效率置信区间的最小值大于峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;最后根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。本发明消除了系统的潜在安全性漏洞,使得协议更加安全。
Description
技术领域
本发明涉及量子密钥分发领域,特别是涉及一种应用于连续变量量子密钥分发系统中的安全密钥速率计算方法及系统。
背景技术
量子密钥分发是量子技术领域的一个重要实际应用,它可以让通信双方通过量子信道和经过认证的经典信道建立一组无条件安全的随机密钥。早期的量子密钥分发系统是基于离散变量,需要单光子源和单光子探测器。与离散变量量子密钥分发系统不同,连续变量量子密钥分发系统采用信号光场的正则分量来编码信息,使用高效率的零差或外差探测器来接收信息,同时该系统与现有的经典光学网络具有良好的兼容性,因此,得到国内外科研人员的广泛关注。
理论上,连续变量量子密钥分发可以被证明是无条件安全的。然而,在实际的系统实现中,由于通道传输效率、系统噪声以及统计误差等因素仍然存在安全性问题。精确地估计系统的各个参数是获得安全密钥的关键。
对连续变量量子密钥分发系统安全密钥的估算主要来自于对量子态的协方差矩阵的估算,而协方差矩阵对通道传输效率、系统噪声比较敏感,我们假设它们是由窃听者Eve引起的,对它们进行准确的估计,就可以估算出Eve 窃取的信息量,除去这部分信息,就可以使得合法的通信双方获得安全的密钥。
由于有限码长的影响,可以得到各个参数的估计区间,也就是置信区间。通过选取合适的参数估计值,可以得到安全密钥速率。传统的参数估计过程,人们认为Eve窃取的信息随着通道传输效率的变化是单调的,我们发现Eve 窃取的信息随着通道传输效率的变化是非单调的(凸函数)。利用传统的参数估计值选取方法会低估Eve窃取信息的上边界,产生安全漏洞。
发明内容
本发明的目的是提供一种安全密钥速率计算方法,提高连续变量量子密钥分发系统的安全性。
为实现上述目的,本发明提供了如下方案:
一种安全密钥速率计算方法,所述方法包括:
对系统噪声进行区间估计,得到系统噪声置信区间;
对通道传输效率进行区间估计,得到通道传输效率置信区间;
选取系统噪声置信区间的上边界为系统噪声估计值;
判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;
若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
可选的,所述安全密钥速率的计算公式为:
其中,N表示一个帧的数据量,β表示逆向数据协调效率,n表示用于提取密钥的数据量,IAB表示发送方与接收方之间的互信息,χBE表示窃取信息, 表示私密放大失败的概率,Tmin表示通道传输效率置信区间的下边界,Tmax表示通道传输效率置信区间的上边界,Tpeak是对应于窃取信息极大值的传输效率,称为峰值传输效率,表示系统噪声置信区间的上边界,表示安全密钥速率的下边界。
Tmin表示通道传输效率置信区间的下边界,Tmax表示通道传输效率置信区间的上边界,T是通道传输效率, erf-1(·)是误差函数,1-δPE为置信度,Vm为调制方差,m为用于参数估计的数据量,为系统噪声的无偏估计,为t的无偏估计。
可选的,所述方法在逆向协调下进行。
可选的,若所述第二判断结果表示所述通道传输效率置信区间的最小值是小于所述峰值传输效率,调整所述调制方差,返回步骤:对通道传输效率进行区间估计,得到通道传输效率置信区间,其中,所述峰值传输效率与所述调整方差相关联。
一种安全密钥速率计算系统,所述系统包括:
第一区间估计模块,用于对系统噪声进行区间估计,得到系统噪声置信区间;
第二区间估计模块,用于对通道传输效率进行区间估计,得到通道传输效率置信区间;
第一选取模块,用于选取系统噪声置信区间的上边界为系统噪声估计值;
第一判断模块,用于判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;
第二选取模块,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
第二判断模块,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
第三选取模块,用于若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
计算模块,用于根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明根据窃取信息随着通道传输效率变化的非单调性,分情况讨论通道传输效率估计值的选取,准确估算窃取信息的上边界,消除了系统的潜在安全性漏洞,使得协议更加安全。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明安全密钥速率计算方法的流程图;
图2为本发明安全密钥速率计算系统的结构连接图;
图3是互信息IAB和窃取信息χBE随通道传输效率T变化的曲线;
图4是有限码长领域的密钥速率随着传输距离变化的曲线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明基于连续变量量子密钥分发系统在逆向协调下窃听方窃取的信息量随通道传输效率变化的非单调性,通过优化参数估计步骤,选取合适的参数估计值,提出了一种更为安全的密钥速率计算方法。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明安全密钥速率计算方法的流程图。如图1所示,所述安全密钥速率计算方法包括:
其中,表示系统噪声置信区间的下边界,表示系统噪声置信区间的上边界,erf-1(·)是误差函数,1-δPE为置信度,m为用于参数估计的数据量,为系统噪声的无偏估计,xi表示发送的数据,yi表示接收的数据。
Tmin表示通道传输效率置信区间的下边界,Tmax表示通道传输效率置信区间的上边界,T是通道传输效率, erf-1(·)是误差函数,1-δPE为置信度,Vm为调制方差,m为用于参数估计的数据量,为系统噪声的无偏估计,为t的无偏估计,xi表示发送的数据,yi表示接收的数据。
步骤13:选取系统噪声置信区间的上边界为系统噪声估计值;
步骤14:判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;峰值传输效率与系统有关,当系统固定了,设定值就固定了,可以通过调整调整方差重新确定峰值传输效率,峰值传输效率即系统中通道传输效率的极大值,如图3所示的Tpeak。
步骤15:若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
步骤16:若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
步骤17:若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
可选的,若所述第二判断结果表示所述通道传输效率置信区间的最小值是小于所述峰值传输效率,调整所述调制方差,重新估计通道传输效率置信区间和计算峰值传输效率。
步骤18:根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
所述安全密钥速率的计算公式为:
其中,N表示一个帧的数据量,β表示逆向数据协调效率,n表示用于提取密钥的数据量,IAB表示发送方与接收方之间的互信息,χBE表示窃取信息, 表示私密放大失败的概率,Tmin表示通道传输效率置信区间的下边界,Tmax表示通道传输效率置信区间的上边界,表示系统噪声置信区间的上边界,表示安全密钥速率的下边界。
其中,上述方法均在在逆向协调下进行。
图2为本发明安全密钥速率计算系统的结构连接图。如图2 所示,所述安全密钥速率计算系统包括:
第一区间估计模块21,用于对系统噪声进行区间估计,得到系统噪声置信区间;
第二区间估计模块22,用于对通道传输效率进行区间估计,得到通道传输效率置信区间;
第一选取模块23,用于选取系统噪声置信区间的上边界为系统噪声估计值;
第一判断模块24,用于判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;
第二选取模块25,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
第二判断模块26,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
第三选取模块27,用于若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
计算模块28,用于根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
我们以GG02协议为例,在量子信号传输之后,Alice和Bob获得一帧的相互关联的数据,Alice端的数据记为x,Bob端的数据记为y,则有以下关系式:
y=tx+z (1)
假设一个帧的数据量为N,其中m个数据用于参数估计,则根据数理统计方法,t和σ2具有如下的无偏估计:
它们是独立的估计值,分别具有如下分布:
χBE随着T的变化并非单调函数,是个凸函数,存在一个极大值点,我们记这个点对应的T为Tpeak,有:
当Tmin和Tmax不满足上述两种情况,可以适当调整调制方差,使得它们满足条件。因此:
当Tmax<Tpeak,为了估计Eve窃取信息量的最大值,我们应该选择T的置信区间的上边界Tmax。
当Tmin>Tpeak,为了估计Eve窃取信息量的最大值,我们应该选择T的置信区间的下边界Tmin。
根据上面参数估计值的选取结果,可以估算出安全密钥速率的下边界。
在有限码长领域,密钥速率:
因此安全密钥速率的下界为:
选取合适的数据进行参数估计,利用在图4中,我们给出了最终的安全密钥速率随着传输距离L变化的仿真结果,从左到右, N=106,N=107,N=108,N=109,N=1010以及N=∞(极限情况),其他参数设定为β=0.97,VM=3,ε=0.01,
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种安全密钥速率计算方法,其特征在于,所述方法包括:
对系统噪声进行区间估计,得到系统噪声置信区间;
对通道传输效率进行区间估计,得到通道传输效率置信区间;
选取系统噪声置信区间的上边界为系统噪声估计值;
判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;
若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
7.根据权利要求1所述的安全密钥速率计算方法,其特征在于,所述方法在逆向协调下进行。
8.根据权利要求3所述的安全密钥速率计算方法,其特征在于,所述安全密钥速率计算方法还包括:
若所述第二判断结果表示所述通道传输效率置信区间的最小值小于所述峰值传输效率,调整所述调制方差,返回步骤:对通道传输效率进行区间估计,得到通道传输效率置信区间;其中,所述峰值传输效率与所述调制方差相关联。
9.一种安全密钥速率计算系统,其特征在于,所述系统包括:
第一区间估计模块,用于对系统噪声进行区间估计,得到系统噪声置信区间;
第二区间估计模块,用于对通道传输效率进行区间估计,得到通道传输效率置信区间;
第一选取模块,用于选取系统噪声置信区间的上边界为系统噪声估计值;
第一判断模块,用于判断所述通道传输效率置信区间的最大值是否小于峰值传输效率,得到第一判断结果;
第二选取模块,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值小于所述峰值传输效率,选取通道传输效率置信区间的上边界为通道传输效率估计值;
第二判断模块,用于若所述第一判断结果表示所述通道传输效率置信区间的最大值大于所述峰值传输效率,判断所述通道传输效率置信区间的最小值是否大于所述峰值传输效率,得到第二判断结果;
第三选取模块,用于若所述第二判断结果表示所述通道传输效率置信区间的最小值大于所述峰值传输效率,选取通道传输效率置信区间的下边界为通道传输效率估计值;
计算模块,用于根据所述系统噪声估计值和所述通道传输效率估计值计算安全密钥速率的下边界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883479.6A CN108599949B (zh) | 2018-08-06 | 2018-08-06 | 一种安全密钥速率计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883479.6A CN108599949B (zh) | 2018-08-06 | 2018-08-06 | 一种安全密钥速率计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108599949A CN108599949A (zh) | 2018-09-28 |
CN108599949B true CN108599949B (zh) | 2020-06-12 |
Family
ID=63623092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810883479.6A Expired - Fee Related CN108599949B (zh) | 2018-08-06 | 2018-08-06 | 一种安全密钥速率计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108599949B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114614975B (zh) * | 2020-12-03 | 2024-01-23 | 科大国盾量子技术股份有限公司 | 量子通信中安全密钥速率的计算方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106100835A (zh) * | 2016-07-19 | 2016-11-09 | 北京邮电大学 | 一种用于连续变量量子密钥分发中的高精度参数估计及高后处理数据利用率方法 |
EP3244566A1 (en) * | 2016-05-11 | 2017-11-15 | Institut Mines-Telecom | Phase reference sharing schemes for continuous-variable quantum cryptography |
CN108199839A (zh) * | 2018-02-02 | 2018-06-22 | 北京邮电大学 | 一种用所有数据进行参数估计和密钥提取的方法 |
CN108306733A (zh) * | 2018-05-04 | 2018-07-20 | 北京邮电大学 | 一种用于离散调制连续变量量子密钥分发中的后处理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105024809B (zh) * | 2015-07-22 | 2018-05-04 | 上海交通大学 | 基于高斯调制相干态的长距离连续变量量子密钥分发方法 |
-
2018
- 2018-08-06 CN CN201810883479.6A patent/CN108599949B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3244566A1 (en) * | 2016-05-11 | 2017-11-15 | Institut Mines-Telecom | Phase reference sharing schemes for continuous-variable quantum cryptography |
CN106100835A (zh) * | 2016-07-19 | 2016-11-09 | 北京邮电大学 | 一种用于连续变量量子密钥分发中的高精度参数估计及高后处理数据利用率方法 |
CN108199839A (zh) * | 2018-02-02 | 2018-06-22 | 北京邮电大学 | 一种用所有数据进行参数估计和密钥提取的方法 |
CN108306733A (zh) * | 2018-05-04 | 2018-07-20 | 北京邮电大学 | 一种用于离散调制连续变量量子密钥分发中的后处理方法 |
Non-Patent Citations (2)
Title |
---|
基于连续变量的量子密钥分发系统中长距离传输技术与安全性研究;王天一;《中国博士学位论文全文数据库 信息科技辑》;20160315;I136-41 * |
连续变量量子密钥分发协议的安全性证明及协议改进;张一辰;《中国博士学位论文全文数据库 信息科技辑》;20180215;I136-103 * |
Also Published As
Publication number | Publication date |
---|---|
CN108599949A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112448814A (zh) | 面向电网广域协调控制的非对称双场量子密钥分发方法 | |
CN109902506B (zh) | 一种多隐私预算的本地差分隐私数据分享方法和系统 | |
CN110798227B (zh) | 模型预测优化方法、装置、设备及可读存储介质 | |
CN111525955B (zh) | 一种基于稀疏贝叶斯学习的可见光通信的均衡方法及系统 | |
CN108599949B (zh) | 一种安全密钥速率计算方法及系统 | |
CN113051628B (zh) | 基于残差学习的芯片侧信道攻击降噪预处理方法 | |
CN114595830B (zh) | 一种面向边缘计算场景下的隐私保护联邦学习方法 | |
CN103259662A (zh) | 一种新的基于整数分解问题的代理签名及验证方法 | |
Banelli | Bayesian estimation of a Gaussian source in Middleton's class-A impulsive noise | |
CN112865891A (zh) | Pucch信号检测方法、装置、通信设备和存储介质 | |
Seike et al. | Fork rate-based analysis of the longest chain growth time interval of a pow blockchain | |
CN109635349A (zh) | 一种噪声增强最小化克拉美罗界的方法 | |
CN105959095B (zh) | 一种混沌加密通信中驱动与响应系统的同步方法 | |
Chen et al. | Blind nonlinear system identification based on a constrained hybrid genetic algorithm | |
CN115834085A (zh) | 一种生成签名的方法和装置、安全认证方法和装置 | |
Kohli et al. | Variable forgetting factor LS algorithm for polynomial channel model | |
CN111970287B (zh) | 往返式连续变量量子密钥分配噪声补偿方法及其系统 | |
Long et al. | Bit error rate improvement for chaos shift keying chaotic communication systems | |
CN111901113B (zh) | 一种基于高斯密集调制的连续变量量子密钥分发方法 | |
Hou et al. | Affine projection algorithms based on sigmoid cost function | |
CN111246463B (zh) | 基于alpha稳定分布随机过程偏斜度参数的隐蔽通信方法 | |
Kazemi et al. | Covert communications through mobile voice channels | |
KR100994848B1 (ko) | 광 부호 분할 다중 접속 시스템에서의 고속 부호 획득 방법및 그 장치 | |
CN113806775A (zh) | 一种基于卷积优化的区块链报文处理方法及装置 | |
CN108768610B (zh) | 一种基于混沌系统的bpsk信号加密方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200612 Termination date: 20210806 |
|
CF01 | Termination of patent right due to non-payment of annual fee |