CN110049507B - 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法 - Google Patents

无线内容分发网络中基于鞅理论的最优缓冲资源分配方法 Download PDF

Info

Publication number
CN110049507B
CN110049507B CN201910367715.3A CN201910367715A CN110049507B CN 110049507 B CN110049507 B CN 110049507B CN 201910367715 A CN201910367715 A CN 201910367715A CN 110049507 B CN110049507 B CN 110049507B
Authority
CN
China
Prior art keywords
buffer
user
intermediate node
buffer capacity
capacity
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
CN201910367715.3A
Other languages
English (en)
Other versions
CN110049507A (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.)
Nanjing Institute of Technology
Original Assignee
Nanjing Institute of Technology
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 Nanjing Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN201910367715.3A priority Critical patent/CN110049507B/zh
Publication of CN110049507A publication Critical patent/CN110049507A/zh
Application granted granted Critical
Publication of CN110049507B publication Critical patent/CN110049507B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/02Arrangements for optimising operational condition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • H04W72/50Allocation or scheduling criteria for wireless resources
    • H04W72/53Allocation or scheduling criteria for wireless resources based on regulatory allocation policies

Abstract

本发明提供了无线内容分发网络中基于鞅理论的最优缓冲资源分配方法。首先确定缓冲中的积压长度、缓冲区容量、数据到达速率和服务速率之间的关系。运用鞅理论,确定缓冲区溢出概率表达式。给定每个服务用户能容忍的最大缓冲区溢出概率,确定每个服务用户需要的最小缓冲区容量。根据中间节点贡献的缓冲区容量与服务的用户需要的最小缓冲区容量关系,确定用户与边缘设备的连接关系,从而构建最小化缓冲区溢出概率和的优化问题,用注水法,确定最优的缓冲资源分配方案。降低了中间节点缓冲区溢出的概率,提高内容分发的效率,进一步提升无线内容分发网络的分发性能。

Description

