WO2006040805A1 - 検査行列生成装置及び通信装置 - Google Patents

検査行列生成装置及び通信装置 Download PDF

Info

Publication number
WO2006040805A1
WO2006040805A1 PCT/JP2004/014958 JP2004014958W WO2006040805A1 WO 2006040805 A1 WO2006040805 A1 WO 2006040805A1 JP 2004014958 W JP2004014958 W JP 2004014958W WO 2006040805 A1 WO2006040805 A1 WO 2006040805A1
Authority
WO
WIPO (PCT)
Prior art keywords
matrix
check matrix
random number
parity check
order
Prior art date
Application number
PCT/JP2004/014958
Other languages
English (en)
French (fr)
Inventor
Wataru Matsumoto
Original Assignee
Mitsubishi Denki Kabushiki Kaisha
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 Mitsubishi Denki Kabushiki Kaisha filed Critical Mitsubishi Denki Kabushiki Kaisha
Priority to US11/664,456 priority Critical patent/US20080184084A1/en
Priority to CNA2004800441456A priority patent/CN101036300A/zh
Priority to EP04792226A priority patent/EP1806849A4/en
Priority to PCT/JP2004/014958 priority patent/WO2006040805A1/ja
Priority to JP2006540790A priority patent/JP4772689B2/ja
Publication of WO2006040805A1 publication Critical patent/WO2006040805A1/ja

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/63—Joint error correction and other techniques
    • H03M13/635—Error control coding in combination with rate matching
    • H03M13/6362—Error control coding in combination with rate matching by puncturing
    • H03M13/6368—Error control coding in combination with rate matching by puncturing using rate compatible puncturing or complementary puncturing
    • H03M13/6393—Rate compatible low-density parity check [LDPC] codes
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/033—Theoretical methods to calculate these checking codes
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1148—Structural properties of the code parity-check or generator matrix
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1148—Structural properties of the code parity-check or generator matrix
    • H03M13/116—Quasi-cyclic LDPC [QC-LDPC] codes, i.e. the parity-check matrix being composed of permutation or circulant sub-matrices
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1148—Structural properties of the code parity-check or generator matrix
    • H03M13/118—Parity check matrix structured for simplifying encoding, e.g. by having a triangular or an approximate triangular structure

Definitions

  • the present invention provides a low-density parity check (LDPC) as an error correction code.
  • LDPC low-density parity check
  • the present invention relates to a parity check matrix generator when a parity-check code is used and a communication device equipped with this parity check matrix generator, and in particular, a test that can search for a deterministic and stable check matrix for LDPC codes.
  • the present invention relates to a matrix generation device and a communication device equipped with this check matrix generation device.
  • the transmission-side communication device includes an encoder and a modulator, while the reception-side device includes a demodulator and a decoder.
  • the flow of encoding and decoding when LDPC codes are used will be described.
  • the encoder on the transmission side generates a check matrix H by a conventional method to be described later. Then, the generator matrix G is obtained based on the following conditions.
  • the encoder receives a message (m m to m) of information length K and generates the above
  • the modulator applies BPSK (Binary Phase Shift) to the generated codeword C.
  • BPSK Binary Phase Shift
  • the demodulator performs digital demodulation such as BP SK, QPSK, and multi-level QAM on the modulated signal received via the communication channel, and the decoder further performs the LDPC code.
  • the demodulated demodulation results are iteratively decoded using the “sum-product algorithm”, and the estimation results (corresponding to the original mm to m) are output.
  • FIG. 1 is a diagram showing a parity check matrix for a conventional LDPC code proposed by Non-Patent Document 1, for example.
  • the matrix shown in Fig. 1 is a binary matrix of "1" and "0", and the "1" part is painted over. All other blank parts are “0”.
  • the number of “1” in one row is 4 (this is expressed as the order of the row), the number of “1” in 1 column (this is expressed as the order of the column) is 3, and all This is generally called “Regular—LDPC code” because the order of columns and rows is uniform.
  • the matrix is divided into three blocks, and random replacement is performed on the second block and the third block.
  • Non-Patent Document 2 proposed a method using a Euclidean geometric code as an LDPC code that can generate a deterministic deterministic without using a computer search and exhibits relatively stable and good characteristics. .
  • This method is described as “Regular—LDPC code” composed of regular ensembles!
  • Non-Patent Document 2 a method for generating a parity check matrix of an LDPC code using Euclidean geometric code EG (2, 2 6), which is a kind of finite geometric code, has been proposed. At 4 points, characteristics approaching 1.45 dB from the Shannon limit are obtained.
  • FIG. 2 is a diagram showing the configuration of the Euclidean geometric code EG (2, 2 2 ), for example, and has a “Regular—LDPC code” structure in which the respective degrees of the rows and columns are 4 and 4.
  • K 2 2 s — 3 s
  • the Euclidean geometric code has a structure in which the arrangement of “1” in each row is cyclically shifted for each row, and the code is easily and deterministically configured.
  • the row and column orders are further changed based on the Euclidean geometric code, and the rows and columns are expanded as necessary. For example, when the order of the column of EG (2, 2 2 ) is separated into 1Z2, in Non-Patent Document 2, four orders of four in one column are separated every other two.
  • Fig. 3 is a diagram showing an example in which the order of a sequence of Euclidean geometric codes is regularly separated into 4 as many as 2 forces.
  • Non-Patent Document 3 reports that the characteristics of “Irregular to LDPC code” are better than the characteristics of “Regular-LDPC code”. It was theoretically analyzed by Non-Patent Document 4 and Non-Patent Document 5. Note that the above “Irregular to LDPC code” represents an LDPC code in which the order of columns and rows is not uniform or either.
  • Non-Patent Document 5 analyzed the LDPC code “Sum—Product algorithm” on the assumption that the log-likelihood ratio (LLR) of the input and output of the iterative decoder can approximate a Gaussian distribution. Find the ensemble of the correct row and column order.
  • LLR log-likelihood ratio
  • Non-Patent Document 1 R.G. Gallager. "Low-Density Parity Check Codes.” IRE Trans. On IT pp21-28 January 1962.
  • Non-Patent Document 2 Y. Kou, S. Lin, and M. P. C. Fossorier, "Low Density Parity Check
  • Non-Patent Document 3 M. G. Luby, M. Mitzenmacher, M. A. Shokrollahi, and D. A.
  • Non-Patent Document 4 TJ Richardson and R. Urbanke, i'he capacity of low-density parity-check codes under message-passing decoding, "IEEE Trans. Inform. Theory, vol.47, No.2, pp.599— 618, Feb. 2001.
  • Non-Patent Document 5 S.- Y. Chung, T. J. Richardson, and R. Urbanke, "Analysis of
  • the minimum value of the check matrix cycle (referred to as the inner diameter) is often about 6 (Figs. 2 and 2).
  • the longer the check matrix cycle the wider the Hamming distance of the codeword, and the general analysis results that the characteristics of Sum-product decoding are good indicate that the error correction capability is There was a problem that the performance capability was insufficient.
  • the present invention has been made to solve the above-described problems, and is deterministic, has stable characteristics, has good error correction capability, and supports various ensembles, code lengths, and coding rates.
  • the purpose of the present invention is to obtain a parity check matrix generator for LDPC codes and a communication device equipped with this parity check matrix generator, which can easily generate a parity check matrix for Irregular-LDPC codes in a short time.
  • deterministic has a common algorithm for both transmission and reception, and derives the same result with limited specifications.
  • limited specifications such as coding rate and code length are used. This means that it is possible to generate a parity check matrix for LDPC codes.
  • the parity check matrix generation apparatus includes an information length 'coding rate' column maximum order input unit that inputs a predetermined information length, a coding rate, and a maximum order of the column; A parameter calculation unit for a pseudo random number replacement matrix that calculates parameters for a pseudo random number replacement matrix, which will be described later, using the information length, the coding rate, and the maximum order of the column, and the calculated pseudo random number replacement matrix Using the pseudo random number substitution matrix generating unit that generates the pseudo random number substitution matrix using the pseudo random number sequence and the Latin square matrix according to the parameters, and the input predetermined information length, coding rate, and maximum order of the column, A degree distribution optimization parameter calculation unit for calculating parameters of combinations of possible degree distributions for optimizing the order distribution of a check matrix generation configurable using the pseudo-random number permutation matrix; and Up It is generated by the order distribution optimization unit that optimizes the order distribution of the parity check matrix generation by the density evolution method, and the optimized order distribution of the parity check matrix generation, with
  • FIG. 1 is a diagram showing a parity check matrix for a conventional LDPC code.
  • FIG. 2 is a diagram showing a configuration of Euclidean geometric code EG (2, 2 2 ).
  • FIG. 3 is a diagram showing an example in which the order of a sequence of Euclidean geometric codes is regularly separated into 4 as many as 2 forces.
  • FIG. 4 is a block diagram illustrating a configuration of an LDPC code Z decoding system.
  • FIG. 5 is a block diagram showing a configuration of a parity check matrix generation device according to Embodiment 1 of the present invention.
  • FIG. 6 is a flowchart showing a process flow of the parity check matrix generation apparatus according to Embodiment 1 of the present invention.
  • FIG. 7 is a block diagram showing a configuration of a parity check matrix generation device according to Embodiment 2 of the present invention.
  • FIG. 8 is a flowchart showing a processing flow of a parity check matrix generation device according to Embodiment 2 of the present invention.
  • FIG. 9 is a diagram showing bit error rate characteristics and frame error rate characteristics with respect to the signal power to noise power ratio per bit of information.
  • FIG. 10 is a block diagram showing a configuration of a parity check matrix generation device according to Embodiment 3 of the present invention.
  • FIG. 11 is a flowchart showing a process flow of the parity check matrix generation device according to the third embodiment of the present invention.
  • FIG. 12 is a block diagram showing a configuration of a parity check matrix generation device according to Embodiment 4 of the present invention. is there.
  • FIG. 13 is a flowchart showing a process flow of a parity check matrix generation device according to Embodiment 4 of the present invention.
  • FIG. 14 is a diagram showing an example of a multi-edge type ensemble.
  • FIG. 15 is a diagram showing an example of an irregular LDPC code.
  • FIG. 16 is a diagram illustrating an example of a parity check matrix configured according to the degree distribution of FIG. 15.
  • FIG. 17 is a diagram showing an example of the degree distribution of an irregular LDPC code used for evaluation.
  • FIG. 18 is a diagram showing a frame error rate characteristic of an irregular LDPC code.
  • FIG. 19 is a block diagram showing a configuration of a parity check matrix generation device according to Embodiment 5 of the present invention.
  • FIG. 20 is a block diagram showing a configuration of a communication apparatus according to Embodiment 6 of the present invention.
  • FIG. 4 is a block diagram showing the configuration of the LDPC code decoding system.
  • the communication device on the transmission side includes an encoder 101 and a modulator 102
  • the communication device on the reception side includes a demodulator 104 and a decoder 105.
  • a description is given of the flow of decoding and decoding when using the generated parity check matrix for LDPC codes.
  • a check matrix H generated by a check matrix generation device for LDPC codes is used as described later.
  • Encoder 101 obtains generator matrix G based on the following conditions.
  • modulator 102 performs digital modulation such as BPSK, QPSK, and multi-level QAM on the generated codeword C and transmits the codeword C.
  • the generator matrix G must be generated by the Gaussian elimination method or the like for the parity check matrix H force, and the parity check matrix H is low density, whereas the generator matrix G
  • the ratio of “1” in the binary generator matrix is large with a relatively high density, there is a problem that a large amount of memory is required to store the generator matrix G.
  • the following literature proposes a method of constructing a matrix having a lower triangular structure (or upper triangular structure) in the check matrix H.
  • the encoder 1 01 When C satisfies / C r _0 and the information message "-(", "" 2 , ..., " ⁇ ) is given, the encoder 1 01 generates a parity element according to the following equation. This procedure makes it possible to easily realize the sign ⁇ without using the generator matrix G. In this invention, it is necessary to construct a parity check matrix H including a lower triangular matrix.
  • the communication device on the receiving side receives information about necessary specifications from the communication device on the transmitting side in advance and generates a check matrix, and the demodulator 104 receives it via the communication path 103.
  • the modulated signal is digitally demodulated such as BPSK, QPSK, multi-level QAM, etc., and the decoder 105 uses the generated parity check matrix to calculate the sum-product for the LDPC encoded demodulation result. Performs iterative decoding using the “algorithm” and outputs the estimation result (corresponding to the original u U... !!).
  • Non-Patent Document 5 a conventional parity check matrix generation method for “Irregular—LDPC code” theoretically analyzed by Non-Patent Document 5 will be described in detail.
  • the decoding characteristics of the sum-product decoding method are calculated by iteratively calculating the probability density function of the log-likelihood ratio (LLR) of the input and output in the decoder 105 according to the decoding process of the sum-product decoding method.
  • LLR log-likelihood ratio
  • the density evolution method which is a parity check matrix generation method for LDPC codes described in Non-Patent Document 5 above, and the log-likelihood ratio (LLR) of the input and output in the decoder 105 can approximate a Gaussian distribution.
  • the check matrix is represented by a Tanner graph, and each column is defined as a variable node and each row is defined as a check node. Also, if there is a “1” at the intersection of a certain column and a certain row, the noble node and the check node are connected by “branches”!
  • LLR message propagation from the check node to the variable node is analyzed.
  • s m is the average value of uO
  • uO is passed through a transmission line containing Gaussian noise with variance ⁇ 2 and uO n
  • ⁇ p i represents the ratio of branches belonging to variable nodes and check nodes of degree i, respectively.
  • D is the order of the maximum noble node, and d is the order of the maximum check node.
  • ⁇ (X) in the above equation (1) is defined as the following equation (4).
  • R represents the set of all real numbers
  • u represents the ensemble average of the LLR output of the variable node.
  • T is the ensemble of the LLR output value of the check node at the first iteration
  • the condition for determining the threshold of SNR where the error can be 0 is t (s) ⁇ ⁇ (represented as R +) when 1 ⁇ ⁇ .
  • the condition for obtaining the SNR threshold where the error can be 0 is r (s) ⁇ 0
  • step ST101 and step ST102 are repeatedly executed, the generation function (X) of the order distribution of the noble node and the generation function ⁇ (X) of the order distribution of the check node Search for better ensembles with linear programming.
  • the signal power is normalized to 1 and the SNR limit (threshold) at which the error can be 0 is obtained from the following equation (10).
  • threshold (dB) -10 * loglOI2 * ⁇ ⁇ )---(l ⁇ )
  • Non-Patent Document 5 the parity check matrix obtained by the maximum value of the code rate becomes fluid, and the code rate that is fixed as a design specification fluctuates. There is a problem. Further, in Non-Patent Document 5, the derivation of the generation function (X) of the order distribution of the noble node and the derivation of the generation function ⁇ (X) of the order distribution of the check node are repeatedly performed for a predetermined number of times. Therefore, if a certain amount of time is required for the search process, it is difficult to easily cope with various problems, various ensembles, code lengths, and code ratios. There is also a title.
  • FIG. 5 is a block diagram showing a configuration of a parity check matrix generation apparatus for LDPC codes according to Embodiment 1 of the present invention.
  • This parity check matrix generation apparatus includes an information length 'coding rate' column maximum order input unit 11, a pseudo random number substitution matrix parameter calculation unit 12, a pseudo random number substitution matrix generation unit 1 3, a degree distribution optimization parameter calculation unit 14, An order distribution optimization unit 15 and a parity check matrix generation unit 16 are provided.
  • This check matrix generation device is mounted in the transmission side and reception side communication devices in FIG. 4, and the generated check matrix may be output to the encoder 101 and the decoder 105, or the transmission side and the reception side.
  • the generated check matrix may be stored in a memory (not shown) in the communication device on the transmission side and the reception side, and output to the encoder 101 and the decoder 105.
  • FIG. 6 is a flowchart showing a processing flow of the parity check matrix generation apparatus according to Embodiment 1 of the present invention.
  • step ST11 information length 'coding rate' sequence maximum order input section
  • step ST12 the pseudorandom number substitution matrix parameter calculation unit 12 calculates a parameter for the pseudorandom number substitution matrix using the predetermined information length, coding rate, and maximum order of the column input in step ST11. To do.
  • step ST13 the pseudo random number substitution matrix generation unit 13 generates a pseudo random number substitution matrix using the pseudo random number sequence and the Latin square matrix according to the parameters for the pseudo random number substitution matrix calculated in step ST12. .
  • step ST14 the order distribution optimization parameter calculation unit 14 can be configured using a pseudo-random permutation matrix using the predetermined information length, coding rate, and maximum order of columns input in step ST11 above. Check Calculates the possible order distribution combination parameters for optimization of matrix generation order distribution.
  • step ST15 the degree distribution optimizing unit 15 uses the combination of possible degree distributions calculated in step ST14 as a constraint to generate a check matrix by the density evolution method. Optimize the order distribution.
  • step ST16 the parity check matrix generation unit 16 generates the parity check matrix for the LDPC code by arranging the pseudorandom substitution matrix generated in step ST13 according to the order distribution of the parity check matrix generation optimized in step ST15. To do.
  • P 2 X p 2 (p + 1 is a prime number) pseudo-random permutation matrix SP parameter p is calculated as follows.
  • step ST13 the pseudo random number substitution matrix generation unit 13 performs the procedure of steps ST13-1 to ST13-6 shown below using the parameter p for the pseudorandom number substitution matrix SP calculated in step ST12. As shown, a pseudorandom permutation matrix SP of p 2 X p (p + 1 is a prime number) is generated using a pseudorandom sequence and a Latin square matrix.
  • V (P) is a continuous power value of (for example,,, ..., (7, '-indicates all non-zero elements of GF (P) if This sequence is a pseudo-random sequence in which the differences between elements are non-uniform and there are no overlapping elements.
  • a substitution pattern S ′) shown below is generated for the pseudo-random number sequence C).
  • () can be used to generate p different Latin square matrices of Ll (q, /) 2 (q), ..., L flesh(?
  • ⁇ ( ⁇ ? i) is a matrix with a row / column for ⁇ , where p is added to each element in the row for each increment of each row.
  • the matrix of ⁇ / columns with each element of for is also a Latin square. This operation also increases the distance between elements in the same column of a row / column matrix generated from flP ,.
  • Equation 17 The basic permutation matrix pattern BP, (q,,.) Is generated as follows.
  • the distance between adjacent elements in this series is increasing. By increasing the distance between adjacent elements, it is possible to suppress the occurrence of short cycles, that is, cycles 4 and 6, which occur with the lower triangular matrix of T described in the first embodiment.
  • Equation 19 A plurality of basic permutation sequences SPj , (a two-element two- random pseudorandom permutation matrix.,-(J defined by the following equation is generated.
  • step ST14 the order distribution optimization parameter calculation unit 14 uses the predetermined information length K, the sign rate rate, and the maximum degree d of the sequence input in step STil to simulate the maximum degree d. Calculates the possible order distribution combination parameters for optimizing the order distribution of the parity check matrix that can be configured using a similar random number permutation matrix.
  • the combinations of possible degree distributions are calculated under the constraints of the maximum degree d of the column and the check matrix including the lower triangular matrix.
  • 5 ⁇ means any ⁇ .
  • the sequence of ⁇ arranged on the diagonal line means that 1 is arranged in a staircase pattern like ⁇ from the top row to the bottom row of H.
  • the constraint conditions determined here are the size of the pseudorandom substitution matrix SP, its maximum number of arrangements, and arrangement locations.
  • a combination of possible degree distributions is calculated under the constraints of the check matrix including the maximum order and the lower triangular matrix.
  • the degree distribution indicates, for example, the ratio of the number of columns of column degree 4 to the total number of columns of the check matrix, and a pseudo-random substitution matrix (column weight 1, row weight) so that the degree of column becomes 4 according to the ratio. Create 4) by combining 4) in the column direction.
  • the ratio that the column weight 4 can be obtained using the pseudorandom permutation matrix with respect to the total number of columns of the check matrix H is 0-4Z8, and the ratio that the column weight 3 can take is 0-5Z8.
  • the combination of these column orders and row orders and possible ratios is output as a parameter for optimization of the order distribution performed in the next step ST15.
  • step ST15 the order distribution optimizing unit 15 optimizes the order distribution of the parity check matrix generation by the density evolution method as shown below with the combination of possible order distributions calculated in step ST14 as a constraint.
  • the optimal degree distribution is calculated from the combinations of the degree distributions obtained in step ST14.
  • the degree distribution indicates, for example, what percentage of the total check matrix the number of columns of the column order 4 is, and the pseudorandom substitution matrix (column weights) so that the column order becomes 4 according to the ratio. 1. Create a combination of four row weights 1) in the column direction. As a result, a check matrix can be generated in the next step ST16.
  • the pseudorandom number in the check matrix H is determined based on the size of the pseudorandom substitution matrix SP in the check matrix H obtained in step ST14, the maximum number of placement, and the constraints on the placement location.
  • Arbitrary pseudorandom permutation matrix SP is placed at the place of permutation matrix SP, or the combination can be selected under the condition that 0 matrix of p 2 X p 2 is placed at that place.
  • an ensemble (order distribution) of the “Irregular—LDPC code” based on the required code rate is obtained using optimization in the Gaussian approximation method of the density evolution method.
  • the generation function of the column order distribution (X) and the generation function p (X) of the row order distribution are treated simultaneously as variables, and the optimal generation function is determined by linear programming so that the Gaussian noise ⁇ is maximized. Search for (X) and the generator function ⁇ ) (see equation (11) below). The constraints in this search are In other words, the following equation (12) is satisfied.
  • the values that can be taken by the variable of the generation function (X) of the column order distribution and the generation function p (X) of the row order distribution are only the values calculated by the constraint conditions of the check matrix H. .
  • the generation function (X) and the generation function ⁇ (X) satisfying a predetermined condition are obtained by a single linear design method as described above, the generation function is set as in Non-Patent Document 5.
  • An ensemble with deterministic and stable characteristics can be generated easily and in a shorter time than the method of repeatedly calculating (X) and calculating the generation function ⁇ (X) and obtaining the optimum values of both.
  • parity check matrix generation section 16 arranges the pseudo-random number permutation matrix SP generated in step ST13 as shown below by the order distribution of the parity check matrix generation optimized in step ST15. To generate a parity check matrix for the LDPC code. In other words, here, based on the degree distribution obtained in step ST15 above, in the check matrix ⁇ Determine the placement of the pseudorandom permutation matrix SP.
  • Irregular to LDP that is deterministic, has stable characteristics, has good error correction capability, and supports various ensembles, code lengths, and code ratios. If the parity check matrix for C code can be easily generated in a short time, the effect can be obtained.
  • the first embodiment is extended to explain how to generate a parity check matrix for RC (Rate comapat3 ⁇ 4le) —LDPC code that is compatible with multiple code rates (only the coding rate changes without changing the information length) To do.
  • FIG. 7 is a block diagram showing a configuration of a parity check matrix generation apparatus for LDPC codes according to Embodiment 2 of the present invention.
  • This parity check matrix generation apparatus includes an information length 'coding rate' column maximum order input unit 21, a pseudo random number substitution matrix parameter calculation unit 22, a pseudo random number substitution matrix generation unit 23, a degree distribution optimization parameter calculation unit 24, An order distribution optimization unit 25 and a check matrix generation unit 26 are provided.
  • This check matrix generation device is mounted in the transmission side and reception side communication devices in FIG. 4, and the generated check matrix may be output to the encoder 101 and the decoder 105, or the transmission side and the reception side.
  • the generated check matrix may be stored in a memory (not shown) in the communication device on the transmission side and the reception side, and output to the encoder 101 and the decoder 105.
  • FIG. 8 is a flowchart showing the flow of processing of the parity check matrix generation apparatus according to Embodiment 2 of the present invention.
  • information length “coding rate” sequence maximum order input section 21 is input with a predetermined information length, a maximum order of a sequence of maximum code rate excluding a plurality of code rate and unsigned rate.
  • the pseudo random number substitution matrix parameter calculation unit 22 uses the predetermined information length input in step ST21, the maximum order of the sequence of the maximum coding rate excluding a plurality of coding rates and uncoded symbols. Calculate the parameters for the pseudorandom substitution matrix.
  • step ST23 the pseudo random number substitution matrix generation unit 23 generates a pseudo random number substitution matrix using the pseudo random number sequence and the Latin square matrix according to the parameters for the pseudo random number substitution matrix calculated in step ST22.
  • step ST24 the order distribution optimization parameter calculation unit 24 uses the predetermined information length input in step ST21, the multiple coding rates, and the maximum order of the maximum coding rate sequence excluding non-coding. Then, the parameters of the possible combinations of order distributions for optimizing the order distribution of the check matrix generation that can be configured using the pseudo random number substitution matrix are calculated.
  • step ST25 the order distribution optimizing unit 25 uses the combination of possible order distributions calculated in step ST24 as a constraint condition to calculate the order distribution of the parity check matrix generation for a plurality of coding rates by the density evolution method.
  • step ST26 the check matrix generation unit 26 arranges the pseudorandom permutation matrix generated in step ST23 according to the order distribution of the check matrix generation for the plurality of coding rates optimized in step ST25, and performs LDPC. Generate a parity check matrix.
  • step ST21 the information length 'coding rate' sequence maximum order input unit 21 receives a predetermined information length K, a plurality of encoding rates R (1), and the maximum order d of the sequence of the maximum encoding rate excluding non-coding. input.
  • 1 1, 2,..., 0 ⁇ R (1) ⁇ R (2),.
  • This step ST21 differs from the step ST11 of the first embodiment in that a plurality of coding rates R (1) and the maximum degree d of the maximum coding rate column excluding uncoding are input. Yes.
  • step ST22 the pseudo random number substitution matrix parameter calculation unit 22 obtains the maximum code sign excluding the predetermined information length K input in step ST21, a plurality of coding rates R (l), and no coding. Using the maximum degree d of the percentage column, calculate the parameters for the pseudorandom permutation matrix as shown below.
  • step ST23 the pseudo random number permutation matrix generation unit 23 performs the same process as the process of step ST13 performed by the pseudo random number permutation matrix generation unit 13 of the first embodiment. That is, the pseudo random number substitution matrix generation unit 23 performs the processing of step ST13-1—step ST13-6.
  • step ST24 the order distribution optimization parameter calculation unit 24 obtains the maximum code excluding the predetermined information length K, the plurality of code rate R (l) and the unsigned value input in step ST21.
  • a combination of possible order distributions for optimizing the order distribution of a parity check matrix that can be configured using a pseudo-random permutation matrix using Parameters are calculated.
  • step ST25 the degree distribution optimization unit 25 uses a combination of the possible degree distributions calculated in step ST24 as a constraint, and a check matrix for a plurality of coding rates by the density evolution method as shown below. Optimize the generation order distribution.
  • Equation 24 Furthermore, let ⁇
  • A be a parity check matrix of AixW for _ ⁇ .
  • H _ 2) be a parity check matrix of (M) to X (Nil) defined by the following equation consisting of a 0 matrix of ⁇ ⁇ — and ⁇ ⁇ ⁇ and an additional parity check matrix of o.
  • Equation 27 the degree distribution of ⁇ ⁇ (, ⁇ ⁇ ⁇ that minimizes ⁇ w) is found.
  • the difference between the SNR of the iteration threshold of H estimated by the Gaussian approximation method and the Shannon limit is expressed as (iB It is expressed by the following formula.
  • parity check matrix generation section 26 uses the pseudorandom permutation matrix generated in step ST23 according to the order distribution of the parity check matrix generation for the plurality of code ratios optimized in step ST25. To generate a check matrix for the LDPC code as shown below.
  • R (3) 1-M
  • FIG. 6 is a diagram showing frame error rate (FER) characteristics.
  • the information length is 1536 bits
  • the communication channel is AWGN
  • the maximum number of repetitions is 200.
  • the decoding method is “Sum-Product algorithm”. This characteristic is represented by the parameters M-M,
  • the each of the difference BER 10- 4 points, about 1.8 (dB), 1.7 (dB ) ⁇ 1.45 (dB), and the information length 1536 RC approaching bits despite the Shannon limit is not so long — LDPC code has been realized.
  • FIG. 10 is a block diagram showing a configuration of a parity check matrix generation apparatus for LDPC codes according to Embodiment 3 of the present invention.
  • This check matrix generation device includes an information length 'coding rate' column maximum order input unit 31, a pseudo random number substitution matrix parameter calculation unit 32, a pseudo random number substitution matrix generation unit 33, an order distribution optimization parameter calculation unit 34, An order distribution optimization unit 35 and a check matrix generation unit 36 are provided.
  • This check matrix generation apparatus is mounted in the transmission side and reception side communication apparatuses in FIG. 4, and the generated check matrix may be output to the encoder 101 and the decoder 105, or the transmission side and the reception side may be output.
  • the generated check matrix may be stored outside the communication device on the transmission side and the reception side, and output to the encoder 101 and the decoder 105.
  • FIG. 11 is a flowchart showing a process flow of the parity check matrix generation device according to Embodiment 3 of the present invention.
  • the information length 'coding rate' sequence maximum order input unit 31 inputs a predetermined information length, code rate, and sequence maximum order.
  • the pseudo-random number substitution matrix parameter calculation unit 32 uses the predetermined information length, coding rate, and maximum order of the columns input in step ST31, and first and second pseudo simulations described later. Calculate the parameters for the random number substitution matrix.
  • step ST33 the pseudo random number substitution matrix generation unit 33 uses the pseudo random number matrix and the Latin square matrix according to the parameters for the first and second pseudo random number substitution matrices calculated in step ST32.
  • a first pseudo random number permutation matrix is generated, and the first pseudo random number permutation matrix is further compressed in the row direction to generate a second pseudo random number permutation matrix.
  • the order distribution optimization parameter calculation unit 34 uses the predetermined information length, coding rate, and maximum order of the columns input in step ST31 to generate the first and second pseudorandom substitution matrices. Is used to calculate possible order distribution combination parameters for optimizing the order distribution of the configurable check matrix.
  • step ST35 the order distribution optimizing unit 35 optimizes the order distribution of the parity check matrix generation by the density evolution method using the combination of possible order distributions calculated in step ST34 as a constraint.
  • step ST36 the parity check matrix generation unit 36 arranges the first and second pseudorandom number permutation matrices generated in step ST33 according to the order distribution of the parity check matrix generation optimized in step ST35, and sets the LDPC code. Check matrix for Is generated.
  • Pseudorandom number replacement matrix size P ⁇ ⁇ M / (Pseudorandom number replacement matrix
  • step ST33 the pseudorandom permutation matrix generation unit 33 calculates the first and second pseudorandom permutation matrixes calculated in step ST32 above.
  • the first pseudorandom permutation matrix is generated using the pseudorandom number sequence and the Latin square matrix with the parameters for Generate. That is, the first pseudorandom substitution matrix from J? And b found in step ST 32 above
  • P is a prime number such that / >> p and is a primitive element of.
  • This sequence C, withdraw ,, (,) is a pseudo-random sequence with non-uniform differences between elements ⁇ .
  • the following formula generates a Latin square matrix i; ( q , l) whose elements are pseudo-random number sequences and their replacement patterns.
  • Equation 35 The following formula generates a matrix 83 ⁇ 4, (, 0 with a constant interval per row of the basic permutation matrix pattern HP (q, i), where / is the interval for per-row permutation. I) If ft and top are disjoint,
  • Equation 36 The basic permutation matrix pattern BP h by the following equation. Generate basic permutation sequence BPS (,) from p (q, i). ii-1,2, .., p.
  • This operation generates a basic permutation matrix for generating a matrix of column degree 1 and row degree 1 of wa xp.fc.
  • the short cycle generated with the lower triangular matrix of T described in Embodiment 1 can be suppressed. Also, by the distance between the elements non-uniform, the generation of short cycles in the above d p. In the basic permutation matrix p .fc can be suppressed. Also, check the basic permutation matrix of p.fcxp.h generated to avoid the occurrence of cycle 4, and if cycle 4 exists, p, b, G. , Init, and hop are varied, and combinations that do not have cycle 4 are searched and regenerated. However, in the procedure of the present embodiment, cycle 4 is rarely present and rarely needs to be regenerated.
  • ii), i ⁇ i), and iv) are the row and column permutations based on SP originate,,-/, all of which are p.bxp.b matrixes with a column degree of 1 and a row degree of 1. By combining these matrices, a matrix with column order ⁇ X row order 1 to column order 4 X row order 4 can be generated.
  • Four different types of pseudorandom permutation matrices can be generated by this row and column permutation.
  • Step ST 3 Generate (pi>) x ( P .fc / v) pseudorandom permutation matrix c with column degree cw and row degree 1 using basic permutation matrix sra _ / obtained in step 7 by the following procedure .
  • the column is compressed to I / cw, that is, compressed in the row direction.
  • c to DPS, ..
  • J is the identity matrix of (p 'fc /).
  • cycle 4 When two or more combinations of are arranged in the column direction, it is known that cycle 4 will occur, so it shall be arranged in one column only in the row direction.
  • step ST34 the order distribution optimization parameter calculation unit 34 uses the predetermined information length, code rate, and maximum order of the sequence input in step ST31 as follows. Then, the parameters of the combination of possible order distributions for optimizing the order distribution of the check matrix generation that can be configured using the first and second pseudorandom substitution matrices are calculated.
  • step ST35 the order distribution optimizing unit 35 uses the density evolution method with the combination of possible order distributions calculated in step ST34 as a constraint, as in step ST15 of the first embodiment. Optimize the order distribution of the check matrix generation.
  • parity check matrix generation section 36 arranges the first and second pseudo-random number permutation matrices generated in step ST33 according to the order distribution of the parity check matrix generation optimized in step ST35. Generate a parity check matrix for LDPC codes.
  • the first pseudorandom permutation matrix ⁇ precede,, _ I, BPS h ⁇ till, —II, BPS, —III, using the optimal combination of degree distributions obtained in step ST 3 5 Generate a parity check matrix H using BPS-IV and the second pseudorandom substitution matrix c and.
  • a good arrangement is obtained by checking the number of cycles by searching.
  • FIG. 12 is a block diagram showing a configuration of a parity check matrix generation apparatus for LDPC codes according to Embodiment 4 of the present invention.
  • This check matrix generation device includes an information length 'coding rate' column maximum order input unit 41, a pseudo random number substitution matrix parameter calculation unit 42, a pseudo random number substitution matrix generation unit 43, an order distribution optimization parameter calculation unit 44, An order distribution optimization unit 45 and a check matrix generation unit 46 are provided.
  • This parity check matrix generator is configured to transmit and receive in FIG.
  • the generated check matrix may be output to the encoder 101 and the decoder 105 installed in the communication device on the transmission side, or installed outside the communication device on the transmission side and the reception side.
  • the data may be stored in a memory (not shown) in the communication device on the transmission side and the reception side and output to the encoder 101 and the decoder 105.
  • FIG. 13 is a flowchart showing a process flow of the parity check matrix generation device according to Embodiment 4 of the present invention.
  • BRS — I, BPS — ⁇ , BPS — III, BPS — IV are used to generate, for example, a matrix of column degree 4 X row degree 4, but in embodiment 4, Generate a matrix of direct column degree 'X row degree TM' using £ IH ⁇ described later.
  • step ST41 the information length 'coding rate' sequence maximum order input unit 31 inputs a predetermined information length, code rate, and sequence maximum order.
  • step ST42 the pseudo-random number substitution matrix parameter calculation unit 42 uses a predetermined information length, code rate, and maximum order of the column input in step ST41, to describe a third pseudo-random number substitution matrix described later. Parameters for the second pseudorandom permutation matrix same as in the second embodiment are calculated.
  • step ST43 the pseudo-random number substitution matrix generation unit 43 performs pseudo-randomization based on the parameters for the third pseudo-random number substitution matrix and the parameters for the second pseudo-random number substitution matrix calculated in step ST42. Generate a first pseudorandom permutation matrix using a random number sequence and a Latin square matrix, and further replace and expand the rows and columns of the generated first pseudorandom permutation matrix to generate a third pseudorandom permutation matrix At the same time, the generated first pseudorandom substitution matrix is further compressed in the row direction to generate a second pseudorandom substitution matrix.
  • the order distribution optimization parameter calculation unit 44 uses the predetermined information length, coding rate, and maximum column order input in step ST41 to optimize the order distribution for parity check matrix generation. The possible order distribution combination parameters are calculated.
  • step ST45 the degree distribution optimizing unit 45 generates a parity check matrix by a density evolution method or a multi-edge type density evolution method with the combination of possible order distributions calculated in step ST44 as a constraint. Optimize the order distribution.
  • Step ST46 ⁇ ⁇ The parity check matrix generation unit 46 arranges the third and second pseudorandom permutation matrices generated in step ST43 according to the order distribution of the parity check matrix generation optimized in step ST45, and uses this for the LDPC code. Generate a parity check matrix.
  • the parameters for the third pseudorandom substitution matrix and the parameters for the second pseudorandom substitution matrix are calculated as follows.
  • Equation 45 After performing the same calculation as step ST 3 2 in Embodiment 3, replace the pseudo random number permutation matrix BPS — BPS —Il'BPS —m, BPS _w ⁇ ⁇ with F step S T4 3 ⁇ 1 to ST 4 Replace with a pseudorandom substitution matrix EL JP.
  • the pseudorandom number permutation matrix generation unit 43 determines the pseudorandom number sequence and the pseudorandom number sequence based on the parameters for the third pseudorandom number permutation matrix calculated in step ST42 and the parameters for the second pseudorandom number permutation matrix.
  • Equation 46 That is, here, the third pseudo-random number permutation matrix £ LH is created from the values obtained in step ST42 through the following steps ST 43-1 to ST43-9.
  • step ST43 the method of generating the second pseudo random number permutation matrix by further compressing the first pseudo random number permutation matrix in the row direction is the same as that in step ST33-1 of the third embodiment. It is the same as the process up to —8.
  • Equation 47 Generate a basic random number sequence c, defined by the following equation. In this index, ,, indicates the value of C, composer, (1), that is, the initial value.
  • Equation 50 The basic permutation matrix pattern of P'xp 'was generated up to step ST43-3. The procedure for multiplying this by b 'for rows and columns is described below.
  • the basic random number sequence (is generated in the same procedure as step ST431 1. This , is a basic random number sequence generated by changing P 'to b' for 0,.
  • This operation generates a basic permutation matrix to generate a matrix of column order ⁇ 'and row order / ⁇ ' of P '' b'xp '-b'.
  • the pseudo permutation matrix with column degree cvv 'and row degree ⁇ ' can be generated by the above operations. Since column order cvv 'and row order rw' of this pseudorandom permutation matrix are both 4 or less and there are many cases where cycle 4 does not exist, column order ⁇ 'and row order ⁇ ' have a maximum value of 4.
  • step ST44 the order distribution optimization parameter calculation unit 44 uses the predetermined information length, code rate, and maximum order of the sequence input in step ST41, as shown below. Calculate the possible combination of order distributions to optimize the order distribution of the check matrix generation.
  • Equation 56 The combination of ⁇ ⁇ ⁇ ⁇ .), C ,, c / and the triangular matrix of the upper right triangle 0 is used to find a feasible combination.
  • step ST45 the degree distribution optimizing unit 45 uses the combination of possible degree distributions calculated in step ST44 as a constraint condition, as shown below, by the density evolution method or the multi-edge type density evolution method. Optimize the order distribution of the check matrix generation. That is, here, the same as step ST15 in the first embodiment, under the constraint conditions of possible order distribution combinations for optimization of the order distribution of the parity check matrix generation calculated in step ST44 above.
  • the best degree distribution may be obtained by optimization of the degree distribution, or the multi-edge type density development method described below may be used.
  • the density evolution method described in Embodiment 1 described above is based on the fact that the Irregular LDPC code ensemble (showing the degree distribution) is expressed as a row degree ratio and a column order ratio.
  • the type ensemble is expressed as a ratio of “combinations of branches (branch types)” (order type).
  • branch types branch types
  • order type order type
  • Fig. 14 is a diagram showing an example of a multi-edge type ensemble.
  • AWGN Gaussian noise
  • Figure 4 is represented as the table shown in Figure 15.
  • FIG. 15 shows an example of the degree distribution of irregular LDPC codes.
  • the multi-edge type density evolution method will be described.
  • the probability density function is divided for each branch type without compiling messages to the node node check node together ("Mult Edge Type LDPC Codes DRAFT ", T. Richardson, R. Urbanke, on the web).
  • Equation 58 Let I be the set of all Imae branch types. The order of branch type i for ie I, the probability density function transmitted to the check node through branch type i, A. , and the probability density function transmitted from the check node through branch type i to / ⁇ To do. Calculate the probability density function for each branch type by calculating the probability density function for each order type, multiplying the ratio of the branch to the branch type i, and adding them.
  • VIER-(A 2) for all branch types /
  • the order for the channel type bl b2 is ⁇ ⁇ d h2 ,... and the initial values are / admir,, ti , Pm .... Calculate the probability density function for each order type, multiply by the ratio of the branch to the branch type i, By combining, the probability density function for each branch type is calculated.
  • the ensemble (order distribution) of the “Irregular—LDPC code” based on the required coding rate is obtained.
  • the ratio of the channel type b and degree type d to the nobleable node code length N is V
  • the ratio of the Eck node to the code length N is treated as a variable at the same time
  • step ST46 the parity check matrix generation unit 46 uses the order distribution of the parity check matrix generation optimized in step ST45 to generate the third and second pseudo-signals generated in step ST43 as shown below.
  • a check matrix for LDPC codes is generated by arranging a random number substitution matrix.
  • Equation 63 using the combination of the optimal degree distribution obtained in step ST 4 3, the third pseudorandom permutation matrix P (cw ., And the second pseudorandom permutation matrix c ,, Generate check matrix H. Regarding the arrangement of each pseudorandom substitution matrix, check the number of cycles by searching, etc., and obtain a good arrangement.
  • FIG. 17 shows an example of the degree distribution of the irregular LDPC code used in the evaluation.
  • FIG. 18 is a diagram showing the frame error rate characteristics of the irregular LDPC code.
  • Embodiment 4 Using LDPC code (ldpcl440 in the figure) with 0 bit and coding rate of 0.5, and Embodiment 4
  • FER frame error rate of LDPC code (ldpc4896 in the figure) with an information length of 4896 bits and a coding rate of 0.5, configured using p 5-24 -53, -3, «, 9, /, 2 Simulation results are shown.
  • turbo code length and coding rate of the turbo code used in the third generation communication are shown in FIG. 18 for comparison (turbo576, turbo 1440). , tu rbo3168, turbo4896).
  • the power that turbo code is dominant in the conventional short V and code length of 10000 bits or less as shown in Fig. 18, the present invention shows almost equivalent performance at information length of 576 bits.
  • the performance of the proposed LDPC code is always superior in the long! And information length, and it can be confirmed that the code construction method shows extremely high performance.
  • turbo codes With conventional LDPC codes, a code length of 10,000 bits or more shows good performance. Conversely, turbo codes are said to be superior for code lengths of 10,000 bits or less. However, in Fig. 18, the information length is 4896 bits (the code length is 0.5 so the code length is 9792 bits), which shows better performance than the turbo code. It also shows a tendency for performance to be better than turbo codes as the code length increases.
  • the conventional LDPC code has a code length of about 10,000 bits and has the same performance as the turbo code. In this invention, the code length is about 1000 bits (in Fig. 18, the information length is 576 bits and the coding rate is 0.5, so the code is It is 1152 bits long) and shows the same performance as the turbo code.
  • the same effect as in the first embodiment can be obtained, and the error correction capability can be further improved by using the multi-edge type density evolution method.
  • FIG. 19 is a block diagram showing a configuration of a parity check matrix generation apparatus for LDPC codes according to Embodiment 5 of the present invention.
  • the parity check matrix generation apparatus according to the fifth embodiment includes a pseudo random number permutation matrix generation unit 13, a check matrix generation unit 16, a pseudo random number substitution matrix parameter storage unit 17, and an order distribution storage unit 18.
  • This check matrix generation apparatus is mounted in the transmission side and reception side communication apparatuses in FIG. 4, and the generated check matrix may be output to the encoder 101 and the decoder 105, or the transmission side and the reception side.
  • the generated check matrix may be stored outside the communication device on the transmission side and the reception side and output to the encoder 101 and the decoder 105 after being installed outside the communication device on the side.
  • the pseudo random number substitution matrix parameter storage unit 17 stores parameters for the pseudo random number substitution matrix calculated using a predetermined information length, coding rate, and maximum degree of the column. is doing.
  • the parameters for this pseudorandom permutation matrix are generated externally in the same manner as in Embodiment 1.
  • the pseudo random number substitution matrix generation unit 13 uses the pseudo random number substitution matrix and the Latin square matrix to determine the pseudo random number substitution matrix according to the parameters for the pseudo random number substitution matrix stored in the pseudo random number substitution matrix parameter storage unit 17. Is generated.
  • the order distribution storage unit 18 also optimizes the order distribution of the parity check matrix that can be configured using a pseudo-random number permutation matrix calculated using a predetermined information length, coding rate, and maximum order of the column.
  • the order distribution of the check matrix generation optimized by the density evolution method is stored as a constraint with the combination of possible order distributions for optimization.
  • the order distribution of this check matrix generation is optimized externally in the same manner as in the first embodiment.
  • the check matrix generation unit 16 arranges the generated pseudo-random number permutation matrix according to the optimized check matrix generation order distribution stored in the degree distribution storage unit 18, and checks the low-density parity check code. Generate a matrix.
  • the pseudo random number substitution matrix generation unit 13 the check matrix generation unit 16, the pseudo random number substitution matrix parameter storage unit 17, and the degree distribution storage unit 18 of the first embodiment are configured.
  • the pseudo random number substitution matrix generation unit 23 and the check matrix generation unit 26 the pseudo random number substitution matrix parameter storage unit 27, and the degree distribution storage unit of the second embodiment.
  • 28 which may be configured by pseudo random number substitution matrix generation unit 33 and check matrix generation unit 36, pseudo random number substitution matrix parameter storage unit 37 and degree distribution storage unit 38 (not shown).
  • the pseudo random number permutation matrix generation unit 43, the check matrix generation unit 46, the pseudo random number permutation matrix parameter storage unit 47, and the degree distribution storage unit 48 (not shown) of the fourth embodiment may be configured. You may do it.
  • FIG. 20 is a block diagram showing a configuration of a communication apparatus according to Embodiment 6 of the present invention.
  • Transmission-side communication apparatus 201 includes check matrix generation apparatus 100 according to any one of Embodiments 1 to 5, encoder 101, and modulator 102.
  • receiving-side communication apparatus 202 includes check matrix generation apparatus 100 according to any of Embodiments 1 to 5, demodulator 104, and decoder 105.
  • parity check matrix generation apparatus 100 generates parity check matrix H by any of the methods in Embodiment 1 to Embodiment 5, and encoder 101 generates Find the generator matrix G using the check matrix H, and receive the message u---u)
  • Modulator 102 performs digital modulation such as BPSK, QPSK, and multilevel QAM on the generated codeword C and transmits the result.
  • parity check matrix generation apparatus 100 In communication apparatus 202 on the receiving side, parity check matrix generation apparatus 100 generates parity check matrix H by any one of Embodiment 1 to Embodiment 5, and demodulator 104 transmits communication path 103.
  • BPSK, QPSK, multi-level QAM, etc. are digitally demodulated to the modulated signal received via the decoder 105, and the decoder 105 uses the generated check matrix to generate the LDPC encoded demodulation result. Iterative decoding is performed using the “sum—product algorithm”, and the estimation result (corresponding to the original u U... !!) is output.
  • the use of the check matrix generation device 100 according to any one of the first to fifth embodiments makes it deterministic and stable in characteristics.
  • a parity check matrix for an Irregula LDPC code that has good error correction capability and supports various ensembles, code lengths, and coding rates.
  • the parity check matrix generation apparatus and the parity check matrix generation method according to the present invention are, for example, deterministic and stable in characteristics, have good error correction capability, and have a variety of ensembles, code lengths, and code signs. It is suitable for easily generating a parity check matrix for LDPC codes corresponding to the rate in a short time.

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)

Abstract

　所定の情報長、符号化率及び列の最大次数を用いて擬似乱数置換行列のためのパラメータを算出し、算出された擬似乱数置換行列のためのパラメータにより擬似乱数系列とラテン方陣行列を使用して擬似乱数置換行列を生成し、所定の情報長、符号化率及び列の最大次数を用いて、擬似乱数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のための取りうる次数分布の組み合わせのパラメータを算出し、算出された取りうる次数分布の組み合わせを拘束条件として検査行列生成の次数分布を最適化し、最適化された検査行列生成の次数分布により、生成された擬似乱数置換行列を配置して低密度パリティ検査符号用の検査行列を生成する。

Description

検査行列生成装置及び通信装置
技術分野
[0001] この発明は、誤り訂正符号として低密度パリティ検査 (LDPC： Low-Density
Parity-Check)符号を採用した場合における検査行列生成装置及びこの検査行列 生成装置を搭載した通信装置に関するものであり、特に、確定的で特性が安定した LDPC符号用の検査行列を探索可能な検査行列生成装置及びこの検査行列生成 装置を搭載した通信装置に関するものである。
背景技術
[0002] 以下、従来の LDPC符号用の検査行列生成方法について説明する。従来の LDP C符号ィ匕 Z復号システムにお 、て、送信側の通信装置は符号化器と変調器とを備え る構成とし、一方、受信側の装置は復調器と復号器とを備える構成とする。ここでは、 従来の LDPC符号用の検査行列生成方法を説明する前に、 LDPC符号を使用した 場合の符号化及び復号の流れにっ 、て説明する。
[0003] まず、送信側の符号化器では、後述する従来の方法で検査行列 Hを生成する。そ して、以下の条件に基づいて生成行列 Gを求める。
G :KX N行歹 lj (Kは†青報長、 Nは符号長）
GHT=0 (Tは転置行列）
[0004] その後、符号化器では、情報長 Kのメッセージ (m m〜m )を受け取り、上記生成
1 2 K
行列 Gを用いて符号語 Cを生成する。
C= (m m "*m ) G
1 2 K
= (c c - --c ) (ただし、 H (c c - --c )T=0)
1 2 N 1 2 N
[0005] そして、変調器では、生成した符号語 Cに対して、 BPSK (Binary Phase Shift
Keying)、 QPSK (Quadrature Phase Shift Keying)、多値 QAM (Quadrature Amplitude Modulation)等のデジタノレ変調を行って送信する。
[0006] 一方、受信側では、復調器が、通信路を介して受け取った変調信号に対して、 BP SK、 QPSK,多値 QAM等のデジタル復調を行い、さらに、復号器が、 LDPC符号 化された復調結果に対して、「sum— productアルゴリズム」による繰り返し復号を実 施し、推定結果 (もとの m m 〜m に対応）を出力する。
1 2 K
[0007] 次に従来の LDPC符号用の検査行列生成方法を具体的に説明する。
第 1図は例えば非特許文献 1により提案された従来の LDPC符号用の検査行列を 示す図である。第 1図に示す行列は、「1」と「0」の 2値の行列で、「1」の部分を塗りつ ぶしている。他の空欄の部分は全て「0」である。この行列は、 1行の「1」の数（これを 行の次数と表現する）が 4で、 1列の「1」の数 (これを列の次数と表現する）が 3であり 、全ての列と行の次数が均一なため、これを一般に「Regular— LDPC符号」と呼んで いる。また、非特許文献 1による符号では、例えば、第 1図に示すように、行列を 3プロ ックに分け、 2ブロック目と 3ブロック目に対してランダム置換を行っている。
[0008] し力しながら、このランダム置換には、所定のルールがないため、より特性の良好な 符号を見つけるためには、計算機による時間の力かる探索を行わなければならなか つた o
[0009] そこで、例えば、計算機探索によらなくても確定的に行列を生成でき、比較的安定 した良好な特性を示す LDPC符号として、ユークリット幾何符号を用いる方法が非特 許文献 2によって提案された。この方法では、規則的な ensemble (アンサンブル)で 構成された「Regular— LDPC符号」につ!/ヽて説明されて!ヽる。
[0010] 非特許文献 2によれば、有限幾何符号の一種であるユークリット幾何符号 EG (2, 2 6)を用いて LDPC符号の検査行列を生成する方法が提案されており、誤り率 10— 4点 において、シャノン限界から 1. 45dBに接近した特性を得ている。
第 2図は例えばユークリット幾何符号 EG (2, 22)の構成を示す図であり、行、列の それぞれの次数が 4, 4の「Regular— LDPC符号」構造をして!/、る。
[0011] したがって、ユークリット幾何符号 EG (m, 2s)の場合、その特性は以下のように規 定される。ここで、 EG (m, 2s)は GF (2s)上の m次元ユークリット幾何符号である。 符号長： N = 22s-1
冗長ビット長： N— K= 3S— 1
情報長： K= 22s— 3s
最小距離： d = 2S+ 1 密度： r= 2sZ(22s—l)
[0012] 第 2図を見ても分力るように、ユークリット幾何符号は、各行の「1」の配置が行毎に 巡回シフトした構造になっており、符号が容易にかつ確定的に構成できる特長がある
[0013] 非特許文献 2による検査行列生成方法では、さらに、上記ユークリット幾何符号に 基づいて行と列の次数を変更し、行、列を必要に応じて拡張している。例えば、 EG ( 2, 22)の列の次数を 1Z2に分離する場合、非特許文献 2では、 1列内に 4つある次 数を 1つ置きに 2個づっ分離する。
第 3図はユークリット幾何符号の列の次数を 4力も 2に規則的に分離した例を示す 図である。
[0014] 一方、上記「Regular - LDPC符号」の特性よりも「Irregular~LDPC符号」の特性 の方が良好であることが、非特許文献 3により報告された。そして、それは非特許文 献 4及び非特許文献 5によって理論的に解析された。なお、上記「Irregular~LDPC 符号」は、列と行の次数がそれぞれ又はどちらか一方が均一でない LDPC符号を表 す。
[0015] 特に、非特許文献 5では、繰り返し復号器における入力と出力の対数尤度比 (LLR )がガウス分布に近似できると仮定して LDPC符号の「Sum— Productアルゴリズム」 を解析し、良好な行と列の次数のアンサンブルを求めて 、る。
[0016] 非特許文献 1 : R.G. Gallager. "Low-Density Parity Check Codes. "IRE Trans. On IT pp21-28 January 1962.
非特許文献 2 : Y. Kou, S. Lin, and M. P. C. Fossorier, "Low Density Parity Check
Codes Based on Finite Geometries: A Rediscovery and New Results, IEEE Trans.
Inform. Theory, vol.47, No.2, pp.2711- 2736, Feb. 2001.
非特許文献 3 : M. G. Luby, M. Mitzenmacher, M. A. Shokrollahi, and D. A.
Spielman, "Improved Low-Density Parity-Check Codes Using Irregular Graphs and
Belief Propagation," Proceedings of 1998 IEEE International Symposium on
Information Theory, pp. 171, Cambridge, Mass., August 16—21, 1998.
非特許文献 4 : T. J. Richardson and R. Urbanke, i'he capacity of low-density parity-check codes under message-passing decoding," IEEE Trans. Inform. Theory, vol.47,No.2, pp.599— 618, Feb. 2001.
非特許文献 5 : S.- Y. Chung, T. J. Richardson, and R. Urbanke, "Analysis of
Sum-Product Decoding of Low-Density Parity-Check Codes Using a Gaussian Approximation," IEEE Trans. Inform. Theory, vol.47, No.2, pp.657- 670, Feb. 2001.
[0017] 従来の LDPC符号用の検査行列生成方法は以上のように行われているので、検査 行列のサイクルの最小値（内径と呼ぶ）が多くの場合 6程度であり（第 2図、第 3図の 場合は 6)、検査行列のサイクルが大きいほうが符号語のハミング距離が広がり、 Su m— product復号法における特性が良いとされている一般的な解析結果からみると、 誤り訂正能力の性能面力 は不十分であるという課題があった。
[0018] この発明は上記のような課題を解決するためになされたもので、確定的で特性が安 定し、誤り訂正能力が良好で、かつ多様なアンサンブル、符号長、符号化率に対応 した Irregular— LDPC符号用の検査行列を短時間で容易に生成できる LDPC符号 用の検査行列生成装置及びこの検査行列生成装置を搭載した通信装置を得ること を目的とする。
ここで、「確定的とは」送受信共に共通のアルゴリズムを有し、限られた諸元で同一 の結果を導き出すもので、この発明の場合、符号化率や符号長等の限られた諸元で 、 LDPC符号用の検査行列を生成することが可能という意味である。
発明の開示
[0019] この発明に係る検査行列生成装置は、所定の情報長、符号化率及び列の最大次 数を入力する情報長'符号化率'列最大次数入力部と、入力された上記所定の情報 長、符号化率及び列の最大次数を用いて後述の擬似乱数置換行列のためのパラメ ータを算出する擬似乱数置換行列用パラメータ算出部と、算出された上記擬似乱数 置換行列のためのパラメータにより擬似乱数系列とラテン方陣行列を使用して上記 擬似乱数置換行列を生成する擬似乱数置換行列生成部と、入力された上記所定の 情報長、符号化率及び列の最大次数を用いて、上記擬似乱数置換行列を使用して 構成可能な検査行列生成の次数分布の最適化のための取りうる次数分布の組み合 わせのパラメータを算出する次数分布最適化用パラメータ算出部と、算出された上 記取りうる次数分布の組み合わせを拘束条件として、密度発展法により検査行列生 成の次数分布を最適化する次数分布最適化部と、最適化された上記検査行列生成 の次数分布により、生成された上記擬似乱数置換行列を配置して低密度パリティ検 查符号用の検査行列を生成する検査行列生成部とを備えたものである。
[0020] この発明により、確定的で特性が安定し、誤り訂正能力が良好で、かつ多様なアン サンブル、符号長、符号化率に対応した Irregular— LDPC符号用の検査行列を短 時間で容易に生成できると 、う効果が得られる。
図面の簡単な説明
[0021] [図 1]従来の LDPC符号用の検査行列を示す図である。
[図 2]ユークリット幾何符号 EG (2, 22)の構成を示す図である。
[図 3]ユークリット幾何符号の列の次数を 4力も 2に規則的に分離した例を示す図であ る。
[図 4]LDPC符号ィ匕 Z復号システムの構成を示すブロック図である。
[図 5]この発明の実施の形態 1による検査行列生成装置の構成を示すブロック図であ る。
[図 6]この発明の実施の形態 1による検査行列生成装置の処理の流れを示すフロー チャートである。
[図 7]この発明の実施の形態 2による検査行列生成装置の構成を示すブロック図であ る。
[図 8]この発明の実施の形態 2による検査行列生成装置の処理の流れを示すフロー チャートである。
[図 9]情報 1ビット当りの信号電力対ノイズ電力比に対するビット誤り率特性とフレーム 誤り率特性を示す図である。
[図 10]この発明の実施の形態 3による検査行列生成装置の構成を示すブロック図で ある。
[図 11]この発明の実施の形態 3による検査行列生成装置の処理の流れを示すフロー チャートである。
[図 12]この発明の実施の形態 4による検査行列生成装置の構成を示すブロック図で ある。
[図 13]この発明の実施の形態 4による検査行列生成装置の処理の流れを示すフロー チャートである。
[図 14]Multi— edge typeのアンサンブルの例を示す図である。
[図 15]irregular LDPC符号の例を示す図である。
[図 16]第 15図の次数分布に従って構成した検査行列の例を示す図である。
[図 17]評価に用いた irregular LDPC符号の次数分布の例を示す図である。
[図 18]irregular LDPC符号のフレーム誤り率特性を示す図である。
[図 19]この発明の実施の形態 5による検査行列生成装置の構成を示すブロック図で ある。
[図 20]この発明の実施の形態 6による通信装置の構成を示すブロック図である。 発明を実施するための最良の形態
[0022] 以下、この発明をより詳細に説明するために、この発明を実施するための最良の形 態について、添付の図面に従って説明する。なお、この実施の形態によりこの発明が 限定されるものではない。
実施の形態 1.
まず、この実施の形態 1の LDPC符号用の検査行列生成方法を説明する前に、生 成された LDPC符号用の検査行列を使用する符号化器及び復号器の位置付け、並 びに「Irregular— LDPC符号」用の従来の検査行列生成方法にっ 、て説明する。
[0023] 第 4図は LDPC符号ィ匕 Z復号システムの構成を示すブロック図である。第 4図にお いて、送信側の通信装置は符号化器 101と変調器 102を含む構成とし、受信側の通 信装置は復調器 104と復号器 105を含む構成とする。ここで、生成された LDPC符 号用の検査行列を使用する場合の符号ィ匕及び復号の流れについて説明する。
[0024] 送信側の通信装置では、後述するように、 LDPC符号用の検査行列生成装置によ り生成された検査行列 Hを使用する。そして、符号化器 101では、以下の条件に基 づ、て生成行列 Gを求める。
G :KX N行歹 lj (Kは†青報長、 Nは符号長）
GHT=0 (Tは転置行列） [0025] その後、符号化器 101では、情報長 Kのメッセージ (u u…！！ ）を受け取り、上記生
1 2 K
成行列 Gを用いて符号語 Cを生成する。
C= (u U…！！ )G
1 2 K
=(c c ---c ) (ただし、 H(c c ---c )T=0)
1 2 N 1 2 N
そして、変調器 102では、生成した符号語 Cに対して、 BPSK、 QPSK、多値 QAM 等のデジタル変調を行って送信する。
[0026] しかし、この符号化器 101の場合、生成行列 Gを検査行列 H力もガウス消去法等で 生成しなければならず、かつ検査行列 Hが低密度であるのに対して、生成行列 Gが 比較的高密度で 2元の生成行列中の" 1 "の比率が大きいため、生成行列 Gを格納す るメモリ量も多く必要であるという課題がある。
この課題を解決するために、検査行列 Hの中に下三角構造 (又は上三角構造)の 行列を含む構成にする手法が下記文献等で提案されている。
M.Rashidpour and S.H. Jamah, Low-Density Parity-Check Codes with Simple Irregular Semi-Random Parity Check Matrix for Finite-Length Applications," pp439-443, in proc. IEEE PIMRC 2003.
[0027] [数 1] 今 Tを以下のような ΛίχΛ ,ΜΕΖ,。の下三角構造の行列とする。 ここで Ζ>„ は正の整 数全体の集合を意味する。
Figure imgf000009_0001
行 j列の要素 で構成される 2元 Λ/ e Ζ.„ ,行列 H - [A,Jを以下のように表す。 ここで//は Λ/xWの行列とする。 また は任意の 2元 χΚ,Κ =Ν-Μ行列である。 組織符号語べクトル Cを (^(«,,^,,.^^ρ,, ,...,^)とすると、 この符号語べクトル
Cは / Cr _0を満たし、 情報メッセージ "- ("，，"2,...，"κ)が与えられた場合、 符号化器 1 01が以下の式によりパリティ要素を生成する。
Figure imgf000009_0002
この手順により、生成行列 Gを用いなくても符号ィ匕が容易に実現できるため、この発 明では、下三角構造の行列を含めた検査行列 Hを構成することを条件とする。
[0028] 一方、受信側の通信装置では、事前に送信側の通信装置から必要な諸元に関す る情報を受け取り検査行列を生成し、復調器 104が、通信路 103を介して受け取つ た変調信号に対して、 BPSK、 QPSK、多値 QAM等のデジタル復調を行い、さらに 、復号器 105が、 LDPC符号化された復調結果に対して、生成された検査行列を用 いて「sum— productアルゴリズム」による繰り返し復号を実施し、推定結果 (もとの u U…！！ に対応）を出力する。
2 K
[0029] 次に上記非特許文献 5によって理論的に解析された、「Irregular— LDPC符号」用 の従来の検査行列生成方法について詳細に説明する。ここでは、 sum— product復 号法の復号プロセスに従って、復号器 105における入力と出力の対数尤度比 (LLR )の確率密度関数を反復計算していくことにより、 sum— product復号法の復号特性 を解析する密度発展法と!ヽぅ手法及びその近似計算法であるガウス近似法を用いて 、良好な行と列の次数のアンサンブルを求めて 、る。
[0030] なお、上記非特許文献 5に記述された LDPC符号用の検査行列生成方法である密 度発展法及び復号器 105における入力と出力の対数尤度比 (LLR)がガウス分布に 近似できるとして計算を簡易化したガウス近似法（Gaussian Approximation)では、前 提として、検査行列をタナ-グラフの表現を用いており各列をバリアブルノードと定義 し、各行をチェックノードと定義する。また、ある列とある行の交点に「1」がある場合に 、ノリアブルノードとチェックノードは「枝」で接続されて!ヽると!/、う。
[0031] まず、チェックノードからバリアブルノードへの LLRメッセージ伝搬を解析する。 0< s<∞と 0≤t<∞という条件において、以下の関数の（1)式を定義する。ただし、 s = m は uOの平均値であり、 uOは分散値 σ 2のガウスノイズを含む伝送路を経由して uO n
受信した信号の対数尤度比 (LLR)のアンサンブル平均であり、 tは所定の繰り返し の時点におけるチェックノードの LLR出力値のアンサンブル平均である。
[数 2] (i - l)t)
Figure imgf000011_0001
なお、 λ p iは、それぞれ次数 iのバリアブルノードとチェックノードに属する枝の比 率を表す。また、 dは最大ノ リアブルノードの次数であり、 dは最大チェックノードの次
1 r
数である。また、上記え（X)及び p (X)は、それぞれバリアブルノード及びチェックノ ードの次数配分 (バリアブルノードとチェックノードの各 1行、各 1列内の「1」の数を次 数と表現する）の生成関数を表し、下記の（2)式及び（3)式のように表すことができる
[数 3]
Figure imgf000011_0002
Figure imgf000011_0003
また、上記（1)式の φ (X)は下記の (4)式のように定義する。
[数 4]
〔4)
Figure imgf000011_0004
ここで、 Rは実数全体の集合を表し、 uはバリアブルノードの LLR出力のアンサンブル 平均を表す。
そして、上記（1)式は、等価的に下記の（5)式と表すことができる。
[数 5]
Figure imgf000012_0001
… ） なお、 tは 1番目の繰り返し時点におけるチェックノードの LLR出力値のアンサンブ
1
ル平均である。
[0034] ここで、誤りが 0となりうる SNRの限界（threshold)を求めるための条件は、 1→∞のと きに t (s)→∞(R+と表現する）となることであり、この条件を満たすためには、以下の
1
条件である（6)式を満たす必要がある。
[数 6] t < f(s,t), 全ての tER+ …(
[0035] 次にノ リアブルノードからチェックノードへの LLRメッセージ伝搬を解析する。 0< s <∞と 0<r≤lという条件において、以下の関数の（7)式を定義する。なお、 rの初期 xは （s)である。
0
[数 7]
Figure imgf000012_0002
h(s,r) = ^ihi(S,r) - (τ )
i-2
[0036] そして、上記（7)式は、等価的に下記の（8)式と表すことができる。
[数 8] = ^-x) … ）
[0037] ここで、誤りが 0となりうる SNRの限界（threshold)を求めるための条件は、 r (s)→0
1 となることであり、この条件を満たすためには、以下の条件である（9)式を満たす必要 がある。
[数 9] r > h(s,r), 全ての Γ Ε(0，φ(3)) [0038] さらに、上記非特許文献 5では、上記式を用いて以下の手順でバリアブルノードと チェックノードの最適な次数を探索して 、る (密度発展法のガウス近似法)。
[ステップ ST101]
ノ リアブルノードの次数配分の生成関数え（X)とガウスノイズ σ が与えられていると 仮定し、チェックノードの次数配分の生成関数 ρ (X)を変数として、符号化率が最大 となる点を探索する。なお、この探索における拘束条件は、 (1) = 1と正規ィ匕するこ とと、上記 (6)式を満たすことである。
[ステップ ST102]
チェックノードの次数配分の生成関数 ρ ）とガウスノイズ σ が与えられて 、ると仮 定し (例えば、 ST101の結果より得られる値）、ノ リアブルノードの次数配分の生成関 数え（X)を変数として、符号ィ匕率が最大となる点を探索する。なお、この探索におけ る拘束条件は、 λ (1) = 1と正規ィ匕することと、上記 (9)式を満たすことである。
[ステップ ST103]
最大の符号ィ匕率を求めるために、上記ステップ ST101と上記ステップ ST102を繰 り返し実行し、ノ リアブルノードの次数配分の生成関数え（X)とチェックノードの次数 配分の生成関数 ρ (X)のより良好なアンサンブルを線形計画法で探索する。
[ステップ ST104]
ガウスノイズ σ より信号電力を 1と正規ィ匕して、誤りが 0となりうる SNRの限界（ threshold)を下記の（10)式より求める。
[数 10] threshold (dB)= -10 * loglOI2 * σ^ ) - - - (l θ)
[0039] し力しながら、上記非特許文献 5では、符号ィ匕率の最大値により得られる検査行列 が流動的になり、設計時の仕様として固定される符号ィ匕率が変動してしまうという課 題がある。また、上記非特許文献 5では、ノ リアブルノードの次数配分の生成関数え (X)の導出と、チェックノードの次数配分の生成関数 ρ (X)の導出とを所定回数にわ たって繰り返し行って 、るために、探索処理にある程度の時間を要すると 、う課題や 、多様なアンサンブル、符号長、符号ィ匕率に容易に対応することができないという課 題もある。
[0040] そこで、この実施の形態 1にお 、ては、確定的で特性が安定し、かつ多様なアンサ ンブル、符号長、符号ィ匕率に対応した「Irregular— LDPC符号」用の検査行列を、短 時間で容易に探索する方法にっ 、て説明する。
[0041] 第 5図はこの発明の実施の形態 1による LDPC符号用の検査行列生成装置の構成 を示すブロック図である。この検査行列生成装置は、情報長'符号化率'列最大次数 入力部 11、擬似乱数置換行列用パラメータ算出部 12、擬似乱数置換行列生成部 1 3、次数分布最適化用パラメータ算出部 14、次数分布最適化部 15及び検査行列生 成部 16を備えている。この検査行列生成装置は、第 4図における送信側及び受信側 の通信装置内に搭載され、生成した検査行列を符号化器 101及び復号器 105に出 力しても良いし、送信側及び受信側の通信装置外に設置され、生成した検査行列を 送信側及び受信側の通信装置内のメモリ（図示せず）に格納し符号化器 101及び復 号器 105に出力しても良い。
[0042] 第 6図はこの発明の実施の形態 1による検査行列生成装置の処理の流れを示すフ ローチャートである。ステップ ST11において、情報長'符号化率'列最大次数入力部
11は所定の情報長、符号ィ匕率及び列の最大次数を入力する。ステップ ST12にお いて、擬似乱数置換行列用パラメータ算出部 12は上記ステップ ST11で入力された 所定の情報長、符号化率及び列の最大次数を用いて擬似乱数置換行列のための ノ ラメータを算出する。
[0043] ステップ ST13において、擬似乱数置換行列生成部 13は、上記ステップ ST12で 算出された擬似乱数置換行列のためのパラメータにより擬似乱数系列とラテン方陣 行列を使用して擬似乱数置換行列を生成する。ステップ ST14において、次数分布 最適化用パラメータ算出部 14は上記ステップ ST11で入力された所定の情報長、符 号化率及び列の最大次数を用いて、擬似乱数置換行列を使用して構成可能な検査 行列生成の次数分布の最適化のための取りうる次数分布の組み合わせのパラメータ を算出する。
[0044] ステップ ST15において、次数分布最適化部 15は、上記ステップ ST14で算出され た取りうる次数分布の組み合わせを拘束条件として、密度発展法により検査行列生 成の次数分布を最適化する。ステップ ST16において、検査行列生成部 16は上記ス テツプ ST15で最適化された検査行列生成の次数分布により、上記ステップ ST13で 生成された擬似乱数置換行列を配置して LDPC符号用の検査行列を生成する。
[0045] 次に第 6図のフローチャートにおけるステップ ST11— ST16の詳細処理について 説明する。
ステップ ST11にお 、て、情報長 ·符号化率 ·列最大次数入力部 11は所定の情報 長、符号化率、列の最大次数として、システムで事前に設定されている情報長 K= l 44、符号化率 rate = 0. 5、列の最大次数 d=4を入力する。
[0046] ステップ ST12において、擬似乱数置換行列用パラメータ算出部 12は上記ステツ プ ST11で入力された情報長 K= 576、符号化率 rate = 0. 5、列の最大次数 d=4を 用いて、 p2 X p2 (p+ 1は素数)の擬似乱数置換行列 SPのためのパラメータ pを以下 のように算出する。
擬似乱数置換行列 SPの行方向の個数 =d=4
検査行列 Hの列数 N (=符号長 N) =K/rate
= 144/0. 5 = 288
検査行列 Hの行数 M = N X ( 1— rate)
= 288 X 0. 5 = 144
擬似乱数置換行列3?のサィズ 2=1^7(1= 14474 = 36
よって p = 6
[0047] ステップ ST13において、擬似乱数置換行列生成部 13は、上記ステップ ST12で 算出された擬似乱数置換行列 SPのためのパラメータ pにより、以下に示すステップ S T13— 1— ST13— 6の手順に示すように、擬似乱数系列とラテン方陣行列を使用して p2 X p (p+ 1は素数)の擬似乱数置換行列 SPを生成する。
LDPC符号を用いた符号ィ匕 Z復号においては、一般的に、タナーグラフにおける 2 部グラフ上にぉ 、て「サイクル 4」及び「サイクル 6」 t 、つた検査行列の小さ!/ヽサイク ルの個数が少ないほど良好な特性を得ることができる。したがって、 LDPC符号とし ては、「サイクル 4」や「サイクル 6」といった小さいサイクルの発生を抑制する構造が望 ましぐ擬似乱数置換行列生成部 13は検査行列の小さいサイクルの発生を抑制する W 200
擬似乱数置換行列 SPを生成する。
[0048] [ステップ ST13—1]
[数 11] 以下の式で定義される擬似乱数系列 を生成する。
C(l) = l, C(/+l)-G„xC(i)mod Ρ, ί ，2,.. P _2,
ここで Pは尸= +1となる素数であり、 G。は の原始元である。 なお、 （V(P)の要 素である の連続するべき乗の値 (例えば ， ，...， (7, ' - が GF(P)の非零の全ての要 素を示す場合、 を原始元とよぶ。 この系列 は要素間の差分が非一様であり重複 する要素がない擬似乱数系列となる。
例として尸- 7,(^ =3の場合を以下に示す。
1 2 3 4 5 6
C(i)： 1 3 2 6 4 5
[0049] [ステップ ST13—2]
[数 12]
擬似乱数系列 C )に対する下式に示す置換パターン S '）を生成する。
·('') (；><'.) mod尸， /' =1,2 1,2,..., p.
この操作により、 i^.(i)を 行 i列の行列で表現した場合、 各行内の各要素 および各列内の各要素が全て異なるラテン方陣を構成することができる。 例を以下に示す。
/: 1 2 3 4 5 6
LU,(i) 1 2 3 4 5 6
LB2(i) 2 4 6 1 3 5
Figure imgf000016_0001
LBt(i) 4 1 5 2 6 3
LB;(i) 5 3 1 6 4 2
LB, i) 6 5 4 3 2 1
[0050] [ステップ ST13—3]
[数 13] 下式により擬似乱数系列及びその置換パターンを要素にもつラテン方陣行列 ((?，り を生成する。
i , (q, LB^dq + _ 2) mod p) + 1)) = C(i),
j -1,2,..., p, i - 1,2，...， p, q- 2,.., p,
この操作により、 ( ）を使って Ll (q, /) 2 (q ),..., L„ (？,りの p個の異なったラテ ン方陣行列を生成することができる。
例を以下に示す。
/; i 2 3 4 5 6
(U) 1 3 2 6 4 5
(2J) 3 2 6 4 5 1 (6，/) 5 1 3 2 6 4
6 1 4 3 5 2
(2，/) 1 4 3 5 2 6
W) 2 6 1 4 3 5
5 4 6 2 3 1
2，り 4 6 2 3 1 5 6 ) 1 5 4 6 2 3
[0051] [数 14]
この例を見ても分かるように、 q=iに対して =1,/)の各要素を下記のように並 ベてみると、 各行内の各要素及び各列内の各要素が異なるラテン方陣行列になってい ることが分かる。
い. 1 2 3 4 5 6
1,(1 1 3 2 6 4 5
/.2 (],() 6 1 4 3 5 2
/.,(!,/) 4 2 1 5 6 3
L,(,i) 3 6 5 1 2 4
Λ5(1,ί) 2 5 3 4 1 6
Lft(l, ) 5 4 6 2 3 1
同様に q = pに対して、 (q,l)の各要素並べた j行 ί列の行列をみてもラテン方陣 になる。
[0052] [ステップ ST13— 4] 下式によりラテン方陣行列 りとパラメ一夕 pを用いて基本置換行列パターン （ を生成する。
BPj(q ) = Lj(q,i) + p x(i-\),
= 12,..., p, i - 1,2,..., p, q - 1,2,..' p.
この操作により .の列次数 1、 行次数 1の擬似乱数置換行列を生成する ための基本置換行列パターンを生成できる。
^ (<? i)は ^りのある に対して 行/列の行列を生成した場合， 各行をひと つインクリメントするごとに pだけその行内の各要素に加算したものである。 例を以下に示す,
[数 16]
i 1 2 3 4 5 6
( ) 3 6 5 1 2 4
P4(2 i) 12 11 7 10 9
17 13 14 16 15 18
Β Λ^) 19 20 22 21 24 23
BP4(5，i) 26 28 27 30 29 25
34 33 36 35 31 32
前ステップ （3) の議論よりある に対する の各要素を並べた 佇/列の 行列もラテン方陣になる。 また、 この操作により、 ある に対する flP,. ( りから 生成する 行/列の行列の同一列の要素間の距離は広がる。
[0054] [ステップ ST13—5]
[数 17] 下式により基本置換行列パターン BP, (q, ,·)の開始点を. 固シフトした複数の基本置換 系列 5尸,., (りを生成する。
P, (0 = BPj (( + Λ-2 )mod p)+i,\i/p]),
j = 1,2,..., p, i l,2,... p2, >5 12".,
この操作により P 2xP 2.の列次数 1、 行次数 1の擬似置換行列を生成するた めの基本置換系列を生成できる。
この系列の近接する要素間の距離は広がっている。 隣接する要素間の距離を 広げることで実施の形態 1で述べた Tの下三角行列との間で発生する短いサ ィクル、 つまりサイクル 4やサイクル 6の発生を抑制できる。
[0055] 例を以下に示す, i： 1 2 3 4 5 6
SJ 3 (i) 17 19 26 34 3 18
7 8 9
('■) 13 20 28
1 14 15
' " (り 14 22 27
19 2U 21
S尸 4J (り 16 21 30 3 133331223 ο 66432521
25 26 27
15 24 29 ί 31 32 33
18 23 25
[ステップ ST13—6]
[数 19] 複数の基本置換系列 SPj、 (りから下式により定義される 2元 2の複数の擬似乱数 置換行列 .， -( Jを生成する。
Figure imgf000019_0001
この操作により、 ρ 2 χ / の列次数 1、 行次数 1の複数の擬似乱数置換行列を 生成できる。 この, ,.、において の値が等しく、 ）の値が異なる複数の, の 組み合わせで構成される行列には、 各行の " 1 "のある列番号の一致するもの が 2個以上の場合がないためサイクル 4は存在せず、 最小サイクル数は 6とな る。 また、 複数の, 、において の値が異なった場合でも、 サイクル 4となら ない組み合わせが存在しており、 この組み合わせは探索等で求めることができ る。
例を以下に示す。
Figure imgf000019_0002
ステップ ST14において、次数分布最適化用パラメータ算出部 14は上記ステップ S Ti lで入力された所定の情報長 K、符号ィ匕率 rate及び列の最大次数 dを用いて、擬 似乱数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のため の取りうる次数分布の組み合わせのパラメータを算出する。ここでは、列の最大次数 dと下三角行列を含む検査行列が持つ制約条件のもとで、取りうる次数分布の組み 合わせを算出する。
[0058] 目標とする列の最大次数 dと情報長 Kと符号ィ匕率 rateから下式に示すような検査行 列 Hの制約条件を求める。
[数 20]
Figure imgf000020_0001
ここで、 5尸は任意の《 を意味する。 また、 対角線上に配置している \の連なり は Hの最 位行から最下位行までを通して τのように階段状に 1を配していることを 意味する。
[0059] ここで定まる制約条件は擬似乱数置換行列 SPのサイズとその最大配置数と配置箇 所である。
例として、情報長 K= 144、符号化率 mte = 0. 5、列の最大次数 d=4の場合を以 下に示す。
擬似乱数置換行列 SPの行方向の個数 =d=4
検査行列 Hの列数 N(=符号長 N) =K/rate
= 144/0. 5 = 288 検査行列 Hの行数 M=NX (1 - ra te)
= 288X0. 5 = 144
擬似乱数置換行列のサイズ p2 = MZd= 144/4 = 36
よって p = 6 情報長部分の擬似乱数置換行列 SPの最大数
=NXrate/p2= 288X0. 5/36=4
これらの結果から、検査行列 Hに対し、左端力 4X4個の擬似乱数置換行列 SPと 検査行列 Hの右側の左下三角部分に 3 + 2+1個の擬似乱数置換行列 SPが配置可 能となる。
[0060] この最大次数と下三角行列を含む検査行列が持つ制約条件の下で、取りうる次数 分布の組み合わせを算出する。次数分布とは、例えば、列次数 4の列数の検査行列 全体の列数に対する比率を示しており、その比率に応じて列次数 4となるように擬似 乱数置換行列（列重み 1、行重み 1)を列方向に 4個組み合わせて作成する。また、 行次数に対しても同様に行う。例えば、上記例では列重み 4が検査行列 H全体の列 数に対して擬似乱数置換行列を使って取りうる比率は 0— 4Z8であり、列重み 3の取 りうる比率は 0— 5Z8となる。これらの列次数及び行次数と取りうる比率の組み合わ せを、次のステップ ST15で行う次数分布の最適化のためのパラメータとして出力す る。
[0061] ステップ ST15において、次数分布最適化部 15は、上記ステップ ST14で算出され た取りうる次数分布の組み合わせを拘束条件として、以下に示すように密度発展法 により検査行列生成の次数分布を最適化する。ここでは、上記ステップ ST14で求め た次数分布の組み合わせのうち最適な次数分布を算出する。ここで、次数分布とは、 例えば、列次数 4の列数が全体の検査行列の何パーセントかを示しており、その比 率に応じて列次数 4となるように擬似乱数置換行列（列重み 1、行重み 1)を列方向に 4個組み合わせて作成する。これにより次のステップ ST16にお 、て検査行列を生成 できる。
[0062] このステップ ST15において、上記ステップ ST14で求めた検査行列 Hにおける擬 似乱数置換行列 SPのサイズとその最大配置数と配置箇所の制約条件のもとに、検 查行列 H内の擬似乱数置換行列 SPの箇所において任意の擬似乱数置換行列 SP を配置する、又はその箇所は p2 X p2の 0行列を配置するという条件で組み合わせを 選択できるようにする。
[0063] つぎに、密度発展法のガウス近似法における最適化を用いて、要求された符号ィ匕 率に基づく「Irregular— LDPC符号」のアンサンブル（次数配分）を求める。
列の次数配分の生成関数え（X)と行の次数配分の生成関数 p (X)を同時に変数と して扱い、ガウスノイズ σ が最大になるように、線形計画法で最適な生成関数え（X) と生成関数 Ρ )を探索する（下記の（11)式を参照)。この探索における拘束条件は 、後述する（12)式を満たすことである。
このとき、列の次数配分の生成関数え（X)と行の次数配分の生成関数 p (X)の変 数のとりうる値は、上記検査行列 Hの制約条件で算出される値のみとする。
[数 21] (5，Γ) = φ ( 8 + (ί— ι) 2ドρ j)φ(ι _ (ι—り h(s，r)
Figure imgf000022_0001
λ(ΐ) = 1,ρ(ΐ) = 1 -"(1 2)
r > h(s,r),全ての Γ Ε(0，φ(5))
[0065] なお、上記 sは、送信信号として { 1, 1 }の 2値信号が出力され、ガウス通信路を通 つて受信した信号の対数尤度比 (LLR)の平均値であり、 s = 2/ σ 2により算出でき る。
[0066] このように、所定の条件を満たす生成関数え（X)と生成関数 ρ (X)を 1回の線形計 画法で求めているため、非特許文献 5のように、生成関数え（X)の導出と生成関数 ρ (X)の算出を繰り返し実行し、双方の最適値を求める方法よりも、容易かつ短時間に 、確定的でかつ特性が安定したアンサンブルを生成できる。
[0067] ステップ ST16において、検査行列生成部 16は、上記ステップ ST15で最適化され た検査行列生成の次数分布により、上記ステップ ST13で生成された擬似乱数置換 行列 SPを下記に示すように配置して LDPC符号用の検査行列を生成する。すなわ ち、ここでは、上記ステップ ST15で求めた次数分布に基づき、検査行列 Ηにおける 擬似乱数置換行列 SPの配置を決定する。
[数 22]
Figure imgf000023_0001
の配置およびパラメータを決める際には、 サイクル 4 , 6等の短いサイクル の発生がない組み合わせを ,のゾ及び. Sを調整して探索することで決定する。 SP,―、は
/及び .、を変更することで異なった擬似乱数置換行列となり、 その組み合わせによるサ ィクル 4や 6発生はもともと起こりにくい構造となっているため、 探索より組み合わ せの導 ! ,は比較的容易である点が特徴として上げられる。
[0068] 以上のように、この実施の形態 1によれば、確定的で特性が安定し、誤り訂正能力 が良好で、かつ多様なアンサンブル、符号長、符号ィ匕率に対応した Irregular~LDP C符号用の検査行列を短時間で容易に生成できると 、う効果が得られる。
[0069] 実施の形態 2.
次に実施の形態 1を拡張し、複数の符号ィ匕率において互換性のある (情報長が変 わらず符号化率のみ変化する） RC(Rate comapat¾le)— LDPC符号の検査行列生成 方法について説明する。
[0070] 第 7図はこの発明の実施の形態 2による LDPC符号用の検査行列生成装置の構成 を示すブロック図である。この検査行列生成装置は、情報長'符号化率'列最大次数 入力部 21、擬似乱数置換行列用パラメータ算出部 22、擬似乱数置換行列生成部 2 3、次数分布最適化用パラメータ算出部 24、次数分布最適化部 25及び検査行列生 成部 26を備えている。この検査行列生成装置は、第 4図における送信側及び受信側 の通信装置内に搭載され、生成した検査行列を符号化器 101及び復号器 105に出 力しても良いし、送信側及び受信側の通信装置外に設置され、生成した検査行列を 送信側及び受信側の通信装置内のメモリ（図示せず）に格納し符号化器 101及び復 号器 105に出力しても良い。
[0071] 第 8図はこの発明の実施の形態 2による検査行列生成装置の処理の流れを示すフ ローチャートである ステップ ST21において、情報長'符号化率'列最大次数入力部 21は所定の情報長、複数の符号ィ匕率及び無符号ィ匕を除く最大符号ィ匕率の列の最 大次数を入力する。ステップ ST22において、擬似乱数置換行列用パラメータ算出 部 22は上記ステップ ST21で入力された所定の情報長、複数の符号化率及び無符 号ィ匕を除く最大符号化率の列の最大次数を用いて擬似乱数置換行列のためのパラ メータを算出する。
[0072] ステップ ST23において、擬似乱数置換行列生成部 23は、上記ステップ ST22で 算出された擬似乱数置換行列のためのパラメータにより擬似乱数系列とラテン方陣 行列を使用して擬似乱数置換行列を生成する。ステップ ST24において、次数分布 最適化用パラメータ算出部 24は上記ステップ ST21で入力された所定の情報長、複 数の符号化率及び無符号化を除く最大符号化率の列の最大次数を用いて、擬似乱 数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のための取り うる次数分布の組み合わせのパラメータを算出する。
[0073] ステップ ST25において、次数分布最適化部 25は、上記ステップ ST24で算出され た取りうる次数分布の組み合わせを拘束条件として、密度発展法により複数の符号 化率に対する検査行列生成の次数分布を最適化する。ステップ ST26において、検 查行列生成部 26は上記ステップ ST25で最適化された複数の符号化率に対する検 查行列生成の次数分布により、上記ステップ ST23で生成された擬似乱数置換行列 を配置して LDPC符号用の検査行列を生成する。
[0074] 次に第 8図のフローチャートにおけるステップ ST21— ST26の詳細処理について 説明する。
ステップ ST21において、情報長'符号化率'列最大次数入力部 21は所定の情報 長 K、複数の符号化率 R (1)及び無符号化を除く最大符号化率の列の最大次数 dを 入力する。ここで、複数の符号ィ匕率 R (l)において、 1= 1, 2, . . . , 0<R(1) <R (2) く. . .く 1とする。このステップ ST21では、実施の形態 1のステップ ST11と比較し て、複数の符号化率 R (1)及び無符号化を除く最大符号化率の列の最大次数 dを入 力する点が異なっている。
入力された符号化率 R (1)の中で、最も大き 、符号化率 R (max— 1)のための検査 行列 H の行数を bとし、列数 =符号長を Nとする。ここで、 R(max) = 1は無符号 化を意味しており、この説明では用いないものとする。
[0075] ステップ ST22において、擬似乱数置換行列用パラメータ算出部 22は、上記ステツ プ ST21で入力された所定の情報長 K、複数の符号化率 R(l)及び無符号化を除く 最大符号ィ匕率の列の最大次数 dを用いて、以下に示すように、擬似乱数置換行列の ためのパラメータを算出する。
検査行列 Hの列数 N(=符号長 N) =K/R(1)
検査行列 Hの行数 M = N-K
次に符号ィ匕率 R (max— 1)の列の最大次数 d (システムで事前に設定されている値） 力 擬似乱数置換行列のサイズ P2を以下より求める。
p2 =M/d
例えば、 K=864、 R(max-l) =R(3) =2/3, d=3とすると、
N=K/R(3) =864/(2/3) =1296
M=N—K= 1296— 864=432
p2 =M/d=432/3 = 144
よって、 p = 12となる。
[0076] ステップ ST23にお 、て、擬似乱数置換行列生成部 23は、実施の形態 1の擬似乱 数置換行列生成部 13が行うステップ ST13の処理と同じ処理を行う。すなわち、擬似 乱数置換行列生成部 23はステップ ST13—1—ステップ ST13— 6の処理を行う。
[0077] ステップ ST24において、次数分布最適化用パラメータ算出部 24は上記ステップ S T21で入力された所定の情報長 K、複数の符号ィ匕率 R(l)及び無符号ィ匕を除く最大 符号ィ匕率の列の最大次数 dを用いて、以下に示すように、擬似乱数置換行列を使用 して構成可能な検査行列生成の次数分布の最適化のための取りうる次数分布の組 み合わせのパラメータを算出する。
[0078] 例えば、上記ステップ ST21で、情報長 K=864、複数の符号化率 R(l) = 1/3, R(2) =1Z2、 R(3) =2/3,無符号ィ匕を除く最大符号ィ匕率 R (3)の列の最大次数 d =3を入力したとする。ここで、複数の符号化率 R(l)のうち、最も高い符号ィ匕率 R(3) = 2Z3に対応した検査行列の行数と列数を Mと Nとする。また、ステップ ST22で
3 3
、M =M、N =Nとし、 p=12が導出されている。符号化率 R(2)=lZ2からその 検査行列の行数と列数を Mと Nとし、情報長 Kは変わらないとすると、
2 2
列数 N =K/R(2) =864/(1/2) =1728,
2
行数 M =N X (1-R(2))=1728X (1/2) =864
2 2
となる。同様に、符号ィ匕率 R(1)=1Z3からその検査行列の行数と列数を Mと Nと
1 1 し、情報長 κは変わらないとすると、
列数 N =K/R(1) =864/(1/3) =2592,
1
行数 M =N X (l-R(l)) =2592X (2/3)
1 1
= 1728
となる。行数 M =1728に対し、 p2 Xp2 =144 X 144の擬似乱数置換行列 SPは行
1
方向に 12個配置可能であり、 K=864より、列方向の情報長部分には、擬似乱数置 換行列 SPは 12個配置可能となる。これらの結果から、 R(l)の符号化率を実行する ための検査行列 Hに対し左端から 12X6個の p2 Xp2 =144 X 144の擬似乱数置換 行列 SPと、検査行列 Hの右側の左下三角部分に 11 + 10 + · · · +1個の擬似乱数 置換行列 SPが配置可能となる。この検査行列は次のステップ ST25にて説明する方 法で R(2)、 R (3)も構成することができる。
[0079] ステップ ST25において、次数分布最適化部 25は、上記ステップ ST24で算出され た取りうる次数分布の組み合わせを拘束条件として、下記に示すように密度発展法 により複数の符号化率に対する検査行列生成の次数分布を最適化する。
[0080] [数 23]
• F2= {0, 1} を 2 元の有限体とする。 符号語 CCF とし、 H- ( ）を Cに対応する M Wのパリティ検査行列とする。 符号語 Cはパリティ検査行列 Hx = 0のパリティ 検査式により定義される Λ/個の 2元パリティチェックビットを満足する長さ Wビ ットの系列 Xからなる。 行列 Hがフルランクであると仮定すると、 情報ビット数 は/ f = JV- Mとなり、 符号化率は となる。
[0081] [数 24] さらに、 ΗΛ|,>を符号化率を可変である LDPC 符号 (： w,のパリティ検査行列とする。 ここで、 ；?(/)は符号化率であり、 Λ ),/一 l,2,...，nm 0 < (1) < /f(2) <… < (max) = Iである (max) =1は無符号化を意味する） 。 また、 を _υのための AixWのパリテ ィ検查行列とする。 そして H _2)を Ηλ— と Μ χ ί の 0行列と o の追加の パリティ検査行列 からなる以下の式で定義される (M )~ X (Nil)のパリティ検 査行列とする。
u
H " /i
2)はともにフルラン^)
N-M
J? Cm ax - 1)- つ、 N- M
N '
[0082] これを図で表現すると以下のようになる。
[数 25]
Figure imgf000027_0001
[0083] これを一般化すると以下のようになる c
[数 26]
H
ここ
Figure imgf000027_0002
の差を示す。
の設計のために密度発展法のガウス近似法により全ての H„,,,,Z-は..., maX-lの次 数配分を同時に最適化する。 最適化問題を以下に示す。
[0084] [数 27] ここでは、 尸 w)を最小化するような ΗΛ(,Ίの次数配分を求める。 ここで、 ガウス近似法で推定する H の反復閾値の SNRとシャノン限界との 差を (iBで表現したものであり、 以下の式で表現する。
讀 E
SNRm 1
一
2び"2
= SNR〃、 ½ (り で眉 "' 、 、 及び βはそれぞれ符号化率 ? (りの際の SNRのシャノン限 界、 ガウス通信路のノイズの分散値、 およびガウス通信路を通って受信した信号の対 数尤度比 （LLR) の平均値及び信号のエネルギーである。
このそれぞれの符号化率に対し （18) 式の条件の下、 G^>w)及び （14) 式を用 い、 λ, ,) (；),p argi ―) 1.. max
Figure imgf000028_0001
を求める。
[0085] 上記ステップ ST24で用いた R(l)=l— M ZN =1/3 R(2)=1—M /N =：
1 1 2 2
/2,R(3)=1-M
3 ZN 3 =2Z3に対し、説明に使用した変数との対応を示すと以 下のようになる。
M =M N =N M =M +t N =N +t
3 3 2 3 2 3
M =M +t, N =N +t
[0086] ステップ ST26にお 、て、検査行列生成部 26は上記ステップ ST25で最適化され た複数の符号ィヒ率に対する検査行列生成の次数分布により、上記ステップ ST23で 生成された擬似乱数置換行列を配置して以下に示すような LDPC符号用の検査行 列を生成する。
例えば、 R(l)=l— M /N =1/3,R(2)=1-M /N = 1/2, R(3) = 1-M
1 1 2 2 S
/N =2Z3に対して、上記ステップ ST25で求めた次数分布に従って、下記のよう
3
な検査行列を生成する。
[0087] [数 28]
Figure imgf000029_0001
[0088] 次に上記で説明した LDPC符号の特性を比較する。
第 9図はこの実施の形態 2による検査行列生成装置で生成された検査行列を使用 した LDPC符号の情報 1ビット当りの信号電力対ノイズ電力比 (EbZNo)に対するビ ット誤り率 (BER)特性とフレーム誤り率 (FER)特性を示す図である。ここでは、情報 長 1536ビット、通信路 AWGN、最大繰り返し数 200回としている。なお、復号法は「 Sum— Productアルゴリズム」である。この特性は、下記に示すパラメータ M— M、
1 3
N— Nを用いてこの実施の形態 2の手順に従って生成した RC— LDPC符号の検査
1 3
行列を使用したものである。この検査行列において、 R(l)=l— M ZN =1/3, R
(2)=1— M /N =1/2, R(3)=l-M /N =2/3である。
[0089] 一方、 R(l) =1Z3、 R(2) =1Z2、 R(3) =2/3におけるシャノン限界の SNRは 、それぞれ— 0.4954 (dB), 0.1870 (dB), 1.0595 (dB)であり、それぞれの差は BER=10— 4点で、約 1.8 (dB) , 1.7(dB)ゝ 1.45(dB)となり、情報長 1536ビットと それほど長くないにもかかわらずシャノン限界に接近する RC— LDPC符号が実現で きている。
[0090] 以上のように、この実施の形態 2によれば、実施の形態 1と同様の効果が得られると 共に、実施の形態 1では符号ィ匕率毎に個別に検査行列を生成する必要があるが、こ の実施の形態 2では複数の符号ィヒ率に対して最終的に一種類の検査行列を用意す るだけで良く、検査行列を容易に生成することができると!/、う効果が得られる。
[0091] 実施の形態 3. 第 10図はこの発明の実施の形態 3による LDPC符号用の検査行列生成装置の構 成を示すブロック図である。この検査行列生成装置は、情報長'符号化率'列最大次 数入力部 31、擬似乱数置換行列用パラメータ算出部 32、擬似乱数置換行列生成 部 33、次数分布最適化用パラメータ算出部 34、次数分布最適化部 35及び検査行 列生成部 36を備えている。この検査行列生成装置は、第 4図における送信側及び受 信側の通信装置内に搭載され、生成した検査行列を符号化器 101及び復号器 105 に出力しても良いし、送信側及び受信側の通信装置外に設置され、生成した検査行 列を送信側及び受信側の通信装置内のメモリ（図示せず）に格納し符号化器 101及 び復号器 105に出力しても良い。
[0092] 第 11図はこの発明の実施の形態 3による検査行列生成装置の処理の流れを示す フローチャートである。ステップ ST31において、情報長'符号化率'列最大次数入力 部 31は所定の情報長、符号ィ匕率及び列の最大次数を入力する。ステップ ST32に おいて、擬似乱数置換行列用パラメータ算出部 32は上記ステップ ST31で入力され た所定の情報長、符号化率及び列の最大次数を用いて、後述の第 1及び第 2の擬 似乱数置換行列のためのパラメータを算出する。
[0093] ステップ ST33において、擬似乱数置換行列生成部 33は、上記ステップ ST32で 算出された第 1及び第 2の擬似乱数置換行列のためのパラメータにより擬似乱数系 列とラテン方陣行列を使用して第 1の擬似乱数置換行列を生成すると共に、第 1の擬 似乱数置換行列をさらに行方向に圧縮して第 2の擬似乱数置換行列を生成する。ス テツプ ST34において、次数分布最適化用パラメータ算出部 34は上記ステップ ST3 1で入力された所定の情報長、符号化率及び列の最大次数を用いて、第 1及び第 2 の擬似乱数置換行列を使用して構成可能な検査行列生成の次数分布の最適化の ための取りうる次数分布の組み合わせのパラメータを算出する。
[0094] ステップ ST35において、次数分布最適化部 35は、上記ステップ ST34で算出され た取りうる次数分布の組み合わせを拘束条件として、密度発展法により検査行列生 成の次数分布を最適化する。ステップ ST36において、検査行列生成部 36は上記ス テツプ ST35で最適化された検査行列生成の次数分布により、上記ステップ ST33で 生成された第 1及び第 2の擬似乱数置換行列を配置して LDPC符号用の検査行列 を生成する。
[0095] 次に第 11図のフローチャートにおけるステップ ST31— ST36の詳細処理について 説明する。
ステップ ST31にお 、て、情報長 ·符号化率 ·列最大次数入力部 31は所定の情報 長、符号化率及び列の最大次数として、システムで事前に設定されている情報長 K = 120、符号化率 rate = 0.5、列の最大次数 d= 16を入力する。
[0096] ステップ ST32において、擬似乱数置換行列用パラメータ算出部 32は、上記ステツ プ ST31で入力された所定の情報長 K=120、符号化率 rate = 0.5及び列の最大 次数 d= 16を用いて、以下に示すように、検査行列 Hの行数 M、列数（=符号化率） N、第 1及び第 2の擬似乱数置換行列のための共通のパラメータ (p'bXp'bの疑似 乱数置換行列の Pと b)を算出する。
[0097] [数 29] 後述の第 1の擬似乱数 s 一 I'BPS J!,BPS _m,BPShv _ινの行方向の最大個数 は 4であるため （ステップ ST33— 7にて説明） 、 残る列次数 4=12を下記ステツ プ S T 34にて述べる第 2の擬似乱数置換行列 , c!i , civ で cw=3によって実 現する。 cw=3の場合、 擬似乱数置換行列 H 'cin の使用個数 =残る列次数 12/cw=4となり cし, cH の 4個の組み合わせを列方向に並べて用いる。
符号 β=検査行列 IIの列数 N=k/rate=120/(l/3)=360、
検査行列 Hの行数 M = N*(1- ratc)=360*(2/3)=240、
擬似 乱数 置 換 行 列 の サ イ ズ P · δ = M / (擬似 乱 数 置 換 行 列
BPS —I BPS^一 11, BPS —!11, BPS _/ の使用個数 +擬似乱数置換行列
w„ / m>cff ， c/ の使用個数） =240/ 30.によって例えば p =6， ύ =5となる。
[0098] ステップ ST33において、擬似乱数置換行列生成部 33は、以下のステップ ST33— 1一 ST33— 8に示すように、上記ステップ ST32で算出された第 1及び第 2の擬似乱 数置換行列のためのパラメータにより擬似乱数系列とラテン方陣行列を使用して第 1 の擬似乱数置換行列を生成すると共に、第 1の擬似乱数置換行列をさらに行方向に 圧縮して第 2の擬似乱数置換行列を生成する。 すなわち、 上記ステップ ST 32で求めた J?と bから第 1の擬似乱数置換行列
BPS _I,BPSllljp—II, BPS —W'BPS _/K及び第 2の擬似乱数置換行列
cl^dl^cIII^cIV^を生成する。
[0099] [ステップ ST33—1]
[数 31] 以下の式で定義される擬似乱数系列 を生成する。 このィンデックスの,',, は
C,„, (1)の値、 つまり初期値を示す。
Cmu{\) = init, C(i' + 1) -G„ xC(i)mod P, i =\,2,...,P -2,
ここで、 Pは/ >>pとなる素数であり、 は の原始元である。 なお、 GF(P)の 要素である (^の連続するべき乗の値 (例えば （7' ,ίϊ2,...,( ρ-' =1)が (?F(P)の非零の全ての 要素を示す場合、 G„を原始元とよぶ。 この系列 C,„,,(,)は要素間の差分が非一様であり ^複する要素がない擬似乱数系列となる。
[0100] [数 32] 例として ^7，^ 3，"ii =3の場合、
： 1 2 3 4 5 6
C )： 3 2 Λ 4 5 I
次に C,.,„ (i)の要素から pより大きい数を間引きする。
例えば/, =5の場合
： 1 2 3 4 5
, ('.)： 3 2 4 5 1
[0101] [ステップ ST33— 2]
[数 33]
下式により擬似乱数系列及びその置換パターンを要素にもつラテン方陣行列 i; (q,l) を生成する。
L(i/,to+i-2)mod ))+l) = C( ,
= 1,2,..., p, q -1,2，..， p.
この操作により、 を使ってラテン方陣行列を生成することができる。
P = 6の場台の例を以下に示す。
1 2 3 4 5 6
L(l,/) 3 2 6 4 5 1
L(2,l) 2 6 4 5 1 3
L(3,l) 6 4 5 1 3 2
L(4 ) 4 5 1 3 2 6
1(5,/) 5 1 3 2 6 4
L(( l) 1 3 2 6 4 5
例えば p =5の場合は以下となる。
1 2 3 4 5
3 2 4 5 1
0 2 4 5 1 3
；) 4 5 1 3 2
/) 5 1 3 2 4
；) 1 3 2 6 4
i) 3 2 6 4 5
[0102] [ステップ ST33— 3]
[数 34] ド式により基本置換行列パターン SP(<7,,)を生成する。
'· =】，2 ,·.·， P, q = 1,2,.., p.
p =6の場合の例を以下に示す。
： 1 2 3 4 5 6
BP(U) 3 2 4 5 1
P(2,i) 8 12 10 11 7 9
P(\i) 18 16 17 13 15 14
BP(4 ) 22 23 19 21 20 24
尸 (5，り 29 25 27 26 30 28
BP(6 ) 31 33 32 36 34 35
[0103] [ステップ ST33— 4]
[数 35] 下式により基本置換行列パターン HP(q，i)の行単位の一定間隔の置換した行列 8¾,( ,0を生成する。 なお、 / には、 行単位の置換を行う際の間隔を設定する。 i) ftと topが互いに素の場合、
ΰ/'„„„ (q,i) = BI {h p mod b) + 1, i)
i i) と topが互いに素でない場合、
ΒΡΙνψ (q,i) ~ BP((hop * (<? - 1) mod b) + 1 + / ύ J, i )
b = 5Mp = 2の場合の例を以下に示す。
t： 1 2 3 4 5 6
3 2 6 4 5 1
IS 16 17 13 15 14
Κ3，/) 29 25 27 26 30 28
8 12 10 11 7 9
Κ5，) 22 23 19 21 20 24
[0104] [ステップ ST33— 5]
[数 36] 下式により基本置換行列パターン BPh。p (q, i)から基本置換系列 BPS (,)を生成する。 i i - 1,2,.., p.
この操作により、 p.bxp.fcの列次数 1 , 行次数 1の行列を生成するための基本 置換行列の " 1 " の位置の行に対する列の位置ベクトルを生成できる。
[0105] 例を以下に示す。
[数 37]
/： 1 2 3 4 5 6
BPS,,,, (0 3 18 29 8 22 2 i： 7 8 9 10 11 12
BPS,,,,,, (り 16 25 12 23 6 17
： 13 14 15 16 17 18
BPS (0 27 10 19 4 13 2β
： 19 20 21 22 23 24
BPS 0) 11 21 5 15 30 7 i： 25 26 27 2« 29 30
BPSiltJ{> (0 20 1 14 28 9 24
[0106] [ステップ ST33— 6] [数 38] 基本置換系列 BPS )から下式により定義される 2元 p A X P .f>行列 BPS '( )を 生成する。
BPSh m)
0, otherwise
この操作により、 わ xp.fcの列次数 1、 行次数 1の行列を生成するための基本置換 行列を生成できる。 ここまでの操作により、 隣接する要素間の距離が広くかつ要素間 の距離が非一様な擬似乱数置換行列を生成できる。
隣接する要素間の距離を広げることで実施の形態 1で述べた Tの下三角行列との間 で発生する短いサイクルを抑制できる。 また、 要素間の距離を非一様にすることで、 上 dp. p.fcの基本置換行列内においても短いサイクルの発生を抑制できる。 また、 サイクル 4の発生を回避するために生成した p.fcxp.hの基本置換行列を確認し、 サイ クル 4が存在していた場合には、 p, b, G。， init, hopを変動させ、 サイクル 4が存在し ない組み合わせを探索し再生成する。 但し、 本実施の形態の手順ではサイクル 4が存 在することは稀であり、 再生成を必要とすることはほとんどない。
例を以下に示す。
Figure imgf000035_0001
[ステップ ST33—7]
[数 39] 基 本 置 換 行 列 spsh„„ か ら 4 種 の 基 本 置 換 行 列
BPS, 'BPS , _II,BPS. —UI,BPS IVを生成" 9る
BPS^, _I{p-b-BPShllll(i)+\,i) = l i =1 , p-b
その他の β尸 5 /の要素はすべて 0である。
i i) BPS— _U{j,i)~BPShop _I ,p b-j+i), 1,2 "..， p b, Hp.b i i i ) BPS _ III {i A) = BPS—― p-b-j + 1), i = 1,2，...， p-b, 卜 1,2，...， p b ί ) BPShv _IV{i,i)-BPShrp _III(i,p-b^j+l), i~ 1 ., P-b, ; =1,2 p-b
この i i), i ί i), iv)は SP „,,—/を元に行と列の置換を行ったのもであり、 全て p.bxp.bの列次数 1 行次数 1の行列である。 これらの行列の組み合わせにより列次 数丄 X行次数 1から列次数 4 X行次数 4までの行列を生成できる。 この行と列の置換により異なった擬似乱数置換行列を 4種類生成できる。
[ステップ ST33—8]
次に列次数 1より大きい基本置換行列を生成する。
[数 40]
ステップ S T 3 3— 7で求めた基本置換行列 sra _ /を用いて列次数 cw、 行次数 1 の （p.i>)x (P.fc / v)擬似乱数置換行列 c を以下の手順により生成する。 このとき、 列 が I/cwに圧縮され、 つまり行方向に圧縮される。 cに DPS,..
Figure imgf000036_0001
」 で、 Jは (p 'fc / ）の単位行列である。
同様の手順により基本置換行列 S/^。„ _ II,BPS^ _ lU,BPShv—IVを用いて列次数^ 行次数 1の基本置換行列 CIし を生成する。 但し、 この c
Figure imgf000036_0002
の組み合わせは列方向に 2個以上並べた場合、 サイクル 4が発生することがわかつて いる為行方向にのみ 1列に並べられるものとする。
[0109] ステップ ST34にお 1、て、次数分布最適化用パラメータ算出部 34は上記ステップ S T31で入力された所定の情報長、符号ィ匕率及び列の最大次数を用いて、以下のよう にして、第 1及び第 2の擬似乱数置換行列を使用して構成可能な検査行列生成の次 数分布の最適化のための取りうる次数分布の組み合わせのパラメータを算出する。
[0110] [数 41]
これら BPS^'—I' BPS —li'BPS一 _ III,BPS^ _ Wと c し の組み合わせと と右卜-三角 0の三角行列を用い、 実現可能な組み合わせを求める。
例えば下記検査行列 I Iに対し左端から 4 X 4個の擬似乱数置換行列
B j, BPSl:,„, _ !I, BPS,,,„, j!I，BPS,,,,,,—IVと 4 X 1個 （CW=3の場合） の擬似乱数置換行 列 t ， c//,„，c///,„, /i^と Hの右側の左下三角部分に 7 + 6十… 1個の擬似乱数置換行 列 !, BPS _ll,BPS l— 'BPS _ IVが配置可能となる。
•BPS - /// BPS —u — I BPS,,,,,, _ IV \ 0 0 0 0 0 0 0
BPS —II BPS _ 1 BPS - - IV _ "1 0 \ 0 0 0 0 0 0
BPS _ l BP'S _ IV BPS - III BPS , II 0 0 \ 0 0 0 0 0
BPS ― IV —Ul BPS —II _ I 0 0 0 \ 0 0 0 0 c ,„. 0 0 0 0 0 0 0 \ 0 0 0
0 0 0 0 0 0 0 0 \ 0 0
0 0 0 0 0 0 0 0 0 \ 0 c/ 0 0 0 0 0 0 0 0 0 0 \
[0111] ステップ ST35において、次数分布最適化部 35は、実施の形態 1のステップ ST15 と同様にして、上記ステップ ST34で算出された取りうる次数分布の組み合わせを拘 束条件として、密度発展法により検査行列生成の次数分布を最適化する。
[0112] ステップ ST36において、検査行列生成部 36は上記ステップ ST35で最適化され た検査行列生成の次数分布により、上記ステップ ST33で生成された第 1及び第 2の 擬似乱数置換行列を配置して LDPC符号用の検査行列を生成する。
[数 42] すなわち、 ステップ S T 3 5で求めた最適な次数分布の組み合わせを用いて、 第 1 の擬似乱数置換行列 θ „,, _ I, BPSh<„,—II, BPS, —III, BPS —IV及び第 2の擬似乱数置換 行列 cし を用いて検査行列 Hを生成する。 各擬似乱数置換行列の配置 に関しては探索等によりサイクル数を確認して良好な配置を求める。
[0113] 以上のように、この実施の形態 3によれば、実施の形態 1と同様の効果が得られる。
[0114] 実施の形態 4.
第 12図はこの発明の実施の形態 4による LDPC符号用の検査行列生成装置の構 成を示すブロック図である。この検査行列生成装置は、情報長'符号化率'列最大次 数入力部 41、擬似乱数置換行列用パラメータ算出部 42、擬似乱数置換行列生成 部 43、次数分布最適化用パラメータ算出部 44、次数分布最適化部 45及び検査行 列生成部 46を備えている。この検査行列生成装置は、第 4図における送信側及び受 信側の通信装置内に搭載され、生成した検査行列を符号化器 101及び復号器 105 に出力しても良いし、送信側及び受信側の通信装置外に設置され、生成した検査行 列を送信側及び受信側の通信装置内のメモリ（図示せず）に格納し符号化器 101及 び復号器 105に出力しても良い。
[0115] 第 13図はこの発明の実施の形態 4による検査行列生成装置の処理の流れを示す フローチャートである。
[数 43]
¾施の形態 3では、 BRS — I , BPS — Π , BPS — III , BPS — IVを使って、 例えば、 列 次数 4 X行次数 4の行列を生成しているが、 実施の形態 4では、 後述する £IH尸 を使って直接列次数 ' X行次数 ™' の行列を生成する。
[0116] ステップ ST41において、情報長'符号化率'列最大次数入力部 31は所定の情報 長、符号ィ匕率及び列の最大次数を入力する。ステップ ST42において、擬似乱数置 換行列用パラメータ算出部 42は上記ステップ ST41で入力された所定の情報長、符 号ィ匕率及び列の最大次数を用いて、後述の第 3の擬似乱数置換行列のためのパラ メータ及び実施の形態 2と同じ第 2の擬似乱数置換行列のためのパラメータを算出す る。
[0117] ステップ ST43において、擬似乱数置換行列生成部 43は、上記ステップ ST42で 算出された第 3の擬似乱数置換行列のためのパラメータ及び第 2の擬似乱数置換行 列のためのパラメータにより、擬似乱数系列とラテン方陣行列を使用して第 1の擬似 乱数置換行列を生成し、生成した第 1の疑似乱数置換行列の行と列をさらに置換し 拡張して第 3の擬似乱数置換行列を生成すると共に、生成した第 1の擬似乱数置換 行列をさらに行方向に圧縮して第 2の擬似乱数置換行列を生成する。ステップ ST44 において、次数分布最適化用パラメータ算出部 44は上記ステップ ST41で入力され た所定の情報長、符号化率及び列の最大次数を用いて、検査行列生成の次数分布 の最適化のための取りうる次数分布の組み合わせのパラメータを算出する。
[0118] ステップ ST45において、次数分布最適化部 45は、上記ステップ ST44で算出され た取りうる次数分布の組み合わせを拘束条件として、密度発展法又は Multi— edge typeの密度発展法により検査行列生成の次数分布を最適化する。ステップ ST46に ぉ 、て、検査行列生成部 46は上記ステップ ST45で最適化された検査行列生成の 次数分布により、上記ステップ ST43で生成された第 3及び第 2の擬似乱数置換行列 を配置して LDPC符号用の検査行列を生成する。
[0119] 次に第 13図のフローチャートにおけるステップ ST41— ST46の詳細処理について 説明する。
ステップ ST41において、情報長'符号化率'列最大次数入力部 41は所定の情報 長、符号化率及び列の最大次数として、システムで事前に設定されている情報長 K = 120、符号化率 rate = 0. 5、列の最大次数 d= 16を入力する。
[0120] ステップ ST42において、擬似乱数置換行列用パラメータ算出部 42は、上記ステツ プ ST41で入力された所定の情報長 K= 120、符号化率 rate = 0. 5及び列の最大 次数 d= 16を用いて、以下に示すように、検査行列 Hの行数 M、列数（=符号化率） N、第 3の擬似乱数置換行列のためのパラメータ及び第 2の擬似乱数置換行列のた めのパラメータを算出する。
すなわち、ここでは、
[数 44] ステップ S T 4 1で入力されたデ一タを用いて、 検査行列 Hの行数 Mを列数 （=符 ^長） N、 ·δ'χ . の列重み CH 、 行重み rw' となる第 3擬似乱数置換行列 EL f .. の と を算出し、 0? .6)x . の第 2の擬似乱数置換行列 cし u civr„の！), を舁出する。
[0121] 例えば、以下のように、第 3の擬似乱数置換行列のためのパラメータ及び第 2の擬 似乱数置換行列のためのパラメータを算出する。
[数 45] 実施の形態 3のステップ ST 3 2の計算と同じ計算を行った後、 擬似乱数置換行列 BPS — BPS —Il'BPS —m,BPS _wの ^ Άを、 — F言己のステップ S T4 3― 1 ~ S T 4 3 ~ 9で説明する擬似乱数置換行列 EL JP . ) に置き換える。
擬似乱数優換行列 — n， _ινの行方向の最大個数は 4で あるため （実施の形態 3のステップ S T 3 3— 7にて説明） 、 残る列次数 4=12を、 以下のステップ S T 4 4にて述ぺる擬似乱数 β換行列 で cw=3によ つて実現する。 cw=3の場合、 擬似乱数置換行列 ^/ /^ の使用個数：残る 列次数 12/cw=4となり C/OT / /„,cj7/„,c/ „の 4個の組み合わせを列方向に並べて用いる。 符号長=検查行列^の列数 =^316=120/(1/3)=360.
検査行列 Hの行数 M=NX (l-rate)=360X (2/3) =240,
擬 似 乱 数 置 換 行 列 の サ イ ズ P · ύ = M / ( 擬 似 乱 数 置 換 行 列 — —n'BPS —m'BPS —rv の 使 用 個 数 + 擬 似 舌し 数 置 換 行 列 w„,C „„,c//J„,c: K„の使用個数） =240/8=30. よって例えば、 第 1及び第 2の擬似乱 数置換行列を生成するためのパラメータは、 p =6, b =5となる。
ここまでは実施の形態 3で求めた ί<と bのィ直であり. 以下のような検査行列 Ι·Ιを想 定している。
Figure imgf000040_0001
■この検査行列 Hの擬似乱数置換行列 BPS. BPS ,BrShap— BPS で される例えば. 列次数 4 X行次数 4で
'
Figure imgf000040_0002
の 4-P'bx4-p-b 部分行列を
以下の検査行列 Hのような列次数 cw' = 4 rw' = の擬似 ¾L数 β換行列 p' b'xp' b' = 4- j bx4 p b の ELHP^ に： ffi Ϊ き換える構造とする。
\ 0 0 0 0 0 0 0
0 \ 0 0 0 0 0 0
0 0 \ o 0 0 0 0
c 0 0 0
0 0 0
H = 0 0 0 \ 0 0 0 0
0 0 0 0
ΰ 0 0 \ 0 Q 0
0 0 0 o 0 0 0
0 0 0 \
0 0 0 0 0 o \ 0
0 0 0 0 0 0 0 \
よって、 例えば、 jp 6 b 5力 ら ^ - P-b -- = 4- 30 = 120 となる。 した;^つ て，例えば、 第 3の擬似乱数置換行列を生成するためのパラメータは， ' = 12, ' = W となる。 ステップ ST43において、擬似乱数置換行列生成部 43は、上記ステップ ST42で 算出された第 3の擬似乱数置換行列のためのパラメータ及び第 2の擬似乱数置換行 列のためのパラメータにより、擬似乱数系列とラテン方陣行列を使用して第 1の疑似 乱数置換行列（第 2の擬似乱数置換行列のためのパラメータと共通のパラメータによ り生成される）を生成し、生成した第 1の擬似乱数置換行列の行と列をさらに置換し 拡張して第 3の擬似乱数置換行列を生成する。
[数 46] すなわち、 ここでは、 ステップ ST42で求めた と から、 下記のステップ ST 4 3— l〜ST43— 9により，第 3の擬似乱数置換行列 £LH を作成する。
[0123] また、このステップ ST43において、第 1の擬似乱数置換行列をさらに行方向に圧 縮して第 2の擬似乱数置換行列を生成する方法は、実施の形態 3のステップ ST33— 1一 ST33—8までの処理と同じである。
[0124] [ステップ ST43— 1]
[数 47] 以下の式で定義される基本乱数系列 c , を生成する。 このインデックスの ,, は C,„, (1)の値つまり初期値を示す。
CM,(l)^i t, CMt(i + l)-G0xCM(i)mod P, i 1,2,...,/»— 2,
ここで Pは尸 > となる素数であり、 G。は (尸）の原始元である。 なお、 (尸） の要索である（7„の連続するべき乗の値 (例えば G' .G1，..., -' =1)が G/'，(尸）の非零の全て の耍素を示す場合、 を原始元とよぶ。 この系列 (りは要素間の差分が非一様であ り重複する要素がない擬似乱数系列となる。 次に C,„„(i)の要素から P より大きい数 を冏引きする。
例として尸 = 7， Gu = 3,init = 3,ρ' = 6の場合、
1 2 3 4 5 6
Figure imgf000041_0001
[0125] [ステップ ST43— 2]
下式により擬似乱数系列及びその置換パターンを要素にもつ行列を生成する [数 48] i = X2,...,p', j 1,2,..，//.
/ = 6の場合の例を以下に示す。
/： 1 2 3 4 5 6
LH( ) 3 2 6 4 5 1
LHil,l) 4 3 5 2 1 6
L (3,l) 5 6 3 1 4 2
LH(4,l) 2 4 1 3 6 5
LH(5,l) 5 1 2 5 3 4
LH(G,l) 1 5 4 6 2 3
[0126] [ステップ ST43— 3]
[数 49] 下式により列次数 CH ,行次数 nv'の基本置換行列パターン £H/^» (ムり を生成 する。
LHP(cw._r„.} ((!■ - 1) *p'lcw' +ムり = LH(i,(j -l)*cw' + k),
ι = 1,2,..., rvc , j = l,2,..,p lew', k =l,2,..,cw'. p' = 6の場合の例を以下に示す。 ： 1 2 3
LHPai) l,i) 3 2 6
(w)(2 ) 4 5 1
LHP^(3,i) 4 3 5
UIPa3)(i,i) 2 1 6
ΙΗΡί })(5, 5 6 3
L PflJ)(6,t) 1 4 2
[0127] [ステップ ST43— 4]
[数 50] ト-記ステップ ST43— 3までで、 P'xp' の基本置換行列パターンを生成した。 次にこれを行及び列に関して b' 倍する手順を説明する。
ステップ ST43 1と同じ手順で基本乱数系列 ,, (りを生成する。 この ,, は ,,,,.,(0に対し P' を b' に変えて生成した基本乱数系列である。
例として 尸 = 7,GD =3，i"ii =3 ' =6 の場合、
；: 1 2 3 4 5 6
'(i): 3 2 6 4 5 1
[0128] [ステップ ST43— 5] [数 51] 下式により擬似乱数系列及びその置換パターンを要素にもつラテン方陣行列 (q ) を生成する。
L" ((1 +【· - 2) mod (ρ')) + 1) = *;, ('.),
2 =1,2,...,^', q=X2,..,p'.
この操作により、 を使ってラテン方陣を生成することができる。
p' = 6の場合の例を以下に示す。
/: 1 2 3 4 5 6
L (l,/) 3 2 6 4 5 1
L''(2,l) 2 6 4 5 1 3
L''{3,1) 6 4 5 1 3 2
£"(4,/) 4 1 3 2 6
£''(5,/) 5 1 3 2 6 4
"(6 ) 1 3 2 6 4 5
[0129] [ステップ ST43—6]
[数 52] ド式により、 列次数 ,行次数 rM の基本置換行列パターン^^ ·. .)( ，りを '倍に拡張するための行列 '^ ,™') (ムり を生成する。
', ((!■― 1) *p'lc + j,k) = lb(W一 1) * cw' + k),
i l,2,...,rw， j = 1,2,.., p lew', k = l,2,..,cw'.
P' = S の場合の例を以下に示す。
'■: 1 2 3
―,) ( ) 3 2 6
.,,(2,0 4 5 1
,,(3,0 2 6 4
(4,'·) 5 1 3
,,(5,0 6 4 5
, 1 3 2
[0130] [ステップ ST43— 7 I
[数 53] 下式により、列方向に拡張する基本置換行列パターン H/^w'
^LHPM (i, j) = LHPM) (i, j) + ((L"P(33) (i, j)― 1) * ρ'),
i = …' *nv'/cw , j = X2,..,cw'.
¾ ¾
以下に P' = 6,6' = 6，cw' = 3，m = 3 の場合の例を示す。
/： 1 2 3
)( ） 15 8 36
2，り 22 29 1
LHPi ,( .） 10 33 23
,(4,0 26 1 18
,(5，り 35 24 27
eL P< .(6, 1 16 8
[0131] [ステップ ST43— 8]
[数 54] ド式により，行方向に拡張する基本置換行列パターン
Figure imgf000044_0001
を生成する。
ELHP(M ((j ~\)-b' + i,k) = ((eLHP^.^ (j, ) + (i - 1) ' - 1) mod p'-b^ + l,
i = 1,2,... ,o , j = 1,2,.., p' -rw lew', k - 1,2,.., cw'. 以下に P' - 6, = 6,cw' = 3,rw' = 3の場台の例を示す。
2 3
,(1. 15 8 36
21 14 6
,,(3,0 27 20 12
Figure imgf000044_0002
£"/尸 („)(36,') 31
[0132] [ステップ ST43—9]
[数 55] 本置換パターン iH尸 ) , )から、 下式により定義される 2元 -Λ' 行列 jipi→w (km jを生成する。
[0, otherwise
この操作により、 P' 'b'xp' -b' の列次数 ν'、 行次数/ τν'の行列を生成するため の基本置換行列を生成できる。 ここまでの操作により列次数 cvv'、 行次数 πν'の擬似 置換行列を生成できる。 この擬似乱数置換行列の列次数 cvv'、 行次数 rw'はいずれも 4以下でサイクル 4が存在しない場合が多いため、 列次数 πν'、 行次数 ν'は最大値 を 4とする。
[0133] ステップ ST44において、次数分布最適化用パラメータ算出部 44は上記ステップ S T41で入力された所定の情報長、符号ィ匕率及び列の最大次数を用いて、下記に示 すように、検査行列生成の次数分布の最適化のための取りうる次数分布の組み合わ せのパラメータを算出する。
[0134] [数 56] これら 尸 ^.)と c , ， c / の組み合わせと右上三角 0の三角行列を用 い， 実現可能な組み合わせを求める。
例えば下記検査行列 Hに対し左端から 1個の擬似乱数置換行列 ELHP と 4 X 1個 （cw=3の場合） の擬似乱数置換行列 c/w，w , c/ / , c /[^が配置可能となる。
\ 0 0 0 0 0 0 0
0 \ 0 0 0 0 0 0
Hi . 4,
0 0 \ 0 0 0 0 0
c 0 0 0
0 0 0 \ 0 0 0 0
H = 0 0 0
0 0 0 0 \ 0 0 0
c" 0 0 0
0 0 0 0 0 \ 0 0
cl 0 0 0
0 0 0 0 0 0 \ 0
0 0 0 0 0 0 0 \ この取りうる次数分布の組み合わせ力 次数分布の最適化の拘束条件となる。 ステップ ST45において、次数分布最適化部 45は、上記ステップ ST44で算出され た取りうる次数分布の組み合わせを拘束条件として、下記に示すように、密度発展法 又は Multi - edge typeの密度発展法により、検査行列生成の次数分布を最適化 する。 [0136] すなわち、ここでは、上記ステップ ST44で算出された検査行列生成の次数分布の 最適化のための取りうる次数分布の組み合わせの拘束条件のもと、実施の形態 1の ステップ ST15と同様の次数分布の最適化により最も良い次数分布を求めても良いし 、又は、以下に述べる Multi— edge typeの密度発展法でも良い。
[0137] 上記実施の形態 1で述べた密度発展法は、 Irregular LDPC符号アンサンブル（ 次数の分布を示すもの）が行次数の割合、列次数の割合で表現されることに対して、 Multi-edge typeのアンサンブルは「枝 (枝タイプ）の組合せ」（次数タイプ）の割合 で表現される。なお、 Multi— edge typeのアンサンブルの枠組みでは、どのビットを パンクチヤする力まで規定したアンサンブルを考えることができ、通信路タイプ別の考 慮が可能である。（"Multi- Edge Type LDPC Codes DRAFT", T.Richardson, R.Urbanke, web上)
[0138] 第 14図は Multi— edge typeのアンサンブルの例を示す図である。通信路タイプ（ 第 14図の例の場合は BEC (Binary Erasure Channel)の通信路（b = (1 0)か、分散 値 σ 2のガウス分布に従うノイズが加えられた通信路、いわゆる加法性白色ガウス雑 音 (AWGN)の通信路 (b = (0 1) )を bとし、次数タイプ dのバリアブルノードの符号 長 Nに対する割合を V 、チヱックノードの符号長 Nに対する割合を; z とすると、第 1 b，d d
4図は第 15図に示す表のように表現される。
第 15図は irregular LDPC符号の次数分布の例を示す図である。
[0139] 以下に具体的な検査行列の例を示す。例えば第 15図の表の Variable (バリアブル ノード:行列の各列がそれぞれバリアブルノードになる） b = (l 0)の V の合計は 1 b，d
である。第 15図の表の Variableの 1番上の比率は V =0. 5であり、 d= (2 0 0 b，d
0)であるが、これは第 15図の次数分布に従つて構成した第 16図に示す検査行列の 15列目力も 24列目を示している。第 15図の表の Variableの 3番目には、 v =0.
b,d
2、 b= (l 0)、d= (0 3 3 0)が示されており、第 16図の検査行列の 11列目から 14列目を示している。この部分は消失確率 1の通信路を意味しており、現実にはこの 部分の符号ビットはパンクチヤされて送信されない。受信側では、この部分の計算は 0と 1の確率が 1Z2という値で計算される。受信値の LLR (対数尤度比）としては 0と なる。 第 15図の表の Constraintはチェックノード（行列の各行がそれぞれチェックノード) の比率を表している。第 15図の表の Constraintの 1番上の比率は =0. 4であり、 d
d= (4 1 0 0)である力 これは第 16図の検査行列の 5行目力も 12行目を示して いる。
[0140] なお、タイプ iの枝のうち、次数タイプ dのノードに接続している枝の割合は、
[数 57] バリアブルノー k' , « ^'Λ, (次数タイプ dの枝タイプ iの数） レ b,d
, レ b,d (枝タイプ iの全数） — ,（U)
チ ノ—ド . (次数タイプ dの枝タイプ ίの数 ά,η,
ノ · ^ヌメ (枝タイプ iの全数) = /I l)" で求められる。 （Standard Irregular LDPCのえ、 ρにあたるもの）
なお、 1は ALL1を示す。
[0141] 次に Multi— edge typeの密度発展法について説明する。 Multi— edge type の密度発展法では、ノ リアブルノード力 チェックノードへのメッセージを一緒にまと めずに、枝のタイプ毎に確率密度関数を分けて計算を進める ("Multト Edge Type LDPC Codes DRAFT", T.Richardson, R.Urbanke, web上)。
行処理（チェックノード力ゝらのメッセージ）
[数 58] 今枝タイプ全ての集合を Iとする。 今 i e Iに対して枝タイプ iの次数を 、 枝 タイプ i を通してチェックノードへ送信される確率密度関数を A.、 枝タイプ i を通 してチェックノードから送信される確率密度関数を/^とする。 次数タイプ毎に確率密度関数を計算し、枝タイプ iに対する枝の割合をかけ、足し 合わせることで、枝タイプ毎の確率密度関数を計算する。
[数 59] pu: = A, W =β Z mh-1 tanh-lanh-) ) ... (Al)
7 (1) ' \ 2ノノ
Vj e I \ /
また、 繰り返し復号における繰り返し数を /とおき、 枝タイプ Πこ対して繰り返し /回目のチェックノードのアンサンブル平均を ί:' とおくと、 以下を満たすことが sum-producl復号における繰り返し復号において漸近的に誤りが限りなく 0に近づく 必要条件である。
全ての枝タイプ/に対して ' W ,VIER --(A 2)
[0142] 列処理（バリアブルノードからのメッセージ）
[数 60]
通信路タイプ bl b2に対する次数を ί Μ dh2, …初期値を /„,,ti、 Pm …とする。 次数タイプ毎に確率密度関数を計算し、枝タイプ iに対する枝の割合をかけ、足し 合わせることで、枝タイプ毎の確率密度関数を計算する。
[数 61]
Figure imgf000048_0001
H® " J (®Ί,Η®" ) ·' … また、 繰り返し復号における繰り返し数を とおき、 枝タイプ;'に対して繰り返し/ 回冃のパリアブルノードのアンサンブル平均を mc とおくと、 以下を満たすことが sum- product復号における繰り返し復号において漸近的に誤りが限りなく 0に近づく 必 5?条件である。
全ての枝タイプ /に対して
>",^-., , ΕΛ÷. ---(A 4)
[0143] さらに、以下を満たすことが sum product復号における繰り返し復号において漸 近的に誤りが限りなく 0に近づく必要十分条件である。
式 (A2)又は式 (A4)の!、すせれかを満足すること。 · · (A5)
[0144] 次に、 Multi— edge typeの密度発展法における最適化を用いて、要求された符 号化率に基づく「Irregular— LDPC符号」のアンサンブル（次数配分）を求める。 通信路タイプ b、次数タイプ dのノリアブルノードの符号長 Nに対する割合を V 、チ エックノードの符号長 Nに対する割合を としそれぞれ同時に変数として扱い、 AW d
GN通信路（通信路タイプ AWGNの bl : =b= (01) )ガウスノイズ σ (下式）が最大 になるように、線形計画法で最適な V と を探索する。この際、この実施の形態で b,d d
は通信路タイプは AWGNの bl :=b=(0 1)と BECの b2 :=b= (l 0)のみとする。
[数 62]
Figure imgf000049_0001
[0145] この探索における拘束条件は、条件式 (A5)を満たすことである。
この際、通信路タイプ b、次数タイプ d、 V μ の変数のとりうる値は実施の形態 1 b，d、 d
力 実施の形態 4で述べた生成方法による検査行列 Hの制約条件で導出される値の みとする。
[0146] ステップ ST46において、検査行列生成部 46は上記ステップ ST45で最適化され た検査行列生成の次数分布により、下記に示すように、上記ステップ ST43で生成さ れた第 3及び第 2の擬似乱数置換行列を配置して LDPC符号用の検査行列を生成 する。
[数 63] ここでは、 ステップ S T 4 3で求めた最適な次数分布の組み合わせを用いて， 第 3 の擬似乱数置換行列 P(cw., 及び第 2の擬似乱数置換行列 c ， ， を 用いて検査行列 Hを生成する。 各擬似乱数置換行列の配置に関しては探索等によりサ イクル数を確認して良好な配置を求める。
[0147] 実施の形態 3と実施の形態 4の構成法による LDPC符号の性能評価結果を示す。
第 17図は評価に用いた irregular LDPC符号の次数分布の例を示す図である。
[0148] 第 18図は irregular LDPC符号のフレーム誤り率特性を示す図である。
[数 64] この第 1 8図は、 第 1 7図に示すこの次数分布を用いて、 実施の形態 3を用い
，尸- 17,i7,, -3，/™'i =9，toP = 2のパラメ一夕を用いて構成した情報長 5 7 6ビッ ト、 符号化率 0. 5の LDPC符号 （図中の Idpc576) と、 実施の形態 4を用い J^,^, のパフメータとして/) =48，わ = 30，尸 = 71，*^, = l,init = 13,cw - 4, w - 4、 l c„, II m,cIIl ,„,ciV,wの パラメ一夕として ρ -30，δ -12,/>-31,ίϊ„ = mit - = 1を用いて構成した情報長 144
0ビット、 符号化率 0. 5の LDPC符号 （図中の ldpcl440) と、 実施の形態 4を用い
ELHPm rw) の パ ラ メ 一 夕 と し て p UZ'fc -S ^-lS 'Gu -S.im'i lS.ov ^TO - A 、 〃 ，"//,,„ /„のパラメ一タとして 33,わ- 34,尸- 53, Gu -3,init - 9, hop -2を用いて 構成した情報長 3 1 68ビット、 符号化率 0. 5の LDPC符号 （図中の ldpc3168) と、 実施の形態 4を用い LW/^. )のパラメータとして
) = 204,b = 24,^ = 211,Ο,, =3,m ( -13,^ - 4,^ -4、 l m, II iw,cill m,dVlit,のノヽ。ラメ—夕として p 5 - 24 -53， -3,«, 9,/, 2を用いて構成した情報長 4896ビッ卜、 符号化 率 0. 5の LDPC符号 （図中の ldpc4896) の FER (フレーム誤り率） シミュレーション 結果を示している。
[0149] なお、第 18図中には比較のため第 3世代通信（3GPP)に用いられている turbo符 号の同一符号長、符号化率のものを比較で示してある（turbo576, turbo 1440, tu rbo3168, turbo4896)。第 18図を見ても分力、るように従来 10000ビット以下の短 V、符号長では turbo符号が優位とされてきた力 本発明により情報長 576ビットでほ ぼ同等性能を示し、それより長!、情報長では常に提案の LDPC符号の性能が優位と なっており、極めて高 、性能を示す符号構成法であることが確認できる。
[0150] 従来の方式の LDPC符号では、 10000ビット以上の符号長で良い性能を示してい る力 逆に符号長 10000ビット以下では turbo符号が優位と言われている。しかし、 第 18図では、情報長 4896ビット (符号ィ匕率 0.5なので符号長 9792ビット）で、 turb o符号より良い性能を示している。また、符号長が長くなるに従って turbo符号より性 能が良くなる傾向も示している。従来の LDPC符号ではこの 10000ビットくらいの符 号長で turbo符号と同等性能を示している力 この発明では、符号長 1000ビット程 度 (第 18図では情報長 576ビット、符号化率 0.5なので符号長 1152ビット)で turbo 符号と同等性能を示している。
[0151] 以上のように、この実施の形態 4によれば、実施の形態 1と同様の効果が得られると 共に、 Multi-edge typeの密度発展法を使用することにより、さらに誤り訂正能力 が良好な LDPC符号の検査行列を生成できるという効果が得られる。
[0152] 実施の形態 5.
第 19図はこの発明の実施の形態 5による LDPC符号用の検査行列生成装置の構 成を示すブロック図である。この実施の形態 5の検査行列生成装置は、擬似乱数置 換行列生成部 13、検査行列生成部 16、擬似乱数置換行列用パラメータ格納部 17 及び次数分布格納部 18を備えている。この検査行列生成装置は、第 4図における送 信側及び受信側の通信装置内に搭載され、生成した検査行列を符号化器 101及び 復号器 105に出力しても良いし、送信側及び受信側の通信装置外に設置され、生成 した検査行列を送信側及び受信側の通信装置内のメモリ（図示せず）に格納し符号 ィ匕器 101及び復号器 105に出力しても良い。
[0153] 第 19図において、擬似乱数置換行列用パラメータ格納部 17は、所定の情報長、 符号化率及び列の最大次数を用いて算出された擬似乱数置換行列のためのパラメ ータを格納している。この擬似乱数置換行列のためのパラメータは、外部で実施の形 態 1と同様にして生成される。擬似乱数置換行列生成部 13は、擬似乱数置換行列 用パラメータ格納部 17に格納されている擬似乱数置換行列のためのパラメータによ り、擬似乱数系列とラテン方陣行列を使用して擬似乱数置換行列を生成する。
[0154] また、次数分布格納部 18は、所定の情報長、符号化率及び列の最大次数を用い て算出された擬似乱数置換行列を使用して構成可能な検査行列生成の次数分布の 最適化のための取りうる次数分布の組み合わせを拘束条件として、密度発展法によ り最適化された検査行列生成の次数分布を格納して!/、る。この検査行列生成の次数 分布は、外部で実施の形態 1と同様にして最適化されている。検査行列生成部 16は 、次数分布格納部 18に格納されている最適化された検査行列生成の次数分布によ り、生成された擬似乱数置換行列を配置して低密度パリティ検査符号用の検査行列 を生成する。
[0155] この実施の形態 5では、実施の形態 1の擬似乱数置換行列生成部 13及び検査行 列生成部 16並びに擬似乱数置換行列用パラメータ格納部 17及び次数分布格納部 18により構成しているが、実施の形態 2の擬似乱数置換行列生成部 23及び検査行 列生成部 26並びに擬似乱数置換行列用パラメータ格納部 27及び次数分布格納部 28 (図示せず）により構成しても良ぐ実施の形態 3の擬似乱数置換行列生成部 33 及び検査行列生成部 36並びに擬似乱数置換行列用パラメータ格納部 37及び次数 分布格納部 38 (図示せず）により構成しても良ぐ実施の形態 4の擬似乱数置換行列 生成部 43及び検査行列生成部 46並びに擬似乱数置換行列用パラメータ格納部 47 及び次数分布格納部 48 (図示せず）により構成しても良い。
[0156] 以上のように、この実施の形態 5によれば、実施の形態 1と同様の効果が得られる。
[0157] 実施の形態 6.
第 20図はこの発明の実施の形態 6による通信装置の構成を示すブロック図である。 送信側の通信装置 201は、実施の形態 1から実施の形態 5のいずれかの検査行列 生成装置 100と、符号化器 101と変調器 102を備えている。また、受信側の通信装 置 202は、実施の形態 1から実施の形態 5のいずれかの検査行列生成装置 100と、 復調器 104と復号器 105を備えて 、る。
[0158] 送信側の通信装置 201において、検査行列生成装置 100は実施の形態 1から実 施の形態 5のいずれかの方法により検査行列 Hを生成し、符号化器 101は、生成さ れた検査行列 Hを使用して生成行列 Gを求め、そして、メッセージお u - --u )を受け
1 2 K 取り、求めた生成行列 Gを用いて符号語 Cを生成する。変調器 102は、生成した符号 語 Cに対して、 BPSK、 QPSK、多値 QAM等のデジタル変調を行って送信する。
[0159] 受信側の通信装置 202において、検査行列生成装置 100は実施の形態 1から実 施の形態 5のいずれかの方法により検査行列 Hを生成し、復調器 104は、通信路 10 3を介して受け取った変調信号に対して、 BPSK、 QPSK、多値 QAM等のデジタル 復調を行い、復号器 105が、 LDPC符号化された復調結果に対して、生成された検 查行列を用いて「sum— productアルゴリズム」による繰り返し復号を実施し、推定結 果 (もとの u U…！！ に対応）を出力する。
1 2 K
[0160] 以上のように、この実施の形態 6によれば、実施の形態 1から実施の形態 5のいず れかの検査行列生成装置 100を使用することにより、確定的で特性が安定し、誤り訂 正能力が良好で、かつ多様なアンサンブル、符号長、符号化率に対応した Irregula r LDPC符号用の検査行列を使用することができるという効果が得られる。
産業上の利用可能性 以上のように、この発明に係る検査行列生成装置及び検査行列生成方法は、例え ば、確定的で特性が安定し、誤り訂正能力が良好で、かつ多様なアンサンブル、符 号長、符号ィ匕率に対応した LDPC符号用の検査行列を短時間で容易に生成するも のに適している。

Claims

請求の範囲
[1] 所定の情報長、符号化率及び列の最大次数を入力する情報長'符号化率'列最大 次数入力部と、
入力された上記所定の情報長、符号化率及び列の最大次数を用いて後述の擬似 乱数置換行列のためのパラメータを算出する擬似乱数置換行列用パラメータ算出部 と、
算出された上記擬似乱数置換行列のためのパラメータにより擬似乱数系列とラテン 方陣行列を使用して上記擬似乱数置換行列を生成する擬似乱数置換行列生成部と 入力された上記所定の情報長、符号化率及び列の最大次数を用いて、上記擬似 乱数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のための 取りうる次数分布の組み合わせのパラメータを算出する次数分布最適化用パラメータ 算出部と、
算出された上記取りうる次数分布の組み合わせを拘束条件として、密度発展法によ り検査行列生成の次数分布を最適化する次数分布最適化部と、
最適化された上記検査行列生成の次数分布により、生成された上記擬似乱数置換 行列を配置して低密度パリティ検査符号用の検査行列を生成する検査行列生成部 とを備えた検査行列生成装置。
[2] 上記情報長'符号化率'列最大次数入力部は所定の情報長、複数の符号化率及 び無符号化を除く最大符号化率の列の最大次数を入力し、
上記擬似乱数置換行列用パラメータ算出部は、入力された上記所定の情報長、複 数の符号化率及び無符号化を除く最大符号化率の列の最大次数を用いて上記擬 似乱数置換行列のためのパラメータを算出し、
上記次数分布最適化用パラメータ算出部は入力された上記所定の情報長、複数 の符号化率及び無符号化を除く最大符号化率の列の最大次数を用いて、上記擬似 乱数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のための 取りうる次数分布の組み合わせのパラメータを算出することを特徴とする請求項 1記 載の検査行列生成装置。
[3] 上記擬似乱数置換行列用パラメータ算出部は、入力された上記所定の情報長、符 号ィ匕率及び列の最大次数を用いて、後述の第 1及び第 2の擬似乱数置換行列のた めのパラメータを算出し、
上記擬似乱数置換行列生成部は、算出された上記第 1及び第 2の擬似乱数置換 行列のためのパラメータにより擬似乱数系列とラテン方陣行列を使用して上記第 1の 擬似乱数置換行列を生成すると共に、上記第 1の擬似乱数置換行列をさらに行方向 に圧縮して上記第 2の擬似乱数置換行列を生成し、
上記次数分布最適化用パラメータ算出部は、入力された上記所定の情報長、符号 化率及び列の最大次数を用いて、上記第 1及び第 2の擬似乱数置換行列を使用し て構成可能な検査行列生成の次数分布の最適化のための取りうる次数分布の組み 合わせのパラメータを算出し、
上記検査行列生成部は、最適化された上記検査行列生成の次数分布により、生成 された上記第 1及び第 2の擬似乱数置換行列を配置して低密度パリティ検査符号用 の検査行列を生成することを特徴とする請求項 1記載の検査行列生成装置。
[4] 上記擬似乱数置換行列用パラメータ算出部は、入力された上記所定の情報長、符 号ィ匕率及び列の最大次数を用いて、後述の第 3の擬似乱数置換行列のためのパラ メータ及び後述の第 2の擬似乱数置換行列のためのパラメータを算出し、
上記擬似乱数置換行列生成部は、算出された上記第 3の擬似乱数置換行列のた めのノ メータ及び上記第 2の擬似乱数置換行列のためのパラメータにより、擬似乱 数系列とラテン方陣行列を使用して生成した第 1の疑似乱数置換行列をさらに行と 列の置換を行って拡張して上記第 3の擬似乱数置換行列を生成すると共に、上記第 1の擬似乱数置換行列をさらに行方向に圧縮して上記第 2の擬似乱数置換行列を生 成し、
上記次数分布最適化用パラメータ算出部は、入力された上記所定の情報長、符号 化率及び列の最大次数を用いて、上記第 3及び第 2の擬似乱数置換行列を使用し て構成可能な検査行列生成の次数分布の最適化のための取りうる次数分布の組み 合わせのパラメータを算出し、
上記次数分布最適化部は、算出された上記取りうる次数分布の組み合わせを拘束 条件として、密度発展法又は Multi - edge typeの密度発展法により検査行列生成 の次数分布を最適化し、
上記検査行列生成部は、最適化された上記検査行列生成の次数分布により、生成 された上記第 3及び第 2の擬似乱数置換行列を配置して低密度パリティ検査符号用 の検査行列を生成することを特徴とする請求項 1記載の検査行列生成装置。
[5] 上記検査行列生成部により上記擬似乱数置換行列を配置することで生成される検 查行列のサイクルは 6以上であることを特徴とする請求項 1記載の検査行列生成装 置。
[6] 所定の情報長、符号化率及び列の最大次数を用いて算出された擬似乱数置換行 列のためのパラメータを入力し、擬似乱数系列とラテン方陣行列を使用して上記擬 似乱数置換行列を生成する擬似乱数置換行列生成部と、
上記所定の情報長、符号化率及び列の最大次数を用いて算出された上記擬似乱 数置換行列を使用して構成可能な検査行列生成の次数分布の最適化のための取り うる次数分布の組み合わせを拘束条件として、密度発展法により最適化された上記 検査行列生成の次数分布を入力し、生成された上記擬似乱数置換行列を配置して 低密度パリティ検査符号用の検査行列を生成する検査行列生成部とを備えた検査 行列生成装置。
[7] 請求項 1記載の検査行列生成装置と、
上記検査行列生成装置により生成された検査行列を使用して生成行列を求めて符 号語を生成する符号化器とを備えた通信装置。
[8] 請求項 1記載の検査行列生成装置と、
上記検査行列生成装置により生成された検査行列を使用して sum— productアル ゴリズムによる繰り返し復号を実施して推定結果を算出する復号器とを備えた通信装 置。
[9] 請求項 6記載の検査行列生成装置と、
上記検査行列生成装置により生成された検査行列を使用して生成行列を求めて符 号語を生成する符号化器とを備えた通信装置。
[10] 請求項 6記載の検査行列生成装置と、 上記検査行列生成装置により生成された検査行列を使用して sum— productアル ゴリズムによる繰り返し復号を実施して推定結果を算出する復号器とを備えた通信装 置。
PCT/JP2004/014958 2004-10-08 2004-10-08 検査行列生成装置及び通信装置 WO2006040805A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US11/664,456 US20080184084A1 (en) 2004-10-08 2004-10-08 Check Matrix Generating Apparatus and Communication Apparatus
CNA2004800441456A CN101036300A (zh) 2004-10-08 2004-10-08 校验矩阵生成装置和通信装置
EP04792226A EP1806849A4 (en) 2004-10-08 2004-10-08 APPARATUS PROVIDING A VERIFICATION MATRIX AND COMMUNICATION APPARATUS
PCT/JP2004/014958 WO2006040805A1 (ja) 2004-10-08 2004-10-08 検査行列生成装置及び通信装置
JP2006540790A JP4772689B2 (ja) 2004-10-08 2004-10-08 検査行列生成方法及び通信方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2004/014958 WO2006040805A1 (ja) 2004-10-08 2004-10-08 検査行列生成装置及び通信装置

Publications (1)

Publication Number Publication Date
WO2006040805A1 true WO2006040805A1 (ja) 2006-04-20

Family

ID=36148107

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2004/014958 WO2006040805A1 (ja) 2004-10-08 2004-10-08 検査行列生成装置及び通信装置

Country Status (5)

Country Link
US (1) US20080184084A1 (ja)
EP (1) EP1806849A4 (ja)
JP (1) JP4772689B2 (ja)
CN (1) CN101036300A (ja)
WO (1) WO2006040805A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1868294A1 (en) * 2006-06-15 2007-12-19 Samsung Electronics Co., Ltd. Apparatus and method of encoding/decoding block low density parity check codes in a communication system
JP2008177654A (ja) * 2007-01-16 2008-07-31 Sony Corp データ処理装置、およびデータ処理方法、並びにコンピュータ・プログラム
EP2053751A1 (en) * 2006-08-04 2009-04-29 Mitsubishi Electric Corporation Inspection matrix generation method, encoding method, communication device, communication system, and encoder

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4743156B2 (ja) * 2007-04-27 2011-08-10 ソニー株式会社 復号装置
US8516352B2 (en) * 2009-07-21 2013-08-20 Ramot At Tel Aviv University Ltd. Compact decoding of punctured block codes
US8516351B2 (en) * 2009-07-21 2013-08-20 Ramot At Tel Aviv University Ltd. Compact decoding of punctured block codes
US8375278B2 (en) * 2009-07-21 2013-02-12 Ramot At Tel Aviv University Ltd. Compact decoding of punctured block codes
US9397699B2 (en) * 2009-07-21 2016-07-19 Ramot At Tel Aviv University Ltd. Compact decoding of punctured codes
CN102025441B (zh) * 2009-09-11 2013-07-31 北京泰美世纪科技有限公司 Ldpc码校验矩阵的构造方法、ldpc码的编码方法和编码装置
US9374184B2 (en) * 2012-03-23 2016-06-21 Nokia Solutions And Networks Oy Controlling of code block to physical layer mapping
US9424307B2 (en) * 2012-10-11 2016-08-23 Scott E. Lilienthal Multivariate data analysis method
US20150046765A1 (en) 2013-02-08 2015-02-12 Sony Corporation Data processing apparatus and data processing method
CN104393959A (zh) * 2014-11-28 2015-03-04 电子科技大学 一种编码调制方法
JP6911540B2 (ja) * 2017-05-31 2021-07-28 ソニーグループ株式会社 送信装置、送信方法、受信装置、及び、受信方法
CN112558924B (zh) * 2020-12-11 2022-08-26 黑龙江大学 伪随机二值序列的局部弱随机现象的快速检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003198383A (ja) * 2001-12-27 2003-07-11 Mitsubishi Electric Corp Ｌｄｐｃ符号用検査行列生成方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7340671B2 (en) * 2003-10-10 2008-03-04 Regents Of The University Of California Decoding low density parity codes
KR100922956B1 (ko) * 2003-10-14 2009-10-22 삼성전자주식회사 저밀도 패리티 검사 코드의 부호화 방법
US20060242530A1 (en) * 2005-03-31 2006-10-26 Nec Laboratories America, Inc. Method for constructing finite-length low density parity check codes

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003198383A (ja) * 2001-12-27 2003-07-11 Mitsubishi Electric Corp Ｌｄｐｃ符号用検査行列生成方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHUNG S.-Y. ET AL: "Analysis of Sum-Product Decoding of Low-Density Parity-Check Codes Using a Gaussian Approximation", IEEE TRANSACTIONS ON INFORMATION TEHEORY, vol. 47, no. 2, February 2001 (2001-02-01), pages 657 - 670, XP002969535 *
MATSUMOTO W. ET AL: "LDPC Fugo Koseiho no Kiso to Oyo", THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS 2003 NEN KISO KYOKAI SOCIETY TAIKAI KOEN RONBUNSHU, 23 September 2003 (2003-09-23), pages SS17 - SS18, XP002998474 *
RICHARDSON T.J. ET AL: "The Capacity of Low-Density Parity-Check Codes Under Message-Passing Decoding", IEEE TRANSACTIONS ON INFORMATION THEORY, vol. 47, no. 2, February 2001 (2001-02-01), pages 599 - 618, XP002228426 *
See also references of EP1806849A4 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1868294A1 (en) * 2006-06-15 2007-12-19 Samsung Electronics Co., Ltd. Apparatus and method of encoding/decoding block low density parity check codes in a communication system
US8020062B2 (en) 2006-06-15 2011-09-13 Samsung Electronics Co., Ltd. Apparatus and method of encoding/decoding block low density parity check codes in a communication system
EP2053751A1 (en) * 2006-08-04 2009-04-29 Mitsubishi Electric Corporation Inspection matrix generation method, encoding method, communication device, communication system, and encoder
EP2053751A4 (en) * 2006-08-04 2011-06-22 Mitsubishi Electric Corp INVESTIGATION MATRIX GENERATION PROCESS, CODING METHOD, COMMUNICATION DEVICE, COMMUNICATION SYSTEM AND CODIER
JP2008177654A (ja) * 2007-01-16 2008-07-31 Sony Corp データ処理装置、およびデータ処理方法、並びにコンピュータ・プログラム

Also Published As

Publication number Publication date
JP4772689B2 (ja) 2011-09-14
CN101036300A (zh) 2007-09-12
US20080184084A1 (en) 2008-07-31
JPWO2006040805A1 (ja) 2008-05-15
EP1806849A1 (en) 2007-07-11
EP1806849A4 (en) 2008-03-05

Similar Documents

Publication Publication Date Title
CN1866751B (zh) 一种低密度奇偶校验码的构造方法及装置
US8185797B2 (en) Basic matrix, coder/encoder and generation method of the low density parity check codes
EP3457575B1 (en) Encoding method and device and decoding method and device for structured ldpc
Pishro-Nik et al. Results on punctured low-density parity-check codes and improved iterative decoding techniques
JP4598085B2 (ja) 検査行列生成方法
Sadeghi et al. Low-density parity-check lattices: Construction and decoding analysis
US8966336B2 (en) Selective merge and partial reuse LDPC (low density parity check) code construction for limited number of layers belief propagation (BP) decoding
WO2006075382A1 (ja) 符号化方法、復号方法及びそれらの装置
JP5436688B2 (ja) 線形ブロック符号を使用する通信システムにおけるパリティ検査行列を生成する方法及び装置とそれを用いる送受信装置及び方法
WO2006040805A1 (ja) 検査行列生成装置及び通信装置
CN102394659B (zh) Ldpc码校验矩阵构造方法及对应矩阵乘法运算装置
JPWO2009060627A1 (ja) 符号化方法および送信装置
JP2004266463A (ja) 検査行列生成方法および検査行列生成装置
WO2008016117A1 (fr) Procédé de génération de matrice d&#39;inspection, procédé d&#39;encodage, dispositif et système de communication et encodeur
CN101373976A (zh) 生成ldpc校验矩阵的方法和设备
US20100050045A1 (en) Method of generating a parity check matrix for ldpc encoding and decoding
Shibata et al. Design of irregular LDPC codes without markers for insertion/deletion channels
CN105871385B (zh) 一种ldpc卷积码构造方法
Wu et al. Design of high-rate LDPC codes based on matroid theory
Malema Low-density parity-check codes: construction and implementation.
KR101145673B1 (ko) 특정 구조를 가지는 ｌｔ코드를 이용한 부호화 방법 및 이를 기반으로 하는 랩터 코드를 이용한 부호화 방법
Hong et al. Optimal puncturing of block-type LDPC codes and their fast convergence decoding
Hussein et al. Comparisons of soft decision decoding algorithms based LDPC wireless communication system
Zikry et al. Performance Analysis of LDPC Decoding Techniques
Baik et al. Irregular low-density parity-check lattices

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 11664456

Country of ref document: US

Ref document number: 1020077007540

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 200480044145.6

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: 2006540790

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2004792226

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2004792226

Country of ref document: EP