(57)【特許請求の範囲】
【請求項１】
移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するためのドップラー修正係数を発生するシステムにおいて、
信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するためのチャンネルインパルス応答決定回路と、
１組の既知の記号及びチャンネルインパルス応答から基準ベクトルを発生するための基準回路と、
上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生するための手段と、
到来信号のノイズ分散を決定するための手段と、
上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を発生するための推定装置と、
を備えたことを特徴とするシステム。
【請求項２】
上記ドップラー修正係数βは、次の式に基づいて発生され、
β＝Ｅ(β)＋ｐH(ｚ−ｙ−Ｅ(ｚ−ｙ))／(δ2／Ｃβ＋ｐHｐ)
但し、Ｅ(β)は、βの予想値であり、ｐHは、ｐをベクトル(ｙi、ｔi)Tとすれば、ｐの複素共役転置マトリクスであり、ｚ−ｙは、エラーベクトルであり、Ｅ(ｚ−ｙ)は、受信サンプルにおける予想エラーであり、δ2は、ノイズ分散であり、そしてＣβは、定数である請求項１に記載のシステム。
【請求項３】
Ｅ(β)及びδ2／Ｃβは、ゼロにセットされる請求項２に記載のシステム。
【請求項４】
ドップラー修正係数βを使用して受信信号サンプルにドップラー修正を適用するための手段を更に備えた請求項１、２又は３に記載のシステム。
【請求項５】
上記信号は、ＴＤＭＡ移動通信システムではタイムスロット内の送信バーストで構成され、そして各タイムスロットに対して新たなドップラー修正係数βが発生される請求項１ないし４のいずれかに記載のシステム。
【請求項６】
以前のタイムスロットからのβに対して発生された値の平均値としてＥ(β)を推定するための手段を備えた請求項５に記載のシステム。
【請求項７】
Ｅ(β)の実数部分は、ゼロにセットされる請求項５に記載のシステム。
【請求項８】
受信した信号サンプルをｅ-j*imag(E(β))tだけデローテーションするための手段を更に備え、但し、Ｅ(β)は、βの予想値であり、次いで、デローテーションした受信信号サンプルを使用して、次の式、
β＝Ｅ(β)＋ｐH(ｚ−ｙ)／(δ2／Ｃβ＋ｐHｐ)
に基づき現在タイムスロットにおけるドップラー修正係数βを発生する請求項１ないし７のいずれかに記載のシステム。
【請求項９】
移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するための方法において、
信号が受信されるチャンネルに対するチャンネルインパルス応答を決定し、
チャンネルインパルス応答における１組の既知の記号から基準ベクトルを発生し、
上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生し、
到来信号のノイズ分散を決定し、
上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を推定し、そして
受信信号サンプルにドップラー修正係数を適用して、ドップラーシフトを補償する、
という段階を含むことを特徴とする方法。
【請求項１０】
上記信号は、ＴＤＭＡ移動通信システムでは一連の送信バーストで構成され、そして各送信バーストに対してドップラー修正係数が発生される請求項９に記載の方法。
【請求項１１】
上記ドップラー修正係数βは次の式に基づいて発生され、
β＝Ｅ(β)＋ｐH(ｚ−ｙ−Ｅ(ｚ−ｙ))／(δ2／Ｃβ＋ｐHｐ)
但し、Ｅ(β)は、βの予想値であり、ｐHは、ｐをベクトル(ｙi、ｔi)Tとすれば、ｐの複素共役転置マトリクスであり、ｚ−ｙは、エラーベクトルであり、Ｅ(ｚ−ｙ)は、受信サンプルにおける予想エラーであり、δ2は、ノイズ分散であり、そしてＣβは、定数である請求項９又は１０に記載の方法。
【請求項１２】
移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するためのドップラー修正係数を発生する装置において、
信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するための手段と、
１組の既知の記号及びチャンネルインパルス応答から基準ベクトルを発生するための手段と、
上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生するための手段と、
到来信号のノイズ分散を決定するための手段と、
上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を発生するための手段と、
を備えたことを特徴とする装置。
【請求項１３】
移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するためのドップラー修正係数を発生するベースステーションコントローラにおいて、
信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するための手段と、
１組の既知の記号及びチャンネルインパルス応答から基準ベクトルを発生するための手段と、
上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生するための手段と、
到来信号のノイズ分散を決定するための手段と、
上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を発生するための手段と、
を備えたことを特徴とするベースステーションコントローラ。
【請求項１４】
移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するためのドップラー修正係数を発生するためのコンピュータプログラムを記録したコンピュータにより読み取り可能な記録媒体であって、上記コンピュータプログラムは、
信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するためのコード部分と、
１組の既知の記号及びチャンネルインパルス応答から基準ベクトルを発生するためのコード部分と、
上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生するためのコード部分と、
到来信号のノイズ分散を決定するためのコード部分と、
上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を発生するためのコード部分と、
を含むことを特徴とする記録媒体。
【発明の詳細な説明】
【０００１】
【技術分野】
本発明は、移動通信システムにおけるドップラーシフト補償の推定に係る。
【０００２】
【背景技術】
移動通信システムでは、ベースステーションに対して移動する移動ステーションから送信された信号が、良く知られたドップラー効果を受け、移動ステーションにより送信された周波数に対しベースステーションに受信される周波数に周波数シフトを生じさせる。この周波数シフトをここでドップラーシフトと称する。ドップラーシフトは、ベースステーションに対する移動ステーションの移動速度及び方向に依存する。従って、ドップラー効果は、ベースステーションに対する移動ステーションの移動方向に基づいて周波数の増加又は減少を与える。ドップラーシフトの大きさは、移動ステーションがベースステーションに対して移動する速度に依存する。
送信ユニット及び受信ユニットの同期が適切でないときには、別の種類の周波数オフセットが生じる。
【０００３】
既存の移動通信設備は、特定チャンネル上の特定信号を選択するベースステーション内の周波数検出回路が信号におけるある量のドップラーシフトを考慮に入れられるという点で、ある種のドップラー補償を与える。
ＡＵ６６４６２６号は、移動ステーションがベースステーションに近づいたり離れたりするときにベースステーションと移動ステーションとの間に伝播する無線信号のドップラーシフトを補償するための方法及び回路構成に関連したものである。伝播時間遅延が減少する状態で移動ステーションがベースステーションに近づくときには、伝播時間遅延が時間と共にいかに変化するかを決定するために伝播時間遅延がある間隔で積分される。この情報は、移動ステーションがベースステーションを越えて移動するときに急激なドップラーシフトを補償するように特定時間における無線信号周波数を変化させるのに使用される。ＴＤＭＡシステムでは、整数の時間フレームに等しい間隔で質問がなされる。ＧＳＭ規格では、時間フレームが８個の連続するタイムスロットより成り、そして１つの送信バーストが特定の移動ステーションとベースステーションとの間を１つの時間フレームにおいて通過する。
【０００４】
従って、このドップラー補償方法では、将来の時間における受信を改善するために過去の信号及び到来する信号に応答して判断がなされ、即ちそれは反作用的なシステムである。従って、到来する信号におけるドップラーシフトをリアルタイムで能動的に補償できるシステムを実施することが要望される。
更に、フェージングが生じるために、ドップラー効果を受ける信号は、振幅及び周波数が変化する。従って、振幅及び周波数の両方においてドップラー効果を推定することが要望される。
【０００５】
【発明の開示】
本発明の１つの特徴によれば、移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するためのドップラー修正係数を発生するシステムにおいて、信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するためのチャンネルインパルス応答決定回路と、１組の既知の記号及びチャンネルインパルス応答から基準ベクトルを発生するための基準回路と、上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生するための手段と、到来信号のノイズ分散を決定するための手段と、上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を発生するための推定装置とを備えたシステムが提供される。
【０００６】
ＧＳＭ規格によれば、ＴＤＭＡ送信バーストは、信号が受信されるチャンネルに対するチャンネルインパルス応答を決定するのに通常使用されるトレーニングシーケンスを備えている。これは、受信したトレーニングシーケンスを記憶されたトレーニングシーケンスでコンボリューションすることによって行われる。従って、トレーニングシーケンスは、基準ベクトルを発生するのに使用する１組の既知の記号を構成することができる。或いは又、受信信号サンプルから推定される記号で１組の既知の記号を構成することもできる。
ＴＤＭＡシステムでは、信号は、各タイムスロットにおける一連の送信バーストより成る。ドップラー修正係数は、各送信バーストに対して発生され、そしてそのバーストの信号サンプルを修正するのに使用できる。
【０００７】
更に、ドップラー効果を修正するために、ドップラー修正係数は、送信ユニットと受信ユニットとの間の誤った周波数同期から生じる周波数オフセットも修正する。
本発明の更に別の特徴によれば、移動通信システムにおいて移動ステーションとベースステーションとの間に送信される信号のドップラーシフトを補償するための方法において、信号が受信されるチャンネルに対するチャンネルインパルス応答を決定し、チャンネルインパルス応答における１組の既知の記号から基準ベクトルを発生し、上記基準ベクトル及び受信した信号サンプルからエラーベクトルを発生し、到来信号のノイズ分散を決定し、上記ノイズ分散、エラーベクトル、受信信号サンプル及びサンプリング時間を使用して、振幅修正を表わす実数部分及び位相修正を表わす虚数部分を有するドップラー修正係数を推定し、そして受信信号サンプルにドップラー修正係数を適用してドップラーシフトを補償するという段階を含む方法が提供される。
【０００８】
【発明を実施するための最良の形態】
本発明を良く理解すると共に、本発明をいかに実施するかを理解するために、添付図面を参照して本発明の好ましい実施形態を詳細に説明する。
図１は、ＧＳＭ規格に基づく移動通信システムにおける通常バーストを示す図である。この図は、ベースステーションに受信されたバーストを表わす。ＧＳＭ規格に基づくＴＤＭＡシステムの場合に、移動ステーションは、ベースステーションコントローラにより割り当てられた周波数チャンネルに変調信号としてバーストを送信する。１つの周波数チャンネルは、８個までのバーストをサポートし、各バーストは各コールに関連され、各コールには、バーストを送信するためのタイムスロットが割り当てられる。ＧＳＭ規格に基づくＴＤＭＡシステムのこれ以上の詳細は、当業者に知られているので、ここでは説明しない。
【０００９】
通常のバーストは、２６ビットのトレーニングシーケンス（ＴＲＳ）を取り巻く５８ビットの２つのパケット（ＤＡＴＡ）で構成される。通常バーストの各端には３つのテールビット（ＴＳ）が追加される。トレーニングシーケンス（ＴＲＳ）は、移動ステーション（ＭＳ）により送信される所定のビットシーケンスであり、これは、ベースステーションコントローラ（ＢＳＣ）に知られている。図１において、受信信号サンプルｙkは、バーストの第ｋビットを表わし、この場合、トレーニングシーケンスのスタートを表わすと仮定する。トレーニングシーケンスは、ベースステーションコントローラにおいて、バーストが送信されるチャンネルのインパルス応答を推定し、少なくとも本発明の実施形態では、ドップラー補償係数を推定するのに使用される。送信される実際の情報は、バーストのデータビット（ＤＡＲＡ）に配置される。
【００１０】
上述したように、移動ステーションからベースステーションへ信号が通過する環境は、とりわけ、移動ステーションとベースステーションとの間の距離、及びエリア内の建物及び他の構造物により生じる干渉によって大幅に変化し得る。その結果、ベースステーションに受信される信号の信号強度及びクオリティも大幅に変化する。更に、移動ステーションが移動するために、ベースステーションに受信される信号がドップラーシフトを受け、これを修正しなければならない。
図２は、ＧＳＭシステムにおいてドップラー補償を実施するのに適した回路１を示す。図２の種々のブロックは、相互接続された個別のエンティティとして示されているが、必ずしも個別の物理的エンティティを表わすものではなく、実施される種々の段階を概略的に表わすことが意図される。ブロックは、回路として実施することもできるし、或いは適当にプログラムされたマイクロプロセッサが、ブロックに個々に指定された機能を実行してもよい。
【００１１】
アンテナ２０は、移動ステーションＭＳから信号１１を受信する。アンテナ２０は、ＲＦ回路２２に接続される。この回路２２は、受信したバーストに対して動作して、周波数を基本帯域周波数にダウンシフトすると共に、バーストをサンプリングして、アナログ信号からデジタルサンプル値ｙiを形成し、但し、ｉはサンプルインデックスである。ＲＦ回路２２は、各バーストに対するサンプル値ｙ0、・・ｙnをサンプル時間ｔ0、・・ｔnと共に、通常、送信信号の予想ビットレートで出力する。ＲＦ回路２２の出力は、ライン２４に沿ってチャンネルインパルス応答（ＣＩＲ）ブロック１０、分散及びエラー計算器１６、フィルタ・イコライザ回路１２、ドップラー修正係数を推定するためのβ推定ブロック３６、及び変換回路４０に供給され、ドップラー修正係数をバーストに適用できるようにする。
メモリ３２は、所定のビットシーケンスであるトレーニングシーケンスｘを保持し、これは、移動ステーションＭＳによりトレーニングシーケンスとして送信されそしてベースステーションで受信されて、サンプルｙTRS＝Ｙk→ｙk+25に変換される。この基準トレーニングシーケンスｘは、基準発生器１４及びチャンネルインパルス応答（ＣＩＲ）ブロック１０に供給される。又、基準発生器１４は、チャンネルインパルス応答ブロック１０から推定されたチャンネルインパルス応答ｈも受信する。
【００１２】
ＣＩＲブロック１０は、受信したトレーニングシーケンスｙTRSを含むバーストを受信し、そして例えば、その受信したトレーニングシーケンスｙTRSと既知のトレーニングシーケンスｘとの間のクロス相関を式１に基づいて計算することにより推定チャンネルインパルス応答ｈを計算する。他の形式の推定方法も使用できる。
ｈ＝ｘｃｏｒｒ（ｘ、ｙTRS）                          （式１）
クロス相関を行う前に、デジタル形態で記憶された既知のトレーニングシーケンスｘは、ＧＳＭ規格に基づき、トレーニングシーケンスが送信のためにＭＳにおいて変調されるやり方と同様に、ｉ、ｑ変調されることが明らかである。クロス相関は、５つのタップ値（ｈ(ｊ)j=σto4）の形態のチャンネルインパルス応答を発生するように既知のやり方で行われる。
【００１３】
良く知られているように、推定されたインパルス応答ｈは、データが同じ平均ノイズを受けたかのように、受信バーストにおけるデータの予想推定値を計算するのに使用される。
ＣＩＲブロックは、タイミング進み情報τも発生し、これは、割り当てられたタイムスロットのどこに受信バーストが位置するか決定するのに使用される。
各バーストに対し、そのバーストの推定チャンネルインパルス応答ｈは、ＣＩＲブロック１０により計算され、そしてフィルタ／イコライザ回路に供給され、該回路は、そのバーストのデータＤＡＴＡ（ｒ）を回復できるようにする。良く知られたように、フィルタ／イコライザ回路１２は、受信バーストに対するチャンネルインパルス応答ｈ及びタイミング情報τを受信し、その信号を復調、フィルタリング及びデコードして、データを既知のやり方で回復することができる。
【００１４】
基準発生器１４は、基準ベクトルｚを発生し、これは、既知のトレーニングシーケンス及びインパルス応答のコンボリューションを使用して計算される。基準発生器１４は、次の計算を行う。
ｚ＝ｈ＊ｘ                                    （式２）
但し、ｘは、使用する変調に基づく既知の記号値を表わす。
ｚkは、ベクトルｚのｋ番目のサンプルを表わし、そしてｚは、ベクトルサンプルｚkの対角マトリクスを表わし、即ちｚ＝ｘｈである。ここで、ｘは、既知のトレーニングシーケンス記号の対角マトリクスである。
基準ベクトルｚは、基準発生器から分散及びエラー計算器１６へ供給される。上述したように、この分散計算器は、受信したトレーニングシーケンスｙTRSを含むサンプルｙiも受信する。この分散計算器は、次の式に基づいて分散ｖａｒ（σ2）を計算する。
【数１】
0004542712.tif 000002
ここで、項ref_lengthは、基準ベクトルｚの長さを表わす定数である。これは、サンプルの数（２２）にビット分離を乗算することにより計算される。
【００１５】
式４において、ｙiの値は、バースト（ｙk→ｙk+25）に対する受信したトレーニングシーケンスのサンプル値である。
又、分散及びエラー計算器は、エラーベクトルｚ−ｙ（即ち、各サンプルに対するエラー値ｚi−ｙi）も発生する。
エラーベクトルｚ−ｙ及び分散σ2は、β推定ブロック３６へ供給される。β推定ブロック３６は、ＲＦ回路２２から信号サンプル及び時間ｙi、ｔiも受信する。ここで、時間ｔiは、値ｙiのサンプリング時間である。β推定ブロックは、プログラムを実行するプロセッサとして実施することもできるし、又は他の適当な手段によって実施することもできる。これは、実数部及び虚数部を有する複素数であるドップラー補償係数βを発生するように動作する。従って、βは、次のように表わされる。
α＋ｊφ
但し、αは、フェージングによる振幅変化を表わし、そしてφは、ドップラー角度である。β推定ブロック３６は、ドップラー修正係数βを発生するに必要な値を保持するためのメモリも備え、これらの値は、Ｃβ、Ｅ(β)及びＥ(ｚ−ｙ)を含む。これらの項は、図３を参照して以下に説明する。
【００１６】
ドップラー修正係数βは、各バーストに対する信号サンプルｙi及びｔiと共に、変換回路４０に供給され、次の式に基づいてドップラー修正が行われる。
ｙn-compensated＝ｙn(１−βｔn)
変換回路４０の別の実施は、振幅及び位相成分を別々の修正ステップに分離することである。位相は、サンプルにｅ-j*imag(β)tnを乗算することにより修正でき、そして振幅は、個別のステップで修正できる。
ｙn-phase-compensated＝ｙn＊ｅ-j*imag(β)tn
ｙn-compensated＝ｙn-phase-compensated（１−ｒｅａｌ(β）ｔn）
（式５）
上記式によれば、ｉｍａｇ(β)＝φ、及びｒｅａｌ(β)＝αである。従って、受信信号は、ｒｅａｌ(β)が負であるか正であるかに基づいてバースト中に増幅又は変更される。これは、１組の補償された信号サンプルｙcompを発生し、これらサンプルは、フィルタ・イコライザ回路１２に供給され、バーストのデータとして改善された値を発生することができる。補償された信号は、補償された位相及び振幅の両者を有するか、又は、位相のみを有する。
【００１７】
トレーニングシーケンスｘを使用して、既知の記号を基準として供給した上記形態とは別に、イコライザ回路３０で推定されるビットを使用することができる。例えば、ビタビイコライザのイコライザ回路３０は、フィルタ・イコライザ回路１２からフィルタされ、復調されそしてイコライズされた信号ＤＡＴＡ（Ｒ）を受け取る。イコライザ回路３０は、バーストのデータシーケンスデータの一部分（図１の推定ブロックから導出された部分）に作用して、移動ステーションから送信された出力ビットを推定する。この出力は、ここでｅｓｔｉｍ＿ｂｉｔｓと称する。イコライザ回路３０は、既知の移動通信システムの場合のようにビットの判断を行うように動作し、従って、ここでは詳細に説明しない。これら推定されたビットは、トレーニングシーケンスに代わって使用され、基準ベクトルｚを計算するための既知の値を供給することができる。ｅｓｔｉｍ＿ｂｉｔｓは、チャンネルインパルス応答ｈと共に基準回路３４に供給されて、エラーベクトルｚ’−ｙが発生される。但し、ｚ’＝ｅｓｔｉｍ＿ｂｉｔｓ＊ｈである。
【００１８】
ドップラー修正係数βの発生の基礎となる理論について以下に説明する。
（純粋な）ドップラー周波数推定の問題は、次のように書くことができる。
Ｘｈ＝Ｙｒ＋ｗ
但し、ｒは回転ベクトル、Ｙは対角サンプルマトリクス、そしてｗはノイズである。
【数２】
0004542712.tif 000003
Ｘｈは、既知の記号（トレーニングシーケンス記号の対角マトリクスｘ）とインパルス応答ｈとの間のコンボリューションに対応する。
従って、ドップラー周波数推定では、ドップラー角φを推定しなければならない。この問題は、非直線的であり、従って、特にリアルタイムで実施する場合に解決されるべき手の込んだビット問題である。直線化ステップは、ｅjΦt0＝１＋ｊφｔ0を仮定して行うことができる。フェージングのために、必要とされるドップラー周波数修正は、振幅が一定であると仮定することができない。
項ｊφｔ0が実数成分も有する場合には、振幅変化をおおよそ推定することができる。従って、項ｊφｔ0は、βを複素数（α＋ｊφ）とすれば、より一般的な仮定１＋βｔ0に置き換えられる。
より一般的には、信号サンプルyiは、eβを乗ずることによってドップラー効果に対して修正されうる。ここで、βの虚数部は、ドップラー角φを表し、実数部は、振幅αを表す。従って、直線近似eβ≒１＋βを使用して、振幅変動が概略的に推定されうる。
【００１９】
直線的な問題は、次のように書くことができる（簡単化のために、ｚ＝Ｘｈ＝［ｚ0,ｚ1,,・・,ｚn-1］T）。
ｚi＝ｙi（１＋β(ｔi)）＋ｗi
これは、マトリクス表示で次のように再フォーマットすることができる。
ｚ−ｙ＝ｐβ＋ｗ                                    （式６）
ここで、ｐ＝［ｙ0ｔ0,ｙ1ｔ1,・・,ｙn-1ｔn-1］T及びβを推定しなければならない。
式６の解は、ＬＭＭＳＥ推定子を使用して見出すことができる。
β＝Ｅ(β)＋Ｃβ(z-y)Ｃ(z-y)(z-y)-1（ｚ−ｙ−Ｅ(ｚ−ｙ)）
但し、Ｅ(β)は、βの予想値であり、Ｃβ(z-y)、Ｃ(z-y)(z-y)は、２つのエレメント間の共分散マトリクス、即ちＥ（β(ｚ−ｙ)H））及びＥ（(ｚ−ｙ)(ｚ−ｙ)H）であり、そしてＥ(ｚ−ｙ)は、受信バーストにおける予想エラーである。
【００２０】
この段階において、ノイズは、δ2の分散を有するホワイトガウスノイズであり、そしてある操作の後に、次のようになると仮定することができる。
β＝Ｅ(β)＋ｐH(ｚ−ｙ−Ｅ(ｚ−ｙ))／(δ2／Ｃβ＋ｐHｐ)    （式７）
但し、Ｃβ＝Ｅ((β−Ｅ(β))(β−Ｅ(β))H)                    （式８）
ｐHは、ｐの複素共役転置マトリクスである。
式７は、β推定ブロック３６に使用されて、βを発生する。β推定ブロックは図３に概略的に示されている。実際に、これは、必要なデータを保持するメモリと共にプロセッサにより実行されるソフトウェアとして実施することができる。β推定回路５０は、各タイムスロットにおいて受け取られたデータ（ｙi、ｔi）に作用してβの値を発生し、これは、変換回路４０により使用される。β累積ブロック５２は、先行するタイムスロットにわたってβ値を累積及び平均化し、Ｅ(β)を発生する。β累積機能は、適応フィルタにより実施することができる。（ｚ−ｙ）累積ブロック５４は、予想エラーＥ(ｚ−ｙ)を発生する。これを行うために、先行するタイムスロットからのｚ−ｙの値、又はβ累積ブロック５２からのＥ(β)の値を使用して、周波数エラーを推定することができる。Ｃβ発生ブロック５６は、Ｃβの値を発生する。式７のパラメータに対して多数の可能性が存在する。
【００２１】
Ｅ(β)がゼロにセットされ、そしてδ2／Ｃβもゼロにセットされる場合には、式７は、通常の直線的回帰を解く。
上述したように、βは、各タイムスロットに対して発生されることが意図される。予想値Ｅ(β)は、以前のタイムスロットから推定β値の平均を計算することにより発生できる。従って、これらの値は、各タイムスロットに対して計算されたときに保持することができ、そしてβ推定ブロック３６においてＥ(β)ブロックに保持された平均値を発生するのに使用できる。
別の簡単化は、Ｅ(β)の実数部がゼロであると仮定することである。これは、典型的なフェージングチャンネル状態では長時間の振幅変化が予想されないので、適切な仮定である。従って、周波数のエラー成分に関してのみ長時間平均をとれば充分である。又、位相変化だけを修正すべき場合にも、実数部がゼロと仮定することができる。
【００２２】
受信したサンプルｙiをｅ-j*imag(β)tだけデローテーション(derotating)することにより予想エラーＥ(ｚ−ｙ)を計算する必要性を排除することができる。このため、βは、以前のタイムスロットからのβの初期推定値でよい。この場合に、式７は、次のように変化する。
β＝Ｅ(β)＋ｐH(ｚ−ｙ)／(δ2／Ｃβ＋ｐHｐ)                （式９）
従って、これは、ベクトルｐを形成する前にサンプルｙiをｅ-j*imag(β)tだけデローテーションするような前処理ステップ（図３の５０ａを参照）を必要とする。
Ｃ(β)の値は、種々の解決策を使用してＣβ発生ブロックにより発生することができる。これは定数値であってもよい。これは、以前のタイムスロットの情報を使用し、例えば、式８を使用して推定することができる。これは、βの予想値及び／又はノイズ分散δ2の関数であってもよい。これは、適応フィルタを使用して推定することもできるし、或いは上記の組合せによって推定することもできる。
この技術は、短い遅延の拡散チャンネル及び純粋な周波数エラーに特に適している。
【図面の簡単な説明】
【図１】  移動通信システムにおける信号バーストを示す図である。
【図２】  ドップラーシフト補償を実施するための回路のブロック図である。
【図３】  β推定ブロックを示す図である。
【図１】
0004542712.tif 000004
【図２】
0004542712.tif 000005
【図３】
0004542712.tif 000006