Interleaving method, interleaving apparatus, turbo encoding method, and turbo encoder

Abstract: An apparatus for interleaving a data sequence comprises means for receiving a data sequence (40) having a plurality of blocks (B 1 , B 2 ...B N ) where each block has a length based on a prime number (P). For example, 664 bits of input data are divided into 8 blocks of 83 bits each. The blocks are then written into a two-dimensional array in a first stage (41) of the processing. At a second stage (42) an intra-permutation is carried out. The intra-permutation of the second stage is performed using an address table derived from a Galois field of characteristic P. Finally, an inter-permutation is performed in which the order of the rows of intra-permutated data is altered. Data can then be read out of the two-dimensional buffer. The invention also relates to an interleaving method and a turbo encoding method and apparatus based on the interleaving method and apparatus respectively.

Description: BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention generally relates to the turbo encoding technique that effectively copes with a burst error. More particularly, the present invention is concerned with an interleaving method, an interleaving apparatus, a turbo encoding method, and a turbo encoder in which pruning is not performed at all or only a small number of bits is pruned away, so that computational complexity can be reduced. The present invention can be applied to fields required to improve reliability of communications using error correction codes, such as digital transmissions and digital recording. The present invention is particularly effective in fields which require flexibility of communications such as multimedia. 2. Description of the Related Art There are a number of documents relating to turbo encoding: Hughes Network Systems; Tdoc SMG2 UMTS-L1 765/98 discloses a Galois-Field Interleaver for turbo codes in which input data is read row by row into a 2 c columns x 2 r rows interleaving matrix, row and column permutations are performed to randomize the data position and the data is read out column by column. The permutations of the data of each row is based on powers of a primitive element in the Galois field GF(2 c ). Canon CRF; Tdoc SMG2 UMTS-L1 571/98 discloses a one-dimensional Algebraic-Interleaver for turbo-codes; Nortel Neworks; Tdoc SMG2 UMTS-L1 51/99 discloses a 2-dimensional Algebraic Interleaver for turbo codes; NTT DoCoMo; Tdoc SMG2 UMTS-L1 027/99 discloses a Multi-stage Interleaver (MIL) for turbo codes, the interleaver using a MIL lookup table. Recently, a turbo encoder has been proposed which utilizes a code having high capability in error correction. Such a turbo encoder is made up of a plurality of encoders, which are coupled together via an interleaver (which is means for performing interleaving processing) in order to reduce the correlation between redundant sequences associated with the respective encoders. The interleaver is a key element which determines the performance of turbo encoding. Figs. 1A and 1B show an example of the turbo encoder. As shown in Fig. 1A, the turbo encoder includes recursive systematic convolutional encoders 12 (RSC1) and 13 (RSC2), and an interleaver 11. As shown in Fig. 1B, each of the encoders 12 and 13 is made up of adders 14 and 15 and delay elements (D) 16 and 17 connected as shown therein. The turbo encoder receives an input data sequence d (K bits) and outputs encoded data sequences X1 - X3. In order to reduce the correlation between the redundant bits X1 and X2, the interleaver 11 is provided at the input side of the encoder (RSC2) 13. As shown in Fig. 1C, a turbo decoder is made up to two decoders 1 and 2, two interleavers 3 and 4, and a deinterleaver 5. In digital systems, permutation in interleaving is performed on a given unit basis of bit or symbol. The permutation is implemented by using a buffer or a pattern for permutation. When the buffer is used, data is written therein and is then read therefrom in a different sequence. When the pattern for permutation is used, data is permuted by referring to the pattern, which describes information concerning a permutation based on interleaving. The pattern described will be referred to as an interleave pattern. A description will now be given of bit-based permutation processing using the interleave pattern. Fig. 2 shows interleaving of a 16-bit sequence. In Fig. 2, a 16-bit sequence 67 is interleaved on the bit unit basis by referring to an interleave pattern table 68, which defines a sequence of interleaving. More particularly, the zeroth to 15th bits of the sequence 67 is written into a two-dimensional buffer by the interleave pattern table 68, and is then read therefrom in the order of 0, 8, 4, 12, ..., as indicated by an arrow in Fig. 2. Thus, a bit sequence after interleaving is obtained as shown in Fig. 2. The interleaver involved in interleaving is required to have the following three capabilities of: (1) handling a variety of frame lengths (for example, thousands to ten thousands); (2) producing the interleaved sequence with a small number of parameters; and (3) reducing computational complexity in creating the interleave pattern. As to the first capability, if the parameters are merely prepared for all of the different frame lengths, a huge number of parameters will be prepared, and a huge memory capacity will be required to store the parameters. Thus, the above is impractical. There is another disadvantage that it takes a long time to compute respective optimal parameters for each of the different frame lengths. The above problems may be solved by designing interleavers with a small number of parameters. This is related to the second capability. Interleavers equal in number to a power of 2 are prepared and pruning of data is performed. However, pruning of data requires an increased number of parameters for optimization, and the interleavers may not have good performance with respect to all of the frame lengths. That is, the interleavers have good performance for some frame lengths but do not operate well for other frame lengths. Reduction in the amount of data to be pruned away is also related to the third capability. In this regard, the present inventors have proposed an improvement in priming and performance (International Application No. PCT/JP98/05027). However, even the proposed improvement has high computational complexity in producing the interleave patterns. SUMMARY OF THE INVENTION It is a general object of the present invention to provide an improved interleaving technique having the above-mentioned three capabilities. A more specific object of the present invention is to provide an interleaving method, an interleaving apparatus, a turbo encoding method and a turbo encoder capable of efficiently achieving randomization of input sequences of a variety of frame lengths with reduced computational complexity. The invention is an interleaving method, an interleaver, a turbo encoding method, and a turbo encoder, as defined respectively in Claims 1, 10, 17 and 20. BRIEF DESCRIPTION OF THE DRAWINGS Other objects, features and advantages of the present invention will more apparent from the following detailed description of preferred embodiments when read in conjunction with the accompanying drawings in which: Figs. 1A and 1B are block diagrams of a conventional turbo encoder; Fig. 1C is a block diagram of a turbo decoder; Fig. 2 is a diagram of a conventional interleaving method of interleaving a 16-bit input data sequence; Fig. 3 is a block diagram of a turbo encoder according to a first embodiment of the present invention; Fig. 4 shows prime numbers up to 200; Fig. 5 is a diagram illustrating a first possible configuration of an interleaver shown in Fig. 3; Fig. 6 shows prime numbers less than 150 and associated primitive roots; Figs. 7A and 7B respectively show examples of sequence permutation tables; Fig. 8 is a diagram illustrating a second possible configuration of the interleaver shown in Fig. 3; Fig. 9 is a diagram illustrating a third possible configuration of the interleaver shown in Fig. 3; Fig. 10 is a flowchart of the turbo encoder according to the present invention; . Fig. 11 is a graph for explaining an error floor in turbo codes; Fig. 12 is a block diagram of a turbo encoder according to a second embodiment of the present invention; Fig. 13 is a diagram illustrating a fourth possible configuration of the interleaver shown in Fig. 3; and Figs. 14A, 14B and 14C respectively show examples of sequence permutation tables used in the fourth embodiment of the present invention. DESCRIPTION OF THE PREFERRED EMBODIMENTS A description will be given, with reference to the accompanying drawings, of embodiments of the present invention. First Embodiment Fig. 3 is a block diagram of a turbo encoder according to the first embodiment of the present invention. In Fig. 3, parts that are the same as those shown in the previously described figures are given the same reference numbers. The turbo encoder shown in Fig. 3 differs from that shown in Fig. 1A in the following: (a) a bit addition process part 21 is newly added; (b) an interleaver of a new configuration is used; and (c) a puncture processing part 23 is newly added. A detailed description will be given of the turbo encoder shown in Fig. 3 by further referring to a flowchart of Fig. 10. Bit Addition Process As preprocessing of interleaving, the bit addition process part 21 adjusts the input sequence so that it has a suitable number of bits for interleaving (steps 101 - 103 shown in Fig. 10). The bit addition process can be implemented by various types of conventional error correction coding. In this case, the bit addition process part 21 is a CRC encoder. Of the various types of conventional error correction coding, it is preferable to use a type of bit repetition in which bits are periodically repeated because it is flexible and is easily implemented. A detailed description will be given of the bit repetition in a case where the number of bits input to the turbo encoder is N IN (which corresponds to K in Fig. 3). The bit repetition includes the following four steps (1) - (4). At step (1), N IN , is divided by 8, and the resultant value n is obtained. The reason why N IN , is divided by 8 will be described later. At step (2), the prime number P that is greater than n and closest to n is obtained. At step (3), the difference between 8 times of P and N IN is calculated, and the resultant value is denoted as a. At step (4), a bits (dummy bits) are added to N IN bits of the input sequence. An example for N IN = 650 will be described below. At step (1), n = 650/8 = 81.25. At step (2), the prime number that is greater than 81.25 (= n) and closest thereto is 83 from a table shown in Fig. 4. That is, P = 83. At step (3), 83*8 = 664, and thus a = 14. That is, the number of dummy bits to be added to N IN is 14. At step (4), 14 dummy bits are added to the input sequence of 650 (=N IN ) bits. For example, the 14 dummy bits are added to the end of the 650-bit input sequence. The number (N IN + a) of bits thus obtained, that is, the number (K + a) of bits in Fig. 3 is divided by 8 without exception, and the quotient is always the prime number. The reason why 8 is used is that the number of rows of a two-dimensional array handled at a first stage of interleaving performed in the interleaver 22 is equal to 8, as will be described in detail later. It is possible to employ an arbitrary number such as 10 or 20 when the two-dimensional array used in the interleaver 22 has 10 or 20 rows. In such a case, 10 or 20 is used at the aforementioned steps (1) - (3). It can be seen from the above description that the process performed by the bit addition process part 21 obtains the number of rows of the two-dimensional array at step 101 shown in Fig. 10, and determines the number of columns thereof which is the prime number, and adds, to the input sequence, the number of dummy bits equal to the difference between the product of the numbers of rows and columns and the number of bits of the input sequence. Besides the bit repetition, it is possible to use block codes or convolutional codes in order to implement the bit addition process. It is also possible to employ a simple method of adding known bits to a known position. Interleaver 22 The interleaver 22 has, for example, one of three different configurations described-below. (First configuration) Fig. 5 shows the first configuration of the interleaver 22. The interleaver 22 having the first configuration includes the first, second and third stages 41, 42 and 43. The first stage 41 corresponds to step 104 shown in Fig. 10. The second stage 42 corresponds to steps 105 - 108 shown in Fig. 10. The third stage 43 corresponds to step 109 shown in Fig. 10. (1) First stage 41: An input sequence 40 (that is the output of the bit addition process part 21 and consists of, for example 664 bits) is divided into N. In Fig. 5, the 664-bit input sequence 40 is divided into eight blocks B 1 - B 8 , which are then written into a two-dimensional array (buffer), which consists of 8 rows and 83 columns. It will be noted that the 664-bit input sequence 40 includes 650 information bits and 14 dummy bits. The input sequence 40 of 664 bits can be divided by 8 and the quotient thereof is the prime number 83. Thus, the number of rows of the two-dimensional buffer is 8, and the number of columns thereof is the prime number 83. (2) Second stage 42: An intra-permutation is performed, as will be described later. In the intra-permutation, the sequence of the bits arranged in each row is permuted. (3) Third stage 43: An inter-permutation is performed in which the order of the rows arranged in the two-dimensional buffer is permuted. The inter-permutation uses, for example, a inter permutation pattern obtained by learning (directed to lengthening the free distance). One row is the unit of the inter-permutation. After the first through third stages are performed, data is read from the two-dimensional buffer in the longitudinal direction (column direction) at step 110 shown in Fig. 10, whereby an interleaved coded sequence 44 can be obtained. A further description will now be given of the second stage 42. The intra-permutation at the second stage 42 uses a table created by the following steps S1 - S7 as an address table, and processes input data written into the two-dimensional buffer by referring to the address table. Step S1: Step S1 is to obtain the primitive root g 0 of the Galois field of the characteristic P (which corresponds to the number of columns and is equal to 83 in the case shown in Fig. 5) at step 105 shown in Fig. 10, and to create a table to described in the order of exponential (power) notation of the primitive root in which the elements of the Galois field are expressed by anti-logarithms, which are arranged in the order of exponential notation. The primitive root g 0 of the Galois field of the characteristic P can be chosen from the table shown in Fig. 6. In other words, step S1 is to compute a mapping sequence c(i) for intra-row permutation defined below: c(i) = (g 0 i )(modP) where i = 0, 1, 2, .... (P-2), and c(P-1) = 0. In the case of P = 83, the primitive root g 0 (which corresponds to q in Fig. 6) of 83 is 2, and the elements of the associated Galois field are 0, 1, 2, ..., 82. Thus, c(0), c(1), c(2), ..., c(82) are as follows: c(0), c(1), c(2), ..., c(82) = 2 0 (mod83), 2 1 (mod83), 2 1 (mod83), ..., 2 82 (mod83) = 1, 2, 4, 8, 16, 32, 64. 45. 7, 14, ..., 42, 0 The table t0 can be formed from the above one-dimensional sequence, as shown in Fig. 7B, in which the combination of the numerals arranged in the transverse-axis and longitudinal-axis directions of the table t0 describes the exponent. For example, the combination of numerals of 1 and 6 respectively in the transverse- and longitudinal-axis directions indicates an exponent of 16. It can be seen from the table t0 that the results of, for example, 2 1 (mod83) and 2 16 (mod83) are respectively 4 and 49. The result of 2 82 (mod83) is set equal to 0. Step S2: The table t0 is defined as a sequence permutation table for the first row of the two-dimensional buffer. Step S2 corresponds to a case when a parameter I indicative of the row number is set equal to 1 at step 106 shown in Fig. 10. That is, the numbers defined in the sequence permutation table t0 indicate the bit allocations after permutation. As shown in Fig. 7B, the sequence permutation table t0 includes the one-dimensional sequence (pattern) starting from the left uppermost position thereof: Table t0: 1, 2, 4, 8, 16, ..., 42, 0 For example, let us assume that the input data mapped into the first row of the two-dimensional buffer is as follows: A0, A1, A2, A3, ..., A82 The bits in the first row of the two-dimensional buffer are permuted by referring to the sequence permutation table t0. For example, bits A0 and A1 which correspond to "1" and "2" in the sequence permutation table t0 are placed in the original positions even after permutation. Bit A2 which corresponds to "4" in the table t0 is permuted so as to be located in the fourth position after permutation. Similarly, bit A3 which corresponds to "8" in the table t0 is permuted so as to be located in the eighth position after permutation. Bit A82 which corresponds to the last number "0" in the table t0 is placed in the original position. The above permutation is performed at step 107 shown in Fig. 10. The data obtained after permutation of sequence (2) has the following sequence: A0, A1, A72, A2, A27, A76. A8, ..., A82 Step S3: Step S3 is to obtain (N - 1) numbers (integers) which have the mutually prime relationship with respect to (P - 1) where N is an integer greater than 2 and denotes the number of rows. In the example being considered, P = 83 and N = 8. Thus, 7 (= 8 - 1) integers which are numbers p1, p2, p3, p4, p5, p6 and p7 that have the mutually prime relationship with respect to 82 (= P - 1 = 83 - 1) are obtained. The 7 integers p1, p2, p3, p4, p5, p6 and p7 are respectively 3, 5, 7, 11, 13, 17 and 19 (1 and 2 are excluded). Step S4: At step 108 of Fig. 10, it is determined whether I is smaller than N (the number of rows). If the answer of step 108 is YES, the process returns to step 107. In the example being considered, I is incremented by 1 and is then equal to 2. Then, a sequence permutation table for permutation of data in the second row is created as follows. Data is cyclically read from the sequence permutation table t0 one by one at intervals of p1, and a sequence thus obtained is denoted as t1. For example, when data is circularly read from the sequence permutation table t0 at intervals of 3, the sequence t1 is obtained: t1: 1, 16, 7, 29, 49, ... The sequence t1 thus obtained is formed into sequence permutation table t1. The above process is also performed by computing the following (mathematically equivalent thereto): c(i) = (g 1 i )(mod83) where g1 is the primitive root obtained from g1 = (g 0 i )(mod83). Step S5: The table t1 obtained at step S4 is defined as a sequence permutation table which is referred to when the order of data arranged in the second row of the two-dimensional buffer is permuted. Step S6: Steps 4 and 5 are repeated by using p2, p3, p4, p5, p6 and p7, so that sequences t2 - t7 can be obtained. The sequences t2 - t7 are formed into tables t2 - t7, which are defined as sequence permutation tables referred to at the time of permutation of the third to eighth rows of the two-dimensional buffer. That is, steps 106 - 108 shown in Fig. 10 can be described as follows. First, prime numbers l i which satisfy the following are obtained where i = 2 - r, and r is the number of rows of the two-dimensional array: (i) (83-1, l i ) = 1 (82 and l i are mutually prime); and (ii) l i > 6. For r = 8, the prime numbers l 2 - l 8 are those to be obtained, and are respectively 7, 11, 13, 17, 19, 23 and 19 from the table shown in Fig. 6. Then, the pieces of data in the table to are cyclically read one by one at intervals of l i , so that the sequence permutation tables t2 - t7 can be created. In this case, the number "0" located at the end of the table to is excluded in the read operation. Step 7: The data arranged in the first through eighth rows of the two-dimensional buffer in which the blocks B 1 - B 8 are written are permuted in accordance with the sequence permutation tables t0 - t7. More particularly, the sequence of the data of the block B 1 is permuted in accordance with the sequence permutation table t0. Similarly, the sequence of the data of the block B 2 is permuted in accordance with the sequence permutation table t1. The same process is carried out for each of the blocks B 2 - B 7 . Finally, the sequence of the data of the block B 8 is permuted in accordance with the sequence permutation table t8. As described above, the sequence permutation tables t0 - t7 are created, and then the permutation process for the blocks B 1 - B 8 is performed. Alternatively, as shown in Fig. 10, the creation of the sequence permutation table and the permutation process are successively carried out for every block. The process of the second stage can be designed so that the sequence permutation tables are prepared and recorded beforehand and the permutation process refers to the tables recorded beforehand. (Second configuration of the interleaver 22) A description will now be given of the second configuration of the interleaver 22, which is the same as that above except for the second stage. The permutation process of the second stage is implemented by referring to tables created by steps S11 - S15 described below. Step S11: Step S11 is the same as the aforementioned step S1. That is, step S11 is to obtain the primitive root g 0 of the Galois field of the characteristic P and to create a table T0 described in the order of exponential notation of the primitive root in which the elements of the Galois field are expressed by anti-logarithms, which are arranged in the order of exponential notation. It will be noted that table T0 is the same as table t0 shown in Fig. 7B. Step S12: Step S12 is to obtain eight integers which have the mutually prime relationship with respect to the primitive root g 0 and are equal in number to the rows of the two-dimensional buffer. Let the eight numbers be denoted as q1, q2, q3, q4, q5, q6, q7 and q8. When the prime number P is equal to 83, for example, the primitive root is equal to 2 and the eight integers which are prime numbers are obtained as follows: 3, 5, 7, 11, 13, 17, 19 and 21 (except for 1 and 2). Step 13: The prime number q1 is added (modP) to each of the items of data of the table T0 obtained at step S12. When the table T0 defines the following sequence: T0: 1, 2, 4, 8, 16, ..., 42. 0, the resultant sequence for q1 = 3 is as follows: 4, 5, 7. 11, 19, ..., 3 Then, the anti-logarithm values thus obtained are converted into values in exponential notation. The following are obtained by a table of Fig. 7A that is the reverse operation of Fig. 7B: = 2, 27, 8, 24, ..., 7, 72 Sequence (7) thus obtained is formed in a table, which is used as a sequence permutation table T1 for the first row of the two-dimensional buffer. Similarly, step S13 is repeated using q2, q3, q4, q5, q6, q7 and q8, so that sequence permutation tables T2 - T8 respectively used for permutation of the second through eight rows of the two-dimensional buffer can be obtained. Step S15: The intra-permutation of the blocks B 1 - B 8 is performed in accordance with the sequence permutation tables T1 - T8. It is possible to prepare and record the sequence permutation tables T1 - T8 beforehand in a memory. (Third configuration) A description will be given, with reference to Fig. 9, of the third configuration of the interleaver 22. Referring to Fig. 9, an input sequence 80 consisting of 1140 bits is written into an interleaver 600 having a two-dimensional buffer of a 72 x 16 array. Then, each row of the 72 x 16 interleaver 600 is read for every 16 bits. The first row of the interleaver 600 is interleaved by using a 4 x 4 interleaver 610, and the second row thereof is interleaved by using a 6 x 3 interleaver 620. Similarly, the third row of the interleaver 600 is interleaved by an 8 x 2 interleaver 630. That is, the respective rows are interleaved by the different interleavers. Alternatively, an identical interleaver may be used to interleave each of the rows of the two-dimensional array. It is also possible to use an identical interleaver for some rows of the two-dimensional array. Then, the data thus interleaved are read in the longitudinal direction (0, 16, 32, 48, ...), so that an output data sequence 90 can be obtained. Since the last row of the buffer consists of only four bits, a 4 x 1 interleaver is applied thereto. However, the 4 x 4 or 2 x 2 interleaver may be used for the last row. The data of the last row can be read therefrom in the order of 1136, 1137, 1138 and 1139. However, in Fig. 9, the data of the last row are read in the reverse order, namely. 1139, 1138, 1137 and 1136. It is also possible to read data from the 72 x 16 interleaver except for the last row and to thereafter read the data of the last row and place them at given intervals. The interleaver 22 can be formed of any of the first through third configurations. After the interleaved data is produced by the interleaver 22, step 110 of Fig. 10, that is, the third state 43 shown in Figs. 5 and 8 is carried out. In the first and second configurations of the interleaver 22, the process at step 110 shown in Fig. 10 can be modified so as to prevent occurrence of a pattern which causes an error floor of the turbo codes. Fig. 11 is a graph showing an error floor. The error floor denotes a phenomenon in which the bit error rate (BER) is not improved as much as the S/N ratio is improved. In the graph of Fig. 11, the error floor starts to occur when the bit error rate is equal to 10 -7 to 10 -8 , and is not much improved even when the bit error rate is further improved. With the above phenomenon in mind, it is preferable to read the data from the two-dimensional array (buffer) in a fixed order but any of a plurality of predetermined orders. That is, it is preferable to determine the order of reading the data from the two-dimensional buffer after the intra-permutation on the basis of the value of the error floor. Thus, it is possible to suppress the occurrence of the error floor in the turbo codes. If the input sequence is divided into 10 blocks B 1 - B 10 , the data is read from the blocks B 10 , B 9 , B 8 , B 7 , B 6 , B 5 , B 4 , B 3 , B 2 and B 1 in this order. If the input sequence is divided into 20 blocks B 1 - B 20 , the data is read from the blocks B 19 , B 9 , B 14 , B 4 , B 0 , B 2 , B 5 , B 7 , B 12 , B 18 , B 16 , B 13 , B 17 , B 15 , B 3 , B 1 , B 6 , B 11 , B 8 , and B 10 in this order. It is also possible to read the data from B 19 , B 9 , B 14 , B 4 , B 0 , B 2 , B 5 , B 7 , B 12 , B 18 , B 10 , B 8 , B 13 , B 17 , B 3 , B 1 , B 16 , B 6 , B 15 and B 11 . One of the predetermined orders of reading the data is selected so that the occurrence of the error floor can be suppressed. As described above, when the input sequence is divided into 10 blocks, data written into the two-dimensional buffer is read therefrom in the reverse direction. This is simple and is implemented easily. (Puncture process) The conventional turbo encoder shown in Fig. 1A receives the K-bit input sequence and outputs the (3*K + T1 + T2)-bit encoded output where T1 is the number of tail bits output from RSC1, and T2 is the number of tail bits output from RSC2. In contrast, in the turbo encoder shown in Fig. 3, the bit addition process part 21 adds the a dummy bits to the N bits of the input sequence, so that the (N + a) bits are applied to the interleaver 22, RSC1 and RSC2. That is, 3a bits are auxiliary bits in total. In order to delete the 3a auxiliary bits, the puncturing process part 23 performs a puncturing process for the 3a auxiliary bits. A method of periodically deleting the redundant bits is generally used for the puncturing process suitable for the turbo codes. The above method can be applied to the puncturing process part 23. Thus, the output of the puncturing process part 23 receives (3K + 3a + T1 + T2) bits and outputs (3K + T1 and T2) to the next stage. Second Embodiment Fig. 12 is a block diagram of a turbo encoder according to a second embodiment of the present invention. In Fig. 12, parts that are the same as those shown in Fig. 3 are given the same reference numbers. The turbo encoder shown in Fig. 12 has the bit addition process part 21 that is placed at the input side of only the interleaver 22. That is, the encoded sequence X1 is the same as the input data sequence from an information source. Further, the RSC1 processes the input data sequence from the information source as per se. In order to delete the dummy bits added to the input data sequence by the bit addition process part 21, a pruning process part 123 is provided between the interleaver 22 and the RSC 2. The interleaver 22 has one of the aforementioned first through third configurations. Now, a fourth configuration of the interleaver 22 will be described. The fourth configuration of the interleaver 22 can be obtained by slightly modifying the first or second configuration. Such a modification is an improvement in obtaining the prime number, namely, the number of columns of the two-dimensional array. This will be described in detail below. The bit addition process part 21 operates as follows. At step (1), N IN is divided by 8, and the resultant value n is obtained. At step (2), the prime number P that is greater than n and closest to n is obtained. Further, (P - 1) and (P + 1) are prepared. Then, one of the numbers P, (P - 1) and (P + 1) that is equal to or greater than n and is closest to n is chosen. At step (3), the difference between 8 times of P and N IN is calculated, and the resultant value is denoted as a. At step (4), 4 bits (dummy bits) are added to N IN bits of the input sequence. An example for N IN = 660 will be described below. At step (1), n = 660/8 = 82.25 (the quotient is 82, and the remainder is 4). At step (2), the prime number P that is greater than 82.25 (= n) and closest thereto is 83 from the table shown in Fig. 4. Further, (P - 1) = 82, and (P + 1) = 84. Since the number 83 is greater than 82.25 and is closed thereto, the number 83 is chosen. At step (3). 83*8 = 664, and thus a = 4. That is, the number of dummy bits to be added to N IN is 4. At step (4), 4 dummy bits are added to the input sequence of 660 (=N IN ) bits. For example, the 4 dummy bits are added to the end of the 660-bit input sequence. The number (N IN + a) of bits thus obtained, that is, the number of (K + a) bits in Fig. 3 is divided by 8 without exception, and the quotient is always any of the prime number P, (P - 1) and (P + 1). If the prime number P is chosen at step (2), the sequence permutation tables can be created by the manner that has been described with reference to Fig. 5. However, if (P - 1) or (P + 1) is chosen, for example, if 82 or 84 is chosen, the two-dimensional array has 82 or 84 columns, and thus the sequence permutation tables associated with P = 83 are not used. The sequence permutation tables suitable for 82 or 84 columns can be created by modifying the aforementioned sequence permutation table t0 for 83 columns, as described below. Fig. 14A shows the sequence permutation table to for the first row of the two-dimensional array having 83 columns, which is the same as that shown in Fig. 7B. The sequence permutation table for the first row of the array having 82 columns (let t0 -1 be that table) is obtained by deleting "0" located at the end of the sequence of the table t0 for 83 columns. That is, the one-dimensional sequence t0 -1 for 82 columns is as follows: t0 -1 : 1, 2, 4, 8, 16, ..., 42. The table t0 -1 ranges from element 1 to element 82, and 1 is subtracted from each of all the elements. The resultant table ranges from element 0 to element 81, and is used as the sequence permutation table t0 -1 for the first row for 82 columns. The sequence permutation table for the first row of the array having 84 columns (let t0 -1 be that table) is obtained by adding the prime number P to the position next "0" located at the end of the sequence for 83 columns. That is, the one-dimensional sequence t0 -1 for 84 columns is as follows; t0 +1 : 1, 2, 4, 8, 16, ..., 42, 0, 83. By executing steps 106 - 108 shown in Fig. 10, sequence permutation tables t1 - t7, t1 -1 - t7 -1 , and t1 +1 - t7 -1 for the second through eight rows of the arrays having 83, 82 and 84 columns can be created, respectively. It is possible to create and record the above tables beforehand. The use of (P - 1) and (P + 1) makes it possible to decrease the difference between the number of the input data sequence and the number of bits (equal to the number of columns) processed by the interleaver 22 and to reduce the number of bits to be deleted by the pruning process. The process of the fourth configuration of the interleaver 22 shown in Fig. 12 can be applied to the interleaver 22 shown in Fig. 3. Another modification can be made. In the first and second embodiments of the present invention, the input data sequence is divided into the fixed number of blocks. This can be modified so that k different numbers of blocks for division can be used where k is an integer equal to or greater than 2. Thus, k interleavers are created, and one of them which brings the best performance is chosen and used. Let us consider a case where k = 10 and 20 and the input data sequence applied to the interleaver 22 consists of 640 bits. For k = 10, the 640-bit input data sequence is divided into 10 blocks, and the sequence permutation tables for the 10 blocks are labeled #1. For k = 20, the 640-bit input data sequence is divided into 20 32-bit blocks, and the sequence permutation tables for the 20 blocks are labeled #2. Then, one of the sets #1 and #2 of interleavers which provides the better bit error rate and/or the frame error rate is selected. The input data sequences consisting of different numbers of bits should be divided into different numbers of blocks in terms of the bit/frame error rate. That is, the number of blocks is adaptively changed taking into consideration the number of bits forming the input data sequence. Thus, it is possible to improve the performance of the turbo encoder. The present invention is not limited to the specifically described embodiments, variations and modifications, and other variations and modifications may be made without departing from the scope of the present invention.

