CN102932311A - 基于四项加权分数傅里叶变换的带宽压缩方法 - Google Patents

基于四项加权分数傅里叶变换的带宽压缩方法 Download PDF

Info

Publication number
CN102932311A
CN102932311A CN2012104868902A CN201210486890A CN102932311A CN 102932311 A CN102932311 A CN 102932311A CN 2012104868902 A CN2012104868902 A CN 2012104868902A CN 201210486890 A CN201210486890 A CN 201210486890A CN 102932311 A CN102932311 A CN 102932311A
Authority
CN
China
Prior art keywords
signal
estimated
fourier transform
score fourier
new
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
CN2012104868902A
Other languages
English (en)
Other versions
CN102932311B (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.)
Harbin Institute of Technology
Original Assignee
Harbin 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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201210486890.2A priority Critical patent/CN102932311B/zh
Publication of CN102932311A publication Critical patent/CN102932311A/zh
Application granted granted Critical
Publication of CN102932311B publication Critical patent/CN102932311B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

基于四项加权分数傅里叶变换的带宽压缩方法,涉及通信系统的带宽资源节约方法。为了达到在多载波以及混合载波体制下节省带宽的目的,在通信系统的发送端,对长度为N的信号X,利用N×N的加权矩阵Wα进行四项加权分数傅里叶变换,得到加权域信号Y;发送长度为N/2的加权域信号Yh;接收端接收到加权域信号Yh,对信号Yh求矩阵广义逆,获得估计信号
Figure DDA00002467220700011
同时对信道进行信噪比SNR估计;任意选择一组随机序列作为估计信号
Figure DDA00002467220700012
后一半的估计值,该估计值记为
Figure DDA00002467220700013
将信号估计值
Figure DDA00002467220700014
和估计信号
Figure DDA00002467220700015
的前一半组成新的信号估计值
Figure DDA00002467220700017
对新的信号估计值进行四项加权分数傅里叶变换,获得新的加权域信号
Figure DDA00002467220700018
再根据搜索条件,搜索最佳信号,重新组合获得原始信号。本发明适用于通信技术领域。

Description

