JP2690196B2 - 最尤系列推定器及び最尤系列推定方法 - Google Patents

Description

【発明の詳細な説明】 技術分野 本発明は、ディジタル通信に使用される受信機等にお
いて、伝送路の歪みを補償する等化器として利用される
最尤系列推定器に関するものである。

背景技術 近年、ディジタル移動通信の開発が急速に行われてい
る。だが、陸上移動通信では、遅延をともなう多数の干
渉波と移動端末が高速に移動することによって、周波数
選択性フェージング（fading）が発生し受信信号波形が
著しく歪むという問題がある。このため等化器によって
この歪みを補償する必要がある。この等化器を構成する
最尤系列推定は、周波数選択性フェージングのように、
伝送路の遅延特性に起因して歪んだ受信信号波形から、
正しい送信データを得るための最も有効な等化方式の一
つである。

先ず、図６及び図７を参照しつつ、ディジタル移動通
信について簡単に説明する。

ディジタル移動通信では、限られた周波数帯域を効率
的に利用するために、また、固定通信網で実施されてい
るISDN（Intergrated Service Digital Network、ディ
ジタル総合ネットワーク）サービスと容易に接続できる
ようにするため、時分割多重方式（Time Division Mult
iple Access;TDMA）が用いられる見込みである。そのTD
MAのフレーム構成例を図６に示す。

図６において、１フレームは、６個のタイムスロット
Slot1〜６から構成されている。また、このタイムスロ
ットの中の一つまたは二つのタイムスロットを一加入者
に割り当てている。各タイムスロットSlot1〜６は、同
期及び等化器のトレーニングのための28ビットのトレー
ニング系列SYNC、12ビットの制御情報の系列SACCH、12
ビットの隣接チャンネル識別用の系列CDVCC、合計260ビ
ットのデータ部分DATA、及び12ビットの未使用領域RSVD
より構成されている。

各加入者に割り当てられたタイムスロットの先頭位置
t slot（ｋ）は、フレーム同期機構によって検出され
る。フレーム同期機構によって検出される検出位置は、
推定による推定先頭位置Et slo（ｋ）（但し、Ｅは推定
の意）である。この推定先頭位置Et slot（ｋ）は、あ
くまでも推定位置なので、必ずしも正確位置ではなく、
フレーム同期機構の補正によって変動する可能性があ
る。

図７はディジタル移動通信の送受信機を示す構成ブロ
ック図である。

この送受信機では、送信機10の出力側に、伝送路20を
介して受信機30が接続されている。

送信機10は、符号化器11、送信ローパスフィルタ（LP
F）12、及び、変調器13等で構成されている。

受信機30は、復調器31、受信ローパスフィルタ（LP
F）32、フレーム同期機構33、サンプリング部34、等化
器35、及び復号器36等で構成されている。

この送受信機の動作を説明する。

送信機10は、符号化器11によって入力データb_mを送信
シンボルｘ＝｛x₁,x₂,…,x_N｝に変換し、送信ローパス
フィルタ12によって撻域制限して送信複素ベースバンド
信号ｓ（ｔ）を生成し、変調器13によって送信複素ベー
スバンド信号ｓ（ｔ）を搬送波で変調信号s_c（ｔ）とし
て伝送路20に送信する。

受信機30は、復調器31によって伝送路20を通った信号
r_c（ｔ）を複素ベースバンド信号ｒ（ｔ）に変換し、受
信ローパスフィルタ32を通して帯域制限された受信複素
ベースバンド信号ｙ（ｔ）を得る。

次に、フレーム同期機構33によって、この信号ｙ
（ｔ）より加入者に割り当てられたタイムスロットの推
定先頭位置（時刻）Et slot（ｋ）を推定し、サンプリ
ング部34によって、信号ｙ（ｔ）を、 ｙ＝〔ｙ（Et slot（ｋ）＋τ_ｌ＋（ｎ−１）Ｔ）｝ （ｎ＝1,2,…,N） ……（１） 但し、Et slot（ｋ）;k個目の割宛アームスロッの推
定先頭時刻（Ｅは推定の意味） τ₁;固定サンプリング位相 T;シンボル間隔時間 でサンプリングする。

そして、等化器35によって、信号ｙ（ｔ）のサンプル
値系列ｙから、周波数選択性フェージングによる伝送路
20の特性を補償し、送信シンボルを推定する。

最後に、復号器36によって、送信シンボル系列の推定
値Ex＝｛Ex_n｝（ｎ＝1,2,…,N、Ｅは推定の意味）を復
号し送信されたデータの推定値Eb_mを得る。

次に、前述した、等化器を構成する最尤系列推定につ
いて説明する。

最尤系列推定については、文献；『J.G.Proakis“Dig
ital Communications"（1983）New York:McGraw−Hill,
p.548〜554.p.610〜648』が開示されている。

前記文献に開示されているように、最尤系列推定は、
ある有限区間での受信信号系列ｙ＝｛y₁,y₂,…,y_N｝が
得られた時に、伝送路のインパルス応答ｈ（ｔ）を既知
として、受信信号系列ｙ＝｛y₁,y₂,…,y_N｝を実現する
確立（尤度）の最も大きい送信シンボル系列ｘ＝｛x₁,x
₂,…,x_N｝を推定するものである。この最尤系列は、伝
送路雑音として白色ガウス雑音を仮定すると、結局、 を最大にするシンボル系列ｘ＝｛x₁,x₂,…,x_N｝を求め
ることによって得られる。

（２）式は、畳み込み符号の復号法として知られるビ
タビ・アルゴリズム（Viterbi algorithm）を用いて効
率的に計算される。

以下、図７の伝送路モデルを示す図８を参照しつつ、
最尤系列推定におけるビタビ・アルゴリズムの原理につ
いて簡単に説明する。

今、伝送路を図８に示すインパルス応答が有限の離散
時間モデルで仮定する。図８において、40はシンボル間
隔Ｔの遅延素子、41は乗算器、42は累積器、43は加算器
を示す。

ｈ＝｛h_j｝（ｊ＝0,…,L）は図７における送・受信ロ
ーパスフィルタ12、32も含めた伝送路20のインパルス応
答ｈ（ｔ）のシンボル間隔Ｔでのサンプル値ｈ（ｔ−j
T）であり、インパルス応答の長さは（Ｌ＋１）Ｔであ
る。またW_nは伝送路の雑音で、加法的な白色ガウス雑音
である。この仮定により（２）式は となる。従って、（３）式においてｋ＝ｎまでの和J
_nは、ｋ＝ｎ−１までの和J_n-1を用いて、 と表される。ここで、J_nはｋ＝１からｎまでの受信信号
系列の対数尤度に比例する量であり、パスメトリック
（経路基準）と呼ばれ、また、（４）式の右辺第２項
は、次に述べる状態遷移についての対数尤度に比例する
量であり、ブランチメトリック（状態遷移基準）と呼ば
れる。

一方、図８に示した伝送路モデルの時刻ｎ−１におけ
る状態は、状態ベクトル S_n-1＝｛x_n-1,x_n-2,…,x_n-L｝ ……（５） によって表される。送信シンボルの種類がＭ個の場合、
伝送路のとり得る状態はM^L個ある。

次に、時刻ｎ−１の状態S_n-1から時刻ｎにおける状態
S_nへの遷移を考えると、M^L個の状態S_nの各々に対して、
Ｍ通りの状態S_n-1からの遷移がある。この時刻間の状態
遷移をブランチと呼び、また、各状態をたどるMⁿ通りの
経路をパスと呼ぶ。各時刻のM^L個のとり得る状態の各々
に対して、一つ前の時刻のＭ個の状態からのパスがあ
る。

ビタビ・アルゴリズムは、各時刻毎に、各状態におけ
るＭ通りの可能性なパスについて（４）式のパスメトリ
ックを計算し、最も値の大きいパスを選択する。従っ
て、各時刻毎に常にM^L通りのパスが残され、過去のパス
は次第に一つに絞られてくる。最終的に一つに収束した
パスより、送信シンボル系列の推定値が得られる。

図９は、従来の最尤系列推定器の構成例を示すブロッ
ク図である。

この最尤系列推定器は、ビタビ・アルゴリズム処理部
50、及び伝送路推定部60で構成されている。

