CN109788479B - 一种最小化保密中断概率的分布式协作干扰功率分配方法 - Google Patents
一种最小化保密中断概率的分布式协作干扰功率分配方法 Download PDFInfo
- Publication number
- CN109788479B CN109788479B CN201910146085.7A CN201910146085A CN109788479B CN 109788479 B CN109788479 B CN 109788479B CN 201910146085 A CN201910146085 A CN 201910146085A CN 109788479 B CN109788479 B CN 109788479B
- Authority
- CN
- China
- Prior art keywords
- antenna
- alice
- rho
- power
- cooperative
- 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
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
Description
技术领域
本发明涉及无线通信领域,为一种最小化保密中断概率的分布式协作干扰功率分配方法。
背景技术
无线通信技术的发展使通信的有效性和可靠性得到进一步增强,人们已经越来越重视无线通信中的信息安全问题。电磁波的广播特性是无线通信中信息泄露问题的根源,因此,增强无线通信保密性的一个关键点是保障信息在物理层的安全传输。鉴于此,作为对传统加密技术有力补充的物理层安全传输技术引起了研究人员的广泛关注。
协作干扰技术是目前物理层安全研究的一个方向,协作节点通过发送人工噪声信号使得窃听方受到的干扰大于期望接收方,以此加强无线通信的物理层安全。
2011年3月在IEEE Transactions on Signal Processing,VOL.59,NO.3上发表的题为“Optimal Cooperative Jamming to Enhance Physical Layer Security UsingRelays”的文章研究了协作干扰的功率分配方案,该方案要求协作干扰用户之间共享干扰信号。2013年7月在IEEE Transactions on Information Forensics and Security,VOL.8,NO.7上发表的“Uncoordinated Cooperative Jamming for SecretCommunications”提出了一种多天线分布式协作干扰方案,该方案要求协作用户必须配置多天线。考虑到现实中各协作用户可能仅配置单天线,并且彼此之间无法共享干扰信号,此时上述两种方案均不能满足实际需求。
2017年5月在IEEE Transactions on Vehicular Technology,VOL.66,NO.5上发表的题为“On the Secrecy Rate Maximization With Uncoordinated CooperativeJamming by Single-Antenna Helpers”的文章提出了一种全新的分布式协作干扰方案,该方案考虑到协作干扰的实际应用场景,只要求协作用户配备单天线并独立发送噪声信号,同时以期望接收方作为分配功率的控制中心。然而,同大多数现有工作一样,上述文章关注于保密中断概率约束下的保密速率最大化问题,所得到的功率分配结果也不能保证最优性。
事实上,在实际的通信过程中,保密速率一般是预先设定的,此时以最小化保密中断概率为目标确定功率分配方案更具实际意义,这也是本发明所要解决的问题。
发明内容
本发明的目的在于针对上述现有技术中的问题,提供一种最小化保密中断概率的分布式协作干扰功率分配方法,能够在给定保密速率的情况下最小化保密中断概率,增强无线通信系统的物理层安全传输性能,计算复杂度较低,效率高,能够满足通信的实时性要求。
为了实现上述目的,本发明最小化保密中断概率的分布式协作干扰功率分配方法为:
信道状态信息包括K个单天线协作节点Alice-1,Alice-2,…,Alice-K与单天线期望接收方Bob之间的信道系数hB,k、K个单天线协作节点Alice-1,Alice-2,…,Alice-K与单天线窃听方Eve之间的平均信道增益ΓE,k、单天线发信方Alice-0与单天线期望接收方Bob之间的信道系数hB,0以及单天线发信方Alice-0与单天线窃听方Eve平均信道增益ΓE,0;为单天线期望接收方Bob接收机的加性高斯白噪声功率;γB,0=|hB,0|2,RS为预先设定的保密速率;
所述的参数ρLB与ρUB分别表示搜索过程中ρop所处区间的左右端点,ρop为参数ρ的最佳取值,其能够使作为目标函数的保密中断概率自然对数F(ρ)最小,这里的参数ρ表示信噪比门限,当窃听方的信噪比高于参数ρ时会发生保密中断事件;
步骤2,利用黄金分割法嵌套二分法的两层一维搜索在区间[ρLB,ρUB]中寻找ρop;
步骤3,依据步骤2确定的ρop计算K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率,并将这些发送功率反馈给各协作节点,完成功率分配。
优选的,本发明最小化保密中断概率的分布式协作干扰功率分配方法,采用的通信系统包括一个单天线发信方Alice-0、K个单天线协作节点Alice-1~Alice-K、一个单天线期望接收方Bob和一个单天线窃听方Eve;单天线发信方Alice-0与单天线期望接收方Bob构成通信对,通信过程中单天线发信方Alice-0以固定功率P0发送承载信息的信号,K个单天线协作节点Alice-1,Alice-2,…,Alice-K发送独立的高斯白噪声信号干扰单天线窃听方Eve;K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率由单天线期望接收方Bob确定并反馈给各发送节点;无线信道慢衰落并能被建模为准静态瑞利信道。
所述步骤2的具体操作步骤如下:
步骤2.1,预设搜索精度ε1>0,令ρl=ρLB,ρr=ρUB;步骤2.2,令ρ1=ρr-0.618(ρr-ρl),ρ2=ρl+0.618(ρr-ρl),分别计算出F(ρ1)与F(ρ2),令y1=F(ρ1),y2=F(ρ2);步骤2.3,若y1≥y2,则令ρl=ρ1,ρ1=ρ2,y1=y2,ρ2=ρl+0.618(ρr-ρl),重新计算F(ρ2),并令y2=F(ρ2);否则令ρr=ρ2,ρ2=ρ1,y2=y1,ρ1=ρr-0.618(ρr-ρl),重新计算F(ρ1),并令y1=F(ρ1);步骤2.4,若|ρr-ρl|>ε1,返回步骤2.3;否则令结束搜索;其中,ρl与ρr分别表示搜索过程中ρop所处区间的左右端点;ε1根据精度需求设定。
所述的步骤2.2中由给定的ρ值计算F(ρ)的具体步骤如下:
步骤2.2.1,预设搜索精度ε2>0,令λl=0,λr=αρ;步骤2.2.2,令由λm计算对于与Pk<0的情况,分别令及Pk=0,然后利用得到的Pk(k=1,2,…,K)计算目标函数步骤2.2.3,若f≥0,则令λl=λm;否则令λr=λm;步骤2.2.4,若|λr-λl|≤ε2,结束搜索并令然后由λm计算对于与Pk<0的情况,分别令及Pk=0;否则返回步骤2.2.2;步骤2.2.5,依据步骤2.2.4记录的Pk计算λl与λr分别表示搜索过程中λm所处区间的左右端点;参数Pk代表K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率;ε2根据精度需求设定;表示第k个单天线协作节点Alice-k的功率上限;为单天线窃听方Eve接收机的加性高斯白噪声功率。
所述的步骤1中参数ρLB与ρUB的具体取值为:
所述的步骤3的具体步骤如下:
步骤3.1,预设搜索精度ε3>0,令λl=0,λr=αρop;步骤3.2,令由λm计算对于与Pk<0的情况,分别令及Pk=0,然后由Pk(k=1,2,…,K)计算目标函数步骤3.3,若f≥0,则令λl=λm;否则令λr=λm;步骤3.4,若|λr-λl|≤ε3,结束搜索并令然后由λm计算对于与Pk<0的情况,分别令及Pk=0;否则返回步骤3.2;步骤3.5,将步骤3.4所记录的Pk(k=1,2,…,K)作为各协作节点的发送功率反馈给第k个单天线协作节点Alice-k,完成功率分配;其中,ε3根据精度需求设定。
与现有技术相比,本发明具有如下的有益效果:首先通过引入信噪比门限ρ,将关于功率分配的K维优化问题转化为对ρop的一维搜索问题,简化了求解过程。其次,本发明在对F(ρ)的求解中采用了一维搜索方法,由此形成的本发明两层一维搜索算法具有很低的复杂度,在实际使用中能够满足通信的实时性要求;最后,利用本发明所得到的功率分配方案能够在给定保密速率的情况下使得保密中断概率最小,即本发明方法能够获得最优的功率分配。
附图说明
图1基于分布式协作干扰的SISOSE窃听信道模型;
图2本发明分配方法的流程图;
图3利用本发明的功率分配方法实施协作干扰与传统单发单收无协作干扰的保密中断概率随发送功率变化的对比曲线图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
如图1所示,本发明所采用的通信系统包括一个单天线发信方Alice-0、K个单天线协作节点Alice-1~Alice-K、一个单天线期望接收方Bob和一个单天线窃听方Eve;Alice-0与Bob构成通信对,Alice-1,Alice-2,…,Alice-K发送独立的高斯白噪声信号干扰窃听方Eve;节点Alice-1~Alice-K的发送功率由Bob确定并反馈给各发送节点。
记hB,k与hE,k分别为Alice-k(k=0,1,…,K)到Bob与Eve的信道系数,假定无线信道慢衰落并可被建模为准静态瑞利信道,hB,k与hE,k相互独立且都服从均值为零的复高斯分布,即hB,k~CN(0,ΓB,k),hE,k~CN(0,ΓE,k),k=0,1,…,K。这里的ΓB,k和ΓE,k表示复高斯分布的方差,代表了信道的平均增益。定义γB,k=|hB,k|2,γE,k=|hE,k|2。
在通信过程中,Bob可以获取hB,k(k=0,1,…,K)与ΓE,k(k=0,1,…,K);Alice-0以固定功率P0发送承载信息的信号,与此同时,Alice-k(k=1,2,…,K)通过发送干扰信号来降低Eve的信噪比,以此增强Alice-0的物理层安全传输性能,其中,xk~CN(0,1),Pk为可调的Alice-k的发射功率。若用CB与CE分别表示主信道(Alice-0到Bob)和窃听信道(Alice-0到Eve)的信道容量,则保密容量CS=max(0,CB-CE)。由于实际中无法获得窃听信道的准确状态,需要将其当作随机变量,此时CS也变成随机变量。当保密容量小于预先设定的保密速率(记为RS,RS≥0)时将会发生保密中断事件,通信中希望保密中断概率越小越好。
给定保密速率时最小化保密中断概率的问题本质上是一个多维优化问题,直接求解该优化问题十分耗时,很难满足通信的实时性要求。利用图2所示的本发明功率分配方法可以快速求解该优化问题并得到最优解,具体包含以下步骤:
步骤2,利用黄金分割法嵌套二分法的两层一维搜索在区间[ρLB,ρUB]中寻找ρop;
步骤3,依据步骤2确定的ρop计算Alice-1~Alice-K的发送功率并将其反馈给各协作节点,完成功率分配;
步骤1中ρLB与ρUB的具体取值为:
步骤2具体包括以下步骤:
步骤2.1,预设搜索精度ε1>0,令ρl=ρLB,ρr=ρUB;
步骤2.2,令ρ1=ρr-0.618(ρr-ρl),ρ2=ρl+0.618(ρr-ρl),分别计算出F(ρ1)与F(ρ2),令y1=F(ρ1),y2=F(ρ2);
步骤2.3,若y1≥y2,则令ρl=ρ1,ρ1=ρ2,y1=y2,ρ2=ρl+0.618(ρr-ρl),依照下文所述的步骤2.2.1~步骤2.2.5重新计算F(ρ2),并令y2=F(ρ2);否则令ρr=ρ2,ρ2=ρ1,y2=y1,ρ1=ρr-0.618(ρr-ρl),依照下文所述的步骤2.2.1~步骤2.2.5重新计算F(ρ1),并令y1=F(ρ1);
其中,ε1为ρ的预设搜索精度,可根据精度需求设定。
进一步地,步骤2.2中由给定的ρ值计算F(ρ)的具体步骤如下:
步骤2.2.1,预设搜索精度ε2>0,令λl=0,λr=αρ;
步骤2.2.3,若f≥0,则令λl=λm;否则令λr=λm;
步骤3具体包括以下步骤:
步骤3.1,预设搜索精度ε3>0,令λl=0,λr=αρop;
步骤3.3,若f≥0,则令λl=λm;否则令λr=λm;
步骤3.5,将步骤3.4所记录的Pk(k=1,2,…,K)作为各协作节点的发送功率反馈给Alice-k(k=1,2,…,K),完成功率分配;
其中,ε3为λ的预设搜索精度,可根据精度需求设定。
通过仿真对比有无协作干扰两种情况下保密中断概率(SOP)随信噪比变化的情况。图3给出了相同保密速率条件下保密中断概率随信噪比变化的曲线,其中共仿真了RS=0.25,0.5两种情况。以无协作干扰的SISO系统为参考对象,仿真过程中,限定协作节点数目K=10,主节点发送功率为P0,噪声功率将信噪比定义为P0/σ2。仿真中设σ2恒为1,协作节点的功率上限均为另设ΓB,k=ΓE,k=1(k=0,1,…,K)。随机生成2000次独立信道进行蒙特卡洛仿真。
由图3可见,不论保密速率为0.25还是0.5,利用本发明的功率分配方法实施协作干扰所达到的保密中断概率均小于无协作干扰的SISO系统保密中断概率。
同时可以看出,相较于无协作干扰的SISO系统,信噪比的增大可以显著提升基于本发明的分布式协作干扰通信系统的保密性能。
综上所述,从技术方案理论分析以及仿真结果都能够验证本发明方法在保障无线通信物理层安全传输方面的有效作用。
以上所述仅为本发明的较佳实施例,并不用以对本发明做任何形式上的限定,本领域技术人员应当理解的是,在不脱离本发明精神和原则的条件下,本发明还可以进行若干修改或简单替换,这些修改及替换也均落入由所提交权利要求划定的范围之内。
Claims (3)
1.一种最小化保密中断概率的分布式协作干扰功率分配方法,其特征在于,包括:
P0为单天线发信方Alice-0的发送功率;
步骤1中参数ρLB与ρUB的具体取值为:
信道状态信息包括K个单天线协作节点Alice-1,Alice-2,…,Alice-K与单天线期望接收方Bob之间的信道系数hB,k、K个单天线协作节点Alice-1,Alice-2,…,Alice-K与单天线窃听方Eve之间的平均信道增益ΓE,k、单天线发信方Alice-0与单天线期望接收方Bob之间的信道系数hB,0以及单天线发信方Alice-0与单天线窃听方Eve平均信道增益ΓE,0;为单天线期望接收方Bob接收机的加性高斯白噪声功率;γB,0=|hB,0|2,RS为预先设定的保密速率;
所述的参数ρLB与ρUB分别表示搜索过程中ρop所处区间的左右端点,ρop为参数ρ的最佳取值,其能够使作为目标函数的保密中断概率自然对数F(ρ)最小,这里的参数ρ表示信噪比门限,当窃听方的信噪比高于参数ρ时会发生保密中断事件;
步骤2,利用黄金分割法嵌套二分法的两层一维搜索在区间[ρLB,ρUB]中寻找ρop;
具体操作步骤如下:
步骤2.1,预设搜索精度ε1>0,令ρl=ρLB,ρr=ρUB;步骤2.2,令ρ1=ρr-0.618(ρr-ρl),ρ2=ρl+0.618(ρr-ρl),分别计算出F(ρ1)与F(ρ2),令y1=F(ρ1),y2=F(ρ2);步骤2.3,若y1≥y2,则令ρl=ρ1,ρ1=ρ2,y1=y2,ρ2=ρl+0.618(ρr-ρl),重新计算F(ρ2),并令y2=F(ρ2);否则令ρr=ρ2,ρ2=ρ1,y2=y1,ρ1=ρr-0.618(ρr-ρl),重新计算F(ρ1),并令y1=F(ρ1);步骤2.4,若|ρr-ρl|>ε1,返回步骤2.3;否则令结束搜索;其中,ρl与ρr分别表示搜索过程中ρop所处区间的左右端点;ε1根据精度需求设定;
步骤3,依据步骤2确定的ρop计算K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率,并将这些发送功率反馈给各协作节点,完成功率分配;
具体步骤如下:
2.根据权利要求1所述的最小化保密中断概率的分布式协作干扰功率分配方法,其特征在于:采用的通信系统包括一个单天线发信方Alice-0、K个单天线协作节点Alice-1~Alice-K、一个单天线期望接收方Bob和一个单天线窃听方Eve;单天线发信方Alice-0与单天线期望接收方Bob构成通信对,通信过程中单天线发信方Alice-0以固定功率P0发送承载信息的信号,K个单天线协作节点Alice-1,Alice-2,…,Alice-K发送独立的高斯白噪声信号干扰单天线窃听方Eve;K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率由单天线期望接收方Bob确定并反馈给各发送节点;无线信道慢衰落并能被建模为准静态瑞利信道。
3.根据权利要求1所述的最小化保密中断概率的分布式协作干扰功率分配方法,其特征在于,步骤2.2中由给定的ρ值计算F(ρ)的具体步骤如下:
步骤2.2.1,预设搜索精度ε2>0,令λl=0,λr=αρ;步骤2.2.2,令由λm计算对于与Pk<0的情况,分别令及Pk=0,然后利用得到的Pk(k=1,2,…,K)计算目标函数步骤2.2.3,若f≥0,则令λl=λm;否则令λr=λm;步骤2.2.4,若|λr-λl|≤ε2,结束搜索并令然后由λm计算对于与Pk<0的情况,分别令及Pk=0;否则返回步骤2.2.2;步骤2.2.5,依据步骤2.2.4记录的Pk计算λl与λr分别表示搜索过程中λm所处区间的左右端点;参数Pk代表K个单天线协作节点Alice-1,Alice-2,…,Alice-K的发送功率;ε2根据精度需求设定;表示第k个单天线协作节点Alice-k的功率上限;为单天线窃听方Eve接收机的加性高斯白噪声功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146085.7A CN109788479B (zh) | 2019-02-27 | 2019-02-27 | 一种最小化保密中断概率的分布式协作干扰功率分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146085.7A CN109788479B (zh) | 2019-02-27 | 2019-02-27 | 一种最小化保密中断概率的分布式协作干扰功率分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109788479A CN109788479A (zh) | 2019-05-21 |
CN109788479B true CN109788479B (zh) | 2020-04-28 |
Family
ID=66487590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910146085.7A Active CN109788479B (zh) | 2019-02-27 | 2019-02-27 | 一种最小化保密中断概率的分布式协作干扰功率分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109788479B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111555836B (zh) * | 2020-04-27 | 2021-12-28 | 西安交通大学 | 一种基于多节点协作干扰阻塞的无线隐蔽通信模型及方法 |
CN112153674B (zh) * | 2020-08-17 | 2022-02-11 | 西安交通大学 | 一种通信系统功率分配方法、装置、终端及可读存储介质 |
CN112260796B (zh) * | 2020-10-21 | 2023-02-21 | 三维通信股份有限公司 | 控制干扰信号发射的方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825677A (zh) * | 2014-02-28 | 2014-05-28 | 北京航空航天大学 | 一种骨干网容量受限时的多基站协作传输方法 |
CN105657698A (zh) * | 2015-12-24 | 2016-06-08 | 电子科技大学 | 多小区网络中基于协作干扰的安全传输方法 |
CN106160897A (zh) * | 2016-07-29 | 2016-11-23 | 浙江工业大学 | 一种基于带宽和功率分配的防窃听抗干扰频谱接入方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103227615B (zh) * | 2012-01-31 | 2015-09-09 | 富士通株式会社 | 自动增益控制装置和方法、功率调整装置和无线发射系统 |
CN103475441B (zh) * | 2013-07-26 | 2016-06-08 | 北京邮电大学 | 无线多跳网络中基于簇的协作干扰传输方法 |
CN104581960B (zh) * | 2015-01-26 | 2018-04-17 | 西安交通大学 | 一种保密中断概率约束下的多用户协作干扰功率分配方法 |
US9661503B2 (en) * | 2015-05-14 | 2017-05-23 | King Fahd University Of Petroleum And Minerals | Spectrum-efficient secondary users grouping method for two-tier cognitive radio networks |
-
2019
- 2019-02-27 CN CN201910146085.7A patent/CN109788479B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103825677A (zh) * | 2014-02-28 | 2014-05-28 | 北京航空航天大学 | 一种骨干网容量受限时的多基站协作传输方法 |
CN105657698A (zh) * | 2015-12-24 | 2016-06-08 | 电子科技大学 | 多小区网络中基于协作干扰的安全传输方法 |
CN106160897A (zh) * | 2016-07-29 | 2016-11-23 | 浙江工业大学 | 一种基于带宽和功率分配的防窃听抗干扰频谱接入方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109788479A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109788479B (zh) | 一种最小化保密中断概率的分布式协作干扰功率分配方法 | |
AU2020100505A4 (en) | Physical Layer Security Method Applicable To Non-Orthogonal Multiple Access System | |
CN108495337B (zh) | 基于noma的无线携能通信系统最大安全速率优化方法 | |
CN109640321B (zh) | 基于最优中继选择的协作干扰物理层安全传输方法 | |
CN104581960B (zh) | 一种保密中断概率约束下的多用户协作干扰功率分配方法 | |
CN108880734A (zh) | 量子回溯搜索优化的CCFD-Massive MIMO系统功率分配方法 | |
CN110492996A (zh) | 一种应用于多用户大规模mimo系统中的密钥生成方法 | |
CN109728865B (zh) | 大规模天线阵中基于人工噪声的窃听编码方法 | |
CN109150855B (zh) | 一种优化功率资源的鲁棒性无线通信安全传输方法 | |
CN105657698B (zh) | 多小区网络中基于协作干扰的安全传输方法 | |
CN106131941A (zh) | 协作信道非理想信道估计条件下基于人工噪声的协作网络中功率分配方法 | |
Cui et al. | Energy-efficient RSMA for multigroup multicast and multibeam satellite communications | |
Ara et al. | A zero-sum power allocation game in the parallel Gaussian wiretap channel with an unfriendly jammer | |
CN110071748B (zh) | 一种多发单收系统的人工噪声功率分配方法 | |
Abdulateef et al. | Performance analyses of channel estimation and precoding for massive MIMO downlink in the TDD system | |
CN108366357A (zh) | 基于统计QoS保障的D2D异构蜂窝安全传输方法 | |
Yu et al. | Energy efficiency tradeoff in interference channels | |
Noam et al. | Competitive spectrum management with incomplete information | |
He et al. | Resource Allocation for Throughput Maximization in Cognitive Radio Network with NOMA. | |
CN112134603B (zh) | Mu-mimo系统安全传输发射信号合成方法 | |
Feng et al. | Random caching design for multi-user multi-antenna HetNets with interference nulling | |
Deraz et al. | A Comparative Study of User Pairing Techniques of MIMO-NOMA System in different channel models | |
CN106792899B (zh) | 基于次用户选择的认知无线网络物理层安全传输方法 | |
CN112153674B (zh) | 一种通信系统功率分配方法、装置、终端及可读存储介质 | |
CN113225825B (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 |