CN108173790A - 一种超奈奎斯特信号的传输方法 - Google Patents
一种超奈奎斯特信号的传输方法 Download PDFInfo
- Publication number
- CN108173790A CN108173790A CN201711294563.6A CN201711294563A CN108173790A CN 108173790 A CN108173790 A CN 108173790A CN 201711294563 A CN201711294563 A CN 201711294563A CN 108173790 A CN108173790 A CN 108173790A
- Authority
- CN
- China
- Prior art keywords
- signal
- matrix
- transmission method
- super nyquist
- coefficient
- 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
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03828—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties
- H04L25/03834—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties using pulse shaping
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/40—Transceivers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03343—Arrangements at the transmitter end
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/06—Dc level restoring means; Bias distortion correction ; Decision circuits providing symbol by symbol detection
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Dc Digital Transmission (AREA)
Abstract
一种超奈奎斯特信号的传输方法,涉及高速光信号的发送与接收,包括步骤:S1.发送端发送训练信号序列,经过星座图映射后的符号序列sn分割成多个数据块,每个数据块的首端及尾端分别添加长度为v的循环前缀,再进行超奈奎斯特滤波后发出;S2.接收端去除循环前缀,基于信道矩阵分解所得到的特征矩阵,分别计算sf与yf,再计算出均衡矩阵;S3.发送端正式发送信号序列,按照S1的方式发出信号;S4.接收端收到正式发送的信号,按照S2的方式得到yf,再采用S2中得到的均衡矩阵进行信号均衡,恢复出发送信号。本发明通过超奈奎斯特整形,压缩信号带宽,释放光器件的带宽压力,避免因带宽不足带来的信号串扰。
Description
技术领域
本发明涉及高速光信号的发送与接收,具体涉及一种超奈奎斯特信号的传输方法。
背景技术
随着宽带应用的不断普及与发展,城域光网络以及数据中心光网络需要向更高的速率演进。但是,现有光器件带宽的不足严重限制了光网络提速的步伐,并且这种现状难以在短时间内从技术和成本方面得到完美的解决。而基于现有的器件带宽进行提速,必然会给信号带来严重的串扰。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种超奈奎斯特信号的传输方法,通过超奈奎斯特整形,压缩信号带宽,释放光器件的带宽压力,避免因带宽不足带来的信号串扰。
为达到以上目的,本发明采取一种超奈奎斯特信号的传输方法,包括步骤:
S1.发送端发送训练信号序列,信号序列经过星座图映射后的符号序列sn分割成长度为N的多个数据块,对每个所述数据块的首端及尾端分别添加长度为v的循环前缀,再进行超奈奎斯特滤波后发出信号;
S2.接收端去除收到数据块首端及尾端的所述循环前缀,得到长度为N的数据块的信道矩阵,基于所述信道矩阵分解所得到的特征矩阵,分别计算sf与yf,再计算出均衡矩阵;其中sf与yf分别代表发送端和接收端的数据块经过特征矩阵映射后的数据块对象;
S3.发送端开始正式发送信号序列,按照所述S1的方式发出信号;
S4.接收端收到发送端正式发送的信号,按照所述S2的方式得到yf,再采用S2中得到的均衡矩阵进行信号均衡,恢复出发送信号。
在上述技术方案的基础上,所述S1中超奈奎斯特滤波后的信号y(t)满足公式
其中,Es表示符号的平均功率,sn表示发送符号序列中索引为n的符号,n表示序号,g(t)表示信号载波,μ(t)表示系统噪音,T表示符号周期,α表示符号加速因子,α<1。
在上述技术方案的基础上,所述信号载波g(t)为升余弦函数载波、奈奎斯特载波或其它正交载波。
在上述技术方案的基础上,所述S2中长度为N的数据块表示为
其中,s=[s0,s1,…,sN-1wT表示发送数据块,nNoise表示噪音分量,G表示信道矩阵。
在上述技术方案的基础上,所述信道矩阵G第k行系数在k≤N-2v时,表示为
gk=[01,…,0k-1,g(-vαT),…,g(vαT),01,…,0N-2v-k],
即,行系数分三部分组成,分别是首端k-1个零;中间2v+1个系数表示因信道造成的符号间串扰,对应g(i*αT),i=[-v,…,0,…,v],T为符号周期,α表示符号加速因子;以及尾端的N-2v-k个零。
在上述技术方案的基础上,所述信道矩阵G第k行系数在k>N-2v时,表示为
在上述技术方案的基础上,所述信道矩阵G内部是循环结构,该矩阵分解为
其中Q是特征矩阵,是对角矩阵,且特征矩阵Q表述为
其中qu,r表示特征矩阵中第u行,第r列的系数。
在上述技术方案的基础上,通过特征矩阵Q,得到sf与yf分别为
sf=Q*s
由于表示为所以有计算均衡矩阵W为
w(i,i)=a(i,i)*/la(i,i)|2
其中,a(i,i)是对角矩阵在第i行与第i列的系数,w(i,i)是均衡矩阵W在第i行与第i列的系数,均衡矩阵W在非对角位置的系数为零。
在上述技术方案的基础上,所述S4中,通过均衡矩阵W与特征矩阵Q恢复发送信号
在上述技术方案的基础上,所述发送端发送信号之前进行并串转化;所述接收端接收信号经同步处理,进行串并转化后,在去除所述循环前缀。
本发明的有益效果在于:通过超奈奎斯特整形,可以将信号带宽压缩至原信号带宽的一半以内,因此能够极大的释放光器件的带宽压力,进而实现高速信号的发送与接收,避免因带宽不足带来的信号串扰。同时,压缩信号带宽可以使得频谱利用率增高,本发明实现方式简单,计算复杂度较低。
附图说明
图1为本发明实施例发送端信号处理流程图;
图2为本发明实施例接收端接收训练信号的处理流程图;
图3为本发明实施例接收端信号处理流程图。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
本发明超奈奎斯特(FTN)信号的传输方法,包括步骤:
S1.发送端发送训练信号序列,如图1所示,将信号序列经过星座图映射,得到映射后符号序列sn;将符号序列sn分割成长度为N的多个数据块,每个数据块即:[s0,s1,s2…sN-1]。针对每个数据块,在首端及尾端分别添加长度为v的循环前缀,首端循环前缀为[sN-v,sN-v+1,sN-v+2…sN-1],尾端循环前缀为[s0,s1,s2…sv-1]。对已添加循环前缀的整个符号序列进行超奈奎斯特滤波,使得输出信号y(t)满足以下公式(1):
其中,Es表示符号的平均功率,sn表示发送符号序列中索引为n的符号,n表示序号,g(t)表示信号载波,μ(t)表示系统噪音,T表示符号周期,α表示符号加速因子,通常情况下,α<1。并且,发送端需要经过并串转化后,再发出信号。其中,所述信号载波g(t)可以为升余弦函数载波、奈奎斯特载波或其它正交载波。
S2.接收端去除收到数据块首端及尾端的所述循环前缀,得到长度为N的数据块的信道矩阵,基于信道矩阵分解所得到的特征矩阵,分别计算sf与yf,再计算出均衡矩阵;其中,sf代表发送端待发出的数据块经过特征矩阵映射后的数据块对象,yf代表接收端收到的数据块经过特征矩阵映射后的数据块对象。
S3.发送端开始正式发送信号序列,按照所述S1中的方式发出信号。
S4.接收端收到发送端正式发送的信号,按照所述S2的方式得到yf,再采用S2中得到的均衡矩阵进行信号均衡,恢复出发送信号。
如图2所示,步骤S2中,接收端接收训练信号具体处理流程如下:
S201.接收端收到来自发送端的信号,首先经过信号同步处理。
S202.对同步处理后的序列信号进行串并转化,得到长度为N+2v的多个数据块。
S203.分别去除每个数据块首端及尾端长度为v的循环前缀,得到长度为N的数据块,每块数据块可以表示为:
其中,s=[s0,s1,…,sN-1]T表示发送数据块,nNoise表示噪音分量,G表示信道矩阵。本申请中,采用加粗字体表示矢量。信道矩阵G中,第k行系数在k≤N-2v时,表示为:
gk=[01,…,0k-1,g(-vαT),…,g(vαT),01,…,0N-2v-k],
即,行系数分三部分组成,分别是首端k-1个零;中间2v+1个系数表示因信道造成的符号间串扰,对应g(i*αT),i=[-v,…,0,…,v],T为符号周期,α表示符号加速因子;以及尾端的N-2v-k个零。
信道矩阵G中,第k行系数在k>N-2v时,表示为:
S204.由于信道矩阵G内部是循环结构,该矩阵可以进行如下分解:
其中,Q是特征矩阵;是对角矩阵。特别的,特征矩阵Q有如下表述形式:
其中,qu,r表示特征矩阵中第u行,第r列的系数。
通过特征矩阵Q,计算得到sf与yf,即:
sf=Q*s (5)
由于可以表示为所以有:
根据公式(7)可以基于之前计算得到的sf与yf求得对角矩阵
S205.计算均衡矩阵W:
w(i,i)=a(i,i)*/|a(i,i)|2 (8)
其中,a(i,i)是对角矩阵在第i行与第i列的系数,w(i,i)是均衡矩阵W在第i行与第i列的系数,均衡矩阵W在非对角位置的系数为零。
上述步骤S204和S205形成训练阶段,通过训练阶段得到频域均衡系数w(i,i)。如图3所示,为步骤S4的具体流程,包括如下步骤:
S401.接收端收到发送端正式发送的信号,首先经过信号同步处理。
S402.对同步处理后的序列信号进行串并转化,得到长度为N+2v的多个数据块。
S403.分别去除每个数据块首端及尾端长度为v的循环前缀,得到长度为N的数据块,每块数据块可以表示为:
其中,s=[s0,s1,…,sN-1]T表示发送数据块,nNoise表示噪音分量,G表示信道矩阵。本申请中,采用加粗字体表示矢量。信道矩阵G中,第k行系数在k≤N-2v时,表示为:
gk=[01,…,0k-1,g(-vαT),…,g(vαT),01,…,0N-2v-k],
即,行系数分三部分组成,分别是首端k-1个零;中间2v+1个系数表示因信道造成的符号间串扰,对应g(i*αT),i=[-v,…,0,…,v],T为符号周期,α表示符号加速因子;以及尾端的N-2v-k个零。
信道矩阵G中,第k行系数在k>N-2v时,表示为:
S404.由于信道矩阵G内部是循环结构,该矩阵可以进行如下分解:
其中,Q是特征矩阵;是对角矩阵。特别的,特征矩阵Q有如下表述形式:
其中,qu,r表示特征矩阵中第u行,第r列的系数。
通过特征矩阵Q,得到sf与yf,即:
sf=Q*s (5)
由于可以表示为所以有:
S405.通过训练阶段得到的均衡矩阵W与特征矩阵Q恢复发送信号:
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (10)
1.一种超奈奎斯特信号的传输方法,其特征在于,包括步骤:
S1.发送端发送训练信号序列,信号序列经过星座图映射后的符号序列sn分割成长度为N的多个数据块,对每个所述数据块的首端及尾端分别添加长度为v的循环前缀,再进行超奈奎斯特滤波后发出信号;
S2.接收端去除收到数据块首端及尾端的所述循环前缀,得到长度为N的数据块的信道矩阵,基于所述信道矩阵分解所得到的特征矩阵,分别计算sf与yf,再计算出均衡矩阵;其中sf与yf分别代表发送端和接收端的数据块经过特征矩阵映射后的数据块对象;
S3.发送端开始正式发送信号序列,按照所述S1的方式发出信号;
S4.接收端收到发送端正式发送的信号,按照所述S2的方式得到yf,再采用S2中得到的均衡矩阵进行信号均衡,恢复出发送信号。
2.如权利要求1所述的超奈奎斯特信号的传输方法,其特征在于:所述S1中超奈奎斯特滤波后的信号y(t)满足公式
其中,Es表示符号的平均功率,sn表示发送符号序列中索引为n的符号,n表示序号,g(t)表示信号载波,μ(t)表示系统噪音,T表示符号周期,α表示符号加速因子,α<1。
3.如权利要求2所述的超奈奎斯特信号的传输方法,其特征在于:所述信号载波g(t)为升余弦函数载波、奈奎斯特载波或其它正交载波。
4.如权利要求2所述的超奈奎斯特信号的传输方法,其特征在于:所述S2中长度为N的数据块表示为
其中,s=[s0,s1,…,sN-1]T表示发送数据块,nNoise表示噪音分量,G表示信道矩阵。
5.如权利要求4所述的超奈奎斯特信号的传输方法,其特征在于:所述信道矩阵G第k行系数在k≤N-2v时,表示为
gk=[01,…,0k-1,g(-vαT),…,g(vαT),01,…,0N-2v-k],
即,行系数分三部分组成,分别是首端k-1个零;中间2v+1个系数表示因信道造成的符号间串扰,对应g(i*αT),i=[-v,…,0,…,v],T为符号周期,α表示符号加速因子;以及尾端的N-2v-k个零。
6.如权利要求5所述的超奈奎斯特信号的传输方法,其特征在于:所述信道矩阵G第k行系数在k>N-2v时,表示为
7.如权利要求6所述的超奈奎斯特信号的传输方法,其特征在于:所述信道矩阵G内部是循环结构,该矩阵分解为
其中Q是特征矩阵,是对角矩阵,且特征矩阵Q表述为
其中qu,r表示特征矩阵中第u行,第r列的系数。
8.如权利要求7所述的超奈奎斯特信号的传输方法,其特征在于:通过特征矩阵Q,得到sf与yf分别为
sf=Q*s
由于表示为所以有计算均衡矩阵W为
w(i,i)=a(i,i)*/|a(i,i)|2
其中,a(i,i)是对角矩阵在第i行与第i列的系数,w(i,i)是均衡矩阵W在第i行与第i列的系数,均衡矩阵W在非对角位置的系数为零。
9.如权利要求8所述的超奈奎斯特信号的传输方法,其特征在于:所述S4中,通过均衡矩阵W与特征矩阵Q恢复发送信号
10.如权利要求1-9任一项所述的超奈奎斯特信号的传输方法,其特征在于:所述发送端发送信号之前进行并串转化;所述接收端接收信号经同步处理,进行串并转化后,在去除所述循环前缀。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711294563.6A CN108173790B (zh) | 2017-12-08 | 2017-12-08 | 一种超奈奎斯特信号的传输方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711294563.6A CN108173790B (zh) | 2017-12-08 | 2017-12-08 | 一种超奈奎斯特信号的传输方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108173790A true CN108173790A (zh) | 2018-06-15 |
CN108173790B CN108173790B (zh) | 2020-01-07 |
Family
ID=62525498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711294563.6A Active CN108173790B (zh) | 2017-12-08 | 2017-12-08 | 一种超奈奎斯特信号的传输方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108173790B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110912651A (zh) * | 2019-12-19 | 2020-03-24 | 西安交通大学 | 一种ftn索引调制信号检测方法 |
WO2023245829A1 (zh) * | 2022-06-21 | 2023-12-28 | 苏州大学 | 数据中心光互连系统及方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104394110A (zh) * | 2014-12-19 | 2015-03-04 | 华南师范大学 | 一种时域超奈奎斯特非正交传输导频设计方法 |
US20150070089A1 (en) * | 2013-09-09 | 2015-03-12 | MagnaCom Ltd. | Adaptive nonlinear model learning |
US20150237407A1 (en) * | 2014-02-17 | 2015-08-20 | Electronics And Telecommunications Research Institute | Signal receiving apparatus based on faster than nyquist and signal decoding method thereof |
CN105099968A (zh) * | 2015-07-10 | 2015-11-25 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种多径信道中的超奈奎斯特速率通信系统 |
CN105933255A (zh) * | 2016-04-01 | 2016-09-07 | 中国人民解放军理工大学 | 低复杂度超奈奎斯特传输方法 |
CN106941465A (zh) * | 2017-05-10 | 2017-07-11 | 中国科学院上海高等研究院 | 超奈奎斯特率块的传输方法、发射机、接收机及系统 |
CN107018102A (zh) * | 2017-03-23 | 2017-08-04 | 西安电子科技大学 | 基于Ungerboeck模型的超奈奎斯特信号检测方法 |
CN107204947A (zh) * | 2017-06-21 | 2017-09-26 | 中国科学院上海高等研究院 | 一种ftn预均衡传输方法、发射机、接收机及系统 |
-
2017
- 2017-12-08 CN CN201711294563.6A patent/CN108173790B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150070089A1 (en) * | 2013-09-09 | 2015-03-12 | MagnaCom Ltd. | Adaptive nonlinear model learning |
US20150237407A1 (en) * | 2014-02-17 | 2015-08-20 | Electronics And Telecommunications Research Institute | Signal receiving apparatus based on faster than nyquist and signal decoding method thereof |
CN104394110A (zh) * | 2014-12-19 | 2015-03-04 | 华南师范大学 | 一种时域超奈奎斯特非正交传输导频设计方法 |
CN105099968A (zh) * | 2015-07-10 | 2015-11-25 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种多径信道中的超奈奎斯特速率通信系统 |
CN105933255A (zh) * | 2016-04-01 | 2016-09-07 | 中国人民解放军理工大学 | 低复杂度超奈奎斯特传输方法 |
CN107018102A (zh) * | 2017-03-23 | 2017-08-04 | 西安电子科技大学 | 基于Ungerboeck模型的超奈奎斯特信号检测方法 |
CN106941465A (zh) * | 2017-05-10 | 2017-07-11 | 中国科学院上海高等研究院 | 超奈奎斯特率块的传输方法、发射机、接收机及系统 |
CN107204947A (zh) * | 2017-06-21 | 2017-09-26 | 中国科学院上海高等研究院 | 一种ftn预均衡传输方法、发射机、接收机及系统 |
Non-Patent Citations (2)
Title |
---|
CHENG XU等: "Sub-Symbol-Rate Sampling of Super-Nyquist Signals", 《ECOC 2016,42ND EUROPEAN CONFERENCE AND EXHIBITION ON OPTICAL COMMUNICATIONS》 * |
刘文超: "FTN传输系统干扰信号消除技术研究", 《中国优秀硕士学位论文全文数据库(信息科技辑)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110912651A (zh) * | 2019-12-19 | 2020-03-24 | 西安交通大学 | 一种ftn索引调制信号检测方法 |
WO2023245829A1 (zh) * | 2022-06-21 | 2023-12-28 | 苏州大学 | 数据中心光互连系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108173790B (zh) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20160070615A1 (en) | Fec coding/decoding data processing method and related apparatus | |
CN103546402B (zh) | 一种发送信号的方法、装置和系统 | |
US9467311B2 (en) | Method for sending and receiving signal, and corresponding device and system | |
CN109004979B (zh) | 一种室内可见光单极性ofdm通信系统的实现方法 | |
US8989577B2 (en) | Methods and systems for implementing time-division duplexing in the physical layer | |
CN107204947A (zh) | 一种ftn预均衡传输方法、发射机、接收机及系统 | |
CN102104396A (zh) | 一种基于cs理论的脉冲超宽带通信系统 | |
CN108173790A (zh) | 一种超奈奎斯特信号的传输方法 | |
CN107483378A (zh) | 基于dft的ftn块传输方法、发射机、接收机及系统 | |
CN108111455B (zh) | 基于压缩感知的随钻测井数据v-ofdm声波传输方法 | |
CN106301593B (zh) | 自适应盲偏振解复用处理方法和装置 | |
US9071358B2 (en) | Repeater fiber-coax units | |
DE102018010503B4 (de) | Hochgeschwindigkeitslösungen für kopplungsverbindungen mit unterstützung für zeitkontinuierliche, bandinterne rückkanalkommunikation und proprietäre merkmale | |
CN112511472B (zh) | 一种基于神经网络的时频二阶均衡方法及通信系统 | |
CN106972891B (zh) | 一种数字调制方法、数字传输装置及设备 | |
WO2011029247A1 (zh) | 一种编码/解码的方法及装置 | |
CN104901906B (zh) | 一种相干光通信信道估计方法和系统 | |
CN105553529B (zh) | 一种少模光纤传输系统及其数字信号恢复方法 | |
CN109495415B (zh) | 基于数字余弦变换和分段量化的数字移动前传方法及链路 | |
CN115882959A (zh) | 一种基于机器学习的非线性补偿方法 | |
CN205378159U (zh) | 一种ocdma和ofdm混合的无源光网络系统 | |
CN107666450A (zh) | 一种提高光通信信道估计精度的方法与系统 | |
CN113824663A (zh) | 面向相干光通信系统的信号调制格式识别方法及识别系统 | |
CN108983191B (zh) | Ofdm雷达通信一体化系统的低速信号处理方法 | |
CN110071888B (zh) | 一种高速数传中的快速时间同步方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |