CN107438038A - 一种fbmc/oqam的导频设计和同步信道估计方法 - Google Patents

一种fbmc/oqam的导频设计和同步信道估计方法 Download PDF

Info

Publication number
CN107438038A
CN107438038A CN201710424321.8A CN201710424321A CN107438038A CN 107438038 A CN107438038 A CN 107438038A CN 201710424321 A CN201710424321 A CN 201710424321A CN 107438038 A CN107438038 A CN 107438038A
Authority
CN
China
Prior art keywords
pilot tone
row
frequency
pilot
data
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
CN201710424321.8A
Other languages
English (en)
Other versions
CN107438038B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201710424321.8A priority Critical patent/CN107438038B/zh
Publication of CN107438038A publication Critical patent/CN107438038A/zh
Application granted granted Critical
Publication of CN107438038B publication Critical patent/CN107438038B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/0202—Channel estimation
    • H04L25/0224—Channel estimation using sounding signals
    • H04L25/0228—Channel estimation using sounding signals with direct estimation from sounding signals
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • 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/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • 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/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2602—Signal structure
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2626—Arrangements specific to the transmitter only
    • H04L27/2627—Modulators
    • H04L27/264—Pulse-shaped multi-carrier, i.e. not using rectangular window
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2668—Details of algorithms
    • H04L27/2669—Details of algorithms characterised by the domain of operation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2668—Details of algorithms
    • H04L27/2673—Details of algorithms characterised by synchronisation parameters
    • H04L27/2675—Pilot or known symbols
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2697—Multicarrier modulation systems in combination with other modulation techniques
    • H04L27/2698—Multicarrier modulation systems in combination with other modulation techniques double density OFDM/OQAM system, e.g. OFDM/OQAM-IOTA system
    • 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
    • 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)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

一种FBMC/OQAM的导频设计和同步信道估计方法，在FBMC/OQAM***的发射端***导频；在***接收端，利用与导频列相邻的数据列的反转对称性，对接收信号做滑动相关以确定起始位置；然后在频域利用同一子载波上导频的对应关系，得到频偏估计值，通过迭代减小与导频列相邻的数据列和信道噪声对估计的干扰；频偏补偿后首先假设与导频列相邻的数据列是零保护列，由此计算出伪导频并在频域进行信道估计，然后解调出与导频列相邻的数据符号；最后利用解调出的与导频列相邻的数据符号重新计算伪导频来估计信道。本发明在FBMC/OQAM***中联合估计了时偏、频偏和信道，不但提高了频谱效率且估计性能也有显著提升。

Description

一种FBMC/OQAM的导频设计和同步信道估计方法
技术领域
本发明涉及无线通信领域，具体涉及一种FBMC/OQAM的导频设计和同步信道估计方法。
背景技术
由于多载波调制***不仅能够提供高速的数据传输，成倍的提升***容量，同时能够有效地对抗信道的频率选择性衰落，因此受到了学术界和工业界的广泛关注和研究。然而，带外泄漏过高一直是多载波调制***的一大弊病。偏置正交幅度调制的滤波器组多载波(FBMC/OQAM，Filter Bank Multicarrier with Offset Quadrature AmplitudeModulation)***通过采用具有良好时频聚焦特性的原型滤波器，可以有效降低信号的带外泄漏。此外， FBMC/OQAM(OFDM/OQAM)***通过多相滤波器和快速傅里叶变换等操作，大大降低了自身的复杂度和运算量，应用前景广阔。
基于数据辅助的时偏、频偏和信道估计一直以来都是FBMC/OQAM***的研究热点。目前，已有大量文献分开讨论了时频同步和信道估计问题。Mattera D,Tanda M在SignalProcessing,2012,92(9):2284–2292上发表的“Data-aided synchronization for OFDM/OQAM systems[J]”提出利用FBMC/OQAM的共轭对称性来进行符号定时，并利用重复导频序列的相互关系计算频偏，该方法需要添加很多同步符号和保护符号才能抵消FBMC/OQAM***的重叠性，频谱效率较低。现在受到广泛认可的信道估计方法是J Du，S Signell在IEEEInternational Conference on Communications,2009:1-6上发表的“Novel preamble-based channel estimation for OFDM/OQAM systems[C]”中提出的干扰近似法(IAM)。同时有一些文献联合估计了频偏和信道，G.Cheng,Y.Xiao,S.Li在IEICE Transaction onCommunications,2012, E95-B(5):1848-1851上发表的“Joint Frequency Offset andChannel Estimation for OFDM/OQAM system[J]”利用同一子载波上的导频在频域的对应关系计算频偏，频偏补偿后基于干扰消除法(ICM)估计信道。然而同时估计时偏、频偏和信道参数时的估计性能和频谱效率仍有待提高。
发明内容
本发明的目的是提供一种FBMC/OQAM的导频设计和同步信道估计方法，在 FBMC/OQAM***中联合估计了时偏、频偏和信道，不但提高了频谱效率且估计性能也有显著提升。
为了实现以上目的，本发明所采用的技术方案如下：
一种FBMC/OQAM的导频设计和同步信道估计方法，在FBMC/OQAM***的发射端***导频：频率时间坐标上，一列幅度减半的数据与导频列交叉重复排列，并在最后再添加该列幅度减半的数据以间隔导频与正常幅度的数据，幅度减半的数据与导频共五列；其中导频按照[1,j,-1,-j]T的规律不断重复；在FBMC/OQAM***的接收端，首先利用与导频列相邻的数据列的反转对称性，对接收信号做滑动相关以确定起始位置；然后在频域利用同一子载波上导频的对应关系，得到频偏估计值，通过迭代减小与导频列相邻的数据列和信道噪声对估计的干扰；频偏补偿后首先假设与导频列相邻的数据列是零保护列，由此计算出伪导频并在频域进行信道估计，然后解调出与导频列相邻的数据符号；最后利用解调出的与导频列相邻的数据符号计算伪导频来估计信道，完成时偏、频偏和信道的联合估计。
本发明进一步的改进在于，具体包括以下步骤：
(1)在发送端对信息比特流进行QPSK调制和串并转换，然后进行OQAM预处理，将复数数据符号转换成实数数据符号；
(2)在发送端数据符号块的头部***导频：在频率时间坐标上，每一列对应一个时间的 M个并行子载波，导频***的要求如下：导频放到数据块的第二列和第四列，每列按照[1,j,-1,-j]T的规律不断循环重复；第一、三、五列为幅度减半的相同的待发送实数数据列，记为：a＝[a0,a1,a2,…,aM/2,…,aM-2,aM-1]T，其中aM-i＝(-1)iai,i＝1,2,…,M/2-1，所以数据列a进一步表示为a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T；第五列之后均为正常幅度的待发送的实数数据列；
(3)对***导频的数据块进行调制，得到发送端待发送的调制信号s[k]为：
其中，am,n为频率时间坐标(m,n)上的数据符号或者导频，是给待发送实数数据符号添加的初始相位；
(4)经过多径信道后，接收端收到的信号与发送端发送的信号之间存在频率偏移，接收信号r[k]表示为
其中，h[k]代表的是频率选择性信道，其长度为Δ+1；n[k]表示均值为0，方差为σ2的加性高斯白噪声；φ表示以子载波间隔归一化的载波频偏，φ∈(-0.5,0.5]；
(5)首先进行接收信号r[k]的符号定时，判断起始位置，确定FBMC/OQAM的接收信号的起始位置为
(6)确定接收信号的起始位置后，进行频偏估计，得到频偏估计值：
(7)根据频偏估计值进行频偏补偿，计算出伪导频并在频域进行信道估计；
(8)在频域进行信道均衡，解调出全部的实数数据符号，进行OQAM后置处理，从实数符号变回到复数符号，进行并串转换和QPSK解调，最终恢复出发送的比特数据流。
本发明进一步的改进在于，步骤(5)的具体过程包括如下步骤：
(5.1)由于多相滤波器组的引入，数据符号实际操作中添加的初始相位为jm+k；因此，数据列a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T添加初始相位、经过IFFT得到的序列 b＝[b0b1bM/2b2]T满足反转对称特性，序列b1的顺序对应元素与序列b2的逆序对应元素相等，表示为即b1,t＝b2,M/2-2-t，t＝0,1,…,M/2-2；依照最小二乘原则，找到最满足上述反转对称性的地方，也就是b1的起始位置：
其中表示b1(τ1)的共轭；
(5.2)将接收信号r[k]代入(5.1)的式中，则找到接收信号中最满足上述反转对称性的地方；根据反转对称性定义基于最小二乘准则的代价函数Φ[τ]
对代价函数Φ(τ)进行相应滤波以减小噪声的干扰：
对接收信号进行一维搜索，通过求解下面的最大化问题：
由此确定FBMC/OQAM的接收信号的起始位置为
本发明进一步的改进在于，步骤(6)的具体包括如下步骤：
(6.1)确定接收信号的起始位置后，接收端经过串并转换、多相滤波、FFT、删除初始相位，得到频时格点(m0,n0)处的数据符号表示为
其中是原型滤波器p[k]的自模糊函数；是信道噪声在频时格点(m0,n0)处的噪声；
(6.2)考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，认为超出该范围的频时格点上的数据对目标频时格点上的数据造成的影响忽略不计；两列导频数据满足假设与导频列相邻的数据列都为零保护列，且忽略噪声，则导频符号表示为：
其中
(6.3)将与相乘，其中表示的共轭，由于因此得到：
所以CFO的估计值由下式得到：
(6.4)根据频偏估计值对接收信号r[k]进行频偏补偿，得到频偏补偿后的接收信号重新经过串并转换、多相滤波、FFT、删除初始相位得到重复计算得到一次迭代之后的频偏估计值。
本发明进一步的改进在于，步骤(6)还包括：
(6.5)重复(6.4)N-1次，得到N次迭代之后的频偏估计值，N为整数，且N≥1；通过迭代，降低与导频列相邻的数据列和信道噪声对估计结果的影响。
本发明进一步的改进在于，步骤(7)的具体包括如下步骤：
(7.1)经过串并转换、多相滤波、FFT、删除初始相位得到经过频偏补偿后的频时格点 (m0,n0)处的数据符号假设频偏完全补偿，则表示为：
(7.2)定义为信道在子载波m上的频率响应，多径信道最大时延Δ＜＜M时，每一个子载波上的信道频域响应与邻近子载波上的信道频域响应近似相等；同时，对于时频域聚焦性良好的原型滤波器而言，在多径时延[0,Δ]范围内，Ap(-τ,0)≈Ap(0,0)；考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，经过频偏补偿后的频时格点(m0,n0)处的数据符号进一步表示为：
其中表示所受到的内在干扰；
记并称为伪导频，信道估计的方法如下：
(7.3)首先假设与导频列相邻的数据列都为零保护列，进行信道参数的初步估计，对两列导频分别得到的信道估计结果求平均值，得到信道参数的初步估计值
本发明进一步的改进在于，步骤(7)还包括：
(7.4)步骤(7)还包括以下步骤：(7.4)利用对与导频列相邻的数据符号进行信道均衡，然后减去已知导频的干扰，对得到的第一、三、五列数据符号求平均值以减小噪声的干扰，并进行取实操作、判决，初步恢复出与导频列相邻的数据符号；
(7.5)对于位于第二列的导频，除了利用步骤(7.4)中解调出的与导频列相邻的数据符号外，还利用位于第四列导频的已知信息，重新计算伪导频进行信道估计得到 作为信道参数的最终估计值。
本发明与前述现有时频同步以及信道估计方法相比，其有益效果体现在：本发明通过在 FBMC/OQAM***的发射端***特定导频：频率时间坐标上，一列幅度减半的数据与导频列交叉重复排列，并在最后再添加该列幅度减半的数据以间隔导频与正常幅度的数据，幅度减半的数据与导频共五列；其中导频按照[1,j,-1,-j]T的规律不断重复。在***接收端，首先利用与导频列相邻的数据列的反转对称性，对接收信号做滑动相关以确定起始位置；然后在频域利用同一子载波上导频的对应关系，得到频偏估计值，通过迭代减小与导频列相邻的数据列和信道噪声对估计的干扰；频偏补偿后首先假设与导频列相邻的数据列是零保护列，由此计算出伪导频并在频域进行信道估计，然后解调出与导频列相邻的数据符号；最后利用解调出的与导频列相邻的数据符号重新计算伪导频来估计信道；完成时偏、频偏和信道的联合估计。本发明在FBMC/OQAM***中联合估计了时偏、频偏和信道，不但提高了频谱效率且估计性能也有显著提升。本发明首次将干扰利用思想运用于频偏估计，使得在相同导频长度的情况下，频偏估计性能优于同类估计方法。并用同一段非常短的导频序列同时完成了符号定时、频偏估计和信道估计三项工作，频谱效率大大提高。
附图说明
图1为本发明发射机端导频结构示意图。
图2为本发明时频同步和信道估计方法流程图。
图3为符号定时的代价函数的在AWGN信道、Eb/N0＝20dB情况下的一次实现过程。
图4为符号定时的均方误差曲线。
图5为频偏估计的均方误差曲线。
图6为信道估计的均方误差曲线。
图7为误比特率性能曲线。
具体实施方式：
下面结合附图对本发明做进一步详细说明。
一种FBMC/OQAM***中的时频同步和信道估计方法流程如图2所示。
1、发射端
(1)在发射端对信息比特流进行QPSK调制和串并转换，然后进行OQAM预处理，将复数数据符号转换成实数数据符号。
(2)在发送端数据符号块的头部***如图1所示的导频。在频率时间坐标上，每一列对应一个时间的M个并行子载波，导频***的要求如下：导频放到数据块的第二列和第四列，每列按照[1,j,-1,-j]T的规律不断循环重复；第一、三、五列为幅度减半的相同的待发送实数数据列，记为：a＝[a0,a1,a2,…,aM/2,…,aM-2,aM-1]T，其中aM-i＝(-1)iai,i＝1,2,…,M/2-1，所以数据列a可进一步表示为a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T；第五列之后均为正常幅度的待发送的实数数据列。
(3)对***导频的数据块进行调制，得到发送端待发送的调制信号为：
其中，am,n为频率时间坐标(m,n)上的数据符号或者导频，是给待发送实数数据符号添加的初始相位。
2、接收端
(1)经过多径信道后，接收端收到的信号与发送端发送的信号之间存在频率偏移，接收信号r[k]表示为
其中，h[k]代表的是频率选择性信道，其长度为Δ+1。n[k]表示均值为0，方差为σ2的加性高斯白噪声。φ表示以子载波间隔归一化的载波频偏，φ∈(-0.5,0.5]。
(2)首先进行接收信号r[k]的符号定时，判断起始位置。由于多相滤波器组的引入，数据符号实际操作中添加的初始相位为jm+k。因此，数据列a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T添加初始相位、经过IFFT得到的序列b＝[b0b1bM/2b2]T满足反转对称特性，序列b1的顺序对应元素与序列b2的逆序对应元素相等，表示为即b1,t＝b2,M/2-2-t，t＝0,1,…,M/2-2。依照最小二乘(LS)原则，可以找到最满足上述反转对称性的地方，也就是b1的起始位置：
其中表示b1(τ1)的共轭。
将接收信号r[k]代入上式中，则可以找到接收信号中最满足上述反转对称性的地方。根据反转对称性定义基于最小二乘准则的代价函数
图3给出了Φ(τ)在AWGN信道、Eb/N0＝20dB情况下的一次实现过程，由于 FBMC/OQAM***的重叠性和a的重复性，满足上述反转对称性的地方有四个，所以图3中Φ(τ)出现四个峰值，因此对Φ(τ)进行相应滤波以减小噪声的干扰：
对接收信号进行一维搜索，通过求解下面的最大化问题：
由此确定FBMC/OQAM的符号起始位置为
(3)如图2所示，确定符号起始位置后，接收端先进行频偏估计。接收端经过串并转换、多相滤波、FFT、删除初始相位，得到频时格点(m0,n0)处的数据符号表示为
其中是原型滤波器p[k]的自模糊函数。是信道噪声在频时格点(m0,n0)处的噪声。
考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，认为超出该范围的频时格点上的数据对目标频时格点上的数据造成的影响可以忽略不计。假设与导频列相邻的数据列都为零保护列，且忽略噪声，考虑到两列导频数据满足 则导频符号可以表示为：
其中
(4)将与相乘，其中表示的共轭，由于因此得到：
所以CFO的估计值可以由下式得到：
根据频偏估计值对接收信号r[k]进行频偏补偿，得到频偏补偿后的接收信号重新经过串并转换、多相滤波、FFT、删除初始相位得到重复计算得到一次迭代之后的频偏估计值。
重复频偏估计值对接收信号r[k]进行频偏补偿N-1次，得到N次迭代之后的频偏估计值，其中，N为整数，且N≥1。通过迭代，可以降低与导频列相邻的数据列和信道噪声对估计结果的影响。
(5)根据N次迭代之后的频偏估计值进行频偏补偿，重新经过串并转换、多相滤波、FFT、删除初始相位得到经过频偏补偿后的频时格点(m0,n0)处的数据符号假设频偏已经完全补偿，则可表示为：
定义为信道在子载波m上的频率响应，多径信道最大时延Δ＜＜M时，每一个子载波上的信道频域响应与邻近子载波上的信道频域响应近似相等。同时，对于时频域聚焦性良好的原型滤波器而言，在多径时延[0,Δ]范围内，其模糊函数Ap(-τ,0)几乎不变，即Ap(-τ,0)≈Ap(0,0)。考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，经过频偏补偿后的频时格点(m0,n0)处的数据符号进一步表示为：
其中表示所受到的内在干扰。
记并称为伪导频。信道估计的方法如下：
(6)首先假设与导频列相邻的数据列都为零保护列，进行信道参数的初步估计，对两列导频分别得到的信道估计结果求平均值，得到信道参数的初步估计值然后利用信道参数的初步估计值对与导频列相邻的数据符号进行信道均衡，然后减去已知导频的干扰，对得到的第一、三、五列数据符号求平均值以减小噪声的干扰，并进行取实操作、判决，初步恢复出与导频列相邻的数据符号。
(7)对于位于第二列的导频，除了利用步骤(6)中解调出的与导频列相邻的数据符号外，还可利用位于第四列导频的已知信息，即考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,2)and(α,β)≠(0,0)}的频时格点。重新计算伪导频进行信道估计得到作为信道参数的最终估计值。
(8)利用进行信道均衡解调出全部的实数数据符号，进行OQAM后置处理，从实数符号变回到复数符号。进行并串转换和QPSK解调，最终恢复出发送的比特数据流。
最后算法结束。
本实施例中使用的FBMC/OQAM***采用QPSK星座调制，***带宽 B＝1/Ts＝11.2MHz，子载波个数为M＝512，原型滤波器为PHYDYAS，重叠因子K＝4，多径衰落信道模型采用ITU Vehicular A信道模型，信道时延为[00.31 0.71 1.09 1.73 2.51](单位为μs)，路径增益为[0 -1 -9 -10 -15 -20](单位为dB)；ITU Vehicular B信道模型，信道时延为[0 0.3 8.9 12.9 17.1 20](单位为μs)，路径增益为[-2.5 0 -12.8 -10 -25.2 -16](单位为dB)。
图4给出了利用本发明方法进行符号定时均方根误差曲线，并与Mattera D,TandaM在 Signal Processing,2012,92(9):2284–2292上发表的“Data-aided synchronizationfor OFDM/OQAM systems[J]”符号定时方法(图中记为DMandMT-2012)进行了比较。图4中横坐标表示比特信噪比Eb/N0，单位分贝(dB)，范围0－30dB。纵坐标表示对应的估计均方根误差。“DMandMT-2012”符号定时方法中每次相关运算利用的信号长度为M/2，而本发明符号定时方法每次相关运算利用的符号长度为M，利用的信息量增多，估计性能得到了提升。
图5给出了利用本发明方法进行频偏估计的均方根误差曲线，并与G.Cheng,Y.Xiao,S.Li 在IEICE Transaction on Communications,2012,E95-B(5):1848-1851上发表的“Joint Frequency Offset and Channel Estimation for OFDM/OQAM system[J]”方法(图中记为Joint-2012)进行了比较。图6给出了利用本发明方法进行信道估计的均方误差曲线，并与“Joint-2012”方法进行了比较。图中纵坐标表示信道估计的均方误差NMSE，单位分贝(dB)，定义“Joint-2012”方法进行频偏和信道估计时，规避消除目标时频格点邻域Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}对目标时频格点的干扰，而本发明方法利用了目标时频格点邻域Ωα,β＝{(α,β)＝(±1,0)}对目标时频格点的干扰，同时通过迭代减小了目标时频格点邻域Ωα,β＝{|(α,β)|≤(1,1)andβ≠0}对目标时频格点的干扰，相比于“Joint-2012”方法估计性能和频谱效率都显著提高。
图7给出了利用本发明方法完成时频同步和信道估计后的误比特率曲线，并与“Joint-2012”方法和理想同步、估计情况(图中记为理想情况)下的结果进行了比较。均衡方式均采用迫零(ZF)均衡。由图可以看出，利用本发明方法后***的误比特率相比于“Joint-2012”方法更接近于理想情况。
本发明方法导频长度小于“Joint-2012”方法的导频长度，频偏估计与信道估计性能都优于“Joint-2012”方法，且本发明方法同时还兼顾了符号定时，符号定时的性能也优于“DMandMT-2012”方法。本发明在FBMC/OQAM***中联合估计了时偏、频偏和信道，在提升估计性能的同时提高了频谱效率。
本发明重点研究导频的设计和时偏、频偏、信道的联合估计问题，并在保证估计性能的前提下，提高频谱效率。

Claims (7)

1.一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，在FBMC/OQAM***的发射端***导频：频率时间坐标上，一列幅度减半的数据与导频列交叉重复排列，并在最后再添加该列幅度减半的数据以间隔导频与正常幅度的数据，幅度减半的数据与导频共五列；其中导频按照[1,j,-1,-j]T的规律不断重复；在FBMC/OQAM***的接收端，首先利用与导频列相邻的数据列的反转对称性，对接收信号做滑动相关以确定起始位置；然后在频域利用同一子载波上导频的对应关系，得到频偏估计值，通过迭代减小与导频列相邻的数据列和信道噪声对估计的干扰；频偏补偿后首先假设与导频列相邻的数据列是零保护列，由此计算出伪导频并在频域进行信道估计，然后解调出与导频列相邻的数据符号；最后利用解调出的与导频列相邻的数据符号计算伪导频来估计信道，完成时偏、频偏和信道的联合估计。
2.根据权利要求1所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，具体包括以下步骤：
(1)在发送端对信息比特流进行QPSK调制和串并转换，然后进行OQAM预处理，将复数数据符号转换成实数数据符号；
(2)在发送端数据符号块的头部***导频：在频率时间坐标上，每一列对应一个时间的M个并行子载波，导频***的要求如下：导频放到数据块的第二列和第四列，每列按照[1,j,-1,-j]T的规律不断循环重复；第一、三、五列为幅度减半的相同的待发送实数数据列，记为：a＝[a0,a1,a2,…,aM/2,…,aM-2,aM-1]T，其中aM-i＝(-1)iai,i＝1,2,…,M/2-1，所以数据列a进一步表示为a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T；第五列之后均为正常幅度的待发送的实数数据列；
(3)对***导频的数据块进行调制，得到发送端待发送的调制信号s[k]为：
其中，am,n为频率时间坐标(m,n)上的数据符号或者导频，是给待发送实数数据符号添加的初始相位；
(4)经过多径信道后，接收端收到的信号与发送端发送的信号之间存在频率偏移，接收信号r[k]表示为
其中，h[k]代表的是频率选择性信道，其长度为Δ+1；n[k]表示均值为0，方差为σ2的加性高斯白噪声；φ表示以子载波间隔归一化的载波频偏，φ∈(-0.5,0.5]；
(5)首先进行接收信号r[k]的符号定时，判断起始位置，确定FBMC/OQAM的接收信号的起始位置为
(6)确定接收信号的起始位置后，进行频偏估计，得到频偏估计值：
(7)根据频偏估计值进行频偏补偿，计算出伪导频并在频域进行信道估计；
(8)在频域进行信道均衡，解调出全部的实数数据符号，进行OQAM后置处理，从实数符号变回到复数符号，进行并串转换和QPSK解调，最终恢复出发送的比特数据流。
3.根据权利要求2所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，步骤(5)的具体过程包括如下步骤：
(5.1)由于多相滤波器组的引入，数据符号实际操作中添加的初始相位为jm+k；因此，数据列a＝[a0,a1,a2,…,aM/2,…,a2,-a1]T添加初始相位、经过IFFT得到的序列b＝[b0 b1 bM/2b2]T满足反转对称特性，序列b1的顺序对应元素与序列b2的逆序对应元素相等，表示为即b1,t＝b2,M/2-2-t，t＝0,1,…,M/2-2；依照最小二乘原则，找到最满足上述反转对称性的地方，也就是b1的起始位置：
其中表示b1(τ1)的共轭；
(5.2)将接收信号r[k]代入(5.1)的式中，则找到接收信号中最满足上述反转对称性的地方；根据反转对称性定义基于最小二乘准则的代价函数Φ[τ]
对代价函数Φ(τ)进行相应滤波以减小噪声的干扰：
对接收信号进行一维搜索，通过求解下面的最大化问题：
由此确定FBMC/OQAM的接收信号的起始位置为
4.根据权利要求2所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，步骤(6)的具体包括如下步骤：
(6.1)确定接收信号的起始位置后，接收端经过串并转换、多相滤波、FFT、删除初始相位，得到频时格点(m0,n0)处的数据符号表示为
其中是原型滤波器p[k]的自模糊函数；是信道噪声在频时格点(m0,n0)处的噪声；
(6.2)考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，认为超出该范围的频时格点上的数据对目标频时格点上的数据造成的影响忽略不计；两列导频数据满足假设与导频列相邻的数据列都为零保护列，且忽略噪声，则导频符号表示为：
其中
(6.3)将与相乘，其中表示的共轭，由于因此得到：
所以CFO的估计值由下式得到：
(6.4)根据频偏估计值对接收信号r[k]进行频偏补偿，得到频偏补偿后的接收信号重新经过串并转换、多相滤波、FFT、删除初始相位得到重复计算得到一次迭代之后的频偏估计值。
5.根据权利要求4所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，步骤(6)还包括：
(6.5)重复(6.4)N-1次，得到N次迭代之后的频偏估计值，N为整数，且N≥1；通过迭代，降低与导频列相邻的数据列和信道噪声对估计结果的影响。
6.根据权利要求2所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，步骤(7)的具体包括如下步骤：
(7.1)经过串并转换、多相滤波、FFT、删除初始相位得到经过频偏补偿后的频时格点(m0,n0)处的数据符号假设频偏完全补偿，则表示为：
(7.2)定义为信道在子载波m上的频率响应，多径信道最大时延Δ＜＜M时，每一个子载波上的信道频域响应与邻近子载波上的信道频域响应近似相等；同时，对于时频域聚焦性良好的原型滤波器而言，在多径时延[0,Δ]范围内，Ap(-τ,0)≈Ap(0,0)；考虑目标频时格点周围Ωα,β＝{|(α,β)|≤(1,1)and(α,β)≠(0,0)}的频时格点，经过频偏补偿后的频时格点(m0,n0)处的数据符号进一步表示为：
其中表示所受到的内在干扰；
记并称为伪导频，信道估计的方法如下：
(7.3)首先假设与导频列相邻的数据列都为零保护列，进行信道参数的初步估计，对两列导频分别得到的信道估计结果求平均值，得到信道参数的初步估计值
7.根据权利要求6所述的一种FBMC/OQAM的导频设计和同步信道估计方法，其特征在于，步骤(7)还包括：
(7.4)步骤(7)还包括以下步骤：(7.4)利用对与导频列相邻的数据符号进行信道均衡，然后减去已知导频的干扰，对得到的第一、三、五列数据符号求平均值以减小噪声的干扰，并进行取实操作、判决，初步恢复出与导频列相邻的数据符号；
(7.5)对于位于第二列的导频，除了利用步骤(7.4)中解调出的与导频列相邻的数据符号外，还利用位于第四列导频的已知信息，重新计算伪导频进行信道估计得到 作为信道参数的最终估计值。
CN201710424321.8A 2017-06-07 2017-06-07 一种fbmc/oqam的导频设计和同步信道估计方法 Active CN107438038B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710424321.8A CN107438038B (zh) 2017-06-07 2017-06-07 一种fbmc/oqam的导频设计和同步信道估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710424321.8A CN107438038B (zh) 2017-06-07 2017-06-07 一种fbmc/oqam的导频设计和同步信道估计方法

Publications (2)

Publication Number Publication Date
CN107438038A true CN107438038A (zh) 2017-12-05
CN107438038B CN107438038B (zh) 2020-03-13

Family

ID=60458734

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710424321.8A Active CN107438038B (zh) 2017-06-07 2017-06-07 一种fbmc/oqam的导频设计和同步信道估计方法

Country Status (1)

Country Link
CN (1) CN107438038B (zh)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109167744A (zh) * 2018-11-06 2019-01-08 上海事凡物联网科技有限公司 一种相位噪声联合估计方法
CN109450827A (zh) * 2018-10-30 2019-03-08 武汉邮电科学研究院有限公司 基于下一代无线波形fbmc-oqam的信道估计方法及***
CN109600332A (zh) * 2018-11-22 2019-04-09 华中科技大学 一种co-fbmc/oqam的调制方法及***
CN109672637A (zh) * 2019-02-28 2019-04-23 重庆邮电大学 一种基于干扰修正迭代的fbmc信道估计方法
CN109726607A (zh) * 2018-12-29 2019-05-07 武汉万集信息技术有限公司 通信信道的切换方法及装置、存储介质、电子装置
CN110086495A (zh) * 2019-01-10 2019-08-02 上海事凡物联网科技有限公司 手持式无线自组网通信设备
CN110266627A (zh) * 2019-05-28 2019-09-20 上海交通大学 基于伪导频和判决反馈的cir和cfo联合估计方法
CN110581813A (zh) * 2019-08-12 2019-12-17 华中科技大学 一种多载波***导频信号的传输方法
CN110691055A (zh) * 2019-10-18 2020-01-14 中国人民解放军空军工程大学 Oqam/ofdm中的时频偏联合估计方法
CN110784425A (zh) * 2019-10-18 2020-02-11 西安烽火电子科技有限责任公司 一种频偏盲消除迭代方法
CN111555993A (zh) * 2020-05-20 2020-08-18 西安电子科技大学 Fbmc***中基于迭代预处理的信道估计方法
CN112054981A (zh) * 2020-08-17 2020-12-08 西安理工大学 一种实正交相移键控调制解调方法
CN112491756A (zh) * 2020-11-23 2021-03-12 西安电子科技大学 一种fbmc***中联合信道均衡的rls信道估计方法
CN113271279A (zh) * 2021-05-14 2021-08-17 成都爱瑞无线科技有限公司 一种窄带物联网随机接入信道的高精度检测方法
CN113300779A (zh) * 2021-04-26 2021-08-24 浙江工业大学 一种导频辅助的co-fbmc/oqam***相位噪声补偿方法
CN113542177A (zh) * 2021-05-31 2021-10-22 上海交通大学 解决脉冲幅度调制信号频偏混叠的方法及***
CN113556305A (zh) * 2021-06-19 2021-10-26 西安电子科技大学 适用于高频率选择性衰落的fbmc迭代信道均衡方法及***
CN114124624A (zh) * 2021-11-19 2022-03-01 成都爱瑞无线科技有限公司 一种基于导频的信道估计方法及***
CN114513396A (zh) * 2022-04-21 2022-05-17 科大天工智能装备技术（天津）有限公司 基于工业5g下的ofdm/oqam的导频方法和***
CN114598392A (zh) * 2022-03-22 2022-06-07 中山大学 多维光调制器支路间高精度同步方法
CN114884792A (zh) * 2022-05-27 2022-08-09 华中科技大学 一种高精度的多载波符号快速恢复方法、装置及***
CN114978820A (zh) * 2022-05-06 2022-08-30 南京信息工程大学 一种基于fbmc雷达通信一体化信号信道估计方法
CN115118558A (zh) * 2022-08-22 2022-09-27 北京科技大学 移动通信中ofdm/oqam的自干扰消除方法及***
US11496350B2 (en) 2018-03-27 2022-11-08 University Of South Carolina Dual-polarization FBMC in wireless communication systems

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101395815A (zh) * 2006-01-17 2009-03-25 株式会社Ntt都科摩 移动台、基站、通信***以及通信方法
CN104995886A (zh) * 2013-12-11 2015-10-21 华为技术有限公司 信息传输方法和装置
EP2913969B1 (en) * 2014-02-26 2016-07-06 NTT DoCoMo, Inc. Reduction of oqam intrinsic interference on pilot for channel estimation
CN106788651A (zh) * 2017-01-22 2017-05-31 西安交通大学 基于非正交多址接入的多地理区域广播***的信息传输方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101395815A (zh) * 2006-01-17 2009-03-25 株式会社Ntt都科摩 移动台、基站、通信***以及通信方法
CN104995886A (zh) * 2013-12-11 2015-10-21 华为技术有限公司 信息传输方法和装置
EP2913969B1 (en) * 2014-02-26 2016-07-06 NTT DoCoMo, Inc. Reduction of oqam intrinsic interference on pilot for channel estimation
CN106788651A (zh) * 2017-01-22 2017-05-31 西安交通大学 基于非正交多址接入的多地理区域广播***的信息传输方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
梁燕，等: "FBMC/OQAM***导频辅助信道估计方法概述", 《重庆邮电大学学报（自然科学版本）》 *
龚帅，等: "基于导频的FBMC/OQAM***的信道估计方法", 《上海师范大学学报（自然科学版）》 *

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11496350B2 (en) 2018-03-27 2022-11-08 University Of South Carolina Dual-polarization FBMC in wireless communication systems
CN109450827A (zh) * 2018-10-30 2019-03-08 武汉邮电科学研究院有限公司 基于下一代无线波形fbmc-oqam的信道估计方法及***
CN109167744B (zh) * 2018-11-06 2021-05-14 上海事凡物联网科技有限公司 一种相位噪声联合估计方法
CN109167744A (zh) * 2018-11-06 2019-01-08 上海事凡物联网科技有限公司 一种相位噪声联合估计方法
CN109600332A (zh) * 2018-11-22 2019-04-09 华中科技大学 一种co-fbmc/oqam的调制方法及***
CN109726607A (zh) * 2018-12-29 2019-05-07 武汉万集信息技术有限公司 通信信道的切换方法及装置、存储介质、电子装置
CN110086495A (zh) * 2019-01-10 2019-08-02 上海事凡物联网科技有限公司 手持式无线自组网通信设备
CN109672637A (zh) * 2019-02-28 2019-04-23 重庆邮电大学 一种基于干扰修正迭代的fbmc信道估计方法
CN109672637B (zh) * 2019-02-28 2021-09-28 重庆邮电大学 一种基于干扰修正迭代的fbmc信道估计方法
CN110266627A (zh) * 2019-05-28 2019-09-20 上海交通大学 基于伪导频和判决反馈的cir和cfo联合估计方法
CN110266627B (zh) * 2019-05-28 2021-07-13 上海交通大学 基于伪导频和判决反馈的cir和cfo联合估计方法
CN110581813A (zh) * 2019-08-12 2019-12-17 华中科技大学 一种多载波***导频信号的传输方法
CN110784425B (zh) * 2019-10-18 2021-12-10 西安烽火电子科技有限责任公司 一种频偏盲消除迭代方法
CN110784425A (zh) * 2019-10-18 2020-02-11 西安烽火电子科技有限责任公司 一种频偏盲消除迭代方法
CN110691055A (zh) * 2019-10-18 2020-01-14 中国人民解放军空军工程大学 Oqam/ofdm中的时频偏联合估计方法
CN111555993A (zh) * 2020-05-20 2020-08-18 西安电子科技大学 Fbmc***中基于迭代预处理的信道估计方法
CN112054981A (zh) * 2020-08-17 2020-12-08 西安理工大学 一种实正交相移键控调制解调方法
CN112054981B (zh) * 2020-08-17 2022-07-15 西安理工大学 一种实正交相移键控调制解调方法
CN112491756A (zh) * 2020-11-23 2021-03-12 西安电子科技大学 一种fbmc***中联合信道均衡的rls信道估计方法
CN113300779A (zh) * 2021-04-26 2021-08-24 浙江工业大学 一种导频辅助的co-fbmc/oqam***相位噪声补偿方法
CN113300779B (zh) * 2021-04-26 2022-04-08 浙江工业大学 一种导频辅助的co-fbmc/oqam***相位噪声补偿方法
CN113271279A (zh) * 2021-05-14 2021-08-17 成都爱瑞无线科技有限公司 一种窄带物联网随机接入信道的高精度检测方法
CN113542177A (zh) * 2021-05-31 2021-10-22 上海交通大学 解决脉冲幅度调制信号频偏混叠的方法及***
CN113542177B (zh) * 2021-05-31 2022-04-08 上海交通大学 解决脉冲幅度调制信号频偏混叠的方法及***
CN113556305A (zh) * 2021-06-19 2021-10-26 西安电子科技大学 适用于高频率选择性衰落的fbmc迭代信道均衡方法及***
CN114124624A (zh) * 2021-11-19 2022-03-01 成都爱瑞无线科技有限公司 一种基于导频的信道估计方法及***
CN114124624B (zh) * 2021-11-19 2023-11-17 成都爱瑞无线科技有限公司 一种基于导频的信道估计方法及***
CN114598392A (zh) * 2022-03-22 2022-06-07 中山大学 多维光调制器支路间高精度同步方法
CN114598392B (zh) * 2022-03-22 2024-03-15 中山大学 多维光调制器支路间高精度同步方法
CN114513396A (zh) * 2022-04-21 2022-05-17 科大天工智能装备技术（天津）有限公司 基于工业5g下的ofdm/oqam的导频方法和***
CN114978820B (zh) * 2022-05-06 2023-04-21 南京信息工程大学 一种基于fbmc雷达通信一体化信号信道估计方法
CN114978820A (zh) * 2022-05-06 2022-08-30 南京信息工程大学 一种基于fbmc雷达通信一体化信号信道估计方法
CN114884792A (zh) * 2022-05-27 2022-08-09 华中科技大学 一种高精度的多载波符号快速恢复方法、装置及***
CN114884792B (zh) * 2022-05-27 2024-05-24 华中科技大学 一种高精度的多载波符号快速恢复方法、装置及***
CN115118558B (zh) * 2022-08-22 2022-11-08 北京科技大学 移动通信中ofdm/oqam的自干扰消除方法及***
CN115118558A (zh) * 2022-08-22 2022-09-27 北京科技大学 移动通信中ofdm/oqam的自干扰消除方法及***

