JP5489552B2 - 復号方法及び復号装置 - Google Patents

復号方法及び復号装置 Download PDF

Info

Publication number
JP5489552B2
JP5489552B2 JP2009146345A JP2009146345A JP5489552B2 JP 5489552 B2 JP5489552 B2 JP 5489552B2 JP 2009146345 A JP2009146345 A JP 2009146345A JP 2009146345 A JP2009146345 A JP 2009146345A JP 5489552 B2 JP5489552 B2 JP 5489552B2
Authority
JP
Japan
Prior art keywords
column
llr
row
value
minimum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2009146345A
Other languages
English (en)
Other versions
JP2011004229A (ja
Inventor
堅也 杉原
好邦 宮田
英夫 吉田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2009146345A priority Critical patent/JP5489552B2/ja
Priority to EP10006102.7A priority patent/EP2264906B1/en
Priority to US12/818,748 priority patent/US8601337B2/en
Priority to CN201010208249.3A priority patent/CN101931418B/zh
Publication of JP2011004229A publication Critical patent/JP2011004229A/ja
Application granted granted Critical
Publication of JP5489552B2 publication Critical patent/JP5489552B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, 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/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error 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/11Error 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/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1111Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1117Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using approximations for check node processing, e.g. an outgoing message is depending on the signs and the minimum over the magnitudes of all incoming messages according to the min-sum rule
    • H03M13/112Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using approximations for check node processing, e.g. an outgoing message is depending on the signs and the minimum over the magnitudes of all incoming messages according to the min-sum rule with correction functions for the min-sum rule, e.g. using an offset or a scaling factor
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, 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/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error 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/11Error 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/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1111Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1117Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using approximations for check node processing, e.g. an outgoing message is depending on the signs and the minimum over the magnitudes of all incoming messages according to the min-sum rule
    • H03M13/1122Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using approximations for check node processing, e.g. an outgoing message is depending on the signs and the minimum over the magnitudes of all incoming messages according to the min-sum rule storing only the first and second minimum values per check node
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, 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/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error 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/11Error 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/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1131Scheduling of bit node or check node processing
    • H03M13/114Shuffled, staggered, layered or turbo decoding schedules
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, 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/65Purpose and implementation aspects
    • H03M13/6577Representation or format of variables, register sizes or word-lengths and quantization
    • H03M13/6583Normalization other than scaling, e.g. by subtraction

Landscapes

  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)

Description

この発明は、誤り訂正符号の復号方法及び復号装置に関するものである。
低密度パリティ検査符号(Low-Density Parity-Check符号;以下、LDPC符号と称す)は、1の数が疎のm×n行列で定義される誤り訂正符号である。この行列を検査行列と呼ぶ。LDPC符号の復号法として、Sum−Product復号法が1960年代に考案されている。
しかしながら、Sum−Product復号法は、復号性能が高い反面、演算量が大きく、回路実装は困難である。これに対して、非特許文献1に記載のOffset BP−based復号法は、Sum−Product復号法を近似計算することにより、演算量を大幅に削減している。
先ず、Sum−Product復号法で行われる具体的な演算について説明する。
この復号法のアルゴリズムは、行演算と列演算と呼ばれる2つの演算からなり、これらの演算を交互に繰り返すことで復号を行う。
行演算は、検査行列の各行に対して行う演算である。
検査行列のr行で要素が1である列番号の集合をN(r)とおく。N(r)に含まれる各列cに対し、下記式(1)により行対数尤度比(Log-Likelihood Ratio;以下、LLRと称す)εr,cを計算する。なお、zr,cは、後に説明する列演算で計算される列LLRであるが、繰り返しの1回目は、受信系列からビットごとに計算されるLLRであるλcが格納される。c’はc以外のN(r)である。
Figure 0005489552
上記(1)式の各演算記号の定義は、下記式(2)〜(5)に示す通りである。
Figure 0005489552

Figure 0005489552


Figure 0005489552
Figure 0005489552



次に列演算について説明する。
行演算で求めた行LLRεr,cから下記式(6)を用い列LLRzr,cを計算する。なお、λcは、上述したLLRである。また、検査行列のc列で要素が1である行番号の集合をM(c)とおく。r’はr以外のM(c)である。
Figure 0005489552
また、送信ビットの推定を下記式(7)を用いて行う。下記式(7)による演算の結果として、zcが0以上であれば、受信系列のcビット目を1とし、zcが0より小さければcビット目を0と推定する。
Figure 0005489552
推定ビット列は、パリティチェックにより符号語であるか否かを判定できる。もし符号語でなければ、列演算で計算した列LLRを用いて行演算を行い、以降これまで説明した演算をやり直す。予め定めた繰り返し回数の上限に達したか、もしくは推定ビットが符号語であれば、推定ビット列を出力し復号は終了となる。Sum−Product復号法では、上記式(1)で表した行演算に膨大な演算量が必要となる。
非特許文献1のOffset BP−based復号法は、上記式(1)の代わりに、下記式(8)による行演算を行う。なお、オフセット値αは、検査行列ごとに求める定数であり、復号性能が最も高くなるαを、密度発展法(Density Evolution)を用いて算出する。Offset BP−based復号法では、最小値を求めてオフセット値αを引くだけでよく、行演算の演算量は大幅に削減される。
Figure 0005489552
Offset BP−based復号法は、Sum−Product復号法の近似演算であることは既に述べたが、次にそのことをより詳しく説明する。行演算中に現れる上記式(4)で定義される2項演算は、下記式(9)と同値であることが知られている。
Figure 0005489552
Offset BP−based復号法では、上記式(9)の第2項以降の対数関数を定数αで近似して得られる。なお、Sum−Product復号法の行演算で用いる上記式(1)には、再帰演算(つまり、上記式(3)で示す演算)が含まれ、この再帰演算による影響をも定数αに包含して近似したものが、Offset BP−based復号法の行演算に用いた上記式(8)である。
J. Chen, A. Dholakia, E. Eleftheriou, M.P.C. Fossorier, and X.-Y. Hu, "Reduced-Complexity Decoding of LDPC Codes", IEEE Trans. Communications, vol. 53, issue 8, pp.1288-1299, Aug. 2005.
上述したように、Offset BP−based復号法では、対数関数及び再帰演算による影響を定数αで近似している。
しかしながら、Sum−Product復号法による実際の数値は、行演算の変数によって大きく変動する。このため、Offset BP−based復号法では、近似精度が低く、復号性能が十分ではないという課題があった。
この発明は、上記のような課題を解決するためになされたもので、Offset BP−based復号法よりSum−Product復号法における行演算を高い精度で近似することができ、誤り訂正符号の復号性能を向上することができる復号方法及び復号装置を得ることを目的とする。
この発明に係る復号方法は、例えば、LDPC符号の検査行列の行重みに対応する行LLRを算出する行演算ステップと、行LLRから検査行列の列重みに対応する列LLRを算出する列演算ステップと、列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定ステップと、推定ビット列が符号語であるか否かを検査する検査ステップとを備え、行演算ステップと列演算ステップを交互に繰り返して符号語を復号するにあたり、行演算ステップにて、列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、列LLRの絶対値の最小値の大きさに応じて変動させたオフセット値を当該列LLRの絶対値の最小値から引いた値を、当該列LLRの列に対応する行LLRとすることを特徴とするものである。
この発明によれば、行演算において、例えば、列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、列LLRの絶対値の最小値の大きさに応じて変動させたオフセット値を当該列LLRの絶対値の最小値から引いた値を、当該列LLRの列に対応する行LLRとするので、Offset BP−based復号法よりも、Sum−Product復号法における行演算を高い精度で近似することができ、誤り訂正符号の復号性能を向上することができるという効果がある。
誤り訂正符号を用いたディジタル通信の手順を概略的に示す図である。 この発明の実施の形態1による誤り訂正復号方法の流れを示すフローチャートである。 g(x,δ)とxの関係を示すグラフである。 実施の形態1による行演算の一例を示すフローチャートである。 実施の形態1によるディジタル通信システムの構成を示すブロック図である。 図5の復号器の構成を示すブロック図である。 図6の行演算器の構成を示すブロック図である。 この発明の実施の形態3による行演算器の構成を示すブロック図である。 この発明の実施の形態4による行演算器の構成を示すブロック図である。
実施の形態1.
図1は、誤り訂正符号を用いたディジタル通信の手順を概略的に示す図である。図1に示すように、ディジタル通信では、先ず送信ビット列を符号化し、符号語とする(手順A)。次に、符号化したビット列を変調する(手順B)。送信された変調信号は、通信路を介して受信されると、受信側で復調されて受信信号値に変換される(手順C)。最後に、受信信号値に基づいて軟判定による復号を行い、推定ビット列を得る(手順D)。なお、この発明では、低密度パリティ検査符号(LDPC符号)を用いて符号化を行う。
図2は、この発明の実施の形態1による誤り訂正復号方法の流れを示すフローチャートである。先ず、図1に示す手順Cの復調によって得られた符号語長のビットのそれぞれの受信信号値から、ビットごとにLLRを計算する(ステップST1)。次に、ステップST1で得られたLLRを用いて行LLRを算出し(ステップST2)、この行LLRから列LLRを計算する(ステップST3)。
続いて、ステップST3で得られた列LLRを基にビットごとの推定を行い、推定ビット列を算出する(ステップST4)。この後、推定ビット列が符号語であるか否かを判定する(ステップST5)。このとき、符号語であれば(ステップST5;YES)、この推定ビット列を出力して(ステップST7)、処理を終了する。
また、推定ビット列が符号語ではないと判定すると(ステップST5;NO)、行演算以降の処理の繰り返し回数が最大繰り返し数に到達したか否かを判定する(ステップST6)。ここで、最大繰り返し数に到達している場合(ステップST6;YES)には、ステップST7に移行して推定ビット列を出力し、処理を終了する。
一方、最大繰り返し数に到達していなければ(ステップST6;NO)、ステップST2の行演算処理に戻り、これ以降の処理を繰り返す。ただし、2回目以降の行演算では、直前の列演算で算出した列LLRに基づいて行演算を行う。
次に、この実施の形態1による行演算の詳細を説明する。
上記式(4)は、x=min{x,y}と仮定し、δ=y−xとすれば、代入操作によって下記式(10)に変換することができる。
Figure 0005489552
さらに、上記式(10)の第2項及び第3項を、下記式(11)のようにg(x,δ)とする。
Figure 0005489552
図3は、上記式(11)で算出したg(x,δ)とxの関係を示すグラフであり、横軸がxで縦軸がg(x,δ)である。図3において、δは固定値とし、δ=2,1,0.5,0.1,0のそれぞれの場合で、上記式(11)を用いてg(x,δ)を算出した結果を図示している。上述のようにOffset BP−based復号法は、g(x,δ)を定数で近似する復号法であるが、図3に示すグラフから明らかなように、xが0に近い場合には、δの値に関わらず、g(x,δ)も0に近い値となる。
そこで、この実施の形態1では、xがある定数γより小さい場合にはg(x,δ)=0と近似し、xがγ以上の場合には、Offset BP−based復号法と同様にg(x,δ)を定数で近似する。以上では、上記式(4)における近似についての説明であるが、上記式(1)の再帰演算(つまり、上記式(3)の演算)の影響を考慮に入れれば、|zr,c'|の最小値がγ以上のときのみ|zr,c'|の最小値から定数のオフセット値βを引く近似が可能である。これを定式化して表したものが、下記式(12)である。下記式(12)中のsgr,cは、下記式(13)で表される。
なお、βとγは、検査行列により異なる値をとる定数である。また、これらの定数は、密度発展法や復号シミュレーション、もしくはその他の手段によって復号性能が最大となる値を見出し、その値に予め設定する。
Figure 0005489552

Figure 0005489552
次に、上記式(12)の行演算について説明する。
図4は、実施の形態1による行演算の一例を示すフローチャートであり、上記式(12)に従った行演算を示している。なお、行演算は、各行それぞれに対して行うが、それを表したのが外側の行ループである。また、検査行列の行列要素が1となる列のそれぞれに対して、行演算では行LLRを計算する。これを表しているのが内側の列ループである。ここでは、r行の行演算における列cに対する行LLRの演算方法を示す。
先ず、行rの行列要素が1の列(列c以外)から列LLRの絶対値の最小値を計算する(ステップST1a)。この値を最小値minとおく。続いて、最小値minがγ以上であるか否かを判定する(ステップST2a)。最小値minがγ以上であると判定された場合(ステップST2a;YES)は、ステップST3aの処理に移行して、当該最小値minから定数βを引いた(min−β)が0より大きいか否かを判定する。
(min−β)が0より大きい場合(ステップST3a;YES)、上記式(13)に従って、sgr,cを計算する(ステップST4a)。この後、ステップST4aで算出したsgr,cに(min−β)を掛け算した値であるsgr,c・(min−β)を出力する(ステップST5a)。
一方、最小値minがγに満たない場合(ステップST2a;NO)は、ステップST4aと同様に、上記式(13)に従って、sgr,cを計算する(ステップST6a)。この後、sgr,cに最小値minを掛け算した値であるsgr,c・minを出力する(ステップST7a)。
また、(min−β)が0以下である場合(ステップST3a;NO)、0を出力する(ステップST8a)。ここまでが、実施の形態1における行演算であるが、行演算結果が上記式(12)と等価であれば、上記演算に限定されるものではない。
例えば、ステップST4aとステップST6aを共通化して、ステップST2aの前段で処理を行ってもよい。また、γ≧βであれば、ステップST3a及びステップST8aは必要なく、行わなくてもよい。
このようにすることで、Sum−Product復号法における行演算を高い精度で近似でき、誤り訂正符号の復号性能を向上できる。また、演算量の小さいOffset BP−based復号法と比べても比較演算が1回増加しただけのほぼ同じ演算量であり、Sum−Product復号法に比べれば演算量を格段に削減できる。
次に、上述の誤り訂正復号方法を利用した誤り訂正符号化及び復号を実行するシステム及び当該復号を実行する復号器について説明する。
図5は、実施の形態1によるディジタル通信システムの構成を示すブロック図である。図5において、実施の形態1によるディジタル通信システム1は、送信側の構成として、符号化器2及び変調器3を備え、受信側の構成として、復調器5及び復号器6を備える。符号化器2は、送信ビット列を符号化し符号語とする。変調器3は、符号化器2で符号化されたビット列を変調信号に変換して通信路4へ送信する。受信側では、復調器5が、通信路4を介して受信された信号を復調し、受信信号値を得る。復号器6では、復調器5によって得られた符号語長のビットのそれぞれの受信信号値から推定ビット列を推定し出力する。
図6は、図5の復号器の構成を示すブロック図である。図6において、復号器6は、入力回路7、LLR演算器8、行演算器(行演算手段)9、列演算器(列演算手段)10、ビット推定器(ビット推定手段)11、パリティ検査器(検査手段)12、繰り返し回数判定回路13及び出力回路14を備える。各構成は、図6に示すように直列に接続されており、さらにパリティ検査器12と出力回路14との接続、及び繰り返し回数判定回路13と行演算器9との接続もなされている。
なお、図6の構成は一例であり、他の構成も可能である。この発明では、行演算器9の内部構成に特徴があり、その他の構成を制約するものではない。例えば、列演算器10とビット推定器11(及びパリティ検査器12)は並列に接続してもよいし、繰り返し回数判定回路13をパリティ検査器12の前段に配置しても構わない。
また、復号器6は、図2に示したフローチャートの復号処理を実行する。
復号器6では、復調器5による復調処理で得られた受信信号値を入力すると、入力回路7が、入力した受信信号値を符号語ごとの固まりに分け、適切なタイミングでLLR演算器8へ出力する。LLR演算器8は、符号語長のビットごとの受信信号値から各LLRを計算する(図2;ステップST1)。
行演算器9は、LLR演算器8で算出されたLLRを用いて行LLRを計算する(図2;ステップST2)。ただし、繰り返し回数判定回路13から処理が渡された場合は、列演算器10が計算した直近の列LLRを用いて行LLRを計算する。列演算器10は、行演算器9が算出した行LLRから列LLRを計算する(図2;ステップST3)。
ビット推定器11は、列演算器10が算出した列LLRから推定ビット列を計算する(図2;ステップST4)。パリティ検査器12は、ビット推定器11で算出された推定ビット列のパリティチェックを行い(図2;ステップST5)、推定ビット列が符号語であれば、出力回路14に推定ビット列を送る。符号語でなければ、その旨を繰り返し回数判定回路13へ通知する。
繰り返し回数判定回路13は、行演算(もしくは列演算)を行った回数を符号語ごとにカウントしており、パリティ検査器12から上記通知を受けると、予め定めた繰り返し回数に到達しているか否かを判定する(図2;ステップST6)。ここで、予め定めた繰り返し回数に到達していれば、繰り返し回数判定回路13は、出力回路14へ推定ビット列を送る。到達していなければ、行演算器9へ処理を渡す。出力回路14は、上述のようにして得られた推定ビット列を復号器6の出力として出力する(図2;ステップST7)。
図7は、図6の行演算器の構成を示すブロック図である。図7において、実施の形態1による行演算器9は、入力回路15、最小2値演算部16、minのオフセット処理部17、2nd minのオフセット処理部18、出力行LLR選択部19、sgr,c演算部20、掛け算回路21及び出力回路22を備える。入力回路15は、行演算器9への入力データが最初に通る回路である。また、出力回路22は、行演算器9の出力に繋がっている。
また、行演算器9は、図4のフローチャートに相当する行演算を実行する。
図7の構成では、行演算を1行ずつ行い、この行から各列に対して出力する行LLRを一度に演算する。
先ず、入力回路15は、検査行列の対応する行において行列要素が1となる全ての列の列LLRを出力する。
最小2値演算部16は、入力回路15から入力した列LLRの絶対値から最小値minとその列LLRの列番号c’、及び2番目に小さい値である2nd minを計算して、出力する(図4;ステップST1a)。
minのオフセット処理部17では、最小値minがγ以上であるか否かを判定(図4;ステップST2a)し、最小値minがγ未満であれば、最小値minをそのまま出力行LLR選択部19へ出力する。また、最小値minがγ以上であると、minのオフセット処理部17は、(min−β)が0よりも大きいか否かを判定する(図4;ステップST3a)。(min−β)が0よりも大きい場合、minのオフセット処理部17は、(min−β)を出力行LLR選択部19へ出力する。一方、(min−β)が0以下であれば、出力行LLR選択部19へ0を出力する。
なお、2nd minのオフセット処理部18においても、minのオフセット処理部17と同様の動作をする。
出力行LLR選択部19では、列c’に対する行LLRについては、2nd minのオフセット処理部18によるオフセット処理後の値(2nd min−β)、それ以外の列に対しては、minのオフセット処理部17によるオフセット処理後の値(min−β)が出力されるよう出力値を選択する。
sgr,c演算部20は、入力回路15から列LLRの符号を入力し、上記式(13)を用いて、sgr,cの演算をそれぞれの列に対して行う(図4;ステップST4a、ステップST6a)。
掛け算回路21は、各列に対して行LLRの絶対値とsgr,c(=+1又は−1)との掛け算を行い(図4;ステップST5a、ステップST7a)、行LLRの正負を決定する。出力回路22は、掛け算回路21による演算結果の行LLRの出力制御を行う。
なお、この実施の形態1による行演算器9は、図7に示した回路構成に限定されるものではない。例えば、図4のフローチャートに忠実に従って行演算を行う回路としてもよいし、演算内容が等価であれば他の手段による回路であってもよい。演算結果さえ同じであれば、等価な復号性能となる。
以上のように、この実施の形態1によれば、図7に示すような回路構成によって、図4の手順で行演算を行うことにより、Sum−Product復号法における行演算を高い精度で近似することができ、誤り訂正符号の復号性能を向上させることができる。
また、Offset BP−based復号法と比べて比較演算が1回増加しただけのほぼ同じ演算量であることから、回路規模の小さいOffset BP−based復号法と同等の回路規模で実現できる。このため、Sum−Product復号法に比べて、演算量を格段に削減できる上、回路規模の小規模化も実現することができる。
なお、上記実施の形態1では、γによって2通りに場合分けしたが、より多くの場合分けをすることも可能である。つまり、|zr,c'|の最小値の大きさによってオフセット値を変動させ、|zr,c'|の最小値からそれぞれの場合に応じたオフセット値を引くという方法が可能である。場合分けを多くすることによって演算量は大きくなるが、より近似精度が高くなり復号性能は向上する。
実施の形態2.
上記実施の形態1で近似を行った関数g(x,δ)は、上記式(11)の通り、xとδを変数としている。この関数g(x,δ)において、上記実施の形態1では、変数xの値の変動を考慮して近似を行った。
しかしながら、上記実施の形態1で示した図3のグラフから明らかなように、変数xがある程度大きい領域では関数g(x,δ)の値が変数δの影響によって大きく変動する。上記実施の形態1は、この変動を無視してオフセット値βを定数として近似している。
そこで、この実施の形態2では、変数δの値によってオフセット値を変動させる方法について述べる。なお、変数δによってオフセット値を変動させる技術としては、例えば、参考文献1に記載のδ−min復号法がある。
(参考文献1)阪井塁、松本渉、吉田英夫、“更新式の近似に基づくLow−Density Parity−Check 符号簡易復号法の検討”,電子情報通信学会論文誌A,vol.J90−A,no.2,pp.83−91,Feb.2007.
参考文献1に記載のδ−min復号法は、変数xによる関数g(x,δ)の変動を無視し、変数δのみからオフセット値を決定する手法であるため、図3のグラフで示したように変数xが0に近い値のときの近似精度が低い。
そこで、この実施の形態2では、xの値がγ未満ときは、g(x,δ)=0と近似し、γ以上の場合は、δ−min復号法で用いられているオフセット値D(δ)を用いる。つまり、下記式(14)のようにする。なお、D(δ)は、下記式(15)に示す関数である。変数δはδ=y−xである。
Figure 0005489552

Figure 0005489552
以上のように、この実施の形態2によれば、上記式(14)で上記式(4)の2項演算を置き換え、g(x,δ)を上述のように近似することで、上記実施の形態1における行演算よりも近似精度が向上し、復号性能を高くすることができる。
実施の形態3.
上記実施の形態2は関数g(x,δ)のみを近似したが、この実施の形態3では、上記実施の形態1のように、行演算の再帰演算(上記式(3)の演算)も近似する。つまり、上記実施の形態1の行演算(上記式(12)の演算)のオフセット値βを、上記実施の形態2におけるD(δ)で置き換えた下記式(16)を利用する。ただし、δは|zr,c'|の最小値と2番目に小さい値との差とする。このように行演算を近似することによって、行演算の近似精度と復号性能が上記実施の形態1よりも高く、演算量が上記実施の形態2よりも小さくなる。
Figure 0005489552
図8は、この発明の実施の形態3による行演算器の構成を示すブロック図であり、既に説明した図6における行演算器の内部を実施の形態3に合わせて図示したものである。図8において、実施の形態3による行演算器(行演算手段)9Aは、入力回路15、sgr,c演算部20、掛け算回路21、出力回路22、最小3値演算部23、オフセット処理部24及び出力行LLR選択部25を備える。入力回路15、sgr,c演算部20、掛け算回路21及び出力回路22の動作は、上記実施の形態1と同じであるので、異なる部分について説明する。
最小3値演算部23は、列LLRの絶対値のうち、最小値minとその列番号c’、2番目に小さい値2nd minとその列番号c’’、そして3番目に小さい値3rd minを演算して出力する。
オフセット処理部24は、上記式(16)に従ってオフセット値を計算してオフセット処理を行うが、次の3つを出力する。
1つは、minと2nd minから計算するminのオフセット処理後の値であり、minがγ未満ならminを出力し、minがγ以上ならδ=2nd min−minとして、min−D(δ)もしくは0を出力する。
残りの2つの出力も同様であり、minと3rd minから計算するminのオフセット処理後の値、2nd minと3rd minから計算する2nd minのオフセット処理後の値をそれぞれ計算し出力する。
出力行LLR選択部25は、列c’に対する行LLRとして2nd minと3rd minとから計算した2nd minのオフセット処理後の値を選択する。また、列c’’に対する行LLRとしては、3rd minとminとから計算したminのオフセット処理後の値を選択し、それ以外の列に対しては、2nd minとminとから計算したminのオフセット処理後の値を選択し、それぞれを出力する。
以上のように、この実施の形態3によれば、図8に示すような回路構成で、上記式(16)を利用した行演算を行うことにより、上記実施の形態1よりも、Sum−Product復号法における行演算を高い精度で近似することができ、誤り訂正符号の復号性能を向上させることができる。
実施の形態4.
上記実施の形態3は復号能力が高いものの、最小3値を演算する必要があり、回路規模が大きくなる。そこで、この実施の形態4では、上記実施の形態1と同様に、最小2値の演算を行いながら、上記実施の形態3と同様のδによるオフセット処理を部分的に行う構成について述べる。
図9は、この発明の実施の形態4による行演算器の構成を示すブロック図である。図9において、実施の形態4による行演算器(行演算手段)9Bは、入力回路15、sgr,c演算部20、掛け算回路21、出力回路22、最小2値演算部16a、オフセット処理部24a及び出力行LLR選択部25aを備える。入力回路15、sgr,c演算部20、掛け算回路21及び出力回路22の動作は、上記実施の形態1及び上記実施の形態3と同じである。最小2値演算部16a、オフセット処理部24a及び出力行LLR選択部25aの動作が、上記実施の形態3と異なる。
最小2値演算部16aは、上記実施の形態1で示した最小2値演算部16と一部異なり、最小値minとその列番号c’、2番目に小さい値2nd minに加えて、2nd minの列番号c’’も出力する。
また、オフセット処理部24aは、3つの出力を行う。
1つは、minに対して、オフセット値を定数βとして上記実施の形態1と同様の処理を行って得られた値(min−β)である。
2つ目は、2nd minに対しても同様に、オフセット値を定数βとして、上記実施の形態1と同様の処理を行って得られた値(2nd min−β)である。
3つ目は、上記実施の形態3と同様の処理で、δ=2nd min−minとし、オフセット値をD(δ)としてオフセット処理を行って得られた値(min−D(δ))である。
出力行LLR選択部25aは、列c’に対しては、2nd minにβによってオフセット処理をした値(2nd min−β)を選択する。また、列c’’に対しては、minにβによってオフセット処理をした値(min−β)を選択し、それ以外の列に対しては、minにD(δ)でオフセット処理をした値(min−D(δ))を選択する。出力行LLR選択部25aの後段における処理は、上記実施の形態1と同様である。
以上のように、この実施の形態4によれば、図9に示すような回路構成で、上記の行演算を行うことで、上記実施の形態1による行演算よりも、列c’、列c’’以外に対する近似精度が高くなり、復号性能を向上できる。また、最小3値の演算を最小2値の演算に抑えることができることから、上記実施の形態3よりも回路規模を小さくできる。
1 ディジタル通信システム、2 符号化器、3 変調器、4 通信路、5 復調器、6 復号器、7 入力回路、8 LLR演算器、9,9A,9B 行演算器(行演算手段)、10 列演算器(列演算手段)、11 ビット推定器(ビット推定手段)、12 パリティ検査器(検査手段)、13 繰り返し回数判定回路、14 出力回路、15 入力回路、16,16a 最小2値演算部、17 minのオフセット処理部、18 2nd minのオフセット処理部、19,25,25a 出力行LLR選択部、20 sgr,c演算部、21 掛け算回路、22 出力回路、23 最小3値演算部、24,24a オフセット処理部。

Claims (8)

  1. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号方法において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算ステップと、
    前記行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算ステップと、
    前記列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定ステップと、
    推定ビット列が前記符号語であるか否かを検査する検査ステップとを備え、
    前記行演算ステップと前記列演算ステップを交互に繰り返して前記符号語を復号するにあたり、
    前記行演算ステップにて、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、前記列LLRの絶対値の最小値の大きさに応じて変動させたオフセット値を当該列LLRの絶対値の最小値から引いた値を、当該列LLRの列に対応する行LLRとすることを特徴とする復号方法。
  2. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号方法において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算ステップと、
    前記行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算ステップと、
    前記列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定ステップと、
    推定ビット列が前記符号語であるか否かを検査する検査ステップとを備え、
    Sum−Product復号法のアルゴリズムを用いて、前記行演算ステップと前記列演算ステップを交互に繰り返して前記符号語を復号するにあたり、
    前記行演算ステップにて、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、前記Sum−Product復号法のアルゴリズムの行演算で実行される2項演算の対象となる2値の最小値の大きさに応じて、前記2値の差を変数δとするδ−min復号法で規定されるオフセット値を変動させ、当該オフセット値を前記2値の最小値から引いた差分値を前記2項演算の結果とすることを特徴とする復号方法。
  3. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号方法において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算ステップと、
    前記行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算ステップと、
    前記列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定ステップと、
    推定ビット列が前記符号語であるか否かを検査する検査ステップとを備え、
    前記行演算ステップと前記列演算ステップを交互に繰り返して前記符号語を復号するにあたり、
    前記行演算ステップにて、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、前記列LLRの絶対値の最小値、2番目に小さい値、及び3番目に小さい値である最小3値の大きさに応じて、前記最小3値のうちの2値の差を変数δとするδ−min復号法で規定されるオフセット値を変動させ、当該オフセット値を前記列LLRの絶対値の最小値から引いた差分値を、当該列LLRの列に対応する行LLRとすることを特徴とする復号方法。
  4. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号方法において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算ステップと、
    前記行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算ステップと、
    前記列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定ステップと、
    推定ビット列が前記符号語であるか否かを検査する検査ステップとを備え、
    前記行演算ステップと前記列演算ステップを交互に繰り返して前記符号語を復号するにあたり、
    前記行演算ステップにて、
    絶対値が最小となる列LLRに対しては、当該最小値の大きさに応じたオフセット値を当該列LLRの絶対値の2番目に小さい値から引いた差分値を、当該列LLRの列に対応する行LLRとし、
    絶対値が2番目に小さい列LLRに対しては、当該列LLRの絶対値の大きさに応じたオフセット値を当該列LLRの絶対値の最小値から引いた差分値を、当該列LLRの列に対応する行LLRとし、
    上記以外の列に対しては、前記列LLRの絶対値の2番目に小さい値から前記列LLRの絶対値の最小値を引いた差を変数δとするδ−min復号法で規定されるオフセット値を当該列LLRの絶対値の最小値から引いた差分値を、当該列に対応する行LLRとすることを特徴とする復号方法。
  5. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号装置において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算手段と、
    前記行演算手段で算出された行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算手段と、
    前記列演算手段で算出された列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定手段と、
    前記ビット推定手段で算出された推定ビット列が前記符号語であるか否かを検査する検査手段とを備え、
    前記行演算手段は、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、前記列LLRの絶対値の最小値の大きさに応じて変動させたオフセット値を当該列LLRの絶対値の最小値から引いた値を、当該列LLRの列に対応する行LLRとすることを特徴とする復号装置。
  6. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号装置において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算手段と、
    前記行演算手段で算出された行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算手段と、
    前記列演算手段で算出された列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定手段と、
    前記ビット推定手段で算出された推定ビット列が前記符号語であるか否かを検査する検査手段とを備え、
    前記行演算手段は、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、Sum−Product復号法のアルゴリズムの行演算で実行される2項演算の対象となる2値の最小値の大きさに応じて、前記2値の差を変数δとするδ−min復号法で規定されるオフセット値を変動させ、当該オフセット値を前記2値の最小値から引いた差分値を前記2項演算の結果とすることを特徴とする復号装置。
  7. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号装置において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算手段と、
    前記行演算手段で算出された行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算手段と、
    前記列演算手段で算出された列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定手段と、
    前記ビット推定手段で算出された推定ビット列が前記符号語であるか否かを検査する検査手段とを備え、
    前記行演算手段は、
    前記列LLRの絶対値の最小値と予め設定された値との大きさ関係を判定し、当該最小値が予め設定された値以上である場合、前記列LLRの絶対値の最小値、2番目に小さい値、及び3番目に小さい値である最小3値の大きさに応じて、前記最小3値のうちの2値の差を変数δとするδ−min復号法で規定されるオフセット値を変動させ、当該オフセット値を前記列LLRの絶対値の最小値から引いた差分値を、当該列LLRの列に対応する行LLRとすることを特徴とする復号装置。
  8. LDPC(Low-Density Parity-Check)符号で符号化された符号語を復号する復号装置において、
    前記LDPC符号の検査行列の行重みに対応する行LLR(Log-Likelihood Ratio;対数尤度比)を算出する行演算手段と、
    前記行演算手段で算出された行LLRから前記検査行列の列重みに対応する列LLRを算出する列演算手段と、
    前記列演算手段で算出された列LLRを基にビットごとの推定を行って推定ビット列を算出するビット推定手段と、
    前記ビット推定手段で算出された推定ビット列が前記符号語であるか否かを検査する検査手段とを備え、
    前記行演算手段は、
    絶対値が最小となる列LLRに対しては、当該最小値の大きさに応じたオフセット値を当該列LLRの絶対値の2番目に小さい値から引いた差分値を、当該列LLRの列に対応する行LLRとし、
    絶対値が2番目に小さい列LLRに対しては、当該列LLRの絶対値の大きさに応じたオフセット値を当該列LLRの絶対値の最小値から引いた差分値を、当該列LLRの列に対応する行LLRとし、
    上記以外の列に対しては、前記列LLRの絶対値の2番目に小さい値から前記列LLRの絶対値の最小値を引いた差を変数δとするδ−min復号法で規定されるオフセット値を当該列LLRの絶対値の最小値から引いた差分値を、当該列に対応する行LLRとすることを特徴とする復号装置。