受信信号ｙ（ｔ）は、ｋ個目の割当タイムスロットを
処理する時、サンプリング部で（１）式にサンプリング
され、ビタビ・アルゴリズム処理部50と伝送路推定部60
に入力される。

伝送路推定部60では、サンプル値直列ｙの図６に示す
タイムスロットのトレーニング系列部分においては、サ
ンプル値系列ｙと、これに対応する、図７の受信機30内
の各構成に対し必要な情報を与える図示せぬ上位装置か
ら、与えられる図６のトレーニング系列φを用いて、伝
送路のインパルス応答をLMS（Least Mean Square）等の
適応アルゴリズムにより推定する。サンプル値系列ｙの
うち図６に示すタイムスロットのトレーニング系列部分
以外においては、サンプル値系列ｙと、送信シンボル系
列の推定値Exを用いて、伝送路のインパルス応答の推定
を続ける。

推定された伝送路のインパルス応答Ehはビタビ・アル
ゴリズム処理部50に入力される。

ビタビ・アルゴリズム処理部50では、受信信号のサン
プル値系列ｙと推定された伝送路のインパルス応答Ehを
用いて先に述べた原理に従って送信シンボルの推定を行
なう。

しかしながら、従来の最尤系列推定器では、予め適用
される条件下でのシミュレーション等を行い、ビット誤
り率をみて決定した妥協点にサンプリング位相を固定し
ているため、必ずしも、全ての条件下で、ビット誤り率
が最小となる最適なサンプリング位相となっていないと
いう問題がある。

さらに、最尤系列推定器においては、干渉波の数や遅
延が異なると、ビット誤り率が最小となる最適なサンプ
リング位相が異なり、しかも、サンプリング位相の最適
なサンプリング位相からのづれはビット誤り率の悪化を
引き起こすという問題があり、それらを解決することが
困難であった。

本発明は、前記従来技術が持っていた課題として、サ
ンプリング位相を固定しているために必ずしも全ての条
件下でビット誤り率が最小となる最適なサンプリング位
相になっていないという点と、干渉波の数や遅延が異な
ると、ビット誤り率が最小となる最適なサンプリング位
相が異なり、しかもサンプリング位相の最適なサンプリ
ング位相からのづれがビット誤り率の悪化を引き起こす
という点について解決し、複数の異なるサンプリング位
相の中から、ビット誤り率が最小になると推定されるサ
ンプリング位相を計算し、そのサンプリング位相で等化
処理を行う最尤系列推定器を提供することを目的とす
る。

発明の開示 本発明の、最尤系列を推定する最尤系列推定器は、サ
ンプリング位相計算部、サンプリング処理部、伝送路推
定部及びビタビ・アルゴリズム処理部とを備えている。

ここで、サンプリング位相計算部は、伝送路からの受
信信号ｙ（ｔ）を、タイムスロットの推定先頭時刻Et s
lok（ｋ）から、サンプリング位相の系列τ＝｛τ
（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位相
で、シンボル間隔Ｔにてサンプリングし、Ｍ通りの異な
る前記受信信号ｙ（ｔ）のサンプル値系列 y_i＝｛ｙ（Et slot（ｋ）＋τ（ｉ）＋（ｎ−１）
Ｔ）｝ （ｉ＝1,2,…,M、ｎ＝1,2,…,Nt） を求め、前記サンプル値系列y_iとそのトレーニング系列
φより、トレーニング系列φの複素相互相関Ｃ＝｛C_i｝
（ｉ＝1,2,…,M）を計算し、このＭ個の複素相互相関C_i
の絶対値を比較し、最もその絶対値が大きいサンプリン
グ位相τ（I_k）を選択し、その選択したサンプリング位
相τ（I_k）をもとにした計算より、所定のサンプリング
位相τ_ｓを求める機能を有している。

サンプリング処理部は、所定のタイミングで入力され
る前記受信信号ｙ（ｔ）を推定先頭時刻Et slot（ｋ）
から前記サンプリング位相τ_ｓ及びシンボル間隔Ｔでサ
ンプリングしてサンプル値系列ｙ＝｛y_n｝（ｎ＝1,2,
…,N）を求める機能を有している。

伝送路推定部は前記サンプル値系列ｙとこれに対応す
る前記トレーニング系列φについて及び、送信シンボル
系列の推定値Exと前記サンプル値系列ｙについて、適応
等化アルゴリズムによって前記伝送路のインパルス応答
の推定値Eh＝｛Eh_j｝（ｊ＝0,1,…,L）を求める機能を
有している。

また、ビダビ・アルゴリズム処理部は、前記推定値Eh
に基づきビタビ・アルゴリズムに従って前記サンプル値
系列ｙから前記送信シンボル系列の推定値Exを求めて出
力する機能を有している。

本発明によれば、以上のように最尤系列推定器を構成
したので、サンプリング位相計算部では、複数の異なる
サンプリング位相でサンプリングした受信信号ｙ（ｔ）
のサンプル値系列y_iの各々について、サンプル値系列y_i
とその系列の等化器のトレーニング系列φに対応する部
分とから、等化器のトレーニング系列φの複素相互相関
Ｃを計算する。そして、該複素相互相関C_iが最大となる
サンプリング位相τ（I_k）を選択した後、その選択した
サンプリング位相τ（I_k）をもとにした計算により、ビ
ット誤り率が最小になると推定されるサンプリング位相
τ_ｓを計算し、それをサンプリング処理部へ与える。

サンプリング処理部では、所定のタイミングで入力さ
れる受信信号ｙ（ｔ）をサンプリング位相τ_ｓ及びシン
ボル間隔Ｔでサンプリングしてサンプル値系列ｙを求
め、それを伝送路推定部およびビタビ・アルゴリズム処
理部へ送る。

伝送路推定部は、サンプリング値系列ｙのうちのトレ
ーニング系列の部分では、サンプル値系列ｙとそれに対
応するトレーニング系列φから伝送路のインパルス応答
の推定値Ehを求め、サンプリング値系列ｙのうちのトレ
ーニング系列以外の部分では、受信信号ｙ（ｔ）のサン
プル値系列ｙと送信シンボル系列の推定値Exを用いて伝
送路のインパルス応答の推定値Ehを求め、求めた推定値
Ehをビタビ・アルゴリズム処理部へ与える。

ビダビ・アルゴリズム処理部では、伝送路のインパル
ス応答の推定値Ehに基づき、送信シンボル系列の推定値
Exを求める。

以上により本発明によれば、複数の異なるサンプリン
グ位相の中から、ビット誤り率が最小になると推測され
るサンプリング位相を簡単かつ的確に計算でき、ビット
誤り率が改善された精度のよい送信シンボル系列の推定
が行える。

従って、従来のように、サンプリング位相が固定され
ているため、必ずしも全ての条件下でビット誤り率が最
小となる最適なサンプリング位相になっていないという
問題点や、干渉波の遅延量が異なると、ビット誤り率が
最小となる最適なサンプリング位相が異なり、しかもサ
ンプリング位相の最適なサンプリング位相からのづれが
ビット誤り率の悪化を引き起こすという問題を解決する
ことができる。

図面の簡単な説明 図１は、本発明の第１の実施例である最尤系列推定器
の構成を示す図であり、図２は、本発明のサンプリング
処理部の動作を示すフローチャートであり、図３は、本
発明の伝送路推定部の動作を示すフローチャートであ
り、図４は、本発明のビタビ・アルゴリズム処理部の動
作を示すフローチャートであり、図５は、本発明の第１
の実施例のサンプリング位相計算部の動作を示すフロー
チャートであり、図６は、TDMAのフレーム構成図であ
り、図７は、ディジタル移動通信の送受信機の構成を示
す図であり、図８は、図７の伝送路モデルを示す図であ
り、図９は、従来の最尤系列推定器の構成を示す図であ
り、図10は、図１の遅延量0.25T（10μ sec）時のサン
プリング位相特性図であり、図11は、図１の遅延量1.00
T（40μ sec）時のサンプリング位相特性図であり、図1
2は、本発明の第２の実施例である最尤系列推定器の構
成を示す図であり、図13は、本発明の第２の実施例のサ
ンプリング位相計算部の動作を示すフローチャートであ
り、図14は、本発明の第３の実施例である最尤系列推定
器の構成を示す図であり、図15は、本発明の第３の実施
例のサンプリング位相決算部の動作を示すフローチャー
トであり、図16は本発明の第４の実施例である最尤系列
推定器の構成を示す図であり、図17は、本発明の第４の
実施例のサンプリング位相計算部の動作を示すフローチ
ャートである。

