JP2009225005A - データ処理装置、データ処理方法、及び、プログラム - Google Patents

Abstract

【課題】生き残りステートの数を少なくしつつ、トレリスを用いたデータ系列の推定を、精度良く行う。
【解決手段】ブランチメトリック計算部２０１は、分岐するステートの数を、制限ステート数Mに制限しながら、ブランチメトリックを計算し、ステートメトリック計算部２０２は、生き残りステートの候補のステートメトリックを、ブランチメトリックを用いて計算する。最小値検出部２０３は、生き残りステートの候補のステートメトリックの最小ステートメトリックを検出する。パスメモリ２０８は、最小ステートメトリックのステートから、打ち切り長L-1だけ遡ったステートである有効ステートから分岐するステートを、生き残りステートとして記憶し、有効ステートに対応するデータの候補を、そのデータの推定値として選択する。本発明は、例えば、シンボル系列の推定に適用できる。
【選択図】図１５

Description

本発明は、データ処理装置、データ処理方法、及び、プログラムに関し、特に、例えば、トレリスにおいて生き残る生き残りステートの数を少なくしつつ、トレリスを用いたデータ系列の推定を、精度良く行うことができるようにするデータ処理装置、データ処理方法、及び、プログラムに関する。

通信や放送等の分野で幅広く使用されている符号として、例えば、畳み込み符号がある。畳み込み符号は、符号語を、最も確からしいデータ系列（符号化の対象の情報の系列）に復号する最尤復号を実現することができるという大きな特長がある。最尤復号を実現する復号方法として、ビタビ復号法がある。

ビタビ復号法では、トレリスを用いて、最も確からしいデータ系列が推定（最尤推定）される。

すなわち、ビタビ復号法では、トレリスの各ステートに対して、そのステートを通る最小距離（最小のメトリック）のパス（ステートの系列）の候補を求めて記憶することが繰り返され、無駄な計算を省きながら、最終的に、受信語と最小距離の符号語を選ぶことができる。

ここで、ビタビ復号法には、復号長（符号長）が増大しても、演算量が線形にしか増えないという特長がある。

一方、日本の地上波ディジタルテレビジョン放送の１つの方式であるISDB(Integrated Services Digital Broadcasting)-T等では、OFDM(Orthogonal Frequency Division Multiplexing)（直交周波数分割多重方式）が採用されている。

OFDMは、各シンボルの後ろの部分を、ガードインターバルとして先頭に付加することにより、遅延波（マルチパスのエコー）に対する耐性を有する。但し、ガードインターバルを超える遅延が存在する場合には、シンボル間干渉とキャリア間干渉とが発生し、受信特性（復号特性）が著しく劣化する。

ガードインターバル長を超える遅延波が存在した場合でも、例えば、ベースバンドのOFDM信号に、ハニング窓を掛け、FFT(Fast Fourier Transform)をすることにより、シンボル間干渉を軽減すると同時に、隣接するキャリア以外から発生するキャリア間干渉を軽減することができる。

いま、OFDM信号のキャリア（サブキャリア）を特定するキャリア番号をnと、OFDM信号に、ハニング窓を掛け、FFTをして得られるOFDM信号（以下、干渉低減信号ともいう）をY_nと、それぞれ表すこととすると、干渉低減信号Y_nは、式（１）で表される。

・・・（１）

ここで、キャリア番号がnのキャリアを、キャリア#nと表すこととすると、X_nは、キャリア#nで送信される送信シンボルを表す。

また、N_nは、干渉低減信号Y_nに含まれる雑音成分を表し、シンボル間干渉、キャリア間干渉、及び、加法性白色雑音の和である。さらに、H_nは、キャリア#nに対する伝送路の特性を表し、以下、伝送路情報ともいう。

なお、式（１）において、X_n-1H_n-1と乗算されている係数-1/2、X_nH_nと乗算されている係数+1、及び、X_n+1H_n+1と乗算されている係数-1/2は、ハニング窓の係数である。

また、以下では、キャリア番号#nのキャリア（キャリア#n）は、キャリアの周波数が小さい順の、n番目のキャリアであることとする。

式（１）によれば、干渉低減信号Y_nは、送信シンボルX_nとX_n-1を内部状態（内部ステート）として記憶するメモリ、つまり、２つのメモリを用いた畳み込み符号化により得られる畳み込み符号と見ることができる。したがって、送信シンボルX_n（送信シンボルX_nからなる送信シンボル系列X₀,X₁,・・・）は、伝送路情報H_nが与えられ、また、キャリア番号#nが0以下の範囲の送信シンボルX_nを、例えば、0とすると、送信シンボルX_nは、トレリスを用いた系列の推定を行うことにより推定することができる。

すなわち、送信シンボルX_nは、例えば、ビタビ復号法によって最尤推定することができる。

送信シンボルX_nを、ビタビ復号法によって最尤推定するビタビ復号器において、トレリスのステート数は、OFDM信号の１次変調の信号点（IQコンスタレーション上の信号点）の数の２乗の値となる。

ISDB-Tでは、一次変調として、64QAM(Quadrature Amplitude Modulation),16QAM,QPSK(Quadrature Phase Shift Keying)、又はBPSK(Binary Phase Shift Keying)が行われる。

例えば、１次変調として、64QAMが行われる場合、トレリスのステート数は、4096(=64×64)となる。

また、ビタビ復号器において、トレリスの１個のステートに対して計算する必要があるブランチメトリックの数は、１次変調の信号点の数となる。

したがって、OFDM信号の１つのキャリアについて計算する必要があるブランチメトリックの数は、１次変調の信号点の数の３乗の値となる。

ここで、図１は、一般なビタビ復号法のトレリスを示している。

なお、図１のトレリスでは（後述するトレリスでも同様）、横方向は、時間を表し、縦方向に並ぶ丸印は、各時刻のステートを表す。また、ステートどうしを結ぶ線分が、枝（ブランチ）を表す。

１次変調として、64QAMが行われる場合、トレリスの各時刻のステート数は、上述したように、4096になる。そして、トレリスの各時刻、つまり、OFDMの各キャリアについて計算する必要があるブランチメトリックの数は、262144(=64×64×64)となり、ビタビ復号器を小規模に実装することは、極めて困難になる。

ビタビ復号器等の、データ系列の推定を行う系列推定器の規模の増加を低減する方法としては、例えば、非特許文献１に記載されている方法（以下、Mアルゴリズムともいう）がある。

Mアルゴリズムでは、トレリスにおいて生き残る生き残りステート、つまり、系列推定器で記憶しておくステートを、以下のように抑えることで、系列推定器の規模の増加を低減する。

すなわち、Mアルゴリズムでは、各時刻において、１個の生き残りステートから分岐するステート（枝）の数を、所定の数（以下、制限ステート数ともいう）Mに制限し、各生き残りステートから分岐するステートのステートメトリックをソートする。さらに、Mアルゴリズムでは、ステートメトリックをソートしたステートの中から、ステートメトリックが小さい順に、ある一定の数Cのステートを、新たな生き残りステートとして選択し、以下、同様の処理が繰り返される。

そして、Mアルゴリズムでは、生き残りステートからステートを遡って、系列の推定が行われる。

図２は、Mアルゴリズムを説明する図である。

図２では、制限ステート数Mが2になっており、最新の時刻において選択する生き残りステートの数（以下、最新生き残りステート数ともいう）Cが4になっている。

図２Ａは、Mアルゴリズムによるトレリスを示している。

図２Ａでは、例えば、時刻t=nの生き残りステートとして、４個のステートが存在する。ここで、図２において、白丸印が、生き残りステートを表し、黒丸印が、生き残りステートに選択されなかったステートを表す。

Mアルゴリズムでは、時刻t=nの４個の生き残りステートそれぞれについて、その生き残りステートから分岐する、時刻t=n+1のステートが、制限ステート数Mである２個のステートに制限される。

すなわち、時刻t=nの生き残りステートから分岐する枝のうちの、ブランチメトリックが小さい上位２(=M)本の枝につながる時刻t=n+1のステートが、生き残りステートの候補として求められる。

図２Ａでは、時刻t=nの生き残りステートが４個であり、各生き残りステートについて、時刻t=n+1のステートの２個が、生き残りステートの候補とされるから、時刻t=n+1の８個のステートが、生き残りステートの候補となる。

その後、Mアルゴリズムでは、時刻t=n+1の８個の生き残りステートの候補それぞれのステートメトリックが求められ、さらに、ステートメトリックがソートされる。

さらに、時刻t=n+1の８個の生き残りステートの候補のうちの、ステートメトリックが小さい順に、４(=C)個の生き残りステートの候補が、時刻t=n+1の生き残りステート（新たな生き残りステート）として選択される。

そして、時刻t=n+1の生き残りステートから分岐する枝のうちの、ブランチメトリックが小さい上位２(=M)本の枝につながる時刻t=n+2のステートが、生き残りステートの候補として求められ、以下、同様の処理が繰り返される。

図２Ｂ及び図２Ｃは、Mアルゴリズムにおいて、図２Ａに示したように、生き残りステートが選択される場合に記憶しておく必要があるステートを示している。

すなわち、図２Ｂは、時刻t=n+1の４個の生き残りステートが選択されたときに、その、例えば、2時刻前の時刻t=n-1のステートの中で記憶しておく必要があるステート（生き残りステート）を示している。

時刻t=n-1のステートの中で、時刻t=n+1の４個の生き残りステートにつながるステートは、系列の推定を行うために記憶しておく必要があり、図２Ｂでは、時刻t=n-1のステートのうちの、図中、点線で囲む３個のステートを記憶しておく必要がある。

図２Ｃは、時刻t=n+2の４個の生き残りステートが選択されたときに、その、例えば、2時刻前の時刻t=nのステートの中で記憶しておく必要があるステートを示している。

時刻t=nのステートの中で、時刻t=n+2の４個の生き残りステートにつながるステートは、系列の推定を行うために記憶しておく必要があり、図２Ｃでは、時刻t=nのステートのうちの、図中、点線で囲む４個のステートを記憶しておく必要がある。

"Source and Channel Coding: An Algorithmic Approach", J.B.Anderson, S.Mohan, Kluwer Academic Press, 1991.

Mアルゴリズムでは、過去の時刻のステートの中で記憶しておく必要があるステートのステート数が、例えば、上述のように、３個や４個等に変化する。

したがって、Mアルゴリズムを採用する系列推定器では、過去の時刻のステートの中で記憶しておく必要があるステートのステート数の最大値分の記憶容量のメモリ（いわゆるパスメモリ）を設ける必要があり、系列推定器の規模が大になる。

ここで、図２に示したように、制限ステート数Mが2であり、最新生き残りステート数Cが4である場合において、Mアルゴリズムで記憶しておく必要がある、例えば、2時刻前のステート数の最大値は、4となる。

本発明は、このような状況に鑑みてなされたものであり、生き残りステートの候補の数を適切に制限すること等によって、生き残りステートの数を少なくしつつ、トレリスを用いたデータ系列の推定を、精度良く行うことができるようにするものである。

本発明の第１の側面のデータ処理装置、又は、プログラムは、データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックを計算するブランチメトリック計算手段と、前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックを、前記ブランチメトリックを用いて計算するステートメトリック計算手段と、前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出する検出手段と、前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートを、前記生き残りステートとして記憶する記憶手段と、前記有効ステートに対応する、前記データ系列を構成するデータの候補を、前記データの推定値として選択する選択手段とを備えるデータ処理装置、又は、データ処理装置として、コンピュータを機能させるためのプログラムである。

本発明の第１の側面のデータ処理方法は、データ系列の推定を行うデータ処理装置が、前記データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックを計算し、前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックを、前記ブランチメトリックを用いて計算し、前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出し、前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートを、前記生き残りステートとして記憶し、前記有効ステートに対応する、前記データ系列を構成するデータの候補を、前記データの推定値として選択するステップを含むデータ処理方法である。

以上のような第１の側面においては、前記データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックが計算され、前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックが、前記ブランチメトリックを用いて計算される。さらに、前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックが検出され、前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートが、前記生き残りステートとして記憶される。そして、前記有効ステートに対応する、前記データ系列を構成するデータの候補が、前記データの推定値として選択される。

本発明の第２の側面のデータ処理装置は、データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限し、前記データ系列の推定を打ち切る単位である打ち切り長L分の前記ステートの系列を用いて、前記データ系列を構成するデータの推定値を出力するときに、前記トレリスにおいて生き残る生き残りステートの候補の数を、M^L個に制限するデータ処理装置である。

以上のような第２の側面においては、データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数が、所定の数Mに制限され、前記データ系列の推定を打ち切る単位である打ち切り長L分の前記ステートの系列を用いて、前記データ系列を構成するデータの推定値が出力される。前記トレリスにおいて生き残る生き残りステートの候補の数は、M^L個に制限される。

なお、データ処理装置は、独立した装置であっても良いし、１つの装置を構成している内部ブロックであっても良い。

また、プログラムは、伝送媒体を介して伝送することにより、又は、記録媒体に記録して、提供することができる。

本発明の第１及び第２の側面によれば、データ系列の推定を行うことができ、特に、例えば、トレリスにおいて生き残る生き残りステートの数を少なくしつつ、トレリスを用いたデータ系列の推定を、精度良く行うことができる。

図３は、本発明を適用した受信装置の一実施の形態の構成例を示すブロック図である。

図３において、受信装置は、例えば、ISDB-Tの規格に準拠するOFDM信号を受信して処理する。