JP2009146345A 2009-06-19 2009-06-19 復号方法及び復号装置 Active JP5489552B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2009146345A JP5489552B2 (ja) 2009-06-19 2009-06-19 復号方法及び復号装置
EP10006102.7A EP2264906B1 (en) 2009-06-19 2010-06-11 Modified Offset BP-based LDPC decoding
US12/818,748 US8601337B2 (en) 2009-06-19 2010-06-18 Decoding method and decoding device
CN201010208249.3A CN101931418B (zh) 2009-06-19 2010-06-18 解码方法以及解码装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009146345A JP5489552B2 (ja) 2009-06-19 2009-06-19 復号方法及び復号装置

Publications (2)

Publication Number Publication Date
JP2011004229A JP2011004229A (ja) 2011-01-06
JP5489552B2 true JP5489552B2 (ja) 2014-05-14

Family

ID=42341634

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009146345A Active JP5489552B2 (ja) 2009-06-19 2009-06-19 復号方法及び復号装置

Country Status (4)

Country Link
US (1) US8601337B2 (ja)
EP (1) EP2264906B1 (ja)
JP (1) JP5489552B2 (ja)
CN (1) CN101931418B (ja)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5634888B2 (ja) 2011-01-12 2014-12-03 富士フイルム株式会社 インク組成物、画像形成方法及び印画物
CN103973315B (zh) * 2013-01-25 2019-01-18 中兴通讯股份有限公司 一种低密度奇偶校验码译码装置及其译码方法
CN103259545B (zh) * 2013-04-26 2017-02-15 西安理工大学 基于振荡的准循环低密度奇偶校验码置信传播译码方法
JP5971670B2 (ja) * 2013-12-09 2016-08-17 三菱電機株式会社 誤り訂正復号装置
EP2892157A1 (en) * 2014-01-02 2015-07-08 Alcatel Lucent Offset Min-Sum decoding of LDPC codes
CN104052499B (zh) * 2014-06-04 2017-04-19 华中科技大学 一种ldpc码的纠删译码方法及系统
JP6511284B2 (ja) * 2015-02-13 2019-05-15 パナソニック株式会社 最小値選択回路、復号器及び最小値選択方法
US9749089B2 (en) * 2015-11-04 2017-08-29 Mitsubishi Electric Research Laboratories, Inc. Fast log-likelihood ratio (LLR) computation for decoding high-order and high-dimensional modulation schemes
CN106708649B (zh) * 2015-11-13 2019-08-06 深圳大心电子科技有限公司 解码方法、存储器储存装置及存储器控制电路单元
CN108933603B (zh) * 2018-02-28 2022-02-22 和芯星通科技(北京)有限公司 一种实现校验节点处理的方法及装置
US10707901B2 (en) * 2019-01-08 2020-07-07 Intel Corporation Die-wise residual bit error rate (RBER) estimation for memories
US12040896B2 (en) 2020-04-02 2024-07-16 Samsung Electronics Co., Ltd. Methods and systems for optimizing computation of log-likelihood ratio (LLR) for decoding modulated symbols
CN114665923B (zh) * 2022-03-09 2024-05-07 深圳市国电科技通信有限公司 用于载波通信的数据处理方法、装置及处理器
FR3135581B1 (fr) * 2022-05-12 2024-05-24 Airbus Defence & Space Sas Optimisation du calcul des messages de contrôle parité dans un processus de décodage LDPC

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1526647B1 (en) * 2002-07-02 2008-10-01 Mitsubishi Electric Corporation Generation of a check matrix for irregular low-density parity-check (LDPC) codes
KR100924189B1 (ko) * 2004-12-02 2009-10-29 미쓰비시덴키 가부시키가이샤 복호 장치 및 통신 장치
US7562279B2 (en) 2005-05-20 2009-07-14 Mitsubishi Electric Research Laboratories, Inc. 2D-normalized min-sum decoding for ECC codes
KR100983692B1 (ko) * 2005-07-13 2010-09-24 미쓰비시덴키 가부시키가이샤 통신 장치 및 복호 방법
JP4341639B2 (ja) * 2006-05-15 2009-10-07 住友電気工業株式会社 復号装置および復号プログラム
CN101267209B (zh) * 2008-04-29 2010-09-08 清华大学 Ldpc译码的循环式分级最小值计算方法及其实现装置
CN101803205B (zh) * 2008-08-15 2013-12-18 Lsi公司 近码字的ram列表解码

Also Published As

Publication number Publication date
EP2264906A3 (en) 2012-11-07
CN101931418B (zh) 2014-03-26
US8601337B2 (en) 2013-12-03
CN101931418A (zh) 2010-12-29
JP2011004229A (ja) 2011-01-06
US20100325514A1 (en) 2010-12-23
EP2264906A2 (en) 2010-12-22
EP2264906B1 (en) 2015-09-09

Similar Documents

Publication Publication Date Title
JP5489552B2 (ja) 復号方法及び復号装置
JP4777261B2 (ja) 復号装置及び通信装置
Sarkis et al. Increasing the throughput of polar decoders
KR102274328B1 (ko) 체크 노드 계산유닛을 관리하는 방법, 장치, 및 상기 방법을 실행하기 위한 소프트웨어
CN108183713B (zh) 基于改进型最小和算法的ldpc译码器及其译码方法
US10498363B2 (en) Low density parity check decoder using binary logarithm and decoding method thereof
KR20060044395A (ko) 저밀도 패리티검사부호화방식에 따라서 실현되는 복호장치및 전처리장치
JP2012151676A (ja) 復号装置および復号方法
KR20070039353A (ko) 저밀도 패러티 검사 복호기에서의 검사 노드 갱신 방법
KR101583139B1 (ko) 높은 처리량과 낮은 복잡성을 갖는 연속 제거 극 부호 복호 장치 및 그 방법
JP2018019401A (ja) リード・ソロモン復号器及び復号方法
KR101882620B1 (ko) 신드롬 기반 연판정 bch 복호 장치 및 그 방법
CN106301388A (zh) 多进制ldpc码译码方法
KR101600759B1 (ko) 연속 제거 극 부호 복호기를 위한 간단한 병합 처리 연산 방법 및 장치
KR20150031568A (ko) 디지털 비디오 방송 시스템에서 LDPC(Low Density Parity Check) 복호기 및 LDPC 복호기의 복호화 방법
CN111164897B (zh) 广义低密度奇偶校验码
US9531577B2 (en) Bit-likelihood calculating apparatus and bit-likelihood calculating method
TWI487291B (zh) 循環碼解碼器及其方法
JP5556570B2 (ja) 信号処理回路及び受信装置
US8607132B1 (en) Matching signal dynamic range for turbo equalization system
JP5653536B2 (ja) Llr算出器及び誤り訂正復号装置
JP6177141B2 (ja) 対数尤度比算出装置、対数尤度比算出方法、及び対数尤度比算出用プログラム
RU2522299C1 (ru) Способ и устройство помехоустойчивого декодирования сигналов, полученных с использованием кода проверки на четность с низкой плотностью
JP4973647B2 (ja) 誤り訂正符号の復号評価装置
JP2011160491A (ja) 復号器

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120119

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130516

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130806

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131004

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140225

R150 Certificate of patent or registration of utility model

Ref document number: 5489552

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250