CN108965174B - 大规模mimo系统上行链路的联合信道估计和数据解调方法 - Google Patents
大规模mimo系统上行链路的联合信道估计和数据解调方法 Download PDFInfo
- Publication number
- CN108965174B CN108965174B CN201810748828.3A CN201810748828A CN108965174B CN 108965174 B CN108965174 B CN 108965174B CN 201810748828 A CN201810748828 A CN 201810748828A CN 108965174 B CN108965174 B CN 108965174B
- Authority
- CN
- China
- Prior art keywords
- matrix
- phase noise
- vector
- antennas
- channel
- 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
- 238000000034 method Methods 0.000 title claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 42
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0452—Multi-user MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0851—Joint weighting using training sequences or error signal
-
- 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/03891—Spatial equalizers
- H04L25/03898—Spatial equalizers codebook-based design
- H04L25/0391—Spatial equalizers codebook-based design construction details of matrices
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
- H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
- H04L27/36—Modulator circuits; Transmitter circuits
- H04L27/362—Modulation using more than one carrier, e.g. with quadrature carriers, separately amplitude modulated
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
- H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
- H04L27/38—Demodulator circuits; Receiver circuits
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L5/00—Arrangements affording multiple use of the transmission path
- H04L5/0001—Arrangements for dividing the transmission path
- H04L5/0014—Three-dimensional division
- H04L5/0023—Time-frequency-space
- H04L5/0025—Spatial division following the spatial signature of the channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L5/00—Arrangements affording multiple use of the transmission path
- H04L5/003—Arrangements for allocating sub-channels of the transmission path
- H04L5/0048—Allocation of pilot signals, i.e. of signals known to the receiver
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
Abstract
本发明属无线通信技术领域,涉及一种大规模MIMO系统上行链路的联合信道估计和数据解调方法。本发明采用了变分贝叶斯推断算法,变分贝叶斯推断算法是一种求解未知随机变量的后验分布的算法,通过不断地迭代,得到样本已知的条件下的隐藏变量的均值与方差本发明能够在相位噪声存在的条件下实现对大规模MIMO系统上行链路实现准确的信道和数据的估计,同时该算法复杂度相对较低,与传统变分贝叶斯推断算法相比,将信道估计的协方差矩阵的求逆运算转化为对角矩阵的求逆,同时利用相对较为简单的ZF合并的方式来判决数据符号,来避免把数据符号当作随机变量时计算协方差矩阵的求逆运算。而相位噪声的协方差矩阵在求逆时由于维度较小,计算复杂度可以忽略。
Description
技术领域
本发明属于无线通信技术领域,涉及在相位噪声存在的情况下,利用基于变分贝叶斯推断算法对大规模MIMO系统上行链路进行联合信道估计和数据解调的快速算法。
背景技术
在现代无线通信系统中,大规模MIMO系统由于其较高的频谱效率和能量效率而被广泛认为是下一代移动通信的核心技术,通常情况下,基站拥有几百根天线,可以在同时同频的条件下为数十个用户服务,从而显著提高频谱效率。随着基站天线数的增加,大规模MIMO的天线增益可以使每个用户的发送信号的功率显著降低,从而提高能量效率。
为了充分发挥大规模MIMO系统的优势,信道状态信息需要提前获得,然后在获得信道信息的情况下进行数据的判决和解调,然而,在大规模MIMO系统中进行信道和数据的估计往往是困难的,目前已经存在的一些算法都具有较高的复杂度,尤其是在相位噪声存在的情况下。大规模MIMO通信系统的信号在传输过程中,除了经历信道的衰落以外,还要受到射频器件非线性因素的影响,这两个因素使在接收端系统的性能降低。通信系统中射频前端的非理想部分主要包括相位噪声,IQ幅度相位不平衡,功率放大器非线性失真等,相位噪声,实际上是对频率源频率稳定度的一种表征。通常情况下,频率稳定度分为长期频率稳定度和短期频率稳定度。所谓短期频率稳定度,是指由随机噪声引起的相位起伏或频率起伏。至于因为温度、老化等引起的频率慢漂移,则称之为长期频率稳定度。通常主要考虑的是短期稳定度问题,可以认为相位噪声就是短期频率稳定度,只不过是一个物理现象的两种不同表示方式。对于振荡器,频率稳定度是它在整个规定的时间范围内产生相同频率的一种量度。如果信号频率存在瞬时的变化,不能保持不变,那么信号源就存在着不稳定性,起因就是相位噪声。在大规模MIMO通信系统中,发送端与接收端都需要产生相应的载波以完成相应的射频与基带间的频谱转换。然而产生载波的晶体振荡器与锁相环存在一定的差异性,造成了载波频率与目标频率存在短时的随机差异,进而造成所产生的正弦波信号发生随机相位跳变,表现为相位噪声。对于正交频分的调制方式,相位噪声会产生公共相位误差和载波间干扰,这将严重影响系统的性能。
发明内容
本发明的目的在于在相位噪声存在的情况下,提供一种针对大规模MIMO-OFDM系统上行链路的信道估计和数据解调的快速算法,具有较低的复杂度,便于在硬件上实现。
本发明采用了变分贝叶斯推断算法,变分贝叶斯推断算法是一种求解未知随机变量的后验分布的算法,通过不断地迭代,得到样本已知的条件下的隐藏变量的均值与方差。
为了便于本领域内技术人员对本发明技术方案的理解,首先对本发明采用的系统模型进行说明。
考虑带有相位噪声的MIMO OFDM系统上行链路的模型,发射端有K个用户,每个用户有1根天线,接收端基站有M根天线,发射端第k个用户和接收端第m根天线之间的时域信道矢量记为其中L为信道矢量的长度。对于每个OFDM符号,接收端第m根天线的时域信号表达式为
其中,是第m根天线上的时域接收信号,N是OFDM子载波的个数,是接收端第m根天线的相位噪声矩阵,是第k个用户到接收端第m根天线之间的Toeplitz信道矩阵,它的第1列为其中01×(N-L)表示元素全为0、长度为N-L的行矢量。F∈CN×N是归一化的FFT矩阵,它的第i行第j个元素为dk=[dk,1,dk,2,…,dk,N]T是第k个用户发送的数据或者导频序列。是时域的复高斯白噪声序列,
其中Hm,k=diag{[Hm,k,1,Hm,k,2,…,Hm,k,N]T},
将(3)改写为
将(5)进一步改写为
其中1表示长度为N的全1列向量,θm=[θm,1,θm,2,…,θm,N]T为实高斯分布的相位噪声矢量,即θm=N(0,Φ)。由于θm的协方差矩阵Φ为实对称矩阵,其特征值是实数,并且可以用正交矩阵进行相似对角化:
Φ=UΛUT (8)
其中Λ=diag{[λ1,λ2,…,λN]T}是对角矩阵,对角元素为Φ的降序排列的特征值,U是正交矩阵,它的每一列是Λ对应列的特征值的特征向量。通过计算可以发现,Λ中的对角元只有前若干项的值较大,其它的元素和前若干项相比很小,因此可以只取前I项来近似,即
Φ≈VΓVT(9)
Γ=diag{[λ1,λ2,…,λI]T}是以Λ中前I个特征值为对角元素的对角矩阵,V∈CN×I是由前U的前I列组成的矩阵。对相位噪声矢量θm作线性变换
θm=Ux'm≈Vxm (10)
由高斯分布的性质可知,xm=N(0,Γ),由于Γ为对角阵,所以xm的各个分量之间是相互独立的。把(10)代入(7)可得
hm服从复高斯的先验分布
在传统的变分贝叶斯推断中,计算信道矢量和数据符号的协方差矩阵时都要涉及大维度矩阵的求逆运算,这将大大增加算法复杂度,不利于在硬件上实现。本发明通过一定的近似手段来进行化简,从而使算法的复杂度显著降低。
本发明通过如下步骤实现:
S1、通过以下步骤的迭代来实现对信道矢量的估计:
S11、假定相位噪声为0,作为迭代的初始值;
S12、计算信道矢量的后验分布的协方差矩阵:
S14、计算相位噪声展开矢量的后验分布的均值和方差:
S15、对信道矢量的先验协方差矩阵D进行更新;
S16、循环步骤S21—S23,在已知接收信号的条件下,信道矢量将收敛于一个稳定的值;
S2、通过以下步骤的迭代来实现对数据符号的估计:
S21、计算相位噪声展开矢量的后验分布的均值和方差
S22、利用估计的相位噪声对接收信号进行补偿,再进行ZF合并,利用最大似然判决得到数据符号的估计值
S23、循环步骤S21—S22,在已知接收信号的条件下,数据符号矢量将收敛于一个稳定的值。
本发明的有益效果为能够在相位噪声存在的条件下实现对大规模MIMO系统上行链路实现准确的信道和数据的估计,同时该算法复杂度相对较低,与传统变分贝叶斯推断算法相比,将信道估计的协方差矩阵的求逆运算转化为对角矩阵的求逆,同时利用相对较为简单的ZF合并的方式来判决数据符号,来避免把数据符号当作随机变量时计算协方差矩阵的求逆运算。而相位噪声的协方差矩阵在求逆时由于维度较小,计算复杂度可以忽略。
附图说明
图1是本发明使用的相位噪声影响下的大规模MIMO系统上行链路示意图;
图2是本发明使用的信道模型图;
图3是本发明实现信道估计算法的流程图;
图4是本发明的方法下的BER性能曲线图。
具体实施方式
下面结合附图对本发明进行详细的描述:
S1、通过以下步骤的迭代来实现对信道矢量的估计:
S11、假定相位噪声为0,作为迭代的初始值;
S12、计算信道矢量的后验分布的协方差矩阵:
S14、计算相位噪声展开矢量的后验分布的均值和方差:
S16、循环步骤S21—S23,在已知接收信号的条件下,信道矢量将收敛于一个稳定的值;
S2、通过以下步骤的迭代来实现对数据符号的估计:
S21、计算相位噪声展开矢量的后验分布的均值和方差
S22、利用估计的相位噪声对接收信号进行补偿,再进行ZF合并,利用最大似然判决得到数据符号的估计值
S23、循环步骤S21—S22,在已知接收信号的条件下,数据符号矢量将收敛于一个稳定的值。
图4是利用本发明的算法进行信道和数据符号估计后的系统BER曲线图,调制方式采用64QAM,相位噪声水平为-85dBc/Hz@1MHz,天线分组数为8,信道长度为64,用于估计信道的导频采用相位均匀分布的复指数符号,基站天线数为64,用户数为5,OFDM子载波个数为512,算法迭代次数为2。从图中可以看出,在没有利用本发明的算法进行相位噪声补偿的情况下,系统性能很差,本发明的算法可以有效抑制相位噪声带来的不利影响。同时,图中也画出了传统变分贝叶斯推断算法的性能曲线,由于本发明的算法属于该算法的近似,该算法的理论性能应优于本发明的算法,但是图中可以看出两条曲线基本重合,本发明的算法的复杂度远远低于传统贝叶斯推断算法,因此本发明的算法更具实用价值。
Claims (1)
1.大规模MIMO系统上行链路的联合信道估计和数据解调方法,设定带有相位噪声的MIMO OFDM系统上行链路中,发射端有K个用户,每个用户有1根天线,接收端基站有M根天线,发射端第k个用户和接收端第m根天线之间的时域信道矢量记为其中L为信道矢量的长度;对于每个OFDM符号,接收端第m根天线的时域信号表达式为
其中,是第m根天线上的时域接收信号,N是OFDM子载波的个数,是接收端第m根天线的相位噪声矩阵,是第k个用户到接收端第m根天线之间的Toeplitz信道矩阵,它的第1列为其中01×(N-L)表示元素全为0、长度为N-L的行矢量;F∈CN×N是归一化的FFT矩阵,它的第i行第j个元素为dk=[dk,1,dk,2,…,dk,N]T是第k个用户发送的数据或者导频序列;是时域的复高斯白噪声序列,
其中Hm,k=diag{[Hm,k,1,Hm,k,2,…,Hm,k,N]T},
将(3)改写为
将(5)进一步改写为
其中1表示长度为N的全1列向量,θm=[θm,1,θm,2,…,θm,N]T为实高斯分布的相位噪声矢量,即θm=N(0,Φ);设定θm的协方差矩阵Φ为实对称矩阵,其特征值是实数,用正交矩阵进行相似对角化:
Φ=UΛUT (8)
其中Λ=diag{[λ1,λ2,…,λN]T}是对角矩阵,对角元素为Φ的降序排列的特征值,U是正交矩阵,它的每一列是Λ对应列的特征值的特征向量;通过计算可以发现,Λ中的对角元只有前若干项的值较大,其它的元素和前若干项相比很小,因此可以只取前I项来近似,即
Φ≈VΓVT (9)
Γ=diag{[λ1,λ2,…,λI]T}是以Λ中前I个特征值为对角元素的对角矩阵,V∈CN×I是由U的前I列组成的矩阵;对相位噪声矢量θm作线性变换
由高斯分布的性质可知,xm=N(0,Γ),由于Γ为对角阵,所以xm的各个分量之间是相互独立的;把(10)代入(7)可得
hm服从复高斯的先验分布
其特征在于,所述方法包括以下步骤:
S1、通过以下步骤的迭代来实现对信道矢量的估计:
S11、假定相位噪声为0,作为迭代的初始值;
S15、对信道矢量的先验协方差矩阵D进行更新;
S16、循环步骤S13—S15,在已知接收信号的条件下,信道矢量将收敛于一个稳定的值;
S2、通过以下步骤的迭代来实现对数据符号的估计:
S21、计算相位噪声展开矢量的后验分布的均值和方差
S22、利用估计的相位噪声对接收信号进行补偿,再进行ZF合并,利用最大似然判决得到数据符号的估计值
S23、循环步骤S21—S22,在已知接收信号的条件下,数据符号矢量将收敛于一个稳定的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810748828.3A CN108965174B (zh) | 2018-07-10 | 2018-07-10 | 大规模mimo系统上行链路的联合信道估计和数据解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810748828.3A CN108965174B (zh) | 2018-07-10 | 2018-07-10 | 大规模mimo系统上行链路的联合信道估计和数据解调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108965174A CN108965174A (zh) | 2018-12-07 |
CN108965174B true CN108965174B (zh) | 2021-06-08 |
Family
ID=64483569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810748828.3A Active CN108965174B (zh) | 2018-07-10 | 2018-07-10 | 大规模mimo系统上行链路的联合信道估计和数据解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108965174B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109639606B (zh) * | 2019-01-22 | 2021-11-26 | 广东工业大学 | 单载波空间复用mimo的本振相位噪声估计方法、装置 |
CN112953864B (zh) * | 2021-04-15 | 2022-10-14 | 东南大学 | 一种保护隐私的无蜂窝混合大规模mimo系统信道估计方法 |
CN113315561B (zh) | 2021-05-25 | 2022-04-08 | 之江实验室 | 一种mimo系统中共参考多通道相位噪声抑制方法 |
CN115766350A (zh) * | 2022-09-28 | 2023-03-07 | 中国传媒大学 | 大规模mimo系统中一种同时信道估计与定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107370693A (zh) * | 2017-08-07 | 2017-11-21 | 电子科技大学 | 大规模mimo系统和dp先验下多用户信道估计方法 |
CN107947839A (zh) * | 2017-11-27 | 2018-04-20 | 电子科技大学 | 用于大规模mimo系统的相位噪声补偿抑制方法 |
CN108111441A (zh) * | 2018-01-12 | 2018-06-01 | 电子科技大学 | 基于变分贝叶斯推断的信道估计方法 |
-
2018
- 2018-07-10 CN CN201810748828.3A patent/CN108965174B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107370693A (zh) * | 2017-08-07 | 2017-11-21 | 电子科技大学 | 大规模mimo系统和dp先验下多用户信道估计方法 |
CN107947839A (zh) * | 2017-11-27 | 2018-04-20 | 电子科技大学 | 用于大规模mimo系统的相位噪声补偿抑制方法 |
CN108111441A (zh) * | 2018-01-12 | 2018-06-01 | 电子科技大学 | 基于变分贝叶斯推断的信道估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108965174A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107947839B (zh) | 用于大规模mimo系统的相位噪声补偿抑制方法 | |
CN108965174B (zh) | 大规模mimo系统上行链路的联合信道估计和数据解调方法 | |
CN108736938B (zh) | 用于大规模mimo上行链路信道估计和数据解调方法 | |
Wang et al. | Bayesian optimal data detector for mmWave OFDM system with low-resolution ADC | |
CN107231216B (zh) | 基于gamp算法的相位噪声补偿抑制方法 | |
EP2453587B1 (en) | Method, apparatus and system for transmitting channel state information | |
CN110166401B (zh) | 大规模mimo正交频分复用系统的相位噪声抑制方法 | |
CN109150260B (zh) | 带有双端相位噪声的大规模mimo系统上行链路数据估计方法 | |
CN108881078B (zh) | 基于变分贝叶斯推断的毫米波系统双端相位噪声抑制方法 | |
CN108924075B (zh) | 基于最大后验准则的毫米波系统双端相位噪声抑制方法 | |
CN109257080B (zh) | 大规模mimo系统下行链路中的多用户相位噪声补偿抑制方法 | |
Han et al. | Eigendecomposition-based partial FFT demodulation for differential OFDM in underwater acoustic communications | |
Domouchtsidis et al. | Constant envelope MIMO-OFDM precoding for low complexity large-scale antenna array systems | |
CN108965172B (zh) | 带有相位噪声的大规模mimo系统上行链路信道估计方法 | |
Kumari et al. | Two choice hard thresholding pursuit (TCHTP) for delay-Doppler channel estimation in OTFS | |
Liu et al. | Research on pilot-based channel estimation algorithms | |
Teng et al. | Joint estimation of channel and I/Q imbalance in massive MIMO: A two-timescale optimization approach | |
Zhang et al. | A Cross-Domain Iterative OTFS Receiver for Sparse Doubly Selective Channels | |
Taner et al. | Alternating projections method for joint precoding and peak-to-average-power ratio reduction | |
CN108965195B (zh) | 大规模mimo系统下行链路中的单用户相位噪声补偿抑制方法 | |
CN108768480B (zh) | 带有相位噪声的大规模mimo系统上行链路数据估计方法 | |
CN108696465A (zh) | 带有双端相位噪声的大规模mimo系统上行链路信道估计方法 | |
Yang et al. | Deep transfer learning-based adaptive beamforming for realistic communication channels | |
Mokhtar et al. | New design of pilot patterns for joint semi-blind estimation of CFO and channel for OFDM systems | |
CN107248901A (zh) | 基于分块和gamp算法融合的相位噪声补偿抑制方法 |
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 |