すなわち、受信装置は、ハニング窓乗算部１１，FFT部１２、送信系列推定部１３、干渉成分推定部１４、干渉成分除去部１５、誤り訂正部１６、及び、伝送路特性推定部１７から構成される。

ハニング窓乗算部１１には、時間領域の信号である、ベースバンドのOFDM信号が供給される。ハニング窓乗算部１１は、そこに供給されるOFDM信号に、例えば、ハニング窓等の窓関数をかけ、FFT部１２に供給する。

FFT部１２は、ハニング窓乗算部１１からのOFDM信号のFFTを行い、その結果得られる、前述の式（１）で表される干渉低減信号Y_nを、送信系列推定部１３、干渉成分推定部１４、及び、干渉成分除去部１５に供給する。

ここで、ハニング窓は、時間領域、及び周波数領域のうちのいずれの領域でかけてもよい。

すなわち、図３では、時間領域の信号であるOFDM信号に、ハニング窓をかけ、その後、FFTを行うことで、式（１）の干渉低減信号Y_nを得るようにしたが、時間領域の信号であるOFDM信号のFFTを行い、その結果得られる周波数領域のOFDM信号に、ハニング窓をかけることにより、式（１）の干渉低減信号Y_nを得ることが可能である。

いま、時間領域の信号であるOFDM信号のFFTを行うことで周波数領域のOFDM信号を、R_nと表すこととすると、OFDM信号R_nは、式（２）で表される。

・・・（２）

なお、H_ICIは、OFDM信号R_nに含まれるキャリア間干渉を表し、N_ISIは、OFDM信号R_nに含まれるシンボル間干渉を表す。N_AWGNは、OFDM信号R_nに含まれる加法性白色雑音を表す。

以下、キャリア間干渉H_ICI、及び、シンボル間干渉N_ISIの両方を含めて、干渉成分ともいう。

前述したように、ガードインターバルを超える遅延の遅延波（マルチパスのエコー）が存在する場合には、干渉成分の影響が強くなり、受信特性が著しく劣化する。

式（２）のOFDM信号R_nに、ハニング窓をかけることにより、式（１）の干渉低減信号Y_nに含まれる干渉成分を、キャリア#nの両側のキャリア#n-1と#n+1のそれぞれからのシンボル間干渉N_ISI、及び、キャリア間干渉N_ICIに限定し、他のキャリアからのシンボル間干渉N_ISI、及び、キャリア間干渉N_ICIを抑制することができる。

ここで、式（２）の周波数領域のOFDM信号R_nに、ハニング窓をかけて、式（１）の干渉低減信号Y_nを得る処理は、OFDM信号R_nと、ハニング窓の係数である-1/2,1,-1/2とを畳み込むことで行うことができる。

すなわち、OFDM信号にハニング窓をかけることは、時間領域では、OFDM信号とハニング窓の係数との乗算で行うことができ、周波数領域では、OFDM信号とハニング窓の係数との畳み込みで行うことができる。

また、式（１）の干渉低減信号Y_nは、前述したように、畳み込み符号と見ることができ、これにより、送信シンボルX_nは、トレリスを用いた系列の推定を行うことによって推定することができる。

したがって、図３の受信装置では、OFDM信号に、ハニング窓をかけることで、トレリスを用いた送信シンボルX_nのシンボル系列（データ系列）の推定のための、式（１）の拘束条件を成立させているということができる。

送信系列推定部１３は、トレリスを用いたデータ系列の推定を行うデータ処理装置であり、FFT部１２からの畳み込み符号としての干渉低減信号Y_nを用いて、送信シンボル（のシンボル系列）を推定し、その送信シンボルの推定値を、干渉成分推定部１４に供給する。

すなわち、送信系列推定部１３には、FFT部１２から干渉低減信号Y_nが供給される他、伝送路特性推定部１７から、伝送路情報H_n+1が供給されるようになっている。

送信系列推定部１３は、FFT部１２からの干渉低減信号Y_nと、伝送路特性推定部１７からの伝送路情報H_n+1とを用いて、トレリスにより、送信シンボル（データ系列としてのシンボル系列を構成する送信シンボル）を推定する。

具体的には、送信系列推定部１３は、干渉低減信号Y_nを表す式（１）を、その雑音成分N_nを無視して、送信シンボルX_n+1について解いた式（３）に従い、干渉低減信号Y_n、及び、伝送路情報H_n+1とを用いて、送信シンボルX_n+1の受信値X'_n+1を求める。

・・・（３）

そして、送信系列推定部１３は、詳細については後述するが、受信値X'_n+1に基づいて、送信シンボルのシンボル系列の推定に用いるトレリスのステートから分岐するステートの数を、制限ステート数Mに制限し、シンボル系列の推定を打ち切る単位である打ち切り長L分のステートの系列を用いて、シンボル系列を構成する送信シンボルの推定値を出力するが、そのときに、トレリスにおいて生き残る生き残りステートの候補の数を、M^L個に制限する。

干渉成分推定部１４は、FFT部１２からの干渉低減信号Y_nと、送信系列推定部１３からの送信シンボル（のシンボル系列）の推定値とを用いて、干渉低減信号Y_nに含まれる干渉成分を推定し、その干渉成分の推定値を、干渉成分除去部１５に供給する。

干渉成分除去部１５は、FFT部１２からの干渉低減信号Y_nから、干渉成分推定部１４からの干渉成分を除去し、その結果得られる干渉成分除去信号を、誤り訂正部１６に供給する。

誤り訂正部１６は、干渉成分除去部１５からの干渉成分除去信号に対して、誤り訂正処理を施し、その訂正結果である、送信シンボルを構成する情報（ビット）の系列（OFDMによって送信されてきた情報の系列）を出力する。

伝送路特性推定部１７には、OFDMのパイロット信号が供給されるようになっている。伝送路特性推定部１７は、そこに供給されるパイロット信号から、伝送路情報H_n+1を推定し、送信系列推定部１３に供給する。

次に、図４ないし図６を参照して、図３の送信系列推定部１３による、送信シンボルの推定の処理について説明する。

送信系列推定部１３は、シンボル系列の推定に用いるトレリスのステートから分岐するステートの数を、制限ステート数Mに制限しながら、その制限ステート数Mのステートへの枝のブランチメトリックを計算する。

さらに、送信系列推定部１３は、ブランチメトリックを用いて、制限ステート数Mのステートのステートメトリックを計算し、打ち切り長L分のステートの系列を用いて、送信シンボルの推定値を出力する。

送信系列推定部１３は、以上の処理を、FFT部１２から、干渉低減信号Y_nが供給されるたびに繰り返す。

すなわち、図４は（図５、及び図６も同様）、制限ステート数Mを4とするとともに、打ち切り長Lを3とした場合のトレリスを示している。

ここで、打ち切り長Lとは、上述したように、シンボル系列の推定を打ち切る単位であり、幾つの受信値に対して、１個の送信シンボルの推定値を出力するか（1個の送信シンボルの推定に、幾つの受信値を用いるか）を表す。したがって、打ち切り長Lが、例えば、1の場合は、１個の受信値を用いて、１個の送信シンボルの推定値を出力することになる。但し、この場合には、系列の推定にはならない。

また、打ち切り長Lが、例えば、3の場合は、3個の受信値を用いて、つまり、打ち切り長L分のステートの系列を用いて、１個の送信シンボルの推定値を出力することになる。

例えば、いま、時刻t=n+1に注目すると、送信系列推定部１３は、1時刻前の時刻t=nの生き残りステートから分岐するステートの数を、制限ステート数Mである4個に制限して、時刻t=nの生き残りステートそれぞれについて、制限ステート数Mである4個の、時刻t=n+1のステートへのブランチメトリックを計算する。

ここで、時刻t=nの生き残りステートから分岐する時刻t=n+1へのステートを、制限ステート数Mである4個のステートに制限する方法については、後述する。

また、時刻t=n+1に注目したときに、その1時刻前の時刻t=nの送信シンボルX_nの推定値のj番目の候補を、X''_n,i,jと表す。

さらに、送信シンボルX_nの推定値の候補（以下、シンボル候補ともいう）X''_n,i,jに対応する時刻t=nの生き残りステート（以下、生き残りステートX''_n,i,jともいう）につながる、1時刻前の時刻t=n-1のi番目のシンボル候補を、X''_n-1,iと表す。

また、時刻t=nの生き残りステートX''_n,i,jから分岐する時刻t=n+1へのステートに対応する、k番目のシンボル候補を、X''_n+1,i,j,kと表す。

トレリスにおいては、時刻t=n+1のステートX''_n+1,i,j,kは、時刻t=nの生き残りステートX''_n+1,i,jにつながっており、さらに、時刻t=nの生き残りステートX''_n+1,i,jは、時刻t=n-1の生き残りステートX''_n+1,iにつながっている。

時刻t=n+1に注目した場合に、1時刻前の時刻t=nの生き残りステートX''_n+1,i,jから分岐する枝のうちの、ブランチメトリックが計算される枝につながる時刻t=n+1のステートX''_n+1,i,j,kは、以降に、生き残りステートとして生き残るステートの候補である。この生き残りステートとして生き残るステートの候補である、時刻t=n+1のステートX''_n+1,i,j,kを、以下、生き残り候補X''_n+1,i,j,kともいう。

なお、サフィックスi，j、及びkは、いずれも、0,1,2,3(=L-1)の4(=L)個の値をとる。

送信系列推定部１３は、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリックを求めると、そのブランチメトリックを、時刻t=nの生き残りステートX''_n+1,i,jのステートメトリックに加算することで、生き残り候補X''_n+1,i,j,kのステートメトリックを求める。

さらに、送信系列推定部１３は、生き残り候補X''_n+1,i,j,kのステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出する。そして、送信系列推定部１３は、トレリスにおいて、最小ステートメトリックの生き残り候補X''_n+1,i,j,kから、打ち切り長L-1、つまり、いまの場合、2(=3-1)時刻だけ遡った、時刻t=n-1の生き残りステート（以下、有効ステートともいう）X''_n+1,iを検出する。

図５は、有効ステートX''_n+1,iの検出の結果を示している。

図５では、時刻t=n+1の生き残り候補X''_n+1,i,j,k（図４）となっているステートのうちの、ステートS_n+1が最小ステートメトリックのステートになっている。

送信系列推定部１３は、時刻t=n+1の2(=L-1)時刻前の時刻t=n-1の生き残りステートX''_n+1,iのうちの、最小ステートメトリックのステートS_n+1から2時刻だけ遡った生き残りステートS_n-1を、有効ステートとして検出する。

そして、送信系列推定部１３は、トレリスの生き残りステートX''_n+1,iとX''_n+1,i,j、及び、生き残り候補X''_n+1,i,j,kの中で、有効ステートS_n-1から分岐するステート（直接分岐するステートと、間接的に分岐するステート）を、新たな生き残りステートとして記憶し（有効ステートS_n-1から、その有効ステートS_n-1から分岐する最新の時刻のステートまでのパスを記憶し）、他のステートをドロップする（生き残りステートから除外する）。

さらに、送信系列推定部１３は、時刻t=n+1の2(=L-1)時刻前の時刻t=n-1の有効ステートS_n-1に対応するシンボル候補X''_n+1,iを、送信シンボルX_n+1の推定値として選択して出力する（送信シンボルX_n+1の推定を、打ち切り長L分のステートの系列（図５では、時刻tがn-1ないしn+1の3(=L)時刻分のステートの系列）で打ち切り、その推定結果を出力する）。

その後、送信系列推定部１３は、時刻t=n+1の生き残りステートの１個から分岐するステートの数を、制限ステート数Mである4個に制限して、時刻t=n+1の生き残りステートそれぞれについて、制限ステート数Mである4個の、時刻t=n+2のステートへのブランチメトリックを計算し、以下、同様の処理を繰り返す。

図６は、時刻t=n+2に注目した場合のトレリスを示している。

図６においては、時刻t=n-1の有効ステートS_n-1から分岐するステートだけが、生き残りステートになっている。そして、最新の時刻、すなわち、注目している時刻t=n+2のステートのうちの、１時刻前の時刻t=n+1の生き残りステートに直接つながるステートだけが、生き残り候補になる。

つまり、時刻t=n-1の有効ステートS_n-1から、3(=L)時刻だけ後の時刻t=n+2のステートから、生き残り候補が選択される。

上述したように、送信系列推定部１３は、1時刻前の時刻の生き残りステートの１個から分岐するステート（生き残り候補）を、制限ステート数Mである4個に制限する。

したがって、ある時刻tの有効ステートからつながる（有効ステートから直接分岐する）、次の時刻t+1の生き残りステートの数（総数）は、制限ステート数Mになる。

そして、時刻t+1のM個の生き残りステートそれぞれからつながる、次の時刻t+2の生き残りステートの数は、制限ステート数Mに制限されるから、次の時刻t+2の生き残りステートの総数は、制限ステート数Mの2乗であるM²になる。

以下、同様にして、時刻t+Tの生き残りステートの総数は、制限ステート数MのT乗であるM^Tになる。

一方、上述したように、生き残り候補の時刻は、有効ステートの時刻tからL時刻だけ後の時刻t+Lになる。

したがって、生き残り候補の総数は、M^L個、すなわち、例えば、図４ないし図６に示したように、制限ステート数Mが4で、打ち切り長Lが3である場合には、64(=4³)個に制限される。

次に、図７を参照して、時刻t=nの１個の生き残りステートから分岐する時刻t=n+1へのステート（生き残り候補）を、制限ステート数Mのステートに制限する方法について説明する。

時刻t=nの１個の生き残りステートから分岐する時刻t=n+1へのステート（時刻t=n+1への枝）は、１次変調のIQコンスタレーション上の信号点の数だけ存在する。

送信系列推定部１３は、時刻t=nの１個の生き残りステートから分岐する時刻t=n+1への生き残り候補を、１次変調のIQコンスタレーション上の信号点の数のステートのうちの、制限ステート数Mのステートに制限する分岐ステート制限を行う。

分岐ステート制限の方法としては、例えば、時刻t=nの１個の生き残りステートから分岐する時刻t=n+1の、１次変調のIQコンスタレーション上の信号点の数のステートへのブランチメトリックを計算し、そのブランチメトリックが小さい上位M個の枝につながるステートを、生き残り候補として選択する第１の方法を採用することができる。

但し、第１の方法では、１次変調のIQコンスタレーション上の信号点の数のブランチメトリックを計算する必要がある。

そこで、分岐ステート制限の方法としては、時刻t=n+1の生き残り候補を、トレリスを用いたシンボル系列の推定のための式（１）の拘束条件に従って計算されるIQコンスタレーション上の点に近いM個の信号点が表すシンボル候補に対応するステートに制限する第２の方法を採用することもできる。

第２の方法では、送信系列推定部１３は、１次変調のIQコンスタレーション上の信号点の中から、式（１）の拘束条件に従って計算されるIQコンスタレーション上の点、つまり、式（３）に従って計算される、送信シンボルX_n+1の受信値X'_n+1に近い上位M個の信号点を検出し、そのM個の信号点それぞれが表すシンボル候補に対応するステートを、生き残り候補とする。

図７は、制限ステート数Mを4にした場合に、第２の方法によって検出される、送信シンボルX_n+1の受信値X'_n+1に近い上位4(=M)個の信号点を示している。

すなわち、図７Ａは、１次変調としてQPSKを採用した場合に検出される4個の信号点を示している。

QPSKでは、IQコンスタレーション上の信号点の数は、4個であるため、その4個の信号点が、受信値X'_n+1に近い上位4個の信号点として検出される。

図７Ｂは、１次変調として16QAMを採用した場合に検出される4個の信号点を示している。

16QAMでは、IQコンスタレーション上の信号点の数は、16個であるため、その16個の信号点のうちの、受信値X'_n+1に近い上位4個の信号点が検出される。

図７Ｃは、１次変調として64QAMを採用した場合に検出される4個の信号点を示している。

64QAMでは、IQコンスタレーション上の信号点の数は、64個であるため、その64個の信号点のうちの、受信値X'_n+1に近い上位4個の信号点が検出される。

なお、第２の方法によって得られるM個の生き残り候補は、ブランチメトリックが小さい上位M個のステートにほぼ一致する。

また、第１及び第２の方法のいずれを採用する場合であっても、信号点の数が多い、例えば、64QAMでは、多くのステートがドロップされる（生き残り候補から除外される）こととなる。

しかしながら、ここでドロップされるステートは、仮に、生き残り候補として残したとしても、有効ステートとして選択される可能性が極めて低いため、生き残り候補をM個のステートに制限して、多くのステートをドロップすることが、受信特性に与える影響は、極めて小さい。

さらに、受信値X'_n+1に近い上位M個の信号点は、受信値X'_n+1の、IQコンスタレーション上の座標(I,Q)（の範囲）に対して、あらかじめ求めておくことができるので、第２の方法によれば、生き残り候補を、ブランチメトリックを計算する必要がある第１の方法と比較して、極めて少ない処理量で、迅速に求めることができる。

ところで、送信系列推定部１３は、上述したように、トレリスの生き残りステート、及び、生き残り候補の中で、有効ステートから分岐するステートを、新たな生き残りステートとして記憶するため、生き残りステートの数を少なくすることができる。

すなわち、図８は、図２のMアルゴリズムの場合と同様に、制限ステート数Mが2であり、かつ、最新生き残りステート数（最新の時刻の生き残りステートの数）Cが4である場合の、送信系列推定部１３による送信シンボルの推定を表すトレリスを示している。

ここで、送信系列推定部１３による送信シンボルの推定のアルゴリズムを、以下、Xアルゴリズムともいう。

制限ステート数Mが2である場合において、最新生き残りステート数Cが4であることは、Xアルゴリズムでは、打ち切り長Lが3であることに相当する。

したがって、Xアルゴリズムでは、例えば、時刻t=nに注目すると、1時刻前の時刻t=n-1の生き残りステートから分岐するステートの数を、制限ステート数Mである2個に制限して、時刻t=n-1の生き残りステートそれぞれについて、制限ステート数Mである2個の、時刻t=nの生き残り候補へのブランチメトリックが計算される。

さらに、Xアルゴリズムでは、時刻t=nの生き残り候補へのブランチメトリックを用いて、生き残り候補のステートメトリックが求められ、最小ステートメトリックが検出される。そして、トレリスにおいて、最小ステートメトリックの生き残り候補S_nから、打ち切り長L-1、つまり、いまの場合、2(=3-1)時刻だけ遡った、時刻t=n-2の生き残りステートS_n-2が、有効ステートとして検出される。

その後、トレリスの生き残りステート、及び、生き残り候補の中で、有効ステートS_n-2から分岐するステートが、新たな生き残りステートとして記憶される。

ここで、制限ステート数Mは2であるから、有効ステートS_n-2につながる時刻t=n-1の生き残りステートの数は、制限ステート数Mである2個になり、さらに、時刻t=n-1の生き残りステートにつながる時刻t=nの生き残りステートの数は、制限ステート数Mの2乗である4個になる。つまり、最新生き残りステート数Cは、図２のMアルゴリズムの場合と同様に、4になる。

次に、時刻t=n+1においては、Xアルゴリズムでは、やはり、1時刻前の時刻t=nの生き残りステートの１個から分岐するステートの数を、制限ステート数Mである2個に制限して、時刻t=nの4個の生き残りステートそれぞれについて、制限ステート数Mである2個の、時刻t=n+1の生き残り候補へのブランチメトリックが計算される。

さらに、Xアルゴリズムでは、時刻t=nの生き残り候補へのブランチメトリックを用いて、生き残り候補のステートメトリックが求められ、最小ステートメトリックが検出される。そして、トレリスにおいて、最小ステートメトリックの生き残り候補S_n+1から、打ち切り長L-1、つまり、いまの場合、2(=3-1)時刻だけ遡った、時刻t=n-1の生き残りステートS_n-1が、有効ステートとして検出される。

その後、トレリスの生き残りステート、及び、生き残り候補の中で、有効ステートS_n-1から分岐するステートが、新たな生き残りステートとして記憶される。

ここで、制限ステート数Mは2であるから、有効ステートS_n-1につながる時刻t=nの生き残りステートの数は、制限ステート数Mである2個になり、さらに、時刻t=nの生き残りステートにつながる時刻t=n+1の生き残りステートの数は、制限ステート数Mの2乗である4個になる。つまり、最新生き残りステート数Cは、図２のMアルゴリズムの場合と同様に、4になる。

以下、Xアルゴリズムでは、同様の処理が繰り返される。

Mアルゴリズムでは、図２で説明したように、過去の時刻のステートの中で記憶しておく必要があるステートの数は変化する。そして、制限ステート数Mが2であり、最新生き残りステート数Cが4である場合においては、Mアルゴリズムで記憶しておく必要がある、例えば、2時刻前のステート数の最大値は、4となる。

一方、Xアルゴリズムでは、過去の時刻のステートの中で記憶しておく必要があるステートの数は、時刻ごとに固定である。例えば、制限ステート数Mが2であり、最新生き残りステート数Cが4である場合において記憶しておく必要がある、例えば、2時刻前のステート数は、2となる。

具体的には、Xアルゴリズムでは、例えば、時刻t=n+1の生き残り候補からは、有効ステートS_n-1に（間接的に）つながる4個のステートが、新たな生き残りステートとして選択されるが、そのときには、1時刻前の時刻t=nで得られた有効ステートS_n-2につながるステートを、生き残りステートとして記憶しておく必要がある。

時刻t=n+1の生き残り候補から、4個のステートが、新たな生き残りステートとして選択されるときの2時刻前の時刻t=n-1のステートの中で、有効ステートS_n-2につながるステートの数は、図８において点線で囲むように、2個であり、この数は、制限ステート数Mが2であり、かつ、打ち切り長Lが3である限り、変わらない。

したがって、Xアルゴリズムでは、生き残りステートの数を少なくすることができる。すなわち、Mアルゴリズムでは、（最大で）４個のステートを記憶しておく必要があるのに対して、Xアルゴリズムでは、２個のステートを記憶しておくだけで済む。

その結果、Xアルゴリズムによれば、ステートを記憶するメモリの記憶容量が少なくて済み、送信系列推定部１３、ひいては、図３の受信装置の規模を、Mアルゴリズムの場合より低減することができる。

また、Xアルゴリズムでは、生き残り候補の数を、制限ステート数Mに制限するが、上述したように、生き残り候補を制限することが受信特性に与える影響は、極めて小さいので、トレリスを用いたシンボル系列の推定を、精度良く行うことができる。

さらに、Xアルゴリズムでは、生き残り候補のステートメトリックの中から、最小ステートメトリックを検出するだけで良く、図２で説明したように、ステートメトリックをソートする必要があるMアルゴリズムに比較して、シンボル系列の推定を、迅速に行うことができる。

次に、図９ないし図１４を参照して、送信系列推定部１３（Xアルゴリズム）による、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理について説明する。

なお、以下では、制限ステート数Mを、例えば、4とするとともに、打ち切り長Lを、例えば、3とする。

また、ISDB-Tの規格に準拠するOFDM信号は、MODE1，MODE2、及びMODE3が規定されているが、ここでは、MODE1，MODE2、及びMODE3のうちの、例えば、MODE3について説明する。

MODE3では、１シンボルのOFDM信号は、5617のキャリア（サブキャリア）#1ないし#5617を有する。

図３の受信装置では、FFT部１２から、送信系列推定部１３に対して、キャリア#0ないし#5616の干渉低減信号Y_nが、例えば、周波数の低い順に供給されるとともに、FFT部１２から、送信系列推定部１３に対して、キャリア#nの干渉低減信号Y_nが供給されるタイミングで、伝送路特性推定部１７から、送信系列推定部１３に対して、伝送路情報H_n+1が供給される。

図９は、FFT部１２から、送信系列推定部１３に対して、干渉低減信号Y₀が供給されるとともに、伝送路特性推定部１７から、送信系列推定部１３に対して、伝送路情報H₁が供給されたときの、送信系列推定部１３の処理を説明する図である。

すなわち、図９Ａは、１番目のキャリア#1が、キャリア番号がn+1のキャリア（式（３）の送信シンボルX_n+1を送信するキャリア）になる場合のキャリアの並びを示している。

送信系列推定部１３では、１番目のキャリア#1(=n+1)による送信シンボルX₁(=X_n+1)の受信値X'₁(=X'_n+1)が、式（３）に従い、干渉低減信号Y₀(=Y_n)、送信シンボルX₀(=X_n)の推定値X''₀(=X''_n)、及び送信シンボルX_-1(=X_n-1)の推定値X''_-1(=X''_n-1)、並びに、伝送路情報H₁(H_n+1)，H₀(=H_n)、及びH_-1(=H_n-1)を用いて求められる。

ここで、１シンボルのOFDM信号に、キャリア番号#nが0以下のキャリア、及び、5618以上のキャリアは存在しないため、キャリア番号#nが0以下のキャリア、及び、5618以上のキャリアの送信シンボルX_n、及び、その推定値X_n''は、0とする。

したがって、推定値X''₀(=X''_n)、及びX''_-1(=X''_n-1)は、0であるから、受信値X'₁(=X'_n+1)は、結局、式（３）のX_n、及びX_n-1を0とした式（４）に従って求められることになる。

・・・（４）

送信系列推定部１３は、受信値X'₁を求めた後、IQコンスタレーション上で、ユークリッド距離が受信値X'₁により近い、制限ステート数Mである4個の信号点を選択する。さらに、送信系列推定部１３は、受信値X'₁により近い4個の信号点それぞれが表す4個のシンボルを、送信シンボルX₁の推定値の候補（シンボル候補）X''_1,i=X''_1,0，X''_1,1，X''_1,2，X''_1,3とする。

その後、送信系列推定部１３は、4個のシンボル候補X''_1,iに対応する4個のステートX''_1,iを、時刻t=1(=n+1)の生き残りステートの候補（生き残り候補）として、1時刻前の時刻t=0(=n)のステートから、時刻t=1(=n+1)の生き残り候補X''_1,iへのブランチメトリックを計算する。

ここで、時刻t=nのステートX''_n（推定値X''_nに対応するステート）から、時刻t=n+1のステートX''_n+1（推定値X''_n+1に対応するステート）へのブランチメトリックをBMで表すこととすると、ブランチメトリックBMは、例えば、式（５）に従って計算される。

・・・（５）

式（５）によれば、ブランチメトリックBMは、干渉低減信号Y_nに含まれるキャリア間干渉、シンボル間干渉、及び加法性白色雑音の和に相当する。送信系列推定部１３によるシンボル系列の推定は、ブランチメトリックBM（の和）を最小にするように行われる。

図９Ｂは、シンボル候補X''_1,iに対応するステートX''_1,iが、生き残り候補とされた状態のトレリスを示している。

いま、時刻t=1(=n+1)の生き残り候補X''_1,iへのブランチメトリックBMを、BM_iと表すこととすると、送信系列推定部１３は、時刻t=1(=n+1)の4個の生き残り候補X''_1,iのステートメトリックsm_iを、ブランチメトリックBM_iとする。

そして、時刻t=1(=n+1)では、その時刻t=1の生き残り候補X''_1,iから、打ち切り長Lである3だけ遡ったステートがないため、送信系列推定部１３は、4個の生き残り候補X''_1,iのすべてを、生き残りステートX''_1,iとして選択する。

図１０は、FFT部１２から、送信系列推定部１３に対して、干渉低減信号Y₁が供給されるとともに、伝送路特性推定部１７から、送信系列推定部１３に対して、伝送路情報H₂が供給されたときの、送信系列推定部１３の処理を説明する図である。

すなわち、図１０Ａは、２番目のキャリア#2が、キャリア番号がn+1のキャリアになる場合のキャリアの並びを示している。

送信系列推定部１３では、２番目のキャリア#2(=n+1)による送信シンボルX₂(=X_n+1)の受信値X'₂(=X'_n+1)が、式（３）に従い、干渉低減信号Y₁(=Y_n)、送信シンボルX₁(=X_n)の推定値X''₁(=X''_n)、及び送信シンボルX₀(=X_n-1)の推定値X''₀(=X''_n-1)、並びに、伝送路情報H₂(H_n+1)，H₁(=H_n)、及びH₀(=H_n-1)を用いて求められる。

ここで、図９で説明したように、キャリア番号#nが0以下のキャリア、及び、5618以上のキャリアの送信シンボルX_nの推定値X_n''は、0であり、したがって、推定値X''₀(=X''_n-1)は、0であるから、受信値X'₂(=X'_n+1)は、結局、式（３）のX_n-1を0とした式（６）に従って求められることになる。

・・・（６）

なお、式（６）において、受信値X'_2,iは、生き残りステートX''_1,iにつながるステートから、生き残り候補を選択するのに用いられる受信値X'₂を表す。

生き残りステートX''_1,iとしては、ステートX''_1,0，X''_1,1，X''_1,2、及び、X''_1,3の４個のステートがあるので、その4個の生き残りステートX''_1,iについて、4個の受信値X'_2,iが、それぞれ求められる。

送信系列推定部１３は、4個の受信値X'_2,iを求めた後、IQコンスタレーション上で、その4個の受信値X'_2,iそれぞれについて、受信値X'_2,iにより近い、制限ステート数Mである4個の信号点を選択する。さらに、送信系列推定部１３は、受信値X'_2,iにより近い4個の信号点それぞれが表す4個のシンボルを、送信シンボルX₂のシンボル候補X''_2,i,j=X''_2,i,0，X''_2,i,1，X''_2,i,2，X''_2,i,3とする。

ここで、4個のシンボル候補X''_2,i,jは、4個の受信値X'_2,iそれぞれについて求められるため、全部で、16個のシンボル候補X''_2,i,jが求められる。

その後、送信系列推定部１３は、16個のシンボル候補X''_2,i,jに対応する16個のステートX''_2,i,jを、時刻t=2(=n+1)の生き残り候補として、1時刻前の時刻t=1(=n)の生き残りステートX''_1,iから、時刻t=2の生き残り候補X''_2,i,jへのブランチメトリックを、上述の式（５）に従って計算する。

図１０Ｂは、16個のシンボル候補X''_2,i,jに対応する16個のステートX''_2,i,jが、生き残り候補とされた状態のトレリスを示している。

時刻t=2(=n+1)では、時刻t=1(=n)の生き残りステートX''_1,iの１個から、制限ステート数Mである4個の生き残り候補X''_2,i,0，X''_2,i,1，X''_2,i,2，X''_2,i,3が分岐している。

いま、時刻t=2(=n+1)の生き残り候補X''_2,i,jへのブランチメトリックBMを、BM_i,jと表すこととすると、送信系列推定部１３は、直前の時刻t=1(=n)の生き残りステートX''_1,iのステートメトリックsm_iに、ブランチメトリックBM_i,jを加算することで、時刻t=2(=n+1)の16個の生き残り候補X''_2,i,jのステートメトリックsm_i,jを求める。

そして、時刻t=2(=n+1)では、その時刻t=2の生き残り候補X''_2,i,jから、打ち切り長Lである3だけ遡ったステートがないため、送信系列推定部１３は、16個の生き残り候補X''_2,i,jのすべてを、生き残りステートX''_2,i,jとして選択する。

図１１は、FFT部１２から、送信系列推定部１３に対して、干渉低減信号Y₂が供給されるとともに、伝送路特性推定部１７から、送信系列推定部１３に対して、伝送路情報H₃が供給されたときの、送信系列推定部１３の処理を説明する図である。

すなわち、図１１Ａは、３番目のキャリア#3が、キャリア番号がn+1のキャリアになる場合のキャリアの並びを示している。

送信系列推定部１３では、３番目のキャリア#3(=n+1)による送信シンボルX₃(=X_n+1)の受信値X'₃(=X'_n+1)が、式（３）に従い、干渉低減信号Y₂(=Y_n)、送信シンボルX₂(=X_n)の推定値X''₂(=X''_n)、及び送信シンボルX₁(=X_n-1)の推定値X''₁(=X''_n-1)、並びに、伝送路情報H₃(H_n+1)，H₂(=H_n)、及びH₁(=H_n-1)を用いて求められる。

すなわち、送信系列推定部１３は、直前の時刻t=2(=n)の16個の生き残りステートX''_2,i,jのそれぞれについて、１個ずつの受信値X'_3,i,j(=X'_n+1,i,j)、つまり、合計で、16個の受信値X'_3,i,jを、式（３）に相当する式（７）に従って求める。

・・・（７）

送信系列推定部１３は、16個の受信値X'_3,i,jを求めた後、IQコンスタレーション上で、その16個の受信値X'_3,i,jそれぞれについて、受信値X'_3,i,jにより近い、制限ステート数Mである4個の信号点を選択する。さらに、送信系列推定部１３は、受信値X'_3,i,jにより近い4個の信号点それぞれが表す4個のシンボルを、送信シンボルX₃のシンボル候補X''_3,i,j,k=X''_3,i,j,0，X''_3,i,j,1，X''_3,i,j,2，X''_3,i,j,3とする。

ここで、4個のシンボル候補X''_3,i,j,kは、16個の受信値X'_3,i,jそれぞれについて求められるため、全部で、64個のシンボル候補X''_3,i,j,kが求められる。

その後、送信系列推定部１３は、64個のシンボル候補X''_3,i,j,kに対応する64個のステートX''_3,i,j,kを、時刻t=3(=n+1)の生き残り候補として、1時刻前の時刻t=2(=n)の生き残りステートX''_2,i,jから、時刻t=3の生き残り候補X''_3,i,j,kへのブランチメトリックを、上述の式（５）に従って計算する。

図１１Ｂは、64個のシンボル候補X''_3,i,j,kに対応する64個のステートX''_3,i,j,kが、生き残り候補とされた状態のトレリスを示している。

時刻t=3(=n+1)では、時刻t=2(=n)の生き残りステートX''_2,i,jの１個から、制限ステート数Mである4個の生き残り候補X''_3,i,j,0，X''_3,i,j,1，X''_3,i,j,2，X''_3,i,j,3が分岐している。

いま、時刻t=3(=n+1)の生き残り候補X''_3,i,j,kへのブランチメトリックBMを、BM_i,j,kと表すこととすると、送信系列推定部１３は、直前の時刻t=2(=n)の生き残りステートX''_2,i,jのステートメトリックsm_i,jに、ブランチメトリックBM_i,j,kを加算することで、時刻t=3(=n+1)の64(=M^L)個の生き残り候補X''_3,i,j,kのステートメトリックsm_i,j,kを求める。

その後、送信系列推定部１３は、64(=M^L)個の生き残り候補X''_3,i,j,kのステートメトリックsm_i,j,kのうちの最小ステートメトリックを検出する。さらに、送信系列推定部１３は、最小ステートメトリックの生き残り候補X''_3,i,j,kから、打ち切り長L-1である2だけ遡った時刻t=1(=(n+1)-(L-1))の生き残りステートX''_1,iを、有効ステートとして検出する。

そして、送信系列推定部１３は、有効ステートX''_1,iに対して、直接及び間接的につながる生き残りステートと、生き残り候補とを、新たな生き残りステートとして選択する。

さらに、送信系列推定部１３は、有効ステートX''_1,iに対応するシンボル候補X''_1,iを、キャリア#1(=n-1=n+1-(L-1))による送信シンボルX₁(=X_n-1)の推定値として選択する。

したがって、図１１Ｂにおいて、例えば、最小ステートメトリックの生き残り候補X''_3,i,j,kが、生き残り候補X''_3,0,0,0である場合には、その生き残り候補X''_3,0,0,0から2(=L-1)だけ遡った時刻t=1の生き残りステートX''_1,0が有効ステートとして選択される。

そして、有効ステートX''_1,0、有効ステートX''_1,0に直接つながる時刻t=2の4個の生き残りステートX''_2,0,j、有効ステートX''_1,0に間接的につながる（時刻t=2の4個の生き残りステートX''_2,0,jに直接つながる）時刻t=3の16個の生き残り候補X''_3,0,j,kが、新たな生き残りステートとして選択される。

さらに、有効ステートX''_1,0に対応するシンボル候補X''_1,0が、キャリア#1(=n-1=n+1-(L-1))による送信シンボルX₁(=X_n-1)の推定値として選択される。

図１２は、FFT部１２から、送信系列推定部１３に対して、キャリア#nの干渉低減信号Y_nが供給されるとともに、伝送路特性推定部１７から、送信系列推定部１３に対して、伝送路情報H_n+1が供給されたときの、送信系列推定部１３の処理を説明する図である。

キャリア#4ないし#5617が、キャリア番号がn+1のキャリアになる場合においては、送信系列推定部１３では、キャリア#n+1による送信シンボルX_n+1の受信値X'_n+1が、式（３）に従い、干渉低減信号Y_n、送信シンボルX_nの推定値X''_n、及び送信シンボルX_n-1の推定値X''_n-1、並びに、伝送路情報H_n+1，H_n、及びH_n-1を用いて求められる。

すなわち、送信系列推定部１３は、直前の時刻t=nの16個の生き残りステートX''_n,i,jのそれぞれについて、１個ずつの受信値X'_n+1,i,j、つまり、合計で、16個の受信値X'_n+1,i,jを求める。

送信系列推定部１３は、16個の受信値X'_n+1,i,jを求めた後、IQコンスタレーション上で、その16個の受信値X'_n+1,i,jそれぞれについて、受信値X'_n+1,i,jにより近い、制限ステート数Mである4個の信号点を選択する。さらに、送信系列推定部１３は、受信値X'_n+1,i,jにより近い4個の信号点それぞれが表す4個のシンボルを、送信シンボルX_n+1のシンボル候補X''_n+1,i,j,k=X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3とする。

ここで、4個のシンボル候補X''_n+1,i,j,kは、16個の受信値X'_n+1,i,jそれぞれについて求められるため、全部で、64個のシンボル候補X''_n+1,i,j,kが求められる。

その後、送信系列推定部１３は、64個のシンボル候補X''_n+1,i,j,kに対応する64個のステートX''_n+1,i,j,kを、時刻t=n+1の生き残り候補として、1時刻前の時刻t=nの生き残りステートX''_n,i,jから、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリックを、上述の式（５）に従って計算する。

図１２は、64個のシンボル候補X''_n+1,i,j,kに対応する64個のステートX''n+1_,i,j,kが、生き残り候補とされた状態のトレリスを示している。

時刻t=n+1では、直前の時刻t=nの生き残りステートX''_n,i,jの１個から、制限ステート数Mである4個の生き残り候補X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3が分岐している。

いま、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリックBMを、BM_i,j,kと表すこととすると、送信系列推定部１３は、直前の時刻t=nの生き残りステートX''_n,i,jのステートメトリックsm_i,jに、ブランチメトリックBM_i,j,kを加算することで、時刻t=n+1の64(=M^L)個の生き残り候補X''_n+1,i,j,kのステートメトリックsm_i,j,kを求める。

その後、送信系列推定部１３は、64(=M^L)個の生き残り候補X''_n+1,i,j,kのステートメトリックsm_i,j,kのうちの最小ステートメトリックを検出する。さらに、送信系列推定部１３は、最小ステートメトリックの生き残り候補X''_n+1,i,j,kから、打ち切り長L-1である2だけ遡った時刻t=n-1(=(n+1)-(L-1))の生き残りステートX''_n-1,iを、有効ステートとして検出する。

そして、送信系列推定部１３は、有効ステートX''_n-1,iに対して、直接及び間接的につながる生き残りステートと、生き残り候補とを、新たな生き残りステートとして選択する。

さらに、送信系列推定部１３は、有効ステートX''_n-1,iに対応するシンボル候補X''_n-1,iを、キャリア#n-1(=n+1-(L-1))による送信シンボルX_n-1の推定値として選択する。

したがって、図１２において、例えば、最小ステートメトリックの生き残り候補X''_n+1,i,j,kが、生き残り候補X''_n+1,0,0,0である場合には、その生き残り候補X''_n+1,0,0,0から2(=L-1)だけ遡った時刻t=n-1の生き残りステートX''_n-1,0が有効ステートとして選択される。

そして、有効ステートX''_n-1,0に直接つながる時刻t=nの4個の生き残りステートX''_n,0,j、有効ステートX''_n-1,0に間接的につながる（時刻t=nの4個の生き残りステートX''_n,0,jに直接つながる）時刻t=n+1の16個の生き残り候補X''_n+1,0,j,kが、新たな生き残りステートとして選択される。

さらに、有効ステートX''_n-1,0に対応するシンボル候補X''_n-1,0が、キャリア#n-1による送信シンボルX_n-1の推定値として選択される。

以上のように、送信系列推定部１３では、トレリスのステート数にかかわらず、生き残り候補の数を、一定（M^L個）に保ちながら、シンボル系列の推定を、容易に行うことができる。

ところで、送信系列推定部１３では、時刻t=n+1において、キャリア#n+1による送信シンボルX_n+1のシンボル候補X''_n+1,i,j,kに対応するM^L個のステートX''_n+1,i,j,kが、生き残り候補として求められた後、その生き残り候補X''_n+1,i,j,kのうちの、最小ステートメトリックの生き残り候補X''_n+1,i,j,kから、打ち切り長L-1である2だけ遡った時刻t=n-1(=(n+1)-(L-1))の生き残りステートX''_n-1,iが、有効ステートとして検出される。

そして、有効ステートX''_n-1,iに対応するシンボル候補X''_n-1,iが、2(=L-1)時刻前の時刻t=n-1のキャリア#n-1による送信シンボルX_n-1の推定値として選択される。

したがって、送信系列推定部１３では、キャリア#n-1による送信シンボルX_n-1の推定値を求めるには、その2時刻後の時刻t=n+1の生き残り候補X''_n+1,i,j,kまで求める必要がある。

そのため、5617のキャリア#1ないし#5617を有する、１シンボルのOFDM信号については、最後のキャリア#5617の２つ後の（存在しない）キャリア#5619が、キャリア番号がn+1のキャリアになるまで、シンボル系列の推定の処理を行う必要がある。

図１３は、１シンボルのOFDM信号の最後のキャリア#5617の１つ後の（存在しない）キャリア#5618が、キャリア番号がn+1のキャリアになる場合の、送信系列推定部１３の処理を説明する図である。

すなわち、図１３Ａは、5618番目のキャリア#5618が、キャリア番号がn+1のキャリアになる場合のキャリアの並びを示している。

送信系列推定部１３では、直前の時刻t=5617(=n)の16個の生き残りステートX''_5617,i,j(=X''_n,i,j)のそれぞれについて、5618番目のキャリア#5618(=n+1)による送信シンボルX₅₆₁₈(=X_n+1)の、１個ずつのシンボル候補X''_5618,i,j,0が求められる。

すなわち、送信系列推定部１３では、16個の生き残りステートX''_5617,i,j(=X''_n,i,j)について、16個のシンボル候補X''_5618,i,j,0が、それぞれ求められる。

ここで、図９で説明したように、キャリア番号#nが5618以上のキャリアの送信シンボルX_nの推定値X_n''は、0であり、したがって、シンボル候補X''_5618,i,j,0としては、0が求められる。

送信系列推定部１３は、16個のシンボル候補X''_5618,i,j,0に対応する16個のステートX''_5618,i,j,0を、時刻t=5618(=n+1)の生き残り候補として、1時刻前の時刻t=5617(=n)の生き残りステートX''_5617,i,jから、時刻t=5618の生き残り候補X''_5618,i,j,0へのブランチメトリックを、上述の式（５）に従って計算する。

図１３Ｂは、16個のシンボル候補X''_5618,i,j,0に対応する16個のステートX''_5618,i,j,0が、生き残り候補とされた状態のトレリスを示している。

時刻t=5618(=n+1)では、時刻t=5617(=n)の生き残りステートX''_5617,i,jの１個から、１個の生き残り候補X''_5618,i,j,0が分岐している。

いま、時刻t=5618(=n+1)の生き残り候補X''_5618,i,j,0へのブランチメトリックBMを、BM_i,j,0と表すこととすると、送信系列推定部１３は、直前の時刻t=5617(=n)の生き残りステートX''_5617,i,jのステートメトリックsm_i,jに、ブランチメトリックBM_i,j,0を加算することで、時刻t=5618(=n+1)の16個の生き残り候補X''_5618,i,j,0のステートメトリックsm_i,j,0を求める。

その後、送信系列推定部１３は、16個の生き残り候補X''_5618,i,j,0のステートメトリックsm_i,j,0のうちの最小ステートメトリックを検出する。さらに、送信系列推定部１３は、最小ステートメトリックの生き残り候補X''_5618,i,j,0から、打ち切り長L-1である2だけ遡った時刻t=5616(=(n+1)-(L-1))の生き残りステートX''_5616,iを、有効ステートとして検出する。

そして、送信系列推定部１３は、有効ステートX''_561,iに対して、直接及び間接的につながる生き残りステートと、生き残り候補とを、新たな生き残りステートとして選択する。

さらに、送信系列推定部１３は、有効ステートX''_5616,iに対応するシンボル候補X''_5616,iを、キャリア#5616(=n-1)による送信シンボルX₅₆₁₆(=X_n-1)の推定値として選択する。

したがって、図１３Ｂにおいて、例えば、最小ステートメトリックの生き残り候補X''_5618,i,j,0が、生き残り候補X''_5618,0,0,0である場合には、その生き残り候補X''_5618,0,0,0から2(=L-1)だけ遡った時刻t=5616の生き残りステートX''_5616,0が有効ステートとして選択される。

そして、有効ステートX''_5616,0に直接つながる時刻t=5617の4個の生き残りステートX''_5617,0,j、有効ステートX''_5616,0に間接的につながる（時刻t=5617の4個の生き残りステートX''_5617,0,jに直接つながる）時刻t=5618の4個の生き残り候補X''_5618,0,j,kが、新たな生き残りステートとして選択される。

さらに、有効ステートX''_5616,0に対応するシンボル候補X''_5616,0が、キャリア#5616(=n-1)による送信シンボルX₅₆₁₆(=X_n-1)の推定値として選択される。

図１４は、１シンボルのOFDM信号の最後のキャリア#5617の２つ後の（存在しない）キャリア#5619が、キャリア番号がn+1のキャリアになる場合の、送信系列推定部１３の処理を説明する図である。

すなわち、図１４Ａは、5619番目のキャリア#5619が、キャリア番号がn+1のキャリアになる場合のキャリアの並びを示している。

送信系列推定部１３では、直前の時刻t=5618(=n)の4個の生き残りステートX''_5618,i,0(=X''_n,i,0)のそれぞれについて、5619番目のキャリア#5619(=n+1)による送信シンボルX₅₆₁₉(=X_n+1)の、１個ずつのシンボル候補X''_5619,i,0,0が求められる。

すなわち、送信系列推定部１３では、4個の生き残りステートX''_5618,i,0(=X''_n,i,0)について、4個のシンボル候補X''_5619,i,0,0が、それぞれ求められる。

ここで、図９で説明したように、キャリア番号#nが5618以上のキャリアの送信シンボルX_nの推定値X_n''は、0であり、したがって、シンボル候補X''_5619,i,0,0としては、0が求められる。

送信系列推定部１３は、4個のシンボル候補X''_5619,i,0,0に対応する4個のステートX''_5619,i,0,0を、時刻t=5619(=n+1)の生き残り候補として、1時刻前の時刻t=5618(=n)の生き残りステートX''_5618,i,0から、時刻t=5619の生き残り候補X''_5619,i,0,0へのブランチメトリックを、上述の式（５）に従って計算する。

図１４Ｂは、4個のシンボル候補X''_5619,i,0,0に対応する4個のステートX''_5619,i,0,0が、生き残り候補とされた状態のトレリスを示している。

時刻t=5619(=n+1)では、時刻t=5618(=n)の生き残りステートX''_5618,i,0の１個から、１個の生き残り候補X''_5619,i,0,0が分岐している。

いま、時刻t=5619(=n+1)の生き残り候補X''_5619,i,0,0へのブランチメトリックBMを、BM_i,0,0と表すこととすると、送信系列推定部１３は、直前の時刻t=5618(=n)の生き残りステートX''_5618,i,0のステートメトリックsm_i,0に、ブランチメトリックBM_i,0,0を加算することで、時刻t=5619(=n+1)の4個の生き残り候補X''_5619,i,0,0のステートメトリックsm_i,0,0を求める。

その後、送信系列推定部１３は、4個の生き残り候補X''_5619,i,0,0のステートメトリックsm_i,0,0のうちの最小ステートメトリックを検出する。さらに、送信系列推定部１３は、最小ステートメトリックの生き残り候補X''_5619,i,0,0から、打ち切り長L-1である2だけ遡った時刻t=5617(=(n+1)-(L-1))の生き残りステートX''_5617,iを、有効ステートとして検出する。

そして、送信系列推定部１３は、有効ステートX''_5617,iに対応するシンボル候補X''_5617,iを、１シンボルのOFDM信号の最後のキャリア#5617(=n-1)による送信シンボルX₅₆₁₇(=X_n-1)の推定値として選択する。

したがって、図１４Ｂにおいて、例えば、最小ステートメトリックの生き残り候補X''_5619,i,0,0が、生き残り候補X''_5619,0,0,0である場合には、その生き残り候補X''_5619,0,0,0から2(=L-1)だけ遡った時刻t=5617の生き残りステートX''_5617,0が有効ステートとして選択される。

そして、有効ステートX''_5617,0に対応するシンボル候補X''_5617,0が、最後のキャリア#5617(=n-1)による送信シンボルX₅₆₁₇(=X_n-1)の推定値として選択される。

図１５は、図３の送信系列推定部１３の構成例を示している。

図１５では、送信系列推定部１３は、ブランチメトリック計算部２０１、ステートメトリック計算部２０２、最小値検出部２０３、セレクタ２０４、レジスタ(REG)２０５ないし２０７、及び、パスメモリ２０８から構成される。

ブランチメトリック計算部２０１には、FFT部１２（図３）から、干渉低減信号Y_n(D206)が供給されるとともに、伝送路特性推定部１７（図３）から、伝送路情報H_n+1が供給される。

さらに、ブランチメトリック計算部２０１には、レジスタ２０６から、伝送路情報H_n(D202)が供給されるとともに、レジスタ２０７から、伝送路情報H_n-1(D201)が供給される。

また、ブランチメトリック計算部２０１には、パスメモリ２０８から、生き残りステート（に対応するシンボル候補）が供給される。

すなわち、上述したように、制限ステート数Mが4で、打ち切り長Lが3である場合には、ブランチメトリック計算部２０１には、パスメモリ２０８から、時刻t=n+1の1時刻前の時刻t=nの16個の生き残りステート（に対応するシンボル候補）X''_n,i,j(D205)と、時刻t=n+1の2(=L-1)時刻前の時刻t=n-1の4個の生き残りステート（に対応するシンボル候補）X''_n-1,i(D204)とが供給される。

ブランチメトリック計算部２０１は、時刻t=n+1の1時刻前の時刻t=nの16(=M^L-1)個の生き残りステートX''_n,i,jの１個から分岐する、時刻t=n+1のステート（生き残り候補）X_n+1,i,j,kを、制限ステート数Mである4に制限しながら、時刻t=n+1のステート（生き残り候補）X_n+1,i,j,kへのブランチメトリックBM_i,j,kを計算し、そのブランチメトリックBM_i,j,k(D208)を、ステートメトリック計算部２０２に供給する。

すなわち、ブランチメトリック計算部２０１は、FFT部１２（図３）からの干渉低減信号Y_n(D206)、伝送路特性推定部１７からの伝送路情報H_n+1(D203)、レジスタ２０６からの伝送路情報H_n(D202)、レジスタ２０７からの伝送路情報H_n-1(D201)、並びに、パスメモリ２０８からの、時刻t=nの16個の生き残りステートに対応するシンボル候補（送信シンボルX_nの推定値となり得るシンボル候補）X''_n,i,j(D205)、及び、時刻t=n-1の4個の生き残りステートに対応するシンボル候補（送信シンボルX_n-1の推定値となり得るシンボル候補）X''_n-1,i(D204)を用い、式（３）に従って、１個のシンボル候補X''_n,i,jについて、1個ずつの受信値X'_n+1,i,j、つまり、合計で、16個の受信値X'_n+1,i,jを求める。

さらに、ブランチメトリック計算部２０１は、IQコンスタレーション上で、16個の受信値X'_n+1,i,jそれぞれについて、受信値X'_n+1,i,jにより近い、制限ステート数Mである4個の信号点を選択し、その4個の信号点それぞれが表す4個のシンボルを、送信シンボルX_n+1のシンボル候補X''_n+1,i,j,k=X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3とする。

ブランチメトリック計算部２０１は、16個の受信値X'_n+1,i,jそれぞれについて、4個ずつのシンボル候補X''_n+1,i,j,k=X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3を求めると、つまり、全部で、64(=M^L)個のシンボル候補X''_n+1,i,j,kを求めると、その64個のシンボル候補X''_n+1,i,j,kに対応するステートを、生き残り候補X''_n+1,i,j,k(D207)として、パスメモリ２０８に供給する。

また、ブランチメトリック計算部２０１は、1時刻前の時刻t=nの16個の生き残りステートX''_n,i,jから、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kへのブランチメトリックBMを、式（５）に従って計算し、そのブランチメトリックBM（D208)を、ステートメトリック計算部２０２に供給する。

すなわち、ブランチメトリック計算部２０１は、1個の生き残りステートX''_n,i,jについて、4(=M)個の生き残り候補X''_n+1,i,j,kへのブランチメトリックBMを、式（５）に従って計算し、そのブランチメトリックBM（D208)を、ステートメトリック計算部２０２に供給する。

ステートメトリック計算部２０２は、トレリスにおいて生き残る生き残りステートの候補、つまり、64個の生き残り候補X''_n+1,i,j,kのそれぞれについて、ブランチメトリック計算部２０１からのブランチメトリックBM(D208)を用いて、ステートメトリックsmを計算する。