基于四项加权分数傅里叶变换的带宽压缩方法
技术领域
本发明涉及一种通信系统的带宽资源节约方法,具体涉及基于四项加权分数傅里叶变换的带宽资源压缩方法。
背景技术
在无线通信中,特别是在移动通信近几年发展中,移动电话用户的增长速度远远超过了预测值,需要传输的业务量越来越多,业务类型也越来越多,而无线通信的带宽资源是一种有限的资源,如何合理利用宝贵的带宽资源,这是无线通信网络系统关注的重要问题。
传统的单载波体制下,为了节约带宽,提出了单边带通信技术,它是将一个载频和一个边带加以抑制,只发送一个边带的通信方式。但是现有的以OFDM为代表的多载波通信体制,以及以加权分数傅里叶变换为代表的单载波和多载波的混合载波体制下,由于频谱之间相互混叠,已经不能利用单边带通信的优势,只传一个边带进行通信了,只能将载波频率和上下两个边带一起传输,这样使得发射功率中的大部分功率被消耗掉,并且还要占用较宽的通信带宽,如果在多载波体制下也能实现类似单边带通信的技术,这样就能节省一半的通信带宽。
发明内容
本发明为了达到在多载波以及混合载波体制下节省带宽的目的,提出了基于四项加权分数傅里叶变换的带宽压缩方法。
本发明所述基于四项加权分数傅里叶变换的带宽压缩方法,该方法的的具体实现步骤为:
步骤一、在通信系统的发送端,对长度为N的信号X,利用N×N的加权矩阵Wα进行四项加权分数傅里叶变换4-WFRFT,得到长度为N的加权域信号Y;
步骤二、发送端发送信号Yh,所述信号Yh与为加权信号Y的前一半,所述信号Yh的长度为N/2;
步骤三、接收端接收到加权域信号Yh,对信号Yh求矩阵广义逆,获得估计信号
Figure GDA00002467220500011
同时利用信噪比公式SNR=10log10(Px/Pn)对信道进行信噪比估计,其中Pn是噪声的平均功率,Px是原始信号的平均功率,获得信道的信噪比。
步骤四、任意选择一组随机序列作为估计信号
Figure GDA00002467220500012
后一半的估计值,该估计值记为
步骤五、将步骤四中信号估计值
Figure GDA00002467220500021
和信号
Figure GDA00002467220500022
为估计信号
Figure GDA00002467220500023
的前一半,组成新的信号估计值
Figure GDA00002467220500024
对新的信号估计值进行四项加权分数傅里叶变换4-WFRFT,获得新的加权域信号
Figure GDA00002467220500025
步骤六、对步骤五获得的新的加权域信号
Figure GDA00002467220500026
的前一半
Figure GDA00002467220500027
和步骤三获得的加权域信号Yh进行均方误差比得Δ1,同时设定一个阈值ξ1,当满足Δ11时,执行步骤七,否则返回步骤四,重新对估计信号
Figure GDA00002467220500028
的后一半进行估计;
步骤七、将信号
Figure GDA00002467220500029
为新的加权域信号
Figure GDA000024672205000210
的前一半和步骤三接收到的信号Yh的方差统计量进行比较获得Δ2,并且根据所需信号精度设置阈值ξ2,满足Δ22时,获得最佳接收序列为估计信号
Figure GDA000024672205000211
的后一半
Figure GDA000024672205000212
执行步骤八,否则返回步骤四,重新对估计信号
Figure GDA000024672205000213
的后一半进行估计;
步骤八、将步骤七获得的最佳接收序列为估计信号
Figure GDA000024672205000214
的后一半
Figure GDA000024672205000215
与估计信号
Figure GDA000024672205000216
的前一半信号
Figure GDA000024672205000217
重新组合,获得发送端原始信号X的估计值。
本发明应用四项加权分数傅里叶变换,它是不同于经典傅里叶变换的一种新的变换方法,本发明中,发送端发送的信号的长度是现有技术中所发送信号的长度的一半,因此有效的减少了系统信道带宽的占用,同时利用了符号全搜索算法,相对于传统的BPSK调制方法,在频谱利用率上提高了一倍。
附图说明
图1为具体实施方式一所述基于四项加权分数傅里叶变换的带宽压缩方法过程图。
图2为利用符号全搜索信号估计的误比特率曲线和BPSK调制的理论的误比特率-信噪比曲线对比图,图中:
曲线1是利用符号全搜索算法的信噪比—误比特率曲线,
曲线2是BPSK调制的理论的信噪比—误比特率—曲线。
图3为BPSK调制理论的信噪比-频谱利用率曲线和利用符号全搜索算法的信噪比-频谱利用效率曲线对比图,图中:
曲线3是α=0.15时利用符号全搜索算法的信噪比-频谱利用率比曲线,
曲线4是BPSK调制方法的信噪比-频谱利用率的理论曲线。
具体实施方式
具体实施方式一、结合图1说明本实施方式,本实施方式所述基于四项加权分数傅里叶变换的带宽压缩方法的具体实现方法如下:
步骤一、在通信系统的发送端,对长度为N的信号X,利用N×N的加权矩阵Wα进行四项加权分数傅里叶变换4-WFRFT,得到长度为N的加权域信号Y;
步骤二、发送端发送信号Yh,所述信号Yh与为加权信号Y的前一半,所述信号Yh的长度为N/2;
步骤三、接收端接收到加权域信号Yh,对信号Yh求矩阵广义逆,获得估计信号
Figure GDA00002467220500031
同时利用信噪比公式SNR=10log10(Px/Pn)对信道进行信噪比估计,其中Pn是噪声的平均功率,Px是原始信号的平均功率,获得信道的信噪比。
步骤四、任意选择一组随机序列作为估计信号后一半的估计值,该估计值记为
Figure GDA00002467220500033
所述随机序列是由+1,-1随机组合的、长度为N/2的序列;
步骤五、将步骤四中信号估计值
Figure GDA00002467220500034
和信号
Figure GDA00002467220500035
为估计信号
Figure GDA00002467220500036
的前一半,组成新的信号估计值
Figure GDA00002467220500037
对新的信号估计值进行四项加权分数傅里叶变换4-WFRFT,获得新的加权域信号
Figure GDA00002467220500038
步骤六、对步骤五获得的新的加权域信号
Figure GDA00002467220500039
的前一半
Figure GDA000024672205000310
和步骤三获得的加权域信号Yh进行均方误差比得Δ1,同时设定一个阈值ξ1,当满足Δ11时,执行步骤七,否则返回步骤四,重新对估计信号
Figure GDA000024672205000311
的后一半进行估计;
步骤七、将信号
Figure GDA000024672205000312
为新的加权域信号的前一半和步骤三接收到的信号Yh的方差统计量进行比较获得Δ2,并且根据所需信号精度设置阈值ξ2,满足Δ22时,获得最佳接收序列为估计信号
Figure GDA000024672205000314
的后一半执行步骤八,否则返回步骤四,重新对估计信号
Figure GDA000024672205000316
的后一半进行估计;
步骤八、将步骤七获得的最佳接收序列为估计信号
Figure GDA000024672205000317
的后一半
Figure GDA000024672205000318
与估计信号
Figure GDA000024672205000319
的前一半信号
Figure GDA000024672205000320
重新组合,获得发送端原始信号X的估计值。
本实施方式中,四项加权分数傅里叶变换阶数α的确定可以通过现有技术实现,例如,通过实验用已知的信号作进行仿真,通过改变α的值观察信号的恢复情况,实验统计出,当α∈[0.1,0.3]时,
Figure GDA00002467220500041
前半部分解调出来错误概率小于1%。因此发生错误主要是在
Figure GDA00002467220500042
的后半部分,所以利用符号全搜索方法求信号的后一半。
本实施方式利用的数学工具之一是四项加权分数傅里叶变换,它是不同于经典傅里叶变换的一种新的变换方法,在已有的四项加权分数傅里叶变换的文献中,还未发现将其应用到混合载波体制下单边带通信上的研究;
四项加权分数傅里叶变换4-WFRFT的定义式:
F s α [ f ] ( t ) = A 0 α f ( t ) + A 1 α F ( t ) + A 2 α f ( - t ) + A 3 α F ( - t ) - - - ( 1 )
其中
Figure GDA00002467220500044
表示函数或者信号f(t)的四项加权分数傅里叶变换,A0 α~A3 α是加权系数;
具体可表示为:
A k α = cos ( ( α - k ) π 4 ) cos ( 2 ( α - k ) π 4 ) exp ( 3 ( α - k ) πi 4 ) , ( k = 0,1,2,3 ) - - - ( 2 )
其中,i为虚数单位,α的取值范围是[0,1],当α取0时,四项加权分数傅里叶变换的结果是原始信号f(t),是恒等变换,而当α取1时,四项加权分数傅里叶变换的结果是F(t),退化f(t)的傅里叶变换;其中f(t)和F(t)是互为傅里叶变换对,而且f(-t)是f(t)的关于原点对称的函数,即f(-t)是f(t)以原点为中心的反转,同样道理,F(-t)是F(t)以原点为中心的反转;傅里叶变换采用下列定义式:
F ( w ) = 1 2 &pi; &Integral; - &infin; + &infin; f ( t ) e - jwt dt , - &infin; < w < &infin; - - - ( 3 )
以上四项加权分数傅里叶变换是利用加权系数A0 α~A3 α的定义连续形式,我们可以利用加权矩阵Wα来定义四项加权分数傅里叶变换离散形式,其中Wα定义是
W &alpha; = A 0 &alpha; I + A 1 &alpha; F + A 2 &alpha; I ~ + A 3 &alpha; F ~ - - - ( 4 )
A0 α~A3 α仍然是以上定义的加权系数,I是单位矩阵,F是离散傅里叶变换矩阵,
Figure GDA00002467220500048
Figure GDA00002467220500049
分别是单位矩阵和傅里叶矩阵的循环移位,
Figure GDA000024672205000410
Figure GDA000024672205000411
是组成傅里叶变换基的元素;对于离散信号X=(x0,x1…xN-1)′,(.)′表示向量或矩阵的转置,其中,xk表示第k个信号点,其中k=0,1,...,N-1,由加权矩阵Wα可以定义X的四项加权分数傅里叶变换:
Fα[X]=WαX                                (5)
离散四项加权分数傅里叶变换为模拟信号提供数字变换实现;具体形式如下:
Figure GDA00002467220500051
Figure GDA00002467220500052
Figure GDA00002467220500054
本实施方式的另一个数学工具是矩阵广义逆,矩阵广义逆的定义:
设B∈□m×n,其中□m×n是复数域的m×n矩阵集合,若有一个矩阵G∈□n×m,□n×m表示复数域的n×m矩阵集合,满足以下四个条件的全部或者一部分,则称G为B的广义逆矩阵,简称广义逆。
(1)BGB=B                                (7)
(2)GBG=G                                (8)
(3)(GB)H=GB                             (9)
(4)(BG)H=BG                            (10)
其中(.)H表示矩阵的共轭转置,如果G满足以上四个条件,则称G为B的Moor-Penose(摩尔-彭诺斯)逆;如果B是行满秩矩阵,则B的广义逆G=BH(BBH)-1
具体实施方式二、本实施方式是对具体实施方式一所述基于四项加权分数傅里叶变换的带宽压缩技术的具体实现方法的进一步说明,步骤三所述获得估计信号
Figure GDA00002467220500055
的具体过程为:首先提取加权矩阵的前N行组成新的矩阵为
Figure GDA00002467220500056
然后求出
Figure GDA00002467220500057
的广义逆矩阵,记为
Figure GDA00002467220500058
W h &alpha; + = ( W h &alpha; ) H [ W h &alpha; ( W h &alpha; ) H ] - 1 得到估计信号值
Figure GDA000024672205000510
具体实施方式三、本实施方式是对具体实施方式一所述基于四项加权分数傅里叶变换的带宽压缩技术的具体实现方法的进一步说明,步骤六所述阈值ξ1是由
Figure GDA00002467220500061
确定的,其中,ni表示第i个噪声信号点,N表示信号长度,
Figure GDA00002467220500062
其中,Pn由步骤三的信噪比公式推到获得Pn=Px/10SNR/10
对信号长度N=16,加权矩阵的角度α=0.15进行仿真,ξ1根据噪声的估计值取定,ξ2=10-4,采用BPSK调制,得到信噪比--误比特率曲线2,如图2所示。按照步骤一到步骤八进行操作,得到加权分数傅里叶变换只传一半带宽的信噪比--误比特率曲线1,在信噪比为9dB的时候,符号全搜索算法的误比特率已经达到10-4,曲线1和BPSK调制理论曲线2相差不到1dB,由图3得,当信噪比达到9dB时,符号全搜索算法的频谱利用率能达到2bit.s-1.HZ-1,如曲线3所示,与QPSK调制的频谱利用率一样,相比较于传统的BPSK调制方法,如曲线4所示,频谱利用率提高了一倍。节省了一半的带宽的全搜索方法,本发明没有加任何编码,加上纠错编码,性能会更好。
本发明不局限于上述实施方式,还可以是上述各实施方式中所述技术特征的合理组合。

Claims (3)

1.基于四项加权分数傅里叶变换的带宽压缩方法,其特征在于,该方法的具体实现步骤:
步骤一、在通信系统的发送端,对长度为N的信号X,利用N×N的加权矩阵Wα进行四项加权分数傅里叶变换4-WFRFT,得到长度为N的加权域信号Y;
步骤二、发送端发送信号Yh,所述信号Yh与为加权信号Y的前一半,所述信号Yh的长度为N/2;
步骤三、接收端接收到加权域信号Yh,对信号Yh求矩阵广义逆,获得估计信号
Figure FDA00002467220400011
同时利用信噪比公式SNR=10log10(Px/Pn)对信道进行信噪比估计,其中Pn是噪声的平均功率,Px是原始信号的平均功率,获得信道的信噪比;
步骤四、任意选择一组随机序列作为估计信号
Figure FDA00002467220400012
后一半的估计值,该估计值记为
Figure FDA00002467220400013
步骤五、根据步骤四中信号估计值
Figure FDA00002467220400014
和信号
Figure FDA00002467220400015
为估计信号的前一半,组成新的信号估计值
Figure FDA00002467220400017
对新的信号估计值进行四项加权分数傅里叶变换4-WFRFT,获得新的加权域信号
Figure FDA00002467220400018
步骤六、对步骤五获得的新的加权域信号
Figure FDA00002467220400019
的前一半
Figure FDA000024672204000110
和步骤三获得的加权域信号Yh进行均方误差比得Δ1,同时设定一个阈值ξ1,当满足Δ11时,执行步骤七,否则返回步骤四,重新对估计信号
Figure FDA000024672204000111
的后一半进行估计;
步骤七、将信号
Figure FDA000024672204000112
为新的加权域信号
Figure FDA000024672204000113
的前一半和步骤三接收到的信号Yh的方差统计量进行比较获得Δ2,并且根据所需信号精度设置阈值ξ2,满足Δ22时,获得最佳接收序列为估计信号的后一半
Figure FDA000024672204000115
执行步骤八,否则返回步骤四,重新对估计信号
Figure FDA000024672204000116
的后一半进行估计;
步骤八、将步骤七获得的最佳接收序列为估计信号
Figure FDA000024672204000117
的后一半
Figure FDA000024672204000118
与估计信号的前一半信号
Figure FDA000024672204000120
重新组合,获得发送端原始信号X的估计值。
2.根据权利要求1所述基于四项加权分数傅里叶变换的带宽压缩方法,其特征在于,步骤三所述获得估计信号
Figure FDA000024672204000121
的具体过程为:首先提取加权矩阵的前N行组成新的矩阵为
Figure FDA000024672204000122
然后求出
Figure FDA000024672204000123
的广义逆矩阵,记为 得到估计信号值 X &OverBar; = W h &alpha; + Y h .
3.根据权利要求1所述基于四项加权分数傅里叶变换的带宽压缩方法,其特征在于,骤六所述阈值ξ1是由
Figure FDA00002467220400022
确定,其中,ni表示第i个噪声信号点,N表示信号长度,
Figure FDA00002467220400023
其中,Pn由步骤三的信噪比公式推到获得Pn=Px/10SNR/10
CN201210486890.2A 2012-11-26 2012-11-26 基于四项加权分数傅里叶变换的带宽压缩方法 Active CN102932311B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210486890.2A CN102932311B (zh) 2012-11-26 2012-11-26 基于四项加权分数傅里叶变换的带宽压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210486890.2A CN102932311B (zh) 2012-11-26 2012-11-26 基于四项加权分数傅里叶变换的带宽压缩方法