Claims: An interleaving method which interleaves K bits of data by using an interleaver having a buffer arranged in a two-dimensional matrix, characterized in comprising the steps of: (a) generating a prime number P based on the number K of bits of data and the number N of rows of the interleaver so that the number of columns of the interleaver is determined by the prime number P; (b) permuting data in each of the N rows of the interleaver in accordance with a sequence permutation table generated specifically for each of the N rows, based on powers of a primitive root of the prime number P; (c) permuting the N rows based on a predetermined inter-permutation pattern; and (d) reading out data from the buffer in the column direction of the two-dimensional matrix. The interleaving method as claimed in claim 1, characterized in that said step (b) comprises the steps of: (1) generating sequence permutation data to in which elements of a Galois field of a characteristic P are arranged in an order of values of exponent parts of a power notation of a primitive root go of the Galois field; (2) generating (N-1) integers p 1 , p 2 , ..., p N-1 which are mutually prime with respect to (P-1); (3) generating first through (N-1)th sequence permutation data t i (i = 1 - (N-1)) by repeating a process for generating the ith sequence permutation data by cyclically reading data in the sequence permutation data to obtained in said step (1), at intervals of the ith integer p i obtained in said step (2); and (4) permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaving method as claimed in claim 1 characterised in that said step (b) comprises the steps of: (1) generating (N-1) integers p 1 , p 2 , ..., p N-1 which are mutually prime with respect (P-1), and generating (N-1) primitive roots g 1 , g 2 , ..., g N-1 each being arranged in a power notation of a primitive root go of a Galois field of a characteristic P by using the (N-1) integers p 1 , p 2 , ..., p N-1 ; (2) generating sequence permutation data t i (i = 0 - (N-1)), each of the sequence permutation data t i containing elements of the Galois field of the characteristic P arranged in an order of values of exponent parts of a power notation of the primitive root g i , for each of the first through Nth rows of the interleaver; and (3) permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaving method as claimed in claim 1 characterized in that said step (b) comprises the steps of: (1) generating sequence permutation data to in which elements of a Galois field of a characteristic P are arranged in an order of values of exponent parts of a power notation of a primitive root of the Galois field; (2) generating N integers q 1 , q 2 , ..., q N-1 which are mutually prime with respect to the primitive root; and (3) generating sequence permutation data t i (i = 1 - (N-1)), each of the sequence permutation data t i containing elements of the Galois field of the characteristic P arranged in an order of values of exponential parts of a power notation of the ith integer q i obtained said step (2), which integer is added to the ith element of the sequence permutation data to obtained in said step (1), and permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaving method as claimed in claim 1 characterized in that said step (a) comprises the steps of: (1) dividing the number K of the data bits by the number N of the rows to obtain a quotient; and (2) generating a prime number P that is greater than and nearest to the quotient so that the number of columns of the interleaver is determined by selecting one of P, (P-1) and (P+1) that is greater than and nearest to the quotient. The interleaving method as claimed in any of claims 1 to 5 characterized in that the number N of the rows of the interleaver is determined by selecting one of a plurality of predetermined integers that provides optimal transmission characteristics of the interleaved data. The interleaving method as claimed in any of claims 1 to 6 characterized in that , in said step (c), the order of rows is determined based on a value of an error floor in turbo codes in which transmission characteristics of the interleaved data are considered. The interleaving method as claimed in claim 7 characterized in that the number N of the rows of the interleaver is equal to 10, and the order of the rows determined in the step (c) is the order of the 10 th row, the 9 th row, the 8 th row, the 7 th row, the 6 th row, the 5 th row, the 4 th row, the 3 rd row, the 2 nd row, and the 1 st row. The interleaving method as claimed in claim 7 characterized in that the number N of rows of the interleaver is equal to 20, and the order of the rows determined in the step (c) is the order of the 20 th row, the 10 th row, the 15 th row, the 5 th row, the 1 st row, the 3 rd row, the 6 th row, the 8 th row, the 13 th row, the 3 rd row, the 19 th row, the 17 th row, the 14 th row, the 18 th row, the 16 th row, the 4 th row, the 2 nd row, the 7 th row, the 12 th row, the 9 th row, and the 11 th row. An interleaver for interleaving K bits of data by using a buffer arranged in a two-dimensional matrix, characterized in comprising: first means for generating a prime number P based on the number K of bits of data and the number N of rows of the interleaver so that the number of columns of the interleaver is determined by the prime number P; second means for permuting data in each of the N rows of the interleaver in accordance with a sequence permutation table generated specifically for each of the N rows, based on powers of a primitive root of the prime number P; third means for permuting the N rows based on a predetermined inter-permutation pattern; and fourth means for reading out data from the buffer in a column direction of the two-dimensional matrix. The interleaver as claimed in claim 10, characterized in that said second means comprises: means for generating sequence permutation data to in which elements of a Galois field of a characteristic P are arranged in an order of values of exponent parts of a power notation of a primitive root go of the Galois field; means for generating (N-1) integers p 1 , p 2 , ..., p N-1 ; which are mutually prime with respect to (P-1), and for generating first through (N-1)th sequence permutation data t 1 (i=1 - (N-1)) by repeating a process for generating ith sequence permutation data by cyclically reading data in the sequence permutation data to, at intervals of the ith integer p i ; and means for permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaver as claimed in claim 10, characterized in that said second means comprises: means for generating (N-1) integers p 1 , p 2 , ..., p N-1 which are mutually prime with respect (P-1), so that (N-1) primitive roots g 1 , g 2 , ..., g N-1 are generated, each of the primitive roots being arranged in a power notation of a primitive root go of a Galois field of a characteristic P by using the (N-1) integers p1, p2, ..., p N-1 , and for generating sequence permutation data t i (i = 0 - (N-1)), each of the sequence permutation data t i containing elements of the Galois field of the characteristic P arranged in an order of values of exponent parts of a power notation of the primitive root g i , for each of the first through Nth rows of the interleaver; and means for permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaver as claimed in claim 10, characterized in that said second means comprises: means for generating sequence permutation data to in which elements of a Galois field of a characteristic P are arranged in an order of values of exponent parts of a power notation of a primitive root of the Galois field; means for generating N integers q 1 , q 2 , ..., q N-1 which are mutually prime with respect to the primitive root, so that sequence permutation data t i (i = 0 - (N-1)) are generated, each of the sequence permutation data t i containing elements of the Galois field of the characteristic P are arranged in an order of values of exponential parts of a power notation of the ith integer q i , which integer is added to the ith element of the sequence permutation data t 0 ; and means for permuting data in each of the first through Nth rows of the interleaver in accordance with the sequence permutation data t i (i = 0 - (N-1)), respectively. The interleaver as claimed in claim 10, characterized in that the first means divides the number K of the data bits by the number N of the rows to obtain a quotient, and generates a prime number P that is greater than and nearest to the quotient so that the number of columns of the interleaver is determined by selecting one of P, (P-1) and (P+1) that is greater than and nearest to the quotient. The interleaver as claimed in any of claims 10 to 14, characterized in that the number N of the rows of the interleaver is determined by selecting one of a plurality of predetermined integers such that optimal transmission characteristics of the interleaved data are obtained with the number N of the rows equal to the selected integer. The interleaver as claimed in any of claims 10 to 15, characterized in that the order of the rows is determined by the third means based on a value of an error floor in turbo codes with respect to transmission characteristics of the interleaved data. A turbo encoding method comprising the steps of: (a) outputting a sequence of input data; (b) performing a convolutional encoding of the input data sequence to output a sequence of the encoded output data; (c) interleaving the input data sequence in accordance with the interleaving method as claimed in any of claims 1-9 to output a sequence of the interleaved data; and (d) performing a convolutional encoding of the interleaved data output by the step (c), to output a sequence of the encoded interleaved data. The turbo encoding method as claimed in claim 17 further comprising the steps of: (e) adding bits to the input data sequence prior to the interleaving step (c) if a number of bits of the input data sequence is less than a given number of bits, so that the number of bits of the input data sequence with the added bits is equal to the given number; and (f) deleting bits from the interleaved data sequence after the interleaving step (c) such that the number of the deleted bits is equal to the number of the added bits. The turbo encoding method as claimed in claim 18 wherein the bit adding step (c) employs dummy bits. A turbo encoder comprising: a first convolutional encoder performing a convolutional encoding of a sequence of input data to output a sequence of the encoded input data; an interleaving apparatus as claimed in any of claims 10 to 16, the interleaving apparatus arranged to interleave the sequence of the input data to output a sequence of the interleaved data; and a second convolutional encoder arranged to perform a convolutional encoding of the sequence of the interleaved data from the interleaving apparatus to output a sequence of the encoded interleaved data. The turbo encoder as claimed in claim 20 further comprising: means for adding bits to the input data sequence applied to the turbo encoder if a number of bits of the input data sequence is less than a given number of bits; and means for deleting bits from the sequence of the interleaved data so that an output data sequence of the turbo encoder has the same number of bits as the input data sequence. The turbo encoder as claimed in claim 21 wherein the means for adding bits employs dummy bits.

Inventors: SHIBUTANI AKIRA, SUDA HIROHITO, SHIBUTANI, AKIRA, SUDA, HIROHITO, Shibutani, Akira, NTT Mobile Communic. Network,Inc, Suda, Hirohito, NTT Mobile Communic. Network, Inc.
Applicants: NIPPON TELEGRAPH & TELEPHONE, NTT MOBILE COMMUNICATIONS NETWORK INC.
Application year: 2000
Main classification: H03M13
IPC codes: H03M13