无线内容分发网络中基于鞅理论的最优缓冲资源分配方法
技术领域
本发明涉及无线内容分发网络、边缘存储和资源分配技术领域,特别涉及无线内容分发网络中最优缓冲资源分配方法。
背景技术
随着智能终端的快速发展,人们对社交和娱乐,尤其是视频资源的消费呈现出了快速增长的势头。无线网络中视频数据占有越来越重要的地位。思科预测到2021年,无线视频数据将占无线数据量的78%。为了满足人们日益增长的视频消费需求,需要对无线网络的架构进行升级或者重新设计,而对现有网络进行升级改造,造价很高。然而,一方面,网络中数量巨大的视频文件,是由很少一部分流行视频文件组成。另一方面,网络中传输很多重复的视频文件,给回程链路带来了很大的流量压力。为了降低改造现有网络架构的代价,缓解回程链路上的流量压力,无线边缘存储技术就被提出了。无线边缘存储技术是将近期比较受欢迎的视频文件,在非高峰时段,存放到边缘设备,比如小基站、家庭机顶盒中。当用户请求视频文件时,不再通过回程链路进行传输,而是直接通过网络的边缘设备传输给用户,有效减少网络上的流量。
由于现在视频资源总量非常巨大,消费者的喜好也各不相同,边缘设备的存储容量有限,将所有的文件都存在网络边缘中,很不现实。有学者提出将一个受欢迎的文件分成很多小文件块,将部分文件块存在边缘网络中,有用户请求时,直接将部分文件块发给用户。剩余的文件块存在远端服务器中,用户请求时,通过中间节点中继的方式,如图1所示,再转发给目标用户,以加快视频文件的请求速度,提高传输效率。这里,在进行文件块中继时就需要用到中间节点的缓冲区资源了。缓冲(Buffer)资源与缓存(Cache)资源不同,缓冲资源有更快的读写速度,本申请的目的是通过合理分配缓冲区资源,提高内容分发网络的效率,主要解决的是最小化缓冲区溢出概率和的问题。通常,一个中间节点,可能同时为多个用户同时提供中继服务,而中间节点贡献出的缓冲区资源是有限的,在缓冲区可能产生排队和溢出的问题。如何有效分配缓冲区的资源,直接决定了能否实现边缘存储技术。
参考文献1“Scheduling Strategy for Multimedia HeterogeneousHigh-SpeedTrain Networks”(Yun Hu;Hongyan Li;Zheng Chang;Zhu Han;IEEE TRANSACTIONS ONVEHICULAR TECHNOLOGY,VOL.66,NO.4,APRIL2017)的第12页中,在公式38下方定义了hs(sn)即“hs(sn)is the corresponding right eigenvector”,在公式40下方对hA(an)进行了定义,“hA(an)is the corresponding right eigenvector”,即表示对应的右特征向量。
发明内容
本发明提供的无线内容分发网络中基于鞅理论的最优缓冲资源分配方法,目的在于填补无线内容分发网络中的最优缓冲资源分配方法的空白,降低了中间节点缓冲区溢出的概率,提高内容分发的效率。
为了实现上述技术目的,本发明的技术方案包括以下步骤:
S1.首先确定缓冲中的积压长度、缓冲区容量、数据到达速率和服务速率之间的关系。运用鞅理论,确定缓冲区溢出概率表达式。
S2.给定每个服务用户能容忍的缓冲区溢出概率,确定每个服务用户需要的最小缓冲区容量。
S3.根据中间节点贡献的缓冲区容量与服务的用户需要的最小缓冲区容量关系,存在三种情况。
S31.第一种情况,当用户需要的最小缓冲区容量,大于中间节点贡献的缓冲区容量,该用户将被建议直接与边缘设备连接。
S32.第二种情况,当用户需要的缓冲区容量,小于中间节点贡献的缓冲区容量,但所有服务用户需要的缓冲区大于中间节点贡献的缓冲区容量,中间节点需要选择一些用户服务,其他用户直接与边缘设备连接。
S33.第三种情况,当所有服务的用户需要的缓冲区小于中间节点贡献的缓冲区容量,建立优化问题,用注水法,确定最优的缓冲资源分配方案。
步骤S1中,运用鞅理论确定缓冲中的积压长度、缓冲区容量、数据到达速率和服务速率之间的关系。缓冲区溢出概率表达式为:
Figure GDA0002963119910000031
其中,用户在缓冲区中的积压用Qi表示,CB表示中间节点贡献出的缓冲区容量大小,αiCB表示分配给用户
Figure GDA00029631199100000311
的缓冲容量,其中αi,0≤αi≤1表示占中间节点贡献的整个缓冲容量的比例,并且满足
Figure GDA0002963119910000032
用户在时间n请求的数据量表示为a(n);在时间n,从中间节点到用户的服务速率表示为s(n)。
Figure GDA0002963119910000033
表示对
Figure GDA0002963119910000034
求期望。
Figure GDA0002963119910000035
Figure GDA0002963119910000036
表示
Figure GDA0002963119910000037
的谱半径,
Figure GDA0002963119910000038
步骤S2中,若给定用户能容忍的最大缓冲溢出概率为ε,带入公式(7),用户最少需要的缓冲容量比例
Figure GDA0002963119910000039
为:
Figure GDA00029631199100000310
其中,
Figure GDA0002963119910000041
表示
Figure GDA0002963119910000042
的谱半径,
Figure GDA0002963119910000043
CB表示中间节点贡献出的缓冲区容量大小。
Figure GDA0002963119910000044
用户能容忍的最大缓冲溢出概率为ε。
步骤S33中,
最优的缓冲资源分配方案,表达式为:
Figure GDA0002963119910000045
其中,
Figure GDA0002963119910000046
表示分配个用户
Figure GDA0002963119910000047
的缓冲区比例。CB表示中间节点贡献出的缓冲区容量大小。
Figure GDA0002963119910000048
Figure GDA0002963119910000049
表示
Figure GDA00029631199100000410
的谱半径,
Figure GDA00029631199100000411
也就是先给每个用户分配了最小比例
Figure GDA00029631199100000412
剩余的份额,即
Figure GDA00029631199100000413
通过调节λ*,进而调节
Figure GDA00029631199100000414
进行分配,当
Figure GDA00029631199100000415
时结束分配。
Figure GDA00029631199100000416
此时的分配份额
Figure GDA00029631199100000417
能够最小化缓冲区溢出概率和,进一步提升无线内容分发网络的分发性能。
有益效果
本发明运用鞅理论确定缓冲中的积压长度、缓冲区容量、数据到达速率和服务速率之间的关系。进而,确定缓冲区溢出概率表达式。获得的这个概率值是非常紧的,适用于数据量具有突发性质的非泊松到达过程;本发明提出的最优缓冲资源分配方法能够最小化缓冲区溢出概率和,进一步提升无线内容分发网络的分发性能。
附图说明
图1为现有采用中间节点中继分配文件的示意图。
图2为无线内容分发网络中基于鞅理论的最优缓冲资源分配方法。
图3是采用本发明提出的最优缓冲分配算法与均匀分配缓冲区方法,关于缓冲溢出概率和的对比图。
图4是本发明操作步骤流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细描述。
如图1所示,给定一个无线异构网络,包含一个服务器,一个基站,多个中间节点和多个用户。一个中间节点周围有多个服务用户,中间节点为周围用户贡献缓冲资源,提供中继服务,进而增强无线网络内容的转发速度。
1.数据到达模型
假设用户
Figure GDA00029631199100000512
在时间k请求的数据量表示为ai(k)。ai(k)服从具有两个状态
Figure GDA0002963119910000051
的马尔科夫调制开关过程。状态
Figure GDA0002963119910000052
表示没有数据到达,即ai(k)=0,状态
Figure GDA0002963119910000053
表示有数据到达,且ai(k)=Ri。ai(k)的转移矩阵表示为:
Figure GDA0002963119910000054
其中pi表示的是从状态
Figure GDA0002963119910000055
到状态
Figure GDA0002963119910000056
的转移概率,qi表示从状态
Figure GDA0002963119910000057
到状态
Figure GDA0002963119910000058
的转移概率。ai(k)相对应的稳态分布可以表示为
Figure GDA0002963119910000059
在时间[m,n]内,累积的到达数据量可以表示为:
Figure GDA00029631199100000510
其中,Ai(m,n)可以看成是一个二元到达过程。若m=0,我们将用
Figure GDA00029631199100000511
表示。
2.中间节点服务模型
中间节点贡献有限的缓冲容量CB服务周围用户。缓冲区用于缓存从基站发送给用户的数据,再利用近距离优势,采用高速传输技术传输给用户,提高内容分发的效率。在时间k,从中间节点到用户的服务速率可以表示为:
Figure GDA0002963119910000061
其中,B是带宽,N为中间节点同时服务的用户个数,Ptr为传输功率,di为中间节点到用户
Figure GDA0002963119910000062
的距离,l为路径损耗参数,
Figure GDA0002963119910000063
为噪声密度。在时间[m,n]内,中间节点对用户
Figure GDA0002963119910000064
服务的数据量可以表示为:
Figure GDA0002963119910000065
若m=0,我们用Si(n)表示Si(0,n)。从公式(2)可以看出,给定了用户的距离di后,服务速率与时间k无关。下面,我们用中间节点对用户
Figure GDA0002963119910000066
服务的数据量si简化表示si(k)。
同时,我们假设
Figure GDA0002963119910000067
公式(4)表示,用户
Figure GDA0002963119910000068
的服务速率si大于请求数据ai(k)的期望,但是小于ai(k)的峰值。这种情况下,数据在被传输之前,需要在中间节点的缓冲区排队,从而使得缓冲区会有数据量的积压。
3.缓冲溢出概率
因为中间节点只贡献出了有限的缓冲容量,缓冲区太多的积压会产生溢出。而且中间节点为服务的每个用户都预留了一块单独的缓冲区,也就是同一个用户的数据将排在一个队列中。这样的话,每个用户的积压Qi可以表示为:
Figure GDA0002963119910000069
用αiCB表示分配给用户
Figure GDA00029631199100000610
的缓冲容量,其中αi,0≤αi≤1表示占中间节点贡献的整个缓冲容量的比例,并且满足
Figure GDA00029631199100000611
缓冲溢出概率表示为:
Pr(Qi≥(αiCB)) (6)
缓冲中的积压具有随机的特性,很难准确的给出数值。下面采用鞅理论,给出缓冲溢出概率的表达式:
Figure GDA0002963119910000071
其中,其中,用户在缓冲区中的积压用Qi表示,CB表示中间节点贡献出的缓冲区容量大小,αiCB表示分配给用户
Figure GDA0002963119910000072
的缓冲容量,其中αi,0≤αi≤1表示占中间节点贡献的整个缓冲容量的比例,并且满足
Figure GDA0002963119910000073
用户在时间n请求的数据量表示为a(n);在时间n,从中间节点到用户的服务速率表示为s(n)。
Figure GDA0002963119910000074
表示对
Figure GDA0002963119910000075
求期望。
Figure GDA0002963119910000076
Figure GDA0002963119910000077
表示
Figure GDA0002963119910000078
的谱半径,
Figure GDA0002963119910000079
若给定用户能容忍的最大缓冲溢出概率为ε,带入公式(7),用户最少需要的缓冲容量比例为:
Figure GDA00029631199100000710
其中,
Figure GDA00029631199100000711
表示
Figure GDA00029631199100000712
的谱半径,
Figure GDA00029631199100000713
CB表示中间节点贡献出的缓冲区容量大小。
Figure GDA00029631199100000714
用户能容忍的最大缓冲溢出概率为ε。
4.三种情况
下面分三种情况讨论用户最少需要的缓冲容量比例与全部缓冲容量1之间的关系。
第一种情况:当用户需要的最小缓冲区容量,大于中间节点贡献的缓冲区容量,即
Figure GDA0002963119910000081
该用户将被建议直接与边缘设备连接。
第二种情况:当用户需要的缓冲区容量,小于中间节点贡献的缓冲区容量,但所有服务用户需要的缓冲区大于中间节点贡献的缓冲区容量,即
Figure GDA0002963119910000082
Figure GDA0002963119910000083
中间节点需要选择一些用户服务,其他用户直接与边缘设备连接。第二种情况的解决办法,不在本专利的申述范围。
第三种情况:当所有服务的用户需要的缓冲区小于中间节点贡献的缓冲区容量,即
Figure GDA0002963119910000084
进一步建立优化问题,用注水法,确定最优的缓冲资源分配方案。
5.最优缓冲资源分配方法
针对第三种情况,即
Figure GDA0002963119910000085
时,构建如下缓冲区分配问题:
Figure GDA0002963119910000086
Figure GDA0002963119910000087
因为
Figure GDA0002963119910000088
Figure GDA0002963119910000089
(10)问题是一个典型的凸优化问题。构建拉格朗日方程
Figure GDA00029631199100000810
其中,{γ1,...,γN}和λ表示拉格朗日乘子。Karush-Kuhn-Tucher(KKT)条件需要满足如下条件:
Figure GDA00029631199100000811
(11)中的第三个条件,也即是
Figure GDA00029631199100000812
可以得到
Figure GDA00029631199100000813
和γi=0
Figure GDA0002963119910000091
其中,λ*是调节量,使结果满足
Figure GDA0002963119910000092
Figure GDA0002963119910000093
所以结果可进一步写成注水形式:
Figure GDA0002963119910000094
其中,
Figure GDA0002963119910000095
Figure GDA0002963119910000096
也就是先给每个用户分配了最小比例
Figure GDA0002963119910000097
剩余的份额,即
Figure GDA0002963119910000098
通过调节λ*,进而调节
Figure GDA0002963119910000099
进行分配,当
Figure GDA00029631199100000910
时结束分配。此时的分配份额,能够最小化缓冲区溢出概率和,进一步提升无线内容分发网络的分发性能。
图2为最优缓冲资源分配方法的示例。图3为采用本专利提出的最优缓冲分配算法与均匀分配缓冲区方法,关于缓冲溢出概率和的对比图。采用本发明的无线内容分发网络中基于鞅理论的最优缓冲资源分配方法的缓冲溢出概率和明显低于均匀分配法的概率和。
以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。

