JP3808769B2 - Ｌｄｐｃ符号用検査行列生成方法 - Google Patents
Ｌｄｐｃ符号用検査行列生成方法 Download PDFInfo
- Publication number
- JP3808769B2 JP3808769B2 JP2001397922A JP2001397922A JP3808769B2 JP 3808769 B2 JP3808769 B2 JP 3808769B2 JP 2001397922 A JP2001397922 A JP 2001397922A JP 2001397922 A JP2001397922 A JP 2001397922A JP 3808769 B2 JP3808769 B2 JP 3808769B2
- Authority
- JP
- Japan
- Prior art keywords
- column
- code
- check matrix
- row
- ldpc code
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000011159 matrix material Substances 0.000 title claims description 83
- 238000000034 method Methods 0.000 title claims description 66
- 238000009826 distribution Methods 0.000 claims description 27
- 238000012545 processing Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 21
- 230000000694 effects Effects 0.000 description 11
- 230000005540 biological transmission Effects 0.000 description 8
- 238000009795 derivation Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000001788 irregular Effects 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
Images
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/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/1177—Regular LDPC codes with parity-check matrices wherein all rows and columns have the same row weight and column weight, respectively
-
- 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/1151—Algebraically constructed LDPC codes, e.g. LDPC codes derived from Euclidean geometries [EG-LDPC codes]
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Error Detection And Correction (AREA)
- Detection And Correction Of Errors (AREA)
- Detection And Prevention Of Errors In Transmission (AREA)
- Complex Calculations (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、誤り訂正符号としてＬＤＰＣ（Low-Density Parity-Check）符号を採用した符号化器におけるＬＤＰＣ符号用検査行列生成方法に関するものである。
【０００２】
【従来の技術】
図１３は、ＬＤＰＣ符号化／復号システムを示す図である。図１３において、１０１は符号化器であり、１０２は変調器であり、１０３は通信路であり、１０４は復調器であり、１０５は復号器である。ここでは、従来のＬＤＰＣ符号用検査行列生成方法を説明する前に、ＬＤＰＣ符号を使用した場合の符号化，復号の流れについて説明する。
【０００３】
まず、送信側の符号化器１０１では、後述する所定の方法で検査行列Ｈを生成する。そして、以下の条件に基づいて生成行列Ｇを求める。
Ｇ：ｋ×ｎ行列（ｋ：情報長，ｎ：符号語長）
ＧＨT＝０（Ｔは転置行列）
【０００４】
その後、符号化器１０１では、情報長ｋのメッセージ（ｍ1ｍ2…ｍk）を受け取り、上記生成行列Ｇを用いて符号語Ｃを生成する。
Ｃ＝（ｍ1ｍ2…ｍk）Ｇ
＝（ｃ1ｃ2…ｃn） （ただし、Ｈ（ｃ1ｃ2…ｃn）T＝０）
【０００５】
そして、変調器１０２では、生成した符号語Ｃに対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル変調を行い、送信する。
【０００６】
一方、受信側では、復調器１０４が、通信路１０３を介して受け取った変調信号に対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル復調を行い、さらに、復号器１０５が、ＬＤＰＣ符号化された復調結果に対して、「ｓｕｍ−ｐｒｏｄｕｃｔアルゴリズム」によるくり返し復号を実施し、推定結果（もとのｍ1ｍ2…ｍkに対応）を出力する。
【０００７】
以下、従来のＬＤＰＣ符号用検査行列生成方法について説明する。ＬＤＰＣ符号用の検査行列としては、たとえば、ＬＤＰＣの発案者Ｇａｌｌａｇｅｒにより以下のような行列が提案されている（図１４参照）。
【０００８】
図１４に示す行列は、「１」と「０」の２値の行列で、「１」の部分を塗りつぶしている。他の部分は全て「０」である。この行列は、１行の「１」の数（これを行の重みと表現する）が４で、１列の「１」の数（これを列の重みと表現する）が３であり、全ての列と行の重みが均一なため、これを一般に「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」と呼んでいる。また、Ｇａｌｌａｇｅｒの符号では、たとえば、図１４に示すように、行列を３ブロックに分け、２ブロック目と３ブロック目に対してランダム置換を行っている。
【０００９】
しかしながら、このランダム置換には、所定のルールが無いため、より特性の良好な符号を見つけるためには、計算機による時間のかかる探索を行わなければならなかった。
【００１０】
そこで、たとえば、計算機探索によらなくても確定的に行列を生成でき、比較的安定した良好な特性を示すＬＤＰＣ符号として、ユークリット幾何符号を用いる方法が、Ｙ．Ｋｏｕ等（Y. Kou, S. Lin, and M. P. C. Fossorier, "Low Density Parity Check Codes Based on Finite Geometries: A Rediscovery," ISIT 2000, pp. 200, Sorrento, Itary, June 25-30, 2000.）によって提案された。この方法では、規則的なｅｎｓｅｍｂｌｅ（アンサンブル）で構成された「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」について説明されている。
【００１１】
ここでは、有限幾何符号の一種であるユークリッド幾何符号ＥＧ（２，２6）を用いてＬＤＰＣ符号の検査行列を生成する方法について提案され、誤り率１０-4点において、シャノン限界から１．４５ｄＢに接近した特性を得ている。図１５は、たとえば、ユークリッド幾何符号ＥＧ（２，２2）の構成を示す図であり、行，列のそれぞれの重みが４，４の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」構造をしている。
【００１２】
したがって、ユークリッド幾何符号ＥＧ（ｍ，２s）の場合、その特性は、以下のように規定される。
符号長： ｎ＝２2s−１
冗長ビット長： ｎ−ｋ＝３s−１
情報長： ｋ＝２2s−３s
最小距離： ｄmin＝２s＋１
密度： ｒ＝２s／（２2s−１）
【００１３】
図１５を見ても分かるように、ユークリッド幾何符号は、各行の「１」の配置が行毎に巡回シフトした構造になっており、符号が容易にかつ確定的に構成できる特長がある。
【００１４】
Ｙ．Ｋｏｕらによる検査行列の生成方法では、さらに、上記ユークリット幾何符号に基づいて行と列の重みを変更し、行，列を必要に応じて拡張している。たとえば、ＥＧ（２，２2）の列の重みを１／２に分離する場合、Ｙ．Ｋｏｕらの論文では、１列内に４つある重みを１つ置きに２個づつ分離する。図１６は、列の重みを４から２に規則的に分離した例を示す図である。
【００１５】
一方、上記「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」の特性よりも「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の特性の方が良好であることが、Ｌｕｄｙ等（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.）により報告された。なお、上記「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」は、列と行の重みがそれぞれあるいはどちらか一方が均一でないＬＤＰＣ符号を表す。
【００１６】
そして、それは、Ｒｉｃｈａｒｄｓｏｎ等（T. J. Richardson and R. Urbanke, "The capacity of low-density parity-check codes under message-passing decoding," IEEE Trans. Inform. Theory, vol.47,No.2, pp.599-618, Feb. 2001.）、あるいはＣｈｕｎｇ等（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.）によって理論的に解析された。
【００１７】
特に、Ｃｈｕｎｇ等は、繰り返し復号器における入力と出力の対数尤度比（ＬＬＲ）がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析し、良好な行と列の重みのアンサンブルを求めている。
【００１８】
【発明が解決しようとする課題】
しかしながら、たとえば、上記Ｃｈｕｎｇ等による従来のＬＤＰＣ符号用検査行列生成方法は、行内の「１」の点の数（後述するバリアブルノードの次数配分に相当）と、列内の「１」の点の数（後述するチェックノードの次数配分に相当）と、の両方を変数として、下記の（１）式（ｒａｔｅ：符号化率）が最大となるバリアブルノードの次数配分およびチェックノードの次数配分を求めている。すなわち、ＳＮＲ（Signal to Noise Ratio）が最小となるアンサンブルを線形計画法により探索している。
【００１９】
【数１】
【００２０】
そのため、上記「ｒａｔｅ」の最大値により得られる検査行列が流動的になり、特性が安定しない、という問題があった。また、従来のＬＤＰＣ符号用検査行列生成方法は、バリアブルノードの次数配分の導出とチェックノードの次数配分の導出とを所定回数にわたって繰り返し行っているため、探索処理にある程度の時間を要する、という問題もあった。
【００２１】
本発明は、上記に鑑みてなされたものであって、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を、短時間で容易に探索可能なＬＤＰＣ符号用検査行列生成方法を得ることを目的とする。
【００２２】
【課題を解決するための手段】
上述した課題を解決し、目的を達成するために、本発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、復号器における入出力データの対数尤度比がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析することによって、誤りが０となるＳＮＲの限界（threshold）を求めることとし、さらに、符号化レートを固定した状態で、かつガウスノイズが最大になるように、行の重みと列の重みの最適なアンサンブル（thresholdが最小となるアンサンブル）を１回の線形計画法で探索し、当該アンサンブルにしたがってＬＤＰＣ符号用の検査行列を生成することを特徴とする。
【００２３】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、前記アンサンブルを探索後、当該探索結果に基づいてユークリット幾何符号の各行または各列からランダムに「１」を抽出し、各行または各列を分割することによって、Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号の検査行列を生成することを特徴とする。
【００２４】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、前記アンサンブルの重み配分を、重み単位の重み総数が整数で、かつ重み単位の重み総数の総和とユークリット幾何符号の「１」の総数とが等しくなるように調整し、調整後のアンサンブルに基づいて分割処理を行うことを特徴とする。
【００２５】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を作成し、ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする。
【００２６】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を作成し、ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする。
【００２７】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を複数個作成し、列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする。
【００２８】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を複数個作成し、列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする。
【００２９】
【発明の実施の形態】
以下に、本発明にかかるＬＤＰＣ符号用検査行列生成方法の実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
【００３０】
実施の形態１．
本実施の形態のＬＤＰＣ符号用検査行列生成方法を説明する前に、本実施の形態のＬＤＰＣ符号用検査行列生成方法を実現可能な符号化器の位置付け、および「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の従来の検査行列生成方法について説明する。なお、ＬＤＰＣ符号化／復号システムの構成については、先に説明した図１３と同様である。
【００３１】
送信側の符号化器１０１では、後述する本実施の形態のＬＤＰＣ符号用検査行列生成方法で検査行列Ｈを生成する。そして、以下の条件に基づいて生成行列Ｇを求める。
Ｇ：ｋ×ｎ行列（ｋ：情報長，ｎ：符号語長）
ＧＨT＝０（Ｔは転置行列）
【００３２】
その後、符号化器１０１では、情報長ｋのメッセージ（ｍ1ｍ2…ｍk）を受け取り、上記生成行列Ｇを用いて符号語Ｃを生成する。
Ｃ＝（ｍ1ｍ2…ｍk）Ｇ
＝（ｃ1ｃ2…ｃn） （ただし、Ｈ（ｃ1ｃ2…ｃn）T＝０）
【００３３】
そして、変調器１０２では、生成した符号語Ｃに対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル変調を行い、送信する。
【００３４】
一方、受信側では、復調器１０４が、通信路１０３を介して受け取った変調信号に対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル復調を行い、さらに、復号器１０５が、ＬＤＰＣ符号化された復調結果に対して、「ｓｕｍ−ｐｒｏｄｕｃｔアルゴリズム」による繰り返し復号を実施し、推定結果（もとのｍ1ｍ2…ｍkに対応）を出力する。
【００３５】
つぎに、Ｃｈｕｎｇ等（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.）によって理論的に解析された、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の従来の検査行列生成方法について詳細に説明する。ここでは、繰り返し復号器における入力と出力の対数尤度比（ＬＬＲ）がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析し、良好な行と列の重みのアンサンブルを求めている。
【００３６】
なお、上記論文に記述されたＬＤＰＣ符号用検査行列生成方法であるガウス近似法（Gaussian Approximation）では、前提として、検査行列における行内の「１」の点をバリアブルノードと定義し、列内の「１」の点をチェックノードと定義する。
【００３７】
まず、チェックノードからバリアブルノードへのＬＬＲメッセージ伝搬を解析する。０＜ｓ＜∞と０≦ｔ＜∞という条件において、以下の関数（２）式を定義する。なお、ｓ＝ｍu0はｕ０の平均値であり、ｕ０は分散値σn 2のガウスノイズを含む伝送路を経由して受信した信号の対数尤度比（ＬＬＲ）であり、ｔは所定の繰り返しの時点におけるチェックノードのＬＬＲ出力値のアンサンブル平均である。
【００３８】
【数２】
【００３９】
なお、上記λ（ｘ）およびρ（ｘ）は、それぞれバリアブルノードおよびチェックノードの次数配分（バリアブルノードとチェックノードの各１行，各１列内の「１」の数を次数と表現する）の生成関数を表し、（３）式および（４）式のように表すことができる。また、λiとρiは、それぞれ次数ｉのバリアブルノードとチェックノードに属するエッジの比率を表す。また、ｄlは最大バリアブルノードの次数であり、ｄrは最大チェックノードの次数である。
【００４０】
【数３】
【００４１】
【数４】
【００４２】
ただし、φ（ｘ）は下記（５）式のように定義する。
【００４３】
【数５】
【００４４】
そして、（２）式は、等価的に下記（６）式と表すことができる。
【００４５】
【数６】
【００４６】
なお、ｔlはｌ番目の繰り返し時点におけるチェックノードのＬＬＲ出力値のアンサンブル平均である。
【００４７】
ここで、誤りが０となりうるＳＮＲの限界（threshold）を求めるための条件は、ｌ→∞のときにｔl（ｓ）→∞（Ｒ+と表現する）となることであり、この条件を満たすためには、以下の条件（７）式を満たす必要がある。
【００４８】
【数７】
【００４９】
つぎに、バリアブルノードからチェックノードへのＬＬＲメッセージ伝搬を解析する。０＜ｓ＜∞と０＜ｒ≦１という条件において、以下の関数（８）式を定義する。なお、ｒの初期値ｒ0はφ（ｓ）である。
【００５０】
【数８】
【００５１】
そして、（８）式は、等価的に下記（９）式と表すことができる。
【００５２】
【数９】
【００５３】
ここで、誤りが０となりうるＳＮＲの限界（threshold）を求めるための条件は、ｒl（ｓ）→０となることであり、この条件を満たすためには、以下の条件（１０）式を満たす必要がある。
【００５４】
【数１０】
【００５５】
さらに、上記Ｃｈｕｎｇ等の論文では、上記式を用いて以下の手順でバリアブルノードとチェックノードの最適な次数を探索している。
（１）生成関数λ（ｘ）とガウスノイズσnが与えられていると仮定し、生成関数ρ（ｘ）を変数として、前述した（１）式が最大となる点を探索する。なお、この探索における拘束条件は、ρ（１）＝１と正規化することと、上記（７）式を満たすことである。
（２）生成関数ρ（ｘ）とガウスノイズσnが与えられていると仮定し（たとえば、（１）の結果より得られる値）、生成関数λ（ｘ）を変数として、（１）式が最大となる点を探索する。なお、この探索における拘束条件は、λ（１）＝１と正規化することと、上記（１０）式を満たすことである。
（３）最大「ｒａｔｅ」を求めるために、上記（１）と上記（２）を繰り返し実行し、生成関数λ（ｘ）と生成関数ρ（ｘ）のより良好なアンサンブルを線形計画法で探索する。
（４）最後に、ガウスノイズσnより信号電力を１と正規化して、ＳＮＲの限界（threshold）を求める。
【００５６】
【数１１】
【００５７】
しかしながら、上記Ｃｈｕｎｇ等の論文では、「ｒａｔｅ（符号化率）」の最大値により得られる検査行列が流動的になり、設計において仕様として固定されるｒａｔｅが変動してしまうため、実設計にむかない、という問題があった。また、上記ガウス近似法では、バリアブルノードの次数配分の導出とチェックノードの次数配分の導出とを所定回数にわたって繰り返し行っているため、探索処理にある程度の時間を要する、という問題もあった。
【００５８】
そこで、本実施の形態においては、確定的でかつ特性が安定した「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の検査行列を、短時間で容易に探索する。図１は、実施の形態１のＬＤＰＣ符号用検査行列生成方法を示すフローチャートである。
（１）「ｒａｔｅ」が与えられているものと仮定する。すなわち、要求「ｒａｔｅ」を固定する。実際の設計では、目標「ｒａｔｅ」が予め指定されている場合が多いためである。
（２）生成関数λ（ｘ）と生成関数ρ（ｘ）を同時に変数として扱い、ガウスノイズσnが最大になるように、線形計画法で最適な生成関数λ（ｘ）と生成関数ρ（ｘ）を探索する。この探索における拘束条件は、λ（１）＝１，ρ（１）＝１と正規化し、さらに上記（１０）式を満たすことである。
【００５９】
このように、本実施の形態では、上記（９）式と上記（１０）式を満たす生成関数λ（ｘ）と生成関数ρ（ｘ）を１回の線形計画法で求めることとしたため、上記論文のように、生成関数λ（ｘ）の導出と生成関数ρ（ｘ）の導出を繰り返し実行し、双方の最適値を求める方法よりも、容易かつ短時間に、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を生成することができる。
【００６０】
実施の形態２．
実施の形態２では、前述の実施の形態１においてユークリッド幾何符号を採用し、１行あるいは１列の「１」の配置を分割することによって、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の検査行列を生成する。
【００６１】
まず、実施の形態１におけるＬＤＰＣ符号用検査行列生成方法により、生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを導出する。図２は、ｒａｔｅ＝０．５とした場合の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。なお、σGAはガウス近似法により導出した「threshold」時のノイズ分散値を表し、ＳＮＲnorm（ＧＡ）はガウス近似法により導出した「threshold」のＳＮＲとシャノン限界のＳＮＲとの差分を表し、ｘは重みを表し、λxおよびρxはそれぞれバリアブルノードとチェックノードの重み配分を表す。
【００６２】
また、基準になるユークリッド幾何符号は、ＥＧ（２，２5）を想定し、ｄl＝３２とした。また、重み配分λxのｘの値および重み配分ρxのｘの値は、組み合わせで３２（ｄl）を構成できる値とする。
【００６３】
図３は、ＥＧ（２，２5）を想定し、ｄl＝３２とした場合の、分割テーブルを示す図である。図３に示すとおり、図２のｘは、組み合わせにより必ず３２となる。たとえば、図示の７ｘ４と２ｘ２の組み合わせは、重みが３２の１列を、重みが７の４列と重みが２の２列に分割可能なことを表している。このように、ＥＧ（２，２5）の符号を基本として、図３のように重みが３２の各行列を適切に分割すると、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の検査行列を構成することができる。
【００６４】
ここで、分割処理を行う前に、図２に示す生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルの重み配分を以下の手順で調整する。図４は、重み配分調整用テーブルを示す図である。なお、ＥＧ（２，２5）のユークリッド幾何符号は、１０２３行×１０２３列で構成される。
【００６５】
（１）ガウス近似法で求めた生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブル（表１参照）をテーブルの２列目と３列目に設定する。
（２）重み配分λxおよびρx（３列目）と、ＥＧ（２，２5）における全行列の「１」の総数ＴＰ＝３２７３６と、を乗算し、重み単位の重み総数を求め、さらに、当該重み単位の重み総数とその総和を４列目に設定する。
（３）重み単位の重み総数（４列目）を対応する重みｘで割り、重み単位の総列数を求め、それを５列目に設定する。
（４）重み単位の総列数が小数点以下を含む場合、丸め処理（四捨五入，切上げ，切捨て等）を行い、その結果を６列目に設定する。
（５）丸め処理後の重み単位の総列数（６列目）と対応する重みｘとを乗算し、丸め処理後の重み単位の重み総数を求め、それを７列目に設定する。そして、各重み総数の総和（７列目の合計の行）が行列内の「１」の総数（ＴＰ＝３２７３６）と等しいかどうかを確認する。
（６）行列内の「１」の総数に等しくない場合、丸め処理後の重み単位の重み総数（７列目）を整数単位で調整し、その結果を８列目に設定する。この場合、８列目の総和が、行列内の「１」の総数（ＴＰ＝３２７３６）に等しくなるように調整する。
（７）調整後の重み単位の重み総数（８列目）を対応する重みｘで割り、調整後の重み単位の総列数を求め、それを９列目に設定する。調整後の各重みの配分（１１列目）は、可能な限りガウス近似法で求めた値（３列目）に近い値にする。
【００６６】
図５は、重み配分後の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【００６７】
つぎに、ユークリッド幾何符号における１行あるいは１列の分割手順について説明する。
【００６８】
たとえば、分割手順に関して、Ｙ．Ｋｏｕ等の論文では、規則的に分割する方法を提示している。図６は、上記論文における分割手順を示す図である。まず、図６に示すように行列のナンバリングを行う。ここでは、列番号を左端から順に１，２，３，…とし、行番号を上から順に１，２，３，…とする。そして、たとえば、３２点×１列を８点×４列に分割する場合、下記（１２）式にしたがって規則的に分割する。
【００６９】
【数１２】
【００７０】
なお、ｍ＝１，２，３，４とし、ｎ＝０，１，２，３，４，５，６，７とし、ｌはＥＧ（２，２5）の列番号を表す。また、Ｂl（ｘ）はＥＧ（２，２5）のｌ列目の「１」の位置を表し、Ｓm（ｎ）は分割後の行列のｍ列目の「１」の位置を表す。
【００７１】
具体的にいうと、ＥＧ（２，２5）における１列中の「１」の位置を示す行番号は、
Ｂl（ｘ）＝｛1 32 114 136 149 223 260 382 402 438 467 507 574 579 588 622 634 637 638 676 717 728 790 851 861 879 947 954 971 977 979 998｝
となり、その結果、分割後の行列における１〜４列目の「１」の位置を示す行番号は、Ｂl（ｘ）から「１」の番号が規則的に抽出され、
Ｓ1（ｎ）＝｛1 149 402 574 634 717 861 971｝
Ｓ2（ｎ）＝｛32 223 438 579 637 728 879 977｝
Ｓ3（ｎ）＝｛114 260 467 588 638 790 947 979｝
Ｓ4（ｎ）＝｛136 382 507 622 676 851 954 998｝
となる。すなわち、３２点×１列が８点×４列に分割される。
【００７２】
一方、本実施の形態におけるユークリッド幾何符号の分割処理は、上記のように規則的に分割するのではなく、Ｂl（ｘ）から「１」の番号をランダムに抽出する。なお、この抽出処理は、ランダム性が保持されるのであればどのような方法を用いてもよい。
【００７３】
したがって、分割後の行列のｍ列目の「１」の位置の一例をＲm（ｎ）とした場合、Ｒm（ｎ）は、
Ｒ1（ｎ）＝｛1 114 574 637 851 879 977 979｝
Ｒ2（ｎ）＝｛32 136 402 467 588 728 861 971｝
Ｒ3（ｎ）＝｛149 260 382 438 579 638 717 998｝
Ｒ4（ｎ）＝｛223 507 622 634 676 790 947 954｝
となる。
【００７４】
上記のような本実施の形態の分割手順をグラフ上で表現すると、以下のように表現することができる。図７は、分割前のＥＧ（２，２5）のグラフを示す図である。なお、両ノードを結ぶ線はエッジと表現する。図７では、分割前の１０２３行×１０２３列（各行列の重みがそれぞれ３２）のユークリット幾何符号を表現している。また、図８は、ＥＧ（２，２5）のエッジをランダムに選択した、分割後のグラフを示す図である。
【００７５】
ここで、上記で説明したＬＤＰＣ符号の特性を比較する。図９は、Ｅｂ／Ｎｏ（情報１ビットあたりの信号電力対ノイズ電力比）と誤り率特性（ＢＥＲ）との関係を示す図である。なお、繰り返し回数は５０回で、復号法は「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」である。
【００７６】
なお、図中"Simple regular extended EG(2,25)"は、Ｙ．Ｋｏｕ等の発案によるＥＧ（２，２5）の規則的な列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」であり、"Random regular extended EG(2,25)"は、本実施の形態によるＥＧ（２，２5）のランダムな列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」である。図１０は、これらの「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【００７７】
また、図中"Simple irregular extended EG(2,25)"は、実施の形態１の方法によって特定されたアンサンブルに対して、Ｙ．Ｋｏｕ等の発案によるＥＧ（２，２5）の規則的な列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であり、"Random irregular extended EG(2,25)"は、実施の形態１の方法によって特定されたアンサンブルに対して、本実施の形態によるＥＧ（２，２5）のランダムな列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」である。図１１は、これらの「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【００７８】
図９からわかるように、同一レートでは、「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」より「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のほうが性能がよい。また、Ｙ．Ｋｏｕ等の論文のような規則的な分割では、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であっても大幅な改善は見込めないが、本実施の形態のランダムな分割を実施すると性能が画期的に改善される。
【００７９】
つぎに、上記ランダム分割の一例を詳細に説明する。ここでは、ランダム分割を行う場合のランダム系列を容易かつ確定的に生成する。この方法による利点は、送信側と受信側が同じランダム系列を生成できることにある。これは、現実のシステムではきわめて重要となる。また、符号特性の条件が正確に規定できる、という利点もある。
【００８０】
（１）基本のランダム系列を作成する。
以下に、ランダム系列作成の一例を記述する。ここでは、説明の便宜上、ユークリッド幾何符号ＥＧ（２，２4）を用いる。ユークリッド幾何符号ＥＧ（２，２4）の場合、１行に存在する「１」の数は２4個である。
【００８１】
ＰをＰ≧２sを満たす最小の素数とした場合、たとえば、２4のときはＰ＝１７となる。ここで、系列長Ｐ−１＝１６の基本のランダム系列Ｃ（ｉ）を（１３）式にしたがって作成する。
Ｃ（１）＝１
Ｃ（ｉ＋１）＝Ｇ0×Ｃ（ｉ） ｍｏｄ Ｐ …（１３）
ただし、ｉ＝１，…，Ｐ−１とし、Ｇ0はガロア体ＧＦ（Ｐ）の原始元である。その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛1 3 9 10 13 5 15 11 16 14 8 7 4 12 2 6｝
となる。なお、２4の場合は、Ｐが２4＋１となり、ランダム系列長が１６になるので問題ないが、２5の場合等は、Ｐが２5＋１以上になりランダム系列長が２5を超えてしまう。そのような場合は、２5を超える数値をランダム系列から削除することで対応する。
【００８２】
そして、ランダム系列Ｃ（ｉ）を巡回シフトして、ランダム系列のラテン方陣を生成する。図１２は、ランダム系列のラテン方陣を示す図である。
【００８３】
（２）ランダム系列のラテン方陣を用いて、ユークリッド幾何符号ＥＧ（２，２4）の１列目の「１」の位置を示す行番号の配列Ｂl（ｘ）から、「１」の番号をランダムに抽出し、ユークリッド幾何符号の第１列目を分割する。なお、上記配列Ｂl（ｘ）は、
Ｂl（ｘ）＝｛1 14 16 19 45 49 55 107 115 126 127 182 210 224 231 247｝
である。
【００８４】
たとえば、ユークリッド幾何符号の第１列目を４分割する場合（各列の重みが４となる）、上記ランダム系列のラテン方陣の１行目から、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ1（ｎ）＝｛1,3,9,10｝
Ｌ2（ｎ）＝｛13,5,15,11｝
Ｌ3（ｎ）＝｛16,14,8,7｝
Ｌ4（ｎ）＝｛4,12,2,6｝
【００８５】
そして、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を用いて、配列Ｂl（ｘ）から、「１」の番号をランダムに抽出する。その結果、Ｒ（Ｌ1（ｎ））〜Ｒ（Ｌ4（ｎ））は、
Ｒ（Ｌ1（ｎ））＝｛1,16,115,126｝
Ｒ（Ｌ2（ｎ））＝｛210,45,231,127｝
Ｒ（Ｌ3（ｎ））＝｛247,224,107,55｝
Ｒ（Ｌ4（ｎ））＝｛19,182,14,49｝
となる。
【００８６】
（３）同様に、ランダム系列のラテン方陣を用いて、ユークリッド幾何符号ＥＧ（２，２4）の２列目の「１」の位置を示す行番号の配列Ｂ2（ｘ）から、「１」の番号をランダムに抽出し、ユークリッド幾何符号の第２列目を分割する。なお、上記配列Ｂ2（ｘ）は、
Ｂ2（ｘ）＝｛2 15 17 20 46 50 56 108 116 127 128 183 211 225 232 248｝
である。
【００８７】
たとえば、ユークリッド幾何符号の第２列目を４分割する場合（各列の重みが４となる）、上記ランダム系列のラテン方陣の２行目から、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ5（ｎ）＝｛3,9,10,13｝
Ｌ6（ｎ）＝｛5,15,11,16｝
Ｌ7（ｎ）＝｛14,8,7,4｝
Ｌ8（ｎ）＝｛12,2,6,1｝
【００８８】
そして、４つのランダム系列Ｌ5（ｎ）〜Ｌ8（ｎ）を用いて、配列Ｂ2（ｘ）から、「１」の番号をランダムに抽出する。その結果、Ｒ（Ｌ5（ｎ））〜Ｒ（Ｌ8（ｎ））は、
Ｒ（Ｌ5（ｎ））＝｛17,116,127,211｝
Ｒ（Ｌ6（ｎ））＝｛46,232,128,248｝
Ｒ（Ｌ7（ｎ））＝｛225,108,56,20｝
Ｒ（Ｌ8（ｎ））＝｛183,15,50,2｝
となる。以降、同様の手順でユークリッド幾何符号の全ての列を分割する。
【００８９】
（４）なお、上記ランダム系列のラテン方陣でランダム系列が不足するような場合は、（１３）式の基本のランダム系列のＧ0に原始元より大きくかつ２s以下の素数を代入して再度基本のランダム系列を作り、同様の手順で分割する。
【００９０】
つぎに、もう一つのランダム系列のラテン方陣の作成法について説明する。ここでは、説明の便宜上、ユークリッド幾何符号ＥＧ（２，２5）を用いる。ユークリッド幾何符号ＥＧ（２，２5）の場合、１行に存在する「１」の数は２5個である。
【００９１】
上記例では、ＥＧ（２，２4）のため１６（行数）×１６（列数）のラテン方陣を作成したが、これでは不足する場合がある。ここでは、ＥＧ（２，２5）を例にとり、３２（行数）×９６０（列数）まで拡張する。行列のサイズは、３２（行数）×３２（列数）のラテン方陣の組を３０組列方向に並べたことにより決められる。
【００９２】
また、上記の例では、ＰをＰ≧２sを満たす最小の素数としたが、ここでは、ＰをＰ≧２sを満たす最大の素数とする。たとえば、２5のときはＰ＝３１となる。
【００９３】
ここで、系列長Ｐ＝３１の基本のランダム系列Ｃ（ｉ）を（１４）式に従って作成する。
Ｃ（１）＝０
Ｃ（ｉ＋１）＝Ｇo×Ｃ（ｉ） ｍｏｄ Ｐ …（１４）
ただし、ｉ＝１，…，Ｐ−１とし、Ｇoはガロア体ＧＦ（Ｐ）の原始元である。その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛0,1,3,9,27,19,26,17,20,29,25,13,8,24,10,30,28,22,4,12,5,15,14,11,2,6,18,23,7,21｝
となる。
【００９４】
つぎに以下の操作をする。
Ｃ（ｉ）＝Ｃ（ｉ）＋１， Ｃ（３２）＝３２
その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛1,2,4,10,28,20,27,17,18,21,30,26,14,9,25,11,31,29,23,5,13,6,16,15,12,3,7,19,24,8,22,32｝
となる。これを図１７の左側の太枠内に示し、基本のランダム系列とする。
【００９５】
つぎに、Ｃ（ｉ）のランダム系列をある間隔Ｓ（ｊ），ｊ＝１，２，…，Ｐ−１で読む方法について説明する。この間隔は、Ｐ−１個生成可能である。ここではＰ−１＝３０である。
Ｓ（ｊ）＝｛1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30｝
【００９６】
ランダム系列を一定の間隔で飛ばし読みする系列をＬＢj（ｉ）とすると、
ＬＢj（ｉ）＝（（Ｓ（ｊ）＊ｉ） ｍｏｄ Ｐ）＋１
ただし、ｊ＝１，２，…，Ｐ−１であり、ｉ＝１，…，Ｐ−１である。
【００９７】
たとえば ｊ＝１の場合、
ＬＢ1（１）＝｛2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31｝
となる。
【００９８】
ここで、ランダム系列の系列数３２に対し、１とＰから２sまでの整数が不足している。この場合は１，３２が不足している。
ＬＢj（ｊ）＝３２，ＬＢj（３２−ｊ）＝１
を挿入する。
【００９９】
この結果、
ＬＢ1（ｉ）＝｛32,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,1,31｝
となる。
【０１００】
同様にｊ＝２の場合、
ＬＢ2（ｉ）＝｛3,32,5,7,9,11,13,15,17,19,21,23,25,27,29,31,2,4,6,8,10,12,14,16,18,20,22,24,26,1,28,30｝
【０１０１】
ここからラテン方陣を作成する。Ｌjq（ｉ）はｊ番目のラテン方陣のｑ列目の系列である。Ｌjq（ｉ）＝ＬＢj（ｉ）とすると、たとえば、
Ｌ11（ｉ）＝Ｃ（ＬＢ1（ｉ））＝｛32,2,4,10,28,20,27,17,18,21,30,26,14,9,25,11,31,29,23,5,13,6,16,15,12,3,7,19,24,8,1,22｝
となる。これは、図１８の１列目がＬ11（ｉ）に相当する。
【０１０２】
このランダム系列Ｌ11（ｉ）を巡回シフトしてＬ1q（ｉ）の３２（行数）×３２（列数）のラテン方陣を作成する。同様の手順で、ラテン方陣をＬ2q（ｉ），Ｌ3q（ｉ），…，Ｌ30q（ｉ）まで作成し、３２（行数）×（３２×３０）（列数）のラテン方陣の組み合わせを作成する。
【０１０３】
ここで、具体例を説明する。ユークリッド幾何符号ＥＧ（２，２5）の１列目の「１」の位置を示す行番号の配列Ｂ1（ｘ）から「１」の番号をランダムに抽出し、ユークリッド幾何符号の第１列目を分解する。なお、上記配列Ｂ1（ｘ）は、
Ｂ1（ｘ）＝｛1 32 114 136 149 223 260 382 402 438 467 507 574 579 588 622 634 637 638 676 717 728 790 851 861 879 947 954 971 977 979 998｝
である。
【０１０４】
たとえば、ユークリッド幾何符号の第１列目を４分割する場合（各列の重みが８となる）、上記ランダム系列のラテン方の陣組み合わせＬjq（ｉ）の１番目のラテン方陣の１列目から４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ1（ｎ）＝｛32,2,4,10,28,20,27,17｝
Ｌ2（ｎ）＝｛18,21,30,26,14,9,25,11｝
Ｌ3（ｎ）＝｛31,29,23,5,13,6,16,15,｝
Ｌ4（ｎ）＝｛12,3,7,19,24,8,1,22｝
【０１０５】
そして、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を用いて、配列Ｂ1（ｘ）から「１」の番号をランダムに抽出する。その結果Ｒ（Ｌ1（ｎ））〜Ｒ（Ｌ4（ｎ））は、
Ｒ（Ｌ1（ｎ））＝｛998,32,136,438,954,676,947,634｝
Ｒ（Ｌ2（ｎ））＝｛637,717,977,879,579,402,861,467｝
Ｒ（Ｌ3（ｎ））＝｛979,971,790,149,574,223,622,588｝
Ｒ（Ｌ4（ｎ））＝｛507,114,260,638,851,382,1,728｝
となる。以降、同様の手順でユークリッド幾何符号の全ての列を分割する。
【０１０６】
【発明の効果】
以上、説明したとおり、本発明によれば、生成関数λ（ｘ）と生成関数ρ（ｘ）を１回の線形計画法で求めることとしたため、上記論文のように、生成関数λ（ｘ）の導出と生成関数ρ（ｘ）の導出を繰り返し実行し、双方の最適値を求める方法よりも、容易かつ短時間に、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を生成することができる、という効果を奏する。
【０１０７】
つぎの発明によれば、同一レートでは「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」よりも良好な特性を得ることができる、という効果を奏する。また、規則的な分割では、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であっても大幅に特性を改善することはできないが、ランダムな分割を実施すると、特性を画期的に改善することができる、という効果を奏する。
【０１０８】
つぎの発明によれば、重み配分を、重み単位の重み総数が整数で、かつ重み単位の重み総数の総和とユークリット幾何符号の「１」の総数とが等しくなるように調整することによって、より高精度な分割処理を実現できる、という効果を奏する。
【０１０９】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１０】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１１】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を複数個作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１２】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を複数個作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【図面の簡単な説明】
【図１】 実施の形態１のＬＤＰＣ符号用検査行列生成方法を示すフローチャートである。
【図２】 ｒａｔｅ＝０．５とした場合の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【図３】 ＥＧ（２，２5）を想定し、ｄl＝３２とした場合の、分割テーブルを示す図である。
【図４】 重み配分調整用テーブルを示す図である。
【図５】 重み配分後の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【図６】 従来の分割手順を示す図である。
【図７】 分割前のＥＧ（２，２5）のグラフを示す図である。
【図８】 分割後のＥＧ（２，２5）のグラフを示す図である。
【図９】 Ｅｂ／Ｎｏと誤り率特性との関係を示す図である。
【図１０】 「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【図１１】 「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【図１２】 ランダム系列のラテン方陣を示す図である。
【図１３】 ＬＤＰＣ符号化／復号システムを示す図である。
【図１４】 従来のＬＤＰＣ符号用の検査行列を示す図である。
【図１５】 ユークリッド幾何符号ＥＧ（２，２2）の構成を示す図である。
【図１６】 列の重みを４から２に規則的に分離した例を示す図である。
【図１７】 ランダム系列のラテン方陣を示す図である。
【図１８】 ランダム系列のラテン方陣を示す図である。
【符号の説明】
１０１ 符号化器、１０２ 変調器、１０３ 通信路、１０４ 復調器、１０５ 復号器。
【発明の属する技術分野】
本発明は、誤り訂正符号としてＬＤＰＣ（Low-Density Parity-Check）符号を採用した符号化器におけるＬＤＰＣ符号用検査行列生成方法に関するものである。
【０００２】
【従来の技術】
図１３は、ＬＤＰＣ符号化／復号システムを示す図である。図１３において、１０１は符号化器であり、１０２は変調器であり、１０３は通信路であり、１０４は復調器であり、１０５は復号器である。ここでは、従来のＬＤＰＣ符号用検査行列生成方法を説明する前に、ＬＤＰＣ符号を使用した場合の符号化，復号の流れについて説明する。
【０００３】
まず、送信側の符号化器１０１では、後述する所定の方法で検査行列Ｈを生成する。そして、以下の条件に基づいて生成行列Ｇを求める。
Ｇ：ｋ×ｎ行列（ｋ：情報長，ｎ：符号語長）
ＧＨT＝０（Ｔは転置行列）
【０００４】
その後、符号化器１０１では、情報長ｋのメッセージ（ｍ1ｍ2…ｍk）を受け取り、上記生成行列Ｇを用いて符号語Ｃを生成する。
Ｃ＝（ｍ1ｍ2…ｍk）Ｇ
＝（ｃ1ｃ2…ｃn） （ただし、Ｈ（ｃ1ｃ2…ｃn）T＝０）
【０００５】
そして、変調器１０２では、生成した符号語Ｃに対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル変調を行い、送信する。
【０００６】
一方、受信側では、復調器１０４が、通信路１０３を介して受け取った変調信号に対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル復調を行い、さらに、復号器１０５が、ＬＤＰＣ符号化された復調結果に対して、「ｓｕｍ−ｐｒｏｄｕｃｔアルゴリズム」によるくり返し復号を実施し、推定結果（もとのｍ1ｍ2…ｍkに対応）を出力する。
【０００７】
以下、従来のＬＤＰＣ符号用検査行列生成方法について説明する。ＬＤＰＣ符号用の検査行列としては、たとえば、ＬＤＰＣの発案者Ｇａｌｌａｇｅｒにより以下のような行列が提案されている（図１４参照）。
【０００８】
図１４に示す行列は、「１」と「０」の２値の行列で、「１」の部分を塗りつぶしている。他の部分は全て「０」である。この行列は、１行の「１」の数（これを行の重みと表現する）が４で、１列の「１」の数（これを列の重みと表現する）が３であり、全ての列と行の重みが均一なため、これを一般に「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」と呼んでいる。また、Ｇａｌｌａｇｅｒの符号では、たとえば、図１４に示すように、行列を３ブロックに分け、２ブロック目と３ブロック目に対してランダム置換を行っている。
【０００９】
しかしながら、このランダム置換には、所定のルールが無いため、より特性の良好な符号を見つけるためには、計算機による時間のかかる探索を行わなければならなかった。
【００１０】
そこで、たとえば、計算機探索によらなくても確定的に行列を生成でき、比較的安定した良好な特性を示すＬＤＰＣ符号として、ユークリット幾何符号を用いる方法が、Ｙ．Ｋｏｕ等（Y. Kou, S. Lin, and M. P. C. Fossorier, "Low Density Parity Check Codes Based on Finite Geometries: A Rediscovery," ISIT 2000, pp. 200, Sorrento, Itary, June 25-30, 2000.）によって提案された。この方法では、規則的なｅｎｓｅｍｂｌｅ（アンサンブル）で構成された「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」について説明されている。
【００１１】
ここでは、有限幾何符号の一種であるユークリッド幾何符号ＥＧ（２，２6）を用いてＬＤＰＣ符号の検査行列を生成する方法について提案され、誤り率１０-4点において、シャノン限界から１．４５ｄＢに接近した特性を得ている。図１５は、たとえば、ユークリッド幾何符号ＥＧ（２，２2）の構成を示す図であり、行，列のそれぞれの重みが４，４の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」構造をしている。
【００１２】
したがって、ユークリッド幾何符号ＥＧ（ｍ，２s）の場合、その特性は、以下のように規定される。
符号長： ｎ＝２2s−１
冗長ビット長： ｎ−ｋ＝３s−１
情報長： ｋ＝２2s−３s
最小距離： ｄmin＝２s＋１
密度： ｒ＝２s／（２2s−１）
【００１３】
図１５を見ても分かるように、ユークリッド幾何符号は、各行の「１」の配置が行毎に巡回シフトした構造になっており、符号が容易にかつ確定的に構成できる特長がある。
【００１４】
Ｙ．Ｋｏｕらによる検査行列の生成方法では、さらに、上記ユークリット幾何符号に基づいて行と列の重みを変更し、行，列を必要に応じて拡張している。たとえば、ＥＧ（２，２2）の列の重みを１／２に分離する場合、Ｙ．Ｋｏｕらの論文では、１列内に４つある重みを１つ置きに２個づつ分離する。図１６は、列の重みを４から２に規則的に分離した例を示す図である。
【００１５】
一方、上記「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」の特性よりも「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の特性の方が良好であることが、Ｌｕｄｙ等（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.）により報告された。なお、上記「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」は、列と行の重みがそれぞれあるいはどちらか一方が均一でないＬＤＰＣ符号を表す。
【００１６】
そして、それは、Ｒｉｃｈａｒｄｓｏｎ等（T. J. Richardson and R. Urbanke, "The capacity of low-density parity-check codes under message-passing decoding," IEEE Trans. Inform. Theory, vol.47,No.2, pp.599-618, Feb. 2001.）、あるいはＣｈｕｎｇ等（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.）によって理論的に解析された。
【００１７】
特に、Ｃｈｕｎｇ等は、繰り返し復号器における入力と出力の対数尤度比（ＬＬＲ）がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析し、良好な行と列の重みのアンサンブルを求めている。
【００１８】
【発明が解決しようとする課題】
しかしながら、たとえば、上記Ｃｈｕｎｇ等による従来のＬＤＰＣ符号用検査行列生成方法は、行内の「１」の点の数（後述するバリアブルノードの次数配分に相当）と、列内の「１」の点の数（後述するチェックノードの次数配分に相当）と、の両方を変数として、下記の（１）式（ｒａｔｅ：符号化率）が最大となるバリアブルノードの次数配分およびチェックノードの次数配分を求めている。すなわち、ＳＮＲ（Signal to Noise Ratio）が最小となるアンサンブルを線形計画法により探索している。
【００１９】
【数１】
そのため、上記「ｒａｔｅ」の最大値により得られる検査行列が流動的になり、特性が安定しない、という問題があった。また、従来のＬＤＰＣ符号用検査行列生成方法は、バリアブルノードの次数配分の導出とチェックノードの次数配分の導出とを所定回数にわたって繰り返し行っているため、探索処理にある程度の時間を要する、という問題もあった。
【００２１】
本発明は、上記に鑑みてなされたものであって、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を、短時間で容易に探索可能なＬＤＰＣ符号用検査行列生成方法を得ることを目的とする。
【００２２】
【課題を解決するための手段】
上述した課題を解決し、目的を達成するために、本発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、復号器における入出力データの対数尤度比がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析することによって、誤りが０となるＳＮＲの限界（threshold）を求めることとし、さらに、符号化レートを固定した状態で、かつガウスノイズが最大になるように、行の重みと列の重みの最適なアンサンブル（thresholdが最小となるアンサンブル）を１回の線形計画法で探索し、当該アンサンブルにしたがってＬＤＰＣ符号用の検査行列を生成することを特徴とする。
【００２３】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、前記アンサンブルを探索後、当該探索結果に基づいてユークリット幾何符号の各行または各列からランダムに「１」を抽出し、各行または各列を分割することによって、Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号の検査行列を生成することを特徴とする。
【００２４】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、前記アンサンブルの重み配分を、重み単位の重み総数が整数で、かつ重み単位の重み総数の総和とユークリット幾何符号の「１」の総数とが等しくなるように調整し、調整後のアンサンブルに基づいて分割処理を行うことを特徴とする。
【００２５】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を作成し、ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする。
【００２６】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を作成し、ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする。
【００２７】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を複数個作成し、列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする。
【００２８】
つぎの発明にかかるＬＤＰＣ符号用検査行列生成方法にあっては、基本のランダム系列のラテン方陣を複数個作成し、列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする。
【００２９】
【発明の実施の形態】
以下に、本発明にかかるＬＤＰＣ符号用検査行列生成方法の実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
【００３０】
実施の形態１．
本実施の形態のＬＤＰＣ符号用検査行列生成方法を説明する前に、本実施の形態のＬＤＰＣ符号用検査行列生成方法を実現可能な符号化器の位置付け、および「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の従来の検査行列生成方法について説明する。なお、ＬＤＰＣ符号化／復号システムの構成については、先に説明した図１３と同様である。
【００３１】
送信側の符号化器１０１では、後述する本実施の形態のＬＤＰＣ符号用検査行列生成方法で検査行列Ｈを生成する。そして、以下の条件に基づいて生成行列Ｇを求める。
Ｇ：ｋ×ｎ行列（ｋ：情報長，ｎ：符号語長）
ＧＨT＝０（Ｔは転置行列）
【００３２】
その後、符号化器１０１では、情報長ｋのメッセージ（ｍ1ｍ2…ｍk）を受け取り、上記生成行列Ｇを用いて符号語Ｃを生成する。
Ｃ＝（ｍ1ｍ2…ｍk）Ｇ
＝（ｃ1ｃ2…ｃn） （ただし、Ｈ（ｃ1ｃ2…ｃn）T＝０）
【００３３】
そして、変調器１０２では、生成した符号語Ｃに対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル変調を行い、送信する。
【００３４】
一方、受信側では、復調器１０４が、通信路１０３を介して受け取った変調信号に対して、ＢＰＳＫ，ＱＰＳＫ，多値ＱＡＭなどのデジタル復調を行い、さらに、復号器１０５が、ＬＤＰＣ符号化された復調結果に対して、「ｓｕｍ−ｐｒｏｄｕｃｔアルゴリズム」による繰り返し復号を実施し、推定結果（もとのｍ1ｍ2…ｍkに対応）を出力する。
【００３５】
つぎに、Ｃｈｕｎｇ等（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.）によって理論的に解析された、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の従来の検査行列生成方法について詳細に説明する。ここでは、繰り返し復号器における入力と出力の対数尤度比（ＬＬＲ）がガウス分布に近似できると仮定してＬＤＰＣ符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析し、良好な行と列の重みのアンサンブルを求めている。
【００３６】
なお、上記論文に記述されたＬＤＰＣ符号用検査行列生成方法であるガウス近似法（Gaussian Approximation）では、前提として、検査行列における行内の「１」の点をバリアブルノードと定義し、列内の「１」の点をチェックノードと定義する。
【００３７】
まず、チェックノードからバリアブルノードへのＬＬＲメッセージ伝搬を解析する。０＜ｓ＜∞と０≦ｔ＜∞という条件において、以下の関数（２）式を定義する。なお、ｓ＝ｍu0はｕ０の平均値であり、ｕ０は分散値σn 2のガウスノイズを含む伝送路を経由して受信した信号の対数尤度比（ＬＬＲ）であり、ｔは所定の繰り返しの時点におけるチェックノードのＬＬＲ出力値のアンサンブル平均である。
【００３８】
【数２】
なお、上記λ（ｘ）およびρ（ｘ）は、それぞれバリアブルノードおよびチェックノードの次数配分（バリアブルノードとチェックノードの各１行，各１列内の「１」の数を次数と表現する）の生成関数を表し、（３）式および（４）式のように表すことができる。また、λiとρiは、それぞれ次数ｉのバリアブルノードとチェックノードに属するエッジの比率を表す。また、ｄlは最大バリアブルノードの次数であり、ｄrは最大チェックノードの次数である。
【００４０】
【数３】
【数４】
ただし、φ（ｘ）は下記（５）式のように定義する。
【００４３】
【数５】
そして、（２）式は、等価的に下記（６）式と表すことができる。
【００４５】
【数６】
なお、ｔlはｌ番目の繰り返し時点におけるチェックノードのＬＬＲ出力値のアンサンブル平均である。
【００４７】
ここで、誤りが０となりうるＳＮＲの限界（threshold）を求めるための条件は、ｌ→∞のときにｔl（ｓ）→∞（Ｒ+と表現する）となることであり、この条件を満たすためには、以下の条件（７）式を満たす必要がある。
【００４８】
【数７】
つぎに、バリアブルノードからチェックノードへのＬＬＲメッセージ伝搬を解析する。０＜ｓ＜∞と０＜ｒ≦１という条件において、以下の関数（８）式を定義する。なお、ｒの初期値ｒ0はφ（ｓ）である。
【００５０】
【数８】
そして、（８）式は、等価的に下記（９）式と表すことができる。
【００５２】
【数９】
ここで、誤りが０となりうるＳＮＲの限界（threshold）を求めるための条件は、ｒl（ｓ）→０となることであり、この条件を満たすためには、以下の条件（１０）式を満たす必要がある。
【００５４】
【数１０】
さらに、上記Ｃｈｕｎｇ等の論文では、上記式を用いて以下の手順でバリアブルノードとチェックノードの最適な次数を探索している。
（１）生成関数λ（ｘ）とガウスノイズσnが与えられていると仮定し、生成関数ρ（ｘ）を変数として、前述した（１）式が最大となる点を探索する。なお、この探索における拘束条件は、ρ（１）＝１と正規化することと、上記（７）式を満たすことである。
（２）生成関数ρ（ｘ）とガウスノイズσnが与えられていると仮定し（たとえば、（１）の結果より得られる値）、生成関数λ（ｘ）を変数として、（１）式が最大となる点を探索する。なお、この探索における拘束条件は、λ（１）＝１と正規化することと、上記（１０）式を満たすことである。
（３）最大「ｒａｔｅ」を求めるために、上記（１）と上記（２）を繰り返し実行し、生成関数λ（ｘ）と生成関数ρ（ｘ）のより良好なアンサンブルを線形計画法で探索する。
（４）最後に、ガウスノイズσnより信号電力を１と正規化して、ＳＮＲの限界（threshold）を求める。
【００５６】
【数１１】
しかしながら、上記Ｃｈｕｎｇ等の論文では、「ｒａｔｅ（符号化率）」の最大値により得られる検査行列が流動的になり、設計において仕様として固定されるｒａｔｅが変動してしまうため、実設計にむかない、という問題があった。また、上記ガウス近似法では、バリアブルノードの次数配分の導出とチェックノードの次数配分の導出とを所定回数にわたって繰り返し行っているため、探索処理にある程度の時間を要する、という問題もあった。
【００５８】
そこで、本実施の形態においては、確定的でかつ特性が安定した「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」用の検査行列を、短時間で容易に探索する。図１は、実施の形態１のＬＤＰＣ符号用検査行列生成方法を示すフローチャートである。
（１）「ｒａｔｅ」が与えられているものと仮定する。すなわち、要求「ｒａｔｅ」を固定する。実際の設計では、目標「ｒａｔｅ」が予め指定されている場合が多いためである。
（２）生成関数λ（ｘ）と生成関数ρ（ｘ）を同時に変数として扱い、ガウスノイズσnが最大になるように、線形計画法で最適な生成関数λ（ｘ）と生成関数ρ（ｘ）を探索する。この探索における拘束条件は、λ（１）＝１，ρ（１）＝１と正規化し、さらに上記（１０）式を満たすことである。
【００５９】
このように、本実施の形態では、上記（９）式と上記（１０）式を満たす生成関数λ（ｘ）と生成関数ρ（ｘ）を１回の線形計画法で求めることとしたため、上記論文のように、生成関数λ（ｘ）の導出と生成関数ρ（ｘ）の導出を繰り返し実行し、双方の最適値を求める方法よりも、容易かつ短時間に、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を生成することができる。
【００６０】
実施の形態２．
実施の形態２では、前述の実施の形態１においてユークリッド幾何符号を採用し、１行あるいは１列の「１」の配置を分割することによって、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の検査行列を生成する。
【００６１】
まず、実施の形態１におけるＬＤＰＣ符号用検査行列生成方法により、生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを導出する。図２は、ｒａｔｅ＝０．５とした場合の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。なお、σGAはガウス近似法により導出した「threshold」時のノイズ分散値を表し、ＳＮＲnorm（ＧＡ）はガウス近似法により導出した「threshold」のＳＮＲとシャノン限界のＳＮＲとの差分を表し、ｘは重みを表し、λxおよびρxはそれぞれバリアブルノードとチェックノードの重み配分を表す。
【００６２】
また、基準になるユークリッド幾何符号は、ＥＧ（２，２5）を想定し、ｄl＝３２とした。また、重み配分λxのｘの値および重み配分ρxのｘの値は、組み合わせで３２（ｄl）を構成できる値とする。
【００６３】
図３は、ＥＧ（２，２5）を想定し、ｄl＝３２とした場合の、分割テーブルを示す図である。図３に示すとおり、図２のｘは、組み合わせにより必ず３２となる。たとえば、図示の７ｘ４と２ｘ２の組み合わせは、重みが３２の１列を、重みが７の４列と重みが２の２列に分割可能なことを表している。このように、ＥＧ（２，２5）の符号を基本として、図３のように重みが３２の各行列を適切に分割すると、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」の検査行列を構成することができる。
【００６４】
ここで、分割処理を行う前に、図２に示す生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルの重み配分を以下の手順で調整する。図４は、重み配分調整用テーブルを示す図である。なお、ＥＧ（２，２5）のユークリッド幾何符号は、１０２３行×１０２３列で構成される。
【００６５】
（１）ガウス近似法で求めた生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブル（表１参照）をテーブルの２列目と３列目に設定する。
（２）重み配分λxおよびρx（３列目）と、ＥＧ（２，２5）における全行列の「１」の総数ＴＰ＝３２７３６と、を乗算し、重み単位の重み総数を求め、さらに、当該重み単位の重み総数とその総和を４列目に設定する。
（３）重み単位の重み総数（４列目）を対応する重みｘで割り、重み単位の総列数を求め、それを５列目に設定する。
（４）重み単位の総列数が小数点以下を含む場合、丸め処理（四捨五入，切上げ，切捨て等）を行い、その結果を６列目に設定する。
（５）丸め処理後の重み単位の総列数（６列目）と対応する重みｘとを乗算し、丸め処理後の重み単位の重み総数を求め、それを７列目に設定する。そして、各重み総数の総和（７列目の合計の行）が行列内の「１」の総数（ＴＰ＝３２７３６）と等しいかどうかを確認する。
（６）行列内の「１」の総数に等しくない場合、丸め処理後の重み単位の重み総数（７列目）を整数単位で調整し、その結果を８列目に設定する。この場合、８列目の総和が、行列内の「１」の総数（ＴＰ＝３２７３６）に等しくなるように調整する。
（７）調整後の重み単位の重み総数（８列目）を対応する重みｘで割り、調整後の重み単位の総列数を求め、それを９列目に設定する。調整後の各重みの配分（１１列目）は、可能な限りガウス近似法で求めた値（３列目）に近い値にする。
【００６６】
図５は、重み配分後の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【００６７】
つぎに、ユークリッド幾何符号における１行あるいは１列の分割手順について説明する。
【００６８】
たとえば、分割手順に関して、Ｙ．Ｋｏｕ等の論文では、規則的に分割する方法を提示している。図６は、上記論文における分割手順を示す図である。まず、図６に示すように行列のナンバリングを行う。ここでは、列番号を左端から順に１，２，３，…とし、行番号を上から順に１，２，３，…とする。そして、たとえば、３２点×１列を８点×４列に分割する場合、下記（１２）式にしたがって規則的に分割する。
【００６９】
【数１２】
なお、ｍ＝１，２，３，４とし、ｎ＝０，１，２，３，４，５，６，７とし、ｌはＥＧ（２，２5）の列番号を表す。また、Ｂl（ｘ）はＥＧ（２，２5）のｌ列目の「１」の位置を表し、Ｓm（ｎ）は分割後の行列のｍ列目の「１」の位置を表す。
【００７１】
具体的にいうと、ＥＧ（２，２5）における１列中の「１」の位置を示す行番号は、
Ｂl（ｘ）＝｛1 32 114 136 149 223 260 382 402 438 467 507 574 579 588 622 634 637 638 676 717 728 790 851 861 879 947 954 971 977 979 998｝
となり、その結果、分割後の行列における１〜４列目の「１」の位置を示す行番号は、Ｂl（ｘ）から「１」の番号が規則的に抽出され、
Ｓ1（ｎ）＝｛1 149 402 574 634 717 861 971｝
Ｓ2（ｎ）＝｛32 223 438 579 637 728 879 977｝
Ｓ3（ｎ）＝｛114 260 467 588 638 790 947 979｝
Ｓ4（ｎ）＝｛136 382 507 622 676 851 954 998｝
となる。すなわち、３２点×１列が８点×４列に分割される。
【００７２】
一方、本実施の形態におけるユークリッド幾何符号の分割処理は、上記のように規則的に分割するのではなく、Ｂl（ｘ）から「１」の番号をランダムに抽出する。なお、この抽出処理は、ランダム性が保持されるのであればどのような方法を用いてもよい。
【００７３】
したがって、分割後の行列のｍ列目の「１」の位置の一例をＲm（ｎ）とした場合、Ｒm（ｎ）は、
Ｒ1（ｎ）＝｛1 114 574 637 851 879 977 979｝
Ｒ2（ｎ）＝｛32 136 402 467 588 728 861 971｝
Ｒ3（ｎ）＝｛149 260 382 438 579 638 717 998｝
Ｒ4（ｎ）＝｛223 507 622 634 676 790 947 954｝
となる。
【００７４】
上記のような本実施の形態の分割手順をグラフ上で表現すると、以下のように表現することができる。図７は、分割前のＥＧ（２，２5）のグラフを示す図である。なお、両ノードを結ぶ線はエッジと表現する。図７では、分割前の１０２３行×１０２３列（各行列の重みがそれぞれ３２）のユークリット幾何符号を表現している。また、図８は、ＥＧ（２，２5）のエッジをランダムに選択した、分割後のグラフを示す図である。
【００７５】
ここで、上記で説明したＬＤＰＣ符号の特性を比較する。図９は、Ｅｂ／Ｎｏ（情報１ビットあたりの信号電力対ノイズ電力比）と誤り率特性（ＢＥＲ）との関係を示す図である。なお、繰り返し回数は５０回で、復号法は「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」である。
【００７６】
なお、図中"Simple regular extended EG(2,25)"は、Ｙ．Ｋｏｕ等の発案によるＥＧ（２，２5）の規則的な列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」であり、"Random regular extended EG(2,25)"は、本実施の形態によるＥＧ（２，２5）のランダムな列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」である。図１０は、これらの「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【００７７】
また、図中"Simple irregular extended EG(2,25)"は、実施の形態１の方法によって特定されたアンサンブルに対して、Ｙ．Ｋｏｕ等の発案によるＥＧ（２，２5）の規則的な列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であり、"Random irregular extended EG(2,25)"は、実施の形態１の方法によって特定されたアンサンブルに対して、本実施の形態によるＥＧ（２，２5）のランダムな列の分割を実施した場合の、ｒａｔｅ＝０．５の「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」である。図１１は、これらの「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【００７８】
図９からわかるように、同一レートでは、「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」より「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のほうが性能がよい。また、Ｙ．Ｋｏｕ等の論文のような規則的な分割では、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であっても大幅な改善は見込めないが、本実施の形態のランダムな分割を実施すると性能が画期的に改善される。
【００７９】
つぎに、上記ランダム分割の一例を詳細に説明する。ここでは、ランダム分割を行う場合のランダム系列を容易かつ確定的に生成する。この方法による利点は、送信側と受信側が同じランダム系列を生成できることにある。これは、現実のシステムではきわめて重要となる。また、符号特性の条件が正確に規定できる、という利点もある。
【００８０】
（１）基本のランダム系列を作成する。
以下に、ランダム系列作成の一例を記述する。ここでは、説明の便宜上、ユークリッド幾何符号ＥＧ（２，２4）を用いる。ユークリッド幾何符号ＥＧ（２，２4）の場合、１行に存在する「１」の数は２4個である。
【００８１】
ＰをＰ≧２sを満たす最小の素数とした場合、たとえば、２4のときはＰ＝１７となる。ここで、系列長Ｐ−１＝１６の基本のランダム系列Ｃ（ｉ）を（１３）式にしたがって作成する。
Ｃ（１）＝１
Ｃ（ｉ＋１）＝Ｇ0×Ｃ（ｉ） ｍｏｄ Ｐ …（１３）
ただし、ｉ＝１，…，Ｐ−１とし、Ｇ0はガロア体ＧＦ（Ｐ）の原始元である。その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛1 3 9 10 13 5 15 11 16 14 8 7 4 12 2 6｝
となる。なお、２4の場合は、Ｐが２4＋１となり、ランダム系列長が１６になるので問題ないが、２5の場合等は、Ｐが２5＋１以上になりランダム系列長が２5を超えてしまう。そのような場合は、２5を超える数値をランダム系列から削除することで対応する。
【００８２】
そして、ランダム系列Ｃ（ｉ）を巡回シフトして、ランダム系列のラテン方陣を生成する。図１２は、ランダム系列のラテン方陣を示す図である。
【００８３】
（２）ランダム系列のラテン方陣を用いて、ユークリッド幾何符号ＥＧ（２，２4）の１列目の「１」の位置を示す行番号の配列Ｂl（ｘ）から、「１」の番号をランダムに抽出し、ユークリッド幾何符号の第１列目を分割する。なお、上記配列Ｂl（ｘ）は、
Ｂl（ｘ）＝｛1 14 16 19 45 49 55 107 115 126 127 182 210 224 231 247｝
である。
【００８４】
たとえば、ユークリッド幾何符号の第１列目を４分割する場合（各列の重みが４となる）、上記ランダム系列のラテン方陣の１行目から、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ1（ｎ）＝｛1,3,9,10｝
Ｌ2（ｎ）＝｛13,5,15,11｝
Ｌ3（ｎ）＝｛16,14,8,7｝
Ｌ4（ｎ）＝｛4,12,2,6｝
【００８５】
そして、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を用いて、配列Ｂl（ｘ）から、「１」の番号をランダムに抽出する。その結果、Ｒ（Ｌ1（ｎ））〜Ｒ（Ｌ4（ｎ））は、
Ｒ（Ｌ1（ｎ））＝｛1,16,115,126｝
Ｒ（Ｌ2（ｎ））＝｛210,45,231,127｝
Ｒ（Ｌ3（ｎ））＝｛247,224,107,55｝
Ｒ（Ｌ4（ｎ））＝｛19,182,14,49｝
となる。
【００８６】
（３）同様に、ランダム系列のラテン方陣を用いて、ユークリッド幾何符号ＥＧ（２，２4）の２列目の「１」の位置を示す行番号の配列Ｂ2（ｘ）から、「１」の番号をランダムに抽出し、ユークリッド幾何符号の第２列目を分割する。なお、上記配列Ｂ2（ｘ）は、
Ｂ2（ｘ）＝｛2 15 17 20 46 50 56 108 116 127 128 183 211 225 232 248｝
である。
【００８７】
たとえば、ユークリッド幾何符号の第２列目を４分割する場合（各列の重みが４となる）、上記ランダム系列のラテン方陣の２行目から、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ5（ｎ）＝｛3,9,10,13｝
Ｌ6（ｎ）＝｛5,15,11,16｝
Ｌ7（ｎ）＝｛14,8,7,4｝
Ｌ8（ｎ）＝｛12,2,6,1｝
【００８８】
そして、４つのランダム系列Ｌ5（ｎ）〜Ｌ8（ｎ）を用いて、配列Ｂ2（ｘ）から、「１」の番号をランダムに抽出する。その結果、Ｒ（Ｌ5（ｎ））〜Ｒ（Ｌ8（ｎ））は、
Ｒ（Ｌ5（ｎ））＝｛17,116,127,211｝
Ｒ（Ｌ6（ｎ））＝｛46,232,128,248｝
Ｒ（Ｌ7（ｎ））＝｛225,108,56,20｝
Ｒ（Ｌ8（ｎ））＝｛183,15,50,2｝
となる。以降、同様の手順でユークリッド幾何符号の全ての列を分割する。
【００８９】
（４）なお、上記ランダム系列のラテン方陣でランダム系列が不足するような場合は、（１３）式の基本のランダム系列のＧ0に原始元より大きくかつ２s以下の素数を代入して再度基本のランダム系列を作り、同様の手順で分割する。
【００９０】
つぎに、もう一つのランダム系列のラテン方陣の作成法について説明する。ここでは、説明の便宜上、ユークリッド幾何符号ＥＧ（２，２5）を用いる。ユークリッド幾何符号ＥＧ（２，２5）の場合、１行に存在する「１」の数は２5個である。
【００９１】
上記例では、ＥＧ（２，２4）のため１６（行数）×１６（列数）のラテン方陣を作成したが、これでは不足する場合がある。ここでは、ＥＧ（２，２5）を例にとり、３２（行数）×９６０（列数）まで拡張する。行列のサイズは、３２（行数）×３２（列数）のラテン方陣の組を３０組列方向に並べたことにより決められる。
【００９２】
また、上記の例では、ＰをＰ≧２sを満たす最小の素数としたが、ここでは、ＰをＰ≧２sを満たす最大の素数とする。たとえば、２5のときはＰ＝３１となる。
【００９３】
ここで、系列長Ｐ＝３１の基本のランダム系列Ｃ（ｉ）を（１４）式に従って作成する。
Ｃ（１）＝０
Ｃ（ｉ＋１）＝Ｇo×Ｃ（ｉ） ｍｏｄ Ｐ …（１４）
ただし、ｉ＝１，…，Ｐ−１とし、Ｇoはガロア体ＧＦ（Ｐ）の原始元である。その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛0,1,3,9,27,19,26,17,20,29,25,13,8,24,10,30,28,22,4,12,5,15,14,11,2,6,18,23,7,21｝
となる。
【００９４】
つぎに以下の操作をする。
Ｃ（ｉ）＝Ｃ（ｉ）＋１， Ｃ（３２）＝３２
その結果、Ｃ（ｉ）は、
Ｃ（ｉ）＝｛1,2,4,10,28,20,27,17,18,21,30,26,14,9,25,11,31,29,23,5,13,6,16,15,12,3,7,19,24,8,22,32｝
となる。これを図１７の左側の太枠内に示し、基本のランダム系列とする。
【００９５】
つぎに、Ｃ（ｉ）のランダム系列をある間隔Ｓ（ｊ），ｊ＝１，２，…，Ｐ−１で読む方法について説明する。この間隔は、Ｐ−１個生成可能である。ここではＰ−１＝３０である。
Ｓ（ｊ）＝｛1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30｝
【００９６】
ランダム系列を一定の間隔で飛ばし読みする系列をＬＢj（ｉ）とすると、
ＬＢj（ｉ）＝（（Ｓ（ｊ）＊ｉ） ｍｏｄ Ｐ）＋１
ただし、ｊ＝１，２，…，Ｐ−１であり、ｉ＝１，…，Ｐ−１である。
【００９７】
たとえば ｊ＝１の場合、
ＬＢ1（１）＝｛2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31｝
となる。
【００９８】
ここで、ランダム系列の系列数３２に対し、１とＰから２sまでの整数が不足している。この場合は１，３２が不足している。
ＬＢj（ｊ）＝３２，ＬＢj（３２−ｊ）＝１
を挿入する。
【００９９】
この結果、
ＬＢ1（ｉ）＝｛32,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,1,31｝
となる。
【０１００】
同様にｊ＝２の場合、
ＬＢ2（ｉ）＝｛3,32,5,7,9,11,13,15,17,19,21,23,25,27,29,31,2,4,6,8,10,12,14,16,18,20,22,24,26,1,28,30｝
【０１０１】
ここからラテン方陣を作成する。Ｌjq（ｉ）はｊ番目のラテン方陣のｑ列目の系列である。Ｌjq（ｉ）＝ＬＢj（ｉ）とすると、たとえば、
Ｌ11（ｉ）＝Ｃ（ＬＢ1（ｉ））＝｛32,2,4,10,28,20,27,17,18,21,30,26,14,9,25,11,31,29,23,5,13,6,16,15,12,3,7,19,24,8,1,22｝
となる。これは、図１８の１列目がＬ11（ｉ）に相当する。
【０１０２】
このランダム系列Ｌ11（ｉ）を巡回シフトしてＬ1q（ｉ）の３２（行数）×３２（列数）のラテン方陣を作成する。同様の手順で、ラテン方陣をＬ2q（ｉ），Ｌ3q（ｉ），…，Ｌ30q（ｉ）まで作成し、３２（行数）×（３２×３０）（列数）のラテン方陣の組み合わせを作成する。
【０１０３】
ここで、具体例を説明する。ユークリッド幾何符号ＥＧ（２，２5）の１列目の「１」の位置を示す行番号の配列Ｂ1（ｘ）から「１」の番号をランダムに抽出し、ユークリッド幾何符号の第１列目を分解する。なお、上記配列Ｂ1（ｘ）は、
Ｂ1（ｘ）＝｛1 32 114 136 149 223 260 382 402 438 467 507 574 579 588 622 634 637 638 676 717 728 790 851 861 879 947 954 971 977 979 998｝
である。
【０１０４】
たとえば、ユークリッド幾何符号の第１列目を４分割する場合（各列の重みが８となる）、上記ランダム系列のラテン方の陣組み合わせＬjq（ｉ）の１番目のラテン方陣の１列目から４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を抽出する。
Ｌ1（ｎ）＝｛32,2,4,10,28,20,27,17｝
Ｌ2（ｎ）＝｛18,21,30,26,14,9,25,11｝
Ｌ3（ｎ）＝｛31,29,23,5,13,6,16,15,｝
Ｌ4（ｎ）＝｛12,3,7,19,24,8,1,22｝
【０１０５】
そして、４つのランダム系列Ｌ1（ｎ）〜Ｌ4（ｎ）を用いて、配列Ｂ1（ｘ）から「１」の番号をランダムに抽出する。その結果Ｒ（Ｌ1（ｎ））〜Ｒ（Ｌ4（ｎ））は、
Ｒ（Ｌ1（ｎ））＝｛998,32,136,438,954,676,947,634｝
Ｒ（Ｌ2（ｎ））＝｛637,717,977,879,579,402,861,467｝
Ｒ（Ｌ3（ｎ））＝｛979,971,790,149,574,223,622,588｝
Ｒ（Ｌ4（ｎ））＝｛507,114,260,638,851,382,1,728｝
となる。以降、同様の手順でユークリッド幾何符号の全ての列を分割する。
【０１０６】
【発明の効果】
以上、説明したとおり、本発明によれば、生成関数λ（ｘ）と生成関数ρ（ｘ）を１回の線形計画法で求めることとしたため、上記論文のように、生成関数λ（ｘ）の導出と生成関数ρ（ｘ）の導出を繰り返し実行し、双方の最適値を求める方法よりも、容易かつ短時間に、確定的でかつ特性が安定したＬＤＰＣ符号用の検査行列を生成することができる、という効果を奏する。
【０１０７】
つぎの発明によれば、同一レートでは「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」よりも良好な特性を得ることができる、という効果を奏する。また、規則的な分割では、「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」であっても大幅に特性を改善することはできないが、ランダムな分割を実施すると、特性を画期的に改善することができる、という効果を奏する。
【０１０８】
つぎの発明によれば、重み配分を、重み単位の重み総数が整数で、かつ重み単位の重み総数の総和とユークリット幾何符号の「１」の総数とが等しくなるように調整することによって、より高精度な分割処理を実現できる、という効果を奏する。
【０１０９】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１０】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１１】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を複数個作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【０１１２】
つぎの発明によれば、送信側と受信側が同じランダム系列を生成することができる、という効果を奏する。また、ランダム系列のラテン方陣を複数個作成することによって、符号特性の条件を正確に規定できる、という効果を奏する。
【図面の簡単な説明】
【図１】 実施の形態１のＬＤＰＣ符号用検査行列生成方法を示すフローチャートである。
【図２】 ｒａｔｅ＝０．５とした場合の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【図３】 ＥＧ（２，２5）を想定し、ｄl＝３２とした場合の、分割テーブルを示す図である。
【図４】 重み配分調整用テーブルを示す図である。
【図５】 重み配分後の生成関数λ（ｘ）と生成関数ρ（ｘ）のアンサンブルを示す図である。
【図６】 従来の分割手順を示す図である。
【図７】 分割前のＥＧ（２，２5）のグラフを示す図である。
【図８】 分割後のＥＧ（２，２5）のグラフを示す図である。
【図９】 Ｅｂ／Ｎｏと誤り率特性との関係を示す図である。
【図１０】 「Ｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【図１１】 「Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号」のアンサンブルを示す図である。
【図１２】 ランダム系列のラテン方陣を示す図である。
【図１３】 ＬＤＰＣ符号化／復号システムを示す図である。
【図１４】 従来のＬＤＰＣ符号用の検査行列を示す図である。
【図１５】 ユークリッド幾何符号ＥＧ（２，２2）の構成を示す図である。
【図１６】 列の重みを４から２に規則的に分離した例を示す図である。
【図１７】 ランダム系列のラテン方陣を示す図である。
【図１８】 ランダム系列のラテン方陣を示す図である。
【符号の説明】
１０１ 符号化器、１０２ 変調器、１０３ 通信路、１０４ 復調器、１０５ 復号器。
Claims (7)
- 復号器における入出力データの対数尤度比がガウス分布に近似できると仮定してＬＤＰＣ（Low-Density Parity-Check）符号の「Ｓｕｍ−Ｐｒｏｄｕｃｔアルゴリズム」を解析することによって、誤りが０となるＳＮＲの限界（threshold）を求めるＬＤＰＣ符号用検査行列生成方法において、
符号化レートを固定した状態で、かつガウスノイズが最大になるように、行の重みと列の重みの最適なアンサンブル（thresholdが最小となるアンサンブル）を１回の線形計画法で探索し、当該アンサンブルにしたがってＬＤＰＣ符号用の検査行列を生成することを特徴とするＬＤＰＣ符号用検査行列生成方法。 - 前記アンサンブルを探索後、当該探索結果に基づいてユークリット幾何符号の各行または各列からランダムに「１」を抽出し、各行または各列を分割することによって、Ｉｒｒｅｇｕｌａｒ−ＬＤＰＣ符号の検査行列を生成することを特徴とする請求項１に記載のＬＤＰＣ符号用検査行列生成方法。
- 前記アンサンブルの重み配分を、重み単位の重み総数が整数で、かつ重み単位の重み総数の総和とユークリット幾何符号の「１」の総数とが等しくなるように調整し、調整後のアンサンブルに基づいて分割処理を行うことを特徴とする請求項２に記載のＬＤＰＣ符号用検査行列生成方法。
- 基本のランダム系列のラテン方陣を作成し、
ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする請求項２または３に記載のＬＤＰＣ符号用検査行列生成方法。 - 基本のランダム系列のラテン方陣を作成し、
ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣の第ｍ行目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする請求項２または３に記載のＬＤＰＣ符号用検査行列生成方法。 - 基本のランダム系列のラテン方陣を複数個作成し、
列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ列目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ列目の「１」を抽出することを特徴とする請求項２または３に記載のＬＤＰＣ符号用検査行列生成方法。 - 基本のランダム系列のラテン方陣を複数個作成し、
列方向に連結したラテン方陣群行列を用いて、ユークリッド幾何符号の第ｍ行目をｎ分割する場合、上記ラテン方陣群行列の第ｍ列目のランダム系列をｎ分割し、当該ｎ分割後の各ランダム系列を用いてユークリッド幾何符号の第ｍ行目の「１」を抽出することを特徴とする請求項２または３に記載のＬＤＰＣ符号用検査行列生成方法。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001397922A JP3808769B2 (ja) | 2001-12-27 | 2001-12-27 | Ｌｄｐｃ符号用検査行列生成方法 |
| EP02791987A EP1460766A4 (en) | 2001-12-27 | 2002-12-25 | PRODUCTION PROCESS FOR AN LDPC CODE INSPECTION MATRIX |
| PCT/JP2002/013518 WO2003056705A1 (fr) | 2001-12-27 | 2002-12-25 | Procede de generation de matrice d'examen de code ldpc |
| US10/482,484 US7516388B2 (en) | 2001-12-27 | 2002-12-25 | LDPC code inspection matrix generation method |
| CNB028259076A CN100539442C (zh) | 2001-12-27 | 2002-12-25 | 低密度奇偶校验码用矩阵生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001397922A JP3808769B2 (ja) | 2001-12-27 | 2001-12-27 | Ｌｄｐｃ符号用検査行列生成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003198383A JP2003198383A (ja) | 2003-07-11 |
| JP3808769B2 true JP3808769B2 (ja) | 2006-08-16 |
Family
ID=19189258
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001397922A Expired - Fee Related JP3808769B2 (ja) | 2001-12-27 | 2001-12-27 | Ｌｄｐｃ符号用検査行列生成方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US7516388B2 (ja) |
| EP (1) | EP1460766A4 (ja) |
| JP (1) | JP3808769B2 (ja) |
| CN (1) | CN100539442C (ja) |
| WO (1) | WO2003056705A1 (ja) |
Families Citing this family (46)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3893383B2 (ja) * | 2002-02-28 | 2007-03-14 | 三菱電機株式会社 | Ｌｄｐｃ符号用検査行列生成方法および検査行列生成装置 |
| JP3946087B2 (ja) * | 2002-05-23 | 2007-07-18 | 三菱電機株式会社 | 通信システム、受信機および通信方法 |
| JP4163023B2 (ja) * | 2003-02-28 | 2008-10-08 | 三菱電機株式会社 | 検査行列生成方法および検査行列生成装置 |
| KR100922956B1 (ko) | 2003-10-14 | 2009-10-22 | 삼성전자주식회사 | 저밀도 패리티 검사 코드의 부호화 방법 |
| KR100632268B1 (ko) * | 2003-12-24 | 2006-10-11 | 한국전자통신연구원 | Ｌｄｐｃ 부호 부호화 및 복호화 방법, 및 ｌｄｐｃ패러티 검사 행렬 형성 방법. |
| WO2005069492A1 (ja) | 2004-01-20 | 2005-07-28 | Nec Corporation | 検査行列生成方法、データ伝送システム、符号化装置、復号装置および検査行列生成プログラム |
| CN1310433C (zh) * | 2004-03-29 | 2007-04-11 | 上海交通大学 | 采用分层低密度校验码的信道编码方法 |
| KR100540663B1 (ko) * | 2004-04-19 | 2006-01-10 | 삼성전자주식회사 | 저밀도 패리티 체크 코딩의 패리티 체크 행렬 변환 방법 |
| WO2005107081A1 (ja) * | 2004-04-28 | 2005-11-10 | Mitsubishi Denki Kabushiki Kaisha | 再送制御方法および通信装置 |
| JP4050726B2 (ja) * | 2004-06-23 | 2008-02-20 | 株式会社東芝 | 復号装置 |
| US7395490B2 (en) * | 2004-07-21 | 2008-07-01 | Qualcomm Incorporated | LDPC decoding methods and apparatus |
| KR100739684B1 (ko) | 2004-08-05 | 2007-07-13 | 삼성전자주식회사 | 저밀도 패리티 체크 행렬 생성 장치 및 방법 |
| JP4672016B2 (ja) * | 2004-08-09 | 2011-04-20 | エルジー エレクトロニクス インコーポレイティド | 低密度パリティ検査行列を用いた符号化及び復号化方法 |
| FI121431B (fi) | 2004-09-13 | 2010-11-15 | Tamfelt Pmc Oy | Paperikoneessa käytettävä kudosrakenne ja menetelmä sen valmistamiseksi |
| WO2006027818A1 (ja) * | 2004-09-03 | 2006-03-16 | Mitsubishi Denki Kabushiki Kaisha | Ｌｄｐｃ符号生成方法および通信装置 |
| WO2006040805A1 (ja) * | 2004-10-08 | 2006-04-20 | Mitsubishi Denki Kabushiki Kaisha | 検査行列生成装置及び通信装置 |
| CN100550655C (zh) * | 2004-11-04 | 2009-10-14 | 中兴通讯股份有限公司 | 一种低密度奇偶校验码的编码器/译码器及其生成方法 |
| CN100490332C (zh) * | 2004-11-05 | 2009-05-20 | 中国科学技术大学 | 一种基于循环三维立方体网格图的低密度码的构造方法 |
| JP2008529448A (ja) * | 2005-02-03 | 2008-07-31 | エージェンシー フォー サイエンス，テクノロジー アンド リサーチ | データの送信方法、データの受信方法、送信機、受信機、並びにコンピュータプログラム製品 |
| CN100414841C (zh) * | 2005-05-11 | 2008-08-27 | 电子科技大学 | 一种ldpc编码的实现方法 |
| EP1881610B1 (en) | 2005-05-13 | 2015-06-03 | NEC Corporation | Weighted LDPC decoding using the MaxLog approximation |
| CN100424999C (zh) * | 2005-06-06 | 2008-10-08 | 重庆市成鑫机电制造有限公司 | 逐条添加边算法的多进制低密度奇偶校验码编译码方法 |
| KR100804793B1 (ko) | 2005-10-07 | 2008-02-20 | 삼성전자주식회사 | 저밀도 패러티 검사 복호기에서의 검사 노드 갱신 방법 |
| CN100544212C (zh) * | 2006-01-23 | 2009-09-23 | 南京大学 | 高速的减少存储需求的低密度校验码解码器 |
| JP4662278B2 (ja) | 2006-04-28 | 2011-03-30 | 富士通株式会社 | エラー訂正装置、符号器、復号器、方法及び情報記憶装置 |
| US7685494B1 (en) | 2006-05-08 | 2010-03-23 | Marvell International, Ltd. | Error correction coding for varying signal-to-noise ratio channels |
| US8209580B1 (en) * | 2006-05-08 | 2012-06-26 | Marvell International Ltd. | Error correction coding for varying signal-to-noise ratio channels |
| WO2007145491A1 (en) * | 2006-06-15 | 2007-12-21 | Samsung Electronics Co., Ltd. | Apparatus and method of encoding/decoding block low density parity check codes in a communication system |
| CN101119178B (zh) * | 2006-08-01 | 2010-08-25 | 华为技术有限公司 | 信号发送、接收方法及信号发送、接收装置 |
| JP2008167435A (ja) * | 2006-12-27 | 2008-07-17 | Nec Lab America Inc | Ｂｃｊｒ等化器及び反復ｌｄｐｃ復号を用いた偏波モード分散補償 |
| US20090034661A1 (en) * | 2007-07-31 | 2009-02-05 | Pascal Olivier Vontobel | Coordinate-ascent method for linear programming decoding |
| WO2010001502A1 (ja) * | 2008-06-30 | 2010-01-07 | 国立大学法人愛媛大学 | 復号化装置、垂直磁気記録再生装置、受信装置、および、復号化方法 |
| JP2010128083A (ja) * | 2008-11-26 | 2010-06-10 | Oki Electric Ind Co Ltd | シンボル線形変換装置、共有鍵生成装置、シンボル線形変換方法、共有鍵生成方法、シンボル線形変換プログラム及び共有鍵生成プログラム |
| JP5197544B2 (ja) * | 2009-10-05 | 2013-05-15 | 株式会社東芝 | メモリシステム |
| JP5719998B2 (ja) * | 2010-02-10 | 2015-05-20 | パナソニックＩｐマネジメント株式会社 | 送信装置、受信装置、送信方法及び受信方法 |
| KR20110096684A (ko) * | 2010-02-23 | 2011-08-31 | 삼성전자주식회사 | 부가 정보의 피드백을 이용하여 통신하는 무선 네트워크 및 상기 무선 네트워크에서 네트워크 코딩을 이용한 통신 방법 |
| CN102122962B (zh) * | 2011-04-01 | 2013-04-17 | 山东大学 | 基于预测-修正原对偶内点法的ldpc码的lp译码器 |
| CN102158233A (zh) * | 2011-05-09 | 2011-08-17 | 山东大学 | 一种ldpc码的线性规划与最小和级联译码方法 |
| CN102148619B (zh) * | 2011-05-09 | 2013-06-05 | 山东大学 | 一种应用于ldpc码的自适应线性规划译码算法 |
| US8566666B2 (en) * | 2011-07-11 | 2013-10-22 | Lsi Corporation | Min-sum based non-binary LDPC decoder |
| CN110535474B (zh) | 2017-05-05 | 2023-06-06 | 华为技术有限公司 | 信息处理的方法、通信装置 |
| WO2018227681A1 (zh) | 2017-06-15 | 2018-12-20 | 华为技术有限公司 | 信息处理的方法和通信装置 |
| CN109150197B (zh) | 2017-06-27 | 2024-05-14 | 华为技术有限公司 | 信息处理的方法、装置和通信设备 |
| CN108462496B (zh) * | 2018-04-24 | 2021-04-02 | 成都吉纬科技有限公司 | 一种基于随机比特流更新的ldpc译码器 |
| CN109347485A (zh) * | 2018-09-29 | 2019-02-15 | 山东存储之翼电子科技有限公司 | 构造ldpc码校验矩阵的方法及ldpc码编译方法 |
| US20210279055A1 (en) * | 2020-03-03 | 2021-09-09 | Nvidia Corporation | Technique for performing bit-linear transformations |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US1526647A (en) * | 1921-12-10 | 1925-02-17 | Richard H Thomas | System for maintaining the vacuum in a heating plant |
| US6307487B1 (en) * | 1998-09-23 | 2001-10-23 | Digital Fountain, Inc. | Information additive code generator and decoder for communication systems |
| FR2799592B1 (fr) * | 1999-10-12 | 2003-09-26 | Thomson Csf | Procede de construction et de codage simple et systematique de codes ldpc |
| US6857097B2 (en) | 2001-05-16 | 2005-02-15 | Mitsubishi Electric Research Laboratories, Inc. | Evaluating and optimizing error-correcting codes using a renormalization group transformation |
| US6567465B2 (en) * | 2001-05-21 | 2003-05-20 | Pc Tel Inc. | DSL modem utilizing low density parity check codes |
| US6789227B2 (en) * | 2001-07-05 | 2004-09-07 | International Business Machines Corporation | System and method for generating low density parity check codes using bit-filling |
| US6757122B1 (en) * | 2002-01-29 | 2004-06-29 | Seagate Technology Llc | Method and decoding apparatus using linear code with parity check matrices composed from circulants |
| JP3893383B2 (ja) | 2002-02-28 | 2007-03-14 | 三菱電機株式会社 | Ｌｄｐｃ符号用検査行列生成方法および検査行列生成装置 |
| JP4005084B2 (ja) | 2002-07-02 | 2007-11-07 | 三菱電機株式会社 | 検査行列生成方法および検査行列生成装置 |
-
2001
- 2001-12-27 JP JP2001397922A patent/JP3808769B2/ja not_active Expired - Fee Related
-
2002
- 2002-12-25 EP EP02791987A patent/EP1460766A4/en not_active Ceased
- 2002-12-25 CN CNB028259076A patent/CN100539442C/zh not_active Expired - Fee Related
- 2002-12-25 WO PCT/JP2002/013518 patent/WO2003056705A1/ja active Application Filing
- 2002-12-25 US US10/482,484 patent/US7516388B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| EP1460766A4 (en) | 2007-03-14 |
| US7516388B2 (en) | 2009-04-07 |
| JP2003198383A (ja) | 2003-07-11 |
| CN100539442C (zh) | 2009-09-09 |
| EP1460766A1 (en) | 2004-09-22 |
| WO2003056705A1 (fr) | 2003-07-10 |
| CN1608347A (zh) | 2005-04-20 |
| US20040186992A1 (en) | 2004-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3808769B2 (ja) | Ｌｄｐｃ符号用検査行列生成方法 | |
| JP4163023B2 (ja) | 検査行列生成方法および検査行列生成装置 | |
| JP3893383B2 (ja) | Ｌｄｐｃ符号用検査行列生成方法および検査行列生成装置 | |
| JP4005084B2 (ja) | 検査行列生成方法および検査行列生成装置 | |
| US20110119554A1 (en) | Method for transmitting non-binary codes and decoding the same | |
| JP2008035524A (ja) | 反復確率伝搬を使用して記号ブロックを復号化する装置および方法 | |
| US20080184084A1 (en) | Check Matrix Generating Apparatus and Communication Apparatus | |
| Guruswami | Iterative decoding of low-density parity check codes (A Survey) | |
| US7406648B2 (en) | Methods for coding and decoding LDPC codes, and method for forming LDPC parity check matrix | |
| Lentmaier et al. | Exact erasure channel density evolution for protograph-based generalized LDPC codes | |
| Zikry et al. | Performance Analysis of LDPC Decoding Techniques | |
| Zhang et al. | Iterative algorithms for decoding a class of two-step majority-logic decodable cyclic codes | |
| Pal et al. | Error-Control Coding Algorithms and Architecture for Modern Applications Powered by LDPC Codes and Belief Propagation | |
| Moon et al. | An introduction to low-density parity-check codes | |
| Mitchell | Mathematical approach to channel codes with a diagonal matrix structure | |
| Castura | Performance analysis and optimization of reduced complexity Low Density Parity Check decoding algorithms | |
| Yazdani | Construction of irregular rate-compatible low-density parity-check codes for hybrid-ARQ schemes at short block lengths | |
| Albertazzi et al. | Forward error correction | |
| Soleymani et al. | Low Density Parity Check Codes | |
| Xie | Design of rate-compatible structured low density parity check codes | |
| Sun | Non-Binary Protograph-Based LDPC Codes: Analysis, Enumerators and Designs | |
| Vijayakumar | FPGA implementation of low density parity check codes decoder | |
| Kumar et al. | CONSTRUCTION OF REGULAR LDPC LIKE CODES BASED ON FULL RANK CODES AND THEIR ITERATIVE DECODING USING A PARITY CHECK TREE | |
| Uzunoğlu | Performance comparison of message passing decoding algorithms for binary and non-binary low density parity check (LDPC) codes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041028 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060516 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060518 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100526 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100526 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110526 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110526 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120526 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120526 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130526 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140526 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |