CN109274623A - 一种基于大点数fft修正载波频偏的实现方法 - Google Patents
一种基于大点数fft修正载波频偏的实现方法 Download PDFInfo
- Publication number
- CN109274623A CN109274623A CN201811286538.8A CN201811286538A CN109274623A CN 109274623 A CN109274623 A CN 109274623A CN 201811286538 A CN201811286538 A CN 201811286538A CN 109274623 A CN109274623 A CN 109274623A
- Authority
- CN
- China
- Prior art keywords
- frequency deviation
- fft
- frequency
- signal
- carrier wave
- 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
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种基于大点数FFT修正载波频偏的实现方法,先将接收到信号经过正交混频,转化为I、Q两路信号,再把I、Q两路的输出信号组合成复信号,利用两次FFT矩阵法估算复信号的频谱特性,然后在频谱特性上找到最大峰值点对应的位置并计算频偏值,最后将频偏值转换为压控振荡器的频率控制字,从而产生与频偏值相同的频率,再通过反馈环路后,经过正交混频补偿I、Q两路信号,并通过低通滤波器消除高频分量,输出消除频偏后的正确基带信号。
Description
技术领域
本发明属于信号处理技术领域,更为具体地讲,涉及一种基于大点数FFT修正载波频偏的实现方法。
背景技术
在无线通信时,如果通信双方中任何一方处于移动状态时,接收到的信号都会产生一个多普勒频移,多普勒频移是产生载波频差的一个重要原因。另外,由于通信双方载波振荡器以及传输过程中有关器件频率的精度所限,在接收信号与本地载波之间总会存在频差。为了实现信息在通信双方的正确传递,我们应该想办法消除掉这些频差。
基于锁相环的闭环载波同步方法,捕获范围很小,所需的捕获时间很长,特别在低信噪比下,这种闭环结构往往不能适应突发方式下快速载波同步,因此许多开环载波提取结构被提出,并越来越多地应用在通信***中。
基于快速傅立叶变换进行载波频偏估计就属于开环载波提取,其思想源于数字信号处理中信号功率谱密度估计的周期图方法。对接收信号中包含周期性剩余频偏、且受到高斯白噪声及其他周期信号干扰的情况,做变换并进行功率谱密度估计,可以在频域上精确地分离出频偏的功率谱,从而能准确地估计频偏。该算法无需数据辅助,当输入信噪比大于一个门限时,估计精度与输入信噪比无关,只由估计参数确定,且非常适合数字编程实现。
由于该方法的频率分辨率由***采样率及FFT点数决定,目前提升频偏估计精度的方法很多都通过抽点来降低采样率。但是,降低采样率就会相应的减小分析带宽,减小频偏估计范围,因此无法分析较大的频偏。对于高阶通信***,如256QAM调制信号***,被测信号在实时分析带宽内的位置通常是不确定的,频偏可能达到几十MHz,因此没有大的估计带宽就不能估计并消除这种频偏。因此增加FFT点数就是一种不损失分析带宽的改善估计精度的较好的方法,
但随着点数的增加,直接计算会消耗大量的硬件资源,当前的硬件是无法满足的,因此需要研究等效替代算法。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于大点数FFT修正载波频偏的实现方法,利用两次FFT矩阵法估计载波频偏,适用于需要大分析带宽和高频率分辨率来精确估计频偏的场合。
为实现上述发明目的,本发明一种基于大点数FFT修正载波频偏的实现方法,其特征在于,包括以下步骤:
(1)、设接收到信号为f0为r(n)的频率;在不考虑定时同步、码间干扰及噪声的影响下,将r(n)与本地压控振荡器NCO产生的正交本振信号进行混频处理,再经低通滤波器滤除高频分量后转化为I、Q两路,分别为:
I(n)=Acos[2πΔfnTs+θ0]
Q(n)=Asin[2πΔfnTs+θ0]
其中,A为I、Q两路信号的幅值,Δf为载波频偏,Ts为***采样周期,θ0为载波初始相位;
(2)、把I、Q两路的输出信号组合成复信号x(n)=Q(n)+j·I(n),利用两次FFT矩阵法对x(n)做N点FFT变换,得到复信号x(n)的频谱特性
(3)、在X(k)中找到最大峰值点对应的位置记为kmax,并根据kmax计算频偏值Δf,如果kmax出现在(0~N/2-1),Δf为正,|Δf|=kmaxFs/N,如果kmax出现在(N/2~N-1),Δf为负,|Δf|=(N-kmax)Fs/N;
(4)、将计算得到的频偏值Δf转换为压控振荡器NCO的频率控制字,从而产生与频偏Δf相同的频率,再通过反馈环路后,经过正交混频补偿I、Q两路信号,并通过低通滤波器消除高频分量,输出消除频偏后的正确基带信号。
本发明的发明目的是这样实现的:
本发明一种基于大点数FFT修正载波频偏的实现方法,先将接收到信号经过正交混频,转化为I、Q两路信号,再把I、Q两路的输出信号组合成复信号,利用两次FFT矩阵法估算复信号的频谱特性,然后在频谱特性上找到最大峰值点对应的位置并计算频偏值,最后将频偏值转换为压控振荡器的频率控制字,从而产生与频偏值相同的频率,再通过反馈环路后,经过正交混频补偿I、Q两路信号,并通过低通滤波器消除高频分量,输出消除频偏后的正确基带信号。
同时,本发明一种基于大点数FFT修正载波频偏的实现方法还具有以下有益效果:
(1)、本发明在不损失***带宽的前提下,通过增加FFT点数,提高了频率分辨率;
(2)、利用两次FFT矩阵法来替代现有方法中的普通点数FFT,打破IP核运算量级的限制,提升频偏估计精度;
(3)、本发明借鉴图像处理中的二维信号处理方法,将一维的大点数FFT算法转换成适用于矩阵应用的两次FFT矩阵方法,在增加运算复杂度的前提下提高运算效率,解决大点数FFT快速运算的问题。
附图说明
图1是一种基于大点数FFT修正载波频偏的实现方法流程图;
图2是传统方法的频偏校正原理图;
图3是两次FFT矩阵法实现原理图;
图4是旋转因子校正示意图;
图5是频偏计算结果示意图;
图6是载波频偏校正前后对比图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是一种基于大点数FFT修正载波频偏的实现方法流程图。
在本实施例中,如图1所示,本发明一种基于大点数FFT修正载波频偏的实现方法,包括以下步骤:
S1、对接收的载波信号进行正交混频转换
以16QAM信号为例,设接收到信号为f0为r(n)的频率;在不考虑定时同步、码间干扰及噪声的影响下,将r(n)与本地压控振荡器NCO产生的正交本振信号进行混频处理,再经低通滤波器滤除高频分量后转化为I、Q两路,分别为:
I(n)=Acos[2πΔfnTs+θ0]
Q(n)=Asin[2πΔfnTs+θ0]
其中,A为I、Q两路信号的幅值,Δf为载波频偏,Ts为***采样周期,θ0为载波初始相位;
S2、估算复信号的频谱特性
把I、Q两路的输出信号组合成复信号x(n)=Q(n)+j·I(n),如图2所示,传统的方法是对复信号x(n)直接做N点FFT变换,得到复信号频谱特性,这样的方法其主要缺点是受FFT点数限制,频偏估计精度低;
在本发明中,如图3所示,主要是利用两次FFT矩阵法对x(n)做N点FFT变换,得到复信号x(n)的频谱特性,下面我们对其进行详细说明:
S2.1、根据点数N,将复信号x(n)分解为N1行N2列的矩阵X,其中,N=N1×N2;
S2.2、对矩阵X按列进行N2次大小为N1的DFT运算,得到N2个N1点长度的序列
在本实施例中,若直接对矩阵X进行DFT运算,则表达式为
设n=n1·N2+n2,n1=0,1,...,N1-1,n2=0,1,...,N2-1,则上式可以转化为:
设k=k2·N1+k1,k1=0,1,...,N1-1,k2=0,1,...,N2-1,上式可以转化为:
化简后,得:
其中,括号[·]内的第一级加权求和可看作对列进行N1个点的DFT变换,由此得出第二步:对矩阵X按列进行N2次大小为N1的DFT运算,得到N2个N1点长度的序列
S2.3、令则简化后的公式变为:
同样地,在得到的公式中外层的第2级加权求和可看作是对第1级FFT运算得到的数据经过旋转因子加权后再进行第2轮的FFT迭代运算,此次FFT迭代运算是按照行数据进行N2点的FFT加权运算。
由此得出第三步,如图4所示,将Y1中数据按行读出,然后分别乘以旋转因子得到矩阵Y2:
S2.4、将矩阵Y2按照行进行N1次大小为N2点的FFT运算,再按列顺序读出,得到复信号x(n)的频谱特性X(k);
S3、计算频偏值
在X(k)中找到最大峰值点对应的位置记为kmax,并根据kmax计算频偏值|Δf|;
频偏估计值Δf与FFT幅度频谱X(k)的最大峰值点kmax值一一对应;载波频偏估计值Δf与采样速率FS成正比例关系,而与采样点数N成反比例关系,由于kmax是整数,并且由此决定了频偏检测值精度的最小分辨率为Fs/N。
在本实施例中,***采样频率FS=1/TS,由于Δf有正负之分,且复信号的离散傅里叶变换是单边频谱。如果kmax出现在(0~N/2-1),Δf为正,|Δf|=kmaxFs/N,如果kmax出现在(N/2~N-1),Δf为负,|Δf|=(N-kmax)Fs/N;正是有了这个特性,频偏Δf的绝对值大小和正负都能被估计出来。
S4、根据频偏值修正载波频偏
将计算得到的频偏值Δf转换为压控振荡器NCO的频率控制字,从而产生与频偏Δf相同的频率,再通过反馈环路后,经过正交混频补偿I、Q两路信号,并通过低通滤波器消除高频分量,输出消除频偏后的正确基带信号。
实例
在本实施例中,以256QAM信号为例,设接收的已解调的基带信号载波频偏为Δf=40MHz,设采样率为Fs=100MHz,符号速率为Rs=5MSps,FFT点数为N=1M。
将正交混频后的I、Q两路信号合成复信号并做两次矩阵FFT分析,得到频谱特性。
其中,两次矩阵FFT步骤如下:
1)、将采样到的1M点,按行顺序排列,得到1k×1k的采样点矩阵,再按列顺序读出。这里主要是对RAM IP核操作,用两个512kRAM来存储1M点。实现行顺序写入,列顺序读出。
2)、对1k×1k采样点矩阵,对读出的每一列数据做1k点的FFT,得到1M点的结果矩阵,并按列顺序存储,再将得到的1M点的结果矩阵按行顺序读出。由于FFT结果是复数,分为实部和虚部,因此需要4个512K的RAM来存储。其中FFT主要利用FPGA中FFT的IP核完成。
3)、将步骤2)中按行顺序读出的数据乘以旋转因子,然后按行顺序对上述乘法结果做1k点的FFT,得到1M点的FFT结果。这里的旋转因子精度利用Matlab已做了验证,这里将旋转因子乘以100后取整。用Matlab计算旋转因子,产生.coe文件,存储在ROM中。
4)、将步骤3)中的FFT结果,按行顺序存入RAM,再按列顺序输出为1M点的向量,此向量即为1M点的FFT频谱并得到峰值点位置。对4块RAM读写操作,读写操作顺序与步骤1)、2)步中的读写操作一致。
最后根据峰值点位置得到频偏值并转化为NCO的频率控制字,并控制NCO生成与频偏频率相同的频率信号并经过正交混频补偿给I、Q两路信号,并通过低通滤波器滤除高频分量,得到正确的I、Q信号。
在本实施例中,计算频偏的结果如图5所示,频偏消除前后的星座图如图6所示,其中,左边为频偏校正前,右边为频偏校正后,通过左右两幅图对比可知,改频偏校正方法有效的去除了大频偏并恢复了星座图,因此,证明了该方法的有效性,可以很好的完成大频偏的高精度校正。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种基于大点数FFT修正载波频偏的实现方法,其特征在于,包括以下步骤:
(1)、设接收到信号为f0为r(n)的频率;在不考虑定时同步、码间干扰及噪声的影响下,将r(n)与本地压控振荡器NCO产生的正交本振信号进行混频处理,再经低通滤波器滤除高频分量后转化为I、Q两路,分别为:
I(n)=Acos[2πΔfnTs+θ0]
Q(n)=Asin[2πΔfnTs+θ0]
其中,A为I、Q两路信号的幅值,Δf为载波频偏,Ts为***采样周期,θ0为载波初始相位;
(2)、把I、Q两路的输出信号组合成复信号x(n)=Q(n)+j·I(n),利用两次FFT矩阵法对x(n)做N点FF变换,得到复信号x(n)的频谱特性
(3)、在X(k)中找到最大峰值点对应的位置记为kmax,并根据kmax计算频偏值Δf,如果kmax出现在(0~N/2-1),Δf为正,|Δf|=kmaxFs/N,如果kmax出现在(N/2~N-1),Δf为负,|Δf|=(N-kmax)Fs/N;
(4)、将计算得到的频偏值Δf转换为压控振荡器NCO的频率控制字,从而产生与频偏Δf相同的频率,再通过反馈环路后,经过正交混频补偿I、Q两路信号,并通过低通滤波器消除高频分量,输出消除消除频偏后的正确基带信号。
2.根据权利要求1所述的基于大点数FFT修正载波频偏的实现方法,其特征在于,所述步骤(2)中,利用两次FFT矩阵法对x(n)做N点FFT变换的方法为:
(2.1)、根据点数N,将复信号x(n)分解为N1行N2列的矩阵X,其中,N=N1×N2;
(2.2)、对矩阵X按列进行N2次大小为N1的DFT运算,得到N2个N1点长度的序列
(2.3)、令将Y1中数据按行读出,然后分别乘以旋转因子得到矩阵Y2:
(2.4)、将矩阵Y2按照进行N1次大小为N2点的FFT运算,再按列顺序读出,得到复信号x(n)的频谱特性X(k);
其中,n=n1·N2+n2,k=k2·N1+k1,n1=0,1,...,N1-1,n2=0,1,...,N2-1,k1=0,1,...,N1-1,k2=0,1,...,N2-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811286538.8A CN109274623B (zh) | 2018-10-31 | 2018-10-31 | 一种基于大点数fft修正载波频偏的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811286538.8A CN109274623B (zh) | 2018-10-31 | 2018-10-31 | 一种基于大点数fft修正载波频偏的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109274623A true CN109274623A (zh) | 2019-01-25 |
CN109274623B CN109274623B (zh) | 2021-07-13 |
Family
ID=65192127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811286538.8A Active CN109274623B (zh) | 2018-10-31 | 2018-10-31 | 一种基于大点数fft修正载波频偏的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109274623B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110518908A (zh) * | 2019-08-01 | 2019-11-29 | 电子科技大学 | 基于施密特正交化来消除iq两路信号相位误差的方法 |
CN110674456A (zh) * | 2019-09-26 | 2020-01-10 | 电子科技大学 | 一种信号采集***的时频转换方法 |
CN111458563A (zh) * | 2020-03-05 | 2020-07-28 | 熊军 | 一种双谱线幅度测量信号频偏的方法及装置 |
CN111597498A (zh) * | 2020-05-18 | 2020-08-28 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
CN112751794A (zh) * | 2020-12-22 | 2021-05-04 | 湖南艾科诺维科技有限公司 | 基于fpga的非连续相位2fsk信号的符号速率估计方法 |
CN114814815A (zh) * | 2022-03-16 | 2022-07-29 | 华域汽车***股份有限公司 | 一种基于插值和单点dft滤波的解决信号处理跨越损失的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1472900A (zh) * | 2002-08-01 | 2004-02-04 | 上海交通大学 | 正交频分复用通信***载波频率偏移估计方法 |
US20110142170A1 (en) * | 2007-07-27 | 2011-06-16 | Mediatek Inc. | Signal Receiver and Method for Estimating Residual Doppler Frequencies of Signals Thereof |
EP2736208A2 (en) * | 2012-11-27 | 2014-05-28 | Gilat Satellite Networks, Ltd. | Frequency offset estimation |
CN105490980A (zh) * | 2015-12-11 | 2016-04-13 | 航天恒星科技有限公司 | 一种载波频偏估计方法和*** |
CN106341358A (zh) * | 2016-09-20 | 2017-01-18 | 湖北三江航天红林探控有限公司 | 高动态无人机数据链的同步装置及方法 |
CN108092929A (zh) * | 2017-12-27 | 2018-05-29 | 北京理工大学 | 一种用于太赫兹通信的同步方法 |
-
2018
- 2018-10-31 CN CN201811286538.8A patent/CN109274623B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1472900A (zh) * | 2002-08-01 | 2004-02-04 | 上海交通大学 | 正交频分复用通信***载波频率偏移估计方法 |
US20110142170A1 (en) * | 2007-07-27 | 2011-06-16 | Mediatek Inc. | Signal Receiver and Method for Estimating Residual Doppler Frequencies of Signals Thereof |
EP2736208A2 (en) * | 2012-11-27 | 2014-05-28 | Gilat Satellite Networks, Ltd. | Frequency offset estimation |
CN105490980A (zh) * | 2015-12-11 | 2016-04-13 | 航天恒星科技有限公司 | 一种载波频偏估计方法和*** |
CN106341358A (zh) * | 2016-09-20 | 2017-01-18 | 湖北三江航天红林探控有限公司 | 高动态无人机数据链的同步装置及方法 |
CN108092929A (zh) * | 2017-12-27 | 2018-05-29 | 北京理工大学 | 一种用于太赫兹通信的同步方法 |
Non-Patent Citations (1)
Title |
---|
MENG JIE等: "Research on the parallel FFT based on digital down conversion and phase difference correction techniques", 《2017 IEEE 13TH INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT & INSTRUMENTS》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110518908A (zh) * | 2019-08-01 | 2019-11-29 | 电子科技大学 | 基于施密特正交化来消除iq两路信号相位误差的方法 |
CN110518908B (zh) * | 2019-08-01 | 2022-06-14 | 电子科技大学 | 基于施密特正交化来消除iq两路信号相位误差的方法 |
CN110674456A (zh) * | 2019-09-26 | 2020-01-10 | 电子科技大学 | 一种信号采集***的时频转换方法 |
CN110674456B (zh) * | 2019-09-26 | 2022-11-22 | 电子科技大学 | 一种信号采集***的时频转换方法 |
CN111458563A (zh) * | 2020-03-05 | 2020-07-28 | 熊军 | 一种双谱线幅度测量信号频偏的方法及装置 |
CN111597498A (zh) * | 2020-05-18 | 2020-08-28 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
CN111597498B (zh) * | 2020-05-18 | 2023-04-18 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
CN112751794A (zh) * | 2020-12-22 | 2021-05-04 | 湖南艾科诺维科技有限公司 | 基于fpga的非连续相位2fsk信号的符号速率估计方法 |
CN112751794B (zh) * | 2020-12-22 | 2022-06-07 | 湖南艾科诺维科技有限公司 | 基于fpga的非连续相位2fsk信号的符号速率估计方法 |
CN114814815A (zh) * | 2022-03-16 | 2022-07-29 | 华域汽车***股份有限公司 | 一种基于插值和单点dft滤波的解决信号处理跨越损失的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109274623B (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109274623A (zh) | 一种基于大点数fft修正载波频偏的实现方法 | |
CN1750531B (zh) | 失真补偿正交调制器和无线发送器 | |
EP3264622B1 (en) | Low complexity, low power and long range radio receiver | |
CN108092929A (zh) | 一种用于太赫兹通信的同步方法 | |
CN100553246C (zh) | 单通道检测输入信号相位差与相对幅度的方法 | |
CN101611554A (zh) | 相位噪声校正设备及其方法 | |
CN110300079B (zh) | 一种msk信号相干解调方法及*** | |
CN106302296A (zh) | 高动态窄带信号频率跟踪方法 | |
CN108243138B (zh) | 一种适用于水声通信***的组合多普勒估计方法 | |
WO2015070820A1 (zh) | 扩频信号的生成方法、生成装置、接收方法和接收装置 | |
CN107247254A (zh) | 一种非线性调频信号去斜处理方法 | |
CN103748846B (zh) | 数字通信***中的频率偏差估计方法和模块 | |
CN102332935B (zh) | 载波补偿***和方法 | |
CN115956356A (zh) | 对调制信号进行解调 | |
CN107370705A (zh) | 高动态弱连续相位调制信号捕获中fft的优化方法 | |
CN105049174A (zh) | 一种用于oqpsk调制的载波与时钟联合同步方法 | |
CN102655491A (zh) | 相干解调频移键控调制信号的频率偏移估计方法及*** | |
CN103873416B (zh) | 一种evm相位估计与补偿方法 | |
CN103023829B (zh) | 一种调制精度的估计方法及装置 | |
CN106789791A (zh) | 基于共轭对称训练序列的移动通信***载波频偏估计方法 | |
CN104901907B (zh) | 一种动态环境下基于数据辅助的稳健信噪比估计方法 | |
CN106101039A (zh) | 一种基于数据辅助精度可调的频偏估计方法 | |
US7751518B2 (en) | System and method for executing preamble detection, symbol timing recovery, and frequency offset estimation | |
US6470056B1 (en) | Offset QPSK modulation analytic system | |
CN105607091A (zh) | 一种改进的基于扩展卡尔曼滤波的载波跟踪环 |
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 |