CN114244670B - Blind channel estimation method and system based on channel coding assistance - Google Patents
Blind channel estimation method and system based on channel coding assistance Download PDFInfo
- Publication number
- CN114244670B CN114244670B CN202111489753.XA CN202111489753A CN114244670B CN 114244670 B CN114244670 B CN 114244670B CN 202111489753 A CN202111489753 A CN 202111489753A CN 114244670 B CN114244670 B CN 114244670B
- Authority
- CN
- China
- Prior art keywords
- channel estimation
- channel
- module
- result
- blind
- 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 78
- 238000004891 communication Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 56
- 230000008569 process Effects 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 14
- 238000010586 diagram Methods 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 230000003321 amplification Effects 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims 1
- 238000003780 insertion Methods 0.000 abstract description 3
- 230000037431 insertion Effects 0.000 abstract description 3
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000012804 iterative process Methods 0.000 description 4
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- 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
-
- 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
-
- 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/0238—Channel estimation using blind estimation
-
- 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/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2634—Inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators in combination with other circuits for modulation
- H04L27/2636—Inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators in combination with other circuits for modulation with FFT or DFT modulators, e.g. standard single-carrier frequency-division multiple access [SC-FDMA] transmitter or DFT spread orthogonal frequency division multiplexing [DFT-SOFDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Error Detection And Correction (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及一种信道编码辅助的盲信道估计方法及系统,属于通信信号处理领域。The invention relates to a channel coding assisted blind channel estimation method and system, belonging to the field of communication signal processing.
背景技术Background Art
正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)技术由于具有高频谱效率和抗多径衰落等突出优势已广泛应用于各种标准,如:802.11,长期演进技术(Long Term Evolution,LTE),WiMAX等。为了抵消多径衰落的影响,OFDM通信系统需要进行信道估计和均衡。传统的信道估计算法主要为基于频域导频辅助的非盲信道估计方法,其包括最小二乘(Least Square,LS)算法、最小均方误差(Minimum Mean Square Error,MMSE)算法、线性最小均方误差(LMMSE)算法和离散傅立叶变换(Discrete FourierTransform,DFT)算法等。上述算法依赖于插入到OFDM信号中的导频符号,而导频符号的插入会导致OFDM系统时频资源的消耗以及传输速率的降低。针对这一不足,可在OFDM系统中采用盲信道估计算法,此类方法无需在发送信号中插入导频符号,因此可以明显提升OFDM系统的频带利用率,其主要分为两大类:统计性方法和确定性方法。前者利用了发送及接收信号的统计特性,如相关函数、相关矩阵等,实现信道估计的功能,其中使用较多的包括子空间方法与基于线性预编码的方法;后者利用了发送调制信号的固有特性,此方法对应的信号处理模块一般放在接收端DFT操作之后,其中使用较多的包括最大似然盲估计方法,基于加性导频的盲估计算法,和基于接收分集的盲估计算法。然而,上述所有盲估计算法均具有较大的计算复杂度,且其信道估计精度不高。Orthogonal Frequency Division Multiplexing (OFDM) technology has been widely used in various standards, such as 802.11, Long Term Evolution (LTE), WiMAX, etc., due to its outstanding advantages such as high spectrum efficiency and resistance to multipath fading. In order to offset the impact of multipath fading, OFDM communication systems need to perform channel estimation and equalization. Traditional channel estimation algorithms are mainly non-blind channel estimation methods based on frequency domain pilot assistance, including the least squares (LS) algorithm, the minimum mean square error (MMSE) algorithm, the linear minimum mean square error (LMMSE) algorithm, and the discrete Fourier transform (DFT) algorithm. The above algorithms rely on the pilot symbols inserted into the OFDM signal, and the insertion of pilot symbols will lead to the consumption of time and frequency resources of the OFDM system and the reduction of the transmission rate. To address this shortcoming, blind channel estimation algorithms can be used in OFDM systems. Such methods do not require the insertion of pilot symbols in the transmitted signal, so they can significantly improve the frequency band utilization of the OFDM system. They are mainly divided into two categories: statistical methods and deterministic methods. The former utilizes the statistical characteristics of the transmitted and received signals, such as correlation functions, correlation matrices, etc., to achieve the channel estimation function. Among them, the more commonly used methods include subspace methods and methods based on linear precoding; the latter utilizes the inherent characteristics of the transmitted modulated signal. The signal processing module corresponding to this method is generally placed after the DFT operation at the receiving end. Among them, the more commonly used methods include maximum likelihood blind estimation methods, blind estimation algorithms based on additive pilots, and blind estimation algorithms based on receiving diversity. However, all of the above blind estimation algorithms have large computational complexity and their channel estimation accuracy is not high.
发明内容Summary of the invention
针对现有OFDM系统信道估计方法存在的下述技术缺陷:(一)传统非盲估计方法中插入导频符号造成频谱利用率下降的问题;(二)传统盲估计方法中计算复杂度较高和估计性能较差的不足。本发明的主要目的是提供一种基于信道编码辅助的盲信道估计方法及系统,将译码器输出的硬判决比特信息迭代反馈给信道估计模块,将所述硬判决比特信息作为发送数据的先验信息,代替原有导频符号的功能,进而使本发明具有如下主要优点:(一)无需在发送端插入导频符号,进而避免插入导频符号造成频谱利用率下降的问题,显著提高OFDM系统的吞吐量;(二)在非盲信道估计方法基础上,通过将所述硬判决比特信息作为发送数据的先验信息,代替原有导频符号的功能,实现盲信道估计,使本发明具有和非盲信道估计方法相近的估计复杂度和估计精度,能够达到近似理论极限的误码率性能。Aiming at the following technical defects in the existing OFDM system channel estimation method: (i) the problem of reduced spectrum utilization caused by inserting pilot symbols in the traditional non-blind estimation method; (ii) the shortcomings of high computational complexity and poor estimation performance in the traditional blind estimation method. The main purpose of the present invention is to provide a blind channel estimation method and system based on channel coding assistance, iteratively feed back the hard decision bit information output by the decoder to the channel estimation module, and use the hard decision bit information as the prior information of the transmitted data to replace the function of the original pilot symbol, so that the present invention has the following main advantages: (i) there is no need to insert pilot symbols at the transmitting end, thereby avoiding the problem of reduced spectrum utilization caused by inserting pilot symbols, and significantly improving the throughput of the OFDM system; (ii) on the basis of the non-blind channel estimation method, by using the hard decision bit information as the prior information of the transmitted data to replace the function of the original pilot symbol, blind channel estimation is realized, so that the present invention has the estimation complexity and estimation accuracy similar to the non-blind channel estimation method, and can achieve the bit error rate performance close to the theoretical limit.
本发明的目的是通过下述技术方案实现的。The purpose of the present invention is achieved through the following technical solutions.
本发明公开的一种基于信道编码辅助的盲信道估计方法,包括如下步骤:The present invention discloses a blind channel estimation method based on channel coding assistance, comprising the following steps:
步骤一,发送具有如下帧格式的信号。Step 1: Send a signal with the following frame format.
所述帧格式为:长度为Nsym个符号的OFDM帧中,第一个符号插入导频信息,后续(Nsym-1)个OFDM符号无需插入额外的导频符号。The frame format is: in an OFDM frame with a length of N sym symbols, pilot information is inserted into the first symbol, and no additional pilot symbols need to be inserted into the subsequent (N sym -1) OFDM symbols.
步骤二,接收步骤一发送的信号,提取第一个OFDM符号中的导频符号,进行频域信道估计。Step 2: Receive the signal sent in
以x(n)表示发送的时域信号,h(n)表示多径信道的抽头,信道抽头数为L,ω(n)表示噪声,y(n)表示接收信号,*表示线性卷积,接收信号如下式所示,Let x(n) represent the transmitted time domain signal, h(n) represent the tap of the multipath channel, the number of channel taps is L, ω(n) represent the noise, y(n) represent the received signal, * represent linear convolution, and the received signal is shown in the following formula:
去除循环前缀后,接收OFDM符号的有效序列如下式所示,After removing the cyclic prefix, the effective sequence of the received OFDM symbol is as follows,
其中表示循环卷积,N为OFDM系统的FFT点数。记Nsc为OFDM子载波个数,X(k)为x(n)的DFT变换,H(k)为h(n)的DFT变换,Y(k)为y(n)的DFT变换,W(k)为ω(n)的DFT变换,接收数据的频域表达式如下式所示,in Represents circular convolution, N is the number of FFT points in the OFDM system. Let Nsc be the number of OFDM subcarriers, X(k) be the DFT transform of x(n), H(k) be the DFT transform of h(n), Y(k) be the DFT transform of y(n), W(k) be the DFT transform of ω(n), and the frequency domain expression of the received data is as follows:
Y(k)=X(k)·H(k)+W(k),k=0,...,Nsc-1 (3)Y(k)=X(k)·H(k)+W(k),k=0,...,N sc -1 (3)
记X=[X(0),...,X(Nsc-1)]T,H=[H(0),...,H(Nsc-1)]T,Y=[Y(0),...,Y(Nsc-1)]T,W=[W(0),...,W(Nsc-1)]T,接收数据频域表达式的矩阵表示为,Remember X=[X(0),...,X(N sc -1)] T , H=[H(0),...,H(N sc -1)] T ,Y=[Y(0),...,Y(N sc -1)] T ,W=[W(0 ),...,W(N sc -1)] T , the matrix of the received data frequency domain expression is expressed as,
当X(k)已知时通过LS估计算法得出H(k)的估计值。记Xp(k)为导频符号,Xp=[Xp(0),...,Xp(Nsc-1)]T,H的LS估计值如下式所示,When X(k) is known, the estimated value of H(k) is obtained by the LS estimation algorithm. Let Xp (k) be the pilot symbol, Xp = [ Xp (0), ..., Xp ( Nsc -1)] T , LS estimate of H As shown in the following formula,
根据上式得,第一个OFDM符号LS信道估计结果中各个元素的表达式为,According to the above formula, the LS channel estimation result of the first OFDM symbol is Each element The expression of is,
角标中的(i)表示第i个OFDM符号。The (i) in the superscript represents the i-th OFDM symbol.
步骤三,使用第一个OFDM符号的信道估计结果,对第二个OFDM符号进行信道均衡,均衡的结果进行后续数字解调和信道译码,得到译码器变量节点输出结果。步骤三的流程不只限于对第二个OFDM符号进行处理,从第二个OFDM符号起,所有OFDM符号进行信道估计时均需要进行相同操作,即从步骤三开始进入本发明所述方法的迭代过程。为了强调一般性,后续记步骤三处理的OFDM符号为第i个OFDM符号,i=2,3,...。Step three, using the channel estimation result of the first OFDM symbol, channel equalize the second OFDM symbol, and perform subsequent digital demodulation and channel decoding on the equalized result to obtain the decoder variable node output result. The process of step three is not limited to processing the second OFDM symbol. From the second OFDM symbol onwards, the same operation needs to be performed when channel estimating all OFDM symbols, that is, the iterative process of the method of the present invention begins from step three. In order to emphasize the generality, the OFDM symbol processed in step three is subsequently recorded as the i-th OFDM symbol, i = 2, 3,...
使用第i-1个OFDM符号的信道估计结果对第i个OFDM符号进行信道均衡,从第二个OFDM符号开始信号中不再插入导频符号,信道均衡后的结果表达式为,Use the channel estimation result of the i-1th OFDM symbol to perform channel equalization on the i-th OFDM symbol. From the second OFDM symbol onwards, no pilot symbol is inserted into the signal. The result after channel equalization is The expression is,
之后进行后续的数字解调过程,得到各个比特的软信息m=0,...,M-1,M为数字解调后的比特数,将比特软信息传递给LDPC译码器,进行LDPC信道译码,得到包括校验位在内所有变量节点的译码输出结果m=0,...,M-1。After that, the subsequent digital demodulation process is carried out to obtain the soft information of each bit. m=0,...,M-1, M is The number of bits after digital demodulation, the bit soft information is passed to the LDPC decoder, LDPC channel decoding is performed, and the decoding output results of all variable nodes including the check bit are obtained m=0,...,M-1.
步骤四,将步骤三译码器变量节点的输出结果进行硬判决得到比特信息,将所述比特信息重新进行数字调制,得到的调制符号后文称为迭代导频符号(Iterative PilotSymbols,IPS)。由于译码器带来的编码增益,此IPS符号相比步骤三信道均衡的结果准确度显著提高,使用此等效的导频符号辅助盲信道估计。由于信号内部实际并未插入导频符号,本发明所述通信系统相比于采用非盲信道估计的通信系统提高频谱利用率。
对LDPC译码输出结果进行硬判决得到序列b(i)(m),m=0,...,M-1,此结果为0、1比特形式,将译码的硬判决输出结果重新进行数字调制,其调制方式应与信号原本数字调制方式相同,调制结果看作等效的导频符号IPS,记为i=2,3,...。Output results of LDPC decoding A hard decision is made to obtain the sequence b (i) (m), m = 0, ..., M-1. This result is in the form of 0 and 1 bits. The decoded hard decision output result is digitally modulated again. Its modulation method should be the same as the original digital modulation method of the signal. The modulation result is regarded as an equivalent pilot symbol IPS, which is recorded as i=2,3,...
步骤五,把步骤三得到的迭代导频符号IPS视为第i个OFDM符号的先验信息,以IPS为导频符号,基于非盲信道估计算法进行第i个OFDM符号的信道估计,实现以非盲估计算法进行盲信道估计的目的,相比于传统盲信道估计降低实现复杂度且提高估计精度。在估计算法流程上,首先对信号进行基于LS算法的信道估计,再对LS算法的估计结果进行基于DFT变换的去噪处理,得到DFT信道估计算法的估计结果。本步骤中进行基于DFT变换的去噪处理有两个目的,其一是提高信道估计的准确性,提升通信系统的误码性能;其二是降低均衡过程中引入的噪声,优化迭代结构的收敛性,保证迭代信道估计结构在较低信噪比下能够收敛。
将IPS反馈给信道估计模块,把IPS视作导频符号,基于LS算法估计出第i个OFDM符号的信道CFR估计值 Feedback IPS to the channel estimation module, treat IPS as a pilot symbol, and estimate the channel CFR estimate of the i-th OFDM symbol based on the LS algorithm
将LS信道估计算法得到的CFR结果进行基于DFT变换的去噪算法处理。LS算法的估计结果先进行IDFT变换,得到信道脉冲响应如下式所示,The CFR result obtained by the LS channel estimation algorithm is processed by a denoising algorithm based on DFT transform. The estimation result of the LS algorithm is first subjected to IDFT transform to obtain the channel impulse response As shown in the following formula,
在信道脉冲响应中保留有效抽头信息,将其余部分的噪声置零,得到去噪后的信道脉冲响应如下式所示,The channel impulse response The effective tap information is retained, and the noise of the rest is set to zero to obtain the denoised channel impulse response As shown in the following formula,
再将进行DFT变换,得到第i个OFDM符号基于DFT信道估计算法的CFR估计结果 Then Perform DFT transformation to obtain the CFR estimation result of the i-th OFDM symbol based on the DFT channel estimation algorithm
步骤六,重复步骤三到步骤五,由于LDPC译码对数据SNR的放大作用以及DFT去噪算法对噪声的抑制作用,每次迭代的信道估计结果CFR将逐渐精确,在接收信噪比足够大的情况下达到和非盲估计方法几乎相同的估计精确度和误码性能。Step 6: Repeat
本发明还公开基于信道编码辅助的盲信道估计系统,用于实现所述基于信道编码辅助的盲信道估计方法。所述基于信道编码辅助的盲信道估计系统主要由信道估计模块、信道均衡模块、数字解调模块、信道译码模块、硬判决模块和数字调制模块组成。信道估计模块包含LS估计模块和基于DFT变换的去噪模块。信道估计模块涉及步骤二和步骤五,其中基于DFT变换的去噪模块降低信道估计结果的噪声方差,提升信道估计性能,优化迭代结构的收敛效果。信道均衡模块和数字解调模块涉及步骤三,信道均衡模块根据信道估计的结果对接收信号进行补偿,而数字解调模块对补偿后的信号进行星座图解映射得到比特软信息,为信道译码模块提供输入。信道译码模块涉及步骤三,通过信道译码提升通信系统的抗噪声性能,并且保证迭代结构的收敛性。硬判决模块和数字调制模块涉及步骤四,硬判决模块将信道译码输出的软信息结果判决为0、1比特,数字调制模块将硬判决输出的结果重新调制为数字调制符号,为信道估计模块提供等效导频符号IPS。The present invention also discloses a blind channel estimation system based on channel coding assistance, which is used to implement the blind channel estimation method based on channel coding assistance. The blind channel estimation system based on channel coding assistance is mainly composed of a channel estimation module, a channel equalization module, a digital demodulation module, a channel decoding module, a hard decision module and a digital modulation module. The channel estimation module includes an LS estimation module and a denoising module based on DFT transformation. The channel estimation module involves
作为优选,在信道编码过程中,使用纠错性能优异的低密度奇偶校验(Low-Density Parity-Check,LDPC)码,增强信道估计及通信系统的抗噪声性能。而基于LDPC编码辅助的迭代式盲信道估计方法简称为LDPC-ICE(LDPC-code-aided iterative channelestimation)。Preferably, in the channel coding process, a low-density parity check (LDPC) code with excellent error correction performance is used to enhance the channel estimation and the anti-noise performance of the communication system. The iterative blind channel estimation method based on LDPC coding assistance is referred to as LDPC-ICE (LDPC-code-aided iterative channel estimation).
有益效果:Beneficial effects:
1、本发明公开的一种基于信道编码辅助的盲信道估计方法及系统,使用非盲信道估计相同的信道估计算法,但没有在信号中插入导频符号,降低盲信道估计的计算复杂度,提升频谱利用率。将译码器变量节点的输出经过硬判决和调制,得到IPS用于辅助估计模块的计算,这种信道估计算法和基于导频符号辅助的非盲估计方法完全相同。相比于非盲信道估计结构,本发明提出的盲信道估计结构只需要增加硬判决器和调制模块。在硬件实现过程中,硬判决器只涉及对于符号位的判断,调制模块也只需要简单的查找表结构,占用硬件资源极少。因此本发明提出的盲信道估计系统可以达到与非盲信道估计系统相近的计算复杂度,远低于传统盲信道估计算法。1. The present invention discloses a blind channel estimation method and system based on channel coding assistance, which uses the same channel estimation algorithm as the non-blind channel estimation, but does not insert pilot symbols into the signal, thereby reducing the computational complexity of the blind channel estimation and improving the spectrum utilization rate. The output of the decoder variable node is subjected to hard decision and modulation to obtain IPS for the calculation of the auxiliary estimation module. This channel estimation algorithm is exactly the same as the non-blind estimation method based on pilot symbol assistance. Compared with the non-blind channel estimation structure, the blind channel estimation structure proposed in the present invention only needs to add a hard decision maker and a modulation module. In the hardware implementation process, the hard decision maker only involves the judgment of the sign bit, and the modulation module only needs a simple lookup table structure, which occupies very few hardware resources. Therefore, the blind channel estimation system proposed in the present invention can achieve a computational complexity similar to that of the non-blind channel estimation system, which is much lower than the traditional blind channel estimation algorithm.
2、本发明公开的一种基于信道编码辅助的盲信道估计方法及系统,在传统OFDM接收机的基础上增加反馈回路,反馈回路中包括硬判决器和调制模块,结构较为简单,其他信号处理流程相比传统流程未做较大改动,易于移植到已有通信系统中。2. The present invention discloses a blind channel estimation method and system based on channel coding assistance, which adds a feedback loop on the basis of a traditional OFDM receiver. The feedback loop includes a hard decision maker and a modulation module. The structure is relatively simple, and other signal processing processes are not significantly changed compared to traditional processes, and can be easily transplanted into existing communication systems.
3、本发明公开的一种基于信道编码辅助的盲信道估计方法及系统,在迭代反馈的回路中引入信道译码模块和基于DFT变换的去噪模块,前者确保迭代结构能够收敛,后者优化迭代结构在较低信噪比的收敛效果。在较高信噪比下,本发明所述迭代式信道估计性能在较高信噪比下都是收敛的。3. The present invention discloses a blind channel estimation method and system based on channel coding assistance, which introduces a channel decoding module and a denoising module based on DFT transformation in the iterative feedback loop. The former ensures that the iterative structure can converge, and the latter optimizes the convergence effect of the iterative structure at a lower signal-to-noise ratio. Under a higher signal-to-noise ratio, the iterative channel estimation performance of the present invention is converged under a higher signal-to-noise ratio.
4、本发明公开的一种基于信道编码辅助的盲信道估计方法及系统,采用基于DFT变换的去噪模块,降低估计结果的噪声方差,提升信道估计性能。4. The present invention discloses a blind channel estimation method and system based on channel coding assistance, which adopts a denoising module based on DFT transformation to reduce the noise variance of the estimation result and improve the channel estimation performance.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1是本发明所述的一种基于信道编码辅助的盲信道估计方法及系统流程示意图;FIG1 is a schematic diagram of a blind channel estimation method and system flow based on channel coding assistance according to the present invention;
图2是本发明所述OFDM通信系统接收机信号处理结构框图;2 is a block diagram of a receiver signal processing structure of an OFDM communication system according to the present invention;
图3是基于DFT变换的去噪模块示意图;FIG3 is a schematic diagram of a denoising module based on DFT transformation;
图4是本发明所述基于信道编码辅助的盲信道估计方法与传统非盲信道估计方法的BER性能曲线对比图;FIG4 is a comparison diagram of BER performance curves of the blind channel estimation method based on channel coding assistance according to the present invention and the traditional non-blind channel estimation method;
图5是本发明所述基于信道编码辅助的盲信道估计方法与传统非盲信道估计方法的NMSE性能曲线对比图。FIG5 is a comparison diagram of NMSE performance curves of the blind channel estimation method based on channel coding assistance according to the present invention and the traditional non-blind channel estimation method.
具体实施方式DETAILED DESCRIPTION
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例子,并且参照附图,对本发明进一步详细说明。In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in combination with specific implementation examples and with reference to the accompanying drawings.
本实施例中系统参数如下表所示:The system parameters in this embodiment are shown in the following table:
本实施例所提出的盲估计方法整体流程图如图1所示,接收机系统框图如图2所示,The overall flow chart of the blind estimation method proposed in this embodiment is shown in FIG1 , and the block diagram of the receiver system is shown in FIG2 .
本实施例公开的一种基于信道编码辅助的盲信道估计方法方法,具体实现步骤如下:This embodiment discloses a blind channel estimation method based on channel coding assistance, and the specific implementation steps are as follows:
步骤一,发送具有如下帧格式的信号。Step 1: Send a signal with the following frame format.
所述帧格式为:长度为120个符号的OFDM帧中,第一个符号插入导频信息,后续119个OFDM符号无需插入额外的导频符号。The frame format is: in an OFDM frame with a length of 120 symbols, the first symbol is inserted with pilot information, and no additional pilot symbols need to be inserted into the subsequent 119 OFDM symbols.
步骤二,接收步骤一发送的信号,提取第一个OFDM符号中的导频符号,进行频域信道估计。Step 2: Receive the signal sent in
以x(n)表示发送的时域信号,h(n)表示多径信道的抽头,信道抽头数为10,ω(n)表示噪声,y(n)表示接收信号,*表示线性卷积,接收信号如下式所示,Let x(n) represent the transmitted time domain signal, h(n) represent the tap of the multipath channel, the number of channel taps is 10, ω(n) represent the noise, y(n) represent the received signal, * represent linear convolution, and the received signal is shown in the following formula:
去除循环前缀后,接收OFDM符号的有效序列如下式表示,After removing the cyclic prefix, the effective sequence of the received OFDM symbol is expressed as follows:
其中表示循环卷积,OFDM系统的FFT点数为1024。实施例中OFDM子载波个数为960,X(k)为x(n)的DFT变换,H(k)为h(n)的DFT变换,Y(k)为y(n)的DFT变换,W(k)为ω(n)的DFT变换,接收数据的频域表达式如下式所示,in represents circular convolution, and the number of FFT points of the OFDM system is 1024. In the embodiment, the number of OFDM subcarriers is 960, X(k) is the DFT transform of x(n), H(k) is the DFT transform of h(n), Y(k) is the DFT transform of y(n), W(k) is the DFT transform of ω(n), and the frequency domain expression of the received data is shown in the following formula:
Y(k)=X(k)·H(k)+W(k),k=0,...,959 (14)Y(k)=X(k)·H(k)+W(k),k=0,...,959 (14)
记X=[X(0),...,X(959)]T,H=[H(0),...,H(959)]T,Y=[Y(0),...,Y(959)]T,W=[W(0),...,W(959)]T,接收数据频域表达式的矩阵表示为,Let X = [X(0),...,X(959)] T , H=[H(0),...,H(959)] T ,Y=[Y(0),...,Y(959)] T ,W=[W(0),...,W(959)] T ,the matrix of the received data frequency domain expression is:
当X(k)已知时通过LS估计算法得出H(k)的估计值。记Xp(k)为导频符号,Xp=[Xp(0),...,Xp(959)]T,H的LS估计值如下式所示,When X(k) is known, the estimated value of H(k) is obtained by the LS estimation algorithm. Let Xp (k) be the pilot symbol, Xp = [ Xp (0),..., Xp (959)] T , LS estimate of H As shown in the following formula,
根据上式得,第一个OFDM符号LS信道估计结果中各个元素的表达式为,According to the above formula, the LS channel estimation result of the first OFDM symbol is Each element The expression of is,
角标中的(i)表示第i个OFDM符号。The (i) in the superscript represents the i-th OFDM symbol.
步骤三,使用第一个OFDM符号的信道估计结果,对第二个OFDM符号进行信道均衡,均衡的结果进行后续数字解调和信道译码,得到译码器变量节点输出结果。步骤三的流程不只限于对第二个OFDM符号进行处理,从第二个OFDM符号起,所有OFDM符号进行信道估计时均需要进行相同操作,即从步骤三开始进入本发明所述方法的迭代过程。为了强调一般性,后续记步骤三处理的OFDM符号为第i个OFDM符号,i=2,3,...120。Step three, using the channel estimation result of the first OFDM symbol, perform channel equalization on the second OFDM symbol, and perform subsequent digital demodulation and channel decoding on the equalization result to obtain the decoder variable node output result. The process of step three is not limited to processing the second OFDM symbol. From the second OFDM symbol onwards, the same operation needs to be performed when channel estimating all OFDM symbols, that is, the iterative process of the method of the present invention begins from step three. In order to emphasize the generality, the OFDM symbol processed in step three is subsequently recorded as the i-th OFDM symbol, i = 2, 3, ... 120.
使用第i-1个OFDM符号的信道估计结果对第i个OFDM符号进行信道均衡,从第二个OFDM符号开始信号中不再插入导频符号,信道均衡后的结果表达式为,Use the channel estimation result of the i-1th OFDM symbol to perform channel equalization on the i-th OFDM symbol. From the second OFDM symbol onwards, no pilot symbol is inserted into the signal. The result after channel equalization is The expression is,
之后进行后续的数字解调过程,得到各个比特的软信息m=0,...,1919,数字解调后的比特数为1920,将比特软信息传递给LDPC译码器,进行LDPC信道译码,得到包括校验位在内所有变量节点的译码输出结果m=0,...,1919。After that, the subsequent digital demodulation process is carried out to obtain the soft information of each bit. m=0,...,1919, The number of bits after digital demodulation is 1920. The bit soft information is passed to the LDPC decoder for LDPC channel decoding to obtain the decoding output results of all variable nodes including the check bit. m=0,...,1919.
步骤四,将步骤三译码器变量节点的输出结果进行硬判决得到比特信息,将所述比特信息重新进行数字调制,得到的调制符号后文称为迭代导频符号(Iterative PilotSymbols,IPS)。由于译码器带来的编码增益,此IPS符号相比步骤三信道均衡的结果准确度显著提高,使用此等效的导频符号辅助盲信道估计。由于信号内部实际并未插入导频符号,本发明所述通信系统相比于采用非盲信道估计的通信系统提高了频谱利用率。
对LDPC译码输出结果进行硬判决得到序列b(i)(m),m=0,...,1919,此结果为0、1比特形式,将译码的硬判决输出结果重新进行数字调制,其调制方式应与信号原本数字调制方式相同,调制结果看作等效的导频符号IPS,记为i=2,3,...120。Output results of LDPC decoding A hard decision is made to obtain the sequence b (i) (m), m = 0, ..., 1919. This result is in the form of 0 and 1 bits. The decoded hard decision output result is re-digitally modulated. Its modulation method should be the same as the original digital modulation method of the signal. The modulation result is regarded as an equivalent pilot symbol IPS, which is recorded as i=2,3,...120.
步骤五,把步骤三得到的迭代导频符号IPS视为第i个OFDM符号的先验信息,以IPS为导频符号,基于非盲信道估计算法进行第i个OFDM符号的信道估计,实现以非盲估计算法进行盲信道估计的目的,相比于传统盲信道估计降低了实现复杂度且提高了估计精度。在估计算法流程上,首先对信号进行基于LS算法的信道估计,再对LS算法的估计结果进行基于DFT变换的去噪处理,得到DFT信道估计算法的估计结果,基于DFT变换的去噪模块示意图如图3所示。本步骤中进行基于DFT变换的去噪处理有两个目的,其一是提高信道估计的准确性,提升通信系统的误码性能;其二是降低均衡过程中引入的噪声,优化迭代结构的收敛性,保证迭代信道估计结构在较低信噪比下能够收敛,如图4图5所示。Step 5: The iterative pilot symbol IPS obtained in
将IPS反馈给信道估计模块,把IPS视作导频符号,基于LS算法估计出第i个OFDM符号的信道CFR估计值 Feedback IPS to the channel estimation module, treat IPS as a pilot symbol, and estimate the channel CFR estimate of the i-th OFDM symbol based on the LS algorithm
将LS信道估计算法得到的CFR结果进行基于DFT变换的去噪算法处理。LS算法的估计结果先进行IDFT变换,得到信道脉冲响应如下式所示,The CFR result obtained by the LS channel estimation algorithm is processed by a denoising algorithm based on DFT transform. The estimation result of the LS algorithm is first subjected to IDFT transform to obtain the channel impulse response As shown in the following formula,
在信道脉冲响应中保留有效抽头信息,将其余部分的噪声置零,得到去噪后的信道脉冲响应如下式所示,The channel impulse response The effective tap information is retained, and the noise of the rest is set to zero to obtain the denoised channel impulse response As shown in the following formula,
再将进行DFT变换,得到第i个OFDM符号基于DFT信道估计算法的CFR估计结果 Then Perform DFT transformation to obtain the CFR estimation result of the i-th OFDM symbol based on the DFT channel estimation algorithm
步骤六,重复步骤三到步骤五,由于LDPC译码对数据SNR的放大作用以及DFT去噪算法对噪声的抑制作用,每次迭代的信道估计结果CFR将逐渐精确,在接收信噪比足够大的情况下达到和非盲估计方法几乎相同的估计精确度和误码性能。Step 6: Repeat steps 3 to 5. Due to the amplification effect of LDPC decoding on data SNR and the noise suppression effect of DFT denoising algorithm, the channel estimation result CFR of each iteration will gradually become more accurate. When the received signal-to-noise ratio is large enough, it can achieve almost the same estimation accuracy and bit error performance as the non-blind estimation method.
综上六个步骤,如图2所示,本实施例还公开基于信道编码辅助的盲信道估计系统,用于实现所述基于信道编码辅助的盲信道估计方法。所述基于信道编码辅助的盲信道估计系统主要由信道估计模块、信道均衡模块、数字解调模块、信道译码模块、硬判决模块和数字调制模块组成。信道估计模块包含LS估计模块和基于DFT变换的去噪模块。信道估计模块涉及步骤二和步骤五,其中基于DFT变换的去噪模块降低信道估计结果的噪声方差,提升信道估计性能,优化迭代结构的收敛效果。信道均衡模块和数字解调模块涉及步骤三,信道均衡模块根据信道估计的结果对接收信号进行补偿,而数字解调模块对补偿后的信号进行星座图解映射得到比特软信息,为信道译码模块提供输入。信道译码模块涉及步骤三,通过信道译码提升通信系统的抗噪声性能,并且保证迭代结构的收敛性。硬判决模块和数字调制模块涉及步骤四,硬判决模块将信道译码输出的软信息结果判决为0、1比特,数字调制模块将硬判决输出的结果重新调制为数字调制符号,为信道估计模块提供等效导频符号IPS。In summary, as shown in FIG2, the present embodiment also discloses a blind channel estimation system based on channel coding assistance, which is used to implement the blind channel estimation method based on channel coding assistance. The blind channel estimation system based on channel coding assistance is mainly composed of a channel estimation module, a channel equalization module, a digital demodulation module, a channel decoding module, a hard decision module and a digital modulation module. The channel estimation module includes an LS estimation module and a denoising module based on DFT transformation. The channel estimation module involves
即图2虚线箭头涉及的迭代流程。随着迭代次数的增加,盲信道估计结果的精确度逐渐向非盲信道估计的性能收敛,使用归一化均方误差(NMSE)作为精确度的衡量标准,如图5所示;使用盲信道估计结果进行信道均衡的误码率性能也逐渐向非盲信道均衡的性能收敛,如图4所示;并且此迭代式盲信道估计方法的收敛速度很快,使用第二次迭代的信道估计结果进行均衡的误码率性能曲线已经几乎与非盲信道均衡的误码率性能曲线相重合,如图4所示。从图4和图5中LDPC-ICE-LS的性能曲线可以看出,使用LS算法进行迭代在Eb/N0>10dB时,误码性能和估计准确度性能也能收敛到和非盲LS信道估计方法相同的性能,但一方面,其性能与使用基于DFT变换的去噪处理的迭代式盲信道估计方法有很大的差距;另一方面,在Eb/N0<10dB时,其误码性能和估计准确度性能是随着迭代次数增加而发散的,无法收敛。这说明在迭代过程中加入基于DFT变换的去噪处理的必要性,提升估计方法的性能,并且优化迭代收敛的效果。That is, the iterative process involved in the dotted arrow in Figure 2. As the number of iterations increases, the blind channel estimation results The accuracy of the channel estimation gradually converges to the performance of non-blind channel estimation, using the normalized mean square error (NMSE) as the accuracy measure, as shown in Figure 5; using the blind channel estimation results The bit error rate performance of channel equalization gradually converges to the performance of non-blind channel equalization, as shown in Figure 4; and the convergence speed of this iterative blind channel estimation method is very fast. The bit error rate performance curve of equalization using the channel estimation result of the second iteration has almost overlapped with the bit error rate performance curve of non-blind channel equalization, as shown in Figure 4. From the performance curves of LDPC-ICE-LS in Figures 4 and 5, it can be seen that when E b /N 0 >10dB is used for iteration using the LS algorithm, the bit error performance and estimation accuracy performance can also converge to the same performance as the non-blind LS channel estimation method, but on the one hand, its performance is far behind that of the iterative blind channel estimation method using DFT-based denoising processing; on the other hand, when E b /N 0 <10dB, its bit error performance and estimation accuracy performance diverge with the increase in the number of iterations and cannot converge. This shows the necessity of adding DFT-based denoising processing in the iterative process to improve the performance of the estimation method and optimize the effect of iterative convergence.
(图中LDPC-ICE-LS表示基于LS算法的LDPC-ICE方法,LDPC-ICE-DFT表示基于DFT算法的LDPC-ICE方法,Pilot-assisted LS表示导频辅助的LS信道估计方法)。(In the figure, LDPC-ICE-LS represents the LDPC-ICE method based on the LS algorithm, LDPC-ICE-DFT represents the LDPC-ICE method based on the DFT algorithm, and Pilot-assisted LS represents the pilot-assisted LS channel estimation method).
以上所述的基于信道编码辅助的盲信道估计方法的搭建为本位发明的实施例,本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或者修改,均落入本发明保护范围。The above-mentioned channel coding-assisted blind channel estimation method is an embodiment of the present invention, and the present invention should not be limited to the contents disclosed in the embodiment and the accompanying drawings. Any equivalent or modification completed without departing from the spirit disclosed in the present invention shall fall within the protection scope of the present invention.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111489753.XA CN114244670B (en) | 2021-12-08 | 2021-12-08 | Blind channel estimation method and system based on channel coding assistance |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111489753.XA CN114244670B (en) | 2021-12-08 | 2021-12-08 | Blind channel estimation method and system based on channel coding assistance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114244670A CN114244670A (en) | 2022-03-25 |
CN114244670B true CN114244670B (en) | 2023-04-18 |
Family
ID=80753844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111489753.XA Active CN114244670B (en) | 2021-12-08 | 2021-12-08 | Blind channel estimation method and system based on channel coding assistance |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114244670B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114389921B (en) * | 2022-01-25 | 2023-12-26 | 山东大学 | Channel estimation method and system based on comb pilot frequency assistance |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107666451A (en) * | 2017-09-15 | 2018-02-06 | 电子科技大学 | Channel estimation methods for LTE system |
CN108156101A (en) * | 2017-12-18 | 2018-06-12 | 中国人民解放军空军工程大学 | A kind of system combined iterative channel estimation of MIMO-SCFDE and iteration equalizing method |
CN108712353A (en) * | 2018-03-29 | 2018-10-26 | 江苏中科羿链通信技术有限公司 | Soft iterative channel estimation method |
CN113395221A (en) * | 2021-04-25 | 2021-09-14 | 北京邮电大学 | Orthogonal time-frequency-space joint-based channel estimation and symbol detection method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100531161C (en) * | 2005-11-07 | 2009-08-19 | 中国人民解放军理工大学 | Associated semi-blind channel estimating and data detecting method based on superimposed pilot and its device |
EP2071786B1 (en) * | 2007-12-14 | 2020-12-23 | Vodafone Holding GmbH | Method and transceiver for data communication |
US8457258B2 (en) * | 2008-02-25 | 2013-06-04 | Indian Institute Of Technology | Optimal training sequence and channel estimation method and system for superimposed training based OFDM systems |
CN102571666B (en) * | 2011-08-12 | 2015-02-25 | 哈尔滨工程大学 | MMSE (Minimum Mean Square Error)-based equalization method of underwater sound OFDM (Orthogonal Frequency Division Multiplexing) judgment iterative channel |
CN104767587B (en) * | 2015-03-10 | 2018-04-27 | 重庆邮电大学 | Based on the compressed sensing channel estimation methods for combining channel decoding under ofdm system |
-
2021
- 2021-12-08 CN CN202111489753.XA patent/CN114244670B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107666451A (en) * | 2017-09-15 | 2018-02-06 | 电子科技大学 | Channel estimation methods for LTE system |
CN108156101A (en) * | 2017-12-18 | 2018-06-12 | 中国人民解放军空军工程大学 | A kind of system combined iterative channel estimation of MIMO-SCFDE and iteration equalizing method |
CN108712353A (en) * | 2018-03-29 | 2018-10-26 | 江苏中科羿链通信技术有限公司 | Soft iterative channel estimation method |
CN113395221A (en) * | 2021-04-25 | 2021-09-14 | 北京邮电大学 | Orthogonal time-frequency-space joint-based channel estimation and symbol detection method |
Non-Patent Citations (2)
Title |
---|
Xuhui Ding.Digital-Analog Hybrid Equalization of Broadband Signals Based on Equivalent Time Sampling.《2020 International Wireless Communications and Mobile Computing (IWCMC)》.2020,全文. * |
白大斌 ; 陈建春 ; .基于编码反馈的OFDM信道估计方法.电子科技.2010,(10),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN114244670A (en) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113381951B (en) | MFTN Joint Channel Estimation and Equalization Method in Time-variable Frequency Selective Fading Channel | |
WO2011082621A1 (en) | Detection method and device for multi-input multi-output single carrier block transmission system | |
CN108712353B (en) | Soft iteration channel estimation method | |
CN110677359A (en) | Signal receiving method, receiving device and storage medium of orthogonal time-frequency space system | |
CN101132388A (en) | Receiving method and device for receiving coded signal assisted by signal channel condition information | |
WO2011111583A1 (en) | Receiving device, receiving method, receiving program, and processor | |
CN110311872A (en) | Method, system, device and storage medium for underwater acoustic sparse channel estimation | |
CN114401172B (en) | Combined estimation and detection method based on Turbo equalization frame and VAMP | |
CN110971307A (en) | A method for compensating for inherent inter-subcarrier interference in SEFDM system | |
CN114244670B (en) | Blind channel estimation method and system based on channel coding assistance | |
CN110324271B (en) | Amplitude limiting F-OFDM system transceiver design method based on compressed sensing | |
CN113411277B (en) | Construction method of multi-carrier super-Nyquist system receiver | |
CN113381954B (en) | Frequency domain equalization method based on generalized approximate message transmission | |
Ketonen et al. | MIMO detection in single carrier systems | |
CN109639301B (en) | A FTN Equalization Method Based on Confidence Estimation | |
CN104135455B (en) | Iterative receiving method for communication system | |
JP2019501582A (en) | Method and system for demodulating higher-order QAM signals | |
CN113556305A (en) | FBMC iterative channel equalization method and system suitable for high frequency selective fading | |
CN101945062B (en) | Receiving processing method in single-carrier wireless communication system | |
CN114978843B (en) | A time-varying channel tracking method for OFDM systems based on decoding assistance | |
CN111884959B (en) | An iterative equalization method and device for a single-carrier broadband wireless communication system | |
CN118089741B (en) | Navigation data processing method based on turbo equalization in delay-Doppler domain | |
CN114726693B (en) | An Iterative Equalization Method for DFT-s-SEFDM Signals | |
CN116961736B (en) | Uplink communication method for low-orbit satellite terminal with limited power | |
CN107634926B (en) | A High-Order Modulation Soft Decision Method for Joint Channel Information |
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 |