JP2551296B2 - 系列推定装置 - Google Patents
系列推定装置Info
- Publication number
- JP2551296B2 JP2551296B2 JP4109961A JP10996192A JP2551296B2 JP 2551296 B2 JP2551296 B2 JP 2551296B2 JP 4109961 A JP4109961 A JP 4109961A JP 10996192 A JP10996192 A JP 10996192A JP 2551296 B2 JP2551296 B2 JP 2551296B2
- Authority
- JP
- Japan
- Prior art keywords
- response
- transmission line
- state
- transmission
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Description
【０００１】
【産業上の利用分野】本発明は、伝送路特性の時間的な
変動に追随して送信信号系列の推定を行う系列推定装置
に関する。
変動に追随して送信信号系列の推定を行う系列推定装置
に関する。
【０００２】
【従来の技術】最尤系列推定装置（ＭＬＳＥ）は等化能
力の最も優れた等化方式として知られている（例えば、
文献１：Ｇ．Ｄ．Ｆｏｒｎｅｙ，”Ｍａｘｉｍｕｍ Ｌ
ｉｋｅｌｉｈｏｏｄ Ｓｅｑｕｅｎｃｅ Ｅｓｔｉｍａ
ｔｉｏｎ ｏｆ Ｄｉｇｉｔａｌ Ｓｅｑｕｅｎｃｅｓ
ｉｎ ｔｈｅ ｐｒｅｓｅｎｃｅ ｏｆ ｉｎｔｅｒ
ｓｙｍｂｏｌ ｉｎｔｅｒｆｅｒｅｎｃｅ，”ＩＥＥＥ
Ｔｒａｎｓａｃｔｉｏｎ ｏｎ Ｉｎｆｏｒｍａｔｉ
ｏｎ Ｔｈｅｏｒｙ，ｖｏｌ．ＩＴ−１８，ｎｏ．３，
Ｍａｙ １９７２）。最尤系列推定装置は一般に単一の
伝送路応答推定器を備えており、伝送路推定既知の系列
を受信する際にこの伝送路応答推定器を用いて行う。
力の最も優れた等化方式として知られている（例えば、
文献１：Ｇ．Ｄ．Ｆｏｒｎｅｙ，”Ｍａｘｉｍｕｍ Ｌ
ｉｋｅｌｉｈｏｏｄ Ｓｅｑｕｅｎｃｅ Ｅｓｔｉｍａ
ｔｉｏｎ ｏｆ Ｄｉｇｉｔａｌ Ｓｅｑｕｅｎｃｅｓ
ｉｎ ｔｈｅ ｐｒｅｓｅｎｃｅ ｏｆ ｉｎｔｅｒ
ｓｙｍｂｏｌ ｉｎｔｅｒｆｅｒｅｎｃｅ，”ＩＥＥＥ
Ｔｒａｎｓａｃｔｉｏｎ ｏｎ Ｉｎｆｏｒｍａｔｉ
ｏｎ Ｔｈｅｏｒｙ，ｖｏｌ．ＩＴ−１８，ｎｏ．３，
Ｍａｙ １９７２）。最尤系列推定装置は一般に単一の
伝送路応答推定器を備えており、伝送路推定既知の系列
を受信する際にこの伝送路応答推定器を用いて行う。
【０００３】また、伝送路の特性が時間的に変動する場
合には、この伝送路の特性の時間的な変動に追従させる
ような適応最尤系列推定装置も提案されている（例え
ば、文献２：Ｇ．Ｕｎｇｅｒｂｏｅｃｋ，”Ａｄａｐｔ
ｉｖｅ Ｍａｘｉｍｕｍ Ｌｉｋｅｌｅｈｏｏｄ Ｒｅ
ｃｅｉｖｅｒ ｆｏｒ Ｃａｒｒｉｅｒ−Ｍｏｄｕｌａ
ｔｅｄ Ｄａｔａ Ｔｒａｎｓｍｉｓｓｉｏｎ Ｓｙｓ
ｔｅｍｓ，”ＩＥＥＥＴｒａｎｓａｃｔｉｏｎ ｏｎ
Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ，ｖｏｌ．ＣＯＭ−２
２，ｎｏ．５，Ｍａｙ １９７４）。適応最尤系列推定
装置は、まず既知系列を受信する際に伝送路応答を求
め、それ以後情報データ系列を受信するときは伝送路推
定器を適応アルゴリズムを用いて動作させ、伝送路応答
を逐次更新していくことで伝送路特性の時間的変動に追
従していくことを特徴としている。しかし、適応最尤系
列推定装置では、伝送路特性が高速変動にする場合その
伝送路推定器の適応動作が追いつかなくなる。
合には、この伝送路の特性の時間的な変動に追従させる
ような適応最尤系列推定装置も提案されている（例え
ば、文献２：Ｇ．Ｕｎｇｅｒｂｏｅｃｋ，”Ａｄａｐｔ
ｉｖｅ Ｍａｘｉｍｕｍ Ｌｉｋｅｌｅｈｏｏｄ Ｒｅ
ｃｅｉｖｅｒ ｆｏｒ Ｃａｒｒｉｅｒ−Ｍｏｄｕｌａ
ｔｅｄ Ｄａｔａ Ｔｒａｎｓｍｉｓｓｉｏｎ Ｓｙｓ
ｔｅｍｓ，”ＩＥＥＥＴｒａｎｓａｃｔｉｏｎ ｏｎ
Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ，ｖｏｌ．ＣＯＭ−２
２，ｎｏ．５，Ｍａｙ １９７４）。適応最尤系列推定
装置は、まず既知系列を受信する際に伝送路応答を求
め、それ以後情報データ系列を受信するときは伝送路推
定器を適応アルゴリズムを用いて動作させ、伝送路応答
を逐次更新していくことで伝送路特性の時間的変動に追
従していくことを特徴としている。しかし、適応最尤系
列推定装置では、伝送路特性が高速変動にする場合その
伝送路推定器の適応動作が追いつかなくなる。
【０００４】これに対して、特願平２−２０３４３６で
は高速に変動する伝送路に対しても追従することが可能
な新しい形の系列推定装置を提案している。この装置
は、送信信号系列のみならず伝送路の特性も未知である
として、系列毎に対応する伝送路応答を推定してビタビ
アルゴリズムを適用することを特徴とする。伝送路応答
の推定は、送信信号系列候補、伝送路応答、受信信号の
三者で定まる伝送路方程式を系列毎に解くことによって
行っている。これは伝送路応答の最適解を逐次求めるこ
とに相当するので、特願平２−２０３４３６の系列推定
装置は高速な伝送路変動に対しても追従できる。以下、
この装置を基本形ブラインドビタビ等化器、あるいは単
に基本形と呼ぶことにする。
は高速に変動する伝送路に対しても追従することが可能
な新しい形の系列推定装置を提案している。この装置
は、送信信号系列のみならず伝送路の特性も未知である
として、系列毎に対応する伝送路応答を推定してビタビ
アルゴリズムを適用することを特徴とする。伝送路応答
の推定は、送信信号系列候補、伝送路応答、受信信号の
三者で定まる伝送路方程式を系列毎に解くことによって
行っている。これは伝送路応答の最適解を逐次求めるこ
とに相当するので、特願平２−２０３４３６の系列推定
装置は高速な伝送路変動に対しても追従できる。以下、
この装置を基本形ブラインドビタビ等化器、あるいは単
に基本形と呼ぶことにする。
【０００５】
【発明が解決しようとする課題】しかしながら、基本形
ブラインドビタビ等化器では、ビタビアルゴリズムを適
用する状態遷移図（トレリス線図）の状態は伝送路応答
計算に必要なシンボル数をその長さとする全ての送信信
号系列候補により与えられるので、伝送路応答計算に必
要なシンボル数が多い場合に状態数が大きくなり、系列
推定装置を実現する回路が複雑になる、あるいはデジタ
ルシグナルプロセッサ等を用いて実現する場合の処理量
がきわめて多くなるという欠点がある。
ブラインドビタビ等化器では、ビタビアルゴリズムを適
用する状態遷移図（トレリス線図）の状態は伝送路応答
計算に必要なシンボル数をその長さとする全ての送信信
号系列候補により与えられるので、伝送路応答計算に必
要なシンボル数が多い場合に状態数が大きくなり、系列
推定装置を実現する回路が複雑になる、あるいはデジタ
ルシグナルプロセッサ等を用いて実現する場合の処理量
がきわめて多くなるという欠点がある。
【０００６】そこで、本発明の目的は、より少ない状態
数で、すなわち、より簡単な回路あるいはより少ない処
理量で、高速に変動する伝送路に追従することが可能な
系列推定装置を提供することにある。
数で、すなわち、より簡単な回路あるいはより少ない処
理量で、高速に変動する伝送路に追従することが可能な
系列推定装置を提供することにある。
【０００７】
【課題を解決するための手段】第１の発明の系列推定装
置は、受信信号のサンプル値を複数個記憶するレジスタ
と、該レジスタから複数個の前記サンプル値を入力し、
送信される可能性のあるＭ（２以上の整数）値信号系列
から指定された複数の信号系列に対してのみ現時刻の伝
送路応答をそれぞれ推定する伝送路応答計算回路と、該
伝送路応答計算回路で求めた前記複数の信号系列に対す
る前記現時刻の伝送路応答の値の正当性を検査し、正当
である場合は前記伝送路応答計算回路の出力を、正当で
ない場合は同じ信号系列に対する前時刻の伝送路応答推
定値を、現時刻の伝送路応答推定値としてそれぞれ出力
する伝送路応答検査回路と、該伝送路応答検査回路が出
力する前記現時刻の伝送路応答推定値を記憶するととも
に、前時刻に記憶した伝送路応答推定値を逆に前記伝送
路応答検査回路に供給する伝送路応答記憶回路と、前記
伝送路応答検査回路が出力した前記現時刻の伝送路応答
推定値と各縮退状態に対応する非縮退状態と現時刻の送
信信号候補との３者に基づいて前記複数の系列のそれぞ
れの系列に対する仮想受信信号点を求め、前記受信信号
のサンプル値との距離を求めるブランチメトリック計算
回路と、各時刻での前記各縮退状態に対応する非縮退状
態を出力しながら前記ブランチメトリック計算回路の出
力を受けてビタビアルゴリズムにより受信信号を判定す
るビタビプロセッサと、前記非縮退状態を記憶し、前記
伝送路応答計算回路に対して伝送路応答を計算する前記
複数の信号系列を指定する非縮退状態記憶回路と、から
構成され、縮退された状態遷移図に基づいて送信信号系
列の推定を行うことを特徴とする。
置は、受信信号のサンプル値を複数個記憶するレジスタ
と、該レジスタから複数個の前記サンプル値を入力し、
送信される可能性のあるＭ（２以上の整数）値信号系列
から指定された複数の信号系列に対してのみ現時刻の伝
送路応答をそれぞれ推定する伝送路応答計算回路と、該
伝送路応答計算回路で求めた前記複数の信号系列に対す
る前記現時刻の伝送路応答の値の正当性を検査し、正当
である場合は前記伝送路応答計算回路の出力を、正当で
ない場合は同じ信号系列に対する前時刻の伝送路応答推
定値を、現時刻の伝送路応答推定値としてそれぞれ出力
する伝送路応答検査回路と、該伝送路応答検査回路が出
力する前記現時刻の伝送路応答推定値を記憶するととも
に、前時刻に記憶した伝送路応答推定値を逆に前記伝送
路応答検査回路に供給する伝送路応答記憶回路と、前記
伝送路応答検査回路が出力した前記現時刻の伝送路応答
推定値と各縮退状態に対応する非縮退状態と現時刻の送
信信号候補との３者に基づいて前記複数の系列のそれぞ
れの系列に対する仮想受信信号点を求め、前記受信信号
のサンプル値との距離を求めるブランチメトリック計算
回路と、各時刻での前記各縮退状態に対応する非縮退状
態を出力しながら前記ブランチメトリック計算回路の出
力を受けてビタビアルゴリズムにより受信信号を判定す
るビタビプロセッサと、前記非縮退状態を記憶し、前記
伝送路応答計算回路に対して伝送路応答を計算する前記
複数の信号系列を指定する非縮退状態記憶回路と、から
構成され、縮退された状態遷移図に基づいて送信信号系
列の推定を行うことを特徴とする。
【０００８】第２の発明の系列推定装置は、受信信号の
サンプル値を複数個記憶するレジスタと、該レジスタか
ら複数個の前記サンプル値を入力し、送信される可能性
のあるＭ（２以上の整数）値信号系列から指定された複
数の信号系列に対してのみ現時刻の伝送路応答をそれぞ
れ推定する伝送路応答計算回路と、該伝送路応答計算回
路で求めた複数の系列に対する前記現時刻の伝送路応答
の値の正当性を検査し、正当である場合は前記伝送路応
答計算回路の出力を、正当でない場合は同じ信号系列に
対する前時刻の伝送路応答推定値を状態遷移図と非縮退
状態の記憶とが定める規則に従って変換した値を、現時
刻の伝送路応答推定値としてそれぞれ出力する伝送路応
答変換回路と、該伝送路応答変換回路が出力する前記現
時刻の伝送路応答推定値を記憶するとともに前時刻に記
憶した伝送路応答推定値を逆に前記伝送路応答変換回路
に供給する伝送路応答記憶回路と、前記伝送路応答変換
回路が出力した前記現時刻の伝送路応答推定値と各縮退
状態に対応する非縮退状態と現時刻の送信信号候補との
３者に基づいて前記複数の系列のそれぞれの系列に対す
る仮想受信信号点を求め、前記受信信号のサンプル値と
の距離を求めるブランチメトリック計算回路と、各時刻
での前記各縮退状態に対応する非縮退状態を出力しなが
ら前記ブランチメトリック計算回路の出力を受けてビタ
ビアルゴリズムにより最尤系列を判定し、最尤系列の最
古状態に対応する非縮退状態の記憶に従い判定値を出力
するビタビプロセッサと、前記非縮退状態を記憶する非
縮退状態記憶回路と、該非縮退状態記憶回路が与える前
記非縮退状態を変換し、変換された状態を基に前記伝送
路応答計算回路に対して伝送路応答を計算する前記複数
の信号系列を指定する状態変換回路と、から構成され、
縮退された状態遷移図に基づいて送信信号系列の推定を
行うことを特徴とする。
サンプル値を複数個記憶するレジスタと、該レジスタか
ら複数個の前記サンプル値を入力し、送信される可能性
のあるＭ（２以上の整数）値信号系列から指定された複
数の信号系列に対してのみ現時刻の伝送路応答をそれぞ
れ推定する伝送路応答計算回路と、該伝送路応答計算回
路で求めた複数の系列に対する前記現時刻の伝送路応答
の値の正当性を検査し、正当である場合は前記伝送路応
答計算回路の出力を、正当でない場合は同じ信号系列に
対する前時刻の伝送路応答推定値を状態遷移図と非縮退
状態の記憶とが定める規則に従って変換した値を、現時
刻の伝送路応答推定値としてそれぞれ出力する伝送路応
答変換回路と、該伝送路応答変換回路が出力する前記現
時刻の伝送路応答推定値を記憶するとともに前時刻に記
憶した伝送路応答推定値を逆に前記伝送路応答変換回路
に供給する伝送路応答記憶回路と、前記伝送路応答変換
回路が出力した前記現時刻の伝送路応答推定値と各縮退
状態に対応する非縮退状態と現時刻の送信信号候補との
３者に基づいて前記複数の系列のそれぞれの系列に対す
る仮想受信信号点を求め、前記受信信号のサンプル値と
の距離を求めるブランチメトリック計算回路と、各時刻
での前記各縮退状態に対応する非縮退状態を出力しなが
ら前記ブランチメトリック計算回路の出力を受けてビタ
ビアルゴリズムにより最尤系列を判定し、最尤系列の最
古状態に対応する非縮退状態の記憶に従い判定値を出力
するビタビプロセッサと、前記非縮退状態を記憶する非
縮退状態記憶回路と、該非縮退状態記憶回路が与える前
記非縮退状態を変換し、変換された状態を基に前記伝送
路応答計算回路に対して伝送路応答を計算する前記複数
の信号系列を指定する状態変換回路と、から構成され、
縮退された状態遷移図に基づいて送信信号系列の推定を
行うことを特徴とする。
【０００９】
【作用】基本形ブラインドビタビ等化器では、ビタビア
ルゴリズムを適用する状態遷移図（図４に示されるトレ
リス線図）の各状態（これを基本形状態、あるいは非縮
退状態と呼ぶ）は、伝送路応答計算に必要なシンボル数
をその長さとする送信信号系列候補により与えられる。
一方、本発明の系列推定装置では、トレリス線図の状態
（これを縮退状態と呼ぶ）をシンボル数が伝送路応答計
算に必要な数より少ない送信信号系列候補で与える。そ
して、各縮退状態に対しては毎時刻それぞれの生残りパ
スの履歴から過去のシンボルを補い、伝送路応答計算に
必要な送信信号系列候補、すなわち該状態に対応する基
本形状態を再現し、それを用いて伝送路推定を行うよう
にする。そのため、本発明では、状態数を基本形に比べ
て大幅に削減しながら、高速に変動する伝送路に追従す
ることが可能な系列推定装置を実現できる。
ルゴリズムを適用する状態遷移図（図４に示されるトレ
リス線図）の各状態（これを基本形状態、あるいは非縮
退状態と呼ぶ）は、伝送路応答計算に必要なシンボル数
をその長さとする送信信号系列候補により与えられる。
一方、本発明の系列推定装置では、トレリス線図の状態
（これを縮退状態と呼ぶ）をシンボル数が伝送路応答計
算に必要な数より少ない送信信号系列候補で与える。そ
して、各縮退状態に対しては毎時刻それぞれの生残りパ
スの履歴から過去のシンボルを補い、伝送路応答計算に
必要な送信信号系列候補、すなわち該状態に対応する基
本形状態を再現し、それを用いて伝送路推定を行うよう
にする。そのため、本発明では、状態数を基本形に比べ
て大幅に削減しながら、高速に変動する伝送路に追従す
ることが可能な系列推定装置を実現できる。
【００１０】
【実施例】以下、図面を参照して、まず、特願平２−２
０３４３６に示される基本形ブラインドビタビ等化器を
説明し、次に、それと比較する形で本発明を説明する。
０３４３６に示される基本形ブラインドビタビ等化器を
説明し、次に、それと比較する形で本発明を説明する。
【００１１】基本形ブランイドビタビ等化器は、図３で
示される。基本形は、時刻ｔでの伝送路応答を送信され
る可能性のある複数の信号系列候補に対してそれぞれ最
小二乗法により求め、それらを時刻（ｔ＋１）における
ブランチメトリック計算において用いる。
示される。基本形は、時刻ｔでの伝送路応答を送信され
る可能性のある複数の信号系列候補に対してそれぞれ最
小二乗法により求め、それらを時刻（ｔ＋１）における
ブランチメトリック計算において用いる。
【００１２】以下では、時刻ｔでの伝送路インパルス応
答は（Ｌ＋１）シンボル分あるとし、これをベクトルｈ
t T ＝［ｈt 0 、ｈt 1 ，…，ｈt L ］で表現する。受
信信号｛ｒt ｝は各時刻でレジスタ３０２に蓄えられて
いく。いま、時刻（ｔ＋１）の受信信号がレジスタ３０
２に入力されたときに注目する。このとき、レジスタ３
０２の遅延素子群には送信信号がＮ時刻分、すなわちベ
クトル［ｒt …ｒt -N + 1 ］の形で記憶されている。
これを時刻ｔでの受信信号ベクトルｒt とする。
答は（Ｌ＋１）シンボル分あるとし、これをベクトルｈ
t T ＝［ｈt 0 、ｈt 1 ，…，ｈt L ］で表現する。受
信信号｛ｒt ｝は各時刻でレジスタ３０２に蓄えられて
いく。いま、時刻（ｔ＋１）の受信信号がレジスタ３０
２に入力されたときに注目する。このとき、レジスタ３
０２の遅延素子群には送信信号がＮ時刻分、すなわちベ
クトル［ｒt …ｒt -N + 1 ］の形で記憶されている。
これを時刻ｔでの受信信号ベクトルｒt とする。
【００１３】 ｒt T ＝［ｒt ，ｒt - 1 ，…，ｒt - N + 1 ］ （１） 時刻（ｔ＋１）では、ベクトルｒt の各成分が伝送路計
算回路３０３に向けて出力される。
算回路３０３に向けて出力される。
【００１４】伝送路計算回路３０３は、ベクトルｒt を
用いて送信される可能性のある複数の信号系列に対して
時刻ｔでの伝送路応答をそれぞれ最小二乗法により求め
る。いま、時刻ｔまでに送信された信号をＬ時刻分まと
めた系列をベトクルＳt T ＝［ｓt ，ｓt - 1 ，…ｓ
t - L ］、送信信号とは独立な観測過程を含めた上での
加法性伝送路雑音をｖt とすると、時刻ｔでの受信信号
ｒt は、式（１）で示されるようにベクトルｈt とベク
トルｓt との畳込みと雑音の和で与えられる。
用いて送信される可能性のある複数の信号系列に対して
時刻ｔでの伝送路応答をそれぞれ最小二乗法により求め
る。いま、時刻ｔまでに送信された信号をＬ時刻分まと
めた系列をベトクルＳt T ＝［ｓt ，ｓt - 1 ，…ｓ
t - L ］、送信信号とは独立な観測過程を含めた上での
加法性伝送路雑音をｖt とすると、時刻ｔでの受信信号
ｒt は、式（１）で示されるようにベクトルｈt とベク
トルｓt との畳込みと雑音の和で与えられる。
【００１５】 ｒt ＝ｓt T ・ｈt ＋ｖt （２） 以下、式（１）を時刻ｔでの伝送路方程式と呼ぶ。伝送
路方程式を時刻ｔ−Ｎ＋１から時刻ｔまでＮ個まとめる
と、式（３）のようにかける。
路方程式を時刻ｔ−Ｎ＋１から時刻ｔまでＮ個まとめる
と、式（３）のようにかける。
【００１６】 ｒt ＝Ｓt T ・ｈt ＋ｖt （３） ここで、送信信号行列Ｓt 、雑音ベクトルｖt は以下で
定義する。
定義する。
【００１７】
【００１８】 ｖt T ＝［ｖt ，ｖt - 1 ，…，ｖt - N + 1 ］ （５） 伝送路特性計算回路３０３は、（条件）Ｎシンボル（Ｎ
≧Ｌ＋１）間での伝送路の変動は無視できる、すなわ
ち、ｈt ＝ｈt - 1 ＝ｈt - N + 1 。という条件のもと
に、式（３）のＮ時刻分の伝送路方程式から伝送路イン
パルス応答ベクトルｈt を最小二乗推定する。具体的に
は、送信される可能性のある複数の信号系列、すなわち
（ｓt ，ｓt - 1 ，…，ｓt - L - N + 1 ）で考えられ
る全ての信号系列に対してそれぞれ送信信号行列Ｓt を
構成し式（６）によりそれぞれ伝送路インパルス応答ベ
クトルｈt , 1 s を求める。
≧Ｌ＋１）間での伝送路の変動は無視できる、すなわ
ち、ｈt ＝ｈt - 1 ＝ｈt - N + 1 。という条件のもと
に、式（３）のＮ時刻分の伝送路方程式から伝送路イン
パルス応答ベクトルｈt を最小二乗推定する。具体的に
は、送信される可能性のある複数の信号系列、すなわち
（ｓt ，ｓt - 1 ，…，ｓt - L - N + 1 ）で考えられ
る全ての信号系列に対してそれぞれ送信信号行列Ｓt を
構成し式（６）によりそれぞれ伝送路インパルス応答ベ
クトルｈt , 1 s を求める。
【００１９】 ｈt , 1 s ＝（Ｓt T ・Ｓt ）- 1 ・Ｓt T ・ｒt （６） 特に、インパルス応答推定に用いる受信信号の数（Ｎ）
が伝送路応答の数（Ｌ＋１）に等しいときは送信信号行
列Ｓt が正方行列となるので、受信信号に単に送信信号
行列Ｓt の逆行列をかけることで最小二乗推定による伝
送路応答推定値が得られる。
が伝送路応答の数（Ｌ＋１）に等しいときは送信信号行
列Ｓt が正方行列となるので、受信信号に単に送信信号
行列Ｓt の逆行列をかけることで最小二乗推定による伝
送路応答推定値が得られる。
【００２０】 ｈt , 1 s ＝Ｓt - 1 ・ｒt （７） これら推定値を基に、ブランチメトリック計算回路３０
４は、送信信号の全ての組み合わせ（ｓt ，ｓt - 1 ，
…，ｓt - L - N + 1 に対して式（８）に示す尤度（ブ
ランチメトリック）を毎時刻計算する。
４は、送信信号の全ての組み合わせ（ｓt ，ｓt - 1 ，
…，ｓt - L - N + 1 に対して式（８）に示す尤度（ブ
ランチメトリック）を毎時刻計算する。
【００２１】 Ｍt （ｓt ，ｓt - 1 ，…，ｓt - L - N + 1 ：ｓt + 1 ） ＝｜ｒt + 1 −ｓt + 1 T ・ｈt , 1 s ｜2 （８） そして、ビタビプロセッサ３０５は、この値の全時刻ｔ
に渡る和で定まる値（パストメリック）を最小にする全
時刻に渡る送信信号系列をビタビアルゴリズムにより求
める。
に渡る和で定まる値（パストメリック）を最小にする全
時刻に渡る送信信号系列をビタビアルゴリズムにより求
める。
【００２２】簡単のため、２値信号、２波モデル（Ｌ＝
１）において２個（Ｎ＝２）の受信信号より伝送路応答
を推定する場合を例にとると、基本形の動作は以下とな
る。伝送路応答はベクトルｈt T ＝［ｈt 0 、ｈt 1 ］
と書け、２時刻分の伝送路方程式は式（９）で書ける。
１）において２個（Ｎ＝２）の受信信号より伝送路応答
を推定する場合を例にとると、基本形の動作は以下とな
る。伝送路応答はベクトルｈt T ＝［ｈt 0 、ｈt 1 ］
と書け、２時刻分の伝送路方程式は式（９）で書ける。
【００２３】
【００２４】基本形では、ビタビアルゴリズムを動作さ
せるトレリス線図の状態は、送信信号行列Ｓt の成分に
現れる３個の送信信号の組み合わせ、すなわち（ｓt ，
ｓt -1 ，ｓt - 2 ）が定めている。２値の送信信号を
｛１、−１｝とし、それらを便宜上それぞれ｛１、０｝
で表現することにすると、送信される可能性のある信号
系列候補（ｓt ｓt - 1 ｓt - 2 ）は、（０００）、
（００１）、（０１０）、（０１１）、（１００）、
（１０１）、（１１０）、（１１１）、の８通りとな
る。したがって、基本形の状態数は８となり、その状態
遷移図は図４のように書くことができる。図４では、各
状態を表す楕円内の記号の組は慣例に従い、時刻の古い
順にｓt + 2 ｓt - 1 ｓt ）を記号表現しており、例え
ば状態（１１０）は、ｓt + 2 ＝１、ｓt - 1 ＝１、ｓ
t ＝０を表している。また、図４において、各状態を結
ぶ線は状態遷移を表し、状態遷移にともなって送られる
信号候補（以下、遷移信号と呼ぶ）は、状態遷移の右側
に置かれた状態（これを遷移先状態といい、状態遷移の
左側に置かれた状態を遷移元状態という）の最下位桁
（一番右の桁）に現れている。
せるトレリス線図の状態は、送信信号行列Ｓt の成分に
現れる３個の送信信号の組み合わせ、すなわち（ｓt ，
ｓt -1 ，ｓt - 2 ）が定めている。２値の送信信号を
｛１、−１｝とし、それらを便宜上それぞれ｛１、０｝
で表現することにすると、送信される可能性のある信号
系列候補（ｓt ｓt - 1 ｓt - 2 ）は、（０００）、
（００１）、（０１０）、（０１１）、（１００）、
（１０１）、（１１０）、（１１１）、の８通りとな
る。したがって、基本形の状態数は８となり、その状態
遷移図は図４のように書くことができる。図４では、各
状態を表す楕円内の記号の組は慣例に従い、時刻の古い
順にｓt + 2 ｓt - 1 ｓt ）を記号表現しており、例え
ば状態（１１０）は、ｓt + 2 ＝１、ｓt - 1 ＝１、ｓ
t ＝０を表している。また、図４において、各状態を結
ぶ線は状態遷移を表し、状態遷移にともなって送られる
信号候補（以下、遷移信号と呼ぶ）は、状態遷移の右側
に置かれた状態（これを遷移先状態といい、状態遷移の
左側に置かれた状態を遷移元状態という）の最下位桁
（一番右の桁）に現れている。
【００２５】伝送路特性計算回路３０３は、それぞれに
対して伝送路応答ベクトルｈt T ＝［ｈt 0 、ｈt １］
の最小二乗解を求め、それをベクトルｈt , 1 s （０，
０，０）、ベクトルｈt , 1 s （０，０，１）、ベクト
ルｈt , 1 s （０，１，０）、ベクトルｈ
t , 1 s （０，１，１）、ベクトルｈt , 1 s （１，
０，０）、ベクトルｈt , 1 s （１，０，１）、ベクト
ルｈt , 1 s （１，１，０）、ベクトルｈ
t , 1 s （１，１，１）とする。ベクトルｈ
t , 1 s （ｓt - 2 ，ｓt - 1，ｓt ）は、ｓt - 2 ，
ｓt - 1 ，ｓt をそれぞれ０、１のいずれかの値とした
ときの伝送路応答ベクトルｈt T ＝［ｈt 0 、ｈt 1 ］
を示す。
対して伝送路応答ベクトルｈt T ＝［ｈt 0 、ｈt １］
の最小二乗解を求め、それをベクトルｈt , 1 s （０，
０，０）、ベクトルｈt , 1 s （０，０，１）、ベクト
ルｈt , 1 s （０，１，０）、ベクトルｈ
t , 1 s （０，１，１）、ベクトルｈt , 1 s （１，
０，０）、ベクトルｈt , 1 s （１，０，１）、ベクト
ルｈt , 1 s （１，１，０）、ベクトルｈ
t , 1 s （１，１，１）とする。ベクトルｈ
t , 1 s （ｓt - 2 ，ｓt - 1，ｓt ）は、ｓt - 2 ，
ｓt - 1 ，ｓt をそれぞれ０、１のいずれかの値とした
ときの伝送路応答ベクトルｈt T ＝［ｈt 0 、ｈt 1 ］
を示す。
【００２６】ブランチメトリック計算回路３０４は、伝
送路特性計算回路３０３で求められたベクトルｈ
t , 1 s （ｓt - 2 ，ｓt - 1 ，ｓt ）と、時刻（ｔ＋
１）での受信信号ｒt + 1 から、以下の１６通りのブラ
ンチメトリックを求める。
送路特性計算回路３０３で求められたベクトルｈ
t , 1 s （ｓt - 2 ，ｓt - 1 ，ｓt ）と、時刻（ｔ＋
１）での受信信号ｒt + 1 から、以下の１６通りのブラ
ンチメトリックを求める。
【００２７】
【００２８】ここで、Ｍt （ｓt - 2 ，ｓt - 1 ，
ｓt ：ｓt + 1 ）は、ｓt ，ｓt - 1 ，ｓt - 2 とｓ
t + 1 をそれぞれ０、１のいずれかの値としたときブラ
ンチメトリック値である。これらは、図７に示す各状態
遷移、７ａ〜７ｐに対応しており、その対応関係は図１
２に詳細に示されている。これらのブランチメトリック
｛Ｍt（ｓt - 2 ，ｓt - 1 ，ｓt ：ｓt + 1 ）｝はビ
タビプロセッサ３０５に出力され、ビタビプロセッサ３
０５はメトリック最小のパスを選択することにより最尤
受信を行って、端子３０６から受信信号を主力する。こ
のビタビプロセッサ３０５の動作は、前述した文献１、
２に記述されているものと全く同一のものであるから、
その詳細な説明は省略する。
ｓt ：ｓt + 1 ）は、ｓt ，ｓt - 1 ，ｓt - 2 とｓ
t + 1 をそれぞれ０、１のいずれかの値としたときブラ
ンチメトリック値である。これらは、図７に示す各状態
遷移、７ａ〜７ｐに対応しており、その対応関係は図１
２に詳細に示されている。これらのブランチメトリック
｛Ｍt（ｓt - 2 ，ｓt - 1 ，ｓt ：ｓt + 1 ）｝はビ
タビプロセッサ３０５に出力され、ビタビプロセッサ３
０５はメトリック最小のパスを選択することにより最尤
受信を行って、端子３０６から受信信号を主力する。こ
のビタビプロセッサ３０５の動作は、前述した文献１、
２に記述されているものと全く同一のものであるから、
その詳細な説明は省略する。
【００２９】なお、ここではメトリツクを式（１０）式
に示すように二乗誤差を求めるとしたが、通常の最尤推
定に用いられるように式（１０）を展開し、文献１、２
に記述されているマッチドフィルタを用いた形式とした
り、さらに同時刻のすべてのブランチメトリックに共通
なｒt + 1 2 の項を省略したり、あるいは式（１０）の
符号を変えて最大メトリックを求めるようにしても同様
の効果が得られることは明らかである。
に示すように二乗誤差を求めるとしたが、通常の最尤推
定に用いられるように式（１０）を展開し、文献１、２
に記述されているマッチドフィルタを用いた形式とした
り、さらに同時刻のすべてのブランチメトリックに共通
なｒt + 1 2 の項を省略したり、あるいは式（１０）の
符号を変えて最大メトリックを求めるようにしても同様
の効果が得られることは明らかである。
【００３０】さて、基本形では、送信信号が２値信号、
伝送路が２波モデルで、観測する受信信号が２個の場合
でも上記のように８状態を要した。しかし、この状態数
は本発明を用いれば、４状態、さらには２状態にも削減
することができる。具体的には、本発明では、ビタビア
ルゴリズムを動作させるトレリス線図の状態を、送信信
号行列Ｓt の成分に現れる３個の送信視信号の組み合わ
せ（ｓt - 2 ｓt - 1ｓt ）でなく、それを１回縮退さ
せ、２個の送信信号の組み合わせ（ｓ1 t - 1ｓ1 t ）
により、あるいは２回縮退させ、１個の送信信号（ｓ2
t ）により各状態を表すようにする。以下、基本形の状
態を縮退系の状態と特に区別する意味で、（ｓt - 2 ｓ
t - 1 ｓt ）を（ｓ0 t - 2 ｓ0 t - 1 ｓ0 t ）とも書
くことにする。
伝送路が２波モデルで、観測する受信信号が２個の場合
でも上記のように８状態を要した。しかし、この状態数
は本発明を用いれば、４状態、さらには２状態にも削減
することができる。具体的には、本発明では、ビタビア
ルゴリズムを動作させるトレリス線図の状態を、送信信
号行列Ｓt の成分に現れる３個の送信視信号の組み合わ
せ（ｓt - 2 ｓt - 1ｓt ）でなく、それを１回縮退さ
せ、２個の送信信号の組み合わせ（ｓ1 t - 1ｓ1 t ）
により、あるいは２回縮退させ、１個の送信信号（ｓ2
t ）により各状態を表すようにする。以下、基本形の状
態を縮退系の状態と特に区別する意味で、（ｓt - 2 ｓ
t - 1 ｓt ）を（ｓ0 t - 2 ｓ0 t - 1 ｓ0 t ）とも書
くことにする。
【００３１】状態の縮退方法は次の２通りの方法があ
る。 １）単純縮退：単純縮退は、状態の記号表現の下位桁の
共通性により状態を縮退させていく方法である。図４の
基本形状態の記号表現の上位桁を１桁ずつ消去していく
ことで実現される。ｋ回縮退すれば、状態の記号表現は
ｋ桁減少する。見方を変えれば、基本形状態を構成する
送信信号候補から複数の最近の送信信号候補を取り出
し、それを縮退状態とする方法である。例えば、基本形
状態（１１０）に対する縮退状態は、１回の縮退（これ
を１次縮退と呼ぶ）では（１１０）→（１０）、２回の
縮退（これを２次縮退と呼ぶ）では（１０）→（０）で
与えられる。逆に、１次縮退状態（００）は基本形状態
の（０００）と（１００）の２つに、２次縮退状態
（０）は４つの基本形状態（０００）、（１００）、
（０１０）、（１１０）に対応する可能性がある。しか
し、各縮退状態は時刻を指定すれば対応する基本形状態
の内の１つにしか対応しない。
る。 １）単純縮退：単純縮退は、状態の記号表現の下位桁の
共通性により状態を縮退させていく方法である。図４の
基本形状態の記号表現の上位桁を１桁ずつ消去していく
ことで実現される。ｋ回縮退すれば、状態の記号表現は
ｋ桁減少する。見方を変えれば、基本形状態を構成する
送信信号候補から複数の最近の送信信号候補を取り出
し、それを縮退状態とする方法である。例えば、基本形
状態（１１０）に対する縮退状態は、１回の縮退（これ
を１次縮退と呼ぶ）では（１１０）→（１０）、２回の
縮退（これを２次縮退と呼ぶ）では（１０）→（０）で
与えられる。逆に、１次縮退状態（００）は基本形状態
の（０００）と（１００）の２つに、２次縮退状態
（０）は４つの基本形状態（０００）、（１００）、
（０１０）、（１１０）に対応する可能性がある。しか
し、各縮退状態は時刻を指定すれば対応する基本形状態
の内の１つにしか対応しない。
【００３２】図５に、単純縮退１回により図４のトレリ
ス線図を縮退した場合のトレリス線図を示す。図５で
は、各縮退状態を表す楕円内の記号の組は慣例に従い時
刻の古い順に（ｓt - 1 ｓt ）を記号表現しており、例
えば状態（１０）は、ｓt - 1＝１、ｓt ＝０を表して
いる。また、各縮退状態の楕円内の括弧内で示されてい
るのは対応する２つの基本形状態である。また、図５の
トレリス線図内の状態遷移と図４の基本形トレリス線図
内の状態遷移との対応関係は、図１０に示される。
ス線図を縮退した場合のトレリス線図を示す。図５で
は、各縮退状態を表す楕円内の記号の組は慣例に従い時
刻の古い順に（ｓt - 1 ｓt ）を記号表現しており、例
えば状態（１０）は、ｓt - 1＝１、ｓt ＝０を表して
いる。また、各縮退状態の楕円内の括弧内で示されてい
るのは対応する２つの基本形状態である。また、図５の
トレリス線図内の状態遷移と図４の基本形トレリス線図
内の状態遷移との対応関係は、図１０に示される。
【００３３】また、単純縮退状態間の遷移信号は対応す
る基本形状態によらず同一であるので、縮退した状態を
記憶しておけば、遷移信号を特定できる。例えば、１次
縮退状態（００）から（００）への遷移は、遷移元状態
に対応する基本形状態が（０００）、（１００）に関わ
らず０である。 ２）回転縮退：基本形状態を表す信号系列（ｓ0 t - 2
ｓ0 t - 1 ｓ0 t ）内での位相回転が同じものどうしを
同一の縮退状態に対応させる方法もある。回転縮退の縮
退状態は漸化的に与えられる。
る基本形状態によらず同一であるので、縮退した状態を
記憶しておけば、遷移信号を特定できる。例えば、１次
縮退状態（００）から（００）への遷移は、遷移元状態
に対応する基本形状態が（０００）、（１００）に関わ
らず０である。 ２）回転縮退：基本形状態を表す信号系列（ｓ0 t - 2
ｓ0 t - 1 ｓ0 t ）内での位相回転が同じものどうしを
同一の縮退状態に対応させる方法もある。回転縮退の縮
退状態は漸化的に与えられる。
【００３４】１次縮退： （ｓ1 t - 1 ｓ1 t ）＝（ｓ
0 t - 2 ’＋ｓ0 t - 1 ｓ0 t - 2’＋ｓ0 t ） ただし、加算はＭ値信号の場合、ｍｏｄ Ｍの和。
0 t - 2 ’＋ｓ0 t - 1 ｓ0 t - 2’＋ｓ0 t ） ただし、加算はＭ値信号の場合、ｍｏｄ Ｍの和。
【００３５】また、ｓ0 t - 2 ’＋ｓ0 t - 2 ＝０
（ｍｏｄ Ｍ） ２次縮退： （ｓ2 t ）＝（ｓ1 t - 1 ’＋ｓ1 t ） ただし、加算はＭ値信号の場合、ｍｏｄ Ｍの和。
（ｍｏｄ Ｍ） ２次縮退： （ｓ2 t ）＝（ｓ1 t - 1 ’＋ｓ1 t ） ただし、加算はＭ値信号の場合、ｍｏｄ Ｍの和。
【００３６】また、ｓ1 t - 1 ’＋ｓ1 t - 1 ＝０
（ｍｏｄ Ｍ） 例えば、基本形状態（１１０）に対する回転縮退状態
は、１次縮退では（１１０）→（０１）に、２次縮退で
は（０１）→（１）で与えられる。逆に、１次縮退状態
（００）は基本形状態の（０００）と（１１１）の２つ
に、２次縮退状態（０）は４つの基本形状態（００
０）、（１１１）、（０１１）、（１００）に対応する
可能性がある。しかし、各縮退状態は時刻を指定すれば
対応する基本形状態の内の１つにしか対応しない。図６
に、回転縮退により図４のトレリス線図を縮退した場合
のトレリス線図を示す。回転縮退状態間の遷移シンボル
は対応する基本形状態間の遷移シンボルにより決まるの
で、回転縮退状態だけでは一意には決まない。例えば、
１次縮退状態（００）から（００）への遷移は、遷移元
状態に対応する基本形状態が（０００）の場合は０であ
るが、（１１１）の場合は１である。そのため、遷移信
号の特定は各縮退状態に対して記憶する基本形状態から
行う。
（ｍｏｄ Ｍ） 例えば、基本形状態（１１０）に対する回転縮退状態
は、１次縮退では（１１０）→（０１）に、２次縮退で
は（０１）→（１）で与えられる。逆に、１次縮退状態
（００）は基本形状態の（０００）と（１１１）の２つ
に、２次縮退状態（０）は４つの基本形状態（００
０）、（１１１）、（０１１）、（１００）に対応する
可能性がある。しかし、各縮退状態は時刻を指定すれば
対応する基本形状態の内の１つにしか対応しない。図６
に、回転縮退により図４のトレリス線図を縮退した場合
のトレリス線図を示す。回転縮退状態間の遷移シンボル
は対応する基本形状態間の遷移シンボルにより決まるの
で、回転縮退状態だけでは一意には決まない。例えば、
１次縮退状態（００）から（００）への遷移は、遷移元
状態に対応する基本形状態が（０００）の場合は０であ
るが、（１１１）の場合は１である。そのため、遷移信
号の特定は各縮退状態に対して記憶する基本形状態から
行う。
【００３７】縮退したトレリス線図は状態数が少ないの
で、トレリス線図全体での遷移の数も少ない。このトレ
リス線図を基にビタビアルゴリズムを動作させれば、ビ
タビプロセッサの処理量を指数関数的に減少させること
ができる。
で、トレリス線図全体での遷移の数も少ない。このトレ
リス線図を基にビタビアルゴリズムを動作させれば、ビ
タビプロセッサの処理量を指数関数的に減少させること
ができる。
【００３８】一方、系列推定装置に高速追従性を可能と
させる伝送路応答計算は、縮退されたトレリス線図上で
は次のように行う。まず、ビタビアルゴリズムにおい
て、生き残り系列の履歴から各時刻の各縮退状態が実際
にどのような基本形状態と対応しているかを記憶するよ
うにする。そして、伝送路応答計算は、この基本形状態
に対応する送信信号の組合せ（ｓt ，ｓt - 1 ，…，ｓ
t - L - N + 1 ）に対して伝送路方程式を解き、得られ
た伝送路応答を当該縮退状態に対する伝送路応答推定値
としてブランチメトリック計算に用いるようにする。時
刻（ｔ＋１）のブランチメトリック計算において、仮想
受信信号点を求めるための送信信号系列候補は、時刻ｔ
の縮退状態に対応する非縮退状態に関する送信信号信号
から最古信号を除いた系列（ｓt ，ｓt - 1 ，…，ｓ
t - L - N ）と時刻（ｔ＋１）の新な送信信号候補との
組み合わせにより与えられる。すなわち、先に求めた伝
送路応答推定値と（ｓt + 1 ，ｓt ，…，
ｓt - L - N ）とから仮想受信信号点を求め、これと実
際の受信信号との距離を計算しブランチメトリックとす
る。この伝送路応答計算とブランチメトリック計算と
が、縮退されたトレリス線図を基に動作させる本発明の
系列推定装置においても直接解による毎時刻の伝送路推
定を可能としている。また、本発明においては、伝送路
応答計算は基本形状態の全てに対してではなく、縮退さ
れた状態の数だけ行えばよく、伝送路応答推定に関わる
計算量も指数関数的に減少することになるという利点が
ある。
させる伝送路応答計算は、縮退されたトレリス線図上で
は次のように行う。まず、ビタビアルゴリズムにおい
て、生き残り系列の履歴から各時刻の各縮退状態が実際
にどのような基本形状態と対応しているかを記憶するよ
うにする。そして、伝送路応答計算は、この基本形状態
に対応する送信信号の組合せ（ｓt ，ｓt - 1 ，…，ｓ
t - L - N + 1 ）に対して伝送路方程式を解き、得られ
た伝送路応答を当該縮退状態に対する伝送路応答推定値
としてブランチメトリック計算に用いるようにする。時
刻（ｔ＋１）のブランチメトリック計算において、仮想
受信信号点を求めるための送信信号系列候補は、時刻ｔ
の縮退状態に対応する非縮退状態に関する送信信号信号
から最古信号を除いた系列（ｓt ，ｓt - 1 ，…，ｓ
t - L - N ）と時刻（ｔ＋１）の新な送信信号候補との
組み合わせにより与えられる。すなわち、先に求めた伝
送路応答推定値と（ｓt + 1 ，ｓt ，…，
ｓt - L - N ）とから仮想受信信号点を求め、これと実
際の受信信号との距離を計算しブランチメトリックとす
る。この伝送路応答計算とブランチメトリック計算と
が、縮退されたトレリス線図を基に動作させる本発明の
系列推定装置においても直接解による毎時刻の伝送路推
定を可能としている。また、本発明においては、伝送路
応答計算は基本形状態の全てに対してではなく、縮退さ
れた状態の数だけ行えばよく、伝送路応答推定に関わる
計算量も指数関数的に減少することになるという利点が
ある。
【００３９】以上により、基本形ブラインドビタビ等化
器に比べより少ない処理量、記憶量で、基本形と同様の
伝送路推定機能を有する系列推定装置が実現できる。
器に比べより少ない処理量、記憶量で、基本形と同様の
伝送路推定機能を有する系列推定装置が実現できる。
【００４０】図１に第１の発明に係る系列推定装置の一
実施例を示す。図１の系列推定装置は、図５に示す単純
縮退されたトレリス線図を基にビタビアルゴリズムを動
作させる。
実施例を示す。図１の系列推定装置は、図５に示す単純
縮退されたトレリス線図を基にビタビアルゴリズムを動
作させる。
【００４１】入力端子１０１に供給された時刻ｔでの受
信器入力ｒt は、レジスタ１０２に記憶されるとともに
ブランチメトリック計算回路１０４に送られる。レジス
タ１０２に記憶された時刻ｔ−Ｎ＋１から時刻ｔまでＮ
個の受信信号信号は伝送路応答計算回路１０３に入力さ
れる。伝送路応答計算回路１０３は、各時刻で図５の各
縮退状態に対応する基本形状態に対してのみ式（６）あ
るいは式（７）にしたがってベクトルｈt , l s を計算
する。例えば、図８のように、時刻ｔでの（００）、
（０１）、（１０）、（１１）の各縮退状態に対応する
基本形状態がそれぞれ（１００）、（１０１）、（０１
０）、（０１１）（下線が引かれた基本形状態）とする
と、この４つの信号系列に対してのみに伝送路応答計算
を実施する。これに対して図３の従来の系列推定装置で
は、（０００）から（１１１）の８つの信号系列（全て
の基本形状態）に対して計算を実行している。
信器入力ｒt は、レジスタ１０２に記憶されるとともに
ブランチメトリック計算回路１０４に送られる。レジス
タ１０２に記憶された時刻ｔ−Ｎ＋１から時刻ｔまでＮ
個の受信信号信号は伝送路応答計算回路１０３に入力さ
れる。伝送路応答計算回路１０３は、各時刻で図５の各
縮退状態に対応する基本形状態に対してのみ式（６）あ
るいは式（７）にしたがってベクトルｈt , l s を計算
する。例えば、図８のように、時刻ｔでの（００）、
（０１）、（１０）、（１１）の各縮退状態に対応する
基本形状態がそれぞれ（１００）、（１０１）、（０１
０）、（０１１）（下線が引かれた基本形状態）とする
と、この４つの信号系列に対してのみに伝送路応答計算
を実施する。これに対して図３の従来の系列推定装置で
は、（０００）から（１１１）の８つの信号系列（全て
の基本形状態）に対して計算を実行している。
【００４２】行列Ｓt T Ｓt あるいは行列Ｓt が特異と
なる送信信号の組み合わせに対しては、伝送路応答計算
回路１０３はあらかじめ定めた値（例えば０）を出力
し、ベクトルｈt , l s が不定である旨を伝送路応答検
査回路１０５に伝える。ここでは、各状態に対して不定
でない場合に得られる伝送路応答推定値を正当な推定値
と呼ぶ。伝送路応答検査回路１０５は、伝送路応答計算
回路１０３が計算した個々の伝送路応答ベクトルｈ
t , l s が不定であるか否かを調べる。不定である状態
に対しては当該状態に遷移している生き残り系列が前時
刻の状態において採用した伝送路応答推定値を持って現
時刻の伝送路応答推定値に代える。例えば、時刻ｔでの
縮退状態（０１）に対応する基本形状態が（１０１）の
とき、信号行列Ｓt T Ｓt あるいは行列Ｓt は特異とな
る。このとき、時刻ｔの縮退状態（０１）に対する伝送
路応答推定値としては、当該縮退状態（０１）に遷移し
ている生残りパスの前時刻（ｔ−１）での縮退状態（１
０）（なぜなら、時刻ｔの基本形状態の記憶値が（１０
１）であるから、時刻（ｔ−１）での縮退状態が（１
０）であることが分かる）が採用した伝送路応答推定値
で代替する。この前時刻の伝送路応答推定値を伝送路応
答記憶回路１０６が供給する。一方、不定でない場合
は、伝送路応答計算回路１０３で計算した現時刻の伝送
路応答をそのまま用いるようにする。こうして得られた
現時刻の伝送路応答推定値は、その値を持つ基本形状態
とともに伝送路応答記憶回路１０６に記憶される。
なる送信信号の組み合わせに対しては、伝送路応答計算
回路１０３はあらかじめ定めた値（例えば０）を出力
し、ベクトルｈt , l s が不定である旨を伝送路応答検
査回路１０５に伝える。ここでは、各状態に対して不定
でない場合に得られる伝送路応答推定値を正当な推定値
と呼ぶ。伝送路応答検査回路１０５は、伝送路応答計算
回路１０３が計算した個々の伝送路応答ベクトルｈ
t , l s が不定であるか否かを調べる。不定である状態
に対しては当該状態に遷移している生き残り系列が前時
刻の状態において採用した伝送路応答推定値を持って現
時刻の伝送路応答推定値に代える。例えば、時刻ｔでの
縮退状態（０１）に対応する基本形状態が（１０１）の
とき、信号行列Ｓt T Ｓt あるいは行列Ｓt は特異とな
る。このとき、時刻ｔの縮退状態（０１）に対する伝送
路応答推定値としては、当該縮退状態（０１）に遷移し
ている生残りパスの前時刻（ｔ−１）での縮退状態（１
０）（なぜなら、時刻ｔの基本形状態の記憶値が（１０
１）であるから、時刻（ｔ−１）での縮退状態が（１
０）であることが分かる）が採用した伝送路応答推定値
で代替する。この前時刻の伝送路応答推定値を伝送路応
答記憶回路１０６が供給する。一方、不定でない場合
は、伝送路応答計算回路１０３で計算した現時刻の伝送
路応答をそのまま用いるようにする。こうして得られた
現時刻の伝送路応答推定値は、その値を持つ基本形状態
とともに伝送路応答記憶回路１０６に記憶される。
【００４３】縮退状態に対して求めた現時刻の伝送路応
答推定値ベクトルｈt ’は、ブランチメトリック計算回
路１０４に供給される。ブランチメトリック計算回路１
０４は、伝送路応答検査回路１０５が与える現時刻の伝
送路応答推定値ベクトルｈt’を基に、縮退状態からの
全ての遷移に対して式（８）で定まるブランチメトリッ
クを個別に計算する。ここで、時刻（ｔ＋１）のブラン
チメトリック計算において、仮想受信信号点を求めるた
めの送信信号系列候補は、時刻ｔの縮退状態に対応する
基本形状態に関する送信信号信号から最古信号を除いた
系列（ｓt ，ｓt - 1 ，…，ｓt - L - N ）と時刻（ｔ
＋１）の新たな送信信号候補との組み合わせにより与え
られる。すなわち、先に求めた伝送路応答推定値と（ｓ
t + 1 ，ｓt ，…，ｓt - L - N ）とから仮想受信信号
点を求め、これと実際の受信信号との距離を計算しブラ
ンチメトリックとする。時刻ｔの縮退状態を（０１）の
ときは、最新信号１と時刻（ｔ＋１）の新たな送信信号
候補０または１との組み合わせ、すなわち系列１０また
は系列１１に対して仮想受信点を求める。
答推定値ベクトルｈt ’は、ブランチメトリック計算回
路１０４に供給される。ブランチメトリック計算回路１
０４は、伝送路応答検査回路１０５が与える現時刻の伝
送路応答推定値ベクトルｈt’を基に、縮退状態からの
全ての遷移に対して式（８）で定まるブランチメトリッ
クを個別に計算する。ここで、時刻（ｔ＋１）のブラン
チメトリック計算において、仮想受信信号点を求めるた
めの送信信号系列候補は、時刻ｔの縮退状態に対応する
基本形状態に関する送信信号信号から最古信号を除いた
系列（ｓt ，ｓt - 1 ，…，ｓt - L - N ）と時刻（ｔ
＋１）の新たな送信信号候補との組み合わせにより与え
られる。すなわち、先に求めた伝送路応答推定値と（ｓ
t + 1 ，ｓt ，…，ｓt - L - N ）とから仮想受信信号
点を求め、これと実際の受信信号との距離を計算しブラ
ンチメトリックとする。時刻ｔの縮退状態を（０１）の
ときは、最新信号１と時刻（ｔ＋１）の新たな送信信号
候補０または１との組み合わせ、すなわち系列１０また
は系列１１に対して仮想受信点を求める。
【００４４】ブランチメトリック計算回路１０４は図１
３に示される遷移に対して計算したブランチメトリック
値８通りをビタビプロセッサ１０７に出力する。ビタビ
プロセッサ１０７は、ビタビアルゴリズムにより式
（８）のメトリックの全ての時刻の和が最小となる系列
を探し、判定出力を出力端子１０９に出力する。ビタビ
プロセッサ１０７は、同時に、各時刻での縮退状態に対
応する基本形状態を生残りパスの履歴より調べ、それを
非縮退状態記憶回路１０８に出力する。例えば、時刻
（ｔ＋１）の縮退状態（１１）が前出の時刻ｔの縮退状
態（０１）から遷移した場合、縮退状態（１１）に対す
る基本形状態は（０１１）となり、この値が縮退状態
（１１）の非縮退状態として非縮退状態記憶回路１０８
に登録される。そして、この基本形状態の記憶値は、時
刻（ｔ＋１）の伝送路応答計算と時刻（ｔ＋２）のブラ
ンチメトリック計算とに利用される。ビタビプロセッサ
１０７の動作は縮退されたトレリス線図を基にしている
ことを除いて文献１、２の系列推定装置と同様であるの
で、詳細は省略する。
３に示される遷移に対して計算したブランチメトリック
値８通りをビタビプロセッサ１０７に出力する。ビタビ
プロセッサ１０７は、ビタビアルゴリズムにより式
（８）のメトリックの全ての時刻の和が最小となる系列
を探し、判定出力を出力端子１０９に出力する。ビタビ
プロセッサ１０７は、同時に、各時刻での縮退状態に対
応する基本形状態を生残りパスの履歴より調べ、それを
非縮退状態記憶回路１０８に出力する。例えば、時刻
（ｔ＋１）の縮退状態（１１）が前出の時刻ｔの縮退状
態（０１）から遷移した場合、縮退状態（１１）に対す
る基本形状態は（０１１）となり、この値が縮退状態
（１１）の非縮退状態として非縮退状態記憶回路１０８
に登録される。そして、この基本形状態の記憶値は、時
刻（ｔ＋１）の伝送路応答計算と時刻（ｔ＋２）のブラ
ンチメトリック計算とに利用される。ビタビプロセッサ
１０７の動作は縮退されたトレリス線図を基にしている
ことを除いて文献１、２の系列推定装置と同様であるの
で、詳細は省略する。
【００４５】図２に第２の発明に係る系列推定装置の一
実施例を示す。２値信号、２波モデル（Ｌ＝１）におい
て２個（Ｎ＝２）の受信信号より伝送路応答を推定する
場合の例では、図２の系列推定装置は、図６に示す縮退
されたトレリス線図を基にビタビアルゴリズムを動作さ
せる。縮退方法は異なるが、ビタビアルゴリズムを縮退
状態に基づくトレリス線図に対して動作させることは第
１の発明と共通である。
実施例を示す。２値信号、２波モデル（Ｌ＝１）におい
て２個（Ｎ＝２）の受信信号より伝送路応答を推定する
場合の例では、図２の系列推定装置は、図６に示す縮退
されたトレリス線図を基にビタビアルゴリズムを動作さ
せる。縮退方法は異なるが、ビタビアルゴリズムを縮退
状態に基づくトレリス線図に対して動作させることは第
１の発明と共通である。
【００４６】入力端子２０１に供給された時刻ｔでの受
信器入力ｒt は、レジスタ２０２に記憶されるとともに
ブランチメトリック計算回路２０４に送られる。レジス
タ２０２に記憶された時刻ｔ−Ｎ＋１から時刻ｔまでＮ
個の受信信号信号は伝送路応答計算回路２０３に入力さ
れる。伝送路応答計算回路２０３は、図１の伝送路応答
計算回路１０３と異なり、あらかじめ定められた１次縮
退状態に対応する基本形状態に対してのみ式（６）ある
いは式（７）にしたがってベクトルｈt , l sを計算す
る。例えば、図９のように、時刻ｔでの（００）、（０
１）、（１０）、（１１）の各縮退状態に対応する基本
形状態がそれぞれ（１１１）、（００１）、（１０
１）、（１００）（下線が引かれた基本形状態）であっ
ても、（０００）、（００１）、（０１０）、（０１
１）というあらかじめ定めた４つの信号系列に対しての
みに伝送路応答計算を実施すればよい。
信器入力ｒt は、レジスタ２０２に記憶されるとともに
ブランチメトリック計算回路２０４に送られる。レジス
タ２０２に記憶された時刻ｔ−Ｎ＋１から時刻ｔまでＮ
個の受信信号信号は伝送路応答計算回路２０３に入力さ
れる。伝送路応答計算回路２０３は、図１の伝送路応答
計算回路１０３と異なり、あらかじめ定められた１次縮
退状態に対応する基本形状態に対してのみ式（６）ある
いは式（７）にしたがってベクトルｈt , l sを計算す
る。例えば、図９のように、時刻ｔでの（００）、（０
１）、（１０）、（１１）の各縮退状態に対応する基本
形状態がそれぞれ（１１１）、（００１）、（１０
１）、（１００）（下線が引かれた基本形状態）であっ
ても、（０００）、（００１）、（０１０）、（０１
１）というあらかじめ定めた４つの信号系列に対しての
みに伝送路応答計算を実施すればよい。
【００４７】これは、以下の理由による。例えば、図６
において、非縮退状態の（０００）と（１１１）は縮退
状態（００）を共有するが、（０００）を用いて計算し
て得られる伝送路応答ベクトルｈt ,l s（０００）と
（１１１）を用いて計算して得られる伝送路応答ベクト
ルｈt , l s （１１１）とは１８０度回転したものとな
る。（００）からの（０１）への遷移において、遷移信
号は基本形状態が（０００）のときは１、（１１１）の
ときは０であるから、１８０度回転した伝送路応答推定
値との間で仮想受信信号点を求めると全く同一の値とな
る。そのため、同じメトリックを計算することになる。
において、非縮退状態の（０００）と（１１１）は縮退
状態（００）を共有するが、（０００）を用いて計算し
て得られる伝送路応答ベクトルｈt ,l s（０００）と
（１１１）を用いて計算して得られる伝送路応答ベクト
ルｈt , l s （１１１）とは１８０度回転したものとな
る。（００）からの（０１）への遷移において、遷移信
号は基本形状態が（０００）のときは１、（１１１）の
ときは０であるから、１８０度回転した伝送路応答推定
値との間で仮想受信信号点を求めると全く同一の値とな
る。そのため、同じメトリックを計算することになる。
【００４８】したがって、回転縮退状態では、伝送路応
答計算のためには各１次縮退状態に対してどれか一つの
基本形状態に対する伝送路応答を固定的に（（００）に
対しては（０００）などというように）求めておけばよ
い。この変換規則例としては、図１６，１７がある。状
態変換回路２０９は、あらかじめ定めた変換規則に従
い、基本形状態を同じ１次縮退状態を有する基本形状態
に変換する回路であり、伝送路応答計算回路２０３はこ
れらあらかじめ定められた基本形状態に対してのみ伝送
路応答を計算すればよい。
答計算のためには各１次縮退状態に対してどれか一つの
基本形状態に対する伝送路応答を固定的に（（００）に
対しては（０００）などというように）求めておけばよ
い。この変換規則例としては、図１６，１７がある。状
態変換回路２０９は、あらかじめ定めた変換規則に従
い、基本形状態を同じ１次縮退状態を有する基本形状態
に変換する回路であり、伝送路応答計算回路２０３はこ
れらあらかじめ定められた基本形状態に対してのみ伝送
路応答を計算すればよい。
【００４９】行列Ｓt T Ｓt あるいは行列Ｓt が特異と
なる送信信号の組み合わせに対しては、伝送路応答計算
回路２０３はあらかじめ定めた値（例えば０）を出力
し、ベクトルｈt , l s が不定である旨を伝送路応答変
換回路２０５に伝える。伝送路応答変換回路２０５は、
伝送路応答計算回路２０３が計算した個々の伝送路応答
ベクトルｈt , l s が不定であるか否かを調べる。不定
でない場合は、伝送路応答計算回路２０３で計算した現
時刻の伝送路応答をそのまま用いるようにする。一方、
伝送路応答計算回路２０３の出力が不当（特異）の場合
の処理は伝送路応答検査回路１０５と異なる。不定であ
る状態に対しては、当該状態に遷移している生き残り系
列が前時刻の状態において採用した伝送路応答推定値を
基に現時刻の伝送路応答推定値を計算する。例えば、図
９において、時刻ｔでの縮退状態（００）に対応する基
本形状態がもし（０００）のとき、信号行列Ｓt T Ｓt
あるいは行列Ｓt は特異となる。図１６の変換規則に従
うとすると、時刻ｔの縮退状態（００）に対する伝送路
応答推定値としては基本形状態（０００）に対する応答
を計算する必要がある。また、当該縮退状態（００）に
遷移する生残りパスの前時刻ｔ−１での縮退状態が（１
１）であったとすると、伝送路応答記憶回路２０６が記
憶している伝送路応答は基本形状態（０１１）に対する
値である。ここで、基本形状態（０１１）から基本形状
態（０００）には遷移は存在せず、存在するのは基本形
状態（０１１）から基本形状態（１１１）である。基本
形状態（１１１）は基本形状態（０１１）の応答をその
まま継承してもよい。基本形状態（０００）の伝送路応
答は基本形状態（１１１）の伝送路応答を１８０度回転
したものであるから、基本形状態（０００）に対する応
答としては、前時刻の基本形状態（０１１）に対する値
を１８０度回転させる必要がある。このように、回転縮
退では、特異時に、遷移に応じて前時刻伝送路応答を回
転させた値を用いる必要がある。こうして得られた現時
刻の伝送路応答推定値は、伝送路応答記憶回路２０６に
記憶される。
なる送信信号の組み合わせに対しては、伝送路応答計算
回路２０３はあらかじめ定めた値（例えば０）を出力
し、ベクトルｈt , l s が不定である旨を伝送路応答変
換回路２０５に伝える。伝送路応答変換回路２０５は、
伝送路応答計算回路２０３が計算した個々の伝送路応答
ベクトルｈt , l s が不定であるか否かを調べる。不定
でない場合は、伝送路応答計算回路２０３で計算した現
時刻の伝送路応答をそのまま用いるようにする。一方、
伝送路応答計算回路２０３の出力が不当（特異）の場合
の処理は伝送路応答検査回路１０５と異なる。不定であ
る状態に対しては、当該状態に遷移している生き残り系
列が前時刻の状態において採用した伝送路応答推定値を
基に現時刻の伝送路応答推定値を計算する。例えば、図
９において、時刻ｔでの縮退状態（００）に対応する基
本形状態がもし（０００）のとき、信号行列Ｓt T Ｓt
あるいは行列Ｓt は特異となる。図１６の変換規則に従
うとすると、時刻ｔの縮退状態（００）に対する伝送路
応答推定値としては基本形状態（０００）に対する応答
を計算する必要がある。また、当該縮退状態（００）に
遷移する生残りパスの前時刻ｔ−１での縮退状態が（１
１）であったとすると、伝送路応答記憶回路２０６が記
憶している伝送路応答は基本形状態（０１１）に対する
値である。ここで、基本形状態（０１１）から基本形状
態（０００）には遷移は存在せず、存在するのは基本形
状態（０１１）から基本形状態（１１１）である。基本
形状態（１１１）は基本形状態（０１１）の応答をその
まま継承してもよい。基本形状態（０００）の伝送路応
答は基本形状態（１１１）の伝送路応答を１８０度回転
したものであるから、基本形状態（０００）に対する応
答としては、前時刻の基本形状態（０１１）に対する値
を１８０度回転させる必要がある。このように、回転縮
退では、特異時に、遷移に応じて前時刻伝送路応答を回
転させた値を用いる必要がある。こうして得られた現時
刻の伝送路応答推定値は、伝送路応答記憶回路２０６に
記憶される。
【００５０】縮退状態に対して求めた現時刻の伝送路応
答推定値ベクトルｈt ’は、ブランチメトリック計算回
路２０４に供給される。ブランチメトリック計算回路２
０４は、伝送路応答変換回路２０５が与える現時刻の伝
送路応答推定値ベクトルｈt’を基に、縮退状態からの
全ての遷移に対して式（８）で定まるブランチメトリッ
クを個別に計算する。ここで、時刻（ｔ＋１）のブラン
チメトリック計算において、仮想受信信号点を求めるた
めの送信信号系列候補は、時刻ｔの縮退状態に対応する
基本形状態に関する送信信号信号から最古信号を除いた
系列（ｓt ，ｓt - 1 ，…，ｓt - L - N ）と時刻（ｔ
＋１）の新たな送信信号候補との組み合わせにより与え
られる。すなわち、先に求めた伝送路応答推定値と（ｓ
t + 1 ，ｓt ，…，ｓt - L - N ）とから仮想受信信号
点を求め、これと実際の受信信号との距離を計算しブラ
ンチメトリックとする。時刻ｔの縮退状態を（１０）に
対しては、最新信号０と時刻（ｔ＋１）の新たな送信信
号候補０または１との組み合わせ、すなわち系列００ま
たは０１に対して仮想受信点を求める。
答推定値ベクトルｈt ’は、ブランチメトリック計算回
路２０４に供給される。ブランチメトリック計算回路２
０４は、伝送路応答変換回路２０５が与える現時刻の伝
送路応答推定値ベクトルｈt’を基に、縮退状態からの
全ての遷移に対して式（８）で定まるブランチメトリッ
クを個別に計算する。ここで、時刻（ｔ＋１）のブラン
チメトリック計算において、仮想受信信号点を求めるた
めの送信信号系列候補は、時刻ｔの縮退状態に対応する
基本形状態に関する送信信号信号から最古信号を除いた
系列（ｓt ，ｓt - 1 ，…，ｓt - L - N ）と時刻（ｔ
＋１）の新たな送信信号候補との組み合わせにより与え
られる。すなわち、先に求めた伝送路応答推定値と（ｓ
t + 1 ，ｓt ，…，ｓt - L - N ）とから仮想受信信号
点を求め、これと実際の受信信号との距離を計算しブラ
ンチメトリックとする。時刻ｔの縮退状態を（１０）に
対しては、最新信号０と時刻（ｔ＋１）の新たな送信信
号候補０または１との組み合わせ、すなわち系列００ま
たは０１に対して仮想受信点を求める。
【００５１】ブランチメトリック計算回路２０４は図１
４に示される遷移に対して計算したブランチメトリック
値８通りをビタビプロセッサ２０７に出力する。ビタビ
プロセッサ２０７は、ビタビアルゴリズムにより式
（８）のメトリックの全ての時刻の和が最小となる系列
を探し、判定出力を出力端子２１０に出力する。
４に示される遷移に対して計算したブランチメトリック
値８通りをビタビプロセッサ２０７に出力する。ビタビ
プロセッサ２０７は、ビタビアルゴリズムにより式
（８）のメトリックの全ての時刻の和が最小となる系列
を探し、判定出力を出力端子２１０に出力する。
【００５２】回転縮退状態では、縮退状態間の遷移シン
ボルは基本形状態間の遷移シンボルに依存して決まる。
遷移信号の特定は、まずメトリックの全ての時刻の和が
最小となる系列に対して現時刻よりある特定時間だけ遡
った時刻の縮退状態を調べ、次に、その過去の縮退状態
に対応する非縮退状態を非縮退状態記憶回路２０８から
検索し、最後にその非縮退状態が有する遷移信号値を出
力端子２１０に出力することにより行う。
ボルは基本形状態間の遷移シンボルに依存して決まる。
遷移信号の特定は、まずメトリックの全ての時刻の和が
最小となる系列に対して現時刻よりある特定時間だけ遡
った時刻の縮退状態を調べ、次に、その過去の縮退状態
に対応する非縮退状態を非縮退状態記憶回路２０８から
検索し、最後にその非縮退状態が有する遷移信号値を出
力端子２１０に出力することにより行う。
【００５３】ビタビプロセッサ２０７は、同時に、各時
刻での縮退状態に対応する基本形状態を生残りパスの履
歴より調べ、それを非縮退状態記憶回路２０８に出力す
る。例えば、時刻（ｔ＋１）の縮退状態（１１）が前出
の時刻ｔの縮退状態（０１）から遷移した場合、縮退状
態（１１）に対する基本形状態は（０１１）となり、こ
の値が縮退状態（１１）の非縮退状態として非縮退状態
記憶回路２０８に登録される。この基本形状態の記憶値
は、時刻（ｔ＋１）の伝送路応答変換と時刻（ｔ＋２）
のブランチメトリック計算とに利用される。ビタビプロ
セッサ２０７の動作は縮退されたトレリス線図を基にし
ていることを除いて文献１、２の系列推定装置と同様で
あるので、詳細は省略する。
刻での縮退状態に対応する基本形状態を生残りパスの履
歴より調べ、それを非縮退状態記憶回路２０８に出力す
る。例えば、時刻（ｔ＋１）の縮退状態（１１）が前出
の時刻ｔの縮退状態（０１）から遷移した場合、縮退状
態（１１）に対する基本形状態は（０１１）となり、こ
の値が縮退状態（１１）の非縮退状態として非縮退状態
記憶回路２０８に登録される。この基本形状態の記憶値
は、時刻（ｔ＋１）の伝送路応答変換と時刻（ｔ＋２）
のブランチメトリック計算とに利用される。ビタビプロ
セッサ２０７の動作は縮退されたトレリス線図を基にし
ていることを除いて文献１、２の系列推定装置と同様で
あるので、詳細は省略する。
【００５４】以上の実施例では、送信信号が２値、２波
レベルの場合の例について述べたが、２値以外の多値信
号の場合、３波以上の複数の応答がある場合にも本発明
の系列推定装置が有効であることも明らかである。ま
た、以上の実施例では、基となる状態遷移図を１回縮退
した状態遷移図の例を説明したが、本発明の系列推定装
置を基となる状態遷移図を２回以上縮退した状態遷移図
に基づいて動作させることも可能である。
レベルの場合の例について述べたが、２値以外の多値信
号の場合、３波以上の複数の応答がある場合にも本発明
の系列推定装置が有効であることも明らかである。ま
た、以上の実施例では、基となる状態遷移図を１回縮退
した状態遷移図の例を説明したが、本発明の系列推定装
置を基となる状態遷移図を２回以上縮退した状態遷移図
に基づいて動作させることも可能である。
【００５５】
【発明の効果】以上に詳しく述べたように、本発明は、
状態の縮退を利用することにより状態数を削減した遷移
図上でビタビアルゴリズムを動作させることを可能とす
る。したがって、従来方式より少ない記憶量と処理量、
あるいは小さな装置規模で、高速に変動する伝送路に追
従する系列推定装置を提供することができる。
状態の縮退を利用することにより状態数を削減した遷移
図上でビタビアルゴリズムを動作させることを可能とす
る。したがって、従来方式より少ない記憶量と処理量、
あるいは小さな装置規模で、高速に変動する伝送路に追
従する系列推定装置を提供することができる。
【図１】第１の発明に係る系列推定装置の一実施例を示
すブロック図
すブロック図
【図２】第２の発明に係る系列推定装置の一実施例を示
すブロック図
すブロック図
【図３】従来の系列推定装置を示すブロック図
【図４】縮退を行う前の基本となる状態遷移図の例
【図５】第１の発明の系列推定装置の動作の基となる縮
退された状態遷移図の例
退された状態遷移図の例
【図６】第２の発明の系列推定装置の動作の基となる縮
退された状態遷移図の例
退された状態遷移図の例
【図７】従来の系列推定装置のブランチメトリック計算
を説明するためのトレリス線図
を説明するためのトレリス線図
【図８】第１の発明の系列推定装置のブランチメトリッ
ク計算を説明するためのトレリス線図
ク計算を説明するためのトレリス線図
【図９】第２の発明の系列推定装置のブランチメトリッ
ク計算を説明するためのトレリス線図
ク計算を説明するためのトレリス線図
【図１０】１次単純縮退の状態変換図表
【図１１】１次回転縮退の状態変換図表
【図１２】従来の系列推定装置のブランチメトリック計
算を説明するための図表
算を説明するための図表
【図１３】第１の発明の系列推定装置のブランチメトリ
ック計算を説明するための図表
ック計算を説明するための図表
【図１４】第２の発明の系列推定装置のブランチメトリ
ック計算を説明するための図表
ック計算を説明するための図表
【図１５】第２の発明の系列推定装置のブランチメトリ
ック計算を説明するための図表
ック計算を説明するための図表
【図１６】状態変換回路２０９の変換規則を説明するた
めの図表である。
めの図表である。
【図１７】状態変換回路２０９の変換規則を説明するた
めの図表である。
めの図表である。
１０１、２０１、３０１ 入力端子 １０２、２０２、３０２ レジスタ １０３、２０３、３０３ 伝送路応答計算回路 １０５、１０６ マトリススイッチ １０４、２０４、３０４ ブランチメトリック計算回路 １０５ 伝送路応答検査回路 ２０５ 伝送路応答変換回路 １０６、２０６ 伝送路応答記憶回路 １０７、２０７、３０５ ビタビプロセッサ １０８、２０８ 非縮退状態記憶回路 ２０９ 状態変換回路 １０９、２１０、３０６ 出力端子
Claims (2)
- 【請求項１】 受信信号のサンプル値を複数個記憶する
レジスタと、該レジスタから複数個の前記サンプル値を
入力し、送信される可能性のあるＭ（２以上の整数）値
信号系列から指定された複数の信号系列に対してのみ現
時刻の伝送路応答をそれぞれ推定する伝送路応答計算回
路と、該伝送路応答計算回路で求めた前記複数の信号系
列に対する前記現時刻の伝送路応答の値の正当性を検査
し、正当である場合は前記伝送路応答計算回路の出力
を、正当でない場合は同じ信号系列に対する前時刻の伝
送路応答推定値を、現時刻の伝送路応答推定値としてそ
れぞれ出力する伝送路応答検査回路と、該伝送路応答検
査回路が出力する前記現時刻の伝送路応答推定値を記憶
するとともに、前時刻に記憶した伝送路応答推定値を逆
に前記伝送路応答検査回路に供給する伝送路応答記憶回
路と、前記伝送路応答検査回路が出力した前記現時刻の
伝送路応答推定値と各縮退状態に対応する非縮退状態と
現時刻の送信信号候補との３者に基づいて前記複数の系
列のそれぞれの系列に対する仮想受信信号点を求め、前
記受信信号のサンプル値との距離を求めるブランチメト
リック計算回路と、各時刻での前記各縮退状態に対応す
る非縮退状態を出力しながら前記ブランチメトリック計
算回路の出力を受けてビタビアルゴリズムにより受信信
号を判定するビタビプロセッサと、前記非縮退状態を記
憶し、前記伝送路応答計算回路に対して伝送路応答を計
算する前記複数の信号系列を指定する非縮退状態記憶回
路と、から構成され、縮退された状態遷移図に基づいて
送信信号系列の推定を行うことを特徴とする系列推定装
置。 - 【請求項２】 受信信号のサンプル値を複数個記憶する
レジスタと、該レジスタから複数個の前記サンプル値を
入力し、送信される可能性のあるＭ（２以上の整数）値
信号系列から指定された複数の信号系列に対してのみ現
時刻の伝送路応答をそれぞれ推定する伝送路応答計算回
路と、該伝送路応答計算回路で求めた複数の系列に対す
る前記現時刻の伝送路応答の値の正当性を検査し、正当
である場合は前記伝送路応答計算回路の出力を、正当で
ない場合は同じ信号系列に対する前時刻の伝送路応答推
定値を状態遷移図と非縮退状態の記憶とが定める規則に
従って変換した値を、現時刻の伝送路応答推定値として
それぞれ出力する伝送路応答変換回路と、該伝送路応答
変換回路が出力する前記現時刻の伝送路応答推定値を記
憶するとともに前時刻に記憶した伝送路応答推定値を逆
に前記伝送路応答変換回路に供給する伝送路応答記憶回
路と、前記伝送路応答変換回路が出力した前記現時刻の
伝送路応答推定値と各縮退状態に対応する非縮退状態と
現時刻の送信信号候補との３者に基づいて前記複数の系
列のそれぞれの系列に対する仮想受信信号求め、前記受
信信号のサンプル値との距離を求めるブランチメトリッ
ク計算回路と、各時刻での前記各縮退状態に対応する非
縮退状態を出力しながら前記ブランチメトリック計算回
路の出力を受けてビタビアルゴリズムにより最尤系列を
判定し最尤系列の最古状態に対応する非縮退状態の記憶
に従い判定値を出力するビタビプロセッサと、前記非縮
退状態を記憶する非縮退状態記憶回路と、該非縮退状態
記憶回路が与える前記非縮退状態を変換し、変換された
状態を基に前記伝送路応答計算回路に対して伝送路応答
を計算する前記複数の信号系列を指定する状態変換回路
と、から構成され、縮退された状態遷移図に基づいて送
信信号系列の推定を行うことを特徴とする系列推定装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4109961A JP2551296B2 (ja) | 1991-04-30 | 1992-04-28 | 系列推定装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3-128494 | 1991-04-30 | ||
| JP12849491 | 1991-04-30 | ||
| JP4109961A JP2551296B2 (ja) | 1991-04-30 | 1992-04-28 | 系列推定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0730437A JPH0730437A (ja) | 1995-01-31 |
| JP2551296B2 true JP2551296B2 (ja) | 1996-11-06 |
Family
ID=26449670
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4109961A Expired - Lifetime JP2551296B2 (ja) | 1991-04-30 | 1992-04-28 | 系列推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2551296B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5127647B2 (ja) * | 2008-09-24 | 2013-01-23 | 三菱電機株式会社 | 受信装置および復調方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2739318B2 (ja) * | 1988-02-05 | 1998-04-15 | 日本電信電話株式会社 | 最尤受信機 |
| GB2215567B (en) * | 1988-03-05 | 1992-11-18 | Plessey Co Plc | Improvements in or relating to equalisers |
-
1992
- 1992-04-28 JP JP4109961A patent/JP2551296B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0730437A (ja) | 1995-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5325402A (en) | Method and arrangement for estimating data sequences transmsitted using Viterbi algorithm | |
| AU695996B2 (en) | A system for decoding digital data using a variable decision depth | |
| JP4555403B2 (ja) | ディジタル伝送装置の受信機のための状態削減シーケンス推定法によるイコライザ | |
| JPH08256184A (ja) | シンボル間干渉の影響を受けた信号受信のための方法及び装置 | |
| KR101005577B1 (ko) | 기저장된 값을 이용한 메트릭 계산 | |
| JP3180761B2 (ja) | 系列推定方法及び系列推定装置 | |
| US6269116B1 (en) | Method and arrangement for demodulating data symbols | |
| WO1996016483A9 (en) | A system for decoding digital data using a variable decision depth | |
| US5905743A (en) | Apparatus, methods and computer program products for sequential maximum likelihood estimating communications signals using whitening path metrics | |
| JP4904276B2 (ja) | ローカル帰還のある低減状態ビタビ検出器内のパイプライン化判定帰還ユニット | |
| US5450445A (en) | Method and arrangement of estimating data sequences transmitted using viterbi algorithm | |
| US20030115061A1 (en) | MPSK equalizer | |
| US6823027B2 (en) | Method for enhancing soft-value information | |
| US6701483B1 (en) | Fast search-based decoding scheme | |
| JP2551296B2 (ja) | 系列推定装置 | |
| US6782060B2 (en) | Method and apparatus for generating reliability information for channel decoding in a radio receiver | |
| JP4024364B2 (ja) | デジタル伝送システム、デジタル信号の受信機及びデジタル信号の受信方法 | |
| US7054392B2 (en) | Process and device for estimating the successive values of digital symbols, in particular for the equalization of an information transmission channel in mobile telephony | |
| JPH1141115A (ja) | ディジタルデータ復調装置 | |
| Cirpan et al. | Maximum-likelihood estimation of FIR channels excited by convolutionally encoded inputs | |
| RU184852U1 (ru) | Устройство формирования мягкого решения демодулятора сигнала четвертичной фазовой манипуляции со сглаживанием и сдвигом | |
| JP3692644B2 (ja) | 軟判定装置 | |
| JPH0758672A (ja) | ディジタル適応等化器 | |
| JP2560911B2 (ja) | 系列推定装置 | |
| JP2894406B2 (ja) | 最尤系列推定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19960702 |