すなわち、ステートメトリック計算部２０２には、ブランチメトリック計算部２０１から、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kへのブランチメトリックBM(D208)が供給される他、レジスタ２０５に記憶された、1時刻前の時刻t=nの16個の生き残りステートX''_n,i,jのステートメトリックsm(D209)が供給される。

ステートメトリック計算部２０２は、レジスタ２０５からの、1時刻前の時刻t=nの生き残りステートX''_n,i,jのステートメトリックsm(D209)に、ブランチメトリック計算部２０１からの、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリックBM(D208)を加算することで、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリックを求める。

そして、ステートメトリック計算部２０２は、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリックsm(D210)を、最小値検出部２０３、及び、セレクタ２０４に供給する。

最小値検出部２０３は、ステートメトリック計算部２０２からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリックsm(D210)の中で、最小ステートメトリックを検出する。

さらに、最小値検出部２０３は、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kのうちの、最小ステートメトリックの生き残り候補X''_n+1,i,j,kから、打ち切り長L-1だけ遡った時刻t=n-1(=(n+1)-(L-1))の生き残りステートX''_n-1,iを、有効ステートとして検出する。そして、最小値検出部２０３は、有効ステートX''_n-1,iを表す有効ステート情報(D211)を、セレクタ２０４、及び、パスメモリ２０８に供給する。

セレクタ２０４は、ステートメトリック計算部２０２からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリック(D210)の中の、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iに（間接的に）つながる16個の生き残り候補、つまり、新たな生き残りステートとして選択される16個の生き残り候補X''_n+1,i,j,kのステートメトリックを選択する。

さらに、セレクタ２０４は、64個の生き残り候補X''_n+1,i,j,kのステートメトリック(D210)の中から選択された、16個の新たな生き残りステートのステートメトリック(D213)を、レジスタ２０５に供給する。

レジスタ２０５は、セレクタ２０４からの16個の新たな生き残りステートのステートメトリック(D213)を記憶する。レジスタ２０５に記憶された16個の新たな生き残りステートのステートメトリックは、次のキャリアの処理時に、1時刻前の時刻t=nの生き残りステートX''_n,i,jのステートメトリックsm(D209)として、ステートメトリック計算部２０２に供給される。

レジスタ２０６には、伝送路特性推定部１７（図３）から、伝送路情報H_n+1(D203)が供給される。

レジスタ２０６は、伝送路特性推定部１７からの伝送路情報H_n+1(D203)を、1時刻分だけ記憶し（遅延し）、伝送路情報H_n(D202)として、ブランチメトリック計算部２０１、及び、レジスタ２０７に供給する。

レジスタ２０７は、レジスタ２０６からの伝送路情報H_n(D202)を、1時刻分だけ記憶し、伝送路情報H_n-1(D201)として、ブランチメトリック計算部２０１に供給する。

パスメモリ２０８は、セレクタ２０９、レジスタ２１０、セレクタ２１１、レジスタ２１２、及び、セレクタ２１３から構成される。

パスメモリ２０８において、レジスタ２１０及び２１２は、トレリスにおいて、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iから（直接及び間接的に）分岐する生き残りステート、及び、生き残り候補を、新たな生き残り候補として記憶する。

また、パスメモリ２０８において、セレクタ２１３は、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iに対応するシンボル候補を、送信シンボルX_n-1(=X_n+1-(L-1))の推定値X''_n-1,i(=X''_n+1-(L-1),i)(D215)として選択して出力する。

すなわち、パスメモリ２０８のセレクタ２０９には、最小値検出部２０３からの有効ステート情報(D211)と、ブランチメトリック計算部２０１からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,k(D207)とが供給される。

セレクタ２０９は、ブランチメトリック計算部２０１からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,k(D207)の中で、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iから（間接的に）分岐する生き残り候補を、次のキャリアの処理時の時刻t=nの16個の新たな生き残りステートX''_n,i,jとして選択する。さらに、セレクタ２０９は、その16個の新たな生き残りステートX''_n,i,j(D212)を、レジスタ２１０に供給する。

レジスタ２１０は、セレクタ２０９からの16個の新たな生き残りステートX''_n,i,j(D212)を記憶する。さらに、レジスタ２１０は、その記憶値である16個の生き残りステートX''_n,i,j(D205)を、ブランチメトリック計算部２０１、及び、セレクタ２１１に供給する。

セレクタ２１１は、レジスタ２１０からの、時刻t=nの16個の生き残りステートX''_n,i,j(D205)の中で、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iから（直接）分岐する生き残りステートを、次のキャリアの処理時の時刻t=n-1の4個の新たな生き残りステートX''_n-1,iとして選択する。さらに、セレクタ２１１は、その4個の新たな生き残りステートX''_n-1,i(D214)を、レジスタ２１２に供給する。

レジスタ２１２は、セレクタ２１１からの4個の新たな生き残りステートX''_n-1,(D214)を記憶する。さらに、レジスタ２１２は、その記憶値である4個の生き残りステートX''_n-1,(D204)を、ブランチメトリック計算部２０１、及び、セレクタ２１３に供給する。

セレクタ２１３は、レジスタ２１２からの、時刻t=n-1の4個の生き残りステートX''_n-1,i(D204)の中の、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iに対応するシンボル候補を、送信シンボルX_n-1(=X_n+1-(L-1))の推定値X''_n-1,i(=X''_n+1-(L-1),i)(D215)として選択して出力する。

セレクタ２１３が出力する送信シンボルX_n-1(=X_n+1-(L-1))の推定値X''_n-1,i(=X''_n+1-(L-1),i)(D215)は、干渉成分推定部１４（図３）に供給される。

次に、図１６は、図１５の送信系列推定部１３によるシンボル系列の推定の処理を説明するフローチャートである。

ブランチメトリック計算部２０１には、FFT部１２（図３）から、干渉低減信号Y_n(D206)が供給されるとともに、伝送路特性推定部１７（図３）から、伝送路情報H_n+1が供給される。

また、伝送路特性推定部１７からの伝送路情報H_n+1(D203)は、レジスタ２０６にも供給され、レジスタ２０６及び２０７に順次記憶されることで、1時刻前の伝送路情報H_nと、2時刻前の伝送路情報H_n-1にされる。

レジスタ２０６に記憶された伝送路情報H_n(D202)、及び、レジスタ２０７に記憶された伝送路情報H_n-1(D201)は、ブランチメトリック計算部２０１に供給される。

また、このとき、パスメモリ２０８では、レジスタ２１０に記憶された、時刻t=n+1の1時刻前の時刻t=nの16個の生き残りステートX''_n,i,j(D205)と、レジスタ２１２に記憶された、時刻t=n+1の2(=L-1)時刻前の時刻t=n-1の4個の生き残りステートX''_n-1,i(D204)とが、ブランチメトリック計算部２０１に供給される。

ブランチメトリック計算部２０１は、ステップＳ１１において、FFT部１２（図３）からの干渉低減信号Yⁿ(D206)、伝送路特性推定部１７からの伝送路情報H_n+1(D203)、レジスタ２０６からの伝送路情報H_n(D202)、レジスタ２０７からの伝送路情報H_n-1(D201)、並びに、パスメモリ２０８からの、時刻t=nの16個の生き残りステートに対応するシンボル候補X''_n,i,j(D205)、及び、時刻t=n-1の4個の生き残りステートに対応するシンボル候補X''_n-1,i(D204)を用い、式（３）に従って、１個のシンボル候補X''_n,i,jについて、1個ずつの受信値X'_n+1,i,j、つまり、合計で、16個の受信値X'_n+1,i,jを求める。

その後、処理は、ステップＳ１１からステップＳ１２に進み、ブランチメトリック計算部２０１は、IQコンスタレーション上で、16個の受信値X'_n+1,i,jそれぞれについて、受信値X'_n+1,i,jにより近い、制限ステート数Mである4個の信号点を選択し、その4個の信号点それぞれが表す4個のシンボルを、送信シンボルX_n+1のシンボル候補X''_n+1,i,j,k=X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3として検出する。

さらに、ステップＳ１２では、ブランチメトリック計算部２０１は、16個の受信値X'_n+1,i,jそれぞれについての、4個ずつのシンボル候補X''_n+1,i,j,k=X''_n+1,i,j,0，X''_n+1,i,j,1，X''_n+1,i,j,2，X''_n+1,i,j,3、つまり、全部で、64(=M_L)個のシンボル候補X''_n+1,i,j,kに対応する時刻t+1のステートを、生き残り候補X''_n+1,i,j,k(D207)として、パスメモリ２０８のセレクタ２０９に供給する。

その後、処理は、ステップＳ１２からステップＳ１３に進み、ブランチメトリック計算部２０１は、1時刻前の時刻t=nの16個の生き残りステートX''_n,i,jから、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kへのブランチメトリックBMを、式（５）に従って計算する。そして、ブランチメトリック計算部２０１は、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kへのブランチメトリックBM（D208)を、ステートメトリック計算部２０２に供給して、処理は、ステップＳ１３からステップＳ１４に進む。

ステップＳ１４では、ステートメトリック計算部２０２が、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kのそれぞれについて、ブランチメトリック計算部２０１からのブランチメトリックBM(D208)と、レジスタ２０５に記憶された、1時刻前の時刻t=nの16個の生き残りステートX''_n,i,jのステートメトリックsm(D209)とを用いて、ステートメトリックsmを計算する。

すなわち、生き残り候補X''_n+1,i,j,kのブランチメトリックBMを、BM_n+1,i,j,kと表すとともに、生き残りステートX''_n,i,jのステートメトリックsmを、sm_n,i,jと表すこととすると、ステートメトリック計算部２０２は、生き残り候補X''_n+1,i,j,kのステートメトリックsm_n+1,i,j,kを、式sm_n+1,i,j,k=BM_n+1,i,j,k＋sm_n,i,jに従って計算する。

そして、ステートメトリック計算部２０２は、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kのそれぞれのステートメトリックsm_n+1,i,j,k(D210)を、最小値検出部２０３、及び、セレクタ２０４に供給して、処理は、ステップＳ１４からステップＳ１５に進む。

ステップＳ１５では、最小値検出部２０３が、ステートメトリック計算部２０２からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリックsm_n+1,i,j,k(D210)の中で、最小ステートメトリックを検出して、処理は、ステップＳ１６に進む。

ステップＳ１６では、最小値検出部２０３は、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kのうちの、最小ステートメトリックの生き残り候補X''_n+1,i,j,kから、打ち切り長L-1だけ遡った時刻t=n-1(=(n+1)-(L-1))の生き残りステートX''_n-1,iを、有効ステートとして検出する。

さらに、ステップＳ１６では、最小値検出部２０３は、有効ステートX''_n-1,iを表す有効ステート情報(D211)を、セレクタ２０４、並びに、パスメモリ２０８のセレクタ２０９，２１１、及び２１３に供給して、処理は、ステップＳ１７に進む。

ステップＳ１７では、セレクタ２０４が、ステートメトリック計算部２０２からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,kそれぞれのステートメトリックsm_n+1,i,j,k(D210)の中の、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iにつながる16個の生き残り候補（新たな生き残りステートとして選択される16個の生き残り候補）X''_n+1,i,j,kのステートメトリックを選択する。さらに、セレクタ２０４は、64個の生き残り候補X''_n+1,i,j,kのステートメトリック(D210)の中から選択した、16個の新たな生き残りステートのステートメトリック(D213)を、レジスタ２０５に供給する。そして、レジスタ２０５は、セレクタ２０４からの16個の新たな生き残りステートのステートメトリック(D213)を記憶する。

レジスタ２０５に記憶された16個の新たな生き残りステートのステートメトリックは、次のキャリアの処理時に、1時刻前の時刻t=nの生き残りステートX''_n,i,jのステートメトリックsm(D209)として、ステートメトリック計算部２０２に供給される。

さらに、ステップＳ１７では、パスメモリ２０８のセレクタ２１３が、レジスタ２１２に記憶された、時刻t=n-1の4個の生き残りステートX''_n-1,i(D204)の中の、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iに対応するシンボル候補を、送信シンボルX_n-1(=X_n+1-(L-1))の推定値X''_n-1,i(=X''_n+1-(L-1),i)(D215)として選択し、干渉成分推定部１４（図３）に出力する。

また、ステップＳ１７では、パスメモリ２０８のセレクタ２１１が、レジスタ２１０に記憶された、時刻t=nの16個の生き残りステートX''_n,i,j(D205)の中で、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iから（直接）分岐する4個の生き残りステートを、次のキャリアの処理時の時刻t=n-1の4個の新たな生き残りステートX''_n-1,iとして選択する。

そして、セレクタ２１１は、その4個の新たな生き残りステートX''_n-1,i(D214)を、レジスタ２１２に供給し、レジスタ２１２は、セレクタ２１１からの4個の新たな生き残りステートX''_n-1,i(D214)を記憶する。

さらに、ステップＳ１７では、パスメモリ２０８のセレクタ２０９が、ブランチメトリック計算部２０１からの、時刻t=n+1の64個の生き残り候補X''_n+1,i,j,k(D207)の中で、最小値検出部２０３からの有効ステート情報(D211)が表す有効ステートX''_n-1,iから（間接的に）分岐する16個の生き残り候補を、次のキャリアの処理時の時刻t=nの16個の新たな生き残りステートX''_n,i,jとして選択する。

そして、セレクタ２０９は、その16個の新たな生き残りステートX''_n,i,j(D212)を、レジスタ２１０に供給し、レジスタ２１０は、セレクタ２０９からの16個の新たな生き残りステートX''_n,i,j(D212)を記憶する。

その後、ブランチメトリック計算部２０１に対して、FFT部１２（図３）から、次のキャリアについての干渉低減信号Y_n(D206)が供給されるとともに、伝送路特性推定部１７（図３）から、次のキャリアについての伝送路情報H_n+1が供給されると、処理は、ステップＳ１７からステップＳ１１に戻り、以下、同様の処理が繰り返される。

次に、図１７は、本発明を適用した受信装置の他の一実施の形態の構成例を示すブロック図である。

なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。

すなわち、図１７の受信装置は、干渉成分推定部１４、干渉成分除去部１５、及び、誤り訂正部１６が設けられている点で、図３の場合と共通する。

但し、図１７の受信装置は、図３のハニング窓乗算部１１，FFT部１２、送信系列推定部１３、及び、伝送路特性推定部１７にそれぞれ代えて、ハニング窓乗算部１１１，FFT部１１２、送信系列推定部１１３、及び、伝送路特性推定部１１７が設けられている点で、図３の場合と相違している。

図３の受信装置では、上述したように、係数が、-1/2，+1，-1/2のハニング窓、つまり、位相が0のハニング窓（以下、基準ハニング窓ともいう）を、OFDM信号にかけることで、トレリスを用いたシンボル系列の推定のための拘束条件（式（１））を成立させ、シンボル系列の推定が行われる。

これに対して、図１７の受信装置では、位相を、適宜調整したハニング窓、つまり、係数が基準ハニング窓と異なるハニング窓（以下、位相調整ハニング窓ともいう）をも、OFDM信号にかけることで、トレリスを用いたシンボル系列の推定のための拘束条件を成立させる。

そして、図１７の受信装置では、係数が異なる２以上の窓関数としての、基準ハニング窓と位相調整ハニング窓のそれぞれに対して計算されるブランチメトリックを合成することにより、生き残り候補のブランチメトリックを計算し、そのブランチメトリックを用いて、シンボル系列の推定が行われる。

ここで、上述のように、係数が異なる２以上の窓関数を用いて、シンボル系列の推定を行う方法では、１つの窓関数だけを用いる場合に比較して、受信特性をさらに改善することができる。この方法を、窓ダイバーシチともいう。

図１７の受信装置では、窓ダイバーシチによって、シンボル系列の推定が行われる。

すなわち、ハニング窓乗算部１１１には、時間領域の信号である、ベースバンドのOFDM信号が供給される。ハニング窓乗算部１１１は、図３のハニング窓乗算部１１１と同様に、そこに供給されるOFDM信号に、基準ハニング窓をかけ、その結果得られる、基準ハニング窓かけ後のOFDM信号を、FFT部１１２に供給する。

さらに、ハニング窓乗算部１１１は、そこに供給されるOFDM信号に、位相調整ハニング窓をかけ、その結果得られる、位相調整ハニング窓かけ後のOFDM信号を、FFT部１１２に供給する。

FFT部１１２は、図３のFFT部１２と同様に、ハニング窓乗算部１１１からの、基準ハニング窓かけ後のOFDM信号のFFTを行い、その結果得られる、前述の式（１）で表される干渉低減信号Y_nを、送信系列推定部１１３、干渉成分推定部１４、及び、干渉成分除去部１５に供給する。

さらに、FFT部１１２は、ハニング窓乗算部１１１からの、位相調整ハニング窓かけ後のOFDM信号のFFTを行い、その結果得られる干渉低減信号Y^* _nを、送信系列推定部１１３に供給する。

ここで、位相調整ハニング窓かけ後のOFDM信号のFFTを行うことによって得られる干渉低減信号Y^* _nは、式（８）で表される。

・・・（８）

なお、W_m，W_c、及びW_pは、位相調整ハニング窓の係数を表す。係数W_m，W_c、及びW_pは、位相調整ハニング窓の位相によって異なる。

式（８）の干渉低減信号Y^* _nは、式（１）の干渉低減信号Y_nと同様に、畳み込み符号と見ることができ、これにより、送信シンボルX_nは、トレリスを用いた系列の推定を行うことによって推定することができる。

したがって、図１７の受信装置では、OFDM信号に、基準ハニング窓と位相調整ハニング窓のそれぞれをかけることで、トレリスを用いた送信シンボルX_nのシンボル系列（データ系列）の推定のための、式（１）と式（８）の拘束条件を成立させているということができる。

送信系列推定部１１３は、FFT部１１２からの畳み込み符号としての干渉低減信号Y_n及びY^* _nを用いて、送信シンボル（のシンボル系列）を推定し、その送信シンボルの推定値を、干渉成分推定部１４に供給する。

すなわち、送信系列推定部１１３には、FFT部１１２から干渉低減信号Y_n及びY^* _nが供給される他、伝送路特性推定部１１７から、伝送路情報H_n+1、さらには、信頼度が供給されるようになっている。

ここで、信頼度とは、例えば、位相調整ハニング窓の位相が、マルチパスのエコーの位置（位相）（遅延量）に一致している程度（信頼性）を表し、例えば、干渉低減信号Y^* _nに含まれるシンボル間干渉の程度の小ささに対応する。

送信系列推定部１１３は、FFT部１１２からの干渉低減信号Y_n及びY^* _nと、伝送路特性推定部１１７からの伝送路情報H_n+1、及び、信頼度とを用いて、トレリスにより、シンボル系列を推定する。

具体的には、送信系列推定部１１３は、図３の送信系列推定部１３と同様に、式（３）に従い、干渉低減信号Y_n、及び、伝送路情報H_n+1とを用いて、送信シンボルX_n+1の受信値X'_n+1の第１の候補（以下、基準ハニング窓に対する受信値X'_n+1の候補ともいう）を求める。

また、送信系列推定部１１３は、干渉低減信号Y^* _nを表す式（８）を、その雑音成分N_nを無視して、送信シンボルX_n+1について解いた式（９）に従い、干渉低減信号Y^* _n、及び、伝送路情報H_n+1とを用いて、送信シンボルX_n+1の受信値X'_n+1の第２の候補（以下、位相調整ハニング窓に対する受信値X'_n+1の候補ともいう）を求める。

・・・（９）

そして、送信系列推定部１１３は、基準ハニング窓に対する受信値X'_n+1の候補Aと、位相調整ハニング窓に対する受信値X'_n+1の候補Bとを、伝送路特性推定部１１７からの信頼度に応じて合成し、その合成結果を、送信シンボルX_n+1の受信値X'_n+1とする。

すなわち、伝送路特性推定部１１７からの信頼度が、例えば、位相調整ハニング窓の位相がマルチパスのエコーの位相に一致している信頼性を、0ないし1の範囲で表す値αであるとすると、送信系列推定部１１３は、例えば、式(1-α)A+αBに従って、基準ハニング窓に対する受信値X'_n+1の候補Aと、位相調整ハニング窓に対する受信値X'_n+1の候補Bとを合成し、その合成結果を、送信シンボルX_n+1の受信値X'_n+1とする。

その後、送信系列推定部１１３は、図３の送信系列推定部１３と同様に、受信値X'_n+1から、64(=M^L)個の生き残り候補を求め、その生き残り候補へのブランチメトリックの第１の候補と第２の候補とを求める。

すなわち、送信系列推定部１１３は、基準ハニング窓に対して、干渉低減信号Y_nを用い、式（５）に従って、ブランチメトリックを計算し、そのブランチメトリックを、生き残り候補へのブランチメトリックの第１の候補とする。

さらに、送信系列推定部１１３は、位相調整ハニング窓に対して、干渉低減信号Y^* _nを用い、式（５）と同様の式（１０）に従って、ブランチメトリックを計算し、そのブランチメトリックを、生き残り候補へのブランチメトリックの第２の候補とする。

・・・（１０）

そして、送信系列推定部１１３は、基準ハニング窓に対して計算されたブランチメトリックの第１の候補BM₁と、位相調整ハニング窓に対して計算されたブランチメトリックの第２の候補BM₂とを、伝送路特性推定部１１７からの信頼度に応じて合成し、その合成結果を、生き残り候補へのブランチメトリックBMとする。

すなわち、伝送路特性推定部１１７からの信頼度が、例えば、上述したように、位相調整ハニング窓の位相がマルチパスのエコーの位相に一致している信頼性を、0ないし1の範囲で表す値αであるとすると、送信系列推定部１１３は、例えば、式(1-α)BM₁+αBM₂に従って、基準ハニング窓に対して計算されたブランチメトリックの第１の候補BM₁と、位相調整ハニング窓に対して計算されたブランチメトリックの第２の候補BM₂とを合成し、その合成結果を、生き残り候補へのブランチメトリックBMとする。

そして、送信系列推定部１１３は、生き残り候補へのブランチメトリックBMを用いて、その生き残り候補のステートメトリックを求め、以下、図３の送信系列推定部１３と同様にして、送信シンボルX_n+1の推定値X''_n+1を求める。

伝送路特性推定部１１７には、OFDMのパイロット信号が供給されるようになっている。伝送路特性推定部１１７は、図３の伝送路特性推定部１７と同様に、そこに供給されるパイロット信号から、伝送路情報H_n+1を推定し、送信系列推定部１１３に供給する。

さらに、伝送路特性推定部１１７は、そこに供給されるパイロット信号から、マルチパスのエコーの位相を検出する。そして、伝送路特性推定部１１７は、位相調整ハニング窓の位相が、マルチパスのエコーの位相に一致している信頼性を表す信頼度を求め、送信系列推定部１１３に供給する。

図１８は、図１７の送信系列推定部１１３の構成例を示している。

なお、図中、図１５の送信系列推定部１３と同様に構成される部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。

すなわち、図１７の送信系列推定部１１３は、ステートメトリック計算部２０２、最小値検出部２０３、セレクタ２０４、レジスタ２０５ないし２０７、及び、パスメモリ２０８を有する点で、図１５の送信系列推定部１３と共通する。

但し、図１７の送信系列推定部１１３は、図１５のブランチメトリック計算部２０１に代えて、ブランチメトリック計算部３０１が設けられている点で、図１５の送信系列推定部１３と相違する。

ブランチメトリック計算部３０１には、FFT部１１２（図１７）から、干渉低減信号Y_n(D206)及びY^* _nが供給されるとともに、伝送路特性推定部１１７（図１７）から、伝送路情報H_n+1(D203)、及び、信頼度が供給される。

さらに、ブランチメトリック計算部３０１には、レジスタ２０６から、伝送路情報H_n(D202)が供給されるとともに、レジスタ２０７から、伝送路情報H_n-1(D201)が供給される。

また、ブランチメトリック計算部３０１には、パスメモリ２０８から、生き残りステート（に対応するシンボル候補）が供給される。

すなわち、上述したように、制限ステート数Mが4で、打ち切り長Lが3である場合には、ブランチメトリック計算部３０１には、パスメモリ２０８から、時刻t=n+1の1時刻前の時刻t=nの16個の生き残りステート（に対応するシンボル候補）X''_n,i,j(D205)と、時刻t=n+1の2(=L-1)時刻前の時刻t=n-1の4個の生き残りステート（に対応するシンボル候補）X''_n-1,i(D204)とが供給される。

ブランチメトリック計算部３０１は、図１５のブランチメトリック計算部２０１と同様に、時刻t=n+1の1時刻前の時刻t=nの16(=M^L-1)個の生き残りステートX''_n,i,jから分岐する、時刻t=n+1のステート（生き残り候補）X_n+1,i,j,kを、制限ステート数Mである4に制限しながら、時刻t=n+1のステート（生き残り候補）X_n+1,i,j,kへのブランチメトリックBM_i,j,kを計算し、そのブランチメトリックBM_i,j,k(D208)を、ステートメトリック計算部２０２に供給する。

但し、ブランチメトリック計算部３０１は、基準ハニング窓に対して計算されるブランチメトリックの第１の候補BM₁と、位相調整ハニング窓に対して計算されるブランチメトリックの第２の候補BM₂とを、伝送路特性推定部１１７からの信頼度に応じて合成することにより、生き残り候補X_n+1,i,j,kへの（枝の）ブランチメトリックBM_i,j,k(D208)を計算する。

すなわち、ブランチメトリック計算部３０１は、FFT部１１２（図１７）からの干渉低減信号Y_n(D206)、伝送路特性推定部１１７からの伝送路情報H_n+1(D203)、レジスタ２０６からの伝送路情報H_n(D202)、レジスタ２０７からの伝送路情報H_n-1(D201)、並びに、パスメモリ２０８からの、時刻t=nの生き残りステートに対応するシンボル候補（送信シンボルX_nの推定値となり得るシンボル候補）X''_n,i,j(D205)、及び、時刻t=n-1の生き残りステートに対応するシンボル候補（送信シンボルX_n-1の推定値となり得るシンボル候補）X''_n-1,i(D204)を用い、式（３）に従って、基準ハニング窓に対する受信値X'_n+1,i,jの候補を求める。

さらに、ブランチメトリック計算部３０１は、FFT部１１２からの干渉低減信号Y^* _n、伝送路特性推定部１１７からの伝送路情報H_n+1(D203)、レジスタ２０６からの伝送路情報H_n(D202)、レジスタ２０７からの伝送路情報H_n-1(D201)、並びに、パスメモリ２０８からの、時刻t=nの生き残りステートに対応するシンボル候補X''_n,i,j(D205)、及び、時刻t=n-1の生き残りステートに対応するシンボル候補X''_n-1,i(D204)を用い、式（３）に従って、位相調整ハニング窓に対する受信値X'_n+1,i,jの候補を求める。

そして、送信系列推定部１１３は、基準ハニング窓に対する受信値X'_n+1の候補と、位相調整ハニング窓に対する受信値X'_n+1の候補とを、伝送路特性推定部１１７からの信頼度に応じて合成し、その合成結果を、送信シンボルX_n+1の受信値X'_n+1,i,jとする。

その後、ブランチメトリック計算部３０１は、図１５のブランチメトリック計算部２０１と同様に、IQコンスタレーション上で、1個の受信値X'_n+1,i,jごとに、その受信値X'_n+1,i,jにより近い、制限ステート数Mである4個の信号点を選択し、その信号点が表すシンボルを、送信シンボルX_n+1のシンボル候補X''_n+1,i,j,kとする。

ブランチメトリック計算部３０１は、図１５のブランチメトリック計算部２０１と同様にして、64(=M^L)個のシンボル候補X''_n+1,i,j,kを求めると、その64個のシンボル候補X''_n+1,i,j,kに対応するステートを、生き残り候補X''_n+1,i,j,k(D207)として、パスメモリ２０８に供給する。

また、ブランチメトリック計算部３０１は、図１５のブランチメトリック計算部２０１と同様に、式（５）を計算し、その計算結果を、1時刻前の時刻t=nの生き残りステートX''_n,i,jから、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリックBMの第１の候補とする。

すなわち、ブランチメトリック計算部３０１は、基準ハニング窓に対して、干渉低減信号Y_nを用い、式（５）に従って、ブランチメトリックを計算し、そのブランチメトリックを、生き残り候補X''_n+1,i,j,kへのブランチメトリック（1時刻前の時刻t=nの生き残りステートX''_n,i,jから、時刻t=n+1の生き残り候補X''_n+1,i,j,kへのブランチメトリック）BMの第１の候補とする。

また、ブランチメトリック計算部３０１は、位相調整ハニング窓に対して、干渉低減信号Y^* _nを用い、式（１０）に従って、ブランチメトリックを計算し、そのブランチメトリックを、生き残り候補X''_n+1,i,j,kへのブランチメトリックBMの第２の候補とする。

その後、ブランチメトリック計算部３０１は、基準ハニング窓に対して計算されたブランチメトリックBMの第１の候補と、位相調整ハニング窓に対して計算されたブランチメトリックBMの第２の候補とを、伝送路特性推定部１１７からの信頼度に応じて合成する。

そして、ブランチメトリック計算部３０１は、基準ハニング窓に対して計算されたブランチメトリックBMの第１の候補と、位相調整ハニング窓に対して計算されたブランチメトリックBMの第２の候補との合成結果を、生き残り候補へのブランチメトリックBM(D208)として、ステートメトリック計算部２０２に供給する。

以上のように、図３及び図１７の受信装置では、干渉低減信号Y_n(Y^* _n)から、送信シンボルX_n+1の受信値X'_n+1が求められ、送信シンボルX_n+1のシンボル候補が、受信値X'_n+1からユークリッド距離がより近い、制限ステート数M（少数）の信号点が表すシンボルに絞り込まれる（制限される）。

さらに、シンボル候補に対応するステートを、生き残りステートの候補である生き残り候補として、生き残り候補へのブランチメトリック（生き残り候補に遷移する枝のブランチメトリック）が計算され、そのブランチメトリックを用いて、生き残り候補のステートメトリックが計算される。

その後、生き残り候補のステートメトリックの最小ステートメトリックが検出され、その最小ステートメトリックの生き残り候補から、打ち切り長L-1だけ遡った生き残りステートである有効ステートに対応するシンボル候補が、送信シンボルX_n+1-(L-1)の推定値として選択されて出力される。

さらに、有効ステートから分岐する生き残りステート、及び、生き残り候補が、新たな生き残りステートとして選択され、他のステートは、メトリック（ステートメトリック、ブランチメトリック）の計算対象から除外される。

したがって、生き残り候補の数を、一定の数であるM^L個に保ちながら、シンボル系列の推定を行うことができ、生き残り候補の数が膨大になって、装置規模が爆発的に増加することを抑制することができる。

また、窓ダイバーシチを利用する場合には、受信特性をより改善することができる。

なお、制限ステート数M、及び、打ち切り長Lを大とする場合には、受信特性の改善を図ることができるが、生き残りステート、及び、生き残り候補の数が増加し、装置規模が増加する。

一方、制限ステート数M、及び、打ち切り長Lを小とする場合には、生き残りステート、及び、生き残り候補の数が減少して、装置規模が抑制されるが、受信特性が多少劣化する。

したがって、制限ステート数M、及び、打ち切り長Lは、受信特性と装置規模とを比較考量して決めることが望ましい。

また、打ち切り長Lは、2以上の値とすることが望ましい。打ち切り長Lが1の場合には、系列の「推定」（トレリスを用いた推定）ではなくなってしまうからである。

次に、上述した一連の処理は、ハードウェアにより行うこともできるし、ソフトウェアにより行うこともできる。一連の処理をソフトウェアによって行う場合には、そのソフトウェアを構成するプログラムが、汎用のコンピュータ等にインストールされる。

そこで、図１９は、上述した一連の処理を実行するプログラムがインストールされるコンピュータの一実施の形態の構成例を示している。

プログラムは、コンピュータに内蔵されている記録媒体としてのハードディスク５０５やROM５０３に予め記録しておくことができる。

あるいはまた、プログラムは、フレキシブルディスク、CD-ROM(Compact Disc Read Only Memory)，MO(Magneto Optical)ディスク，DVD(Digital Versatile Disc)、磁気ディスク、半導体メモリなどのリムーバブル記録媒体５１１に、一時的あるいは永続的に格納（記録）しておくことができる。このようなリムーバブル記録媒体５１１は、いわゆるパッケージソフトウエアとして提供することができる。

なお、プログラムは、上述したようなリムーバブル記録媒体５１１からコンピュータにインストールする他、ダウンロードサイトから、ディジタル衛星放送用の人工衛星を介して、コンピュータに無線で転送したり、LAN(Local Area Network)、インターネットといったネットワークを介して、コンピュータに有線で転送し、コンピュータでは、そのようにして転送されてくるプログラムを、通信部５０８で受信し、内蔵するハードディスク５０５にインストールすることができる。

コンピュータは、CPU(Central Processing Unit)５０２を内蔵している。CPU５０２には、バス５０１を介して、入出力インタフェース５１０が接続されており、CPU５０２は、入出力インタフェース５１０を介して、ユーザによって、キーボードや、マウス、マイク等で構成される入力部５０７が操作等されることにより指令が入力されると、それに従って、ROM(Read Only Memory)５０３に格納されているプログラムを実行する。あるいは、また、CPU５０２は、ハードディスク５０５に格納されているプログラム、衛星若しくはネットワークから転送され、通信部５０８で受信されてハードディスク５０５にインストールされたプログラム、またはドライブ５０９に装着されたリムーバブル記録媒体５１１から読み出されてハードディスク５０５にインストールされたプログラムを、RAM(Random Access Memory)５０４にロードして実行する。これにより、CPU５０２は、上述したフローチャートにしたがった処理、あるいは上述したブロック図の構成により行われる処理を行う。そして、CPU５０２は、その処理結果を、必要に応じて、例えば、入出力インタフェース５１０を介して、LCD(Liquid Crystal Display)やスピーカ等で構成される出力部５０６から出力、あるいは、通信部５０８から送信、さらには、ハードディスク５０５に記録等させる。

ここで、本明細書において、コンピュータに各種の処理を行わせるためのプログラムを記述する処理ステップは、必ずしもフローチャートとして記載された順序に沿って時系列に処理する必要はなく、並列的あるいは個別に実行される処理（例えば、並列処理あるいはオブジェクトによる処理）も含むものである。

また、プログラムは、１のコンピュータにより処理されるものであっても良いし、複数のコンピュータによって分散処理されるものであっても良い。さらに、プログラムは、遠方のコンピュータに転送されて実行されるものであっても良い。

なお、本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々の変更が可能である。

ビタビ復号法のトレリスを示す図である。 Mアルゴリズムのトレリスを示す図である。 本発明を適用した受信装置の一実施の形態の構成例を示すブロック図である。 送信系列推定部１３による、送信シンボルの推定の処理を説明するためのトレリスを示す図である。 送信系列推定部１３による、送信シンボルの推定の処理を説明するためのトレリスを示す図である。 送信系列推定部１３による、送信シンボルの推定の処理を説明するためのトレリスを示す図である。 生き残りステートから分岐するステート（枝）の数を、制限ステート数Mに制限する方法を説明する図である。 制限ステート数Mが2であり、かつ、最新生き残りステート数Cが4である場合の、Xアルゴリズムのトレリスを示す図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 Xアルゴリズムによる、ISDB-Tの規格に準拠するOFDM信号からのシンボル系列の推定の処理を説明する図である。 送信系列推定部１３の構成例を示すブロック図である。 送信系列推定部１３の処理を説明するフローチャートである。 本発明を適用した受信装置の他の一実施の形態の構成例を示すブロック図である。 送信系列推定部１１３の構成例を示すブロック図である。 本発明を適用したコンピュータの一実施の形態の構成例を示すブロック図である。

符号の説明

１１ ハニング窓乗算部， １２ FFT部， １３ 送信系列推定部， １４ 干渉成分推定部， １５ 干渉成分除去部， １６ 誤り訂正部， １７ 伝送路特性推定部， １１１ ハニング窓乗算部， １１２ FFT部， １１３ 送信系列推定部， １１７ 伝送路特性推定部， ２０１ ブランチメトリック計算部， ２０２ ステートメトリック計算部， ２０３ 最小値検出部， ２０４ セレクタ， ２０５ないし２０７ レジスタ， ２０８ パスメモリ， ２０９ セレクタ， ２１０ レジスタ， ２１１ セレクタ， ２１２ レジスタ， ２１３ セレクタ， ３０１ ブランチメトリック計算部， ５０１ バス， ５０２ CPU， ５０３ ROM， ５０４ RAM， ５０５ ハードディスク， ５０６ 出力部， ５０７ 入力部， ５０８ 通信部， ５０９ ドライブ， ５１０ 入出力インタフェース， ５１１ リムーバブル記録媒体

Claims (11)

1. データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックを計算するブランチメトリック計算手段と、
   前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックを、前記ブランチメトリックを用いて計算するステートメトリック計算手段と、
   前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出する検出手段と、
   前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートを、前記生き残りステートとして記憶する記憶手段と、
   前記有効ステートに対応する、前記データ系列を構成するデータの候補を、前記データの推定値として選択する選択手段と
   を備えるデータ処理装置。
2. 前記生き残りステートの候補の数が、M^L個に制限される
   請求項１に記載のデータ処理装置。
3. 前記ブランチメトリック計算手段は、前記トレリスの前記ステートから分岐する前記ステートを、前記ブランチメトリックが小さい前記M個の枝につながる前記ステートに制限する
   請求項１に記載のデータ処理装置。
4. 前記ブランチメトリック計算手段は、前記トレリスの前記ステートから分岐する前記ステートを、前記トレリスを用いた前記データ系列の推定のための拘束条件に従って計算されるIQコンスタレーション上の点に近い前記M個の信号点が表す前記データの候補に対応するステートに制限する
   請求項１に記載のデータ処理装置。
5. 信号に、窓関数をかけることにより、前記トレリスを用いた前記データ系列の推定のための拘束条件が成立させられる
   請求項１に記載のデータ処理装置。
6. 前記窓関数は、ハニング窓である
   請求項５に記載のデータ処理装置。
7. 前記窓関数は、時間領域、又は周波数領域でかけられる
   請求項５に記載のデータ処理装置。
8. 信号に、窓関数をかけることにより、前記トレリスを用いた前記データ系列の推定のための拘束条件が成立させられ、
   前記ブランチメトリック計算手段は、係数が異なる２以上の前記窓関数それぞれに対して計算される前記ブランチメトリックを合成することにより、前記トレリスの前記ステートから分岐する前記ステートにつながる枝の前記ブランチメトリックを計算する
   請求項１に記載のデータ処理装置。
9. データ系列の推定を行うデータ処理装置が、
   前記データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックを計算し、
   前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックを、前記ブランチメトリックを用いて計算し、
   前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出し、
   前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートを、前記生き残りステートとして記憶し、
   前記有効ステートに対応する、前記データ系列を構成するデータの候補を、前記データの推定値として選択する
   ステップを含むデータ処理方法。
10. データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限しながら、ブランチメトリックを計算するブランチメトリック計算手段と、
    前記トレリスにおいて生き残る生き残りステートの候補のステートメトリックを、前記ブランチメトリックを用いて計算するステートメトリック計算手段と、
    前記生き残りステートの候補のステートメトリックの中で、最小のステートメトリックである最小ステートメトリックを検出する検出手段と、
    前記トレリスにおいて、前記最小ステートメトリックの前記ステートから、前記データ系列の推定を打ち切る単位である打ち切り長L-1だけ遡った前記ステートである有効ステートから分岐するステートを、前記生き残りステートとして記憶する記憶手段と、
    前記有効ステートに対応する、前記データ系列を構成するデータの候補を、前記データの推定値として選択する選択手段と
    して、コンピュータを機能させるためのプログラム。
11. データ系列の推定に用いるトレリスのステートから分岐する前記ステートの数を、所定の数Mに制限し、前記データ系列の推定を打ち切る単位である打ち切り長L分の前記ステートの系列を用いて、前記データ系列を構成するデータの推定値を出力するときに、前記トレリスにおいて生き残る生き残りステートの候補の数を、M^L個に制限する
    データ処理装置。

Images