発明を実施するための最良の形態 第１の実施例 第１の実施例では、例えば図７の受信機30内の等化器
35及びサンプリング部34を構成するための最尤系列推定
器について説明する。ここで、単調非増加または単調非
減少なサンプリング位相の系列τ＝｛τ（ｉ）｝（ｉ＝
1,2,…,M）の中から、ビット誤り率が最小になると推測
されるサンプリング位相τ_ｓを計算する場合について、
実施例の説明を行う。

図１は、本発明の実施例である最尤系列推定器の構成
を示す図であり、図１（ａ）は、最尤系列推定器の全体
の機能ブロック図、図１（ｂ）は図１（ａ）中のサンプ
リング位相計算部の機能ブロック図である。

図１（ａ）において、最尤系列推定器は、集積回路等
の個別回路あるいはプロセッサを用いたプログラム制御
などで構成されるものである。また、最尤系列推定器
は、サンプリング位相計算部100、遅延処理部110、サン
プリング処理部120、ビタビ・アルゴリズム処理部130及
び伝送路推定部140を備えている。

サンプリング位相計算部100及び遅延処理部110の出力
側には、サンプリング処理部120が接続されている。サ
ンプリング処理部120の出力側には、ビタビ・アルゴリ
ズム処理部130及び伝送路推定部140が接続されている。

サンプリング位相計算部100は、伝送路からの受信信
号ｙ（ｔ）をサンプリングし、それらのサンプル値系列
y_i｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）と、そのト
レーニング系列SYNCφ＝｛φ_ｍ｝（ｎ＝1,2,…,Nt）に
対する演算を行って所定のサンプリング位相τ_ｓを計算
しサンプリング処理部120に与える機能を有している。

遅延処理部110は、所定時間遅らせて受信信号ｙ
（ｔ）をサンプリング処理部120に与える機能を有して
いる。

サンプリング処理部120は、ｋ個目の割当タイムスロ
ットを処理する時、遅延処理部110からの受信信号ｙ
（ｔ）を、タイムスロットの推定先頭時刻Et slo（ｋ）
から、サンプリング位相τ_ｓ、シンボル間隔Ｔでサンプ
リングしてサンプルｙ＝｛y_n｝（ｎ＝1,2,…,N）の求
め、それをビタビ・アルゴリズム処理部130及び伝送路
推定部140に与える機能を有している。

ビタビ・アルゴリズム処理部130は、サンプル値系列
ｙの入力と、インパルス応答の推定値Eh＝｛Eh_j｝（ｊ
＝0,1,…,L）に基づき、従来から行われているビタビ・
アルゴリズムの原理に従って送信シンボル系列の推定値
Ex＝｛Ex_n｝（ｎ＝1,2,…,N）を求める機能を有し、ブ
ランチメトリック計算手段、及びパスメトリック計算手
段等で構成されている。

伝送路推定部140は、サンプル値系列ｙのうち図６に
示すタイムスロットのトレーニング系列部分において
は、サンプル値系列ｙとこれに対応するトレーニング系
列φとから、伝送路のインパルス応答の推定値Ehを求
め、サンプル値系列ｙのうち図６に示すタイムスロット
のトレーニング系列以外の部分においては、受信信号ｙ
（ｔ）のサンプル値系列ｙと、送信シンボル系列の推定
値Exを用いて、インパルス応答の推定値Ehを求め、更
に、インパルス応答の推定値Ehをビタビ・アルゴリズム
処理部130へ与える機能を有している。

スイッチ手段150は、伝送路推定部140が、サンプル値
系列ｙのうち図６に示すタイムスロットのトレーニング
系列部分を処理する場合にトレーニング系列φ側に切り
替え、トレーニング系列以外の部分を処理する場合に送
信シンボル系列の推定値Ex側に切り替え、それぞれを伝
送路推定部140に与える機能を有している。

図２はサフンプリング処理部の動作を示すフローチャ
ートであり、図３はビタビ・アルゴリズム処理部の動作
を示すフローチャートであり、図４は伝送路推定部の動
作を示すフローチャートである。

ここで、図２、図３、図４及び図１（ａ）を用いて本
発明の最尤系列推定器の動作を説明する。

受信信号ｙ（ｔ）は、割当タイムスロット毎に、サン
プリング位相計算部100と遅延処理部110とに入力され
る。

サンプリング位相計算部100では、ｋ個目の割当タイ
ムスロットを処理する時、受信信号ｙ（ｔ）を、タイム
スロットの推定先頭時刻Et slt（ｋ）から、Ｍ通りの異
なるサンプリング位相τ（ｉ）と、シンボル間隔Ｔでサ
ンプリングする。次いで、それらのサンプル値系列の
内、等化器のトレーニング系列φに対応する部分 ｛ｙ（Et slot（ｋ）＋τ（ｉ）＋（ｎ−１）Ｔ）｝ （ｉ＝1,2,…,M、ｎ＝1,2,…,Nt） と、等化器のトレーニング系列φとの複素相互相関を各
ｉについて計算する。そして、その複素相互相関の絶対
値が最大となるサンプリング位相を選択し、それをもと
に決算したサンプリング位相τ_ｓをサンプリング処理部
120へ与える。

一方、遅延処理部110に入力された受信信号ｙ（ｔ）
は、サンプリング位相計算部100でサンプリング位相τ
_ｓが選択されるまで、遅延されてから、サンプリング処
理部120に入力される。

サンプリング処理部120は、受信信号ｙ（ｔ）及びサ
ンプリング位相τ_ｓを受信し（ST201）、ｋ個目の割当
タイムスロットを処理する時、受信信号ｙ（ｔ）を、 ｙ＝｛ｙ（Et slot（ｋ）＋τ_ｓ＋（ｎ−１）Ｔ）｝ （ｎ＝1,2,…,N） ……（６） 但し、Et slot（ｋ）:k個目の割当タイムスロットの
推定先頭時刻 τ_s;サンプリング位相 T;シンボル間隔時間 にサンプリングし（ST202）、そのサンプル値系列ｙを
ビタビ・アルゴリズム処理部130と伝送路推定部140に与
える（ST203）。

ビタビ・アルゴリズム処理部130では、サンプリング
処理部120からのサンプル値系列ｙと、伝送路推定部140
からの推定値Eh（ST301）とから、先に述べた原理に従
って送信シンボル系列の推定値Exを求め（ST302）出力
する（ST303）。

伝送路推定部140は、受信したサンプル値系列ｙのう
ち図６に示すタイムスロットのトレーニング系列SYNC部
分においては（ST402）、サンプル値系列ｙと、これに
対応する、スイッチ手段150から切り換え入力されるト
レーニング系列φ（ST403）とから、LMS等の適応等化ア
ルゴリズムによって伝送路のインパルス応答の推定値Eh
を求め（ST404）ビタビ・アルゴリズム処理部130へ与え
る（ST407）。更に、伝送路推定部140は、サンプル値系
列ｙのうち図６に示すタイムスロットのトレーニング系
列部分以外においては（ST402）、受信信号ｙ（ｔ）の
サンプル値系列ｙと、スイッチ手段150から切り換え入
力される送信シンボル系列の推定値Ex（S405）を用い
て、伝送路のインパルス応答の推定値Ehを求め（ST40
6）、その推定値Ehをビタビ・アルゴリズム処理部130へ
与える（ST407）。

次に、本発明の特徴である図１（ｂ）のサンプリング
位相計算部100の構成を説明する。

このサンプリング位相計算部100は、受信信号ｙ
（ｔ）を、タイムスロットの推定先頭時刻Et slot
（ｋ）から、前記サンプリング位相の系列τの各々のサ
ンプリング位相で、シンボル間隔Ｔにてサンプリング
し、受信信号ｙ（ｔ）のＭ通りの異なるサンプル値系列
y_iを求めるサンプリング部101−１〜101−Ｍを備えてい
る。サンプリング部101−１〜101−Ｍの各出力側には、
相関値計算器102−１〜102−Ｍを介して最大相関インデ
ックス選択器103が接続されている。

相関値計算器102−１〜102−Ｍは、サンプル値系列y_i
と、そのトレーニング系列φとから複素相互相関C_i（ｉ
＝1,2,…,M）を計算する機能を有している。

最大相関インデックス選択器103は、相関値計算器102
−１〜102−Ｍで計算されたＭ個の複素相互相関C_iの絶
対値を比較し、最もその絶対値が大きいサンプリング位
相τ（I_k）を選択する。そして、選択したサンプリング
位相τ（I_k）のインデックスI_kを、ｋ回目のサンプリン
グ位相選択時に選んだサンプリング位相のインデックス
I_kとして出力する機能を有し、その出力側に、サンプリ
ング位相平滑器104が接続されている。

サンプリング位相平滑器104は、最大相関インデック
ス選択器103で選択したインデックスI_kの区間荷重平均
値、区間最大値及び区間最小値を求め、該区間最大値か
ら該区間最小値を減算して定数倍した値を該区間荷重平
均値に加算したインデックスに対応するサンプリング位
相τ_ｓを出力する機能を有している。

図５は第１の実施例のサンプリング位相計算部の動作
を示すフローチャートである。

ここで、図５及び図１（ｂ）を用いてサンプリング位
相計算部100の動作を説明する。

受信信号ｙ（ｔ）は、割当タイムスロットごとに入力
され、ｋ個目の割当スロットを処理する時、有限区間毎
に、サンプリング部101−１〜101−Ｍで、 y_i＝｛ｙ（Et slot（ｋ）＋τ（ｉ）＋（ｎ−１）
ｔ）｝ （ｉ＝1,2,…,M、ｎ＝1,2,…,Nt） ……（７） にサンプリングされ、相関値計算器102−１〜102−Ｍに
それぞれ入力される（ST501）。

相関値計算器102−１〜102−Ｍでは、受信信号ｙ
（ｔ）のサンプル値系列y_iと、これに対応するトレーニ
ング系列φとの複素相互相関C_iを次式（８）により計算
し、最大相関インデックス選択器103へ与える（ST50
2）。

但し、＊は複素共役を取ることを意味する 最大相関インデックス選択器103では、Ｍ個の複素相
互相関C_iの絶対値が最大となるサンプリング位相τ
（I_k）を選択し、そのインデックスI_kをｋ回目のサンプ
リング位相選択時に選んだサンプリング位相のインデッ
クスI_kとして、サンプリング位相平滑器104に与える（S
T503）。

サンプリング位相平滑器104では、ｋ回目のサンプリ
ング位相計算時に、次式（９）〜（14）によって得られ
るサンプリング位相τ_ｓを出力する（ST504〜ST507）。

但し、w_a（ｊ）（ｊ＝1,2,…,L_s）：荷重係数 L_s;最大計算長 I_s＝［Iavg（ｋ）＋w_s＋w_s（Imax（ｋ）−Imin（ｋ））］ ……（12） 但し、［ ］；ガウス記号（記号内の値を越えない最
大の整数を返す） w_s:シフト係数 I_s;選択サンプリング位相インデックス 図10と図11に、本実施例を適用した場合のビット誤り
率のシミュレーション結果を示す。このシミュレーショ
ンでは、フレーム構成を図６に示す北米標準ディジタル
セルラ方式として、伝送路モデルを２波無相関レイリー
フェージングモデルとしている。

図10は図１の遅延量0.25T（10μ sec）時のサンプリ
ング位相特性図、及び、図11は図１の遅延量1.00T（40
μ sec）の時のサンプリング位相特性図である。この図
10と図11の横軸は、直接波のロールオフ波形の中心から
のサンプリング位相づれを表し、縦軸はビット誤り率BE
Rを表す。この図10と図11より、本実施例の最尤系列推
定器では、従来の最尤系列推定器に比べ、最適な点から
サンプリング位相がづれた場合のビット誤り率BERが小
さくなっていることがわかる。

このように、本実施例では、サンプリング部101−１
〜101−Ｍにより、複数の異なるサンプリング位相τ
（ｉ）の各々について、その位相でサンプリングした場
合の受信信号ｙ（ｔ）のサンプル値系列y_iのトレーニン
グ系列φに対応する部分｛ｙ（Et slot（ｋ）＋τ
（ｉ）＋（ｎ−１）Ｔ）｝（ｎ＝1,2,…,Nt）の、等化
器のトレーニング系列φに対する複素相互相関C_iを相関
値計算器102−１〜102−Ｍで計算する。そして、最大相
関インデックス選択器103及びサンプリング位相平滑器1
04により、複素相互相関C_iが最大となるサンプリング位
相τ（ｉ）から、サンプリング位相τ_ｓを計算するよう
にしている。そのため、従来に比べてビット誤り率（BE
R）を改善できる。

第２の実施例 第２の実施例では、例えば図７の受信機30内の等化器
35及びサンプリング部34を構成するための最尤系列推定
器について説明する。ここで、単調非増加または単調非
減少なサンプリング位相の系列τ＝｛τ（ｉ）｝（ｉ＝
1,2,…,M）の中から、ビット誤り率が最小になると推測
されるサンプリング位相τ_ｓを計算する場合について、
実施例の説明を行う。

図12は、本発明の第２の実施例である最尤系列推定器
の構成を示す図であり、図12（ａ）は、最尤系列推定器
の全体の機能ブロック図、図12（ｂ）は図12（ａ）中の
サンプリング位相計算部の機能ブロック図である。ここ
で、第１の実施例と同一の構成には、同一の符号を付し
ている。

第２の実施例は、第１の実施例に対して、サンプリン
グ位相計算部内のサンプリング位相平滑器の構成のみが
異なる。図12（ｂ）を用いて、第２の実施例のサンプリ
ング位相計算部200内のサンプリング位相平滑器201の構
成を説明する。

サンプリング位相平滑器201は、最大相関インデック
ス選択器103で選択したインデックスI_kの荷重平均値、
最大値及び最小値を求め、該最大値から該最小値を減算
して定数α倍した値を該荷重平均値に加算したインデッ
クスI_kに対応するサンプリング位相τ_ｓを出力する機能
を有している。

第２の実施例を示す図12の他の構成は、第１の実施例
と同一であるため説明は省略する。

図13は本発明の第２の実施例のサンプリング位相計算
部の動作を示すフローチャートである。

次に、図13及び図12（ｂ）を用いてサンプリング位相
計算部200の動作を説明する。

この第２の実施例のサンプリング位相計算部200で
は、第１の実施例と同様に、サンプリング部101−１〜1
01−Ｍ及び相関値計算器102−１〜102−Ｍによって、受
信信号ｙ（ｔ）及びトレーニング系列φから複素相互相
関C_i（ｉ＝1,2,…,M）を求め、最大相関インデックス選
択器103に与えている（ST1301、ST1302）。

そして、最大相関インデックス選択器103がインデッ
クスI_kをサンプリング位相平滑器201へ与えると（S130
3）、該サンプリング位相平滑器201では、ｋ回目のサン
プリング位相選択時に、次式（14）〜（17）と前記した
式（12），（13）とによって得られるサンプリング位相
τ_ｓを出力する（ST1304〜ST1307）。

Imax（ｋ）＝max I_j ……（14） １≦ｊ≦ｋ Imin（ｋ）＝min I_j ……（15） １≦ｊ≦ｋ Iavg（ｋ）＝（１−α）Iavg（ｋ−１）＋αI_k ……（16） 但し、α：定数 Iavg（０）＝０ ……（17） 以上説明した第２実施例では、第１実施例とほぼ同様
の効果が得られるうえに、サンプリング位相平滑器201
での演算量が少ないため、推定処理を高速化することが
できる。

第３の実施例 第３の実施例では、例えば図７の受信機30内の等化器
35及びサンプリング部34を構成するための最尤系列推定
器について説明する。ここで、任意のサンプリング位相
の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の中から、ビ
ット誤り率が最小になると推定されるサンプリング位相
τ_ｓを計算するについて、実施例の説明を行う。

図14は、本発明の第３の実施例である最尤系列推定器
の構成を示す図であり、図14（ａ）は、最尤系列推定器
の全体の機能ブロック図、図14（ｂ）は図14（ａ）中の
サンプリング位相計算部の機能ブロック図である。ここ
で、第１の実施例と同一の構成には、同一の符号を付し
ている。

第３の実施例は、第１の実施例に対してサンプリング
位相計算内の構成が異なる。図14（ｂ）を用いてサンプ
リング位相計算部300の構成を説明する。図14（ｂ）に
おいて、サンプリング位相計算部300は、受信信号ｙ
（ｔ）を、タイムスロットの推定先頭時刻Et slot
（ｔ）から、サンプリング位相の系列τ＝｛τ（ｉ）｝
（ｉ＝1,2,…,M）の各々のサンプリング位相で、シンボ
ル間隔Ｔをシンボル時間にしてサンプリングし、受信信
号ｙ（ｔ）のＭ通りの異なるサンプル値系列ｙい｛_ini
（ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求めるサンプリン
グ部101−１〜101−Ｍを備える。サンプリング部101−
１〜101−Ｍの各出力側には、相関値計算器102−１〜10
2−Ｍを介して最大相関サンプリング位相選択器302が接
続されている。

相関値計算器102−１〜102−Ｍは、サンプル値系列y_i
とトレーニング系列SYNCφ＝｛φ_ｎ｝（ｎ＝1,2,…,N
t）の複素相互相関C_i（ｉ＝1,2,…,M）を計算する機能
を有している。

最大相関サンプリング位相選択器302は、相関値計算
器102−１〜102−Ｍで計算されたＭ個の複素相互相関C_i
の絶対値を比較し、最もその絶対値が大きいサンプリン
グ位相τsel（ｋ）を選択し、出力する機能を有し、そ
の出力側に、サンプリング位相平滑器301が接続されて
いる。

サンプリング位相平滑器301は、フレーム同期機構33
（図７参照）からの推定先頭時刻Et slot（ｋ）に基づ
いて、フレーム同期が確立した時点での受信信号の入力
タイミングを基準にした座標に変換するための座標変換
項tc（ｋ）を求める。そして、サンプリング位相平滑器
301は、最大相関サンプリング位相選択器302で選択した
サンプリング位相τsel（ｋ）に、前記座標変換項tc
（ｋ）を加えてから、その区間荷重平均値、区間最大値
及び区間最小値を求め、該区間最大値から該区間最小値
と前記座標変換項tc（ｋ）の絶対値を減算して定数倍し
た値を、該区間荷重平均値に加算し前記座標変換項tc
（ｋ）を減算したサンプリング位相τ_ｓを出力する機能
を有している。

第３の実施例を示す図14の他の構成は、第１の実施例
と同一であるため説明は省略する。

図15は本発明の第３の実施例のサンプリング位相計算
部の動作を示すフローチャートである。

次に、図15及び図14（ｂ）を用いてサンプリング位相
計算部300の動作を説明する。

受信信号ｙ（ｔ）は、割当タイムスロットごとに入力
され、ｋ個目の割当スロットを処理する時、サンプリン
グ部101−１〜101−Ｍで、式（７）にてサンプリングさ
れ、相関値計算器102−１〜102−Ｍにそれぞれ入力され
る（ST1501）。

相関値計算器102−１〜102−Ｍの各部では、受信信号
ｙ（ｔ）のサンプル値系列y_iと、これに対応するトレー
ニング系列φとの複素相互相関（C_i）を式（８）により
計算し、最大相関サンプリング位相選択器302へ与える
（ST1502）。

最大相関サンプリング位相選択器302では、Ｍ個の相
互相関C_iの絶対値が最導となるサンプリング位相τsel
（ｋ）を選択し、そのサンプリング位相を、ｋ回目のサ
ンプリング位相選択時に選んだサンプリング位相τsel
（ｋ）として、サンプリング位相平滑器301に与える（S
T1504）。

サンプリング位相平滑器301では、ｋ回目のサンプリ
ング位相計算時に、次式（18）〜（23）によって得られ
るサンプリング位相τ_ｓを出力する（ST1504〜ST150
7）。

Δtc（ｋ）＝Et slot（ｋ）−Et slot（ｋ−１） Ｋ＞１ ……（18） τ_ｓ＝τavg（ｋ）−tc（ｋ）＋w_s（τmax（ｋ） −τmin（ｋ）−|tc（ｋ）｜） ……（23） 但し、w_s;シフト係数 w_a（ｊ）（ｊ＝1,2,…,L_s）；荷重係数 L_s;最大計算長 第４の実施例 第４の実施例では、例えば図７の受信機30内の等化器
35及びサンプリング部34を構成するための最尤系列推定
器について説明する。ここで、任意のサンプリング位相
の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の中から、ビ
ット誤り率が最小になると推測されるサンプリング位相
τ_ｓを計算する場合について、実施例の説明を行う。

図16は、本発明の第４の実施例である最尤系列推定器
の構成を示す図であり、図16（ａ）は、最尤系列推定器
の全体の機能ブロック図、図16（ｂ）は、図16（ａ）中
のサンプリング位相計算部の機能ブロック図である。こ
こで、第３の実施例と同一の構成には、同一の符号を付
している。

第４の実施例は、第３の実施例に対してサンプリング
位相計算部内のサンプリング位相平滑器の構成が異な
る。図16（ｂ）を用いて、第４の実施例のサンプリング
位相計算部400内のサンプリング位相平滑器401の構成を
説明する。

サンプリング位相平滑器401は、フレーム同期機構33
（第７参照）からの推定先頭時刻Et slot（ｋ）に基づ
いて、現時刻のフレーム同期における受信信号の入力タ
イミングを基準にした座標に変換するための第１の座標
変換項tc1（ｋ）及び第２の座標変換項tc2（j,k）を求
める。そして、サンプリング位相平滑器401は、最大相
関サンプリング位相選択器302で選択したサンプリング
位相τsel（ｋ）に、前記第２の座標変換項tc2（j,k）
を加えてからその区間荷重平均値、区間最大値及び区間
最小値を求め、該区間最大値から該区間最小値と第１の
座標変換項tc1（ｋ）の絶対値を減算して定数倍した値
を、該区間荷重平均値に加算したサンプリング位相τ_ｓ
を出力する機能を有している。

第４の実施例を示す図16の他の構成は、第３の実施例
と同一であるため説明は省略する。

図17は本発明の第４の実施例のサンプリング位相計算
部の動作を示すフローチャートである。

次に、図17及び図16（ｂ）を用いてサンプリング位相
計算部400の動作を説明する。

第４の実施例のサンプリング位相計算部400では、第
１の実施例と同様に、サンプリング部101−１〜101−Ｍ
及び相関値計算器101−１〜102−Ｍによって、受信信号
ｙ（ｔ）及びトレーニング系列φから複素相互相関Ｃ＝
｛C_i｝（ｉ＝1,2,…,M）を求め最大相関サンプリング位
相選択機302に与えている（ST1701、ST1702）。

そして、最大相関サンプリング位相選択器302が、相
互相関C_iの絶対値が最大となるサンプリング位相τsel
（ｋ）を選択し、そのサンプリング位相を、ｋ回目のサ
ンプリング位相選択時に選んだサンプリング位相τsel
（ｋ）として、サンプリング位相平滑器401に与えると
（ST1703）、該サンプリング位相平滑器401では、ｋ回
目のサンプリング位相計算時に、次式（24）〜（30）に
よって得られるサンプリング位相τ_ｓを出力する（ST17
04〜ST1707）。

τ_ｓ＝τAVG（ｋ）＋w_s（τMAX（ｋ） −τMIN（ｋ）−|tc1（ｋ）｜） ……（30） 但し、w_s;シフト係数 w_a（ｊ）（ｊ＝1,2,…,l_s）；荷重係数 L_s;最大計算長 フレーム同期機構によるタイミング補正によって、各
タイムスロットで選択されるサンプリング位相の基準点
は、タイムスロット毎に異なる。このことを考慮してい
ない第１の実施例による計算では、サンプリング位相τ
_ｓを得るための計算に誤り率が悪化する等の悪影響を与
える場合があるが、第３の実施例及び第４の実施例で
は、τ_ｓを得るための計算への悪影響を防止することが
できる。