Also Published As

Publication number Publication date
CN107438038B (zh) 2020-03-13

Similar Documents

Publication Publication Date Title
CN107438038A (zh) 一种fbmc/oqam的导频设计和同步信道估计方法
CN102185822A (zh) 一种ofdm/oqam***及其时频同步方法
CN102215205B (zh) 一种ofdm/oqam***及其时频同步方法
KR20160022278A (ko) 기준 신호 송수신 방법 및 장치
CN107086974A (zh) 一种高动态环境下的ofdm同步方法
CN107302419B (zh) 一种用于mimo-ofdm***的低复杂度检测方法
CN102215204A (zh) 基于反馈迭代的ofdm/oqam***及其时频同步方法
CN109861939B (zh) 一种oqpsk频域均衡无线数据传输方法
KR20040110038A (ko) 주파수 영역에 숨겨진 파일럿 신호를 이용하여신호처리하는 ｔｄｓ-ｏｆｄｍ 수신 시스템 및 그의 신호처리 방법
WO2010080361A1 (en) Method and apparatus for tracking fft window position
CN106788935A (zh) Fbmc***中高频谱效率的导频设计和信道估计方法
CN101119350A (zh) 正交频分复用***、快速同步的方法和发送端设备
CN101001235A (zh) 一种弱能量并行pn序列时间同步与频率同步方法
KR20180052003A (ko) Ofdm 시스템에서 단일 탭 등화기를 이용한 부반송파의 왜곡 보상 방법 및 그를 위한 장치
CN101369882B (zh) 一种基于发送端预纠正的频率同步方法
CN101729479A (zh) 一种基于ofdm信号循环平稳特性的盲信道估计方法
CN104253772A (zh) 正交频分复用***的信道估计方法
TH et al. Improvement in Data Transmission Efficiency in Mobile 5G New Radio System using Filter Bank Multicarrier Signals
EP3384650A1 (en) Wireless data communication based on discrete cosine transformation
CN101783782B (zh) 一种自动适应信道特性变化的上行多址接入方法
CN101087285A (zh) 一种新一代无线移动通信***中的同步方法
CN107204953B (zh) 一种cp-fbmc通信***中的盲频偏估计方法
Puska et al. Matched filter time and frequency synchronization method for OFDM systems using PN-sequence preambles
CN103188196B (zh) 正交频分复用***的同步方法
Nguyen Robust time and frequency synchronization in 802.11 a communication wireless system

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