(57)【要約】【課題】受信信号系列から送信シンボルの判定を行う際にシンボル多値度が増加した場合であっても、演算量の増大を抑制すること。【解決手段】伝送路から取り込む受信信号系列に対して伝送路の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、生成したシンボル系列に対して仮判定を行い、仮判定した仮判定シンボルと仮判定シンボルの近傍のシンボルの範囲で伝送路状態を示す複数のシンボル系列を生成し、生成した伝送路状態を示す複数のシンボル系列と、伝送路の推定伝達関数とに基づいて、伝送路状態ごとの推定受信シンボル系列を生成し、受信信号系列から得られるシンボル系列と、推定受信シンボル系列の各々とのメトリックを算出し、算出したメトリックと、仮判定シンボルと、仮判定シンボルの近傍のシンボルとに基づいて最も尤もらしい推定受信シンボル系列を選択して送信シンボル系列の判定を行う。【選択図】図２
2020184696.tif 000002
【特許請求の範囲】
【請求項１】
伝送路から取り込む受信信号系列に対して前記伝送路の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、前記シンボル系列に対して仮判定を行う仮判定部と、
伝送路状態を示す複数のシンボル系列と、前記伝送路の推定伝達関数とに基づいて、前記伝送路状態ごとの推定受信シンボル系列を生成する伝送路推定部と、
前記受信信号系列から得られるシンボル系列と、前記推定受信シンボル系列の各々とのメトリックを算出し、算出した前記メトリックと、前記仮判定部が仮判定した仮判定シンボルと、前記仮判定シンボルの近傍のシンボルとに基づいて所定の推定アルゴリズムにより最も尤もらしい前記推定受信シンボル系列を選択する系列推定アルゴリズム処理部と、
前記最も尤もらしい推定受信シンボル系列に基づいて、トレリスのパスを遡って送信シンボル系列の判定を行うパス遡り判定部と、を備え、
前記系列推定アルゴリズム処理部は、前記仮判定シンボルと前記仮判定シンボルの近傍のシンボルの範囲で前記伝送路状態を示す複数のシンボル系列を生成して前記伝送路推定部に出力するシンボル判定装置。
【請求項２】
前記近傍のシンボルとは、前記仮判定シンボルを中心とした前後のシンボルであり、前記仮判定シンボルと前記近傍のシンボルの総和が（２ｎ＋１）個（ｎは、２以上の整数）になるシンボルである、請求項１に記載のシンボル判定装置。
【請求項３】
前記受信信号系列から得られるシンボル系列は、前記受信信号系列に含まれるシンボル系列であるか、または、前記仮判定部が前記受信信号系列の一部に対して前記推定逆伝達関数による適応等化を行って生成したシンボル系列であるか、または、更に備える適応フィルタ部によって生成されるシンボル系列であり、
前記適応フィルタ部は、入力信号系列の一部に対してパルス幅を圧縮するフィルタリング処理を行って前記シンボル系列を生成し、前記入力信号系列として、前記受信信号系列、または、前記仮判定部が前記受信信号系列の一部に対して前記伝送路の推定逆伝達関数による適応等化を行って生成したシンボル系列を含む信号系列が与えられる、請求項１又は２に記載のシンボル判定装置。
【請求項４】
前記パス遡り判定部が出力する前記シンボル系列と、前記受信信号系列から得られるシンボル系列とに基づいて、前記推定伝達関数を算出する更新処理部、を備える請求項１から３のいずれか一項に記載のシンボル判定装置。
【請求項５】
前記適応フィルタ部を備える場合、前記パス遡り判定部が出力する前記シンボル系列と、前記受信信号系列から得られるシンボル系列とに基づいて、前記推定伝達関数、及び前記適応フィルタ部が行う前記フィルタリング処理のフィルタ係数値を算出する更新処理部を備える、請求項３に記載のシンボル判定装置。
【請求項６】
前記所定の推定アルゴリズムは、ビタビ・アルゴリズムであるか、または、ＢＣＪＲアルゴリズムであり、前記ＢＣＪＲアルゴリズムが適用される場合、前記系列推定アルゴリズム処理部は、前記仮判定シンボルと、前記仮判定シンボルの近傍のシンボルからなるシンボル系列に含まれるバイナリ値の対数尤度比が無限大になる場合、当該バイナリ値の前記対数尤度比として予め定められる有限値を用いる、請求項１から５のいずれか一項に記載のシンボル判定装置。
【請求項７】
前記仮判定部は、前記伝送路の逆伝達関数による適応等化を行って生成した前記シンボル系列と、前記仮判定の判定結果として得られた前記シンボル系列とに基づいて、前記推定逆伝達関数を算出する、請求項１から６のいずれか一項に記載のシンボル判定装置。
【請求項８】
伝送路から取り込む受信信号系列に対して前記伝送路の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、前記シンボル系列に対して行う仮判定を行い、
仮判定した仮判定シンボルと前記仮判定シンボルの近傍のシンボルの範囲で前記伝送路状態を示す複数のシンボル系列を生成し、
生成した前記複数のシンボル系列と、前記伝送路の推定伝達関数とに基づいて、前記伝送路状態ごとの推定受信シンボル系列を生成し、
前記受信信号系列から得られるシンボル系列と、前記推定受信シンボル系列の各々とのメトリックを算出し、
算出した前記メトリックと、前記仮判定シンボルと、前記仮判定シンボルの近傍のシンボルとに基づいて所定の推定アルゴリズムにより最も尤もらしい前記推定受信シンボル系列を選択し、
前記最も尤もらしい推定受信シンボル系列に基づいて、トレリスのパスを遡って送信シンボル系列の判定を行うシンボル判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シンボル判定装置、及びシンボル判定方法に関する。
【背景技術】
【０００２】
近年のスマートフォン、タブレットの急速な普及や高精細な動画配信サービスの様なリッチコンテンツの増加等により、インターネットのバックボーンネットワークが転送するトラフィックは増え続けている。また、企業におけるクラウドサービスの活用も進んでいる。これらのことから、データセンタ（以下「ＤＣ」(Data Center)という。）内、ＤＣ間のネットワークのトラフィックが年率約１．３倍の割合で増加することが予測されている。
【０００３】
現在ＤＣ内やＤＣ間の接続方式には主にイーサネット（登録商標）が導入されている。通信トラフィックの増大に伴い、単一拠点におけるＤＣの大規模化が困難になることが予想されている。そのため、今後は今まで以上にＤＣ間連携の必要性が高まり、ＤＣ間で送受信されるトラフィックの更なる増大が考えられる。このような状況に対応するためには、低コストかつ大容量の短距離光伝送技術の確立が求められる。
【０００４】
現行のイーサネット（登録商標）規格では、１０ＧｂＥ(Gigabit Ethernet（登録商標）)−ＺＲを除き４０ｋｍまでの伝送路に光ファイバ通信が適用されている。また、１００ＧｂＥまでは光のｏｎ、ｏｆｆに２値情報を割り当てる強度変調方式が用いられている。受信側は受光器のみで構成され、長距離伝送で用いるコヒーレント受信方式よりも安価な構成となっている。
【０００５】
１００ＧｂＥでは、変調スピードが２５ＧＢｄ(GigaBaud)、シンボルあたりの情報量が１ｂｉｔ／ｓｙｍｂｏｌのＮＲＺ（Non-return-to-zero）信号を４波多重することで、１００Ｇｂ／ｓの伝送容量を実現している。
【０００６】
１００ＧｂＥの次の世代にあたる４００ＧｂＥの標準化においては、１００ＧｂＥで用いられた経済的なデバイス構成の維持と、信号の帯域利用効率を考慮し、初めて２ｂｉｔｓ／ｓｙｍｂｏｌのＰＡＭ４（4-level pulse-amplitude-modulation）が採用されている。これにより、４００ＧｂＥでは、５０Ｇｂ／ｓの信号を８波多重することで４００Ｇｂ／ｓの伝送容量を実現している。４００ＧｂＥの規格として、例えば、４００ＧＢＡＳＥ−ＦＲ８，ＬＲ８などがある。
【０００７】
今後の更なるトラフィック増大に向け、２０２０年以降、８００ＧｂＥ、１．６ＴｂＥの標準化が予定されている。これらの通信速度は、例えば、図１５に示すように、２００Ｇｂ／ｓのＰＡＭ４を採用して変調スピードを１００ＧＢａｕｄとした信号を４〜８波長多重して実現するか、または、２００Ｇｂ／ｓのＰＡＭ８を採用して変調スピードを７５ＧＢａｕｄとした信号を４〜８波長多重して実現することが予定されている。
【０００８】
更なる大容量化に向けた課題として、伝送容量の増大に伴う信号品質劣化、すなわち、デバイスの帯域制限や波長分散の影響が顕在化する。例えば、図１６に示すように、伝送容量が増大して利用帯域が増加すると、デバイスの帯域制限により、周波数領域５０１（斜め線のハッチング領域）が失われてしまう問題が発生する。また、図１７に示すように、伝送容量が増大すると、波長分散の影響が大きくなり、干渉する領域５０２が増大したりする。
【０００９】
このような問題を解決する手法として、高速な通信速度に対応したＤＡＣ(Digital to Analog Converter)やＡＤＣ（Analog to Digital Converter）を用いたり、波長の分散を補償する分散補償モジュール等を用いたりする手法がある。しかし、こういった機器は高価であり、機器に要するコストが高くなるため、経済的な観点からは、回避したい手法である。経済的な観点において、望まれている手法は、従来の送受信器の構成を維持し、多値化、帯域制限耐力、波長分散耐力を向上させて、低コストな狭帯域のデバイスを活用する手法である。しかし、低コストな狭帯域のデバイスを用いる場合、上述したような通信速度の高速化に伴う帯域制限や波長分散による符号間干渉などによって生じる問題を解決する必要がある。
【００１０】
符号間干渉などの問題によって歪んだ受信信号波形から正しい送信データを得るための最も有効な等化方式として、最尤系列推定（以下「ＭＬＳＥ」（Maximum Likelihood Sequential Estimation）という。）が知られている（例えば、非特許文献１、２参照）。
【００１１】
例えば、図１８は、上述した低コストな狭帯域デバイスを用いて構成された従来の通信システム１００を示すブロック図である。通信システム１００は、送信側の信号生成装置１、伝送路２、受信側のシンボル判定装置９０を備える。信号生成装置１は、外部から与えられるｍ値データ信号を取り込み、電気信号の送信信号系列｛ｓｔ｝を生成する。ここで、ｍは、シンボル多値度であり、２以上の整数である。また、ｔは、送信信号系列を識別する識別番号であり、送信信号系列｛ｓｔ｝に含まれるシンボル数がＮ個であるとすると、例えば、１，２，３，…，Ｎといった整数値が割り当てられる。
【００１２】
伝送路２の強度変調器２−２は、信号生成装置１が出力する電気信号の送信信号系列｛ｓｔ｝を取り込み、取り込んだ電気信号の送信信号系列｛ｓｔ｝によって光源２−１が出射する光を変調し、光信号の送信信号系列｛ｓｔ｝を生成する。光ファイバ２−３は、強度変調器２−２が生成した光信号の送信信号系列｛ｓｔ｝を伝送する。受光器２−４は、光ファイバ２−３が伝送する光信号の送信信号系列｛ｓｔ｝を光信号の受信信号系列｛ｒｔ｝として受光し、電気信号の受信信号系列｛ｒｔ｝に変換して出力する。
【００１３】
このとき、伝送路２を等化回路によって示すと、図１９に示すような構成となる。図１９において、遅延器８２−１〜８２−２Ｌは、入力シンボルを取り込んで記憶し、「Ｔ」の時間経過後に記憶した入力シンボルを出力する。ここで、「Ｔ」は、シンボル間隔であり、シンボルごとの演算のタイミングは、「ｔＴ」となる。
【００１４】
遅延器８１は、入力シンボルを取り込んで記憶し、「−ＬＴ」の時間経過後に記憶した入力シンボルを出力する。なお、遅延量にマイナス符号が付いているため、遅延器８１は、「ＬＴ」の負遅延を与えていることになる。ここでは、伝送路２において、時刻ｔにおける符号の前後Ｌシンボル分の符号間干渉が発生するものとして、遅延器８１により光信号の送信信号系列｛ｓｔ｝の要素ｓｔの前後Ｌシンボル分を伝達関数部８３に与えている。
【００１５】
伝達関数部８３は、遅延器８１，８２−１〜８２−２Ｌが出力するシンボル系列に伝達関数（Ｈ）を適用する。加算器８５は、伝達関数部８３の出力値に対して、雑音成分であるωｔを加算して、受信信号系列｛ｒｔ｝を生成する。ωｔは、平均０、分散δ２の互いに独立なガウスランダム系列である。
【００１６】
図１９の等化回路が生成する受信信号系列｛ｒｔ｝を式で示すと次式（１）となる。式（１）において、ｔ＝１，２，…，Ｎである。
【００１７】
【数１】
2020184696.tif 000003
【００１８】
図２０は、伝送路２が出力する受信信号系列｛ｒｔ｝に基づいてＭＬＳＥによって送信シンボル系列の識別を行う識別回路であるシンボル判定装置９０の内部構成を示した図である。伝送路２における伝達関数（Ｈ）は、未知の関数である。そのため、シンボル判定装置９０は、伝送路２の伝達関数（Ｈ）を推定し、推定した伝達関数（Ｈ’）（以下「推定伝達関数（Ｈ’）」という。）を用いて受信シンボル系列のレプリカを生成する。以下、推定伝達関数（Ｈ’）によって複製した受信シンボル系列を推定受信シンボル系列という。シンボル判定装置９０は、生成した推定受信シンボル系列と、受信信号系列｛ｒｔ｝から得られるシンボル系列とを比較し、最も尤もらしい推定受信シンボル系列を判定結果とする。
【００１９】
ＭＬＳＥでは、条件付き結合確率密度関数ｐＮ（｛ｒＮ｝｛ｓ’Ｎ｝）を最大にする送信信号系列｛ｓ’ｔ｝を探索することで、シンボル判定を行う。条件付き結合確率密度関数ｐＮ（｛ｒＮ｝｛ｓ’Ｎ｝）は、伝送路２を通じて、ｍ値データから生成される系列長Ｎの送信信号系列｛ｓ’ｔ｝が送信された場合に、受信信号系列｛ｒｔ｝が受信される確率であり、次式（２）によって表される。
【００２０】
【数２】
2020184696.tif 000004
【００２１】
条件付き結合確率密度関数ｐＮ（｛ｒＮ｝｛ｓ’Ｎ｝）を最大にすることは、次式（３）で示される距離関数ｄＮを最小にすることと等価である。なお、式（３）において（ｐ−１）／２＝Ｌである。
【００２２】
【数３】
2020184696.tif 000005
【００２３】
式（３）における（ｓ’ｔ−（ｐ−１）／２，…，ｓ’ｔ，…，ｓ’ｔ＋（ｐ−１）／２）は、時刻ｔにおける伝送路２の状態（ステート）μｔ（以下「伝送路状態μｔ」という。）を示している。シンボル系列長が「ｐ」の場合、変調シンボルＩ＝［ｉ１，ｉ２，…，ｉｍ］の全ての組み合わせの数は、「ｍｐ」となる。この場合、伝送路２は、ｍｐ個の有限の伝送路状態を有する有限状態機械とみなすことができる。有限状態機械とみなすことができることから、例えば、ビタビ・アルゴリズムのような最尤系列推定の手法を用いて、受信信号系列｛ｒｔ｝ごとに、逐次計算を行って距離関数ｄＮを算出することができる。
【００２４】
時刻ｔにおいて、伝送路状態μｔに到達する距離関数ｄｔ（｛μｔ｝）は、時刻ｔ−１における距離関数ｄｔ−１（｛μｔ−１｝）と、時刻ｔにおける状態遷移に伴う尤度、すなわちメトリックｂ（ｒｔ；μｔ−１→μｔ）とを用いて次式（４）によって表される。
【００２５】
【数４】
2020184696.tif 000006
【００２６】
メトリックｂ（ｒｔ；μｔ−１→μｔ）は、推定伝達関数（Ｈ’）を用いて次式（５）として表される。
【００２７】
【数５】
2020184696.tif 000007
【００２８】
時刻ｔにおけるメトリックｂは、ｔ−１からｔへの状態遷移にのみ依存し、それ以前の状態遷移には依存しない。ここで、伝送路状態μｔに到達する距離関数の最小値ｄ＿ｍｉｎｔ−１（μｔ−１）と、これに対応する全状態遷移が、時刻ｔ−１における全ての伝送路状態μｔ−１において既知であると仮定する。
【００２９】
この仮定の下で、伝送路状態μｔに到達する距離関数ｄｔ（｛μｔ｝）の最小値を求める場合、全ての状態遷移に対応する距離関数ｄｔ（｛μｔ｝）を求める必要はない。伝送路状態μｔに遷移する可能性のある全ての伝送路状態｛μｔ−１｝について、ｄ＿ｍｉｎｔ−１（μｔ−１）＋ｂ（ｒｔ；μｔ−１→μｔ）を算出し、この中の最小値を求めれば、その値が伝送路状態μｔに到達する全ての距離関数ｄｔ（｛μｔ｝）の最小値であるｄ＿ｍｉｎｔ（μｔ）になる。これを式で示すと次式（６）で表される。
【００３０】
【数６】
2020184696.tif 000008
【００３１】
図２０に戻り、シンボル判定装置９０において、遅延器９２−１〜９２−（ｐ−１）は、入力シンボルを取り込んで記憶し、「Ｔ」の時間経過後に記憶した入力シンボルを出力する。そのため、推定伝達関数部９３には、伝送路２の時刻ｔにおける伝送路状態μｔを示す（ｓ’ｔ−（ｐ−１）／２，…，ｓ’ｔ，…，ｓ’ｔ＋（ｐ−１）／２）のシンボル系列が与えられる。推定伝達関数部９３には、このシンボル系列に対して、推定伝達関数（Ｈ’）を適用する。減算器９４は、受信信号系列｛ｒｔ｝から推定伝達関数部９３の出力値を減算する。絶対値器９５は、減算器９４の出力値の絶対値を算出し、算出した絶対値が式（５）で示されるメトリックｂとなる。
【００３２】
加算比較選択部９１は、伝送路状態μｔに遷移する可能性のある全ての伝送路状態｛μｔ−１｝について、式（６）におけるｄ＿ｍｉｎｔ−１（μｔ−１）＋ｂ（ｒｔ；μｔ−１→μｔ）を算出し、算出した値の中の最小値を距離関数ｄｔ（｛μｔ｝）の最小値であるｄ＿ｍｉｎｔ（μｔ）とする。
【００３３】
パス遡り判定部９６は、加算比較選択部９１が算出した距離関数ｄｔ（｛μｔ｝）の最小値に基づいて、ビタビ・アルゴリズムのトレリスのパスを遡り、信号生成装置１が取り込んだｍ値データの推定値を求め、求めた推定値を判定結果として出力する。
【００３４】
上記のように、ビタビ・アルゴリズムでは、伝送路状態μｔに到達する距離関数ｄｔ（｛μｔ｝）の最小値を求める場合、全ての状態遷移に対応する距離関数ｄｔ（｛μｔ｝）を求めずに、伝送路状態μｔに遷移する可能性のある全ての伝送路状態｛μｔ−１｝について、ｄ＿ｍｉｎｔ−１（μｔ−１）＋ｂ（ｒｔ；μｔ−１→μｔ）を算出するようにしている。これにより、ビタビ・アルゴリズムでは、系列長に対し指数的に増大する演算量を線形的な増大に抑えることができる。そのため、ＭＬＳＥでは、ビタビ・アルゴリズムを用いることで、演算量を抑えた最尤系列の推定を可能としている。
【００３５】
また、ＭＬＳＥでは、送信信号波形が受ける符号間干渉を受信器側のデジタル信号処理によって推定して再現することで、高い等化性能を実現する。したがって、推定精度が高いほど符号間干渉による符号誤りを抑制して、符号間干渉によって歪んだ受信信号波形から正しい送信データを得ることが可能になる。ＭＬＳＥを用いた信号品質劣化抑制技術は、上記のイーサネット（登録商標）の大容量化においても検討されている。
【先行技術文献】
【非特許文献】
【００３６】
【非特許文献１】D. D. FALCONER and F. R. MAGEE, JR. “Adaptive Channel Memory Truncation for Maximum Likelihood Sequence Estimation”, The Bell System Technical Journal Vol. 52, No. 9, pp. 1541-1562、November, 1973.
【非特許文献２】Qiang Zhang, Nebojsa Stojanovic, Liang Zhang, Tianjian Zuo, Changsong Xie, Enbo Zhou, “Up to 190-Gb/s OOK Signal Generation using a Coding and Cutting Technique with a 92 GSa/s DAC”, OFC2017, Th3D.1.
【発明の概要】
【発明が解決しようとする課題】
【００３７】
上述したように、ビタビ・アルゴリズムでは、上記の式（６）に示すような演算が行われる。記憶長（拘束長）を「３」として、式（６）の演算対象の組み合わせをトレリス図によって示すと図２１に示すような組み合わせになる。例えば、一番上の枝の場合、［ｉ１，ｉ１，ｉ１］が、式（５）における（ｓ’ｔ−（ｐ−１）／２，…，ｓ’ｔ，…，ｓ’ｔ＋（ｐ−１）／２）に対応する。図２１に示すようにｍ２個の伝送路状態μｔの各々からｍ本の枝が生じるため、枝の合計はｍ３本となる。
【００３８】
例えば、記憶長「３」で、ｍ＝８とした場合、図２２に示すように、伝送路状態μｔの数は、「６４」となり、枝の数は、「５１２」本となる。ビタビ・アルゴリズムでは、伝送路状態μｔごと、すなわちトレリスの枝ごとにメトリックの演算を行う必要がある。そのため、ｍの数、すなわちシンボル多値度が増えて状態数が増加すると、それに伴ってトレリスの枝数が多くなり、ビタビ・アルゴリズムの演算量が多くなるという問題がある。
【００３９】
上記事情に鑑み、本発明は、受信信号系列から送信シンボルの判定を行う際にシンボル多値度が増加した場合であっても、演算量の増大を抑制することができる技術の提供を目的としている。
【課題を解決するための手段】
【００４０】
本発明の一態様は、伝送路から取り込む受信信号系列に対して前記伝送路の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、前記シンボル系列に対して仮判定を行う仮判定部と、伝送路状態を示す複数のシンボル系列と、前記伝送路の推定伝達関数とに基づいて、前記伝送路状態ごとの推定受信シンボル系列を生成する伝送路推定部と、前記受信信号系列から得られるシンボル系列と、前記推定受信シンボル系列の各々とのメトリックを算出し、算出した前記メトリックと、前記仮判定部が仮判定した仮判定シンボルと、前記仮判定シンボルの近傍のシンボルとに基づいて所定の推定アルゴリズムにより最も尤もらしい前記推定受信シンボル系列を選択する系列推定アルゴリズム処理部と、前記最も尤もらしい推定受信シンボル系列に基づいて、トレリスのパスを遡って送信シンボル系列の判定を行うパス遡り判定部と、を備え、前記系列推定アルゴリズム処理部は、前記仮判定シンボルと前記仮判定シンボルの近傍のシンボルの範囲で前記伝送路状態を示す複数のシンボル系列を生成して前記伝送路推定部に出力するシンボル判定装置である。
【００４１】
本発明の一態様は、上記のシンボル判定装置であって、前記近傍のシンボルとは、前記仮判定シンボルを中心とした前後のシンボルであり、前記仮判定シンボルと前記近傍のシンボルの総和が（２ｎ＋１）個（ｎは、２以上の整数）になるシンボルである。
【００４２】
本発明の一態様は、上記のシンボル判定装置であって、前記受信信号系列から得られるシンボル系列は、前記受信信号系列に含まれるシンボル系列であるか、または、前記仮判定部が前記受信信号系列の一部に対して前記推定逆伝達関数による適応等化を行って生成したシンボル系列であるか、または、更に備える適応フィルタ部によって生成されるシンボル系列であり、前記適応フィルタ部は、入力信号系列の一部に対してパルス幅を圧縮するフィルタリング処理を行って前記シンボル系列を生成し、前記入力信号系列として、前記受信信号系列、または、前記仮判定部が前記受信信号系列の一部に対して前記伝送路の推定逆伝達関数による適応等化を行って生成したシンボル系列を含む信号系列が与えられる。
【００４３】
本発明の一態様は、上記のシンボル判定装置であって、前記パス遡り判定部が出力する前記シンボル系列と、前記受信信号系列から得られるシンボル系列とに基づいて、前記推定伝達関数を算出する更新処理部を備える。
【００４４】
本発明の一態様は、上記のシンボル判定装置であって、前記適応フィルタ部を備える場合、前記パス遡り判定部が出力する前記シンボル系列と、前記受信信号系列から得られるシンボル系列とに基づいて、前記推定伝達関数、及び前記適応フィルタ部が行う前記フィルタリング処理のフィルタ係数値を算出する更新処理部を備える。
【００４５】
本発明の一態様は、上記のシンボル判定装置であって、前記所定の推定アルゴリズムは、ビタビ・アルゴリズムであるか、または、ＢＣＪＲアルゴリズムであり、前記ＢＣＪＲアルゴリズムが適用される場合、前記系列推定アルゴリズム処理部は、前記仮判定シンボルと、前記仮判定シンボルの近傍のシンボルからなるシンボル系列に含まれるバイナリ値の対数尤度比が無限大になる場合、当該バイナリ値の前記対数尤度比として予め定められる有限値を用いる。
【００４６】
本発明の一態様は、上記のシンボル判定装置であって、前記仮判定部は、前記伝送路の逆伝達関数による適応等化を行って生成した前記シンボル系列と、前記仮判定の判定結果として得られた前記シンボル系列とに基づいて、前記推定逆伝達関数を算出する。
【００４７】
本発明の一態様は、伝送路から取り込む受信信号系列に対して前記伝送路の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、前記シンボル系列に対して仮判定を行い、仮判定した仮判定シンボルと前記仮判定シンボルの近傍のシンボルの範囲で伝送路状態を示す複数のシンボル系列を生成し、生成した前記複数のシンボル系列と、前記伝送路の推定伝達関数とに基づいて、前記伝送路状態ごとの推定受信シンボル系列を生成し、前記受信信号系列から得られるシンボル系列と、前記推定受信シンボル系列の各々とのメトリックを算出し、算出した前記メトリックと、前記仮判定シンボルと、前記仮判定シンボルの近傍のシンボルとに基づいて所定の推定アルゴリズムにより最も尤もらしい前記推定受信シンボル系列を選択し、前記最も尤もらしい推定受信シンボル系列に基づいて、トレリスのパスを遡って送信シンボル系列の判定を行うシンボル判定方法である。
【発明の効果】
【００４８】
本発明により、受信信号系列から送信シンボルの判定を行う際にシンボル多値度が増加した場合であっても、演算量の増大を抑制することができる。
【図面の簡単な説明】
【００４９】
【図１】第１の実施形態における通信システムの構成を示すブロック図である。
【図２】第１の実施形態におけるシンボル判定装置の内部構成を示すブロック図である。
【図３】第１の実施形態におけるシンボル判定装置の内部の詳細構成を示すブロック図である。
【図４】第１の実施形態における適応フィルタ部が取り込む系列を説明する図である。
【図５】第１の実施形態における系列推定部の適応フィルタ部が行うパルス幅の圧縮を説明する図である。
【図６】第１の実施形態におけるシンボル判定装置の仮判定部による処理の流れを示すフローチャートである。
【図７】第１の実施形態におけるシンボル判定装置の系列推定部による処理の流れを示すフローチャートである。
【図８】第１の実施形態における系列推定部の近傍加算比較選択部が演算対象とする状態を示すトレリス図である。
【図９】図８においてシンボル多値度を８とした場合のトレリス図である。
【図１０】第２の実施形態におけるシンボル判定装置の内部構成を示すブロック図である。
【図１１】第２の実施形態におけるシンボル判定装置の内部の詳細構成を示すブロック図である。
【図１２】シンボル多値度が８の場合のｇｒａｙ符号の一例を示す図である。
【図１３】ＢＣＪＲアルゴリズムを適用する場合の課題を説明する図である。
【図１４】ＢＣＪＲアルゴリズムを適用する場合にシンボル判定装置が備えるテーブルの一例を示す図である。
【図１５】伝送容量と波長の多重度との関係を示すグラフである。
【図１６】伝送容量が増大した場合のデバイスの帯域制限の影響を示すグラフである。
【図１７】伝送容量が増大した場合の波長分散の影響を示すグラフである。
【図１８】従来の通信システムの構成を示すブロック図である。
【図１９】伝送路の等化回路のブロック図である。
【図２０】従来のシンボル判定装置の内部構成を示すブロック図である。
【図２１】ビタビ・アルゴリズムの演算対象となる状態を示すトレリス図である。
【図２２】図２１においてシンボル多値度を８とした場合のトレリス図である。
【発明を実施するための形態】
【００５０】
（第１の実施形態）
以下、本発明の実施形態について図面を参照して説明する。
図１は、第１の実施形態における通信システムＳの構成を示すブロック図である。通信システムＳは、信号生成装置１、伝送路２、及びシンボル判定装置３を備える。
【００５１】
信号生成装置１及び伝送路２は、図１８に示した従来の通信システム１００が備える信号生成装置１及び伝送路２と同一構成である。信号生成装置１は、外部から与えられるｍ値データ信号を取り込む。ｍは、シンボル多値度であり、例えば２以上の整数である。各シンボルは、数字や記号で表される。例えば、ｍ＝８の場合、各シンボルは、［０，１，２，３，４，５，６，７］の数字で表される。
信号生成装置１は、取り込んだｍ値データ信号を送信シンボル系列とし、送信シンボル系列を含む電気信号の送信信号系列｛ｓｔ｝を生成する。ｔは、送信信号系列を識別する識別番号である。送信信号系列｛ｓｔ｝に含まれるシンボル数がＮ個であるとすると、例えばｔ＝１，２，３，…，Ｎといった整数値が割り当てられる。
【００５２】
伝送路２において、強度変調器２−２は、信号生成装置１が出力する電気信号の送信信号系列｛ｓｔ｝を取り込む。強度変調器２−２は、電気信号の送信信号系列｛ｓｔ｝に含まれているｍ値のシンボル系列に基づいて光源２−１が出射する光を変調する。これにより、強度変調器２−２は、ｍ値のシンボル系列を表す光信号の送信信号系列｛ｓｔ｝を生成する。
【００５３】
光ファイバ２−３は、強度変調器２−２が生成した光信号の送信信号系列｛ｓｔ｝を伝送する。
受光器２−４は、光ファイバ２−３が伝送する光信号の送信信号系列｛ｓｔ｝を光信号の受信信号系列｛ｒｔ｝として受光する。受光器２−４は、受光した光信号の受信信号系列｛ｒｔ｝を電気信号の受信信号系列｛ｒｔ｝に変換して出力する。受光器２−４は、例えばフォトダイオードである。
【００５４】
シンボル判定装置３は、伝送路２が出力する受信信号系列｛ｒｔ｝に基づいて、送信シンボル系列の識別を行う識別回路である。シンボル判定装置３は、図２に示す内部構成を備えている。シンボル判定装置３は、仮判定部３０、系列推定部４０を備える。仮判定部３０には、例えば、ＦＦＥ（Feed Forward Equalizer）が適用される。仮判定部３０は、逆伝達関数を推定した関数（以下「推定逆伝達関数」という。）により受信信号系列｛ｒｔ｝を適応等化して硬判定を行い、送信シンボル系列の仮判定を行う。
【００５５】
系列推定部４０は、伝送路状態を示すシンボル系列に対して推定伝達関数（Ｈ’）を適用して推定受信シンボル系列を生成する。系列推定部４０は、生成した推定受信シンボル系列と、パルス幅を圧縮した受信信号系列｛ｒｔ｝とに基づいてメトリックを算出する。また、系列推定部４０は、算出したメトリックを用いて、仮判定部３０が仮判定したシンボルの各々を中心とする近傍のシンボルの範囲でビタビ・アルゴリズムを実行する。これにより、系列推定部４０は、送信シンボル系列の推定値、すなわち信号生成装置１が取り込んだｍ値データ信号の推定値を求める。
【００５６】
仮判定部３０は、適応フィルタ部３０１、判定処理部３０２、及び更新処理部３０３を備える。適応フィルタ部３０１は、例えば、図３に示すように線形トランスバーサルフィルタである。適応フィルタ部３０１は、伝送路２の伝達関数（Ｈ）の逆伝達関数を推定した推定逆伝達関数により入力信号を適応等化する。
【００５７】
適応フィルタ部３０１は、図３に示すように、遅延器３１，３２−１〜３２−（ｕ−１）、タップ３３−１〜３３−ｕ、及び加算器３４を備える。遅延器３１は、図４に示すように、Ｎ個からなる受信信号系列｛ｒｔ｝の一部であるｕ個のシンボル系列を取り込む。遅延器３１は、取り込んだシンボル系列を記憶し、「（ｕ−１）Ｔ／２」の時間経過後、すなわち「（ｕ−１）／２」シンボル分、遅延させて取り込んだシンボル系列を出力する。
【００５８】
タップ３３−１には、遅延器３１が出力する「（ｕ−１）／２」シンボル分、遅延したシンボルであるｒｔ＋（ｕ−１）／２が与えられる。
【００５９】
遅延器３２−１〜３２−（ｕ−１）の各々は、１つのシンボルを取り込んで記憶し、「Ｔ」の時間経過後、すなわち１シンボル分、遅延させて取り込んだシンボルを出力する。例えば、最初の遅延器３２−１は、「（ｕ−３）／２」シンボル分、遅延したｒｔ＋（ｕ−３）／２のシンボルを出力する。最後の遅延器３２−（ｕ−１）は、「（ｕ−１）／２」シンボル分、遅延したｒｔ＋（ｕ−１）／２のシンボルを出力する。これにより、タップ３３−１〜３３−ｕには、次式（７）で示される系列長ｕのシンボル系列を含んだ信号が与えられることになる。
【００６０】
【数７】
2020184696.tif 000009
【００６１】
タップ３３−１〜３３−ｕの各々には、いわゆるフィルタ係数値であるｆ１，ｆ２，…，ｆ（ｕ＋１）／２，…，ｆｕのタップ利得値が設定されている。このタップ利得値ｆ１〜ｆｕが、推定逆伝達関数を表していることになる。
【００６２】
タップ３３−１〜３３−ｕは、各々に与えられるシンボルに対して各々のタップ利得値ｆ１〜ｆｕを乗算して出力する。加算器３４は、タップ３３−１〜３３−ｕの出力値を合計して出力する。式（７）は、「（ｕ＋１）／２」番目の要素であるｒｔを中心とした系列ということができるため、加算器３４の出力値を次式（８）のように表すことができる。
【００６３】
【数８】
2020184696.tif 000010
【００６４】
適応フィルタ部３０１では、遅延器３１により「（ｕ−１）Ｔ／２」の遅延を与えている。そのため、適応フィルタ部３０１による線形デジタルフィルタリングの演算タイミングｔＴの出力値に対応する入力情報であるｕ個のシンボル系列は、「（ｕ−１）Ｔ／２」分の遅延が発生していることになる。
【００６５】
判定処理部３０２は、適応フィルタ部３０１の出力値に対して硬判定による仮判定を行い、出力値に対応する送信シンボルの推定値を求める。判定処理部３０２は、求めた推定値である仮判定シンボルＡ’を仮判定結果として出力する。
【００６６】
更新処理部３０３は、適応フィルタ部３０１の出力値の目標値を判定処理部３０２が出力する仮判定シンボルＡ’として、適応フィルタ部３０１のタップ３３−１〜３３−ｕの各々のタップ利得値ｆ１〜ｆｕの更新値を算出する。例えば、更新処理部３０３は、タップ利得値ｆ１〜ｆｕの更新値、すなわち推定逆伝達関数を所定の更新アルゴリズムにより算出する。
【００６７】
更新処理部３０３は、図３に示すように、フィルタ更新アルゴリズム処理部３５と、減算器３６とを備える。更新処理部３０３において、減算器３６は、判定処理部３０２が出力する仮判定シンボルＡ’から適応フィルタ部３０１の出力値を減算して得られる減算値を誤差としてフィルタ更新アルゴリズム処理部３５に出力する。
【００６８】
フィルタ更新アルゴリズム処理部３５は、減算器３６が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムによりタップ利得値ｆ１〜ｆｕの更新値を算出する。フィルタ更新アルゴリズム処理部３５は、算出したタップ利得値ｆ１〜ｆｕを、タップ３３−１〜３３−ｕに設定して、タップ利得値ｆ１〜ｆｕの更新を行う。
【００６９】
系列推定部４０は、適応フィルタ部４０１、系列推定アルゴリズム処理部４０２、伝送路推定部４０３、更新処理部４０４、及びパス遡り判定部４０５を備える。適応フィルタ部４０１は、例えば、図３に示すように線形トランスバーサルフィルタであり、伝送路推定部４０３の記憶長を削減するために、受信信号系列｛ｒｔ｝のインパルスレスポンスを圧縮する。ここで、インパルスレスポンスの圧縮とは、図５に示すように、帯域制限や波長分散のために時間的に広がった受信信号系列｛ｒｔ｝のパルス幅を圧縮することであり、圧縮によりシンボル間の干渉を削減することができる。
【００７０】
適応フィルタ部４０１は、図３に示すように、遅延器４１，４２−１〜４２−（ｖ−１）、タップ４３−１〜４３−ｖ、及び加算器４４を備える。遅延器４１は、遅延器３１と同様に、図４に示すように、Ｎ個からなる受信信号系列｛ｒｔ｝の一部であるｖ個のシンボル系列を取り込む。ここで、ｖは、ｕと同一値であってもよいし、異なる値であってもよい。
【００７１】
遅延器４１は、取り込んだシンボル系列を記憶し、「（ｖ−１）Ｔ／２」の時間経過後、すなわち「（ｖ−１）／２」シンボル分、遅延させて取り込んだシンボル系列を出力する。タップ４３−１には、遅延器４１が出力するシンボルであるｒｔ＋（ｖ−１）／２が与えられる。
【００７２】
遅延器４２−１〜４２−（ｖ−１）の各々は、１つのシンボルを取り込んで記憶し、「Ｔ」の時間経過後、すなわち１シンボル分、遅延させて取り込んだシンボルを出力する。例えば、最初の遅延器４２−１は、「（ｖ−３）／２」シンボル分、遅延したｒｔ＋（ｖ−３）／２のシンボルを出力する。最後の遅延器４２−（ｖ−１）は、「（ｖ−１）／２」シンボル分、遅延したｒｔ＋（ｖ−１）／２のシンボルを出力する。これにより、タップ４３−１〜４３−ｖには、次式（９）で示される系列長ｖのシンボル系列を含んだ信号が与えられることになる。
【００７３】
【数９】
2020184696.tif 000011
【００７４】
タップ４３−１〜４３−ｖの各々には、いわゆるフィルタ係数値であるｃ１，ｃ２，…，ｃ（ｖ＋１）／２，…，ｃｖのタップ利得値が設定されている。タップ４３−１〜４３−ｖは、各々に与えられるシンボルに対して各々のタップ利得値を乗算して出力する。加算器４４は、タップ４３−１〜４３−ｖの出力値を合計して出力する。式（９）は、「（ｖ＋１）／２」番目の要素であるｒｔを中心とした系列ということができるため、加算器４４の出力値を式によって表すと、次式（１０）になる。
【００７５】
【数１０】
2020184696.tif 000012
【００７６】
式（１０）から分かるように、適応フィルタ部４０１は、タップ利得値ｃ１，ｃ２，…，ｃ（ｖ＋１）／２，…，ｃｖにより影響度合いが調整されているが、ｖ個分のシンボル系列の情報量を圧縮した１つの出力シンボルを出力していることになる。
【００７７】
ＭＬＳＥの演算量は、パルスの広がり幅に対して指数的に増大することが知られているが、適応フィルタ部４０１によってパルス幅を圧縮することにより演算量の増大を抑えることができる。
【００７８】
適応フィルタ部４０１では、遅延器４１により「（ｖ−１）Ｔ／２」の遅延を与えている。そのため、適応フィルタ部４０１による線形デジタルフィルタリングの演算タイミングｔＴの出力値に対応する入力情報であるｖ個のシンボル系列は、「（ｖ−１）Ｔ／２」分の遅延が発生していることになる。
【００７９】
伝送路推定部４０３には、時刻ｔにおける伝送路２の伝送路状態μｔを表すシンボル系列のうち、仮判定部３０が仮判定したシンボルの各々を中心とする近傍のシンボルの範囲で表されるシンボル系列が与えられる。伝送路推定部４０３は、近傍加算比較選択部５１から与えられるシンボル系列の各々に対して推定伝達関数（Ｈ’）を適用してシンボル系列ごとの推定受信シンボル系列を生成する。
【００８０】
伝送路推定部４０３は、例えば、図３に示すように線形トランスバーサルフィルタである。伝送路推定部４０３は、遅延器６２−１〜６２−（ｘ−１）、タップ６１−１〜６１−ｘ、及び加算器６３を備える。遅延器６２−１〜６２−（ｘ−１）の各々は、１つのシンボルを取り込んで記憶し、「Ｔ」の時間経過後、すなわち１シンボル分、遅延させて取り込んだシンボルを出力する。近傍加算比較選択部５１から与えられる時刻ｔにおける伝送路２の伝送路状態μｔは、次式（１１）に示すシンボル系列｛ｓ’ｔ｝として表すことができる。なお、ｘは、記憶長の値に一致し、記憶長が「３」である場合、ｘ＝３となる。
【００８１】
【数１１】
2020184696.tif 000013
【００８２】
タップ６１−１〜６１−ｘの各々には、式（１１）のシンボル系列｛ｓ’ｔ｝に含まれるシンボルの各々が与えられる。タップ６１−１〜６１−ｘの各々には、いわゆるフィルタ係数値である推定伝達関数（Ｈ’）の係数値ｈ１，ｈ２，…，ｈ（ｘ＋１）／２，…，ｈｘが設定されている。例えば、タップ６１−１は、ｈ１×ｓ’ｔ−（ｘ−１）／２の演算を行う。タップ６１−２〜６１−ｘも、同様に、各々の係数値と、与えられるシンボルとを乗算し、乗算結果を加算器６３に出力する。加算器６３は、乗算結果の合計値を出力する。加算器６３の出力値が、推定受信シンボル系列を構成するシンボルとなり、次式（１２）として表される。
【００８３】
【数１２】
2020184696.tif 000014
【００８４】
系列推定アルゴリズム処理部４０２は、伝送路状態μｔごとのメトリックを算出する。系列推定アルゴリズム処理部４０２は、算出した伝送路状態μｔごとのメトリックを用いて、仮判定部３０が仮判定したシンボルの各々を中心とする近傍のシンボルの範囲でビタビ・アルゴリズムを実行する。
【００８５】
系列推定アルゴリズム処理部４０２は、減算器５４、絶対値器５３、及び近傍加算比較選択部５１を備える。減算器５４は、式（１０）に示す適応フィルタ部４０１の出力値から、式（１２）に示す伝送路推定部４０３の出力値を減算する。減算器５４は、減算により得られる減算値を絶対値器５３に出力する。絶対値器５３は、減算器５４から受けた減算値の絶対値を算出する。絶対値器５３が算出した絶対値が、メトリックであり次式（１３）として表される。
【００８６】
【数１３】
2020184696.tif 000015
【００８７】
近傍加算比較選択部５１は、仮判定部３０が仮判定結果として出力する複数の仮判定シンボルＡ’を取りこむ。近傍加算比較選択部５１は、複数の仮判定シンボルＡ’を用いて、予め定められる記憶長の長さを系列長とする仮判定シンボル系列｛Ａ’ｔ｝を生成する。近傍加算比較選択部５１は、時刻ｔにおける伝送路２の伝送路状態μｔを表すシンボル系列のうち、仮判定シンボル系列｛Ａ’ｔ｝に含まれるシンボルの各々を中心とする近傍のシンボルの範囲で表される複数のシンボル系列｛ｓ’ｔ｝を生成する。近傍加算比較選択部５１は、生成した複数のシンボル系列｛ｓ’ｔ｝を伝送路推定部４０３に出力する。
【００８８】
近傍加算比較選択部５１は、絶対値器５３が出力する伝送路状態μｔごとのメトリックを用いて、仮判定シンボル系列｛Ａ’ｔ｝に含まれるシンボルの各々を中心とする近傍のシンボルの範囲でビタビ・アルゴリズムを実行する。近傍加算比較選択部５１は、ビタビ・アルゴリズムを実行することにより、推定受信シンボル系列の尤もらしさを示す距離関数ｄｔ（｛μｔ｝）を算出し、算出した距離関数ｄｔ（｛μｔ｝）の最小値を検出する。この最小値に対応する推定受信シンボル系列が、最も尤もらしい推定受信シンボル系列になる。
【００８９】
パス遡り判定部４０５は、近傍加算比較選択部５１が検出した距離関数ｄｔ（｛μｔ｝）の最小値に基づいて、トレリスのパスを遡って、送信シンボルの判定を行う。なお、遡るパスの始点は、時刻ｔにおいて伝送路状態μｔに到達する際の距離関数ｄｔ（｛μｔ｝）が最小値になる伝送路状態である。また、パス遡り判定部４０５がパスを遡る際の遡り数「ｗ」は予め定められており、遡り数「ｗ」を固定値にすることでパスを判定する演算量を削減することができる。なお、伝送路推定部４０３の記憶長の数倍程度遡ることで、パスが収束することが知られている。
【００９０】
パス遡り判定部４０５が、シンボル判定することにより得られる各々のシンボルを判定シンボルＡとする。パス遡り判定部４０５は、判定シンボルＡを判定結果として出力する。パス遡り判定部４０５が順次判定する判定シンボルＡを系列に並べた判定シンボル系列｛Ａｔ｝が、送信シンボル系列の推定値、すなわち信号生成装置１が取り込んだｍ値データの推定値となる。
【００９１】
更新処理部４０４は、適応フィルタ部４０１のタップ利得値ｃ１〜ｃｖと、フィルタ更新用伝送路推定部７０の推定伝達関数の係数値ｈ１〜ｈｘを算出する。
【００９２】
更新処理部４０４は、適応フィルタ部４０１のタップ利得値ｃ１〜ｃｖを算出する際、適応フィルタ部４０１の出力値の目標値を、判定結果である判定シンボル系列｛Ａｔ｝を入力情報とするフィルタ更新用伝送路推定部７０の出力値とする。更新処理部４０４は、目標値になるように、適応フィルタ部４０１のタップ４３−１〜４３−ｖの各々のタップ利得値ｃ１〜ｃｖの更新値を所定の更新アルゴリズムにより算出する。
【００９３】
また、更新処理部４０４は、フィルタ更新用伝送路推定部７０の推定伝達関数の係数値ｈ１〜ｈｘを算出する際、フィルタ更新用伝送路推定部７０の出力値の目標値を、適応フィルタ部４０１の出力値とする。更新処理部４０４は、目標値になるように、フィルタ更新用伝送路推定部７０のタップ７１−１〜７１−ｘの各々の係数値ｈ１〜ｈｘの更新値を所定の更新アルゴリズムにより算出する。更新処理部４０４が算出した係数値ｈ１〜ｈｘの更新値は、伝送路推定部４０３のタップ６１−１〜６１−ｘにも適用される。
【００９４】
更新処理部４０４は、図３に示すように、フィルタ更新用伝送路推定部７０、フィルタ更新アルゴリズム処理部７５、遅延器７６、及び減算器７７を備える。フィルタ更新用伝送路推定部７０の構成は、伝送路推定部４０３の構成に対応しており、タップ６１−１〜６１−ｘが、タップ７１−〜７１−ｘに対応し、遅延器６２−１〜６２−（ｘ−１）が、遅延器７２−１〜７２−（ｘ−１）に対応し、加算器６３が、加算器７３に対応する。
【００９５】
遅延器７６は、適応フィルタ部４０１の出力値を「−ｗＴ」の時間、遅延させて減算器７７に出力する。「−ｗＴ」の時間、遅延させるのは、パス遡り判定部４０５における処理により、「−ｗＴ」の遅延が生じるためである。遅延器７６による「−ｗＴ」の遅延により、フィルタ更新用伝送路推定部７０の出力値と、適応フィルタ部４０１の出力値のタイミングが一致することになる。
【００９６】
減算器７７は、フィルタ更新用伝送路推定部７０の出力値から、「−ｗＴ」時間遅延した適応フィルタ部４０１の出力値を減算し、減算により得られた誤差をフィルタ更新アルゴリズム処理部７５に出力する。
【００９７】
フィルタ更新アルゴリズム処理部７５は、減算器７７が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムによりタップ利得値ｃ１〜ｃｖの更新値を算出する。また、フィルタ更新アルゴリズム処理部７５は、減算器３６が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムにより係数値ｈ１〜ｈｘの更新値を算出する。
【００９８】
フィルタ更新アルゴリズム処理部７５は、算出したタップ利得値ｃ１〜ｃｖを、タップ４３−１〜４３−ｖに設定して、タップ利得値ｃ１〜ｃｖの更新を行う。フィルタ更新アルゴリズム処理部７５は、算出した係数値ｈ１〜ｈｘを、タップ７１−１〜７１−ｘ、及び伝送路推定部４０３のタップ６１−１〜６１−ｘに設定して、係数値ｈ１〜ｈｘの更新を行う。
【００９９】
（第１の実施形態におけるシンボル判定装置による処理）
次に、図６から図９を参照しつつ、第１の実施形態におけるシンボル判定装置３による処理について説明する。
【０１００】
（第１の実施形態における仮判定部による処理）
図６は、シンボル判定装置３の仮判定部３０による処理の流れを示すフローチャートである。
【０１０１】
適応フィルタ部３０１の遅延器３１が、受信信号系列｛ｒｔ｝から系列長ｕのシンボル系列を取り込んで記憶する（ステップＳａ１）。遅延器３１は、取り込んだシンボル系列を、「（ｕ−１）Ｔ／２」時間遅延させて出力する。遅延器３２−１〜３２−（ｕ−１）の各々は、遅延器３１が順次出力するシンボルを取り込んで記憶し、記憶したシンボルを「Ｔ」時間経過後に出力する。
【０１０２】
これにより、式（７）で示される受信信号系列｛ｒｔ｝のシンボル系列が、タップ３３−１〜３３−ｕに与えられる。タップ３３−１〜３３−ｕは、各々に与えられるシンボルｒｔ−（ｕ−１）／２〜ｒｔ＋（ｕ−１）／２と、各々に設定されているタップ利得値ｆ１〜ｆｕとを乗算し、乗算した結果を加算器３４に出力する。加算器３４は、乗算結果を合計して式（８）で示される出力値を算出して出力する。この出力値が、推定逆伝達関数により適応等化された受信信号系列｛ｒｔ｝のシンボルとなる（ステップＳａ２）。
【０１０３】
判定処理部３０２は、適応フィルタ部３０１の出力値に対して硬判定による仮判定を行い、送信シンボルの推定値を求める。判定処理部３０２は、求めた推定値である仮判定シンボルＡ’を仮判定結果として出力する（ステップＳａ３）。
【０１０４】
減算器３６は、判定処理部３０２が出力する仮判定シンボルＡ’から適応フィルタ部３０１の出力値を減算して得られる減算値を誤差としてフィルタ更新アルゴリズム処理部３５に出力する。
【０１０５】
フィルタ更新アルゴリズム処理部３５は、減算器３６が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムによりタップ利得値ｆ１〜ｆｕの更新値、すなわち推定逆伝達関数を算出する。フィルタ更新アルゴリズム処理部３５は、算出したタップ利得値ｆ１〜ｆｕを、タップ３３−１〜３３−ｕに設定して、タップ利得値ｆ１〜ｆｕの更新を行う（ステップＳａ４）。
【０１０６】
適応フィルタ部３０１の遅延器３１が、前回のステップＳａ１において取り込んだ系列長ｕのシンボル系列の先頭から１シンボルずらしたシンボルを先頭として、受信信号系列｛ｒｔ｝から系列長ｕのシンボル系列を取り込める場合（ステップＳａ５、Ｙｅｓ）、ステップＳａ１の処理を行う。一方、遅延器３１が、前回のステップＳａ１において取り込んだ系列長ｕのシンボル系列の先頭から１シンボルずらしたシンボルを先頭として、受信信号系列｛ｒｔ｝から系列長ｕのシンボル系列を取り込めない場合（ステップＳａ５、Ｎｏ）、処理を終了する。
【０１０７】
（第１の実施形態における系列推定部による処理）
図７は、シンボル判定装置３の系列推定部４０による処理の流れを示すフローチャートである。適応フィルタ部４０１の遅延器４１が、受信信号系列｛ｒｔ｝から系列長ｖのシンボル系列を取り込んで記憶する（ステップＳｂ１−１）。遅延器４１は、取り込んだシンボル系列を「（ｖ−１）Ｔ／２」時間遅延させて出力する。遅延器４２−１〜４２−（ｖ−１）の各々は、遅延器４１が順次出力するシンボルを取り込んで記憶し、記憶したシンボルを「Ｔ」時間経過後に出力する。
【０１０８】
これにより、式（９）で示される受信信号系列｛ｒｔ｝のシンボル系列が、タップ４３−１〜４３−ｖに与えられる。タップ４３−１〜４３−ｖは、各々に与えられるシンボルｒｔ−（ｖ−１）／２〜ｒｔ＋（ｖ−１）／２と、各々に設定されているタップ利得値ｃ１〜ｃｖとを乗算し、乗算した結果を加算器４４に出力する。加算器４４は、乗算結果を合計して式（１０）により示される出力値を算出し、算出した出力値を受信信号系列｛ｒｔ｝の情報量を圧縮したシンボルとして出力する（ステップＳｂ１−２）。
【０１０９】
ステップＳｂ１−１，Ｓｂ１−２の処理と並列に近傍加算比較選択部５１は、仮判定部３０が仮判定結果として出力する複数の仮判定シンボルＡ’を取り込む（ステップＳｂ２−１）。近傍加算比較選択部５１は、予め定められる記憶長の長さを系列長とする仮判定シンボル系列｛Ａ’ｔ｝を生成する。例えば、近傍加算比較選択部５１が取り込んだシンボルを並べた、シンボル系列が［Ａ’ｔ−（ｐ−１）／２，…，Ａ’ｔ−１，Ａ’ｔ，Ａ’ｔ＋１，…，Ａ’ｔ＋（ｐ−１）／２］であるとする。
【０１１０】
このとき、記憶長が「３」、近傍のシンボル範囲として隣接±１シンボルとすることが予め定められているとする。近傍加算比較選択部５１は、記憶長が「３」であることから、例えば、時刻ｔにおいて「Ａ’ｔ」と、「Ａ’ｔ」の前後１シンボルを含んだ［Ａ’ｔ−１，Ａ’ｔ，Ａ’ｔ＋１］を仮判定シンボル系列｛Ａ’ｔ｝として選択する。
【０１１１】
近傍加算比較選択部５１は、図８に示すように、トレリス図において、中央付近に実線で示す仮判定シンボル系列［Ａ’ｔ−１，Ａ’ｔ，Ａ’ｔ＋１］に含まれるシンボルの各々の隣接±１のシンボルの範囲を時刻ｔにおける伝送路２の伝送路状態μｔを示すシンボル系列｛ｓ’ｔ｝とする。
【０１１２】
すなわち、（１）［Ａ’ｔ−１］と、［Ａ’ｔ−１］の隣接±１である［Ａ’ｔ−１＋１］と［Ａ’ｔ−１−１］の３シンボルと、（２）［Ａ’ｔ］と、［Ａ’ｔ］の隣接±１である［Ａ’ｔ＋１］と［Ａ’ｔ−１］の３シンボルと、（３）［Ａ’ｔ＋１］と、［Ａ’ｔ＋１］の隣接±１である［Ａ’ｔ＋１＋１］と［Ａ’ｔ＋１−１］の３シンボルによって示されるトレリスの枝、すなわち、図８の破線と１本の実線で示されるトレリスの枝を伝送路状態μｔの範囲とする。
【０１１３】
上述したように、シンボル多値度が「ｍ」であり、シンボル系列長が「ｐ」である場合、トレリスの枝は、ｍｐ本になる。仮に、記憶長を「３」にしたとしても枝数は、ｍ３本となる。これに対して、上記のように、隣接±１シンボルにすることで、ｍ＝３にすることができるので、図８に示すように、枝数を３３＝２７本に絞り込むことができる。
【０１１４】
近傍加算比較選択部５１は、仮判定シンボルの各々を中心とした近傍のシンボルの範囲に絞った伝送路状態μｔの各々を表すシンボル系列を伝送路推定部４０３に出力する。このシンボル系列は、式（１１）で表され、記憶長が「３」の場合、［ｓ’ｔ−１，ｓ’ｔ，ｓ’ｔ＋１］になる。
【０１１５】
伝送路推定部４０３の遅延器６２−１〜６２−（ｘ−１）の各々は、取り込んだシンボルを記憶し、「Ｔ」時間経過後に取り込んだシンボルを出力する。タップ６１−１〜６１−ｘの各々には、シンボル系列｛ｓ’ｔ｝に含まれるシンボルの各々が与えられる。
【０１１６】
タップ６１−１〜６１−ｘは、各々の係数値ｈ１〜ｈｘと、与えられたシンボルとを乗算し、乗算結果を加算器６３に出力する。加算器６３は、乗算結果を合計して式（１２）により示される出力値を算出して出力する。加算器６３が順次出力する出力値が、推定受信シンボル系列を構成するシンボルになる（ステップＳｂ２−２）。
【０１１７】
減算器５４は、式（１０）により示される適応フィルタ部４０１の出力値から、式（１２）により示される伝送路推定部４０３の出力値を減算し、減算により得られる減算値を絶対値器５３に出力する。絶対値器５３は、減算器５４から受けた減算値の絶対値を算出する。絶対値器５３が算出した絶対値が、式（１３）として表されるメトリックとなる（ステップＳｂ３）。
【０１１８】
例えば、近傍加算比較選択部５１が、トレリスの枝数を２７本に絞った場合、伝送路状態μｔの各々とは、図８に示す最初の枝［Ａ’ｔ−１−１，Ａ’ｔ−１，Ａ’ｔ＋１−１］から最後の枝の［Ａ’ｔ−１＋１，Ａ’ｔ＋１，Ａ’ｔ＋１＋１］までの２７個のシンボル系列｛ｓ’ｔ｝になる。伝送路推定部４０３は、２７個のシンボル系列｛ｓ’ｔ｝の各々に対して推定伝達関数（Ｈ’）を適用し、式（１２）によって示される２７個の出力値を算出する。したがって、絶対値器５３も、２７個のメトリックを算出することになる。
【０１１９】
近傍加算比較選択部５１は、図８に示すトレリス図において、中央付近に実線で示す仮判定シンボル系列［Ａ’ｔ−１，Ａ’ｔ，Ａ’ｔ＋１］に含まれるシンボルの各々の隣接±１のシンボルの範囲でビタビ・アルゴリズムを実行する。近傍加算比較選択部５１は、ビタビ・アルゴリズムを実行して、距離関数ｄｔ（｛μｔ｝）を算出し、算出した距離関数ｄｔ（｛μｔ｝）の最小値を検出する（ステップＳｂ４）。
【０１２０】
パス遡り判定部４０５は、近傍加算比較選択部５１が算出した距離関数ｄｔ（｛μｔ｝）の最小値に基づいて、トレリスのパスを遡って、シンボル判定を行い、判定シンボルＡを求める。パス遡り判定部４０５は、判定シンボルＡを判定結果として出力する（ステップＳｂ５）。パス遡り判定部４０５が順次判定する判定シンボルＡを系列に並べた判定シンボル系列｛Ａｔ｝が、送信シンボルの推定値、すなわち信号生成装置１が取り込んだｍ値データの推定値となる。
【０１２１】
フィルタ更新用伝送路推定部７０の遅延器７２−１〜７２−（ｘ−１）の各々は、パス遡り判定部４０５が順次出力する判定シンボルＡを取り込んで記憶し、「Ｔ」時間経過後に記憶した判定シンボルＡを出力する。タップ７１−１〜７１−ｘの各々には、判定シンボル系列｛Ａｔ｝に含まれるシンボルの各々が与えられる。タップ７１−１〜７１−ｘは、各々の係数値ｈ１〜ｈｘと、与えられたシンボルとを乗算し、乗算結果を加算器７３に出力する。加算器７３は、乗算結果の合計値を算出して出力する。
【０１２２】
遅延器７６は、適応フィルタ部４０１の出力値を「−ｗＴ」の時間、すなわち「−ｗ」シンボル分、遅延させて減算器７７に出力する。減算器７７は、フィルタ更新用伝送路推定部７０の出力値から、「−ｗＴ」時間遅延した適応フィルタ部４０１の出力値を減算し、減算により得られた誤差をフィルタ更新アルゴリズム処理部７５に出力する。
【０１２３】
フィルタ更新アルゴリズム処理部７５は、減算器７７が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムにより係数値ｈ１〜ｈｘの更新値を算出する。フィルタ更新アルゴリズム処理部７５は、算出した係数値ｈ１〜ｈｘを、各々に対応するタップ７１−１〜７１−ｘ、及びする伝送路推定部４０３のタップ６１−１〜６１−ｘに設定して、係数値ｈ１〜ｈｘの更新を行う（ステップＳｂ６）。
【０１２４】
フィルタ更新アルゴリズム処理部７５は、減算器７７が出力する誤差に基づいて、誤差を小さくするように所定の更新アルゴリズムによりタップ利得値ｃ１〜ｃｖの更新値を算出する。フィルタ更新アルゴリズム処理部７５は、算出したタップ利得値ｃ１〜ｃｖを、各々に対応するタップ４３−１〜４３−ｖに設定して、タップ利得値ｃ１〜ｃｖの更新を行う（ステップＳｂ７）。
【０１２５】
適応フィルタ部４０１の遅延器４１が、前回のステップＳｂ１−１において取り込んだ系列長ｖのシンボル系列の先頭から１シンボルずらしたシンボルを先頭として、受信信号系列｛ｒｔ｝から系列長ｖのシンボル系列を取り込める場合（ステップＳｂ８、Ｙｅｓ）、ステップＳｂ１−１、Ｓｂ２−１の処理を行う。一方、遅延器４１が、前回のステップＳｂ１−１において取り込んだ系列長ｖのシンボル系列の先頭から１シンボルずらしたシンボルを先頭として、受信信号系列｛ｒｔ｝から系列長ｖのシンボル系列を取り込めない場合（ステップＳｂ８、Ｎｏ）、処理を終了する。
【０１２６】
なお、上記の図７に示す処理において、ステップＳｂ６とステップＳｂ７の処理の順番が逆になっていてもよい。
【０１２７】
上記の第１の実施形態の構成において、仮判定部３０は、伝送路２から取り込む受信信号系列に対して伝送路２の推定逆伝達関数による適応等化を行ってシンボル系列を生成し、生成したシンボル系列に対して仮判定を行う。伝送路推定部４０３は、伝送路状態を示す複数のシンボル系列と、伝送路２の推定伝達関数とに基づいて、伝送路状態ごとの推定受信シンボル系列を生成する。系列推定アルゴリズム処理部４０２は、受信信号系列から得られるシンボル系列と、推定受信シンボル系列の各々とのメトリックを算出し、算出したメトリックと、仮判定部３０が仮判定した仮判定シンボルと、仮判定シンボルの近傍のシンボルとに基づいて所定の推定アルゴリズムにより最も尤もらしい推定受信シンボル系列を選択する。また、系列推定アルゴリズム処理部４０２は、仮判定部３０が、仮判定した仮判定シンボルと仮判定シンボルの近傍のシンボルの範囲で伝送路状態を示す複数のシンボル系列を生成する。パス遡り判定部４０５は、最も尤もらしい推定受信シンボル系列に基づいて、トレリスのパスを遡って送信シンボル系列の判定を行う。そのため、系列推定アルゴリズム処理部４０２は、仮判定部３０が仮判定した仮判定シンボルと、仮判定シンボルの近傍のシンボルの範囲に絞って、すなわちトレリスの枝数を削減して演算を行えばよくなる。したがって、受信信号系列から送信シンボルの判定を行う際、シンボル多値度が増加した場合であっても、トレリスの枝数を増加させずに、演算量の増大を防ぐことができる。
【０１２８】
ビタビ・アルゴリズムでは、上述したように、伝送路状態μｔに到達する距離関数の最小値ｄ＿ｍｉｎｔ−１（μｔ−１）と、これに対応する全状態遷移が、時刻ｔ−１における全ての伝送路状態μｔ−１において既知であると仮定する。この仮定により、伝送路状態μｔに到達する距離関数ｄｔ（｛μｔ｝）の最小値を求める場合、全ての状態遷移に対応する距離関数ｄｔ（｛μｔ｝）を求める必要がなく、伝送路状態μｔに遷移する可能性のある全ての伝送路状態｛μｔ−１｝について、ｄ＿ｍｉｎｔ−１（μｔ−１）＋ｂ（ｒｔ；μｔ−１→μｔ）を算出すればよいことになる。このように、演算量が少ないビタビ・アルゴリズムであるが、ビタビ・アルゴリズムの演算量は、トレリスの枝数が支配的であり、高多値度符号や、広域な符号間干渉の推定の際には演算量が多くなってしまう。これに対して、上記の第１の実施形態の構成では、近傍加算比較選択部５１は、更に、仮判定部３０が仮判定した仮判定シンボル系列｛Ａ’ｔ｝に含まれるシンボルの各々を中心とする近傍のシンボルの範囲に絞って伝送路状態μｔに到達する距離関数ｄｔ（｛μｔ｝）を算出するため、更に演算量の削減を行うことが可能になる。
【０１２９】
従来のＭＬＳＥでは、例えば、記憶長「３」の場合、トレリスの枝数は、ｍ３本となる。伝送路推定部４０３には、ｍ３個のシンボル系列が与えられるため、伝送路推定部４０３は、ｍ３回の演算を行う必要があり、また、絶対値器５３が算出するメトリックの数もｍ３になる。これに対して、上記の第１の実施形態の構成では、近傍加算比較選択部５１は、仮判定シンボル系列｛Ａ’ｔ｝に含まれるシンボルの各々を中心とする近傍のシンボルの範囲に絞って伝送路状態μｔを表すシンボル系列を生成している。したがって、シンボル多値度ｍが大きくなった場合でも、仮判定シンボル系列｛Ａ’ｔ｝に含まれるシンボルの各々を中心とする近傍のシンボルの範囲を予め制限しておくことで、シンボル多値度の増加による演算量の増加を抑えることができる。
【０１３０】
仮判定部３０には、例えば、上述したように、ＦＦＥが適用されるが、ＦＦＥは、演算量の少ない適応等化技術である。そのため、第１の実施形態におけるシンボル判定装置３のように、最初に演算量の少ないＦＦＥによって、シンボルの仮判定を行い、仮判定結果として得られた仮判定シンボルの各々を中心とする近傍のシンボル範囲に絞る手法の方が、シンボルの推定を全てＭＬＳＥによって行う手法よりも少ない演算量でシンボル推定を行うことができる。
【０１３１】
隣接±１よりも広い符号間干渉（ISI:Inter-Symbol Interference)を推定する場合、近傍のシンボルの範囲として、隣接±ｎ（ただし、ｎは、２以上の整数）にすることになり、この場合、記憶長を「３」とするとトレリスの枝の数は、（２ｎ＋１）３となり、記憶長を「ｑ」とすると、（２ｎ＋１）ｑという一般式で表すことができる。現時点での信号処理チップの製造コストや消費電力の観点からｎ＝１，２程度が現実的な値であるが、今後の半導体加工技術の性能、すなわち信号処理性能の向上によっては、ｎの値を更に大きな値することが可能である。
【０１３２】
なお、図８において、Ａ’ｔ＋１やＡ’ｔ−１などに相当するシンボルが存在しない場合、例えば、［ｉ１，ｉ２，…，ｉｍ］をシフトしたシンボルを用いる。ｉ１−１の場合、ｉｍとなり、ｉｍ＋１の場合ｉ１となる。これらのシンボルは、仮判定シンボルの近傍ではないが、距離関数ｄｔ（｛μｔ｝）が大きくなるため、最終的な判定結果として選択されない。８値データ信号の場合であって、記憶長が「３」、仮判定シンボルとして［５，３，７］が得られた場合を例に具体的に示すと、図９のようなトレリス図となる。
【０１３３】
なお、上記の第１の実施形態の構成において、系列推定部４０が適応フィルタ部４０１を備えない構成であってもよい。この場合、系列推定アルゴリズム処理部４０２の減算器５４は、適応フィルタ部４０１の出力値の代わりに受信信号系列｛ｒｔ｝を取り込むことになる。
【０１３４】
（第２の実施形態）
図１０は、第２の実施形態におけるシンボル判定装置３ａの内部構成を示すブロック図であり、図１１は、第２の実施形態におけるシンボル判定装置３ａの詳細な内部構成を示すブロック図である。第２の実施形態において、第１の実施形態と同一の構成については、同一の符号を付し、以下、異なる構成について説明する。
【０１３５】
図１０に示すように、シンボル判定装置３ａは、仮判定部３０と系列推定部４０ａを備える。系列推定部４０ａは、適応フィルタ部４０１ａ、系列推定アルゴリズム処理部４０２、伝送路推定部４０３、更新処理部４０４、及びパス遡り判定部４０５を備える。
【０１３６】
適応フィルタ部４０１ａの入力端は、仮判定部３０の適応フィルタ部３０１の出力端に接続されている。より詳細には、図１１に示すように、適応フィルタ部３０１の加算器３４の出力端と、適応フィルタ部４０１ａの遅延器４１の入力端とが接続されている。
【０１３７】
第１の実施形態におけるシンボル判定装置３において説明したように、適応フィルタ部３０１は、推定逆伝達関数による適応等化により、受信信号系列｛ｒｔ｝を送信信号系列に近づけることを目的としている。そのため、推定逆伝達関数による適応等化の際に高周波領域の雑音を増幅させることになる。系列推定部４０ａの適応フィルタ部４０１ａでは、適応フィルタ部３０１が増幅させてしまった高周波領域の雑音を抑える程度の少ないタップ数のフィルタリング処理を行うだけで、高周波領域の雑音を抑えるとともにパルス幅の圧縮処理を行うことができる。すなわち、第２の実施形態では、適応フィルタ部３０１の出力信号を利用することで、適応フィルタ部４０１ａにおけるパルス幅の圧縮処理の効率化を図ることができる。
【０１３８】
なお、第２の実施形態においても、第１の実施形態と同様に、系列推定部４０ａが適応フィルタ部４０１ａを備えない構成としてもよい。この場合、系列推定アルゴリズム処理部４０２の減算器５４は、適応フィルタ部４０１ａの出力値の代わりに適応フィルタ部３０１の出力値を取り込むことになる。
【０１３９】
（ＢＣＪＲアルゴリズムを用いる他の構成例）
上記の第１及び第２の実施形態におけるシンボル判定装置３，３ａが備える近傍加算比較選択部５１は、送信シンボル系列を推定する推定アルゴリズムとして、ビタビ・アルゴリズムを実行していたが、本発明の構成は、当該実施の形態に限られない。例えば、近傍加算比較選択部５１が、シンボルごとに異なる事前確率が存在する場合に有効な最大事後確率（Maximum A Posteriori probability：ＭＡＰ）復号法であるＢＣＪＲアルゴリズムを実行するようにしてもよい。なお、ビタビ・アルゴリズムによる最尤復号法は、ＭＡＰ復号法において事前確率が全てのシンボルにおいて等確率にした場合と等価な手法である。
【０１４０】
ＬＤＰＣ(Low Density Parity Check)符号のＢＣＪＲアルゴリズムでは、ＭＬＳＥのビタビ・アルゴリズムと同様に、トレリス図が用いられる。そのため、ＢＣＪＲアルゴリズムを用いる場合、シンボル判定装置３，３ａの構成を変更する必要はない。近傍加算比較選択部５１が実行するアルゴリズムを、ビタビ・アルゴリズムから、ＢＣＪＲアルゴリズムに置き換えることで、ＢＣＪＲアルゴリズムを用いた送信シンボル系列の推定を行うことができる。
【０１４１】
近傍加算比較選択部５１は、仮判定部３０が仮判定した仮判定シンボルと、仮判定シンボルの各々を中心とする近傍のシンボルの範囲でＢＣＪＲアルゴリズムを実行することになる。このとき、シンボルを表すバイナリ値の「０」か「１」のいずれか一方の値の条件付き確率が０％になることがある。ＢＣＪＲアルゴリズムでは、対数尤度比を軟判定出力としており、条件付き確率が０％になると、対数尤度比が無限大になるという問題がある。
【０１４２】
このＢＣＪＲアルゴリズムの問題点を、具体例を用いて説明する。図１２は、シンボル多値度ｍ＝８の場合のｇｒａｙ符号を示した図である。ｇｒａｙ符号では、誤って隣接シンボルとして判定してしまった場合、バイナリ値の誤りを最小限に抑える符号になっている。例えば、「３」と判定すべきところを「２」と判定してしまった場合、バイナリ値では１つ目と２つ目が同じ「０１」を表しているので、バイナリ値の誤りを３つ目の値のみに抑えることができる。
【０１４３】
例えば、記憶長「３」の場合、ＢＣＪＲアルゴリズムの演算対象の伝送路状態として、図１３に示す９通りの状態が存在しているとする。この場合、シンボル「２」，「３」，「４」についてみると、２番目の位置のバイナリ値が必ず「１」になるため、２番目の位置のバイナリ値が「１」になる条件付き確率が１００％になり、「０」になる条件付き確率が０％になる。この場合、対数尤度比Ｒを算出すると、次式（１４）に示すように無限大になる。
【０１４４】
【数１４】
2020184696.tif 000016
【０１４５】
なお、式（１４）において、Ｐ（ｕ＝１）が「１」になる確率であり、Ｐ（ｕ＝０）が「０」になる確率である。対数尤度比Ｒが、無限大ならず、かつ無限大に近い値を示すことができれば、この問題を解決することができる。そこで、対数尤度比Ｒが無限大になる場合、対数尤度比Ｒとして、予め定められる十分大きな値を適用するという手法が考えられる。例えば、Ｐ（ｕ＝１）＝０．９９９９９９にした場合の対数尤度比Ｒは、Ｒ≒６になるため、Ｒ＝１０程度の値を適用する。なお、バイナリ値の「０」と「１」は、単なる表記上の相違であり、「０」と「１」が逆の場合についても同様である。
【０１４６】
バイナリ値の条件付き確率は、仮判定シンボルと、近傍のシンボルによって決まるため、近傍加算比較選択部５１が、仮判定部３０から仮判定シンボルを受けた時点で、対数尤度比Ｒが無限大になるバイナリ値の位置を特定することができる。
【０１４７】
そこで、近傍加算比較選択部５１の内部の記憶領域に、例えば、図１４に示すようなテーブルを予め記憶させておく。図１４に示すテーブルは、「仮判定シンボル値」、「対象シンボル値」、「バイナリ値の位置」の項目を有している。「仮判定シンボル値」の項目には、仮判定シンボルを示す０〜７の値が書き込まれる。「対象シンボル系列」には、記憶長「３」の場合における近傍シンボルを含んだシンボル系列が書き込まれる。
【０１４８】
「バイナリ値の位置」には、「０」か「１」の一方のバイナリ値の条件付き確率が０％になる位置を示す情報が書き込まれる。例えば、仮判定シンボルが「３」の場合、すなわち「２」，「３」，「４」の系列の場合、２番目のバイナリ値が「０」になる条件付き確率が０％になる。これに対して、仮判定シンボルが「１」の場合、すなわち「０」，「１」，「２」の系列の場合、１番目のバイナリ値が「１」になる条件付き確率が０％になる。
【０１４９】
このようなテーブルを備えておくことで、近傍加算比較選択部５１は、仮判定部３０から仮判定シンボルを受けた際、テーブルを参照することにより、対数尤度比Ｒが無限大になるバイナリ値の位置を特定することができる。そして、特定した位置について、対数尤度比Ｒとして、予め定められる十分大きな値、例えば「１０」程度の値を用いて演算を進める。これにより、近傍加算比較選択部５１は、ＢＣＪＲアルゴリズムの演算を問題なく行うことができる。
【０１５０】
なお、上記の第１及び第２の実施形態の構成において、適応フィルタ部３０１、適応フィルタ部４０１，４０１ａ、伝送路推定部４０３、及びフィルタ更新用伝送路推定部７０に、線形トランスバーサルフィルタが適用されるとしているが、本発明の構成は、当該実施の形態に限られない。適応フィルタ部３０１、適応フィルタ部４０１，４０１ａ、伝送路推定部４０３、フィルタ更新用伝送路推定部７０には、他の線形フィルタや非線形フィルタなど、線形トランスバーサルフィルタ以外のフィルタを適用してもよい。
【０１５１】
また、上記の第１及び第２の実施形態の構成において、フィルタ更新アルゴリズム処理部３５，７５が、所定の更新アルゴリズムにより、誤差を小さくするように、タップ利得値ｆ１〜ｆｕ、タップ利得値ｃ１〜ｃｖ、係数値ｈ１〜ｈｘなどのフィルタ係数値の更新値を算出している。ここで、所定の更新アルゴリズムとして、例えば、ＬＭＳ（Least Mean Square）アルゴリズムやＲＬＳ（Recursive Least Square）アルゴリズムなどの反復近似解法のアルゴリズムが適用される。
【０１５２】
ＬＭＳアルゴリズムが適用される場合、各々のタップに設定されるフィルタ係数値の更新値は、例えば、次式（１５）によって算出される。
【０１５３】
Ｈ（ｎ）ｋ＋１＝Ｈ（ｎ）ｋ＋φ・Ｅ・Ｕ（ｎ）ｋ・・・（１５）
【０１５４】
式（１５）において、ｎは、複数あるタップの識別子であり、ｋは、更新の回数を示す値である。また、Ｈ（ｎ）ｋ＋１が、フィルタ係数値の更新値であり、Ｈ（ｎ）ｋが、更新前のフィルタ係数値であり、Ｕ（ｎ）ｋが、更新値を算出する直前にｎ番目のタップに与えられた入力信号である。Ｅが、誤差であり、φが適宜定められる収束定数である。フィルタ更新アルゴリズム処理部３５，７５は、新たなフィルタ係数値を算出するごとに、例えば、算出したフィルタ係数値を内部の記憶領域に書き込んでおき、次回の更新時に、内部の記憶領域から読み出して更新前のフィルタ係数値Ｈ（ｎ）ｋとして用いる。
【０１５５】
上述した実施形態におけるシンボル判定装置３，３ａをコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。
【０１５６】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【０１５７】
３…シンボル判定装置，３０…仮判定部，４０…系列推定部，３０１…適応フィルタ部，３０２…判定処理部，３０３…更新処理部，４０１…適応フィルタ部，４０２…系列推定アルゴリズム処理部，４０３…伝送路推定部，４０４…更新処理部，４０５…パス遡り判定部
【図１】
2020184696.tif 000017
【図２】
2020184696.tif 000018
【図３】
2020184696.tif 000019
【図４】
2020184696.tif 000020
【図５】
2020184696.tif 000021
【図６】
2020184696.tif 000022
【図７】
2020184696.tif 000023
【図８】
2020184696.tif 000024
【図９】
2020184696.tif 000025
【図１０】
2020184696.tif 000026
【図１１】
2020184696.tif 000027
【図１２】
2020184696.tif 000028
【図１３】
2020184696.tif 000029
【図１４】
2020184696.tif 000030
【図１５】
2020184696.tif 000031
【図１６】
2020184696.tif 000032
【図１７】
2020184696.tif 000033
【図１８】
2020184696.tif 000034
【図１９】
2020184696.tif 000035
【図２０】
2020184696.tif 000036
【図２１】
2020184696.tif 000037
【図２２】
2020184696.tif 000038