Nothing Special   »   [go: up one dir, main page]

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 PDF

Info

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
Application number
CN202111489753.XA
Other languages
Chinese (zh)
Other versions
CN114244670A (en
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202111489753.XA priority Critical patent/CN114244670B/en
Publication of CN114244670A publication Critical patent/CN114244670A/en
Application granted granted Critical
Publication of CN114244670B publication Critical patent/CN114244670B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0238Channel estimation using blind estimation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2626Arrangements specific to the transmitter only
    • H04L27/2627Modulators
    • H04L27/2634Inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators in combination with other circuits for modulation
    • H04L27/2636Inverse 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]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing 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

The invention discloses a blind channel estimation method and a blind channel estimation system based on channel coding assistance, and belongs to the field of communication signal processing. The invention feeds back the hard decision bit information output by the decoder to the channel estimation module in an iterative way, and takes the hard decision bit information as prior information of sending data to replace the function of the original pilot frequency symbol, thereby leading the invention to have the following main advantages: pilot symbols do not need to be inserted at a transmitting end, so that the problem of frequency spectrum utilization rate reduction caused by the insertion of the pilot symbols is avoided, and the throughput of the OFDM system is remarkably improved; on the basis of the non-blind channel estimation method, the blind channel estimation is realized by taking the hard decision bit information as the prior information of the transmitted data to replace the function of the original pilot frequency symbol, so that the method has the estimation complexity and the estimation precision which are similar to those of the non-blind channel estimation method, and can achieve the error rate performance which is approximate to the theoretical limit. The invention adds a feedback loop on the basis of the OFDM receiver, has simple structure and is easy to transplant.

Description

一种基于信道编码辅助的盲信道估计方法及系统A blind channel estimation method and system based on channel coding assistance

技术领域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 step 1, extract the pilot symbol in the first OFDM symbol, and perform frequency domain channel estimation.

以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:

Figure BDA0003398750110000021
Figure BDA0003398750110000021

去除循环前缀后,接收OFDM符号的有效序列如下式所示,After removing the cyclic prefix, the effective sequence of the received OFDM symbol is as follows,

Figure BDA0003398750110000022
Figure BDA0003398750110000022

其中

Figure BDA00033987501100000211
表示循环卷积,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
Figure BDA00033987501100000211
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

Figure BDA0003398750110000023
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 ,
Figure BDA0003398750110000023
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,

Figure BDA0003398750110000024
Figure BDA0003398750110000024

当X(k)已知时通过LS估计算法得出H(k)的估计值。记Xp(k)为导频符号,Xp=[Xp(0),...,Xp(Nsc-1)]T

Figure BDA0003398750110000025
H的LS估计值
Figure BDA0003398750110000026
如下式所示,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 ,
Figure BDA0003398750110000025
LS estimate of H
Figure BDA0003398750110000026
As shown in the following formula,

Figure BDA0003398750110000027
Figure BDA0003398750110000027

根据上式得,第一个OFDM符号LS信道估计结果

Figure BDA0003398750110000028
中各个元素
Figure BDA0003398750110000029
的表达式为,According to the above formula, the LS channel estimation result of the first OFDM symbol is
Figure BDA0003398750110000028
Each element
Figure BDA0003398750110000029
The expression of is,

Figure BDA00033987501100000210
Figure BDA00033987501100000210

角标中的(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符号开始信号中不再插入导频符号,信道均衡后的结果

Figure BDA0003398750110000031
表达式为,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
Figure BDA0003398750110000031
The expression is,

Figure BDA0003398750110000032
Figure BDA0003398750110000032

之后进行后续的数字解调过程,得到各个比特的软信息

Figure BDA0003398750110000033
m=0,...,M-1,M为
Figure BDA0003398750110000034
数字解调后的比特数,将比特软信息传递给LDPC译码器,进行LDPC信道译码,得到包括校验位在内所有变量节点的译码输出结果
Figure BDA0003398750110000035
m=0,...,M-1。After that, the subsequent digital demodulation process is carried out to obtain the soft information of each bit.
Figure BDA0003398750110000033
m=0,...,M-1, M is
Figure BDA0003398750110000034
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
Figure BDA0003398750110000035
m=0,...,M-1.

步骤四,将步骤三译码器变量节点的输出结果进行硬判决得到比特信息,将所述比特信息重新进行数字调制,得到的调制符号后文称为迭代导频符号(Iterative PilotSymbols,IPS)。由于译码器带来的编码增益,此IPS符号相比步骤三信道均衡的结果准确度显著提高,使用此等效的导频符号辅助盲信道估计。由于信号内部实际并未插入导频符号,本发明所述通信系统相比于采用非盲信道估计的通信系统提高频谱利用率。Step 4, hard decision is performed on the output result of the decoder variable node in step 3 to obtain bit information, and the bit information is re-digitally modulated, and the obtained modulation symbol is hereinafter referred to as iterative pilot symbol (IPS). Due to the coding gain brought by the decoder, the accuracy of the result of the IPS symbol is significantly improved compared with the channel equalization in step 3, and this equivalent pilot symbol is used to assist blind channel estimation. Since the pilot symbol is not actually inserted into the signal, the communication system of the present invention improves the spectrum utilization compared with the communication system using non-blind channel estimation.

对LDPC译码输出结果

Figure BDA0003398750110000036
进行硬判决得到序列b(i)(m),m=0,...,M-1,此结果为0、1比特形式,将译码的硬判决输出结果重新进行数字调制,其调制方式应与信号原本数字调制方式相同,调制结果看作等效的导频符号IPS,记为
Figure BDA0003398750110000037
i=2,3,...。Output results of LDPC decoding
Figure BDA0003398750110000036
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
Figure BDA0003398750110000037
i=2,3,...

步骤五,把步骤三得到的迭代导频符号IPS视为第i个OFDM符号的先验信息,以IPS为导频符号,基于非盲信道估计算法进行第i个OFDM符号的信道估计,实现以非盲估计算法进行盲信道估计的目的,相比于传统盲信道估计降低实现复杂度且提高估计精度。在估计算法流程上,首先对信号进行基于LS算法的信道估计,再对LS算法的估计结果进行基于DFT变换的去噪处理,得到DFT信道估计算法的估计结果。本步骤中进行基于DFT变换的去噪处理有两个目的,其一是提高信道估计的准确性,提升通信系统的误码性能;其二是降低均衡过程中引入的噪声,优化迭代结构的收敛性,保证迭代信道估计结构在较低信噪比下能够收敛。Step 5, regard the iterative pilot symbol IPS obtained in step 3 as the prior information of the i-th OFDM symbol, use IPS as the pilot symbol, and perform channel estimation of the i-th OFDM symbol based on the non-blind channel estimation algorithm, so as to achieve the purpose of blind channel estimation with the non-blind estimation algorithm. Compared with the traditional blind channel estimation, it reduces the implementation complexity and improves the estimation accuracy. In the estimation algorithm process, firstly, the signal is estimated based on the LS algorithm, and then the estimation result of the LS algorithm is denoised based on the DFT transform to obtain the estimation result of the DFT channel estimation algorithm. There are two purposes for the denoising based on the DFT transform in this step. One is to improve the accuracy of channel estimation and improve the error performance of the communication system; the other is to reduce the noise introduced in the equalization process, optimize the convergence of the iterative structure, and ensure that the iterative channel estimation structure can converge at a lower signal-to-noise ratio.

将IPS反馈给信道估计模块,把IPS视作导频符号,基于LS算法估计出第i个OFDM符号的信道CFR估计值

Figure BDA0003398750110000038
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
Figure BDA0003398750110000038

Figure BDA0003398750110000039
Figure BDA0003398750110000039

将LS信道估计算法得到的CFR结果进行基于DFT变换的去噪算法处理。LS算法的估计结果先进行IDFT变换,得到信道脉冲响应

Figure BDA00033987501100000310
如下式所示,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
Figure BDA00033987501100000310
As shown in the following formula,

Figure BDA00033987501100000311
Figure BDA00033987501100000311

在信道脉冲响应

Figure BDA00033987501100000312
中保留有效抽头信息,将其余部分的噪声置零,得到去噪后的信道脉冲响应
Figure BDA0003398750110000041
如下式所示,The channel impulse response
Figure BDA00033987501100000312
The effective tap information is retained, and the noise of the rest is set to zero to obtain the denoised channel impulse response
Figure BDA0003398750110000041
As shown in the following formula,

Figure BDA0003398750110000042
Figure BDA0003398750110000042

再将

Figure BDA0003398750110000043
进行DFT变换,得到第i个OFDM符号基于DFT信道估计算法的CFR估计结果
Figure BDA0003398750110000044
Then
Figure BDA0003398750110000043
Perform DFT transformation to obtain the CFR estimation result of the i-th OFDM symbol based on the DFT channel estimation algorithm
Figure BDA0003398750110000044

Figure BDA0003398750110000045
Figure BDA0003398750110000045

步骤六,重复步骤三到步骤五,由于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.

本发明还公开基于信道编码辅助的盲信道估计系统,用于实现所述基于信道编码辅助的盲信道估计方法。所述基于信道编码辅助的盲信道估计系统主要由信道估计模块、信道均衡模块、数字解调模块、信道译码模块、硬判决模块和数字调制模块组成。信道估计模块包含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 steps 2 and 5, wherein the denoising module based on DFT transformation reduces the noise variance of the channel estimation result, improves the channel estimation performance, and optimizes the convergence effect of the iterative structure. The channel equalization module and the digital demodulation module involve step 3, the channel equalization module compensates the received signal according to the result of the channel estimation, and the digital demodulation module performs constellation diagram demapping on the compensated signal to obtain bit soft information, which provides input for the channel decoding module. The channel decoding module involves step 3, which improves the anti-noise performance of the communication system through channel decoding, and ensures the convergence of the iterative structure. The hard decision module and the digital modulation module involve step four. The hard decision module judges the soft information result output by the channel decoding as 0 or 1 bit, and the digital modulation module remodulates the result output by the hard decision into a digital modulation symbol to provide an equivalent pilot symbol IPS for the channel estimation module.

作为优选,在信道编码过程中,使用纠错性能优异的低密度奇偶校验(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:

Figure BDA0003398750110000051
Figure BDA0003398750110000051

本实施例所提出的盲估计方法整体流程图如图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 step 1, extract the pilot symbol in the first OFDM symbol, and perform frequency domain channel estimation.

以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:

Figure BDA0003398750110000061
Figure BDA0003398750110000061

去除循环前缀后,接收OFDM符号的有效序列如下式表示,After removing the cyclic prefix, the effective sequence of the received OFDM symbol is expressed as follows:

Figure BDA0003398750110000062
Figure BDA0003398750110000062

其中

Figure BDA0003398750110000063
表示循环卷积,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
Figure BDA0003398750110000063
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

Figure BDA0003398750110000064
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 ,
Figure BDA0003398750110000064
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:

Figure BDA0003398750110000065
Figure BDA0003398750110000065

当X(k)已知时通过LS估计算法得出H(k)的估计值。记Xp(k)为导频符号,Xp=[Xp(0),...,Xp(959)]T

Figure BDA0003398750110000066
H的LS估计值
Figure BDA0003398750110000067
如下式所示,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 ,
Figure BDA0003398750110000066
LS estimate of H
Figure BDA0003398750110000067
As shown in the following formula,

Figure BDA0003398750110000068
Figure BDA0003398750110000068

根据上式得,第一个OFDM符号LS信道估计结果

Figure BDA0003398750110000069
中各个元素
Figure BDA00033987501100000610
的表达式为,According to the above formula, the LS channel estimation result of the first OFDM symbol is
Figure BDA0003398750110000069
Each element
Figure BDA00033987501100000610
The expression of is,

Figure BDA00033987501100000611
Figure BDA00033987501100000611

角标中的(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符号开始信号中不再插入导频符号,信道均衡后的结果

Figure BDA0003398750110000071
表达式为,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
Figure BDA0003398750110000071
The expression is,

Figure BDA0003398750110000072
Figure BDA0003398750110000072

之后进行后续的数字解调过程,得到各个比特的软信息

Figure BDA0003398750110000073
m=0,...,1919,
Figure BDA0003398750110000074
数字解调后的比特数为1920,将比特软信息传递给LDPC译码器,进行LDPC信道译码,得到包括校验位在内所有变量节点的译码输出结果
Figure BDA0003398750110000075
m=0,...,1919。After that, the subsequent digital demodulation process is carried out to obtain the soft information of each bit.
Figure BDA0003398750110000073
m=0,...,1919,
Figure BDA0003398750110000074
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.
Figure BDA0003398750110000075
m=0,...,1919.

步骤四,将步骤三译码器变量节点的输出结果进行硬判决得到比特信息,将所述比特信息重新进行数字调制,得到的调制符号后文称为迭代导频符号(Iterative PilotSymbols,IPS)。由于译码器带来的编码增益,此IPS符号相比步骤三信道均衡的结果准确度显著提高,使用此等效的导频符号辅助盲信道估计。由于信号内部实际并未插入导频符号,本发明所述通信系统相比于采用非盲信道估计的通信系统提高了频谱利用率。Step 4, hard decision is performed on the output result of the decoder variable node in step 3 to obtain bit information, and the bit information is re-digitally modulated, and the obtained modulation symbol is hereinafter referred to as iterative pilot symbol (IPS). Due to the coding gain brought by the decoder, the accuracy of the result of the IPS symbol is significantly improved compared with the channel equalization in step 3, and this equivalent pilot symbol is used to assist blind channel estimation. Since the pilot symbol is not actually inserted into the signal, the communication system of the present invention improves the spectrum utilization compared with the communication system using non-blind channel estimation.

对LDPC译码输出结果

Figure BDA0003398750110000076
进行硬判决得到序列b(i)(m),m=0,...,1919,此结果为0、1比特形式,将译码的硬判决输出结果重新进行数字调制,其调制方式应与信号原本数字调制方式相同,调制结果看作等效的导频符号IPS,记为
Figure BDA0003398750110000077
i=2,3,...120。Output results of LDPC decoding
Figure BDA0003398750110000076
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
Figure BDA0003398750110000077
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 step 3 is regarded as the prior information of the i-th OFDM symbol. With IPS as the pilot symbol, the channel estimation of the i-th OFDM symbol is performed based on the non-blind channel estimation algorithm, so as to achieve the purpose of blind channel estimation with the non-blind estimation algorithm. Compared with the traditional blind channel estimation, the implementation complexity is reduced and the estimation accuracy is improved. In the estimation algorithm process, the signal is firstly subjected to channel estimation based on the LS algorithm, and then the estimation result of the LS algorithm is subjected to denoising based on the DFT transform to obtain the estimation result of the DFT channel estimation algorithm. The schematic diagram of the denoising module based on the DFT transform is shown in FIG3. There are two purposes for the denoising based on the DFT transform in this step. One is to improve the accuracy of the channel estimation and improve the error performance of the communication system; the other is to reduce the noise introduced in the equalization process, optimize the convergence of the iterative structure, and ensure that the iterative channel estimation structure can converge at a lower signal-to-noise ratio, as shown in FIG4 and FIG5.

将IPS反馈给信道估计模块,把IPS视作导频符号,基于LS算法估计出第i个OFDM符号的信道CFR估计值

Figure BDA0003398750110000078
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
Figure BDA0003398750110000078

Figure BDA0003398750110000079
Figure BDA0003398750110000079

将LS信道估计算法得到的CFR结果进行基于DFT变换的去噪算法处理。LS算法的估计结果先进行IDFT变换,得到信道脉冲响应

Figure BDA00033987501100000710
如下式所示,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
Figure BDA00033987501100000710
As shown in the following formula,

Figure BDA00033987501100000711
Figure BDA00033987501100000711

在信道脉冲响应

Figure BDA0003398750110000081
中保留有效抽头信息,将其余部分的噪声置零,得到去噪后的信道脉冲响应
Figure BDA0003398750110000082
如下式所示,The channel impulse response
Figure BDA0003398750110000081
The effective tap information is retained, and the noise of the rest is set to zero to obtain the denoised channel impulse response
Figure BDA0003398750110000082
As shown in the following formula,

Figure BDA0003398750110000083
Figure BDA0003398750110000083

再将

Figure BDA0003398750110000084
进行DFT变换,得到第i个OFDM符号基于DFT信道估计算法的CFR估计结果
Figure BDA0003398750110000085
Then
Figure BDA0003398750110000084
Perform DFT transformation to obtain the CFR estimation result of the i-th OFDM symbol based on the DFT channel estimation algorithm
Figure BDA0003398750110000085

Figure BDA0003398750110000086
Figure BDA0003398750110000086

步骤六,重复步骤三到步骤五,由于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 steps 2 and 5, wherein the denoising module based on DFT transformation reduces the noise variance of the channel estimation result, improves the channel estimation performance, and optimizes the convergence effect of the iterative structure. The channel equalization module and the digital demodulation module involve step 3, the channel equalization module compensates the received signal according to the result of the channel estimation, and the digital demodulation module performs constellation diagram demapping on the compensated signal to obtain bit soft information, which provides input for the channel decoding module. The channel decoding module involves step 3, which improves the anti-noise performance of the communication system through channel decoding and ensures the convergence of the iterative structure. The hard decision module and the digital modulation module involve step four. The hard decision module judges the soft information result output by the channel decoding as 0 or 1 bit, and the digital modulation module remodulates the result output by the hard decision into a digital modulation symbol to provide an equivalent pilot symbol IPS for the channel estimation module.

即图2虚线箭头涉及的迭代流程。随着迭代次数的增加,盲信道估计结果

Figure BDA0003398750110000087
的精确度逐渐向非盲信道估计的性能收敛,使用归一化均方误差(NMSE)作为精确度的衡量标准,如图5所示;使用盲信道估计结果
Figure BDA0003398750110000088
进行信道均衡的误码率性能也逐渐向非盲信道均衡的性能收敛,如图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
Figure BDA0003398750110000087
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
Figure BDA0003398750110000088
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)

1.一种基于信道编码辅助的盲信道估计方法,其特征在于:包括如下步骤,1. A blind channel estimation method based on channel coding assistance, characterized in that it comprises the following steps: 步骤一,发送具有如下帧格式的信号;Step 1, sending a signal having the following frame format; 所述帧格式为:长度为Nsym个符号的OFDM帧中,第一个符号插入导频信息,后续(Nsym-1)个OFDM符号无需插入额外的导频符号;The frame format is: in an OFDM frame with a length of N sym symbols, the first symbol is inserted with pilot information, and no additional pilot symbols are inserted into the subsequent (N sym -1) OFDM symbols; 步骤二,接收步骤一发送的信号,提取第一个OFDM符号中的导频符号,进行频域信道估计;Step 2, receiving the signal sent in step 1, extracting the pilot symbol in the first OFDM symbol, and performing frequency domain channel estimation; 步骤三,使用第一个OFDM符号的信道估计结果,对第二个OFDM符号进行信道均衡,均衡的结果进行后续数字解调和信道译码,得到译码器变量节点输出结果;步骤三的流程不只限于对第二个OFDM符号进行处理,从第二个OFDM符号起,所有OFDM符号进行信道估计时均需要进行相同操作,即从步骤三开始迭代;Step 3: Use the channel estimation result of the first OFDM symbol to perform channel equalization on the second OFDM symbol, and perform subsequent digital demodulation and channel decoding on the equalization result to obtain the output result of the decoder variable node; the process of step 3 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, iterate from step 3; 步骤三实现方法为,The implementation method of step three is: 使用第i-1个OFDM符号的信道估计结果对第i个OFDM符号进行信道均衡,从第二个OFDM符号开始信号中不再插入导频符号,信道均衡后的结果
Figure FDA0004070777840000011
表达式为,
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
Figure FDA0004070777840000011
The expression is,
Figure FDA0004070777840000012
Figure FDA0004070777840000012
之后进行后续的数字解调过程,得到各个比特的软信息
Figure FDA0004070777840000013
M为
Figure FDA0004070777840000014
数字解调后的比特数,将比特软信息传递给LDPC译码器,进行LDPC信道译码,得到包括校验位在内所有变量节点的译码输出结果
Figure FDA0004070777840000015
After that, the subsequent digital demodulation process is carried out to obtain the soft information of each bit.
Figure FDA0004070777840000013
M is
Figure FDA0004070777840000014
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.
Figure FDA0004070777840000015
步骤四,将步骤三译码器变量节点的输出结果进行硬判决得到比特信息,将所述比特信息重新进行数字调制,得到的调制符号为迭代导频符号(Iterative Pilot Symbols,IPS);由于译码器带来的编码增益,此IPS符号相比步骤三信道均衡的结果准确度显著提高,使用此等效的导频符号辅助盲信道估计,提高频谱利用率;Step 4: hard decision is performed on the output result of the decoder variable node in step 3 to obtain bit information, and the bit information is re-digitally modulated to obtain a modulation symbol, which is an iterative pilot symbol (IPS). Due to the coding gain brought by the decoder, the accuracy of the IPS symbol is significantly improved compared with the result of the channel equalization in step 3. This equivalent pilot symbol is used to assist blind channel estimation and improve spectrum utilization. 步骤四中,对LDPC译码输出结果
Figure FDA0004070777840000016
进行硬判决得到序列b(i)(m),m=0,…,M-1,此结果为0、1比特形式,将译码的硬判决输出结果重新进行数字调制,其调制方式应与信号原本数字调制方式相同,调制结果为等效的导频符号IPS,记为
Figure FDA0004070777840000017
In step 4, the LDPC decoding output result
Figure FDA0004070777840000016
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 an equivalent pilot symbol IPS, which is recorded as
Figure FDA0004070777840000017
步骤五,把步骤三得到的迭代导频符号IPS为第i个OFDM符号的先验信息,以IPS为导频符号,基于非盲信道估计算法进行第i个OFDM符号的信道估计,实现以非盲信道估计算法进行盲信道估计的目的,相比于传统盲信道估计降低实现复杂度且提高估计精度;在信道估计算法流程上,首先对信号进行基于LS信道估计算法的信道估计,再对LS信道估计算法的估计结果进行基于DFT变换的去噪处理,得到DFT信道估计算法的估计结果;本步骤中进行基于DFT变换的去噪处理有两个目的,其一是提高信道估计的准确性,提升通信系统的误码性能;其二是降低均衡过程中引入的噪声,优化迭代结构的收敛性;Step 5: The iterative pilot symbol IPS obtained in step 3 is used as the prior information of the i-th OFDM symbol, and IPS is used as the pilot symbol. Channel estimation of the i-th OFDM symbol is performed based on a non-blind channel estimation algorithm, so as to achieve the purpose of blind channel estimation using a non-blind channel estimation algorithm. Compared with the traditional blind channel estimation, the implementation complexity is reduced and the estimation accuracy is improved. In the channel estimation algorithm process, the signal is firstly subjected to channel estimation based on the LS channel estimation algorithm, and then the estimation result of the LS channel estimation algorithm is subjected to denoising processing based on the DFT transform to obtain the estimation result of the DFT channel estimation algorithm. The denoising processing based on the DFT transform in this step has two purposes: one is to improve the accuracy of channel estimation and enhance the error performance of the communication system; the other is to reduce the noise introduced in the equalization process and optimize the convergence of the iterative structure. 步骤五中,In step five, 将IPS反馈给信道估计模块,把IPS为导频符号,基于LS信道估计算法估计出第i个OFDM符号的信道估计结果
Figure FDA0004070777840000021
Feedback IPS to the channel estimation module, take IPS as the pilot symbol, and estimate the channel estimation result of the i-th OFDM symbol based on the LS channel estimation algorithm
Figure FDA0004070777840000021
Figure FDA0004070777840000022
Figure FDA0004070777840000022
将LS信道估计算法得到的信道估计结果进行基于DFT变换的DFT去噪算法处理;LS信道估计算法的估计结果先进行IDFT变换,得到信道脉冲响应
Figure FDA0004070777840000023
如下式所示,
The channel estimation result obtained by the LS channel estimation algorithm is processed by the DFT denoising algorithm based on DFT transformation; the estimation result of the LS channel estimation algorithm is first subjected to IDFT transformation to obtain the channel impulse response
Figure FDA0004070777840000023
As shown in the following formula,
Figure FDA0004070777840000024
Figure FDA0004070777840000024
在信道脉冲响应
Figure FDA0004070777840000025
中保留有效抽头信息,将其余部分的噪声置零,得到去噪后的信道脉冲响应
Figure FDA0004070777840000026
如下式所示,
The channel impulse response
Figure FDA0004070777840000025
The effective tap information is retained, and the noise of the rest is set to zero to obtain the denoised channel impulse response
Figure FDA0004070777840000026
As shown in the following formula,
Figure FDA0004070777840000027
Figure FDA0004070777840000027
再将
Figure FDA0004070777840000028
进行DFT变换,得到第i个OFDM符号基于DFT信道估计算法的信道估计结果
Figure FDA0004070777840000029
Then
Figure FDA0004070777840000028
Perform DFT transformation to obtain the channel estimation result of the i-th OFDM symbol based on the DFT channel estimation algorithm
Figure FDA0004070777840000029
Figure FDA00040707778400000210
Figure FDA00040707778400000210
步骤六,重复步骤三到步骤五,由于LDPC译码对数据SNR的放大作用以及DFT去噪算法对噪声的抑制作用,每次迭代的信道估计结果将逐渐精确,在接收信噪比足够大的情况下达到和非盲信道估计算法相同的估计精确度和误码性能。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 of each iteration will gradually become more accurate, and when the received signal-to-noise ratio is large enough, it will achieve the same estimation accuracy and error performance as the non-blind channel estimation algorithm.
2.如权利要求1所述的一种基于信道编码辅助的盲信道估计方法,其特征在于:步骤二实现方法为,2. A blind channel estimation method based on channel coding assistance as claimed in claim 1, characterized in that: the implementation method of step 2 is: 以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:
Figure FDA00040707778400000211
Figure FDA00040707778400000211
去除循环前缀后,接收OFDM符号的有效序列如下式所示,After removing the cyclic prefix, the effective sequence of the received OFDM symbol is as follows,
Figure FDA00040707778400000212
Figure FDA00040707778400000212
其中
Figure FDA00040707778400000213
表示循环卷积,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
Figure FDA00040707778400000213
represents circular convolution, N is the number of FFT points in the OFDM system; N sc is the number of OFDM subcarriers, 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,…,Nsc-1 (3)Y(k)=X(k)·H(k)+W(k),k=0,…,N sc -1 (3)
Figure FDA0004070777840000039
Figure FDA00040707778400000310
接收数据频域表达式的矩阵表示为,
remember
Figure FDA0004070777840000039
Figure FDA00040707778400000310
The matrix representation of the received data frequency domain expression is:
Figure FDA0004070777840000032
Figure FDA0004070777840000032
当X(k)已知时通过LS信道估计算法得出H(k)的估计值;记Xp(k)为导频符号,
Figure FDA00040707778400000311
H的LS信道估计结果
Figure FDA0004070777840000034
如下式所示,
When X(k) is known, the estimated value of H(k) is obtained by the LS channel estimation algorithm; Xp (k) is denoted as the pilot symbol,
Figure FDA00040707778400000311
LS channel estimation results of H
Figure FDA0004070777840000034
As shown in the following formula,
Figure FDA0004070777840000035
Figure FDA0004070777840000035
根据上式得,第一个OFDM符号LS信道估计结果
Figure FDA0004070777840000036
中各个元素
Figure FDA0004070777840000037
的表达式为,
According to the above formula, the LS channel estimation result of the first OFDM symbol is
Figure FDA0004070777840000036
Each element
Figure FDA0004070777840000037
The expression of is,
Figure FDA0004070777840000038
Figure FDA0004070777840000038
角标中的(i)表示第i个OFDM符号。The (i) in the superscript represents the i-th OFDM symbol.
3.基于信道编码辅助的盲信道估计系统,用于实现如权利要求1或2所述的一种基于信道编码辅助的盲信道估计方法,其特征在于:主要由信道估计模块、信道均衡模块、数字解调模块、信道译码模块、硬判决模块和数字调制模块组成;信道估计模块包含LS估计模块和基于DFT变换的去噪模块;信道估计模块涉及步骤二和步骤五,其中基于DFT变换的去噪模块降低信道估计结果的噪声方差,提升信道估计性能,优化迭代结构的收敛效果;信道均衡模块和数字解调模块涉及步骤三,信道均衡模块根据信道估计的结果对接收信号进行补偿,而数字解调模块对补偿后的信号进行星座图解映射得到比特软信息,为信道译码模块提供输入;信道译码模块涉及步骤三,通过信道译码提升通信系统的抗噪声性能,并且保证迭代结构的收敛性;硬判决模块和数字调制模块涉及步骤四,硬判决模块将信道译码输出的软信息结果判决为0、1比特,数字调制模块将硬判决输出的结果重新调制为数字调制符号,为信道估计模块提供等效导频符号IPS。3. A blind channel estimation system based on channel coding assistance, used to implement a blind channel estimation method based on channel coding assistance as described in claim 1 or 2, characterized in that: it 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 steps 2 and 5, wherein the denoising module based on DFT transformation reduces the noise variance of the channel estimation result, improves the channel estimation performance, and optimizes the convergence effect of the iterative structure; the channel equalization module and the digital demodulation module The module involves step three, the channel equalization module compensates the received signal according to the result of channel estimation, and the digital demodulation module performs constellation diagram de-mapping on the compensated signal to obtain bit soft information, which provides input for the channel decoding module; the channel decoding module involves step three, which improves the anti-noise performance of the communication system through channel decoding and ensures the convergence of the iterative structure; the hard decision module and the digital modulation module involve step four, the hard decision module judges the soft information result of the channel decoding output as 0 or 1 bit, and the digital modulation module re-modulates the result of the hard decision output into a digital modulation symbol, which provides an equivalent pilot symbol IPS for the channel estimation module. 4.如权利要求3所述的基于信道编码辅助的盲信道估计系统,其特征在于:在信道编码过程中,使用纠错性能优异的低密度奇偶校验(Low-Density Parity-Check,LDPC)码,增强信道估计及通信系统的抗噪声性能;而基于LDPC编码辅助的迭代式盲信道估计方法简称为LDPC-ICE(LDPC-code-aided iterative channel estimation)。4. The blind channel estimation system based on channel coding assistance as described in claim 3 is characterized in that: 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; and the iterative blind channel estimation method based on LDPC coding assistance is referred to as LDPC-ICE (LDPC-code-aided iterative channel estimation).
CN202111489753.XA 2021-12-08 2021-12-08 Blind channel estimation method and system based on channel coding assistance Active CN114244670B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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