(57)【要約】  無線通信システムにおけるチャネル推定のための方法が、以下のステップを含む。チャネル伝達関数が、連続パイロット・セルと散乱パイロット・セルにおける送信信号と受信信号を使用して、連続パイロット・セルと散乱パイロット・セルにおいて計算される。連続パイロット・セルと散乱パイロット・セルにおいて計算されたチャネル伝達関数を使用して時間ドメイン適応補間が実行されて、散乱パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数が得られる。連続パイロット・セルと散乱パイロット・セルにおいて計算されたチャネル伝達関数を使用して周波数ドメイン適応補間が実行されて、非パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数が得られる。
2008543186.tif 000002
【特許請求の範囲】
【請求項１】
無線通信システムにおけるチャネル推定のための方法であって、
（Ｉ）連続パイロット・セルと散乱パイロット・セルにおけるチャネル伝達関数を、前記連続パイロット・セルと前記散乱パイロット・セルにおける、送信信号と受信信号を使用して計算するステップと、
（II）ステップ（Ｉ）からの前記計算されたチャネル伝達関数を使用して時間ドメイン適応補間を実行して、散乱パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数を得るステップと、
（III）ステップ（Ｉ）からの前記計算されたチャネル伝達関数を使用して周波数ドメイン適応補間を実行して、非パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数を得るステップとを含む方法。
【請求項２】
ステップ（II）は、
（Ａ）（１）現在の補間フィルタ係数を使用して、連続パイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数をステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、連続パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップと
を含む補間フィルタ係数を更新するステップ、
（Ｂ）前記更新された補間フィルタ係数を使用して、散乱パイロット・トーンにおいて補間を実行して、前記散乱パイロット・トーンの前記非パイロット・セルにおける伝達関数推定値を得るステップをさらに含む請求項１に記載の方法。
【請求項３】
ステップ（III）は、
（Ａ）（１）現在の補間フィルタ係数を使用して、非パイロット・トーンの非パイロット・セルにおいて補間を実行するステップと、
（２）ステップ（１）からの補間結果を使用して、送信信号に対するチャネルひずみを補正するステップと、
（３）ステップ（２）からの前記補正された受信信号に基づく判定を介して、前記送信信号を推定するステップと、
（４）前記判定の前後に差を計算して、前記差の大きさが、事前設定された閾値より大きい場合、ステップ（５）〜（７）をスキップするステップと、
（５）受信信号と、前記送信信号の前記判定とに基づき、非パイロット・トーンにおけるチャネル伝達関数を計算するステップと、
（６）ステップ（５）からの前記計算されたチャネル伝達関数を、ステップ（１）からの前記補間結果と比較することにより、推定誤差を計算するステップと、
（７）ステップ（６）からの前記計算された推定結果を使用して、非パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップとを含む補間フィルタ係数を更新するステップをさらに含む請求項１に記載の方法。
【請求項４】
ステップ（Ａ）は、ＬＭＳアルゴリズムに基づいて実行される請求項２に記載の方法。
【請求項５】
ステップ（５）〜（７）は、ＬＭＳアルゴリズムに基づいて実行される請求項３に記載の方法。
【請求項６】
前記補間フィルタ係数は、複素非対称係数である請求項２に記載の方法。
【請求項７】
前記補間フィルタ係数は、複素非対称係数である請求項３に記載の方法。
【請求項８】
ステップ（III）は、
（Ａ）（１）現在の補間フィルタ係数を使用して、非パイロット・トーンにおいて補間を実行するステップと、
（２）ステップ（１）からの補間結果を使用して、送信信号に対するチャネルひずみを補正するステップと、
（３）ステップ（２）からの前記補正された受信信号に基づく判定を介して、前記送信信号を推定するステップと、
（４）判定誤差を計算するステップと、
（５）前記判定誤差の大きさが、事前設定された閾値より小さい場合、
（ａ）受信信号と、前記送信信号の前記判定とに基づき、非パイロット・トーンにおけるチャネル伝達関数を計算するステップと、
（ｂ）ステップ（ａ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの前記補間結果と比較することにより、推定誤差を計算するステップと、
（ｃ）ステップ（ｂ）からの前記計算された推定誤差を使用して、非パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップとを含む、前記補間フィルタ係数を更新するステップとを含む補間フィルタ係数を更新するステップをさらに含む請求項１に記載の方法。
【請求項９】
ステップ（II）は、
（Ａ）（１）先行するシンボルからの補間フィルタ係数を使用して、所定の数の連続パイロット・トーンの第１のパイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、前記所定の数の連続パイロット・トーンの前記第１のパイロット・トーンにおいて補間フィルタ係数を更新ステップと、
（４）前記所定の数の連続パイロット・トーンの残りの連続パイロット・トーンのそれぞれに関して、ステップ（１）〜（３）を繰り返すステップと
を含む前記所定の数の連続パイロット・トーンにおいて補間フィルタ係数を更新するステップ、及び
（Ｂ）前記更新された補間フィルタ係数を使用して、散乱パイロット・トーンにおいて補間を実行して、前記散乱パイロット・トーンの前記非パイロット・セルにおける伝達関数推定値を得るステップをさらに含む請求項１に記載の方法。
【請求項１０】
ステップ（II）は、
（Ａ）ＬＭＳアルゴリズムに基づいて補間フィルタ係数を更新するステップ、及び
（Ｂ）前記更新された補間フィルタ係数を使用して、散乱パイロット・トーンにおいて補間を実行して、前記散乱パイロット・トーンの前記非パイロット・セルにおける伝達関数推定値を得るステップをさらに含む請求項１に記載の方法。
【請求項１１】
ステップ（Ａ）は、
（１）現在の補間フィルタ係数を使用して、連続パイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、連続パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップとをさらに含む請求項１０に記載の方法。
【請求項１２】
ステップ（Ａ）は、
（１）先行するシンボルからの更新された補間フィルタ係数を使用して、所定の数の連続パイロット・トーンの第１のパイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、前記所定の数の連続パイロット・トーンの前記第１のパイロット・トーンにおける補間フィルタ係数を更新するステップと、
（４）前記所定の数の連続パイロット・トーンの残りの連続パイロット・トーンのそれぞれに関して、ステップ（１）〜（３）を繰り返すステップをさらに含む請求項１０に記載の方法。
【請求項１３】
ステップ（III）は、
（Ａ）ＬＭＳアルゴリズムに基づいて補間フィルタ係数を更新するステップ、及び
（Ｂ）前記更新された補間フィルタ係数を使用して、非パイロット・トーンにおいて補間を実行するステップをさらに含む請求項１に記載の方法。
【請求項１４】
ステップ（Ａ）は、
（１）現在の補間フィルタ係数を使用して、連続パイロット・シンボルにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、連続パイロット・シンボルの少なくともサブセットにおける補間フィルタ係数を更新するステップとをさらに含む請求項１３に記載の方法。
【請求項１５】
ステップ（Ａ）は、
（１）先行するシンボルからの更新された補間フィルタ係数を使用して、所定の数の連続パイロット・トーンの第１のパイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、前記所定の数の連続パイロット・トーンの前記第１のパイロット・トーンにおける補間フィルタ係数を更新するステップと、
（４）前記所定の数の連続パイロット・トーンの残りの連続パイロット・トーンのそれぞれに関して、ステップ（１）〜（３）を繰り返すステップとをさらに含む請求項１３に記載の方法。
【請求項１６】
無線通信システムにおけるチャネル推定のための方法であって、
（Ｉ）連続パイロット・セルと散乱パイロット・セルにおけるチャネル伝達関数を、前記連続パイロット・セル及び前記散乱パイロット・セルにおける送信信号と受信信号を使用して計算するステップ、
（II）（Ａ）（１）現在の補間フィルタ係数を使用して、連続パイロット・トーンにおいて補間を実行するステップと、
（２）ステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、連続パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップと
を含む補間フィルタ係数を更新するステップ、
（Ｂ）前記更新された補間フィルタ係数を使用して、散乱パイロット・トーンにおいて補間を実行して、前記散乱パイロット・トーンの前記非パイロット・セルにおける伝達関数推定値を得るステップとを含む、時間ドメイン適応補間を使用して、散乱パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数を計算するステップ、及び
（III）（Ａ）（１）現在の補間フィルタ係数を使用して、非パイロット・トーンにおいて補間を実行するステップと、
（２）ステップ（III）−（Ａ）−（１）からの補間結果を使用して、送信信号に対するチャネルひずみを補正するステップと、
（３）ステップ（III）−（Ａ）−（２）からの前記補正された受信信号に基づく判定を介して、前記送信信号を推定するステップと、
（４）判定誤差を計算するステップと、
（５）前記判定誤差の大きさが、事前設定された閾値より小さい場合、
（ａ）受信信号と、前記送信信号の前記判定とに基づき、非パイロット・トーンにおけるチャネル伝達関数を計算するステップと、
（ｂ）ステップ（ａ）からの前記計算されたチャネル伝達関数を、ステップ（III）−（Ａ）−（１）からの前記補間結果と比較することにより、推定誤差を計算するステップと、
（ｃ）ステップ（ｂ）からの前記計算された推定誤差を使用して、非パイロット・トーンの少なくともサブセットにおける補間フィルタ係数を更新するステップとを含む、前記補間フィルタ係数を更新するステップと
を含む前記補間フィルタ係数を更新するステップを含む、周波数ドメイン適応補間を使用して、非パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数を計算するステップを含む方法。
【請求項１７】
無線通信システムにおけるチャネル推定のための方法であって、
（Ｉ）連続パイロット・セルと散乱パイロット・セルにおけるチャネル伝達関数を、前記連続パイロット・セル及び前記散乱パイロット・セルにおける送信信号と受信信号を使用して計算するステップ、
（II）周波数ドメイン適応補間を使用して、非パイロット・セルにおけるチャネル伝達関数を計算するステップ、及び
（III）時間ドメイン適応補間を使用して、非パイロット・シンボルの非パイロット・セルにおけるチャネル伝達関数を計算するステップを含む方法。
【請求項１８】
ステップ（II）は、
（Ａ）ＬＭＳアルゴリズムに基づいて補間フィルタ係数を更新するステップ、及び
（Ｂ）前記更新された補間フィルタ係数を使用して、非パイロット・トーンにおいて補間を実行するステップをさらに含む請求項１７に記載の方法。
【請求項１９】
ステップ（Ａ）は、
（１）現在の補間フィルタ係数を使用して、連続パイロット・シンボルにおいて補間を実行するステップと、
（２）請求項１３におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、連続パイロット・シンボルの少なくともサブセットにおける補間フィルタ係数を更新するステップとをさらに含む請求項１８に記載の方法。
【請求項２０】
ステップ（Ａ）は、
（１）先行するシンボルからの更新された補間フィルタ係数を使用して、所定の数の連続パイロット・トーンの第１のパイロット・トーンにおいて補間を実行するステップと、
（２）請求項１におけるステップ（Ｉ）からの前記計算されたチャネル伝達関数を、ステップ（１）からの補間結果と比較することにより、推定誤差を計算するステップと、
（３）前記計算された推定誤差を使用して、前記所定の数の連続パイロット・トーンの前記第１のパイロット・トーンにおける補間フィルタ係数を更新するステップと、
（４）前記所定の数の連続パイロット・トーンの残りの連続パイロット・トーンのそれぞれに関して、ステップ（１）〜（３）を繰り返すステップとをさらに含む請求項１８に記載の方法。
【発明の詳細な説明】
【関連出願】
【０００１】
（関連出願の相互参照）
本出願は、引用によって全体が組み込まれている、２００５年５月２７日に出願した米国仮出願第６０／６８５，７０４号の利益を主張する。
【技術分野】
【０００２】
本発明は、無線通信システムに関し、より詳細には、ＯＦＤＭ通信システムのための改良されたチャネル推定技術に関する。
【背景技術】
【０００３】
通信システムにおいて、情報を運ぶ信号は、信号ひずみを生じさせる通信チャネルを介して、送信元から宛先に伝送される。受信機において、通信チャネルによって生じさせられた信号ひずみは、送信元からの伝送された信号を正確に回復させるように、適切に補正されなければならない。チャネルひずみを補正する通常の例が、等化器である。等化器は、通常、いくつかのトレーニング信号に基づき、何らかの最適な設定にトレーニングされる。この種の適応等化器は、変化のないチャネル、又はゆっくり変化するチャネルに関してうまく働く。無線通信などの時間につれ急速に変化するチャネルの場合、チャネルの変化は、しばしば、非常に速い。その結果、等化器は、急速に変化するチャネル特性に追従するために、非常に頻繁にトレーニングされなければならない。
【０００４】
ＯＦＤＭ（直交周波数分割多重）が、ブロードバンド通信において好評を得ている。ＯＦＤＭシステムにおいて、データ信号は、多くの等間隔の互いに直交な副搬送波に分散される。ＯＦＤＭ変調は、通常、送信機におけるＩＦＦＴ（逆高速フーリエ変換）や、受信機におけるＦＦＴ（高速フーリエ変換）を介して実施される。図１は、ＯＦＤＭベースの無線受信機のブロック図である。
【０００５】
図１において、無線周波数信号が、アンテナを介してＲＦチューナ１００によって受信される。所望される信号が、チューナ１００によって選択され、ダウンコンバータ／フィルタ・ブロック１１０を介してダウンコンバートされて、濾波される。ブロック１１０の出力は、アナログ・ベースバンド（又は元の無線周波数より相当に低い周波数の通過帯域）信号であり、この信号は、Ａ−Ｄブロック１２０によってデジタル信号に変換される。デジタル信号はシンボルにグループ化され、そのシンボル境界がシンボル同期ブロック１３０において適切に識別され、ガード区間（通常、サイクリック・プレフィックス）がブロック１４０において除去されてから、ＦＦＴブロック１５０に供給される。ＦＦＴ１５０の後、信号は周波数ドメインにあり、等化器ＦＥＱ１７０がその周波数ドメインにおいて使用される。直交性を有して、副搬送波は互いに干渉せず、したがって、ＦＥＱ（周波数ドメイン等化器）は、各副搬送波（ときとして、ビン又は搬送波とも呼ばれる）に対して別々に実装される。シンボルは、いくつかのガード区間（サイクリック・プレフィックス）によって分離されるので、ＩＳＩ（シンボル間干渉）が回避される。したがって、そのようなＦＥＱは単純に１タップ複素スケーリングとなる。この複素タップ係数はトレーニングによって適応的に決定されることが可能であり、データ伝送中に更新される。無線通信チャネルなどの、急速に変化するチャネルの場合、それらの係数は頻繁にトレーニングされなければならない。
【０００６】
それらのＦＥＱ係数をトレーニングする１つの仕方は、ブロック１６０によって実行されるチャネル推定を介してである。副搬送波ｋで、時間ｎＴにおいて（Ｔはシンボル区間）、送信機から伝送される信号がＸ（ｎ，ｋ）であり、チャネル伝達関数がＨ（ｎ，ｋ）であり、受信機で受信される信号が、Ｙ（ｎ，ｋ）と仮定する。Ｙ（ｎ，ｋ）＝Ｘ（ｎ，ｋ）Ｈ（ｎ，ｋ）を得る。Ｈ（ｎ，ｋ）が既知の場合、ＦＥＱ係数を１／Ｈ（ｎ，ｋ）に設定することができ、すると、Ｙ（ｎ，ｋ）／Ｈ（ｎ，ｋ）＝Ｘ（ｎ，ｋ）である。チャネル伝達関数は、受信機には未知である。したがって、この場合の課題は、Ｈ（ｎ，ｋ）を推定することである。
【０００７】
Ｘ（ｎ，ｋ）とＹ（ｎ，ｋ）が知られている場合、Ｈ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｘ（ｎ，ｋ）を推定することが可能である。Ｙ（ｎ，ｋ）は受信機において入手可能である。受信機がＸ（ｎ，ｋ）を知るために、通常、いくつかの予め決めたトレーニング信号が、いくつかの特定の時間／周波数において送信機から伝送される。変化のないチャネル、又はゆっくり変化するチャネルの場合、それらのトレーニング信号は、通常、データ伝送が始まる前の初期トレーニング段階において伝送される。その後、Ｘ（ｎ，ｋ）が、典型的には、受信機判定を介して、又は伝送される基準信号を介して得られる。急速に変化するチャネルの場合、基準信号は、受信機がチャネル伝達関数を十分に頻繁に推定してチャネルの変化に追従できるように、予め決めた時間と周波数で送信機から頻繁に伝送されなければならない。基準信号の伝送は、いくらかのチャネル帯域幅を消費して、データ伝送速度の低下をもたらす。したがって、基準信号を過度に伝送することはできない。通常、基準信号は、時間／周波数の小さいパーセンテージだけ伝送される。受信機はそれらのスナップショット・トレーニング信号を利用して、それらの特定の時間／周波数スナップショットにおいてチャネル伝達関数を計算し、次に、他のすべての時間／周波数におけるチャネル伝達関数を推定する。すべての時間／周波数に関するチャネル伝達関数推定値、Ｈ（ｎ，ｋ）を得た後、１／Ｈ（ｎ，ｋ）が、時間ｎＴにおける第ｋ番の副搬送波に関するＦＥＱ係数として使用される。最後に、伝送された信号の推定値が、Ｙ（ｎ，ｋ）／Ｈ（ｎ，ｋ）として得られる。
【０００８】
典型的な例を、ＤＶＢ−Ｔ仕様において見ることができる。ＤＶＢ−Ｔは、２ｋ又は８ｋの副搬送波でＯＦＤＭ変調を使用する。２ｋモードの場合、４５の副搬送波が、連続パイロット・トーンとして使用される。８ｋモードの場合、１７７の副搬送波が、連続パイロット・トーンとして使用される。ＤＶＢ−Ｈ仕様は、ＤＶＢ−Ｔに基づくが、移動体／ハンドヘルド・アプリケーションに合わせられている。ＤＶＢ−Ｈにおいて、さらなる４ｋモードが、決められている。図２は、ＤＶＢ−ＴとＤＶＢ−Ｈにおけるパイロット挿入パターンを示す。図２は、本開示全体で使用される用語を定義するのに使用される。図２において、横の次元は、周波数ドメインを表し、縦の次元は、時間ドメインを表す。それぞれの黒い円は、「パイロット・セル」と呼ばれ、それぞれの白い円は、「データ・セル」又は「非パイロット・セル」と呼ばれる。図２における各行は、個別の「シンボル」に対応し、各列は、「トーン」と呼ばれる。パイロット・セルだけを有する列（左端の列と右端の列などの）は、「連続パイロット・トーン」と呼ばれ、パイロット・セルだけを有する行は、「連続パイロット・シンボル」と呼ばれる。連続パイロット・トーン又は連続パイロット・シンボルの中の各パイロット・セルは、「連続パイロット・セル」と呼ばれる。パイロット・セルとデータ・セルの両方を有する各列は、「散乱パイロット・セル」と呼ばれ、パイロット・セルとデータ・セルの両方を有する各行は、「散乱パイロット・シンボル」と呼ばれる。散乱パイロット・トーン又は散乱パイロット・シンボルの中の各パイロット・セルは、「散乱パイロット・セル」と呼ばれる。データ・セルだけを有する列は、「非パイロット・トーン」と呼ばれ、データ・セルだけを有する行は、「非パイロット・シンボル」と呼ばれる。図２では、「連続パイロット・シンボル」も、「非パイロット・シンボル」も存在しないことに留意されたい。
【０００９】
図２において、すべてのシンボルの中で、いくつかの副搬送波が、散乱パイロット・セルとして使用される。散乱パイロット・セルは、周波数において１２搬送波分離れており、搬送波位置は、シンボルごとに３つずつシフトしている。その結果、散乱パイロット・セルは、時間的に４シンボル離れている。連続パイロット・トーンを除く残りの時間／周波数において、データ信号が伝送される。パイロット信号は、受信機に知られているので、受信機によって、それらの特定の時間／周波数におけるチャネル伝達関数を計算するのに使用される。パイロット信号は、次に、受信機によってチャネルひずみを補正し、データを適切に検出するために、他のすべての時間／周波数における推定チャネル伝達関数Ｈ（ｎ，ｋ）を計算する（補間する）のに使用される。補間器は、時間と周波数において２次元である。この場合の課題は、どのようにしてＨ（ｎ，ｋ）を正確に、効率的に推定するかである。
【００１０】
通常、２次元補間は、２つの個別の１次元補間で実現可能である。補間はまず、すべての散乱パイロット・トーンにおいて時間ドメインで行われる。或る特定の散乱パイロット・トーンｋｉにおいて、パイロット・セルは、時間的に４シンボル離れて送信されるので、Ｘ（ｎ＋４ｍ，ｋｉ）、ｍ＝０，±１，±２．．．、Ｘ（ｎ＋４ｍ，ｋｉ）とＹ（ｎ＋４ｍ，ｋｉ）から、Ｈ（ｎ＋４ｍ，ｋｉ）が得られ、次に、推定に必要なＨ（ｎ＋４ｍ＋１，ｋｉ）、Ｈ（ｎ＋４ｍ＋２，ｋｉ）、Ｈ（ｎ＋４ｍ＋３，ｋｉ）が得られる。これは、時間ドメイン補間である。周波数ドメインにおいて、散乱パイロット・トーンは、３トーン離れている。時間ｎにおける周波数ドメイン補間は、Ｈ（ｎ，ｋ＋３ｊ）、ｊ＝０，±１，±２．．．を使用して、すべての非パイロット搬送波におけるＨ（ｎ，ｋ＋３ｊ＋１）、Ｈ（ｎ，ｋ＋３ｊ＋２）を推定する。
【００１１】
補間オペレーションのいずれも、ＦＩＲ（有限インパルス応答）フィルタを使用して実施可能である。そのようなＦＩＲフィルタは、単に、ローパス・フィルタである補間フィルタであってもよい。固定補間フィルタが使用される場合、ローパス・フィルタの帯域幅は、最悪ケースのチャネル変動を範囲に含まなければならない。ＤＶＢ−Ｔ／ＤＶＢ−Ｈの場合、時間ドメイン補間フィルタは、通過帯域が最悪ケースのドップラー周波数を範囲に含まなければならない１／４通過帯域ローパス・フィルタで、周波数ドメイン補間フィルタは、通過帯域が最悪ケースのマルチパス遅延拡散を範囲に含まなければならない１／３通過帯域ローパス・フィルタである。ローパス補間フィルタは、しばしば、実数の対称係数を使用する。例えば、時間ドメイン補間の場合、フィルタは、Ｈ（４（ｍ−Ｍ／２＋１），ｋｉ）、Ｈ（４（ｍ−Ｍ／２＋２），ｋｉ）、Ｈ（４（ｍ−Ｍ／２＋３），ｋｉ）、，，，Ｈ（４（ｍ−１），ｋｉ）、Ｈ（４（ｍ），ｋｉ）、Ｈ（４（ｍ＋１），ｋｉ）、．．．Ｈ（４（ｍ＋Ｍ／２），ｋｉ）を入力として、さらに３セットのＭ個の実係数をそれぞれ使用して、Ｈ（４ｍ＋１，ｋｉ）、Ｈ（４ｍ＋２，ｋｉ）、Ｈ（４ｍ＋３，ｋｉ）をそれぞれ推定する。Ｈ（４（ｍ＋１），ｋｉ）、．．．Ｈ（４（ｍ＋Ｍ／２），ｋｉ）は、時間４ｍ＋１、４ｍ＋２、４ｍ＋３においては入手できない将来のチャネル伝達関数であることに留意されたい。明らかに、Ｈ（４ｍ＋１，ｋｉ）、Ｈ（４ｍ＋２，ｋｉ）、Ｈ（４ｍ＋３，ｋｉ）は、Ｈ（４（ｍ＋１），ｋｉ）、．．．Ｈ（４（ｍ＋Ｍ／２），ｋｉ）が入手可能になるまで、計算することができない。これらのフィルタは、中心のまわりに基本的に対称的な係数を有し、したがって、ほぼ同数のＨ値が、両方の側で使用される必要がある。このことは、パイロットにおけるいくつかの過去のＨ値に加え、受信された信号の約Ｍ個の「将来の」シンボルを格納するメモリ空間が必要とされることを意味する。例えば、ＤＶＢ−Ｔ  ８ｋモードの場合、Ｍ＝８である場合、受信機は、１６個の「将来の」シンボルを格納する必要があり、これらのシンボルのそれぞれは、約８ｋの複素信号を有し、これは、相当な量のメモリ空間を要する。より小さいＭを使用することにより、メモリ要件と複雑さを減らすことができるが、パフォーマンスが低下する可能性がある。周波数ドメイン補間において、問題は信号の格納ではなく、境界の外でＨ値が利用できないことによる、両方の境界におけるパフォーマンス低下である。実係数を使用することにより、フィルタ特性も制限される。多くの特定のチャネル条件に関して、複素係数を使用することにより、はるかに大きい柔軟性と、パフォーマンス向上の可能性がもたらされる。最適なフィルタ係数セットは、しばしば、最悪ケースの条件に関して、オフラインで事前に計算され、ＲＯＭの中に格納されて、オペレーション・モードに応じて選択される。
【００１２】
より単純な補間方法が、線型補間である。例えば、時間ドメイン補間において、Ｈ（４ｍ＋１，ｋｉ）、Ｈ（４ｍ＋２，ｋｉ）、Ｈ（４ｍ＋３，ｋｉ）が、Ｈ（４ｍ，ｋｉ）とＨ（４ｍ＋４，ｋｉ）に基づいて計算される。線型補間は、はるかに単純であり、周波数ドメイン補間において、境界問題はまったく存在しない。しかし、線型補間のパフォーマンスは、通常、ローパス補間フィルタと比べてはるかに劣る。
【００１３】
時間ドメインにおける補間は、将来の信号を格納するメモリを要求するので、将来の信号をまったく要求しない予測を使用することが、ときとして、望ましい。しかし、通常、補間は、予測よりも良好なパフォーマンスを提供する。
【００１４】
通信チャネルは、雑音を生じさせる。すなわち、Ｙ（ｎ，ｋ）＝Ｘ（ｎ，ｋ）Ｈ（ｎ，ｋ）＋ｗ（ｎ，ｋ）であり、ただし、ｗ（ｎ，ｋ）は、チャネルによって生じる付加的な雑音である。その結果、パイロットで計算されるＨ値の精度は雑音による影響を受ける。推定に対する雑音の影響を最小限に抑える、いくつかの複雑なアプローチが開発されている。そのようなアプローチの１つが、「ウィーナー・フィルタ」又は「最小平均２乗誤差フィルタ」と呼ばれる。チャネル伝達関数相関行列が計算され、雑音を含む急速に変化するチャネル特性の知識が使用されなければならず（ただし、入手可能ではない）、行列反転がその計算に関与する。このアプローチの実施は、極めて複雑である。
【発明の開示】
【発明が解決しようとする課題】
【００１５】
したがって、最適な補間フィルタ係数をオンラインで効果的かつ効率的に計算するための技術が必要である。さらに、補間フィルタが、一方の側で、他方の側よりも少ないタップを使用して（又は、予測のためのタップをまったく使用しないで）非対称的であり、時間ドメイン補間に関する信号格納要件と、周波数ドメイン補間に関する境界効果をともに減らすことが望ましい。また、補間フィルタが、様々なチャネル条件に対して最適化されるように、より大きい柔軟性を有することも望ましい。
【課題を解決するための手段】
【００１６】
無線通信システムにおけるチャネル推定のための方法が、以下のステップを含む。チャネル伝達関数が、連続パイロット・セルと散乱パイロット・セルにおける送信信号と受信信号を使用して、連続パイロット・セルと散乱パイロット・セルにおいて計算される。連続パイロット・セルと散乱パイロット・セルにおいて計算されたチャネル伝達関数を使用して時間ドメイン適応補間が実行され、散乱パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数が得られる。連続パイロット・セルと散乱パイロット・セルにおいて計算されたチャネル伝達関数を使用して周波数ドメイン適応補間が実行されて、非パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数が得られる。
【発明を実施するための最良の形態】
【００１７】
以下の詳細な説明、及び添付の図面は、本発明の性質及び利点のより良好な理解をもたらす。
【００１８】
本発明によれば、補間フィルタ係数がオンラインで適応的に調整され、推定精度を最適化し、平均２乗誤差を最小限に抑える。一実施形態では、低い複雑度のＬＭＳ（最小平均２乗）適応アルゴリズムが、係数適合のために使用される。ＬＭＳアルゴリズムを使用すると、雑音を含むチャネル特性の予備知識は必要とされない。補間フィルタは、最適な設定に自動的に収束して、平均２乗誤差を最小限に抑える。ＬＭＳアルゴリズムにおける計算は簡単であり、複雑度は低い。チャネル条件が変化するにつれ、フィルタ係数は新たな最適設定に自動的に再収束する。したがって、補間フィルタは、可能な限り少ない数のタップを使用して、パフォーマンスを向上させる。
【００１９】
複素係数が、補間フィルタのために使用されて、パフォーマンスを向上させ、フィルタ・スパンを短くし、フィルタ・タップの非対称の配置を可能にする。この非対称性は、時間ドメイン補間において、より少ない将来の信号を使用すること（又は時間ドメイン予測において将来の信号をまったく使用しないこと）を可能にして、格納要件を減らすことに役立ち、さらに、周波数ドメイン補間に関する境界効果を大幅に低減する。
【００２０】
Ｉ．パイロット・トーンにおいてチャネル伝達関数を計算する
パイロット・トーンにおいて、送信信号Ｘ（ｎ，ｋ）が受信機によって知られている。それらのパイロット・トーンにおける受信信号Ｙ（ｎ，ｋ）は、受信機において入手可能である。チャネル伝達関数Ｈ（ｎ，ｋ）は、単に、Ｈ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｘ（ｎ，ｋ）として計算される。
【００２１】
II．時間ドメイン補間
補間はまず、時間ドメインで実行される。時間ドメイン補間が必要とされるのは、一部のパイロット・トーン（例えば、ＤＶＢ−Ｔ／ＤＶＢ−Ｈにおいて散乱パイロット・トーンと呼ばれる）において、パイロットが、一部の時間だけにしか伝送されないからである。ＤＶＢ−Ｔ／ＤＶＢ−Ｈにおいて、散乱パイロット・セルは４つのシンボルおきに伝送される。時間ドメイン補間が実行されて、散乱パイロット・トーンの非パイロット・セルにおけるチャネル伝達関数が推定される。
【００２２】
Ａ．係数更新
最初のステップは、補間フィルタ係数適合である。一実施形態では、ＤＶＢ−Ｔ／ＤＶＢ−Ｈシステムなどのケースにおいて、適合は、連続パイロットに基づく。適合は、以下の３つのサブステップを行う。
【００２３】
ステップ（１）：現在の補間フィルタを使用して、連続パイロット・トーンにおいて補間を実行する。副搬送波ｋ1、ｋ2、．．．ｋPにおいてＰ個の連続パイロット・トーンが存在するものと仮定する。すべてのシンボル（又はシンボルのサブセット）の中で、各連続パイロット・トーンｋi（ｉ＝１，２，．．．Ｐ）において、補間が、現在の補間フィルタ係数セットを使用して実行される。ＤＶＢ−Ｔ／ＤＶＢ−Ｈの場合、３つの時間ドメイン補間フィルタ係数セットが存在する。それら３つのフィルタは、ＨのＭ１個の過去の値、及びＭ２個の将来の値とともに、Ｈ（４（ｍ−Ｍ１＋１），ｋi），Ｈ（４（ｍ−Ｍ１＋２），ｋi），．．Ｈ（４（ｍ−１），ｋi）、Ｈ（４ｍ，ｋi）、Ｈ（４（ｍ＋１），ｋi）、．．Ｈ（４（ｍ＋Ｍ２），ｋi）を使用して、Ｈ（４ｍ＋１，ｋｉ）、Ｈ（４ｍ＋２，ｋｉ）、Ｈ（４ｍ＋３，ｋｉ）をそれぞれ推定する。シンボルｎにおける連続パイロット・トーンにおいて、３つの補間フィルタが実行されて、以下のとおり、３つのＨ（ｎ，ｋｉ）の推定値を得る。すなわち、
H1(n,ki)＝ΣW1(n-1,m)H(n-1+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
H2(n,ki)＝ΣW2(n-1,m)H(n-2+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
H3(n,ki)＝ΣW3(n-1,m)H(n-3+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
である。
【００２４】
ただし、ＷＬ（ｎ−１，ｍ）は、先行するシンボルｎ−１において更新された第Ｌ番（Ｌ＝１，２，３）補間フィルタの第ｍ番の係数である。各フィルタは、Ｍ１＋Ｍ２個の係数を有する。係数の数は、異なるフィルタに関して同一であるか、近いことが可能であるが、異なる数の係数も可能である。さらに、Ｍ２は、一部のフィルタ又はすべてのフィルタに関して０（補間ではなく、予測）であってもよい。
【００２５】
ステップ（２）：推定誤差を計算する。連続パイロット・トーンにおけるチャネル伝達関数Ｈ（ｎ，ｋｉ）は、前段のセクションＩで計算された。推定誤差Ｅ１（ｎ，ｋｉ）、Ｅ２（ｎ，ｋｉ）、Ｅ３（ｎ，ｋｉ）は、単に、Ｈ（ｎ，ｋｉ）と、前述のステップ（１）で得られた補間結果Ｈ１（ｎ，ｋｉ）、Ｈ２（ｎ，ｋｉ）、Ｈ３（ｎ，ｋｉ）との差として計算される。すなわち、
E1(n,ki)=H(n,ki)-H1(n,ki)
E2(n,ki)=H(n,ki)-H2(n,ki)
E3(n,ki)=H(n,ki)-H3(n,ki)
である。
【００２６】
ステップ（３）：補間フィルタ係数を更新する。補間フィルタ係数は、ステップ（２）において計算された誤差を使用して更新される。それら３つの誤差を使用して、以下のとおり、３つの補間フィルタが、それぞれ更新される。すなわち、
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W1(n,m)=W1(n-1,m)+μE1*(n,ki)H(n-1+4m,ki)
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W2(n,m)=W2(n-1,m)+μE2*(n,ki)H(n-2+4m,ki)
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W3(n,m)=W3(n-1,m)+μE3*(n,ki)H(n-3+4m,ki)
である。
【００２７】
*は、複素共役を表し、μは、更新ステップサイズであることに留意されたい。フィルタ係数更新は、以下のとおり、すべての連続パイロットにおいて、又はすべての連続パイロットのサブセットにおいて実行可能である。すなわち、
ｋｉのすべて、又はサブセットにわたる合計に関して、W1(n,m)=W1(n-1,m)+μΣE1*(n,ki)H(n-1+4m,ki)
ｋｉのすべて、又はサブセットにわたる合計に関して、W2(n,m)=W2(n-1,m)+μΣE2*(n,ki)H(n-2+4m,ki)
ｋｉのすべて、又はサブセットにわたる合計に関して、W3(n,m)=W3(n-1,m)+μΣE3*(n,ki)H(n-3+4m,ki)
である。
【００２８】
この方法は、「ブロック更新」と呼ばれる。
【００２９】
Ｂ．散乱パイロット・トーンにおける補間
散乱パイロット・トーンｋｊにおいて、パイロット・セルがシンボルｎにおいて伝送されない場合、パイロット・セルを伝送するｎより前の最も近いシンボルは、シンボルｎ−１又はシンボルｎ−２又はシンボルｎ−３である。したがって、補間フィルタ（Ｗ１、Ｗ２、又はＷ３）の１つを使用して、推定値Ｈ（ｎ，ｋｊ）が計算される。シンボルｎ−１が、散乱パイロットｋｉにおいてパイロットを伝送する場合、以下のとおりである。すなわち、
H(n,ki)＝ΣW1(n,m)(n-1+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2にわたる合計）
である。
【００３０】
シンボルｎ−２が、散乱パイロットｋｉにおいてパイロットを伝送する場合、以下のとおりである。すなわち、
H(n,ki)＝ΣW2(n,m)(n-2+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2にわたる合計）
である。
【００３１】
シンボルｎ−３が、散乱パイロットｋｉにおいてパイロットを伝送する場合、以下のとおりである。すなわち、
H(n,ki)＝ΣW3(n,m)(n-3+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2にわたる合計）
である。
【００３２】
時間ドメイン補間が、すべての散乱パイロット・トーンにおいて実行されて、散乱パイロット・トーンの非パイロット・セルにおける伝達関数推定値が得られる。
【００３３】
Ｃ．更新のための代替の方法
１つの代替の実施形態では、連続パイロットにわたって掃引するにつれ、更新された係数が推定のために使用可能である。連続パイロット・トーンが、搬送波ｋ1、ｋ2、．．．ｋpに存在するものと仮定する。最初に、以下のとおり、ｎ−１として表される前のシンボルの終わりにおける係数を引き継ぐ。すなわち、
W1(n,k1,m)=W1(n-1,kp,m)
W2(n,k1,m)=W2(n-1,kp,m)
W3(n,k1,m)=W3(n-1,kp,m)
である。
【００３４】
ｋi＝ｋ2から始めて、以下のループを実行する。すなわち、
以下のとおり、推定値を計算し、
H1(n,ki)＝ΣW1(n,ki-1,m)H(n-1+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
H2(n,ki)＝ΣW2(n,ki-1,m)H(n-2+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
H3(n,ki)＝ΣW3(n,ki-1,m)H(n-3+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)、
さらに、以下のとおり、誤差を計算し、
E1(n,ki)=H(n,ki)-H1(n,ki)
E2(n,ki)=H(n,ki)-H2(n,ki)
E3(n,ki)=H(n,ki)-H3(n,ki)、
次に、以下のとおり、係数を更新する。すなわち、
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W1(n,ki,j)=W1(n,ki-1,j)+μE1*(n,ki)H(n-1+4m,ki)
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W2(n,ki,j)=W2(n,ki-1,j)+μE2*(n,ki)H(n-2+4m,ki)
m=-M1+1,-M1+2,...-1,0,1,..M2に関して、W3(n,ki,j)=W3(n,ki-1,j)+μE3*(n,ki)H(n-3+4m,ki)
である。ｋiをｋi+1に移し、最後の連続パイロット・トーンｋpまで、このループを繰り返す。
【００３５】
別の代替の実施形態では、ＬＭＳ更新の頻度が、減ら可能である。例えば、すべてのシンボルにおいて、フィルタの１つだけが、更新されることが可能であり、３つのフィルタが順に更新される。もちろん、減らされた更新では、追従パフォーマンスは、低下する。
【００３６】
さらに別の代替の実施形態では、ＬＭＳ更新が、同一のシンボルの中で１回、又は複数回繰り返されて、収束が速められることが可能である。具体的には、１つのシンボルの中の連続パイロット信号に基づき、ＬＭＳ更新がまず、セクションII−Ａ又はセクションII−Ｃで前述したとおり、実行される。更新の後、新たな係数セットが獲得される。ＬＭＳ更新は、その新たな係数セットから始めて、繰り返され、より新しい係数セットがもたらされる。同一のプロセスが、さらなる回数、繰り返されてもよい。これは、獲得時間を短縮することに役立つ可能性がある。
【００３７】
III．周波数ドメイン補間
時間ドメイン補間の後、すべての連続パイロット・トーンと散乱パイロット・トーンにおいてＨ（ｎ，ｋ）を得ている。その結果、すべてのシンボルの中で、３つのトーンおきに、Ｈ（ｎ，ｋ）が、利用できる。一般性を失うことなしに、最初のそのようなトーンが、３の倍数であると仮定して、それらの値を、ｓ＝Ｓ０，Ｓ０＋１，．．．Ｓ０＋Ｓ−１に関して、Ｈ（ｎ，ｓ３）と表す。周波数ドメイン補間において、ｓ＝Ｓ０，Ｓ０＋１，．．．Ｓ０＋Ｓ−１に関して、Ｈ（ｎ，ｓ３＋１）とＨ（ｎ，ｓ３＋２）を計算する。（最初の搬送波、及び最後の搬送波が、パイロット・トーンであるものと想定して。）
【００３８】
Ａ．周波数ドメイン補間フィルタを更新する
１．パイロット・シンボルを伴なって
すべての副搬送波がパイロットを送信するパイロット・シンボルが、存在する（すなわち、連続パイロット・シンボルが、存在する）場合、補間フィルタ更新は、それらのパイロット・シンボルの中で容易に実行可能である。シンボルｎが、搬送波ｋにおける連続パイロット・シンボルである場合、以下のとおり、２回、チャネル伝達関数Ｈ（ｎ，ｋ）を推定する。すなわち、
H1(n,k)＝ΣV1(n-1,j)H(n,k-1+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2)
H2(n,k)＝ΣV1(n-1,j)H(n,k-2+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2)
である。
【００３９】
次に、以下のとおり、推定誤差が、計算される。すなわち、
E1(n,k)=H(n,k)-H1(n,k)
E2(n,k)=H(n,k)-H2(n,k)
である。
【００４０】
これらの誤差を使用して、以下のとおり、補間フィルタ係数が更新される。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V1(n,j)=V1(n-1,j)+μE1*(n,k)H(n,k-1+3j)
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V2(n,j)=V2(n-1,j)+μE2*(n,k)H(n,k-2+3j)
である。
【００４１】
更新は、以下のとおり、すべての搬送波、又はすべての搬送波のサブセットにおいて実行可能である（ブロック更新）。
V1(n,j)=V1(n-1,j)+μΣE1*(n,k)H(n,k-1+3j)（ｋのすべて、又はサブセットにわたる合計）
V2(n,j)=V2(n-1,j)+μΣE2*(n,k)H(n,k-2+3j)（ｋのすべて、又はサブセットにわたる合計）
である。
【００４２】
代替として、搬送波にわたって掃引するにつれ、更新された係数を推定のために使用することができる。最初に、以下のとおり、ｎ−Ｎとして表される前のパイロット・シンボルの更新の終わりにおける係数を引き継ぐ。すなわち、
V1(n,kmin,j)=V1(n-N,kmax,j)
V2(n,kmin,j)=V2(n-N,kmax,j)
である。
【００４３】
ｋ＝ｋｍｉｎ＋１から始めて、以下のループを実行する。すなわち、
以下のとおり、推定値を計算し、
H1(n,k)＝ΣV1(n,k-1,j)H(n,k-1+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2)
H2(n,k)＝ΣV2(n,k-1,j)H(n,k-2+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2)、
さらに、以下のとおり、誤差を計算し、
E1(n,k)=H(n,k)-H1(n,k)
E2(n,k)=H(n,k)-H2(n,k)
次に、以下のとおり、係数を更新する。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V1(n,k,j)=V1(n,k-1,j)+μE1*(n,k)H(n,k-1+3j)
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V2(n,k,j)=V2(n,k-1,j)+μE2*(n,k)H(n,k-2+3j)
である。ｋを１（又は何らかの整数）だけ増やし、最後のｋ、ｋｍａｘまで、このループを繰り返す。
【００４４】
２．パイロット・シンボルを伴なわずに
すべての副搬送波がパイロットを送信するパイロット・シンボルが、まったく存在しない（すなわち、連続パイロット・シンボルが、まったく存在しない）場合、誤差計算において使用されるＨ（ｎ，ｋ）は、利用できない。このことは、ＤＶＢ−Ｔ／ＤＶＢ−Ｈに関して該当する。その場合、最初に、前のシンボルの中で更新されたフィルタ係数に基づき、すべての非パイロット搬送波に関する周波数ドメイン補間を実行し、その後に、フィルタ係数更新を実行することができる。一実施形態では、フィルタ係数更新は、判定に基づくことが可能である。係数更新に関する誤差信号が、以下のとおり、判定に基づいて計算される。すなわち、
ｋ＝３Ｋ＋１である場合、（Ｋは、整数）フィルタ１を使用して、チャネル伝達関数推定値Ｈ（ｎ，ｋ）を、Ｈ（ｎ，ｋ）＝ΣＶ１（ｎ−１，ｊ）Ｈ（ｎ，ｋ−１＋３ｊ）（ｊ＝−Ｊ１＋１，−Ｊ１＋２，．．．−１，０，１，．．Ｊ２にわたる合計）として計算する。
ｋ＝３Ｋ＋２である場合、フィルタ２を使用して、チャネル伝達関数推定値Ｈ（ｎ，ｋ）を、Ｈ（ｎ，ｋ）＝ΣＶ２（ｎ−１，ｊ）Ｈ（ｎ，ｋ−２＋３ｊ）（ｊ＝−Ｊ１＋１，−Ｊ１＋２，．．．−１，０，１，．．Ｊ２）として計算する。
【００４５】
Ｈ（ｎ，ｋ）を計算した後、Ｘ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｈ（ｎ，ｋ）を計算することにより、Ｘ（ｎ，ｋ）の推定値を計算する。Ｘ（ｎ，ｋ）から、送信信号Ｘ（ｎ，ｋ）の推定値は、判定を介して得られるＸ（ｎ，ｋ）である。誤差は、以下のとおり、計算される。すなわち、
Ｅ（ｎ，ｋ）＝Ｈ（ｎ，ｋ）−Ｈ（ｎ，ｋ）
である。
【００４６】
ただし、Ｈ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｘ（ｎ，ｋ）である。Ｅ（ｎ，ｋ）を使用して、以下のとおり、２つのフィルタ係数セット（Ｖ１又はＶ２）の１つが更新される。すなわち、
ｋ＝３Ｋ＋１である場合、Ｖ１が、以下のとおり、更新される。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V1(n,j)=V1(n-1,j)+μE*(n,k)H(n,k-1+3j)
である。
ｋ＝３Ｋ＋２である場合、Ｖ２が、以下のとおり、更新される。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V2(n,j)=V2(n-1,j)+μE*(n,k)H(n,k-1+3j)
である。
【００４７】
可能なすべての搬送波を使用して、以下のとおり、Ｖ１及びＶ２を更新することができる。
ｋ＝３Ｋ＋１である場合、Ｖ１が、以下のとおり、ブロック更新される。すなわち、
V1(n,j)=V1(n-1,j)+μΣE*(n,k)H(n,k-1+3j)（ｋのすべて、又はサブセットにわたる合計）
である。
ｋ＝３Ｋ＋２である場合、Ｖ２が、以下のとおり、ブロック更新される。すなわち、
V2(n,j)=V2(n-1,j)+μΣE*(n,k)H(n,k-2+3j)（ｋのすべて、又はサブセットにわたる合計）
である。
【００４８】
判定誤差は、フィルタ係数更新に影響を与える可能性があり、フィルタが、発散することさえ生じさせる可能性がある。図１８におけるブロック図を使用して、この潜在的な問題に対処が行われる実施形態を説明する。図１８は、減算ブロック６１０、比較ブロック６２０、選択ブロック６３０が追加されていることを除き、図１におけるのと同様であるＯＦＤＭベースの無線受信機のブロック図である。図１８において、誤差Ｅ（ｎ，ｋ）は減算ブロック６１０によって計算され、次に、比較ブロック６２０を介して事前設定された閾値と比較される。例えば、Ｅ（ｎ，ｋ）の実部と虚部の絶対値（つまり、大きさ）が、計算され、事前設定された閾値と比較されて、そのＥ（ｎ，ｋ）が、更新のために使用されるかどうかが決定される。計算された誤差が、閾値より大きい場合、比較ブロック６２０は、選択ブロック６３０を介して０（ゼロ）を選択し、これにより、その搬送波に関する更新をスキップする。計算された誤差が、閾値より小さい場合、比較ブロック６２０は、その搬送波を更新するために、選択ブロック６３０を介して、計算された誤差Ｅ（ｎ，ｋ）を選択する。別の実施形態では、補間は、以下のとおり、最近に更新された係数に基づいて実行される。すなわち、
ｋ＝３Ｋ＋１である場合、（Ｋは、整数）フィルタ１を使用して、チャネル伝達関数推定値Ｈ（ｎ，ｋ）を、H(n,k)＝ΣV1(n,K-1,j)H(n,k-1+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2にわたる合計）として計算する。
ｋ＝３Ｋ＋２である場合、フィルタ２を使用して、チャネル伝達関数推定値Ｈ（ｎ，ｋ）を、H(n,k)＝ΣV2(n,K-1,j)H(n,k-1+3j)(j=-J1+1,-J1+2,...-1,0,1,..J2）として計算する。
【００４９】
次に、新たに計算された推定値Ｈ（ｎ，ｋ）を使用して、Ｘ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｈ（ｎ，ｋ）を計算し、次に、判定を介してＸ（ｎ，ｋ）からＸ（ｎ，ｋ）を計算し、次に、Ｈ（ｎ，ｋ）＝Ｙ（ｎ，ｋ）／Ｘ（ｎ，ｋ）を計算する。誤差Ｅ（ｎ，ｋ）＝Ｈ（ｎ，ｋ）−Ｈ（ｎ，ｋ）を使用して、以下のとおり、Ｖ１又はＶ２が更新される。すなわち、
ｋ＝３Ｋ＋１である場合、Ｖ１が、以下のとおり、更新される。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V1(n,K,j)=V1(n,K-1,j)+μE*(n,k)H(n,k-1+3j)
である。
ｋ＝３Ｋ＋２である場合、Ｖ２が、以下のとおり、更新される。すなわち、
j=-J1+1,-J1+2,...-1,0,1,..J2に関して、V2(n,K,j)=V2(n,K-1,j)+μE*(n,k)H(n,k-1+3j)
である。
【００５０】
第１の搬送波において、係数値は、前のシンボルの終わりからとられる。セクションII−Ｃで説明した代替の時間ドメイン補間更新方法と同様に、周波数ドメイン補間器更新に関するＬＭＳ更新が、搬送波のサブセットに対して実行されて、複雑度を減らす、又は収束を繰り返し速めることが可能である。
【００５１】
図３は、時間ドメイン補間及び周波数ドメイン補間を含む、本発明の実施形態による補間プロセスの流れ図を示す。チャネル推定器１６０は、ブロック１６１、１６２、１６５を含む。ブロック１６１は、連続パイロット・セルと散乱パイロット・セルを含め、パイロット・セルにおけるチャネル伝達関数を計算する。時間ドメイン補間器１６２は、ブロック１６３と１６４を含む。ブロック１６３は、ＬＭＳアルゴリズムに基づいて時間ドメイン補間器係数更新を計算する。ブロック１６４において、チャネル伝達関数は、すべての散乱パイロット・トーンの非パイロット・セルにおける時間ドメイン補間を使用して計算される。時間ドメイン補間の後、すべての散乱パイロット・トーンにおけるチャネル伝達関数の、時間ドメイン補間を介する推定値又は計算された値が、利用可能である。
【００５２】
周波数ドメイン補間器１６５が、すべての非パイロット・トーンにおけるチャネル伝達関数推定値を計算する。ブロック１６５は、ブロック１６６と１６７を含む。ブロック１６６は、ＬＭＳアルゴリズムに基づいて周波数ドメイン補間器係数更新を計算する。ブロック１６７において、チャネル伝達関数が、すべての非パイロット・トーンにおける周波数ドメイン補間を使用して計算される。周波数ドメイン補間の後、すべてのデータ搬送波におけるすべての伝達関数の推定値が、すべての非パイロット・トーンに関する周波数ドメイン補間、又はすべての散乱パイロット・トーンの非パイロット・セルに関する時間ドメイン補間器を介して利用可能である。それらの結果は、ＦＥＱ１７０に送信されて、ＦＥＱ係数がセットアップされる。
【００５３】
IV．補間フィルタ構造
通常、補間フィルタは、ローパス・デジタル・フィルタとして設計される。通過帯域の帯域幅は、補間比に依存する。例えば、ＤＶＢ−Ｔ／ＤＶＢ−Ｈシステムにおける時間ドメイン補間の場合、補間比は、４であり、ローパス・フィルタの通過帯域は、ナイキスト周波数の１／４である。ＤＶＢ−Ｔ／ＤＶＢ−Ｈシステムにおける周波数ドメイン補間の場合、補間比は、３であり、ローパス・フィルタの通過帯域は、ナイキスト周波数の１／３である。ローパス・フィルタは、通常、中央タップのまわりに対称的な実係数を使用して実施される。したがって、セクションIIで説明した時間ドメイン補間フィルタの場合、Ｍ１＝Ｍ２であり、セクションIIIにおける周波数ドメイン補間の場合、Ｊ１＝Ｊ２である。
【００５４】
しかし、時間ドメイン補間において、以下のとおりのフィルタ・オペレーションを有する（第１のフィルタを例にとる）（セクションII参照）。すなわち、
H(n,ki)＝ΣW1(n,m)H(n-1+4m,ki)(m=-M1+1,-M1+2,...-1,0,1,..M2)
である。
【００５５】
０＜ｍ≦Ｍ２に関して、ｎ−１＋４ｍ＞である。したがって、Ｈ（ｎ，ｋｉ）を計算するのに、「将来の」チャネル伝達関数を必要とする。これは、４（Ｍ１＋Ｍ２）個のシンボルに関する連続パイロット・トーン及び散乱パイロット・トーンにおけるすべてのＨ値、及びすべてのデータ搬送波に関する４Ｍ２個のシンボルに関するすべての受信信号Ｙ（ｎ，ｋ）を格納することを要求する。搬送波の数は、極めて大きいので（ＤＶＢ−Ｔにおいて８ｋまで）、これは、相当に大きいメモリ空間を要求する。メモリ要件を小さくするため、Ｍ２を減らすことが望ましい。
【００５６】
一実施形態では、Ｍ２＜Ｍ１である複素非対称補間フィルタが使用される。複素フィルタは、タイプ当たりの計算を２倍にするものの、パフォーマンスを向上させて、タップの数を減らすのに役立ち、特にＭ２を減らすのに役立つ。さらにメモリ要件を減らすのに、Ｍ２は、フィルタのいくつかに関して０であることさえ可能である。例えば、Ｍ２が、第１のフィルタに関して０であるが、その他の２個のフィルタに関してＭ２＝１である場合、すべてのデータ搬送波に関してわずか３つのシンボルを格納すればよい。
【００５７】
周波数ドメイン補間において、補間フィルタへのすべての入力は、同一のシンボル内にある。「将来の格納」問題は、まったく存在しない。主な問題は、境界にある。Ｊ１＝Ｊ２である場合、計算される搬送波の各側で３Ｊ１の搬送波と３Ｊ２の搬送波を必要とする。このことは、補間フィルタに利用可能な搬送波が十分に存在しない境界以外は問題ではない。その結果、Ｈの推定が両端で劣化する。この問題を解決するのに、小さいＪ１、又は小さいＪ２（１ほどに小さい）を有する非対称フィルタを使用することが望ましい。これは、複素フィルタを使用することによって可能にされる。
【００５８】
補間フィルタは、最悪ケースの条件に基づいて事前設計可能である。例えば、時間ドメイン補間の場合、フィルタは、最高のドップラー周波数を範囲に含まなければならず、周波数ドメイン補間の場合、フィルタは、最長の遅延拡散を範囲に含まなければならない。補間フィルタの阻止帯域は、良好な補間精度と雑音除去を確実にするのに十分なだけの減衰を有さなければならない。雑音は、通常、白色であり、すべての周波数を範囲に含むことに留意されたい。フィルタの通過帯域が、全帯域幅の１／３である場合、雑音の２／３は、阻止帯域内にあり、このため、除去可能である。このため、雑音を除去するのに、補間フィルタ通過帯域は、ドップラー周波数又は遅延拡散が範囲に含まれることが可能である限り、できるだけ狭くなければならない。他方、深い阻止帯域減衰は、補間フィルタ・リソースに重い負担をかける。このため、阻止帯域除去は深すぎてはならないが、パフォーマンス低下が無視できるほどであることを確実にするのに十分なだけちょうどの深さでなければならない。
【００５９】
移動環境の場合、ドップラー周波数、遅延拡散、チャネル雑音条件は、頻繁に変化する。事前設計された補間フィルタは、急速に変化する環境に対して、決して最適化されていない。本発明による適応補間フィルタが、あらゆるチャネル条件の下で、最適なパフォーマンスのための最適なフィルタ設定の理想的な解決を提供する。
【００６０】
実際の無線通信環境において、ドップラー周波数、遅延拡散、チャネル雑音条件は、必ずしもすべて最悪の条件にあるわけではない。補間フィルタは、主にフィルタ長（複雑度）制約のため、限られたリソースを有する。本発明による適応補間は、補間フィルタが補間フィルタの設定を最適化して、任意の特定の急速に変化するチャネル条件下でシステム・パフォーマンスを最適化することを可能にする。
【００６１】
複素フィルタ係数の使用により、一部のチャネル条件に関して最適である可能性がある非対称的なフィルタ応答が可能になる。例えば、マルチパス遅延拡散に関して、第１のパスを始まりとして選択して、その他のパスからの信号が、サイクリック・プレフィックスによって扱われることが可能であるようにしなければならない。このことは、他のすべてのパスがより長い遅延を有し、そのため、遅延拡散が片側であることを意味する。その場合、最適な補間フィルタ応答は、片側でなければならない。これは、複素係数を使用し、係数をオンラインで調整することによって可能になる。片側の応答は、雑音の影響を減らして、両側の応答を有する補間フィルタが範囲に含むことができる最大遅延拡散の２倍の大きさの遅延拡散の補正を可能にすることに役立つ（図４及び図５参照）。
【００６２】
図６、図７は、ＬＭＳ適応アルゴリズムと、最悪ケースのドップラー周波数のために設計された固定の係数とを使用する時間ドメイン補間器のパフォーマンスを比較する。図６では５０Ｈｚという比較的低いドップラー周波数が使用され、図７では１５０Ｈｚという比較的高いドップラー周波数が使用される。補間器の長さは、１７、２５、又は３３である。
【００６３】
図８、図９は、ＬＭＳ適応アルゴリズムと、最悪ケースの遅延拡散のために設計された固定の係数とを使用する周波数ドメイン補間器のパフォーマンスを比較する。図８では、［０  ０．２  ０．５  １．６  ２．３  ５．０］μＳという比較的小さい遅延拡散が、使用され、図９では、｛［０  ０．２  ０．５  １．６  ２．３  ５．０］μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋５０μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋１００μＳ｝という大きい遅延拡散が、使用される。補間器の長さは、２１、２７、又は３３である。
【００６４】
図１０は、１２ｄＢ入力ＳＮＲと２１のタップにおける片側の遅延拡散｛［０  ０．２  ０．５  １．６  ２．３  ５．０］μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋５０μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋１００μＳ｝を有する周波数ドメインＬＭＳ適応補間器の応答を示す。また、補間器への入力のＦＦＴも、青で示されている。右から最初の３つのピークが、３つの遅延グループを表す。ＬＭＳ補間器は、それら２つのピークが通過することを許す。その他すべてのピークは、ＬＭＳ適応補間器の無効に入っており、そのため、除去される。
【００６５】
図１１は、時間ドメインＬＭＳ適応補間器、固定補間器、フィルタ入力のＦＦＴを示す。５０Ｈｚという適度のドップラー周波数が、１２ｄＢの入力信号ＳＮＲで使用される。２５のフィルタ・タップが使用される。通過帯域は、適応補間器において、ドップラーがより低いので、固定の元のフィルタよりも狭いことに留意されたい。このより狭い通過帯域は、より少ない雑音が通過することを許し、そのため、パフォーマンスを向上させる。適応補間器における阻止帯域減衰は、このケースに関して、ちょうど十分な減衰を得るように緩和される。元の補間器におけるより深い減衰は必要ではなく、鋭いフィルタは、しばしば、他の問題を生じさせ、パフォーマンスを低下させる。明らかに、適応フィルタは、利用可能なフィルタ・リソースを最適に使用して、出力平均２乗誤差を最小限に抑えることができる。
【００６６】
図１２は、５０Ｈｚのドップラー周波数の下で非対称時間ドメイン複素フィルタを示す。この非対称フィルタでは、わずか１つの「将来の」Ｈ値しか必要とされず、そのため、「将来の」シンボルのためのストレージがはるかに縮小される。図１１は、そのようなフィルタ、及びそのようなフィルタの入力のＦＦＴを示す。
【００６７】
図１４、図１５は、非対称周波数ドメイン補間器を示す。遅延拡散｛［０  ０．２  ０．５  １．６  ２．３  ５．０］μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋５０μＳ，［０  ０．２  ０．５  １．６  ２．３  ５．０］＋１００μＳ｝が、１２ｄＢの入力ＳＮＲ、及び２１のタップにおいて使用される。非対称周波数ドメイン補間器は、図１６に示される境界効果を減らすことに役立つ。この適応補間の収束が、図１７に示される。
【００６８】
以上は、本発明の様々な実施形態の詳細な説明を提供したが、多くの代替形態、変形形態、及び等価形態が可能である。例えば、適応チャネル推定のいくつかの実施形態を、ＤＶＢ−Ｔ及びＤＶＢ−Ｈの文脈で説明したが、本発明による適応チャネル推定は、他の無線標準の実施形態において使用されることも可能である。したがって、本発明の範囲は、説明した実施形態に限定されるべきではなく、代わりに、添付の特許請求の範囲によって定義される。
【図面の簡単な説明】
【００６９】
【図１】ＯＦＤＭベースの無線受信機のブロック図である。
【図２】ＤＶＢ−Ｔ及びＤＶＢ−Ｈにおけるパイロット挿入パターンを示す図である。
【図３】補間器ベースのチャネル推定器を示すブロック図である。
【図４】両側信号の３による補間を示す図である。
【図５】片側信号の３による補間を示す図である。
【図６】時間ドメイン補間器パフォーマンス（５０Ｈｚドップラー）を示すプロットである。
【図７】時間ドメイン補間器パフォーマンス（１５０Ｈｚドップラー）を示すプロットである。
【図８】周波数ドメイン補間器パフォーマンス（小さい拡散）を示すプロットである。
【図９】周波数ドメイン補間器パフォーマンス（大きい拡散）を示すプロットである。
【図１０】周波数ドメインＬＭＳ適応補間器係数のＦＦＴ、及び入力を示すプロットである。
【図１１】時間ドメインＬＭＳ適応補間器係数のＦＦＴ、固定補間器係数、及びフィルタ入力を示すプロットである。
【図１２】非対称時間ドメイン適応補間器係数を示すプロットである。
【図１３】非対称時間ドメイン適応補間器係数のＦＦＴ、及び入力を示すプロットである。
【図１４】非対称周波数ドメイン適応補間器係数を示すプロットである。
【図１５】非対称周波数ドメイン適応補間器係数のＦＦＴ、及び入力を示すプロットである。
【図１６】境界効果を減らすのに役立つ非対称周波数ドメイン適応補間器を示すプロットである。
【図１７】非対称周波数ドメイン補間の収束を示すプロットである。
【図１８】本発明の実施形態によるＯＦＤＭベースの無線受信機を示すブロック図である。
【図１】
2008543186.tif 000003
【図２】
2008543186.tif 000004
【図３】
2008543186.tif 000005
【図４】
2008543186.tif 000006
【図５】
2008543186.tif 000007
【図６】
2008543186.tif 000008
【図７】
2008543186.tif 000009
【図８】
2008543186.tif 000010
【図９】
2008543186.tif 000011
【図１０】
2008543186.tif 000012
【図１１】
2008543186.tif 000013
【図１２】
2008543186.tif 000014
【図１３】
2008543186.tif 000015
【図１４】
2008543186.tif 000016
【図１５】
2008543186.tif 000017
【図１６】
2008543186.tif 000018
【図１７】
2008543186.tif 000019
【図１８】
2008543186.tif 000020
【国際調査報告】
2008543186.tif 000021