(57)【要約】  本発明は、送信元と受信機との間のチャネルのチャネル推定を実行する方法に関し、周波数依存分散を有する雑音が、送信元からチャネルを通じて受信機によって受信されるシンボルのストリームに加えられ、受信機は、周波数領域において復調を実行する。この方法は、  受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出し、  パイロット位置における雑音分散を推定し、  パイロット位置において抽出されたシンボルから粗チャネル推定値を求め、  推定された雑音分散から重み付け係数を求め、  求められた重み付け係数によって粗チャネル推定値を重み付けし、  重み付き粗チャネル推定値を所定の係数を用いてフィルタリングして、なお、少なくとも１つの所定の係数は他の所定の係数と異なり、所定の係数は定雑音分散について求められ、  求められた重み付け係数及び所定のフィルター係数から正規化係数を求め、  フィルタリングされた重み付き粗チャネル推定値を求められた正規化係数を用いて正規化する。
2019511169.tif 000002
【特許請求の範囲】
【請求項１】
送信元と受信機との間のチャネルのチャネル推定を実行する方法であって、周波数依存分散を有する雑音が、前記送信元から前記チャネルを通じて前記受信機によって受信されるシンボルのストリームに加えられ、前記受信機は、周波数領域において復調を実行し、該方法は、
前記受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出するステップと、
パイロット位置における雑音分散を推定するステップと、
パイロット位置において抽出された前記シンボルから粗チャネル推定値を求めるステップと、
前記推定された雑音分散から重み付け係数を求めるステップと、
前記求められた重み付け係数によって前記粗チャネル推定値を重み付けするステップと、
前記重み付き粗チャネル推定値を所定の係数を用いてフィルタリングするステップであって、少なくとも１つの所定の係数は他の所定の係数と異なり、前記所定の係数は定雑音分散について求められる、ステップと、
前記求められた重み付け係数及び前記所定のフィルター係数から正規化係数を求めるステップと、
前記フィルタリングされた重み付き粗チャネル推定値を前記求められた正規化係数を用いて正規化するステップと、
を含むことを特徴とする、方法。
【請求項２】
前記方法は、復調されたシンボルと、雑音電力推定値と、前記正規化された、フィルタリングされた重み付き粗チャネル推定値とを用いて、前記チャネルを等化するステップを更に含むことを特徴とする、請求項１に記載の方法。
【請求項３】
前記重み付け係数は、
【数１】
2019511169.tif 000050
に等しく、ｋはパイロット位置インデックスであり、σｋ２は、パイロット位置インデックスｋの雑音プラス干渉体分散であり、σｂ２は、定雑音分散であることを特徴とする、請求項１又は２に記載の方法。
【請求項４】
前記所定の係数ａｉｋは、所定の行列Ａｂａｓｉｃの形態でグループ化され、該行列は、以下のように求められ、
【数２】
2019511169.tif 000051
ここで、Σ（・）は、チャネルベクトルの相関行列を表し、Ｉは単位行列であることを特徴とする、請求項１から３のいずれか１項に記載の方法。
【請求項５】
前記正規化係数λｉは、以下のように求められ、
【数３】
2019511169.tif 000052
ここで、ｉは、前記フィルタリングの第ｉの出力のインデックスであり、μｉは、所定の値であるか又は前記係数ａｉｋの和であることを特徴とする、請求項１から４のいずれか１項に記載の方法。
【請求項６】
複数の所定の行列が記憶され、前記方法は、前記推定された雑音分散及び／又はチャネルベクトルの相関行列Σｈに関連したパラメーターの推定値に従って、前記複数の所定の行列のうちの１つを選択するステップを更に含むことを特徴とする、請求項１から５のいずれか１項に記載の方法。
【請求項７】
前記復調は、ＯＦＤＭ復調又はＳＣ−ＦＤＭＡ復調又はＳＣ−ＯＦＤＭ復調であることを特徴とする、請求項１から６のいずれか１項に記載の方法。
【請求項８】
送信元と受信機との間のチャネルのチャネル推定を実行するデバイスであって、周波数依存分散を有する雑音が、前記送信元から前記チャネルを通じて前記受信機によって受信されるシンボルのストリームに加えられ、前記受信機は、周波数領域において復調を実行し、該デバイスは、
前記受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出する手段と、
パイロット位置における雑音分散を推定する手段と、
パイロット位置において抽出された前記シンボルから粗チャネル推定値を求める手段と、
前記推定された雑音分散から重み付け係数を求める手段と、
前記求められた重み付け係数によって前記粗チャネル推定値を重み付けする手段と、
前記重み付き粗チャネル推定値を所定の係数を用いてフィルタリングする手段であって、少なくとも１つの所定の係数は他の所定の係数と異なり、前記所定の係数は定雑音分散について求められる、手段と、
前記求められた重み付け係数及び前記所定のフィルター係数から正規化係数を求める手段と、
前記フィルタリングされた重み付き粗チャネル推定値を前記求められた正規化係数を用いて正規化する手段と、
を備えることを特徴とする、デバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、包括的には、送信元と受信機との間のチャネルのチャネル推定を実行する方法及びデバイスに関する。
【背景技術】
【０００２】
無線電気通信では、地上通信又は衛星通信のいずれにおいても、従来から、チャネル推定が、パイロットシンボル、すなわち、幾つかの状況において異なる信号対雑音比を受ける受信パイロットから行われている。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
ＯＦＤＭ、ＳＣ−ＯＦＤＭ及びＳＣ−ＦＤＭＡのように、復調が周波数領域において行われるシステムでは、送信は、有色雑音、すなわち、周波数依存分散を有する雑音を受けるおそれがある。
【０００４】
雑音のこの周波数依存性は、狭帯域干渉体若しくは隣接する干渉体、又は任意の種類の干渉体に起因している可能性がある。
【０００５】
従来から、そのようなチャネル推定は、一般に受信パイロットに基づいており、パイロット位置に対するチャネル相関を考慮に入れた雑音低減及び補間フィルターを含む。雑音低減及び補間フィルターは、受信パイロットシンボルに影響を及ぼす信号対雑音比に依存する。そのようなフィルターを計算する、例えばＭＭＳＥ／ウィナーフィルターのような解決策が存在する。しかしながら、その結果もたらされるフィルター係数計算は、かなり複雑であり、有色雑音分散を事前に知ることができないことから、リアルタイムで行われることになる。
【課題を解決するための手段】
【０００６】
本発明は、良好なチャネル推定性能を維持しながら、そのようなフィルター係数計算の複雑度を低減することを可能にする方法及びデバイスを提供することを目的とする。
【０００７】
このために、本発明は、送信元と受信機との間のチャネルのチャネル推定を実行する方法であって、周波数依存分散を有する雑音が、送信元からチャネルを通じて受信機によって受信されるシンボルのストリームに加えられ、受信機は、周波数領域において復調を実行し、この方法は、
受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出するステップと、
パイロット位置における雑音分散を推定するステップと、
パイロット位置において抽出されたシンボルから粗チャネル推定値を求めるステップと、
推定された雑音分散から重み付け係数を求めるステップと、
求められた重み付け係数によって粗チャネル推定値を重み付けするステップと、
重み付き粗チャネル推定値を所定の係数を用いてフィルタリングするステップであって、少なくとも１つの所定の係数は他の所定の係数と異なり、所定の係数は定雑音分散について求められる、ステップと、
求められた重み付け係数及び所定のフィルター係数から正規化係数を求めるステップと、
フィルタリングされた重み付き粗チャネル推定値を求められた正規化係数を用いて正規化するステップと、
を含むことを特徴とする、方法に関する。
【０００８】
また、本発明は、送信元と受信機との間のチャネルのチャネル推定を実行するデバイスであって、周波数依存分散を有する雑音が、送信元からチャネルを通じて受信機によって受信されるシンボルのストリームに加えられ、受信機は、周波数領域において復調を実行し、このデバイスは、
受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出する手段と、
パイロット位置における雑音分散を推定する手段と、
パイロット位置において抽出されたシンボルから粗チャネル推定値を求める手段と、
推定された雑音分散から重み付け係数を求める手段と、
求められた重み付け係数によって粗チャネル推定値を重み付けする手段と、
重み付き粗チャネル推定値を所定の係数を用いてフィルタリングする手段であって、少なくとも１つの所定の係数は他の所定の係数と異なり、所定の係数は定雑音分散について求められる、手段と、
求められた重み付け係数及び所定のフィルター係数から正規化係数を求める手段と、
フィルタリングされた重み付き粗チャネル推定値を求められた正規化係数を用いて正規化する手段と、
を備えることを特徴とする、デバイスに関する。
【０００９】
したがって、チャネル推定は、定雑音分散のみの所定の係数を用いた簡単なフィルタリングに関して、最適なフィルターのリアルタイム計算と比較して低複雑度を有するものに改善される。
【００１０】
特定の特徴によれば、方法は、復調されたシンボルと、雑音電力推定値と、正規化された、フィルタリングされた重み付き粗チャネル推定値とを用いて、チャネルを等化するステップを更に含む。
【００１１】
したがって、チャネル推定が改善され、全体性能が改善されるので、復調は改善される。
【００１２】
特定の特徴によれば、重み付け係数は、
【数１】
2019511169.tif 000003
に等しく、ｋはパイロット位置インデックスであり、σｋ２は、パイロット位置インデックスｋの雑音プラス干渉体分散であり、σｂ２は、定雑音分散である。
【００１３】
したがって、重み付き係数は簡単に求められ、大きな影響を受けるパイロットが他の位置におけるチャネル推定にあまり影響を及ぼさないことが保証される。
【００１４】
特定の特徴によれば、所定の係数ａｉｋは、所定の行列Ａｂａｓｉｃの形態でグループ化され、この行列は、以下のように求められる。
【数２】
2019511169.tif 000004
【００１５】
ここで、Σ（・）は、チャネルベクトルの相関行列を表し、Ｉは単位行列である。
【００１６】
したがって、基本行列は、非有色雑音について最適化され、パイロット推定方式は、雑音が有色でない場合、例えば、干渉体が存在しないときに最適である。
【００１７】
特定の特徴によれば、正規化係数λｉは、以下のように求められる。
【数３】
2019511169.tif 000005
【００１８】
ここで、ｉは、フィルタリングの第ｉの出力のインデックスであり、μｉは、所定の値であるか又は係数ａｉｋの和である。
【００１９】
したがって、正規化係数は簡単に求められ、その結果得られるチャネル推定値が最適なものに近い平均エネルギーを有することが保証される。
【００２０】
特定の特徴によれば、複数の所定の行列が記憶され、方法は、推定された雑音分散及び／又はチャネルベクトルの相関行列Σｈに関連したパラメーターの推定値に従って、複数の所定の行列のうちの１つを選択するステップを更に含む。
【００２１】
したがって、チャネル推定プロセスは更に改善され、所定の基本フィルターは最適なものにより近くなる。
【００２２】
特定の特徴によれば、復調は、ＯＦＤＭ復調又はＳＣ−ＦＤＭＡ復調又はＳＣ−ＯＦＤＭ復調である。
【００２３】
したがって、復調は、通常、周波数領域において行われ、信号を周波数領域に変換することによる複雑度は追加されない。
【００２４】
更に別の態様によれば、本発明は、プログラマブルデバイス内に直接ロード可能とすることができるコンピュータープログラムであって、このコンピュータープログラムがプログラマブルデバイスにおいて実行されると本発明による方法のステップを実施する命令又はコード部を含む、コンピュータープログラムに関する。
【００２５】
このコンピュータープログラムに関する特徴及び利点は、本発明による方法及び装置に関連して上記で述べたものと同じであるので、ここでは繰り返さないことにする。
【００２６】
本発明の特徴は、一例示の実施形態の以下の説明を読むことによってより明らかになる。この説明は、添付図面に関して作成されたものである。
【図面の簡単な説明】
【００２７】
【図１】本発明が実施される無線システムを表す図である。
【図２】本発明が実施される受信機のアーキテクチャを表す図である。
【図３】受信機の無線インターフェースの構成要素のブロック図である。
【図４】本発明による受信機によって実行されるアルゴリズムの一例を開示する図である。
【図５】本発明が受信機において実施されたときのパケットエラーレートの進展の一例を開示する図である。
【発明を実施するための形態】
【００２８】
図１は、本発明が実施される無線システムを表している。
【００２９】
本発明は、送信元Ｓｒｃによって転送された信号が受信機Ｒｅｃに送信される一例において開示される。送信元Ｓｒｃは、衛星Ｓｔ又は地上送信機Ｔｔに備えることができる。
【００３０】
簡潔にするために、図１には、１つの衛星Ｓｔ及び１つの地上送信機Ｓｔしか示されていないが、無線リンクは、これよりも多くの数の衛星Ｓｔ及び／又は地上送信機Ｓｔを備えることができる。
【００３１】
衛星Ｓｔは、信号をブロードキャストすることもできるし、デジタルビデオブロードキャスティング次世代ハンドヘルドブロードキャスティングシステム（ＤＶＢ  ＮＧＨ：Ｄｉｇｉｔａｌ  Ｖｉｄｅｏ  Ｂｒｏａｄｃａｓｔｉｎｇ  Ｎｅｘｔ  Ｇｅｎｅｒａｔｉｏｎ  Ｂｒｏａｄｃａｓｔｉｎｇ  ｓｙｓｔｅｍ  ｔｏ  ｈａｎｄｈｅｌｄ）物理層仕様の２０１２年１１月に公表されたＤＶＢ文書Ａ１６０に開示されているような衛星通信を行うこともできる。
【００３２】
地上送信機Ｔｔは、３ＧＰＰ  ＴＳＧ−ＲＡＮ，ＴＲ２５．８１４：「Ｐｈｙｓｉｃａｌ  Ｌａｙｅｒ  Ａｓｐｅｃｔｓ  ｆｏｒ  Ｅｖｏｌｖｅｄ  ＵＴＲＡ」，バージョン７．１．０，２００６−０９の仕様に準拠した無線セルラー電気通信ネットワークの基地局とすることができる。
【００３３】
本発明は、ブロードキャスト又はこれまでの通信のために、ＳＣ−ＦＤＭＡ、ＳＣ−ＯＦＤＭ及びＯＦＤＭのような変調に適用することができる。
【００３４】
簡潔にするために、図１には、１つの受信機Ｒｅｃしか示されていないが、信号は、これよりも多くの数の受信機Ｒｅｃにブロードキャスト又は送信される。
【００３５】
受信機Ｒｅｃは、モバイル端末であってもよい。
本発明によれば、受信機Ｒｅｃは、
受信されたシンボルのストリームから、パイロット位置においてシンボルを抽出し、
パイロット位置における雑音分散を推定し、
パイロット位置において抽出されたシンボルから粗チャネル推定値を求め、
推定された雑音分散から重み付け係数を求め、
求められた重み付け係数によって粗チャネル推定値を重み付けし、
重み付き粗チャネル推定値を所定の係数を用いてフィルタリングして、なお、少なくとも１つの所定の係数は他の所定の係数と異なり、所定の係数は定雑音分散について求められ、
求められた重み付け係数及び所定のフィルター係数から正規化係数を求め、
フィルタリングされた重み付き粗チャネル推定値を求められた正規化係数を用いて正規化する。
【００３６】
図２は、本発明が実施される受信機のアーキテクチャを表す図である。
【００３７】
受信機Ｒｅｃは、例えば、バス２０１によって互いに接続された構成要素と、図４に開示されるようなプログラムによって制御されるプロセッサ２００とに基づくアーキテクチャを有する。
【００３８】
ここで、受信機Ｒｅｃは、専用集積回路に基づくアーキテクチャを有することができることに留意すべきである。
【００３９】
バス２０１は、プロセッサ２００を、読み出し専用メモリＲＯＭ２０２、ランダムアクセスメモリＲＡＭ２０３及び無線インターフェース２０５にリンクする。
【００４０】
メモリ２０３は、変数と、図４に開示されるようなアルゴリズムに関連したプログラムの命令とを収容するように意図されたレジスタを含む。
【００４１】
プロセッサ２００は、無線インターフェース２０５の動作を制御する。
【００４２】
読み出し専用メモリ２０２は、図４に開示されるようなアルゴリズムに関連したプログラムの命令を含む。これらの命令は、受信機Ｒｅｃに電源が投入されると、ランダムアクセスメモリ２０３に転送される。
【００４３】
無線インターフェース２０５は、受信パイロットが種々の信号対雑音比を受けているパイロットから行われる信頼性のあるチャネル推定を実行する手段を備える。
【００４４】
例えば、本発明は、ＯＦＤＭ及びＳＣ−ＦＤＭＡ（ＳＣ−ＯＦＤＭ）のように、復調が周波数領域において行われるシステムに関する。送信は、有色雑音、すなわち周波数依存分散を有する雑音を受ける。
【００４５】
雑音のこの周波数依存性は、狭帯域干渉体若しくは隣接する干渉体、又は任意の種類の干渉体に起因する可能性がある。
【００４６】
従来、Ｓｔｅｆａｎ  ＫａｉｓｅｒのＰｈＤレポート「Ｍｕｌｔｉ−Ｃａｒｒｉｅｒ  ＣＤＭＡ  Ｍｏｂｉｌｅ  Ｒａｄｉｏ  Ｓｙｓｔｅｍｓ  −  Ａｎａｌｙｓｉｓ  ａｎｄ  Ｏｐｔｉｍｉｚａｔｉｏｎ  ｏｆ  Ｄｅｔｅｃｔｉｏｎ，  Ｄｅｃｏｄｉｎｇ  ａｎｄ  Ｃｈａｎｎｅｌ  Ｅｓｔｉｍａｔｉｏｎ」（Ｊａｎｕａｒｙ  １９９８，  Ｍｕｎｉｃｈ）に開示されているように、幾つかのチャネル推定方式は、雑音分散が一定であると仮定している。この簡単な手法によって、通常の固定された雑音分散を前提として、フィルター係数を事前に計算することが可能になり、したがって、リアルタイム計算が回避される。
【００４７】
例えば、Ｃｈａｎｇ−Ｙｉ  Ｙａｎｇ及びＢｏｒ−ｓｅｎ  Ｃｈｅｎの論文「Ａｄａｐｔｉｖｅ  Ｃｈａｎｎｅｌ−Ｔｒａｃｋｉｎｇ  Ｍｅｔｈｏｄ  ａｎｄ  Ｅｑｕａｌｉｚａｔｉｏｎ  ｆｏｒ  ＭＣ−ＣＤＭＡ  Ｓｙｓｔｅｍｓ  ｏｖｅｒ  Ｒａｐｉｄｌｙ  Ｆａｄｉｎｇ  Ｃｈａｎｎｅｌ  ｕｎｄｅｒ  Ｃｏｌｏｒｅｄ  Ｎｏｉｓｅ」（Ｅｕｒａｓｉｐ  Ｊｏｕｒｎａｌ  ｏｎ  Ａｄｖａｎｃｅｄ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，  Ｖｏｌｕｍｅ  ２０１０）に開示されるような別の手法は、最小平均二乗誤差（ＭＭＳＥ）基準を最小にするフィルターを直接設計することである。この手法は、雑音が有色であり得ること、すなわち、分散がインデックスｋに依存することを考慮に入れる利点を有する。
【００４８】
フィルタリングをチャネルベクトルｈ＝（ｈｋ）の直接推定によって表す場合、フィルタリング（ｈの推定）は、行列Ａによる乗算に対応する。
【数４】
2019511169.tif 000006
ここで、ｇは、雑音フィルタリング前のパイロット位置において取得される粗チャネル推定値である。
【００４９】
次に、ＭＭＳＥ基準に従ってＡを定義することは、以下のようにＡを選択することに対応する。
【数５】
2019511169.tif 000007
【００５０】
対応する一般的なＭＭＳＥ解は、以下のものであることが知られている。
【数６】
2019511169.tif 000008
【００５１】
ここで、Σ（・）は、相関行列Σｘ＝Ｅ（ｘｘＨ）を表し、ｘＨはベクトルｘの複素共役を表す。
【００５２】
雑音及びチャネルｈが独立であると仮定することによって、以下の式が得られる。
【数７】
2019511169.tif 000009
【００５３】
チャネル相関行列Σｈは、システム設計、特にパイロット位置と、例えば遅延スプレッドのような幾つかのチャネル特性とに依存する。以下では、このチャネル相関行列が既知であると仮定する。例えば、従来から、チャネルの遅延スプレッドは、ＯＦＤＭシステム、ＳＣ−ＯＦＤＭシステム又はＳＣ−ＦＤＭＡシステムのサイクリックプレフィックス継続時間に等しいと仮定される。
【００５４】
異なるパイロットに影響を及ぼす雑音成分は無相関であると仮定することによって、雑音相関行列は、以下の式として得られる。
【数８】
2019511169.tif 000010
【００５５】
この手法は欠点を有する。すなわち、新たな雑音相関行列Σｎごとにフィルター係数を再計算しなければならず、特にフィルターサイズが大きい場合に、関連した複雑度が高い。フィルター効率は、このサイズに強く依存することを想起することにする。分散分布は、比較的高速に変化する可能性があり、過度に高速に変化すると、低複雑度でフィルター係数のリアルタイム計算が可能でなくなる可能性があることが起こり得る。本発明によれば、受信機Ｒｅｃは、定雑音分散σｂ２について、すなわち、雑音分散の変動の知識を有しないで計算された登録済みのフィルター係数を有し、少なくとも１つのフィルター係数は、他のフィルター係数と異なる。例えば、ＭＳＥ（平均二乗誤差）基準が用いられる場合、フィルター係数に対応する行列は以下に等しい。
【数９】
2019511169.tif 000011
【００５６】
ここで、Ｉは単位行列を表す。
【００５７】
σｂ２は、干渉プラス熱雑音に対応する全てのσｋ２よりも小さいものと仮定することにする。これは、例えば、σｂ２が熱雑音分散に対応する場合に当てはまる。その場合、分散は、以下のように記述することができる。
【数１０】
2019511169.tif 000012
【００５８】
σｂ２が熱雑音分散に対応する場合、
【数１１】
2019511169.tif 000013
値は干渉分散に対応する。そうでない場合、
【数１２】
2019511169.tif 000014
は、全体の分散σｋ２と、所定のフィルター行列Ａｂａｓｉｃの計算に用いられる所定の値σｂ２との間の差に単に対応する。
【００５９】
最適な行列Ａの代わりに行列Ａｂａｓｉｃを用いることは、Ａ＝Σｈ（Σｈ＋Σｎ）−１において、ΣｎについてＤｉａｇ（σｋ２）をＤｉａｇ（σｂ２）に置換することに対応する。その結果得られる係数には大きな相違が生じ、したがって、性能に大きな相違が生じる。
【００６０】
行列Ａは、以下のように書き換えることができる。
【数１３】
2019511169.tif 000015
ここで、
【数１４】
2019511169.tif 000016
である。
【００６１】
ここで、
【数１５】
2019511169.tif 000017
が非常に大きいとき、近似
【数１６】
2019511169.tif 000018
は、近似
【数１７】
2019511169.tif 000019
よりも決定的でないことに留意しなければならない。
【００６２】
その上、基本行列は、以下のように書き換えることができる。
【数１８】
2019511169.tif 000020
【００６３】
したがって、定雑音分散σｂ２について計算された登録済みのフィルター係数と、σｂ２及びσｋ２に従って計算された係数ωｋとを組み合わせたものである近似的に最適なフィルターＡａｐは、以下のように表すことができる。
【数１９】
2019511169.tif 000021
【００６４】
ここで、雑音相関行列Σｎは、上記数式の第１の部分Ａｂａｓｉｃでは近似されており、第２の部分Ｄｉａｇ（ωｋ）では近似されていないことに留意しなければならない。
【００６５】
次に、たとえ部分的にしかすぎない場合であっても、分散が変動していることが考慮される。これは、ＡをＡｂａｓｉｃに単に置換するときには該当しなかったものである。この定式化は、フィルターの入力であるベクトルｇを重み係数
【数２０】
2019511169.tif 000022
によって重み付けし、次に、Ａｂａｓｉｃフィルターを適用することによって解釈することができる。
【００６６】
重み係数ωｋを計算する変形形態が幾つか存在する。上記導出において、σｂ２は、Ａｂａｓｉｃフィルターの計算に用いられるパラメーターである。これは、通常の信号対雑音比の逆数、又は小さいがゼロではない正則化パラメーター、又はσｋ２の平均に比例するパラメーターのような他の任意のパラメーターとすることができる。熱雑音分散σｍ２及び干渉体分散
【数２１】
2019511169.tif 000023
が推定される場合、重み
【数２２】
2019511169.tif 000024
を用いることができる。
【００６７】
基本フィルター係数をａｉｋで示すと、これは、
【数２３】
2019511169.tif 000025
に等しい新たな重み付きフィルター係数を適用することと等価である。
【００６８】
無線インターフェース２０５は、図３に開示するような構成要素を備える。
【００６９】
図３は、受信機の無線インターフェースの構成要素のブロック図を開示している。
【００７０】
受信機Ｒｅｃによって受信された信号は、復調モジュール３００によって復調される。例えば、復調モジュール３００は、受信信号の離散フーリエ変換（ＤＦＴ）を実行する。
【００７１】
復調後、すなわち周波数領域では、この周波数領域における受信信号は、以下の式によって表すことができる。
【数２４】
2019511169.tif 000026
【００７２】
ここで、ｈｋは、インデックスｋのキャリアのチャネル応答であり、νｋは、同じ周波数における付加雑音である。例えば、項νｋは、ＡＷＧＮ（熱）雑音及び干渉体を加えるものである。この干渉体のために、付加雑音の分散は周波数に依存し、σｋ２によって示される。
【００７３】
復調信号は、パイロット抽出モジュール３０１に提供される。
【００７４】
パイロット抽出モジュール３０１は、パイロット位置において、対応する受信シンボルを抽出し、それらの受信シンボルを粗チャネル推定モジュール３０５に提供する。
【００７５】
粗チャネル推定モジュール３０５は、各パイロット位置ｋにおけるチャネル応答ｈｋの粗推定値ｇｋを、例えば以下のように求める。
【数２５】
2019511169.tif 000027
【００７６】
ここで、ｋは、パイロット位置インデックスであり、ｙｋは、対応する受信シンボルであり、ｐｋは、第ｋの位置における送信パイロット値である。
【００７７】
ここで、この推定は、雑音（熱雑音及び／又は干渉体）による影響を受け、パイロット位置における推定値しか提供しないことに留意しなければならない。フィルタリングの目的は、雑音を低減することと、パイロット位置以外の位置におけるチャネル推定の補間を行うこととである。提案された方式は、主として、雑音プラス干渉低減問題に対処する。
【００７８】
パイロット抽出モジュール３０１は、対応する受信シンボルを雑音電力推定モジュール３０２に提供する。
【００７９】
雑音電力推定モジュール３０２は、付加雑音の分散σｋ２を推定する。提示を簡潔にするために、この推定は完全であると仮定する。
【００８０】
粗チャネル推定は、以下のようにモデル化される。
【数２６】
2019511169.tif 000028
【００８１】
付加雑音
【数２７】
2019511169.tif 000029
は、例えば、熱雑音及び干渉を加えることである。以下では、一般性を失うことなく、Ｅ（｜ｐｋ｜２）＝１であると仮定する。
【００８２】
したがって、付加雑音ｎｋは、ｖｋと同じ分散σｋ２を有し、中心であると仮定される。すなわち、Ｅ（｜ｎｋ｜２）＝σｋ２であり、Ｅ（ｎｋ）＝０である。
【００８３】
雑音電力推定値は、重み係数計算モジュール３０３に提供される。
【００８４】
重み係数モジュール３０３は、既に開示しているように重み係数ωｋを求める。
【００８５】
重み係数は、重み付けモジュール３０６に転送され、この重み付けモジュールは、粗チャネル推定モジュール３０６の出力ｇｋを重み係数ωｋによって重み付けする。
【００８６】
重み付けモジュール３０６の出力は、フィルタリングモジュール３０７に提供され、このフィルタリングモジュールは、重み付き出力ωｋｇｋをフィルタリングするのにＡｂａｓｉｃ行列を用いる。
【００８７】
ここで、Ａｂａｓｉｃ行列は、補間を提供することができる、すなわち、より多くのサンプルを提供するために、２つのシンボル又はシンボルブロックの間で補間を行うことに留意しなければならない。
【００８８】
本発明の特定の実現態様によれば、フィルタリングモジュール３０７の出力は、補正モジュール３０８に提供される。
【００８９】
幾つかの分散σｋ２が大きい場合、幾つかの重み付け係数ωｋは非常に小さい場合があり、フィルタリングモジュール３０７によって提供される推定値は、推定される元の信号ｈｋよりもはるかに低いエネルギーを有する。補正モジュール３０８は、更なる補正ステップによってこれを補償し、これは重み付け係数ωｋに依存する。
【００９０】
正規化係数計算モジュール３０４は、重み付け係数ωｋを用いて、補正係数を求める。
【００９１】
例えば、インデックスｉの各出力に以下の式を乗算することができる。ここで、ｉは、Ａｂａｓｉｃ行列の出力のインデックスである。
【数２８】
2019511169.tif 000030
ここで、μｉは所定の値である。例えば、μｉは、固定値μに設定することができ、例えばμ＝１である。代替的に、μｉは、既知の基本フィルターの対応する値に等しく設定することができる。すなわち、
【数２９】
2019511169.tif 000031
である。
【００９２】
ここで、補間がフィルタリングモジュール３０７によって行われない場合、ｉはｋと同等であることに留意しなければならない。
【００９３】
最後に、最適なフィルターの以下の近似が得られる。
【数３０】
2019511169.tif 000032
【００９４】
補正モジュール３０８の出力
【数３１】
2019511169.tif 000033
は、チャネル値ｈｉの最終推定値に対応する。説明を簡単にするために、以下では、
【数３２】
2019511169.tif 000034
を    
【数３３】
2019511169.tif 000035
に表記変更する。
【００９５】
復調モジュール３００の出力ｙｋ、雑音電力推定モジュール３０２の出力σｋ２及び補正モジュール３０８の出力
【数３４】
2019511169.tif 000036
は、等化モジュール３０９に提供され、この等化モジュールは、復調信号の等化を行う。例えば、等化モジュール３０９は、以下のＭＭＳＥ等化を実施する。
【数３５】
2019511169.tif 000037
【００９６】
球面復号器として、又はターボ等化の場合に事前情報を有する等化器として、より高性能な等化をここで行うことができる。
【００９７】
本発明の特定の実現態様によれば、幾つかの基本フィルター行列は、事前に計算されて記憶される。これらの基本フィルターのうちの１つが、例えば、推定された雑音分散、例えばこれらの分散の平均、若しくはチャネル相関行列Σｈに関連した重要なパラメーターの推定値、例えばチャネル遅延スプレッド、又はそれらの双方に従って選択される。
【００９８】
フィルター選択モジュール３１０は、雑音電力推定モジュール３０２によって提供された分散から及び／又はチャネル相関行列Σｈに関連した重要なパラメーターの推定値、例えばチャネル遅延スプレッドに従って、フィルター係数モジュール３１１に記憶された複数の所定の基本行列の中から基本行列Ａｂａｓｉｃを選択する。
【００９９】
選択された基本行列は、正規化係数モジュール３０４及び基本フィルターモジュール３０７に提供される。
【０１００】
図４は、本発明による受信機によって実行されるアルゴリズムの一例を開示している。
【０１０１】
ステップＳ４００において、受信機Ｒｅｃによって受信された信号は、復調モジュール３００によって復調される。例えば、復調は、受信信号の離散フーリエ変換（ＤＦＴ）を実行する。
【０１０２】
復調後、すなわち周波数領域では、この周波数領域における受信信号は、以下の式によって表すことができる。
【数３６】
2019511169.tif 000038
【０１０３】
ここで、ｈｋは、インデックスｋのキャリアのチャネル応答であり、νｋは、同じ周波数における付加雑音である。例えば、項νｋは、ＡＷＧＮ（熱）雑音及び干渉体を加えるものである。この干渉体のために、付加雑音の分散は周波数に依存し、σｋ２によって示される。
【０１０４】
ステップＳ４０１において、パイロット抽出が行われる。
【０１０５】
このパイロット抽出によって、復調信号のパイロット位置において、対応する受信パイロットシンボルが抽出される。
【０１０６】
ステップＳ４０２において、抽出されたパイロットシンボルから、粗チャネル推定が行われる。
【０１０７】
粗チャネル推定は、各パイロット位置ｋにおけるチャネル応答ｈｋの粗推定値ｇｋを、例えば以下のように求める。
【数３７】
2019511169.tif 000039
【０１０８】
ここで、ｋは、パイロットインデックスであり、ｙｋは、対応する受信シンボルであり、ｐｋは、第ｋの位置における送信パイロット値である。
【０１０９】
ここで、この推定は、雑音（熱雑音及び／又は干渉体）による影響を受け、パイロット位置における推定値しか提供しないことに留意しなければならない。
【０１１０】
ステップＳ４０３において、雑音電力推定が、受信パイロットシンボルに対して行われる。
【０１１１】
雑音電力推定は、付加雑音の分散σｋ２を推定する。
【０１１２】
粗チャネル推定は、以下のようにモデル化される。
【数３８】
2019511169.tif 000040
【０１１３】
付加雑音
【数３９】
2019511169.tif 000041
は、例えば、熱雑音及び干渉を加えることである。以下では、一般性を失うことなく、Ｅ（｜ｐｋ｜２）＝１であると仮定する。
【０１１４】
したがって、付加雑音ｎｋは、ｖｋと同じ分散σｋ２を有し、中心であると仮定される。すなわち、Ｅ（｜ｎｋ｜２）＝σｋ２であり、Ｅ（ｎｋ）＝０である。
【０１１５】
ステップＳ４０４において、重み係数計算が雑音電力推定に基づいて行われる。
【０１１６】
重み係数ωｋは、既に開示したように求められる。
【０１１７】
ステップＳ４０６において、求められた重み係数は、粗チャネル推定の出力ｇｋを重み付けする。
【０１１８】
ステップＳ４０７において、フィルタリングが、粗チャネル推定の重み付き出力ωｋｇｋに対して実行される。このフィルタリングは、重み付き出力ωｋｇｋをフィルタリングするのにＡｂａｓｉｃ行列を用いる。
【０１１９】
ここで、Ａｂａｓｉｃ行列は、補間を提供することができる、すなわち、より多くのサンプルを提供するために、２つのシンボル又はシンボルブロックの間で補間を行うことに留意しなければならない。
【０１２０】
ステップＳ４０９において、補正係数が、重み付け係数ωｋを用いて求められる。
【０１２１】
幾つかの分散σｋ２が大きい場合、幾つかの重み付け係数ωｋは非常に小さい場合があり、フィルタリングモジュール３０７によって提供される推定値は、推定される元の信号ｈｋよりもはるかに低いエネルギーを有する。補正モジュール３０８は、更なる補正ステップによってこれを補償し、これは重み付け係数に依存する。
【０１２２】
例えば、インデックスｉの各出力に以下の式を乗算することができる。ここで、ｉは、Ａｂａｓｉｃ行列の出力のインデックスである。
【数４０】
2019511169.tif 000042
ここで、μｉは所定の値である。例えば、μｉは、固定値μに設定することができ、例えばμ＝１である。代替的に、μｉは、既知の基本フィルターの対応する値に等しく設定することができる。すなわち、
【数４１】
2019511169.tif 000043
である。
【０１２３】
ここで、補間がフィルタリングによって行われない場合、ｉはｋと同等であることに留意しなければならない。
【０１２４】
ステップＳ４０９において、補正が、フィルタリングモジュール３０７の出力に対して行われる。
【０１２５】
最後に、最適なフィルターの以下の近似が得られる。
【数４２】
2019511169.tif 000044
【０１２６】
補正モジュール３０８の出力
【数４３】
2019511169.tif 000045
は、チャネル値ｈｉの最終推定値に対応する。説明を簡単にするために、以下では、
【数４４】
2019511169.tif 000046
を    
【数４５】
2019511169.tif 000047
に表記変更する。
【０１２７】
ステップＳ４１０において、等化が、復調ステップの出力ｙｋと、雑音電力推定ステップの出力σｋ２と、補正ステップの出力
【数４６】
2019511169.tif 000048
とを用いて行われる。例えば、等化は、以下のＭＭＳＥ等化を実施する。
【数４７】
2019511169.tif 000049
【０１２８】
球面復号器として、又はターボ等化の場合に事前情報を有する等化器として、より高性能な等化をここで行うことができる。
【０１２９】
本発明の特定の実現態様によれば、幾つかの基本フィルター行列は、事前に計算されて記憶される。これらの基本フィルターのうちの１つが、ステップＳ４０５において、例えば、推定された雑音分散、例えばこれらの分散の平均、若しくはチャネル相関行列Σｈに関連した重要なパラメーターの推定値、例えばチャネル遅延スプレッド、又はそれらの双方に従って選択される。
【０１３０】
フィルター選択ステップＳ４０５は、雑音電力推定によって提供された分散から及び／又はチャネル相関行列Σｈに関連した重要なパラメーターの推定値、例えばチャネル遅延スプレッドに従って、フィルター係数モジュール３１１に記憶された複数の所定の基本行列の中から基本行列Ａｂａｓｉｃを選択する。
【０１３１】
図５は、本発明が受信機において実施されたときのパケットエラーレートの進展の一例を開示している。
【０１３２】
横軸はキャリア対雑音比を表し、縦軸はパケットエラーレートを表している。
【０１３３】
５０で示す曲線は、干渉体が存在しない場合に最適なチャネル推定が行われた一例を表している。
【０１３４】
５１で示す曲線は、干渉体が存在する場合に最適なチャネル推定が行われた一例を表している。
【０１３５】
５２で示す曲線は、重み付けがフィルタリングステップ４０７の前に行われていない場合の結果を表している。
【０１３６】
５３で示す曲線は、重み付けがフィルタリングステップ４０７の前に行われるとともに正規化ステップが実行された場合の結果を表している。
【０１３７】
その上、利得は、シミュレーション条件に依存することになる。１つのトーン及び−１０ｄＢのＣ／Ｉの場合、提案された方式を適用すると、これらのチャネル状態において、現時点における最先端技術に対して約０．５ｄＢの利得が提供される。その上、完全なチャネル推定と比較した損失は０．１ｄＢ未満であるので、これは、最適なフィルター計算と比較した損失が更に低くなることを意味している。
【０１３８】
本発明は、フィルター入力の重み付けが、ＭＳＥ基準について最適化された基本フィルターを考慮することによって行われる一例において開示されている。フィルター入力の重み付けは、これらの基本フィルターが、固定分散を明示的又は暗黙的に前提とすることによって求められる限り、上記方法が基本フィルターを計算するのに用いたものがいかなるものであっても適用することができる。例えば、周波数応答に基づくフィルター設計は、固定分散を暗黙的に前提としている。
【０１３９】
本発明は、単数又は複数の基本行列が事前に知られている一例において開示されていることに留意しなければならない。本発明は、基本行列がその時々において再計算され、そのレートが、全体の複雑度がこの計算によって大きな影響を受けないほど十分に低い場合にも適用可能であり、また、干渉体分散がこの低いレートよりも高速に変化すると予想される場合にも適用可能である。
【０１４０】
当然のことながら、本発明の範囲から逸脱することなく、上記で説明した本発明の実施形態に対して多くの変更を行うことができる。
【図１】
2019511169.tif 000053
【図２】
2019511169.tif 000054
【図３】
2019511169.tif 000055
【図４】
2019511169.tif 000056
【図５】
2019511169.tif 000057
【国際調査報告】
2019511169.tif 000058
2019511169.tif 000059
2019511169.tif 000060