Claims (1)

1.无线内容分发网络中基于鞅理论的最优缓冲资源分配方法,其特征在于,包括以下步骤:
S1.首先确定缓冲中的积压长度、缓冲区容量、数据到达速率和服务速率之间的关系;运用鞅理论,确定缓冲区溢出概率表达式;
缓冲区溢出概率表达式为:
Figure FDA0002963119900000011
其中,用户在缓冲区中的积压用Qi表示,CB表示中间节点贡献出的缓冲区容量大小,αiCB表示分配给用户ui的缓冲容量,其中αi,0≤αi≤1表示占中间节点贡献的整个缓冲容量的比例,并且满足
Figure FDA0002963119900000012
N表示用户的总个数,其中
Figure FDA0002963119900000013
表示对
Figure FDA0002963119900000014
求期望;a(0)表示时间为0时的请求的初始数据量;
Figure FDA0002963119900000015
用户ui在时间n请求的数据到达速率表示为ai(n);在时间n,从中间节点到用户ui的服务速率表示为si(n);
Figure FDA0002963119900000016
定义为
Figure FDA0002963119900000017
Figure FDA0002963119900000018
表示ai(n)的转移矩阵,用户ui在时间j请求的数据量表示为a(j),
Figure FDA0002963119900000019
表示
Figure FDA00029631199000000110
的谱半径,
Figure FDA00029631199000000111
Figure FDA00029631199000000112
表示对
Figure FDA00029631199000000113
求期望;
S2.给定每个服务用户能容忍的最大缓冲区溢出概率,确定每个服务用户需要的最小缓冲区容量;
给定用户能容忍的最大缓冲溢出概率为ε,带入公式(7),用户最少需要的缓冲容量比例
Figure FDA00029631199000000114
为:
Figure FDA00029631199000000115
其中,
Figure FDA0002963119900000021
S3.根据中间节点贡献的缓冲区容量与服务的用户需要的最小缓冲区容量关系,存在三种情况:
S31.第一种情况,当用户需要的最小缓冲区容量,大于中间节点贡献的缓冲区容量,该用户将被建议与边缘设备直接连接;
S32.第二种情况,当用户需要的缓冲区容量,小于中间节点贡献的缓冲区容量,但所有服务用户需要的缓冲区容量大于中间节点贡献的缓冲区容量,中间节点需要选择一些用户服务,其他用户与边缘设备直接连接;
S33.第三种情况,当所有服务的用户需要的缓冲区容量小于中间节点贡献的缓冲区容量,建立优化问题,用注水法,确定最优的缓冲资源分配方案;
最优的缓冲资源分配方案,表达式为:
Figure FDA0002963119900000022
其中,
Figure FDA0002963119900000023
表示分配用户ui的缓冲区比例;
Figure FDA0002963119900000024
也就是先给每个用户分配了最小比例
Figure FDA0002963119900000025
剩余的份额,即
Figure FDA0002963119900000026
通过调节变量λ*,进而调节
Figure FDA0002963119900000027
进行分配,当
Figure FDA0002963119900000028
时结束分配。
CN201910367715.3A 2019-05-05 2019-05-05 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法 Active CN110049507B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910367715.3A CN110049507B (zh) 2019-05-05 2019-05-05 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910367715.3A CN110049507B (zh) 2019-05-05 2019-05-05 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法

Publications (2)

Publication Number Publication Date
CN110049507A CN110049507A (zh) 2019-07-23
CN110049507B true CN110049507B (zh) 2021-07-20

Family

ID=67280740

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910367715.3A Active CN110049507B (zh) 2019-05-05 2019-05-05 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法

Country Status (1)

Country Link
CN (1) CN110049507B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110557838B (zh) * 2019-09-02 2020-08-18 南京工程学院 无线内容分发网络中基于契约理论的缓冲资源分配方法
CN111885641B (zh) * 2020-07-23 2023-09-12 南京工程学院 一种缓冲辅助中继网络中的缓冲资源分配方法
CN112929217B (zh) * 2021-02-05 2022-07-01 吉林化工学院 一种基于鞅理论的差异化网络流量带宽需求估计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104284201A (zh) * 2014-09-26 2015-01-14 北京奇艺世纪科技有限公司 一种视频内容的处理方法和装置
CN107889160B (zh) * 2017-11-15 2020-03-17 北京邮电大学 一种考虑用户时延的小小区网络边缘部分缓存方法

Also Published As

Publication number Publication date
CN110049507A (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN110267338B (zh) 一种d2d通信中联合资源分配和功率控制方法
CN107682135B (zh) 一种基于noma的网络切片自适应虚拟资源分配方法
CN105191209B (zh) 一种用于调度相同频带内的视频点播流和尽力而为流的方法和装置
CN110049507B (zh) 无线内容分发网络中基于鞅理论的最优缓冲资源分配方法
CN109905918B (zh) 一种基于能效的noma蜂窝车联网动态资源调度方法
US5970062A (en) Method and apparatus for providing wireless access to an ATM network
KR100581085B1 (ko) Mf-tdma 광대역 위성 시스템에서의 자원 할당 장치및 그 방법
CN112954385B (zh) 一种基于控制论和数据驱动的自适应分流决策方法
CN110290507B (zh) 一种d2d通信辅助边缘缓存系统的缓存策略及频谱分配方法
CN111654712A (zh) 适用于移动边缘计算场景的动态自适应流媒体组播方法
Khorov et al. SEBRA: SAND-enabled bitrate and resource allocation algorithm for network-assisted video streaming
CN110557838A (zh) 无线内容分发网络中基于契约理论的缓冲资源分配方法
CN111447511B (zh) 一种带有用户感知体验质量的带宽分配方法
CN113079577A (zh) 基于embb和urllc共存场景下的资源分配方法
CN101711033A (zh) 适合于感知无线电网络的动态信道分配方法及系统
CN101645932A (zh) 一种p2p流媒体启发式数据调度方法
CN111885641B (zh) 一种缓冲辅助中继网络中的缓冲资源分配方法
CN111510882B (zh) 基于用户体验质量的车联网频谱资源分配方法及系统
CN114521027A (zh) 电网资源动态调度的方法、装置、电子设备及存储介质
CN108429919B (zh) 多速率视频在无线网络中的缓存和传输优化方法
KR100960109B1 (ko) 통합 시스템의 자원 할당 방법
CN111492643B (zh) 用于编码缓存的装置和方法
CN110446211B (zh) 一种小蜂窝协作边缘缓存方法及系统
Jiang et al. Call admission control vs. bandwidth reservation: reducing handoff call dropping rate and providing bandwidth efficiency in mobile networks
CN113542124B (zh) 一种d2d缓存网络中信用驱动的协作传输方法

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