なお、本発明は、上記各実施例に限定されず、種々の
変形が可能である。その変形例としては、例えば次の様
なものがある。

（ｉ）最尤系列推定器として、例えば、図1,図12,図14,
図16に示した構成に、整合フィルタや白色化整合フィル
タ等を加え、該白色化フィルタの出力をビタビ・アルゴ
リズム処理部130に入力することも可能である。これら
の構成では、パスメトリックの計算式が異なるだけであ
り、従って、これらの構成の最尤系列推定器にも、上記
実施例を同様に適用できる。

（ii）各実施例で説明したサンプリング位相選択部は、
図１（ｂ），図12（ｂ），図14（ｂ），図16（ｂ）以外
の回路で構成することも可能である。さらに、各実施例
において、図６のTDMAフレーム構成を用いて送信シンボ
ルの推定値を求めるようにしたが、他のデータ系列を用
いて上記実施例を適用することも可能である。

産業上の利用可能性 以上の説明した本発明は、コードレス電話、自動車電
話、携帯電話、ポケットベル、簡易陸上無線電話、テレ
ターミナルシステムなどの陸上移動通信や、船舶電話、
海事衛星通信などの海上移動体通信、又は、航空機公衆
電話などの航空移動体通信に用いるのに適している。

また、各種の公共機関で用いられている無線移動体通
信や、MCA陸上移動無線通信システム、タクシー無線な
どの業務用自衛通信、特定小電力無線局などに用いるの
にも適している。

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁶ 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｈ０４Ｌ 27/38 Ｈ０４Ｌ 27/00 Ｇ

Claims (12)

(57)【特許請求の範囲】
1. 【請求項１】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングし、前記受信信号
   のサンプル値系列y_i＝｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,
   2,…,Nt）を求め出力するサンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
   レーニング系列φの複素相互相関｛C_i｝（ｉ＝1,2,…,
   M）を求め出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め該サンプリ
   ング位相τ（I_k）のインデックスI_kを出力する最大相関
   インデックス選択手段と、 前記インデックスI_kの区間荷重平均値、区間最大値及び
   区間最小値を求め、該区間最大値から該区間最小値を減
   算して定数倍した値を、該区間荷重平均値に加えたイン
   デックスI_sを求め、該インデックスI_sに対応する前記サ
   ンプリング位相τ_ｓを求め出力するサンプリング位相平
   滑手段と、 所定のタイミングで入力される前記受信信号をシンボル
   間隔Ｔ及び前記サンプリング位相τ_ｓでサンプリングし
   サンプル値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）を求め出力
   するサンプリング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列φ以外の部分の処理では該
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
   部と、 前記推定値Ehに基づきビダビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
2. 【請求項２】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングし、前記受信信号
   ｙ（ｔ）のサンプル値系列y_i＝｛y_in｝（ｉ＝1,2,…,
   M、ｎ＝1,2,…,Nt）を求め相関値計算手段に出力するサ
   ンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）を基に該トレ
   ーニング系列φの複素相互関係C_iを （但し、ｉ＝1,2,…,M、＊は複素共役を取ることを意味
   する。）によって求め最大相関インデックス選択手段に
   出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め該サンプリ
   ング位相τ（I_k）のインデックスI_kをサンプリング位相
   平滑手段に出力する最大相関インデックス選択手段と、 ｋ回目のサンプリング位相選択時に前記インデックスI_k
   の区間荷重平均値Iavg（ｋ）、区間最大値Imax（ｋ）及
   び区間最小値Imin（ｋ）を （但し、W_a（ｊ）（ｊ＝1,2,…,L_s）は荷重係数、L_sは
   最大計算長を表す） によって求め、選択サンプリング位相インデックスI_sを I_s＝［Iavg（ｋ）＋W_s（Imax（ｋ）−Imin（ｋ））］ （但し、［ ］は記号内の値を越えない最大の整数を返
   すガウス記号、W_sはシフト係数を表す。） によって求め、前記インデックスI_sに対応するサンプリ
   ング位相τ_ｓを によって求め、該サンプリング位相τ_ｓをサンプリング
   処理部に出力するサンプリング位相平滑手段と、 伝送路から入力される前記受信信号を所定時間遅らせて
   からサンプリング処理部に出力する遅延処理部と、 前記遅延処理部からの前記受信信号をシンボル間隔Ｔ及
   び前記サンプリング位相τ_ｓでサンプリングしサンプル
   値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）を求め伝送路推定部
   及びビタビ・アルゴリズム処理部に出力するサンプリン
   グ処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列φ以外の部分の処理では該
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て推定値Ehを求め、該推定値Ehを出力する伝送路推定部
   と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
3. 【請求項３】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングしてサンプル値系
   列y_i＝｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求
   め、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
   レーニング系列φの複素相互相関｛C_i｝（ｉ＝1,2,…,
   M）を求め、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め、 前記サンプリング位相τ（I_k）のインデックスI_kを求
   め、 前記インデックスI_kの区間荷重平均値、区間最大値及び
   区間最小値を求め、 該区間最大値から該区間最小値を減算して定数倍した値
   を、該区間荷重平均値に加えたインデックスI_sを求め、 前記インデックスI_sに対応する前記サンプリング位相τ
   _ｓを求め、 所定のタイミングで入力される前記受信信号ｙ（ｔ）を
   シンボル間隔Ｔ及び前記サンプリング位相τ_ｓでサンプ
   リングしてサンプル値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）
   を求め、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、 前記サンプル値系列ｙにおける該トレーニング系列のφ
   以外の部分の処理では該サンプル値系列ｙと送信シンボ
   ル系列の推定値Ex＝｛Ex_n｝（ｎ＝1,2,…,N）から適応
   等化アルゴリズムによって該推定値Ehを求め、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ことを特徴とする最尤系列推定方法。
4. 【請求項４】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングし、前記受信信号
   のサンプル値系列y_i＝｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,
   2,…,Nt）を求め出力するサンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
   レーニング系列φの複素相互相関Ｃ＝｛C_i｝（ｉ＝1,2,
   …,M）を求め出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め該サンプリ
   ング位相τ（I_k）のインデックスI_kを出力する最大相関
   インデックス選択手段と、 前記インデックスI_kの荷重平均値、最大値及び最小値を
   求め、該最大値から該最小値を減算して定数倍した値
   を、該荷重平均値に加えたインデックスI_sを求め、該イ
   ンデックスI_sに対応する前記サンプリング位相τ_ｓを求
   め出力するサンプリング位相平滑手段と、 所定のタイミングで入力される前記受信信号をシンボル
   間隔Ｔ及び前記サンプリング位相τ_ｓでサンプリングし
   サンプル値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）を求め出力
   するサンプリング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列のφ以外の部分の処理では
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
   部と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
5. 【請求項５】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングし、前記受信信号
   ｙ（ｔ）のサンプル値系列y_i＝｛y_in｝（ｉ＝1,2,…,
   M、ｎ＝1,2,…,Nt）を求め相関値計算手段に出力するサ
   ンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）を基に該トレ
   ーニング系列φの複素相互相関C_iを （但し、ｉ＝1,2,…,M、＊は複素共役を取ることを意味
   する。）によって求め最大相関インデックス選択手段に
   出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め、該サンプ
   リング位相τ（I_k）のインデックスI_kをサンプリング位
   相平滑手段に出力する最大相関インデックス選択手段
   と、 ｋ回目のサンプリング位相選択時に前記インデックスI_k
   の荷重平均値Iavg（ｋ）、最大値Imax（ｋ）及び最小値
   Imin（ｋ）を Imax（ｋ）＝max I_j １≦ｊ≦ｋ Imin（ｋ）＝min I_j １≦ｊ≦ｋ Iavg（ｋ）＝（１−α）Iavg（ｋ−１）＋αI_k（α：定
   数） Iavg（０）＝０ によって求め、選択サンプリング位相インデックスI_sを I_s＝［Iavg（ｋ）＋W₆（Imax（ｋ）−Imin（ｋ）］ （但し、［ ］は記号内の値を越えない最大の整数を返
   すガウス記号、W_sはシフト係数を表す。） によって求め、前記インデックスI_sに対応するサンプリ
   ング位相τ_ｓを によって求め、該サンプリング位相τ_ｓをサンプリング
   処理部に出力するサンプリング位相平滑手段と、 伝送路から入力される前記受信信号を所定時間遅らせて
   からサンプリング処理部に出力する遅延処理部と、 前記遅延所処理部からの前記受信信号をシンボル間隔Ｔ
   及び前記サンプリング位相τ_ｓでサンプリングしサンプ
   ル値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）を求め伝送路推定
   部及びビタビ・アルゴリズム処理部に出力するサンプリ
   ング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列φ以外の部分の処理では該
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
   部と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
6. 【請求項６】伝送路から入力される受信信号を、単調非
   増加または単調非減少なサンプリング位相の系列τ＝
   ｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサンプリング位
   相でシンボル間隔Ｔにてサンプリングしてサンプル値系
   列y_i＝｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求
   め、 前記サンプリング値系列y_iと該サンプル値系列y_iのトレ
   ーニング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって
   該トレーニング系列φの複素相互相関｛C_i｝（ｉ＝1,
   …,M）を求め、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τ（I_k）を求め、 前記サンプリング位相τ（I_k）のインデックスI_kを求
   め、 前記インデックスI_kの荷重平均値、最大値及び最小値を
   求め、 該最大値から該最小値を減算して定数倍した値を該荷重
   平均値に加えたインデックスI_sを求め、 前記インデックスI_sに対応する前記サンプリング位相τ
   _ｓを求め、 所定のタイミングで入力される前記受信信号ｙ（ｔ）を
   シンボル間隔Ｔ及び前記サンプリング位相τ_ｓでサンプ
   リングしてサンプル値系列ｙ＝｛y_n｝（ｎ＝1,2,…,N）
   を求め、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、 前記サンプル値系列ｙにおける該トレーニング系列φ以
   外の部分の処理では該サンプル値系列ｙと送信シンボル
   系列の推定値Ex＝｛Ex_n｝（ｎ＝1,2,…,N）から適応等
   化アルゴリズムによって該推定値Ehを求め、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ことを特徴とする際尤系列推定方法。
7. 【請求項７】伝送路から入力される受信信号を、サンプ
   リング位相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の
   各々のサンプリング位相でフレーム同期機構によって定
   められる先頭推定時刻からシンボル間隔Ｔにてサンプリ
   ングし、前記受信信号のサンプル値系列y_i＝｛y_in｝
   （ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求め出力するサン
   プリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
   レーニング系列φの複素相互相関｛C_i｝（ｉ＝1,…,M）
   を求め出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τsel（ｋ）を求め出力す
   る最大相関サンプリング位相選択手段と、 前記先頭推定時刻を基に座標変換項を求め、前記サンプ
   リング位相τsel（ｋ）に該座標変換項を加えてから該
   記サンプリング位相τsel（ｋ）の区間荷重平均値、区
   間最大値及び区間最小値を求め、該区間最大値から前記
   座標変換項の絶対値と該区間最小値を減算して定数倍し
   た値を、該区間荷重平均値に加算し前記座標変換項を減
   算したサンプリング位相τ_ｓを求め出力するサンプリン
   グ位相平滑手段と、 所定のタイミングで入力される前記受信信号を前記先頭
   推定時刻からシンボル間隔Ｔ及び前記サンプリング位相
   τ_ｓでサンプリングしサンプル値系列ｙ＝｛y_n｝（ｎ＝
   1,2,…,N）を求め出力するサンプリング出力部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列φ以外の部分の処理では該
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
   部と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
8. 【請求項８】ｋ個目の割当タイムスロットを処理する場
   合、伝送路から入力される受信信号を、サンプリング位
   相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサ
   ンプリング位相でフレーム同期機構によって定められる
   先頭推定時刻Et slot（ｋ）からシンボル間隔Ｔにてサ
   ンプリングし、前記受信信号ｙ（ｔ）のサンプル値系列
   y_i＝｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,2,…,Mt）を求め相
   関値計算手段に出力するサンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｉ｝（ｎ＝1,2,…,Nt）を基に該トレ
   ーニング系列φの複素相互相関C_iを （但し、ｉ＝1,2,…,M、＊は複素共役を取ることを意味
   する。）によって求め最大相関インデックス選択手段に
   出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τsel（ｋ）を求め該サン
   プリング位相τsel（ｋ）をサンプリング位相平滑手段
   に出力する最大相関サンプリング位相選択手段と、 座標変換項tc（ｋ）を Δtc（ｋ）＝Et slot（ｋ）−Et slot（ｋ−１）ｋ＞１ にて求め、前記サンプリング位相τsel（ｋ）の区間荷
   重平均値τavg（ｋ）、区間最大値τmax（ｋ）及び区間
   最小値τmin（ｋ）を （但し、W_a（ｊ）（ｊ＝1,2,…,L_s）は荷重係数、L_sは
   最大計算長を表す） によって求め、選択サンプリング位相τ_ｓを τ_ｓ＝τavg（ｋ）−tc（ｋ） ＋W_s（τmax（ｋ）−τmin（ｋ）−|tc（ｋ）｜） （但し、W_sはシフト係数、W_a（ｊ）（ｊ＝1,2,…,L_s）
   は荷重係数、L_sは最大計算長を表す。） によって求め、該サンプリング位相τ_ｓをサンプリング
   処理部に出力するサンプリング位相平滑手段と、 伝送路から入力される前記受信信号を所定時間遅らせて
   からサンプリング処理部に出力する遅延処理部と、 前記遅延処理部からの前記受信信号を前記先頭推定時刻
   et slot（ｋ）からシンボル間隔Ｔ及び前記サンプリン
   グ位相τ_ｓでサンプリングしサンプル値系列ｙ＝｛y_n｝
   （ｎ＝1,2,…,N）を求め伝送推定部及びビタビ・アルゴ
   リズム処理部に出力するサンプリング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列に対応する該トレーニング系列φから適応等化アルゴ
   リズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
   における該トレーニング系列φ以外の部分の処理では該
   サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
   _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
   て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
   部と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ビタビ・アルゴリズム処理部と、から構成したことを特
   徴とする最尤系列推定器。
9. 【請求項９】伝送路から入力される受信信号を、サンプ
   リング位相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の
   各々のサンプリング位相でフレーム同期機構によって定
   められる先頭推定時刻からシンボル間隔Ｔにてサンプリ
   ングしてサンプル値系列y_i＝｛y_in｝（ｉ＝1,2,…,M、
   ｎ＝1,2,…,Nt）を求め、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
   ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
   レーニング系列φの複素相互相関｛C_i｝（ｉ＝1,…,M）
   を求め、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
   最大となるサンプリング位相τsel（ｋ）を求め、 前記先頭推定時刻を基に座標変換項を求め、 前記サンプリング位相τsel（ｋ）に該座標変換項を加
   えてから該記サンプリング位相τsel（ｋ）の区間荷重
   平均値、区間最大値及び区間最小値を求め、 該区間最大値から前記座標変換項の絶対値と該区間最小
   値を減算して定数倍した値を、該区間荷重平均値に加算
   し前記座標変換項を減算したサンプリング位相τ_ｓを求
   め、 所定のタイミングで入力される前記受信信号ｙ（ｔ）を
   前記先頭推定時刻からシンボル間隔Ｔ及び前記サンプリ
   ング位相τ_ｓでサンプリングしてサンプル値系列ｙ＝
   ｛y_n｝（ｎ＝1,2,…,N）を求め、 前記サンプル値系列ｙにおける前記トレーニング系列φ
   の部分の処理では該サンプル値系列ｙと該サンプル値系
   列ｙに対応する該トレーニング系列φから適応等化アル
   ゴリズムによって伝送路のインパルス応答の推定値Eh＝
   ｛Eh_j｝（ｊ＝0,1,…,L）を求め、 前記サンプル値系列ｙにおける該トレーニング系列φ以
   外の部分の原理では該サンプル値系列ｙと送信シンボル
   系列の推定値Ex＝｛Ex_n｝（ｎ＝1,2,…,N）から適応等
   化アルゴリズムによって該推定値Ehを求め、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
   記サンプル値系列ｙから前記推定値Exを求めて出力する
   ことを特徴とする最尤系列推定方法。
10. 【請求項１０】伝送路から入力される受信信号を、フレ
    ーム同期機構によって定められる先頭推定時刻からサン
    プリング位相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）
    の各々のサンプリング位相でフレーム同期機構によって
    定められる先頭推定時刻からシンボル間隔Ｔにてサンプ
    リングし、前記受信信号のサンプル値系列y_i＝｛y_in｝
    （ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求め出力するサン
    プリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
    ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
    レーニング系列φの複素相互相関｛C_i｝（ｉ＝1,…,M）
    を求め出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
    最大となるサンプリング位相τsel（ｋ）を求め出力す
    る最大相関サンプリング位相選択手段と、 前記先頭推定時刻を基に第１の座標変換項及び第２の座
    標変換項を求め、前記サンプリング位相τsel（ｋ）に
    該第２の座標変換項を加えてから該記サンプリング位相
    τsel（ｋ）の区間荷重平均値、区間最大値及び区間最
    小値を求め、該区間最大値から前記第１の座標変換項の
    絶対値と該区間最小値を減算して定数倍した値を、該区
    間荷重平均値に加算したサンプリング位相τ_ｓを求め出
    力するサンプリング位相平滑手段と、 所定のタイミングで入力される前記受信信号を前記先頭
    推定時刻からシンボル間隔Ｔ及び前記サンプリング位相
    τ_ｓでサンプリングしサンプル値系列ｙ＝｛y_n｝（ｎ＝
    1,2,…,N）を求め出力するサンプリング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
    の部分の処理では該サンプル値系列ｙと該サンプル値系
    列ｙに対応する該トレーニング系列φから適応等化アル
    ゴリズムによって伝送路のインパルス応答の推定値Eh＝
    ｛Eh_j｝（ｉ＝0,1,…,L）を求め、該サンプル値系列ｙ
    における該トレーニング系列φ以外の部分の処理では該
    サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
    _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
    て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
    部と、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
    記サンプル値系列ｙから前記推定値Exを求めて出力する
    ビタビ・アルゴリズム処理部と、から構成したことを特
    徴とする最尤系列推定器。
11. 【請求項１１】ｋ個目の割当タイムスロットを処理する
    場合伝送路から入力される受信信号を、サンプリング位
    相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）の各々のサ
    ンプリング位相でフレーム同期機構によって定められる
    先頭推定時刻Et slot（ｋ）からシンボル間隔Ｔにてサ
    ンプリングし、前記受信信号ｙ（ｔ）のサンプル値系列
    y_i｛y_in｝（ｉ＝1,2,…,M、ｎ＝1,2,…,Nt）を求め相関
    値計算手段に出力するサンプリング手段と、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
    ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）を基に該トレ
    ーニング系列φの複素相互相関C_iを （但し、ｉ＝1,2,…,M、＊は複素共役を取ることを意味
    する。）によって求め最大相関インデックス選択手段に
    出力する相関値計算手段と、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
    最大となるサンプリング位相τsel（ｋ）を求め該サン
    プリング位相τsel（ｋ）をサンプリング位相平滑手段
    に出力する最大相関サンプリング位相選択手段と、 ｋ回目のサンプリング位相選択時に第１の座標変換項tc
    l（ｋ）及び第２の座標変換項tc2（j,k）を にて求め、前記サンプリング位相τ（I_k）の区間荷重平
    均値τAVG（ｋ）、区間最大値τMAX（ｋ）及び区間最小
    値τMIN（ｋ）を （但し、W_a（ｊ）（ｊ＝1,2,…,L_s）は荷重係数、L_sは
    最大計算長を表す） によって求め、選択サンプリング位相τ_ｓを τ_ｓ＝τAVG（ｋ）＋W_s（τMAX（ｋ）−τMIN（ｋ）−|
    tc1（ｋ）｜） （但し、W_sはシフト係数、W_a（ｊ）（ｊ＝1,2,…,L_s）
    は荷重係数、L_sは最大計算長を表す。） によって求め、該サンプリング位相τ_ｓをサンプリング
    処理部に出力するサンプリング位相平滑手段と、 伝送路から入力される前記受信信号を所定時間遅らせて
    からサンプリング処理部に出力する遅延処理部と、 前記遅延処理部からの前記受信信号を前記先頭推定時刻
    Et slot（ｋ）からシンボル間隔Ｔ及び前記サンプリン
    グ位相τ_ｓでサンプリングしサンプル値系列ｙ＝｛y_n｝
    （ｎ＝1,2,…,N）を求め伝送路推定部及びビタビ・アル
    ゴリズム処理部に出力するサンプリング処理部と、 前記サンプル値系列ｙにおける前記トレーニング系列φ
    の部分の処理では該サンプル値系列ｙと該サンプル値系
    列ｙに対応する該トレーニング系列φから適応等化アル
    ゴリズムによって伝送路のインパルス応答の推定値Eh＝
    ｛Eh_j｝（ｊ＝0,1,…,L）を求め、該サンプル値系列ｙ
    における該トレーニング系列φ以外の部分の処理では該
    サンプル値系列ｙと送信シンボル系列の推定値Ex＝｛Ex
    _n｝（ｎ＝1,2,…,N）から適応等化アルゴリズムによっ
    て該推定値Ehを求め、該推定値Ehを出力する伝送路推定
    部と、 前記推定値Ehに求づきビタビ・アルゴリズムに従って前
    記サンプル値系列ｙから前記推定値Exを求めて出力する
    ビタビ・アルゴリズム処理部と、から構成したことを特
    徴とする最尤系列推定器。
12. 【請求項１２】伝送路から入力される受信信号を、サン
    プリング位相の系列τ＝｛τ（ｉ）｝（ｉ＝1,2,…,M）
    の各々のサンプリング位相でフレーム同期機構によって
    定められる先頭推定時刻からシンボル間隔Ｔにてサンプ
    リングしてサンプル値系列y_i｛y_in｝（ｉ＝1,2,…,M、
    ｎ＝1,2,…,Nt）を求め、 前記サンプル値系列y_iと該サンプル値系列y_iのトレーニ
    ング系列φ＝｛φ_ｎ｝（ｎ＝1,2,…,Nt）によって該ト
    レーニング系列φの複素相互相関Ｃ＝｛C_i｝（ｉ＝1,
    …,M）を求め、 前記Ｍ個の複素相互相関C_iの絶対値を比較し該絶対値が
    最大となるサンプリング位相τsel（ｋ）を求め、 前記先頭推定時刻を基に第１の座標変換項及び第２の座
    標変換項を求め、 前記サンプリング位相τsel（ｋ）に該第２の座標変換
    項を加えてから該記サンプリング位相τsel（ｋ）の区
    間荷重平均値、区間最大値及び区間最小値を求め、 該区間最大値から前記第１の座標変換項の絶対値と該区
    間最小値を減算して定数倍した値を、該区間荷重平均値
    に加算したサンプリング位相τ_ｓを求め、 所定のタイミングで入力される前記受信信号ｙ（ｔ）を
    前記先頭推定時刻からシンボル間隔Ｔ及び前記サンプリ
    ング位相τ_ｓでサンプリングしてサンプル値系列ｙ＝
    ｛y_n｝（ｎ＝1,2,…,N）を求め、 前記サンプル値系列ｙにおける前記トレーニング系列φ
    の部分の処理では該サンプル値系列ｙと該サンプル値系
    列ｙに対応する該トレーニング系列φから適応化等アル
    ゴリズムによって伝送路のインパルス応答の推定値Eh＝
    ｛Eh_j｝（ｊ＝0,1,…,L）を求め、 前記サンプル値系列ｙにおける該トレーニング系列φ以
    外の部分の処理では該サンプル値系列ｙと送信シンボル
    系列の推定値Ex＝｛Ex_n｝（ｎ＝1,2,…,N）から適応等
    化アルゴリズムによって該推定値Ehを求め、 前記推定値Ehに基づきビタビ・アルゴリズムに従って前
    記サンプル値系列ｙから前記推定値Exを求めて出力する
    ことを特徴とする最尤系列推定方法。