[go: up one dir, main page]

CN1362789A - Turbo code decoding method and decoder - Google Patents

Turbo code decoding method and decoder Download PDF

Info

Publication number
CN1362789A
CN1362789A CN 01101350 CN01101350A CN1362789A CN 1362789 A CN1362789 A CN 1362789A CN 01101350 CN01101350 CN 01101350 CN 01101350 A CN01101350 A CN 01101350A CN 1362789 A CN1362789 A CN 1362789A
Authority
CN
China
Prior art keywords
state
metric
calculation
decoding
calculation module
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 01101350
Other languages
Chinese (zh)
Other versions
CN1145266C (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.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies Co Ltd
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 Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Priority to CNB011013508A priority Critical patent/CN1145266C/en
Publication of CN1362789A publication Critical patent/CN1362789A/en
Application granted granted Critical
Publication of CN1145266C publication Critical patent/CN1145266C/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Error Detection And Correction (AREA)

Abstract

本发明公开了一种特博码解码方法及其解码器,该方法在最大后验概率解码方法的基础上,增加了对通过状态网格图对初始几个时刻状态度量进行计算控制;该解码器中增设了α和β初始状态度量计算控制模块,对前项状态度量和初始后项状态度量进行控制以消除多余状态度量的引入。本发明与比没有初始状态度量计算控制的解码性能提高大约有0.1db的性能增益。

The invention discloses a turbo code decoding method and a decoder thereof. On the basis of the maximum a posteriori probability decoding method, the method adds the calculation and control of the state metrics at the initial several moments through the state grid graph; the decoding The α and β initial state measure calculation control modules are added in the controller to control the former state measure and the initial subsequent state measure to eliminate the introduction of redundant state measure. The present invention has a performance gain of about 0.1db compared to the decoding performance improvement without initial state metric calculation control.

Description

特博码解码方法及其解码器Turbine Code Decoding Method and Its Decoder

本发明涉及移动通信领域中的信道编解码技术,更具体地指一种可提高解码性能的特博码解码方法及其解码器。The present invention relates to channel coding and decoding technology in the field of mobile communication, more specifically to a turbo code decoding method and a decoder thereof which can improve decoding performance.

在无线通信系统中,由于传输介质的不均匀性和不稳定性,传输的信号会受到时间扩散、衰落等干扰影响,造成接收的比特有随机性的差错。为了防止信道噪声的干扰影响,需要采用一定的方式来提高信息的传送可靠性和有效性。通过增加冗余度来降低误码率的纠错编码方法,被时间证明是一类有效可靠的重要手段。特别在移动通信和卫星通信系统中。纠错码得到广泛的技术应用。In a wireless communication system, due to the inhomogeneity and instability of the transmission medium, the transmitted signal will be affected by interference such as time diffusion and fading, resulting in random errors in the received bits. In order to prevent the interference effect of channel noise, it is necessary to adopt a certain method to improve the reliability and effectiveness of information transmission. The error correction coding method, which reduces the bit error rate by increasing redundancy, has been proved to be an effective and reliable important means by time. Especially in mobile communication and satellite communication systems. Error-correcting codes are widely used in technology.

1993年Berrou等人提出了特博(Turbo,下同)码,它是一种接近伪随机编码的可解码。特博是一种回归系统卷积码(RSC),这种编码方式的优越性来源于其独特的码结构和迭代解码算法。迭代解码突破了最小码距的设计思想,更接近香农随机码的概念。In 1993, Berrou et al. proposed a Turbo (Turbo, the same below) code, which is a decodable code close to a pseudo-random code. Tebo is a regression systematic convolutional code (RSC). The superiority of this coding method comes from its unique code structure and iterative decoding algorithm. Iterative decoding breaks through the design idea of minimum code distance, and is closer to the concept of Shannon random code.

特博码被证明是一种纠错能力很强的码。它的编码器是由两个或多个子编码器通过串联或并联的级联方式构成的,通常较普遍的特博码编码器是由两个卷积码编码器并联而成,输入信息位一路直接送入子编码器1,另一路经过交织器后送入子编码器2编码,编完后的数据再经过打孔器打孔调制到合适的码率输出。如图1所示。Turbine code has been proved to be a code with strong error correction ability. Its encoder is composed of two or more sub-encoders connected in series or in parallel. Usually, the more common turbo code encoder is composed of two convolutional code encoders in parallel, and the input information bits are all the way It is directly sent to the sub-encoder 1, and the other channel is sent to the sub-encoder 2 for encoding after passing through the interleaver. After encoding, the encoded data is then punctured by a puncher and modulated to a suitable code rate for output. As shown in Figure 1.

图1是cdma2000和WCDMA提案中的特博码编码器的结构,交织器的作用是对输入数据的顺序进行重新编排,目的是调整权重的分布,使得子编码器2输入比特流的权重分布与子编码器1的不同。打孔器对两个子编码器输出的六路比特进行打孔抽样和并串转换。Figure 1 shows the structure of the turbo code encoder in the cdma2000 and WCDMA proposals. The function of the interleaver is to rearrange the order of the input data, and the purpose is to adjust the weight distribution so that the weight distribution of the sub-encoder 2 input bit stream is consistent with Sub-encoder 1 is different. The puncturer performs puncture sampling and parallel-to-serial conversion on the six bits output by the two sub-encoders.

卷积码编码器通常用(n0,k0,m)来表示特征,n0是输出比特,k0是输入比特,m是寄存器个数。K表示约束长度,即为卷积码内部移位寄存器的个数m加1。cdma2000提案中特博码子(3,1,4)回归系统卷积码编码器如图2所示。A convolutional code encoder usually uses (n0, k0, m) to represent features, where n0 is the output bit, k0 is the input bit, and m is the number of registers. K represents the constraint length, that is, the number m of shift registers in the convolutional code plus 1. In the cdma2000 proposal, the convolutional code encoder of the turbo code sub-(3, 1, 4) regression system is shown in Fig. 2 .

图2所示的是一个R=1/3码率的回归系统卷积码编码器(RSC)。20是移位寄存器,一共有三个移位寄存器,所以m=3,K=4。22是模2加法器。24是尾比特控制结构,当一帧数据输入完毕后,需要对20寄存器清零,这是将尾比特控制器开关切换到下方,通过三个节拍,将三个寄存器20内的比特作为输入依次清零。What Fig. 2 shows is a R = 1/3 code rate regression systematic convolution code encoder (RSC). 20 is a shift register, and there are three shift registers in total, so m=3, K=4. 22 is a modulo 2 adder. 24 is the tail bit control structure. After a frame of data input is completed, the 20 register needs to be cleared. This is to switch the tail bit controller switch to the bottom. Through three beats, the bits in the three registers 20 are used as input in sequence cleared.

宽带码分多址(WCDMA)提案中的特博码编码器的RSC子编码器与图2的相似,只是没有Y1的输出端,也就是(2,1,3)的回归系统卷积码编码器(RSC)。The RSC sub-encoder of the turbo code encoder in the wideband code division multiple access (WCDMA) proposal is similar to that in Figure 2, except that there is no output terminal of Y1, that is, the regression system convolutional code encoding of (2, 1, 3) controller (RSC).

特博码的解码器采用迭代递推的方法,通过多次迭代来不断提高解码精度。图3是特博码解码器的结构,其中的33,34是指软输入软输出(SISO)解码器。31是解打孔装置,对应于编码器中打孔器的逆操作,32是解交织器,对应于编码器中交织器的逆操作,还原交织前的顺序,35是符号判决器,当输入数据大于等于0时,输出1;当输入数据小于0时,输出0。The decoder of the Turbo code adopts an iterative and recursive method to continuously improve the decoding accuracy through multiple iterations. Fig. 3 is a structure of a turbo code decoder, in which 33 and 34 refer to a soft-input soft-output (SISO) decoder. 31 is a depuncturing device, corresponding to the inverse operation of the puncturer in the encoder, 32 is a deinterleaver, corresponding to the inverse operation of the interleaver in the encoder, and restores the order before interleaving, 35 is a symbol decision device, when the input When the data is greater than or equal to 0, output 1; when the input data is less than 0, output 0.

软输入软输出(SISO)解码器主要分为最大后验概率解码(MAP)和最大似然解码(SOVA)两种方法。Soft-Input Soft-Output (SISO) decoders are mainly divided into two methods: maximum a posteriori probability decoding (MAP) and maximum likelihood decoding (SOVA).

最大后验概率解码(MAP)方法是基于回归卷积编码网格图,通过前反向递归求得每一解码比特的最大后验概率对数比(LLR)。 The Maximum A Posteriori Probability Decoding (MAP) method is based on the regressive convolutional coding trellis graph, and obtains the Maximum Posteriori Probability Log Ratio (LLR) of each decoded bit through forward and backward recursion.

                 dk=1 if^(dk)>0  (2)d k =1 if^(d k )>0 (2)

                 dk=0 if^(dk)≤0其中dk对应编码输入数据;R1 N是解码器接收到的N组解码数据。引入定义: a k i ( m ) = P r ( d k = i , S k = m , R 1 k ) P r ( R 1 k ) = P r ( d k = i , S k = m / R 1 k ) - - - - ( 3 ) β k ( m ) = P r ( R k + 1 N / S k = m ) P r ( R k + 1 N / R 1 N ) - - - ( 4 ) d k =0 if^(d k )≤0 where d k corresponds to encoded input data; R 1 N is the N sets of decoded data received by the decoder. import definition: a k i ( m ) = P r ( d k = i , S k = m , R 1 k ) P r ( R 1 k ) = P r ( d k = i , S k = m / R 1 k ) - - - - ( 3 ) β k ( m ) = P r ( R k + 1 N / S k = m ) P r ( R k + 1 N / R 1 N ) - - - ( 4 )

       γi(Rk,m′,m)=Pr(dk=i,Rk,Sk=m/Sk-1=m′)          (5)γ i (R k , m', m)=P r (d k =i, R k , S k =m/S k-1 =m') (5)

其中:Sk=m表示第K时刻的第m个状态;α、β、γ分别对应状态转移网格图中的前项状态度量、反向状态度量和分支度量。此时解码后验概率对数逻辑似然比(LLR)可表示为:

Figure A0110135000061
α和β满足如下递推关系 a k i ( m ) = Σ m ′ Σ j = 0 1 γ i ( R k , m ′ , m ) a k - 1 j ( m ′ ) Σ m Σ m ′ Σ i = 0 1 Σ j = 0 1 γ i ( R k , m ′ , m ) a k - 1 j ( m ′ ) - - - ( 7 ) β k ( m ) = Σ m ′ Σ j = 0 1 γ i ( R k + 1 , m ′ , m ) β k + 1 ( m ′ ) Σ m Σ m ′ Σ i = 0 1 Σ j = 0 1 γ i ( R k + 1 , m ′ , m ) a k j ( m ′ ) - - - ( 8 ) 根据MAP解码理论,最大后验概率解码(MAP)方法由以下几步完成:步骤1:初始化状态度量α、β a 0 i ( 0 ) = 1 , a 0 i ( m ) = 0 - - ∀ m ≠ 0 - - i = 0,1 - - - - ( 9 ) Where: S k =m represents the mth state at the Kth moment; α, β, and γ correspond to the previous state metric, reverse state metric, and branch metric in the state transition grid diagram, respectively. At this time, the logarithmic likelihood ratio (LLR) of the decoding posterior probability can be expressed as:
Figure A0110135000061
α and β satisfy the following recurrence relation a k i ( m ) = Σ m ′ Σ j = 0 1 γ i ( R k , m ′ , m ) a k - 1 j ( m ′ ) Σ m Σ m ′ Σ i = 0 1 Σ j = 0 1 γ i ( R k , m ′ , m ) a k - 1 j ( m ′ ) - - - ( 7 ) β k ( m ) = Σ m ′ Σ j = 0 1 γ i ( R k + 1 , m ′ , m ) β k + 1 ( m ′ ) Σ m Σ m ′ Σ i = 0 1 Σ j = 0 1 γ i ( R k + 1 , m ′ , m ) a k j ( m ′ ) - - - ( 8 ) According to the MAP decoding theory, the maximum a posteriori probability decoding (MAP) method is completed by the following steps: Step 1: Initialize the state metrics α, β a 0 i ( 0 ) = 1 , a 0 i ( m ) = 0 - - ∀ m ≠ 0 - - i = 0,1 - - - - ( 9 )

         βN(0)=1,βN(m)=0  m≠0(10)β N (0)=1, β N (m)=0 m≠0(10)

步骤2:对接收到的每一码元,计算分支度量γ和前项状态度量αStep 2: For each symbol received, calculate the branch metric γ and the previous state metric α

步骤3:当接收序列的α计算完成后,计算后项状态度量β,并得到每一比特的后验概率对数逻辑似然比LLR,Step 3: After the calculation of α of the received sequence is completed, calculate the subsequent state measure β, and obtain the logarithmic likelihood ratio LLR of the posterior probability of each bit,

通常MAP解码方法主要按照上面三步进行,整个解码过程主要是分支度量γ、状态度量α、β和LLR的计算。Usually, the MAP decoding method is mainly carried out according to the above three steps, and the whole decoding process is mainly the calculation of branch metric γ, state metric α, β and LLR.

由于标准MAP解码中由很多指数和乘法运算,因此MAP解码的实现常在对数域进行,这就是Log MAP。这样标准MAP解码中的乘法运算就变成了加法,而加法运算可采用Jacobian公式简化,用查找表的方式实现。Jacobian公式如下:Since there are many exponential and multiplication operations in standard MAP decoding, the implementation of MAP decoding is often performed in the logarithmic domain, which is Log MAP. In this way, the multiplication operation in the standard MAP decoding becomes addition, and the addition operation can be simplified by the Jacobian formula and realized by a look-up table. The Jacobian formula is as follows:

ln(ex+ey)=max(x,y)+ln(1-exp{-|x-y|})max(x,y)+fc(|x-y|)    (11)l n (e x +e y )=max(x,y)+ln(1-exp{-|xy|})max(x,y)+f c (|xy|) (11)

从标准MAP转移到Log MAP,从算法上讲是等效的。在初始时刻状态度量的处理上,通常的方法是增大初始零状态度量的权值,从而使得从初始零状态发出的状态转移的状态度量值比其他状态度量值大。这样在Log MAP解码的实现过程中,只需要按照MAP解码的实现步骤,把整个解码器的实现转移到对数域即可。Moving from Standard MAP to Log MAP is algorithmically equivalent. In the processing of state metrics at the initial moment, the usual method is to increase the weight of the initial zero state metric, so that the state metric value of the state transition from the initial zero state is larger than other state metric values. In this way, in the implementation process of Log MAP decoding, it is only necessary to transfer the implementation of the entire decoder to the logarithmic domain according to the implementation steps of MAP decoding.

根据MAP解码方法步骤,在MAP解码的过程中,首先要进行前项状态度量α和后项状态度量的初始化,其中,状态度量按照(9)、(10)式进行初始化,由于在标准MAP解码中,在状态转移路径上,当前时刻的状态度量由上一时刻的状态度量和分支度量相乘得到,因此按照(9)、(10)式的初始化不会影响解码结果。当标准MAP解码转移到对数域时,标准MAP解码中状态度量计算中的相乘关系就变成了相加关系,这时无论进行怎样的初始化,在初始状态度量的计算过程中都会引入多余的状态转移路径,如图4中虚线所示。According to the steps of the MAP decoding method, in the process of MAP decoding, the initialization of the previous state metric α and the subsequent state metric must be performed first, wherein the state metrics are initialized according to formulas (9) and (10), because in the standard MAP decoding In , on the state transition path, the state metric at the current moment is multiplied by the state metric at the previous moment and the branch metric, so the initialization according to (9) and (10) will not affect the decoding result. When the standard MAP decoding is transferred to the logarithmic domain, the multiplication relationship in the state metric calculation in the standard MAP decoding becomes an additive relationship. At this time, no matter what kind of initialization is performed, redundant The state transition path of , as shown by the dotted line in Figure 4.

图4中,图4是前项状态度量初始状态转移网格图,图5是后项状态度量初始状态转移网格;43是状态转移路径;44代表某一时刻编码器的某一状态。In Fig. 4, Fig. 4 is the grid diagram of the initial state transition of the previous state metric, and Fig. 5 is the initial state transition grid of the subsequent state metric; 43 is the state transition path; 44 represents a certain state of the encoder at a certain moment.

在对初始状态度量处理时,通常的方法是零状态的状态度量取一个较大的值,而其他状态的状态度量值取一个相对较小的的值,这样尽量使得从零状态转移的状态度量值占较大的权重。但即使采用加大初始零状态状态度量的方式,也不能完全消除附加多余转移路径的状态度量对后级状态度量的计算,这样就相当于在状态度量的计算中引入了人为噪声,从而影响解码性能。When dealing with the initial state metrics, the usual method is to take a larger value for the state metrics of the zero state, and take a relatively small value for the state metrics of other states, so that the state metrics transferred from the zero state can be as far as possible The value takes a larger weight. However, even if the method of increasing the initial zero-state state metric is adopted, the calculation of the subsequent state metric by the state metric of the additional transition path cannot be completely eliminated, which is equivalent to introducing artificial noise in the calculation of the state metric, thereby affecting decoding. performance.

为此,本发明的目的是针对上述特博码解码方法因存在不能消除附加路径状态度量对后级状态度量的计算而引入了人为噪声以及影响解码性能的缺点,提出一种属于最大后验概率特博码解码方法及其解码器。For this reason, the purpose of the present invention is to address the shortcomings of the above-mentioned Turbo code decoding method that cannot eliminate the calculation of the additional path state metric to the subsequent stage state metric, which introduces artificial noise and affects the decoding performance, and proposes a method that belongs to the maximum a posteriori probability A turbo code decoding method and a decoder thereof.

为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solutions:

该特博码解码方法基于最大后验概率解码方法,首先对输入的解码软数据信息的前项状态度量α、反向状态度量β进行初始化,该方法还包括以下步骤:The turbo code decoding method is based on the maximum a posteriori probability decoding method, and first initializes the antecedent state metric α and the reverse state metric β of the input decoded soft data information, and the method also includes the following steps:

a,通过状态网格图对初始几个时刻状态度量进行计算控制,以消除多余状态度量的引入;a, Calculate and control the state metrics at the initial few moments through the state grid diagram to eliminate the introduction of redundant state metrics;

b,对接收到的每一输入解码软数据信息的码元计算出分支度量γ和前项状态度量α;b. Calculate the branch metric γ and the previous state metric α for each received symbol of the decoded soft data information;

c,当上述的码元前项状态度量α计算完成后,计算后项状态度量β,并得到每一比特的后验概率对数逻辑似然比LLR;c. After the calculation of the previous item state metric α of the above-mentioned symbol is completed, the subsequent item state metric β is calculated, and the posterior probability logarithmic logic likelihood ratio LLR of each bit is obtained;

d,经过上面步骤多次迭代后,输出软输出解码信息。d, After multiple iterations of the above steps, output the soft output decoding information.

一种特博码解码器,该解码器基于最大后验概率的状态转移网格图解码器,它包括用于计算状态转移路径中的分支度量值的分支度量γ计算模块、用于计算状态转移中的状态度量值的前项状态度量α计算模块及存储器、用于计算状态转移中的状态度量值的后项状态度量β计算模块及存储器、对前项状态度量α和后项状态度量β的初始化的状态度量初始化模块、查找表、逻辑似然比计算模块、产生软输入软输出解码的控制信号的总控制模块,分支度量γ计算模块接收解码输入数据,其输出接到前项状态度量α计算模块和后项状态度量β计算模块,前项状态度量α计算模块和后项状态度量β计算模块的计算结果分别存储到α、β存储器,查找表提供前项状态度量α计算模块、后项状态度量β计算模块、逻辑似然比计算模块查表,α、β存储器分别输出到逻辑似然比计算模块,该解码器还包括α初始状态度量计算控制模块和β初始状态度量计算控制模块,α初始状态度量计算控制模块与状态度量初始化模块双向接通,α初始状态度量计算控制模块和β初始状态度量计算控制模块分别对前项状态度量α计算模块和后项状态度量β计算模块进行状态转移路径控制,总控制模块提供整个解码器的控制信号,该信号分别接到γ计算模块、α初始状态度量计算控制模块和α计算模块、β初始状态度量计算控制模块和β计算模块。A Tebow code decoder, the decoder is based on the maximum a posteriori probability state transition trellis graph decoder, which includes a branch metric gamma calculation module for calculating the branch metric value in the state transition path, for calculating the state transition The former state measure α calculation module and memory of the state measure value in the state transition, the subsequent state measure β calculation module and memory used to calculate the state measure value in the state transition, the former state measure α and the latter state measure β Initialized state metric initialization module, lookup table, logic likelihood ratio calculation module, general control module that generates soft input and soft output decoding control signals, branch metric γ calculation module receives decoding input data, and its output is connected to the previous state metric α The calculation module and the subsequent state measure β calculation module, the calculation results of the previous state measure α calculation module and the subsequent state measure β calculation module are stored in the α and β memories respectively, and the lookup table provides the former state measure α calculation module, the latter The state metric β calculation module and the logic likelihood ratio calculation module look up the table, and the α and β memories are respectively output to the logic likelihood ratio calculation module. The decoder also includes an α initial state metric calculation control module and a β initial state metric calculation control module. The α initial state metric calculation control module and the state metric initialization module are bidirectionally connected, and the α initial state metric calculation control module and the β initial state metric calculation control module perform state measurement on the previous state metric α calculation module and the subsequent state metric β calculation module respectively. For transfer path control, the general control module provides control signals for the entire decoder, and the signals are respectively connected to the γ calculation module, the α initial state metric calculation control module and the α calculation module, the β initial state metric calculation control module, and the β calculation module.

由于本发明的方法在最大后验概率解码方法的基础上,增加了对通过状态网格图对初始几个时刻状态度量进行计算控制,消除了多余状态度量的引入;以及依本发明的方法而设计的解码器中,增设了α初始状态度量计算控制模块和β初始状态度量计算控制模块,对前项状态度量和初始后项状态度量进行状态转移路径控制以消除多余状态度量的引入,可以提高状态度量的精度和提高译码的性能。因此,经过初始状态度量计算控制处理的Log MAP解码性能比没有初始状态度量计算控制的解码性能提高大约有0.1db的性能增益。也就是说在保证相同误码率的情况下,有初始状态度量计算控制处理译码器可以容忍更差的输入信号质量(输入信号信噪比更小),在无线通信中,就可以减少发射机的功率,减少干扰,提高通信质量。Because the method of the present invention is on the basis of the maximum a posteriori probability decoding method, it increases the calculation and control of the state metrics at the initial few moments through the state grid diagram, eliminating the introduction of redundant state metrics; and according to the method of the present invention In the designed decoder, the α initial state metric calculation control module and the β initial state metric calculation control module are added, and the state transition path control is performed on the previous state metric and the initial subsequent state metric to eliminate the introduction of redundant state metrics, which can improve Accuracy of state metrics and improved decoding performance. Therefore, the decoding performance of the Log MAP after the initial state metric calculation control is improved by about 0.1db compared to the decoding performance without the initial state metric calculation control. That is to say, in the case of ensuring the same bit error rate, the initial state metric calculation control processing decoder can tolerate worse input signal quality (the input signal signal-to-noise ratio is smaller), and in wireless communication, the emission can be reduced Machine power, reduce interference, improve communication quality.

下面结合附图和实施例,对本发明的解码方法作一详细地说明:Below in conjunction with accompanying drawing and embodiment, the decoding method of the present invention is described in detail:

图1为传统的特博码编码器结构原理示意图。FIG. 1 is a schematic diagram of the structure and principle of a traditional turbo code encoder.

图2为图1中编码器的子编码器结构原理示意图。FIG. 2 is a schematic diagram of the structural principle of a sub-encoder of the encoder in FIG. 1 .

图3为对应于图1编码器的传统的解码器结构原理示意图。FIG. 3 is a schematic diagram of the structure and principle of a traditional decoder corresponding to the encoder in FIG. 1 .

图4为现有的解码方法中的前项状态度量初始状态转移网格图。FIG. 4 is a grid diagram of an initial state transition of antecedent state metrics in an existing decoding method.

图5为现有的解码方法中的后项状态度量初始状态转移网格图。Fig. 5 is a grid diagram of the initial state transition of the latter state metric in the existing decoding method.

图6为本发明的解码器结构原理示意图。Fig. 6 is a schematic diagram of the structural principle of the decoder of the present invention.

为了尽可能提高解码的性能,在Log MAP解码过程中必须消除多余状态转移路径的引入。在标准的MAP理论中,在计算前项状态度量α和后项状态度量β过程中,分支度量γ与前一时刻的状态度量α是相乘关系数,那么在状态度量计算时,只需要对初始时刻(前项状态度量α对应第一个时刻,后项状态度量β对应最后一个时刻按上述的公式(9)、(10)进行初始化,这样初始化为零的状态度量与相应状态转移度量γ的乘积仍然为零,因此应不会引入多余的状态度量对后面状态度量的计算产生影响,但是在Log MAP运算中,为了减少MAP解码的运算量,所有的运算都转移到对数域进行,这样原来的相乘就变成对数域的相加,此时如果还按照公式(9)、(10)进行状态度量的初始化,由于状态转移度量不为零,在计算下一个状态度量时就会引入多余的状态度量值,这个多余的状态度量值会影响后面的状态度量值计算准确性,因此对译码性能产生影响。In order to improve the decoding performance as much as possible, the introduction of redundant state transition paths must be eliminated during the Log MAP decoding process. In the standard MAP theory, in the process of calculating the previous state metric α and the subsequent state metric β, the branch metric γ and the previous state metric α are multiplicative coefficients, so when calculating the state metric, only need to The initial moment (the former state metric α corresponds to the first moment, and the latter state metric β corresponds to the last moment) is initialized according to the above formulas (9) and (10), so that the state metric initialized to zero and the corresponding state transition metric γ The product of is still zero, so it should not introduce redundant state metrics to affect the calculation of subsequent state metrics. However, in the Log MAP operation, in order to reduce the amount of MAP decoding operations, all operations are transferred to the logarithmic domain. In this way, the original multiplication becomes the addition of the logarithmic domain. At this time, if the state metric is initialized according to the formulas (9) and (10), since the state transition metric is not zero, when calculating the next state metric, it will be A redundant state metric value will be introduced, and this redundant state metric value will affect the calculation accuracy of the subsequent state metric value, thus affecting the decoding performance.

基于上述的考虑,本发明的方法是通过状态转移网格图对初始几个时刻状态度量计算进行控制,以消除多余状态度量引入,具体解码方法仍然是基于最大后验概率解码方法,首先对输入的解码软数据信息的前项状态度量α、反向状态度量β进行初始化,该方法还包括以下步骤:Based on the above considerations, the method of the present invention controls the calculation of the state metrics at the initial few moments through the state transition grid diagram to eliminate the introduction of redundant state metrics. The specific decoding method is still based on the maximum a posteriori probability decoding method. First, the input The previous item state measure α and the reverse state measure β of the decoded soft data information are initialized, and the method also includes the following steps:

a,通过状态网格图对初始几个时刻状态度量进行计算控制,以消除多余状态度量的引入;a, Calculate and control the state metrics at the initial few moments through the state grid diagram to eliminate the introduction of redundant state metrics;

b,对接收到的每一输入解码软数据信息的码元计算出分支度量γ和前项状态度量α;b. Calculate the branch metric γ and the previous state metric α for each received symbol of the decoded soft data information;

c,当上述的码元前项状态度量α计算完成后,计算后项状态度量β,并得到每一比特的后验概率对数逻辑似然比LLR;c. After the calculation of the previous item state metric α of the above-mentioned symbol is completed, the subsequent item state metric β is calculated, and the posterior probability logarithmic logic likelihood ratio LLR of each bit is obtained;

d,经过上面步骤多次迭代后,后输出软输出解码信息。d. After multiple iterations of the above steps, the soft output decoding information is finally output.

所述的步骤b中的分支度量γ和前项状态度量α均是转移概率,分支度量γ由信道模型和状态转移图决定,前项状态度量α由公式(7)计算得到。Both the branch metric γ and the previous state metric α in the step b are transition probabilities, the branch metric γ is determined by the channel model and the state transition diagram, and the previous state metric α is calculated by formula (7).

所述的步骤c中,当码元前项状态度量α计算完成后,就按公式(8)计算码元的后项状态度量,在计算后项状态度量β的同时就对每个码元对应的后验概率对数逻辑似然比LLR进行计算。In the described step c, after the calculation of the previous item state measure α of the symbol is completed, the subsequent item state measure of the symbol is calculated according to formula (8), and when calculating the latter item state measure β, it corresponds to each symbol The logarithmic likelihood ratio LLR of the posterior probability is calculated.

所述的步骤d中,所述的多次迭代是指每次迭代解码的对数似然比LLR输出经过解交织后作为下一次迭代解码的输入,经过反复循环,最后输出解码信息。In said step d, said multiple iterations means that the LLR output of each iterative decoding is deinterleaved and used as the input of the next iterative decoding, and after repeated cycles, the decoded information is finally output.

在本发明的上述方法中,所谓初始几个时刻的状态度量含义是,按照状态转移网格图(结合图4、5所示),在前项状态度量和后项状态度量的初始时刻并不是所有状态进行状态转移,因此需要根据实际状态转移,对初始状态度量的计算进行控制,正如在背景技术中所描述的那样,在对初始时刻状态度量处理时,通常的方法是零状态度量取一个较大的值,而其它状态度量值取一个相对较小的值以尽量使从零状态转移的状态度量值占较大的权重。In the above-mentioned method of the present invention, the so-called state measurement meaning of several initial moment is, according to the state transition grid diagram (shown in conjunction with Fig. 4, 5), the initial moment at the former state measurement and the subsequent state measurement is not All states undergo state transition, so it is necessary to control the calculation of the initial state metric according to the actual state transition. As described in the background technology, when processing the state metric at the initial moment, the usual method is to take a zero state metric and take one A larger value, while the other state metrics take a relatively small value to try to make the state metrics transitioning from the zero state account for a larger weight.

由于在状态度量计算的初始时刻,并不是所有状态进行状态转移,因此,如果不对初始计算进行控制,就会对所有状态进行状态度量的计算,这样会引入多余状态度量,这些多余状态度量会对以后的状态度量计算产生影响,从而降低译码性能。对初始时刻状态度量进行计算控制,就会消除多余状态度量产生的影响,从而能提高译码的性能。Since not all states undergo state transition at the initial moment of state metric calculation, if the initial calculation is not controlled, state metric calculation will be performed for all states, which will introduce redundant state metrics, which will affect Subsequent state metric calculations have an impact, thereby degrading decoding performance. Calculating and controlling the state metrics at the initial moment will eliminate the influence of redundant state metrics, thereby improving the decoding performance.

在传统特博码解码器的软输入输出(SISO)解码器的基础上,根据本发明的上述方法,该解码器基于最大后验概率的状态转移网格图解码器,它仍然包括用于计算状态转移路径中的分支度量值的分支度量γ计算模块、用于计算状态转移中的状态度量值的前项状态度量α计算模块及存储器、用于计算状态转移中的状态度量值的后项状态度量β计算模块及存储器、对前项状态度量α和后项状态度量β的初始化的状态度量初始化模块、查找表、逻辑似然比计算模块、产生软输入软输出解码的控制信号的总控制模块,分支度量γ计算模块接收解码输入数据,其输出接到前项状态度量α计算模块和后项状态度量β计算模块,前项状态度量α计算模块和后项状态度量β计算模块的计算结果分别存储到α、β存储器,查找表提供前项状态度量α计算模块、后项状态度量β计算模块、逻辑似然比计算模块查表,α、β存储器分别输出到逻辑似然比计算模块。该解码器还包括α初始状态度量计算控制模块和β初始状态度量计算控制模块,α初始状态度量计算控制模块与状态度量初始化模块双向接通,α初始状态度量计算控制模块和β初始状态度量计算控制模块分别对前项状态度量α计算模块和后项状态度量β计算模块进行状态转移路径控制,总控制模块提供整个解码器的控制信号,该信号分别接到分支度量γ计算模块、α初始状态度量计算控制模块和前项状态度量α计算模块、β初始状态度量计算控制模块和后项状态度量β计算模块。On the basis of the soft-input-output (SISO) decoder of the traditional Turbo code decoder, according to the above-mentioned method of the present invention, the decoder is based on the state transition trellis graph decoder of the maximum a posteriori probability, and it still includes for calculating The branch metric gamma calculation module of the branch metric value in the state transition path, the previous state metric α calculation module and memory for calculating the state metric value in the state transition, and the subsequent item state for calculating the state metric value in the state transition Metric β calculation module and memory, a state metric initialization module for initializing the previous state metric α and the subsequent state metric β, a lookup table, a logic likelihood ratio calculation module, and a total control module for generating control signals for soft input and soft output decoding , the branch metric γ calculation module receives the decoded input data, and its output is connected to the previous state metric α calculation module and the subsequent state metric β calculation module, and the calculation results of the previous state metric α calculation module and the subsequent state metric β calculation module are respectively Stored in the α and β memories, the lookup table provides the previous state measure α calculation module, the subsequent state measure β calculation module, and the logical likelihood ratio calculation module look-up table, and the α and β memories are respectively output to the logical likelihood ratio calculation module. The decoder also includes an α initial state metric calculation control module and a β initial state metric calculation control module, the α initial state metric calculation control module is bidirectionally connected with the state metric initialization module, the α initial state metric calculation control module and the β initial state metric calculation The control module performs state transition path control on the previous state metric α calculation module and the subsequent state metric β calculation module, and the general control module provides the control signal of the entire decoder, which is respectively connected to the branch metric γ calculation module, α initial state The metric calculation control module and the previous state metric α calculation module, the β initial state metric calculation control module and the subsequent state metric β calculation module.

请结合图6所示,在该图中,单元51是分支度量γ(Gamma)计算模块,用于计算状态转移路径中的分支度量值;单元52是前项状态度量α(Alpha)计算模块,用于计算状态转移中的前项状态度量值;单元53是后项状态度量β(Beta)计算模块,用于计算状态转移中的后项状态度量值;单元54是逻辑似然比LLR计算模块,根据状态度量值计算后验概率对数比;单元55是状态度量初始化模块,完成对前项状态度量α和后项状态度量β的初始化;单元56是一个八值查找表;单元57是α存储器,用于存储计算得到的前项状态度量值;单元58是β存储器,用于存储计算得到的后项状态度量值;单元59是总控制模块,产生软输入软输出解码的控制信号;单元5a是初始前项状态度量α计算控制模块,使得前项状态度量的计算完全按照实际的状态转移路径进行(如图4、4b中的实线所示);单元5b是初始后项状态度量β计算控制模块,使得后项状态度量的计算完全按照实际的状态转移路径进行。Please in conjunction with shown in Fig. 6, in this figure, unit 51 is a branch metric γ (Gamma) calculation module, is used for calculating the branch metric value in the state transfer path; Unit 52 is a preceding item state metric α (Alpha) calculation module, For calculating the previous item state metric value in the state transition; Unit 53 is the latter item state metric β (Beta) calculation module, for calculating the latter item state metric value in the state transition; Unit 54 is the logical likelihood ratio LLR calculation module , calculate the posterior probability logarithmic ratio according to the state metric value; unit 55 is a state metric initialization module, and completes the initialization of the former item state metric α and the subsequent item state metric β; unit 56 is an eight-valued lookup table; unit 57 is α The memory is used to store the previous item state metric value that is calculated; the unit 58 is a β memory, and is used to store the subsequent item state metric value that is calculated; the unit 59 is a total control module that generates soft input and soft output decoding control signals; the unit 5a is the calculation control module of the initial antecedent state metric α, so that the calculation of the antecedent state metric is completely carried out according to the actual state transition path (as shown by the solid lines in Figure 4 and 4b); unit 5b is the initial subsequent state metric β The calculation control module makes the calculation of the latter state measure completely follow the actual state transition path.

本发明的方法和解码器经过了解码定点仿真,经过初始控制处理的Log MAP解码性能比没有初始控制的解码性能提高大约有0.1dB的性能增益。The method and the decoder of the present invention have undergone decoding fixed-point simulation, and the Log MAP decoding performance after initial control processing is improved by about 0.1dB performance gain compared with the decoding performance without initial control.

仿真条件:Simulation conditions:

解码码块长度5114Decoded code block length 5114

解码输入数据6Bits量化Decoding input data 6Bits quantization

没有初始状态度量计算控制模块时,初始化如下: ln [ a 0 i ( 0 ) ] = 63 , ln [ a 0 i ( m ) ] = 0 - - ∀ m ≠ 0 - - i = 0,1 - - - ( 12 ) When there is no initial state metric calculation control module, the initialization is as follows: ln [ a 0 i ( 0 ) ] = 63 , ln [ a 0 i ( m ) ] = 0 - - ∀ m ≠ 0 - - i = 0,1 - - - ( 12 )

               ln[βN(0)]=63,ln[βN(m)]=0  m≠0  (13)ln[β N (0)]=63, ln[β N (m)]=0 m≠0 (13)

有初始状态度量计算控制模块时,仿真性能不受状态度量初始值的影响。在保证相同误码率的情况下,有初始状态度量计算控制处理译码器可以容忍更差的输入信号质量(输入信号信噪比更小),在无线通信中,就可以减少发射机的功率,减少干扰,提高通信质量。在CDMA2000和WCDMA系统中,则可以减少基站的发射功率,提高小区的容量。When there is an initial state metric calculation control module, the simulation performance is not affected by the initial value of the state metric. In the case of ensuring the same bit error rate, the initial state metric calculation control processing decoder can tolerate worse input signal quality (the input signal signal-to-noise ratio is smaller), and in wireless communication, the power of the transmitter can be reduced , Reduce interference and improve communication quality. In CDMA2000 and WCDMA systems, the transmit power of the base station can be reduced and the capacity of the cell can be increased.

Claims (5)

1.一种特博码解码方法,该解码方法基于最大后验概率的状态转移网格图解码的方法,首先对输入的解码软数据信息码元的前项状态度量α、反向状态度量β进行初始化,其特征在于,该方法还包括以下步骤:1. A turbo code decoding method, the decoding method is based on the method of state transition lattice diagram decoding of maximum a posteriori probability, at first to the preceding item state measure α, reverse state measure β of the decoding soft data information symbol of input Carry out initialization, it is characterized in that, this method also comprises the following steps: a,通过状态网格图对初始几个时刻状态度量进行计算控制,以消除多余状态度量的引入;a, Calculate and control the state metrics at the initial few moments through the state grid diagram to eliminate the introduction of redundant state metrics; b,对接收到的每一输入解码软数据信息的码元计算出分支度量γ和前项状态度量α;b. Calculate the branch metric γ and the previous state metric α for each received symbol of the decoded soft data information; c,当上述的码元前项状态度量α计算完成后,计算后项状态度量β,并得到每一比特的后验概率对数逻辑似然比LLR;c. After the calculation of the previous item state metric α of the above-mentioned symbol is completed, the subsequent item state metric β is calculated, and the posterior probability logarithmic logic likelihood ratio LLR of each bit is obtained; d,经过上面步骤多次迭代后,输出软输出解码信息。d, After multiple iterations of the above steps, output the soft output decoding information. 2、如权利要求1所述的特博码解码方法,其特征在于:所述的步骤b中的分支度量γ和前项状态度量α均是转移概率,分支度量γ由信道模型和状态转移图决定。2. The turbo code decoding method according to claim 1, characterized in that: both the branch metric γ and the previous state metric α in the step b are transition probabilities, and the branch metric γ is determined by the channel model and the state transition diagram Decide. 3、如权利要求1所述的特博码解码方法,其特征在于:所述的步骤c中,在计算后项状态度量β的同时就对每个码元对应的后验概率对数逻辑似然比LLR进行计算。3. The turbo code decoding method as claimed in claim 1, characterized in that: in the step c, when calculating the posterior state measure β, the logarithmic logarithm of the posterior probability corresponding to each symbol is Then the ratio LLR is calculated. 4、如权利要求1所述的特博码解码方法,其特征在于:所述的步骤d中,所述的多次迭代是指每次迭代解码的对数似然比LLR输出经过解交织后作为下一次迭代解码的输入,经过反复循环,最后输出解码信息。4. The turbo code decoding method according to claim 1, characterized in that: in the step d, the multiple iterations refer to the log likelihood ratio LLR output of each iteration decoding after deinterleaving As the input of the next iterative decoding, after repeated loops, the decoding information is finally output. 5、一种特博码解码器,该解码器基于最大后验概率的状态转移网格图解码器,它包括用于计算状态转移路径中的分支度量值的分支度量γ计算模块、用于计算状态转移中的状态度量值的前项状态度量α计算模块及存储器、用于计算状态转移中的状态度量值的后项状态度量β计算模块及存储器、对前项状态度量α和后项状态度量β的初始化的状态度量初始化模块、查找表、逻辑似然比计算模块、产生软输入软输出解码的控制信号的总控制模块,分支度量γ计算模块接收解码输入数据,其输出接到前项状态度量α计算模块和后项状态度量β计算模块,前项状态度量α计算模块和后项状态度量β计算模块的计算结果分别存储到α、β存储器,查找表提供前项状态度量α计算模块、后项状态度量β计算模块、逻辑似然比计算模块查表,α、β存储器分别输出到逻辑似然比计算模块,其特征在于:该解码器还包括α初始状态度量计算控制模块和β初始状态度量计算控制模块,α初始状态度量计算控制模块与状态度量初始化模块双向接通,α初始状态度量计算控制模块和β初始状态度量计算控制模块分别对前项状态度量α计算模块和后项状态度量β计算模块进行状态转移路径控制,总控制模块提供整个解码器的控制信号,该信号分别接到γ计算模块、α初始状态度量计算控制模块和α计算模块、β初始状态度量计算控制模块和β计算模块。5. A turbo code decoder, the decoder is based on the state transition lattice graph decoder of maximum a posteriori probability, which includes a branch metric gamma calculation module for calculating the branch metric value in the state transition path, for calculating The previous state measure α calculation module and memory of the state measure value in the state transition, the subsequent state measure β calculation module and memory for calculating the state measure value in the state transition, and the former state measure α and the latter state measure β initialization state metric initialization module, lookup table, logic likelihood ratio calculation module, general control module that generates soft input and soft output decoding control signals, branch metric γ calculation module receives decoding input data, and its output is connected to the previous state The metric α calculation module and the subsequent state metric β calculation module, the calculation results of the previous state metric α calculation module and the subsequent state metric β calculation module are stored in α and β memories respectively, and the lookup table provides the previous state metric α calculation module, The subsequent state measure β calculation module and the logical likelihood ratio calculation module look up the table, and the α and β memories are respectively output to the logical likelihood ratio calculation module. It is characterized in that: the decoder also includes an α initial state metric calculation control module and a β initial The state metric calculation control module, the α initial state metric calculation control module and the state metric initialization module are bidirectionally connected, and the α initial state metric calculation control module and the β initial state metric calculation control module respectively calculate the previous state metric α calculation module and the subsequent state The metric β calculation module performs state transfer path control, and the general control module provides control signals for the entire decoder, which are respectively connected to the γ calculation module, the α initial state metric calculation control module, and the α calculation module, β initial state metric calculation control module and β calculation module.
CNB011013508A 2001-01-04 2001-01-04 Turbo code decoding method and decoder Expired - Lifetime CN1145266C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB011013508A CN1145266C (en) 2001-01-04 2001-01-04 Turbo code decoding method and decoder

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB011013508A CN1145266C (en) 2001-01-04 2001-01-04 Turbo code decoding method and decoder

Publications (2)

Publication Number Publication Date
CN1362789A true CN1362789A (en) 2002-08-07
CN1145266C CN1145266C (en) 2004-04-07

Family

ID=4652032

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB011013508A Expired - Lifetime CN1145266C (en) 2001-01-04 2001-01-04 Turbo code decoding method and decoder

Country Status (1)

Country Link
CN (1) CN1145266C (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101079642A (en) * 2002-08-20 2007-11-28 恩益禧电子股份有限公司 Method for decoding error correcting code, its program and its device
CN102611461A (en) * 2012-03-28 2012-07-25 哈尔滨工业大学 Decoder of AR4JA code with code rate of 1/2 and information bit length of 1024 bits
CN104468023A (en) * 2013-09-22 2015-03-25 中兴通讯股份有限公司 Decoding method and decoding device
CN105790775A (en) * 2016-05-19 2016-07-20 电子科技大学 Probability calculation unit based on probability Turbo encoder

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101079642A (en) * 2002-08-20 2007-11-28 恩益禧电子股份有限公司 Method for decoding error correcting code, its program and its device
CN102611461A (en) * 2012-03-28 2012-07-25 哈尔滨工业大学 Decoder of AR4JA code with code rate of 1/2 and information bit length of 1024 bits
CN104468023A (en) * 2013-09-22 2015-03-25 中兴通讯股份有限公司 Decoding method and decoding device
CN104468023B (en) * 2013-09-22 2019-11-05 南京中兴软件有限责任公司 Coding/decoding method and device
CN105790775A (en) * 2016-05-19 2016-07-20 电子科技大学 Probability calculation unit based on probability Turbo encoder
CN105790775B (en) * 2016-05-19 2019-01-29 电子科技大学 A kind of probability calculation unit based on probability Turbo decoder

Also Published As

Publication number Publication date
CN1145266C (en) 2004-04-07

Similar Documents

Publication Publication Date Title
CN1154236C (en) Error correction coding type digital transmission method
US6813742B2 (en) High speed turbo codes decoder for 3G using pipelined SISO log-map decoders architecture
CN1258885C (en) Turbo Decoding Method with Error Message Recoding and Feedback
CN1306713C (en) Error correction in CDMA mobile communication system using turbo codes
CN1189935A (en) Parallel-linked truncated convolutional codes and their decoders
CN1761160A (en) Decoding method and device
CN114826284A (en) Iterative decoding method based on extended Turbo code and continuous phase modulation
US7886209B2 (en) Decoding device, decoding method, and receiving apparatus
US7027521B2 (en) Digital transmission method of the error correcting coding type
CN1157854C (en) A High Speed Turbo Code Decoder
CN1183687C (en) TURBO code encoder and encoding method
CN1157883C (en) Maximal posterior probability algorithm of parallel slide windows and its high-speed decoder of Turbo code
CN109831217B (en) A turbo code decoder, a component decoder for a turbo code and a component decoding method
CN1211931C (en) Memory architecture for MAP decoder
CN1145266C (en) Turbo code decoding method and decoder
CN1148006C (en) Method and decoder for decoding turbo code
CN1147169C (en) Decoding method and decoder for Turbo code
CN1129257C (en) Maximum-likelihood decode method f serial backtracking and decoder using said method
CN1773867A (en) Method for decoding Turbo code
CN115529048A (en) A Turbo Code Decoding Method Based on Linear Approximation and Sliding Window
CN1204693C (en) Stochastic system Turbo code coding and decoding method
CN1398047A (en) Turbine coder-decoder with parallel slide windows and its implementation method
CN2506034Y (en) Turbo decoder
CN1710815A (en) A High-Speed Maximum A Posteriori Turbo Decoding Method in Non-logarithmic Domain
CN1883120A (en) Decoder apparatus and decoding method

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
CX01 Expiry of patent term
CX01 Expiry of patent term

Granted publication date: 20040407