Publications (2)

Publication Number Publication Date
CN102932311A true CN102932311A (zh) 2013-02-13
CN102932311B CN102932311B (zh) 2014-12-03

Family

ID=47647014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210486890.2A Active CN102932311B (zh) 2012-11-26 2012-11-26 基于四项加权分数傅里叶变换的带宽压缩方法

Country Status (1)

Country Link
CN (1) CN102932311B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108833326A (zh) * 2018-06-25 2018-11-16 哈尔滨工业大学 多分量功率平均的广义混合载波传输方法
WO2022127811A1 (zh) * 2020-12-15 2022-06-23 华为技术有限公司 一种信号传输的方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101232488A (zh) * 2008-02-25 2008-07-30 哈尔滨工业大学 能降低正交频分复用系统峰均功率比的通信方法
CN101340418A (zh) * 2008-08-15 2009-01-07 哈尔滨工业大学 码分多址系统中降低码序列之间干扰及码序列复用的传输方法
US20090285317A1 (en) * 2008-05-19 2009-11-19 Qualcomm Incorporated Methods and systems for effective channel estimation in ofdm systems
CN102340476A (zh) * 2011-07-14 2012-02-01 哈尔滨工业大学 基于加权分数傅里叶变换扩展的ofdm通信系统
CN102624660A (zh) * 2012-04-23 2012-08-01 哈尔滨工业大学 基于四项加权分数傅里叶变换的窄带干扰抑制的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101232488A (zh) * 2008-02-25 2008-07-30 哈尔滨工业大学 能降低正交频分复用系统峰均功率比的通信方法
US20090285317A1 (en) * 2008-05-19 2009-11-19 Qualcomm Incorporated Methods and systems for effective channel estimation in ofdm systems
CN101340418A (zh) * 2008-08-15 2009-01-07 哈尔滨工业大学 码分多址系统中降低码序列之间干扰及码序列复用的传输方法
CN102340476A (zh) * 2011-07-14 2012-02-01 哈尔滨工业大学 基于加权分数傅里叶变换扩展的ofdm通信系统
CN102624660A (zh) * 2012-04-23 2012-08-01 哈尔滨工业大学 基于四项加权分数傅里叶变换的窄带干扰抑制的方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108833326A (zh) * 2018-06-25 2018-11-16 哈尔滨工业大学 多分量功率平均的广义混合载波传输方法
CN108833326B (zh) * 2018-06-25 2021-04-02 哈尔滨工业大学 多分量功率平均的广义混合载波传输方法
WO2022127811A1 (zh) * 2020-12-15 2022-06-23 华为技术有限公司 一种信号传输的方法及装置

Also Published As

Publication number Publication date
CN102932311B (zh) 2014-12-03

Similar Documents

Publication Publication Date Title
CN104601513B (zh) 用于低功率无线网络中通信的方法和布置
CN103117970B (zh) Mimo系统中全双工天线的选择方法
CN101494528B (zh) 发射分集块传输系统的训练序列设计及其信道估计方法
CN103685096B (zh) 一种基于最优导频的mimo-ofdm系统信道估计方法
CN1787413B (zh) 低峰均比正交频分复用传输技术方法
CN103973629B (zh) 一种降低mimo‑ofdm系统papr的slm方法
CN103166685A (zh) Lte中一种基于联合功率分配的干扰对齐方法
CN101651649A (zh) 基于四项加权分数傅立叶变换与分布式发送天线的载波协同通信方法
CN108990142A (zh) 一种低信令开销的多载波非正交多址接入系统的传输方法
CN102215047A (zh) 一种高速数据传输的软扩频多路复用方法
CN103312405B (zh) 一种时频编码分集mt-cdma系统发射与接收方法
CN110460556B (zh) 正交多载波系统无线数据与能量一体化传输信号设计方法
CN104796371B (zh) 一种降低无线信道pmd效应影响的ofdm自适应极化调制方法
US20180054242A1 (en) Mimo communication method, transmitting device, and receiving device
CN102932311B (zh) 基于四项加权分数傅里叶变换的带宽压缩方法
CN112788725B (zh) 无人机通信中一种基于空间调制的非正交多址接入能效优化方法
CN102611650B (zh) 一种广义多载波系统频域信道估计方法及装置
CN101917208B (zh) 一种抑制相邻小区干扰的方法
CN1758576A (zh) 一种空时分组码mc-ds-cdma系统上行链路发射与接收方法
CN100492956C (zh) 一种空时分组码mt-cdma系统上行链路发射与接收方法
CN104639472B (zh) 上行多用户mimo信道估计方法
CN102833209B (zh) 码分多址通信方法
CN105119861A (zh) 一种交织载波索引调制n阶连续ofdm系统ber改善方法
CN102196543A (zh) 一种基于二元功率分配的移动通信基站能量效率优化方法
CN108768471A (zh) 联合选择映射的mimo-ofdm-cdma扩频方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant