CN106332290B - 基于可持续充电水声多跳通信系统的资源分配方法 - Google Patents

基于可持续充电水声多跳通信系统的资源分配方法 Download PDF

Info

Publication number
CN106332290B
CN106332290B CN201610753691.1A CN201610753691A CN106332290B CN 106332290 B CN106332290 B CN 106332290B CN 201610753691 A CN201610753691 A CN 201610753691A CN 106332290 B CN106332290 B CN 106332290B
Authority
CN
China
Prior art keywords
base station
relay
time slot
traversal
underwater
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
CN201610753691.1A
Other languages
English (en)
Other versions
CN106332290A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201610753691.1A priority Critical patent/CN106332290B/zh
Publication of CN106332290A publication Critical patent/CN106332290A/zh
Application granted granted Critical
Publication of CN106332290B publication Critical patent/CN106332290B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W72/00—Local resource management
    • H04W72/04—Wireless resource allocation
    • H04W72/044—Wireless resource allocation based on the type of the allocated resource
    • H04W72/0446—Resources in time domain, e.g. slots or frames
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种针对时分复用的水声中继放大转发通信系统的资源分配方法，其中声呐中继配备无线充电设备，可从船载基站获取能量。从信道遍历和容量上界最大的角度出发，设计出一种时隙资源的分配方案。仿真结果表明，利用本发明中的分配策略，可以显著提高系统遍历和容量上界，从而提高了系统吞吐量以及系统的工作效率。

Description

基于可持续充电水声多跳通信系统的资源分配方法
技术领域
本发明涉及水声通信系统的中继能源供应技术领域，尤其是一种可显著提高系统的可达遍历和容量上界的方法。
背景技术
水声通信技术是当前研究的热点问题之一。在水声通信系统中，信号在传输过程中受到窄带、高噪声、长时延传输等因素影响，发送节点发射的信号到达目的节点时，信号的大幅度衰减和畸变使得信息传输的可靠性下降。为了解决上述问题，中继放大传输技术应运而生，中继放大传输作为一种很有应用前景的中继策略，可降低基站到水下传感器的传输距离，能够降低传输损耗带来的不利影响，并且简单、易于配置和扩展。与此同时，能量供应问题也是目前的挑战。
在传统的水声传感器网络中，传感器节点通常采用电池供电。有限的能量会制约网络的工作寿命，并且带来很高的运营成本(例如更换数以百计的节点的电池)；另一方面，节点电池携带的能量非常有限，电池存储能量的多少决定了整个节点的寿命。因此，能量获取成为水声通信网络中需要优先考虑的一个问题，也是系统设计中最大的挑战。在此背景下，本文考虑了一种具有能量收集能力的无线供电的水声通信系统。
引入一种新的无线充电技术，作为水声通信系统的中继能源供应方式。中继可以通过微波无线电力传输装置从基站远程补充能量。相比传统的能源供电方式，该技术可以减少频繁手动更换电池的麻烦，并有更高的吞吐量，更长的元件寿命以及更低的网络运营成本。此外，无线充电还可以根据不同的环境和服务需求控制它的发射功率、波形、占用时间和频率尺寸等。这些明显的优势使得无线供电成为一个有发展前景的新模式。
发明内容
本发明要解决的技术问题为：引入一种新的无线充电技术，作为水声通信系统的中继能源供应方式，相比传统的能源供电方式，该技术可以减少频繁手动更换电池的麻烦，并有更高的吞吐量，更长的元件寿命以及更低的网络运营成本。本发明提出了水声中继系统模型，其中中继采用无线充电技术，并基于系统在可达遍历和容量上边界最大，提出了一种时隙分配方案。该方案与其它方案相比，显著提高了系统的可达遍历和容量上界
本发明采取的技术方案具体为：一种基于充电的水声中继系统基于充电的水声中继系统由一个船载基站、K个水下传感器以及每个传感器对应的水面的声呐中继组成。另外，系统的基站、中继和传感器都是单天线的，中继有一个放大器用来放大从传感器接收到的信号，并且中继端配有可充电电源，可以从基站发送的信号中获取能量。
上述系统的通信方法，其特征在于整个通信过程可分为两个阶段，第一阶段为下行链路传输阶段，即基站给中继充电的过程；第二阶段为上行链路传输阶段，即水下传感器通过中继放大转发将信号发送至基站的过程。基站S到第i个中继Bi的下行链路信道功率增益满足
Figure BDA0001097065400000021
其中复随机变量
Figure BDA0001097065400000022
表示下行链路信道信息；第i个传感器Ui到第i个中继Bi以及第i个中继Bi到基站S的上行链路信道功率增益分别满足和
Figure BDA0001097065400000024
其中复随机变量
Figure BDA0001097065400000025
表示水下传感器到中继的信道信息，复随机变量
Figure BDA0001097065400000026
表示中继到基站的信道信息；采用时分复用的通信方式，避免了各个传感器间的干扰；τ0为下行链路基站给中继的归一化充电时间，τi为上行链路中第i个传感器Ui发送信号到基站的归一化时间，且有
Figure BDA0001097065400000027
上述通信方法，其特征在于在第一个阶段，即下行链路能量传输过程中，控制充电时间长度为τ0。第i个中继接收到的能量为
Figure BDA0001097065400000028
其中
Figure BDA0001097065400000029
表示充电效率，我们令
Figure BDA00010970654000000210
xS为基站在第一个时隙发送的信号，xS为复随机信号并满足E[|xS|2]＝PS，而pS表示基站的发射功率。
上述通信方法，其特征在于在第二阶段，即上行链路传输中，第i个传感器的通信时间长度为τi，yBi表示中继端接收信号，nBi表示中继端接收噪声，且
Figure BDA0001097065400000031
其中
Figure BDA0001097065400000032
为此噪声的方差，则有
Figure BDA0001097065400000033
其中xi表示传感器Ui端发射信号，pUi为第i个传感器Ui的发射功率，且xi～CN(0,pUi)，中继以放大因子αi对接收到信号放大并转发给基站，ySi表示基站接收信号，nSi表示基站端接收噪声，且
Figure BDA0001097065400000034
其中
Figure BDA0001097065400000035
为此噪声的方差，则基站接收到的信号为
其中
Figure BDA0001097065400000037
PBi为第i个中继的发射功率，xU即为传感器Ui端发射信号xi，由Ei的表达式可以得到
Figure BDA0001097065400000038
由γSi和pBi的表达式得到水下传感器到基站通信速率表达式为
Figure BDA0001097065400000039
其中τ＝[τ0,τ1,…,τK]，
Figure BDA00010970654000000311
Figure BDA00010970654000000312
ci＝glinSi 2pUi+nBi 2nSi 2
对上述系统进行最大系统遍历和容量上界最大化的方案设计。系统遍历和容量上界定义为系统和速率的期望能够达到的最大值，系统和速率表达式为
由于水声信道的多变性使得其难以估计，对水声信道进行统计平均，即对g1i取均值，且Rsum(τ)关于g1i为凹函数，根据詹森不等式，遍历和容量满足
Figure BDA00010970654000000410
其中E(ai),E(bi),E(ci)分别为
Figure BDA00010970654000000411
Figure BDA00010970654000000412
E(ci)＝nSi 2pUiE(g1i)+nBi 2nSi 2
这之中，
Figure BDA00010970654000000413
表示第i个中继的充电效率，令
Figure BDA00010970654000000414
xS为基站在第一个时隙发送的信号，xS为复随机信号并满足E[|xS|2]＝PS，而pS表示基站的发射功率；下行链路信道功率增益满足
Figure BDA0001097065400000041
复随机变量
Figure BDA0001097065400000042
表示下行链路信道信息；第i个传感器Ui到第i个中继Bi以及中继Bi到基站S的上行链路信道功率增益分别满足和
Figure BDA0001097065400000044
复随机变量
Figure BDA0001097065400000045
表示水下传感器到中继的信道信息，复随机变量表示中继到基站的信道信息；nBi表示中继端接收噪声；nSi表示基站端接收噪声；pUi为第i个传感器Ui的发射功率。
由于各个水下传感器配备的悬浮到水面的中继的距离近似相等，在考虑大尺度衰落下不同中继接收信噪比取同一个值，即E(bi)/E(ai)为同一值。
令
Figure BDA0001097065400000047
那么这个海洋中继通信系统遍历和容量上界优化问题可以描述为
Figure BDA0001097065400000048
Figure BDA0001097065400000049
τi≥0,i＝1,…,K
由解决此问题得到，在遍历和容量上界最大化下的时间分配为
Figure BDA0001097065400000051
其中，
Figure BDA0001097065400000052
x*是f(x)的唯一解，
Figure BDA0001097065400000054
即为实际时间，其中T表示一个通信周期。
本发明的实行流程：
第一步，基站给中继进行充电；
第二步，水下传感器通过中继放大转发将信号发送至基站。
通过基站与中继上的中央控制单位进行时隙控制，保证按照权利要求5中的时隙分配方案
Figure BDA0001097065400000055
对上述两步骤进行时隙分配，其中，充电过程时隙为
Figure BDA0001097065400000056
通信过程时隙为
Figure BDA0001097065400000057
附图说明
图1为本发明水声中继系统的模型；
图2为本发明TDMA上行链路与下行链路的时间分配；
图3为基于路径损耗变化的不同方案对比图；
具体实施方式
以下结合附图和具体实施例进一步描述。
结合图1所示的水声中继系统，该系统由一个船载基站、K个水下传感器以及每个传感器对应的水面的声呐中继组成，其中K为正整数。另外，系统的基站、中继和传感器都是单天线的，中继有一个放大器用来放大从传感器接收到的信号，并且中继端配有可充电电源，可以从基站发送的信号中获取能量。
整个通信过程可分为两个阶段，第一阶段为下行链路传输阶段，即基站给中继充电的过程；第二阶段为上行链路传输阶段，即水下传感器通过中继放大转发将信号发送至基站的过程。如图1所示，第i个传感器Ui到第i个中继Bi以及第i个中继Bi到基站S的上行链路信道功率增益分别满足
Figure BDA0001097065400000061
和
Figure BDA0001097065400000062
其中复随机变量
Figure BDA0001097065400000063
表示水下传感器到中继的信道信息，复随机变量
Figure BDA0001097065400000064
表示中继到基站的信道信息。
结合图2，这里采用时分复用的通信方式，避免了各个传感器间的干扰；τ0为下行链路基站给中继的归一化充电时间，τi为上行链路中第i个传感器Ui发送信号到基站的归一化时间，且有
Figure BDA0001097065400000065
在第一个阶段，即下行链路能量传输过程中，控制充电时间长度为τ0。第i个中继接收到的能量为
Figure BDA0001097065400000069
其中
Figure BDA00010970654000000610
表示充电效率，我们令xS为基站在第一个时隙发送的信号，xS为复随机信号并满足E[|xS|2]＝PS，而ps表示基站的发射功率。
在第二个阶段，即上行链路传输中，第i个传感器的通信时间长度为τi，yBi表示中继端接收信号，nBi表示中继端接收噪声，且nBi～CN(0,σ1 2)，其中
Figure BDA0001097065400000066
为此噪声的方差，则有
Figure BDA0001097065400000067
其中xi表示传感器Ui端发射信号，pUi为第i个传感器Ui的发射功率，且xi～CN(0,pUi)，中继以放大因子αi对接收到信号放大并转发给基站，ySi表示基站接收信号，nSi表示基站端接收噪声，且nSi～CN(0,σ2 2)，其中
Figure BDA0001097065400000068
为此噪声的方差，则基站接收到的信号为
其中PBi为第i个中继的发射功率，xU即为传感器Ui端发射信号xi，由Ei的表达式可以得到
由γSi和pBi的表达式得到水下传感器到基站通信速率表达式为
其中τ＝[τ0,τ1,…,τK]，
Figure BDA0001097065400000077
Figure BDA0001097065400000078
ci＝g1inSi 2pUi+nBi 2nSi 2
结合图1系统，进行最大系统遍历和容量上界最大化的方案设计。系统遍历和容量上界定义为系统和速率的期望能够达到的最大值，系统和速率表达式为
由于水声信道的多变性使得其难以估计，对水声信道进行统计平均，即对g1i取均值，且Rsum(τ)关于g1i为凹函数，根据詹森不等式，遍历和容量满足
Figure BDA0001097065400000076
其中E(ai),E(bi),E(ci)分别为
Figure BDA0001097065400000079
E(ci)＝nSi 2pUiE(g1i)+nBi 2nSi 2
这之中，
Figure BDA00010970654000000813
表示第i个中继的充电效率，令
Figure BDA00010970654000000814
xS为基站在第一个时隙发送的信号，xS为复随机信号并满足E[|xS|2]＝PS，而pS表示基站的发射功率；下行链路信道功率增益满足
Figure BDA0001097065400000081
复随机变量
Figure BDA0001097065400000082
表示下行链路信道信息；第i个传感器Ui到第i个中继Bi以及中继Bi到基站S的上行链路信道功率增益分别满足
Figure BDA0001097065400000083
和
Figure BDA0001097065400000084
复随机变量
Figure BDA0001097065400000085
表示水下传感器到中继的信道信息，复随机变量
Figure BDA0001097065400000086
表示中继到基站的信道信息；nBi表示中继端接收噪声；nSi表示基站端接收噪声；pUi为第i个传感器Ui的发射功率。
由于各个水下传感器配备的悬浮到水面的中继的距离近似相等，在考虑大尺度衰落下不同中继接收信噪比取同一个值，即E(bi)/E(ai)为同一值。
令
Figure BDA0001097065400000087
那么这个海洋中继通信系统遍历和容量上界优化问题可以描述为
Figure BDA0001097065400000088
τi≥0,i＝1,…,K
由解决此问题得到，在遍历和容量上界最大化下的时间分配为
Figure BDA00010970654000000810
其中，
Figure BDA00010970654000000811
x*是f(x)的唯一解，
Figure BDA00010970654000000812
Figure BDA0001097065400000091
即为实际时间，其中T表示一个通信周期。
根据图1所示系统，结合上述时隙分配方案，得到本发明的实行流程：
第一步，基站给中继进行充电；
第二步，水下传感器通过中继放大转发将信号发送至基站。
通过基站与中继上的中央控制单位进行时隙控制，保证按照权利要求5中的时隙分配方案对上述两步骤进行时隙分配，其中，充电过程时隙为
Figure BDA0001097065400000093
通信过程时隙为
Figure BDA0001097065400000094
根据上述时隙分配方案进行仿真比较，如图3所示，证明了本发明提出的分配策略相比于其他传统方案显著提高了系统遍历和容量上界，从而提高了系统吞吐量以及系统的工作效率。

Claims (2)

1.一种基于充电的水声中继系统的时隙分配方法，其特征在于：进行最大系统遍历和容量上界最大化的方案设计，系统遍历和容量上界定义为系统和速率的期望能够达到的最大值，系统和速率表达式为
Figure FDA0002280475860000011
其中τ＝[τ0,τ1,…,τK]；τ0为下行链路基站给中继的归一化充电时间，τi为上行链路中第i个传感器发送信号到基站的归一化时间，且有
对水声信道进行统计平均，即对g1i取均值，且Rsum(τ)关于g1i为凹函数，根据詹森不等式，遍历和容量满足
Figure FDA0002280475860000013
其中E(ai),E(bi),E(ci)分别为
Figure FDA0002280475860000014
E(ci)＝nSi 2pUiE(g1i)+nBi 2nSi 2
Figure FDA0002280475860000016
Figure FDA0002280475860000017
ci＝g1inSi 2pUi+nBi 2nSi 2
这之中，E[·]和E(·)均表示数学期望，
Figure FDA0002280475860000018
表示第i个中继的充电效率，令
Figure FDA0002280475860000019
xS为基站在第一个时隙发送的信号，xS为复随机信号并满足E[|xS|2]＝PS，而pS表示基站的发射功率；下行链路信道功率增益满足
Figure FDA00022804758600000110
复随机变量
Figure FDA00022804758600000111
表示下行链路信道信息；第i个传感器Ui到第i个中继Bi以及中继Bi到基站S的上行链路信道功率增益分别满足和
Figure FDA00022804758600000113
复随机变量表示水下传感器到中继的信道信息，复随机变量
Figure FDA00022804758600000115
表示中继到基站的信道信息；nBi表示中继端接收噪声；nSi表示基站端接收噪声；pUi为第i个传感器Ui的发射功率；
由于各个水下传感器配备的悬浮到水面的中继的距离近似相等，在考虑大尺度衰落下不同中继接收信噪比取同一个值，即E(bi)/E(ai)为同一值；
令
Figure FDA0002280475860000021
那么这个海洋中继通信系统遍历和容量上界优化问题可以描述为
Figure FDA0002280475860000022
Figure FDA0002280475860000023
τi≥0,i＝1,…,K
由解决此问题得到，在遍历和容量上界最大化下的时间分配为
Figure FDA0002280475860000024
其中，
Figure FDA0002280475860000025
x*是f(x)的唯一解，
Figure FDA0002280475860000026
Figure FDA0002280475860000027
即为实际时间，其中T表示一个通信周期。
2.一种基于充电的水声中继系统传输方法，其特征在于包括如下步骤：
第一步，基站给中继进行充电；
第二步，水下传感器通过中继放大转发将信号发送至基站；
通过基站与中继上的中央控制单位进行时隙控制，按照权利要求1所述时隙分配方法对上述两步骤进行时隙分配，
Figure FDA0002280475860000028
其中，充电过程时隙为
Figure FDA0002280475860000029
通信过程时隙为
Figure FDA00022804758600000210
CN201610753691.1A 2016-08-29 2016-08-29 基于可持续充电水声多跳通信系统的资源分配方法 Active CN106332290B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610753691.1A CN106332290B (zh) 2016-08-29 2016-08-29 基于可持续充电水声多跳通信系统的资源分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610753691.1A CN106332290B (zh) 2016-08-29 2016-08-29 基于可持续充电水声多跳通信系统的资源分配方法

Publications (2)

Publication Number Publication Date
CN106332290A CN106332290A (zh) 2017-01-11
CN106332290B true CN106332290B (zh) 2020-01-24

Family

ID=57788242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610753691.1A Active CN106332290B (zh) 2016-08-29 2016-08-29 基于可持续充电水声多跳通信系统的资源分配方法

Country Status (1)

Country Link
CN (1) CN106332290B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108738025B (zh) * 2018-05-09 2023-04-07 大连声文科技发展有限公司 海上客滚船移动基站带宽动态优化分配方法及系统
CN109361430A (zh) * 2018-09-29 2019-02-19 天津大学 一种多用户水声传感网络大容量协作传输方法
CN110932794B (zh) * 2019-12-05 2021-09-21 哈尔滨工程大学 一种基于全深海auv的无冲突双向通信方法
CN111556429B (zh) * 2020-04-02 2021-08-31 清华大学 一种基于契约的水下声学传感网信息收集方法
CN115131952B (zh) * 2022-06-15 2023-10-24 武汉拓宝科技股份有限公司 一种窄带通信方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104539374A (zh) * 2014-12-09 2015-04-22 桂林电子科技大学 基于水面传感网-水声传感网的混合组网方法
CN104811313A (zh) * 2015-04-03 2015-07-29 浙江大学 基于无线能量传输的最佳能量波束和时间分配设计方法
CN105792366A (zh) * 2016-01-29 2016-07-20 西安交通大学 同道干扰下无线充能中继网络最优传输时隙时间分配方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101071421B1 (ko) * 2009-09-01 2011-10-10 인하대학교 산학협력단 부분 중계기 협력 다중화 기술에서 중계기의 위치에 따른 동적 최적 시간 할당 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104539374A (zh) * 2014-12-09 2015-04-22 桂林电子科技大学 基于水面传感网-水声传感网的混合组网方法
CN104811313A (zh) * 2015-04-03 2015-07-29 浙江大学 基于无线能量传输的最佳能量波束和时间分配设计方法
CN105792366A (zh) * 2016-01-29 2016-07-20 西安交通大学 同道干扰下无线充能中继网络最优传输时隙时间分配方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
无线充电的水声通信系统资源分配方案;曹欢欢 等;《东南大学学报》;20160930;第46卷(第5期);第923-927页 *

Also Published As

Publication number Publication date
CN106332290A (zh) 2017-01-11

Similar Documents

Publication Publication Date Title
CN106332290B (zh) 基于可持续充电水声多跳通信系统的资源分配方法
KR20220027946A (ko) 리피터 방법 및 장치
CN104301984B (zh) D2d蜂窝网络中基于时域半双工中继的功率控制方法
CN102812647A (zh) 利用中继进行收发的宽频局域网的无线通信方法及装置
KR101936019B1 (ko) 협력적 비직교 다중접속 시스템에서 동시 무선 정보 전송 및 전력 공급을 위한 협력적 직접 전송 및 중계 전송 방법
CN109474317B (zh) Mr预处理下硬件损伤大规模mimo双向中继系统功率分配方法
CN102333366A (zh) 在移动通信系统中控制上行链路发送功率的装置和方法
Thai et al. Coordinated direct and relay transmission with interference cancelation in wireless systems
KR101938113B1 (ko) 협력적 비직교다중접속 상향링크에서 조직적 직접 및 중계전송 연계기술
CN111740791A (zh) 一种通信控制方法、装置、通信设备及存储介质
KR102253092B1 (ko) 무선 전력 통신 시스템 및 방법
Zhong et al. User cooperation for enhanced throughput fairness in wireless powered communication networks
Melki et al. System performance of two-way decode-and-forward relaying assisted D2D communication underlaying cellular networks
CN108599911B (zh) 一种无线供电的用户协作系统资源分配方法
Lee et al. The realization of full duplex relay and sum rate analysis in multiuser MIMO relay channel
CN102045822B (zh) 无线中继系统中的下行功率控制方法
Cirik et al. Resource allocation in full-duplex OFDMA systems with partial channel state information
Wei et al. Wireless information and power transfer: Spectral efficiency optimization for asymmetric full-duplex relay systems
Hoc et al. Outage and bit error probability analysis in energy harvesting wireless cooperative networks
CN106411369B (zh) 一种基于能耗和节点分类的中继选择方法
CN106330344B (zh) 水声中继通信系统和速率最大准则资源分配方法
Tölli et al. Mode selection and transceiver design for rate maximization in underlay D2D MIMO systems
Tran et al. On the performance of regenerative relaying for SWIPT in NOMA systems
Li et al. On full-duplex link performance under consideration of error vector magnitude
KR101197613B1 (ko) 공유 중계기 및 이동통신 시스템

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant