(57)【要約】本発明は、２次元チャネル推定の概念を、複数の（ＮＴ）送信アンテナから受信側に作用するマルチキャリア変調送信信号を使用するマルチキャリア式多重入力多重出力通信システムに拡大する。送信信号は、パイロットシンボルが挿入された２次元データ系列を搬送すると仮定される。２段階チャネル推定処理の基礎として、複数の送信アンテナが、分割送信アンテナサブセットに分割される（小文字ｌ）。チャネル推定の第１段階において、作用するパイロット系列が、第１段階チャネル推定を行うことにより送信アンテナサブセット（小文字ｌ）に関連して分割され、送信の第１次元におけるチャネル応答の仮推定を得る。チャネル推定の第２段階において、作用するパイロット系列が、各送信アンテナサブセット（小文字ｌ）の各アンテナに対して第２段階チャネル推定を行うことにより送信アンテナサブセット（小文字ｌ）のアンテナに関連して分割され、チャネル応答の推定を得る。
2006515481.tif 000002
【特許請求の範囲】
【請求項１】
以下のステップを備える、複数（ＮＴ）の送信アンテナから作用し、パイロットシンボルが挿入された２次元データ系列を搬送する、マルチキャリア変調送信信号を使用する多重入力多重出力送信システムのための、２次元チャネル推定方法：
2006515481.tif 000116
作用するパイロット系列が、第１段階チャネル推定を行うことにより送信アンテナサブ
2006515481.tif 000117
るステップ、および
2006515481.tif 000118
連して分割され、チャネル応答の推定を得るステップ。
【請求項２】
第１段階チャネル推定はパイロットシンボルの２次元グリッドとして配置されたパイロット系列を使用して行われ、第１段階チャネル推定に使用されるパイロットシンボルは送信の第１次元にのみ依存し、第２段階チャネル推定に使用されるパイロットシンボルは送信の第２次元にのみ依存することを特徴とする請求項１に記載の方法。
【請求項３】
送信の第１次元と送信の第２次元においてパイロット系列が分離可能な最終形式として表現されることを特徴とする請求項２に記載の方法。
【請求項４】
パイロット系列が直交デザインから選択されることを特徴とする請求項３に記載の方法。
【請求項５】
直交デザインがウォルシュ系列および位相シフト系列を含むグループから選択されることを特徴とする請求項４に記載の方法。
【請求項６】
第１次元が周波数方向であることを特徴とする請求項１乃至５のいずれかに記載の方法。
【請求項７】
パイロット系列が以下の最終形式で表現されることを特徴とする請求項６に記載の方法：
2006515481.tif 000119
ここで、
2006515481.tif 000120
し、
2006515481.tif 000121
示し、
2006515481.tif 000122
示し、
2006515481.tif 000123
に含まれるアンテナの数であり、送信アンテナサブセットの数と送信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しいＮＴ＝ＮＴ１・ＮＴ２。
【請求項８】
パイロットの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値を与えるために、補間を用いない方法により、推定の次元におけるパイロットシンボルのグリッド位置の仮推定を行うことを特徴とする請求項６または７に記載の方法。
【請求項９】
パイロットの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値よりも大きな値を与えるために、補間を用いる方法により、推定の次元におけるデータ系列全ての仮推定を行うことを特徴とする請求項６または７に記載の方法。
【請求項１０】
第１段階チャネル推定は以下の式に従い、パイロットシンボルおよび受信したパイロットシンボルの値を用いて行われることを特徴とする請求項７乃至９のいずれかに記載の方法：
2006515481.tif 000124
ここで、
2006515481.tif 000125
【請求項１１】
第１段階チャネル推定が周波数領域で行われることを特徴とする請求項６乃至１０のいずれかに記載の方法。
【請求項１２】
第１段階チャネル推定が時間領域で行われることを特徴とする請求項６乃至１０のいずれかに記載の方法。
【請求項１３】
第２段階チャネル推定が周波数領域で行われることを特徴とする請求項１１または１２に記載の方法。
【請求項１４】
第２段階チャネル推定が時間領域で行われることを特徴とする請求項１１または１２に記載の方法。
【請求項１５】
第１次元が時間方向であることを特徴とする請求項１乃至５のいずれかに記載の方法。
【請求項１６】
パイロット系列が以下の最終形式で表現されることを特徴とする請求項１５に記載の方法：
2006515481.tif 000126
ここで、
2006515481.tif 000127
し、
2006515481.tif 000128
示し、
2006515481.tif 000129
示し、
2006515481.tif 000130
に含まれるアンテナの数であり、送信アンテナサブセットの数と送信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しいＮＴ＝ＮＴ１・ＮＴ２。
【請求項１７】
パイロットの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値を与えるために、補間を用いない方法により、推定の次元におけるパイロットシンボルのグリッド位置の仮推定を行うことを特徴とする請求項１５または１６に記載の方法。
【請求項１８】
パイロットの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値よりも大きな値を与えるために、補間を用いる方法により、推定の次元におけるデータ系列全ての仮推定を行うことを特徴とする請求項１５または１６に記載の方法。
【請求項１９】
第１段階チャネル推定は以下の式に従い、パイロットシンボルおよび受信したパイロットシンボルの値を用いて行われることを特徴とする請求項１６乃至１８のいずれかに記載の方法：
2006515481.tif 000131
ここで、
2006515481.tif 000132
【請求項２０】
第１段階チャネル推定が周波数領域で行われることを特徴とする請求項１５乃至１９のいずれかに記載の方法。
【請求項２１】
第１段階チャネル推定が時間領域で行われることを特徴とする請求項１５乃至１９のいずれかに記載の方法。
【請求項２２】
第２段階チャネル推定が周波数領域で行われることを特徴とする請求項２０または２１に記載の方法。
【請求項２３】
第２段階チャネル推定が時間領域で行われることを特徴とする請求項２０または２１に記載の方法。
【請求項２４】
周波数領域における第１段階チャネル推定が、受信したパイロット
2006515481.tif 000133
2006515481.tif 000134
ることを特徴とする請求項１１または２０に記載の方法。
【請求項２５】
周波数領域における第２段階チャネル推定が、送信アンテナサブセ
2006515481.tif 000135
または２２に記載の方法。
【請求項２６】
周波数領域における第２段階チャネル推定が、平滑化フィルタリングのための時間遅延をもたらす有限インパルス応答フィルタリングにより実行されることを特徴とする請求項２４または２５に記載の方法。
【請求項２７】
周波数領域における第２段階チャネル推定が、予測フィルタリングのための時間遅延をもたらさない有限インパルス応答フィルタリングにより実行されることを特徴とする請求項２４または２５に記載の方法。
【請求項２８】
有限インパルス応答フィルタリングがウィーナ型フィルタリングであることを特徴とする請求項２４乃至２７のいずれかに記載の方法。
【請求項２９】
時間領域における第１段階チャネル推定が、受信したパイロットシンボルを逆フーリエ変換し、逆フーリエ変換の結果を有限インパルス応答フィルタリングし、ゼロパディングし、ゼロパディングの結果をフーリエ変換し、送信アンテナサブセッ
2006515481.tif 000136
ることを特徴とする請求項１２または２１に記載の方法。
【請求項３０】
時間領域における第２段階チャネル推定が、送信アンテナサブセッ
2006515481.tif 000137
ンとして、逆フーリエ変換の結果を有限インパルス応答フィルタリングし、ゼロパディン
2006515481.tif 000138
することにより実行されることを特徴とする請求項１４または２３に記載の方法。
【請求項３１】
フーリエ変換の前に有限インパルスフィルタリングの結果を補間し、２次元データ系列のサブキャリアの数に対応するように時間領域の仮推定の数を調整するステップをさらに備えることを特徴とする請求項２９または３０に記載の方法。
【請求項３２】
受信したパイロットシンボルの数に対する時間領域における仮推定の数の比率が、送信の第１次元におけるパイロットの間隔に対応することを特徴とする請求項２９または３０に記載の方法。
【請求項３３】
第２段階チャネル推定が有限インパルス応答フィルタリングの結果に基づき実行されることを特徴とする請求項２９乃至３２に記載の方法。
【請求項３４】
基地局とそれに対応するアンテナアレイが複数の送信アンテナを構
2006515481.tif 000139
アンテナに関係付けて定義される周波数再利用ファクタを伴う携帯電話通信システムにおいて適用されることを特徴とする請求項１乃至３３のいずれかに記載の方法。
【請求項３５】
複数の送信アンテナ（ＮＴ）から作用し、パイロットシンボルが挿入された２次元データ系列を搬送する、マルチキャリア変調送信信号を使用する多重入力多重出力送信システムのための、以下の要素を備えるチャネル推定装置：
2006515481.tif 000140
するメモリ、
作用するパイロット系列を、第１段階チャネル推定を行うことにより送信アンテナサブ
2006515481.tif 000141
ための少なくとも１つの第１段階推定手段、および
2006515481.tif 000142
連して分割し、チャネル応答の推定を得るための少なくとも１つの第２段階推定手段。
【請求項３６】
第１段階推定手段は、パイロットシンボルの２次元グリッドとして配置されたパイロット系列を使用して第１段階チャネル推定を行い、第１段階チャネル推定に使用されるパイロットシンボルは送信の第１次元にのみ依存し、第２段階チャネル推定に使用されるパイロットシンボルは送信の第２次元にのみ依存することを特徴とする請求項１に記載のチャネル推定装置。
【請求項３７】
第１段階推定手段は、パイロット系列を用いる第１段階チャネル推定を行い、当該パイロット系列は送信の第１次元と送信の第２次元においてパイロット系列が分離可能な最終形式として表現されることを特徴とする請求項３６に記載のチャネル推定装置。
【請求項３８】
パイロット系列が直交デザインから選択されることを特徴とする請求項３７に記載のチャネル推定装置。
【請求項３９】
第１段階推定手段は、直交デザインがウォルシュ系列もしくは位相シフト系列を用いる第１段階チャネル推定を行うことを特徴とする請求項３８に記載のチャネル推定装置。
【請求項４０】
第１段階推定手段は、第１次元が周波数方向である第１段階チャネル推定を行うことを特徴とする請求項３５乃至３９のいずれかに記載のチャネル推定装置。
【請求項４１】
第１段階推定手段は、以下の最終形式で表現されるパイロット系列を用いる第１段階チャネル推定を行うことを特徴とする請求項４０に記載のチャネル推定装置：
2006515481.tif 000143
ここで、
2006515481.tif 000144
し、
2006515481.tif 000145
示し、
2006515481.tif 000146
示し、
2006515481.tif 000147
に含まれるアンテナの数であり、送信アンテナサブセットの数と送信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しいＮＴ＝ＮＴ１・ＮＴ２。
【請求項４２】
第１段階推定手段および／もしくは第２段階推定手段は、補間を用いない方法により、推定の次元におけるパイロットシンボルのグリッド位置の仮推定を行い、パイロットの間に値を与えるチャネル推定を行うことを特徴とする請求項４０または４１に記載のチャネル推定装置。
【請求項４３】
第１段階推定手段および／もしくは第２段階推定手段は、補間を用いる方法により、推定の次元におけるデータ系列全ての仮推定を行い、パイロットの間に推定の値よりも大きな値を与えるチャネル推定を行うことを特徴とする請求項４０または４１に記載のチャネル推定装置。
【請求項４４】
第１段階推定手段は、パイロットシンボルおよび受信したパイロットシンボルの値を用いて以下の式に従い第１段階チャネル推定を行うことを特徴とする請求項４１乃至４３のいずれかに記載のチャネル推定装置：
2006515481.tif 000148
ここで、
2006515481.tif 000149
【請求項４５】
第１段階推定手段は、周波数領域で第１段階チャネル推定を行うことを特徴とする請求項４０乃至４４のいずれかに記載のチャネル推定装置。
【請求項４６】
第１段階推定手段は、時間領域で第１段階チャネル推定を行うことを特徴とする請求項４０乃至４４のいずれかに記載のチャネル推定装置。
【請求項４７】
第２段階推定手段は、周波数領域で第２段階チャネル推定を行うことを特徴とする請求項４５または４６に記載のチャネル推定装置。
【請求項４８】
第２段階推定手段は、時間領域で第２段階チャネル推定を行うことを特徴とする請求項４５または４６に記載のチャネル推定装置。
【請求項４９】
第１段階推定手段は、時間方向で第１段階チャネル推定を行うことを特徴とする請求項３５乃至３９のいずれかに記載のチャネル推定装置。
【請求項５０】
第１段階推定手段は、以下の最終形式で表現されるパイロット系列を用いて第１段階チャネル推定を行うことを特徴とする請求項４９に記載のチャネル推定装置：
2006515481.tif 000150
ここで、
2006515481.tif 000151
し、
2006515481.tif 000152
示し、
2006515481.tif 000153
示し、
2006515481.tif 000154
に含まれるアンテナの数であり、送信アンテナサブセットの数と送信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しいＮＴ＝ＮＴ１・ＮＴ２。
【請求項５１】
パイロットの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値を与えるために、補間を用いない方法により、推定の次元におけるパイロットシンボルのグリッド位置の仮推定を行うことを特徴とする請求項４９または５０に記載の方法。
【請求項５２】
パイロントの間に第１段階チャネル推定および／もしくは第２段階チャネル推定の値よりも大きな値を与えるために、補間を用いる方法により、推定の次元におけるデータ系列全ての仮推定を行うことを特徴とする請求項４９または５０に記載の方法。
【請求項５３】
第１段階推定手段は、以下の式に従い、パイロットシンボルおよび受信したパイロットシンボルの値を用いて第１段階チャネル推定を行うことを特徴とする請求項５０乃至５２のいずれかに記載のチャネル推定装置：
2006515481.tif 000155
ここで、
2006515481.tif 000156
【請求項５４】
第１段階推定手段は、周波数領域で第１段階チャネル推定を行うことを特徴とする請求項４９乃至５３のいずれかに記載のチャネル推定装置。
【請求項５５】
第１段階推定手段は、時間領域で第１段階チャネル推定を行うことを特徴とする請求項４９乃至５３のいずれかに記載のチャネル推定装置。
【請求項５６】
第２段階推定手段は、周波数領域で第２段階チャネル推定を行うことを特徴とする請求項５４または５５に記載のチャネル推定装置。
【請求項５７】
第２段階推定手段は、時間領域で第２段階チャネル推定を行うことを特徴とする請求項５４または５５に記載のチャネル推定装置。
【請求項５８】
第１段階推定手段が、周波数領域において受信したパイロットシン
2006515481.tif 000157
応答フィルタを備えることを特徴とする請求項４５または５４に記載のチャネル推定装置。
【請求項５９】
第２段階推定手段が、周波数領域において送信アンテナサブセット
2006515481.tif 000158
する請求項４７または５６に記載のチャネル推定装置。
【請求項６０】
第２段階推定手段が、周波数領域において平滑化フィルタリングのための時間遅延をもたらす有限インパルス応答フィルタリングを行うためのバッファを備えることを特徴とする請求項５８または５９に記載のチャネル推定装置。
【請求項６１】
第１有限インパルス応答フィルタおよび／もしくは第２有限インパルス応答フィルタがウィーナ型フィルタであることを特徴とする請求項５８乃至６０のいずれかに記載のチャネル推定装置。
【請求項６２】
第１段階推定手段が以下の要素を備えることを特徴とする請求項４６または５５に記載のチャネル推定装置：
受信したパイロットシンボルを逆フーリエ変換する第１逆フーリエ変換部、
オプションとして、逆フーリエ変換の結果を有限インパルス応答フィルタリングする第３有限インパルス応答フィルタ、
ゼロパディング部、そして
2006515481.tif 000159
【請求項６３】
第２段階推定手段が以下の要素を備えることを特徴とする請求項４８または５７に記載のチャネル推定装置：
2006515481.tif 000160
ーリエ変換する第２逆フーリエ変換部、
オプションとして、逆フーリエ変換の結果を有限インパルス応答フィルタリングする第４有限インパルス応答フィルタ、
ゼロパディング部、そして
ゼロパディングの結果をフーリエ変換し、時間領域におけるチャネル応答の推定値
2006515481.tif 000161
【請求項６４】
フーリエ変換の前に有限インパルスフィルタリングの結果を補間し、２次元データ系列のサブキャリアの数に対応するように時間領域の仮推定の数を調整する補間フィルタをさらに備えることを特徴とする請求項６２に記載のチャネル推定装置。
【請求項６５】
受信したパイロットシンボルの数に対する時間領域における仮推定の数の比率が、送信の第１次元におけるパイロットの間隔に対応することを特徴とする請求項６４に記載のチャネル推定装置。
【請求項６６】
第２段階推定手段は第３有限インパルス応答フィルタに接続され、有限インパルス応答フィルタリングの結果に基づき第２段階チャネル推定を行うことを特徴とする請求項６２乃至６５に記載のチャネル推定装置。
【請求項６７】
基地局とそれに対応するアンテナアレイが複数の送信アンテナを構
2006515481.tif 000162
アンテナに関係付けて定義される周波数再利用ファクタを伴う携帯電話通信システムにおいて適用されることを特徴とする請求項３５乃至６６のいずれかに記載のチャネル推定装置。
【請求項６８】
プロダクトが、多重入力多重出力送信チャネルを２次元において推定するためのチャネル推定装置のプロセッサにおいて実行された時に、請求項１乃至３４のいずれかに記載の各ステップを実行させるソフトウェアコード部分を備える、多重入力多重出力送信チャネルを２次元において推定するためのチャネル推定装置の内部メモリに直接的にロード可能なコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、多重入力多重出力通信システムのためのチャネル推定に関し、特にマルチキャリア変調スキームを用いた多重入力多重出力通信システムのための方法および関連するチャネル推定装置に関する。
【背景技術】
【０００２】  無線通信システムにおけるコヒーレント送信技術の使用には、チャネル推定として公知である移動無線チャネルのトラッキングが必要である。多重送信アンテナから送信された信号は相互干渉として観測されるため、多重入力多重出力ＭＩＭＯ通信システムのためのチャネル推定は、単一送信アンテナの場合とは異なる。マルチキャリア変調スキームを用いた多重入力多重出力ＭＩＭＯ通信システムにおいては、マルチキャリア復調後の受信信号は、時間および周波数などの２次元において相関がとられることが典型的である。以下に、マルチキャリア変調スキームの典型例の１つとして、直交周波数分割多重ＯＦＤＭを挙げる。なぜならば、直交周波数分割多重ＯＦＤＭおよびその変形物はマルチキャリア変調スキームにおいて最もよく普及しているためである。【０００３】  多重入力多重出力ＭＩＭＯ通信システムとして公知の多重送受信アンテナを採用した通信システムを、直交周波数分割多重ＯＦＤＭと併用して、移動無線システムの通信容量および品質を向上することが可能である。ナギブ（Ａ．Ｎａｇｕｉｂ）、シシャドリ（Ｎ．Ｓｅｓｈａｄｒｉ）およびカルターバンク（Ａ．Ｃａｌｄｅｒｂａｎｋ）の『無線高速データ通信のための時空間符号化および信号処理』（ＩＥＥＥ信号処理誌、７６〜９２頁、２０００年５月）において記載されている時空間符号などの多重送信アンテナを用いた直交周波数分割多重ＯＦＤＭ通信システム、または特殊な多重方式においては、異なる信号が異なる送信アンテナから同時に送信される。従って、受信信号はこれらの信号が重畳されたものとなり、これがチャネル推定における課題を意味する。チャネルパラメータは、時空間符号を使用した場合にはダイバーシチの結合のために必要となり、また、空間多重方式を使用した場合には重畳信号の分割のために必要となる。【０００４】  このようなチャネル推定のための手法は、リー（Ｙ．Ｌｉ）、シシャドリ（Ｎ．Ｓｅｓｈａｄｒｉ）およびアリャウィシタクル（Ｓ．Ａｒｉｙａｖｉｓｉｔａｋｕｌ）の『移動無線チャネルにおける送信ダイバーシチを用いたＯＦＤＭシステムのためのチャネル推定』（ＩＥＥＥ通信選択エリアジャーナル、１７巻、４６１〜４７０頁、１９９９年３月）、およびリー（Ｙ．Ｌｉ）の『多重送信アンテナを用いたＯＦＤＭシステムのための簡略化したチャネル推定』（ＩＥＥＥ無線通信議事録、１巻、６７〜７５頁、２００２年１月）の、特に離散フーリエ変換ＤＦＴに基づく多重送信アンテナを用いた直交周波数分割多重ＯＦＤＭのためのチャネル推定スキームにおいて記載されている。【０００５】  さらに、ＯＦＤＭ−ＭＩＭＯシステムのための最小二乗ＬＳおよび最小平均二乗誤差ＭＭＳＥ規準に基づく推定装置は、ゴング（Ｙ．Ｇｏｎｇ）およびラティエフ（Ｋ．Ｌｅｔａｉｅｆ）の『時空間符号化広帯域ＯＦＤＭシステムのための低ランクチャネル推定』（ＩＥＥＥ車両技術会議議事録（ＶＴＣ２００１年秋）、アトランチックシティー、米国、７２２〜７７６頁、２００１年）において体系的に演繹されている。関連する解決法においては、公知のパイロットＯＦＤＭシンボルに、Ｌデータを含むＯＦＤＭシンボルが続く１次元手法が取り上げられている。無線���内情報通信網ＷＬＡＮなどの構内システムのように、ＬＯＦＤＭシンボル中にチャネルが大きく変動しない準静的環境において、このスキームは適用可能である。【０００６】  ある程度の移動性を確保するために、リー（Ｙ．Ｌｉ）、シシャドリ（Ｎ．Ｓｅｓｈａｄｒｉ）およびアリャウィシタクル（Ｓ．Ａｒｉｙａｖｉｓｉｔａｋｕｌ）の『移動無線チャネルにおける送信ダイバーシチを用いたＯＦＤＭシステムのためのチャネル推定』（ＩＥＥＥ通信選択エリアジャーナル、１７巻、４６１〜４７０頁、１９９９年３月）、およびリー（Ｙ．Ｌｉ）の『多重送信アンテナを用いたＯＦＤＭシステムのための簡略化したチャネル推定』（ＩＥＥＥ無線通信議事録、１巻、６７〜７５頁、２００２年１月）において提案されているように、Ｌデータを含むＯＦＤＭシンボルの受信中に、受信側で判定指向型チャネル推定に切り替えてもよい。しかしながら、先行判定データシンボルをパイロットシンボルとして使用する判定指向型チャネル推定は、パイロットのみに依存するチャネル推定スキームよりも大幅に複雑となる。【０００７】  さらに、１つの送信アンテナを用いるＯＦＤＭに基づくシステムの場合、時間および周波数においてサンプリング定理を満たす分散パイロットグリッドを利用した２次元チャネル推定を採用可能である。パイロットシンボル援用チャネル推定ＰＡＣＥにおいては、公知のパイロットシンボルをデータストリームに多重化する。補間を使用して情報伝達シンボルのチャネル推定を得る。単一キャリアシステムのためのＰＡＣＥは、カバーズ（Ｊ．Ｋ．Ｃａｖｅｒｓ）の『レイリーフェージングチャネルに対するパイロットシンボル援用変調の分析』（ＩＥＥＥ車両技術議事録、ＶＴ−４０巻、６８６〜６９３頁、１９９１年１１月）において紹介された。ニルソン（Ｒ．Ｎｉｌｓｏｎ）、エドフォース（Ｏ．Ｅｄｆｏｒｓ）、サンデル（Ｍ．Ｓａｎｄｅｌｌ）およびボージェソン（Ｐ．Ｂｏｅｒｊｅｓｓｏｎ）の『ＯＦＤＭに対する２次元パイロットシンボル援用変調の分析』（ＩＥＥＥ個人無線通信国際会議議事録（ＩＣＰＷＣ’９７）、ムンバイ（ボンベイ）、インド、７１〜７４頁、１９９７年）、およびヘーアー（Ｐ．Ｈｏｅｈｅｒ）、カイサー（Ｓ．Ｋａｉｓｅｒ）およびロバートソン（Ｐ．Ｒｏｂｅｒｔｓｏｎ）の『ウィーナフィルタリングによる２次元パイロットシンボル援用チャネル推定』（ＩＥＥＥ音響、音声、信号処理国際会議議事録（ＩＣＡＳＳＰ’９７）、ミュンヘン、ドイツ１８４５〜１８４８頁、１９９７年）において、２次元２Ｄフィルタリングアルゴリズムがパイロットシンボル援用チャネル推定ＰＡＣＥに対し提案されている。しかしながら、一般的に、このような２Ｄ推定装置構造は実際の実施にあたってはあまりに複雑である。【０００８】  複雑さを軽減するために、ヘーアー（Ｐ．Ｈｏｅｈｅｒ）、カイサー（Ｓ．Ｋａｉｓｅｒ）およびロバートソン（Ｐ  Ｒｏｂｅｒｔｓｏｎ）の『時間および周波数におけるパイロットシンボル援用チャネル推定』（ＩＥＥＥ世界電気通信会議（ＧＬＯＢＥＣＯＭ’９７）と連携した通信理論小会議（ＣＴＭＣ）議事録、フェニックス、米国、９０〜９６頁、１９９７年）において、時間および周波数相関の使用を��割することが提案された。２重１次元（２×１Ｄ）パイロットシンボル援用チャネル推定ＰＡＣＥと呼ばれるこの併用スキームでは、周波数方向において１つおよび時間方向において１つの別々のウィーナフィルタが使用される。【０００９】  計算複雑性を軽減する別の手法は、チャネル出力を僅かな変換係数に集中させる変換に基づいている。離散フーリエ変換ＤＦＴに基づく推定装置の場合、ＦＦＴの形式をとる計算上効率的な変換が存在する、またＤＦＴに基づいた補間が簡単であるという利点を有する。リー（Ｙ．Ｌｉ）の『無線システムにおけるＯＦＤＭに対するパイロットシンボル援用チャネル推定』（ＩＥＥＥ車両技術議事録、４９巻、１２０７〜１２１５頁、２０００年７月）において、単一アンテナの場合には２次元ＦＦＴを使用することで、離散フーリエ変換ＤＦＴに基づくパイロットシンボル援用チャネル推定ＰＡＣＥによる手法を、２次元パイロットシンボル援用チャネル推定ＰＡＣＥに拡大している。【００１０】  しかしながら、分散パイロットグリッドを利用した２次元チャネル推定の手法を多重入力多重出力ＭＩＭＯ通信システムに拡大する際の主たる問題として、一定数のパイロットに分割可能な送信アンテナ数の限界がある。それぞれＱ個のタップを有するＮＴ個のチャネルインパルス応答（ＣＩＲ）を推定するのに必要とされるパイロットシンボル最小数
2006515481.tif 000003
号化広帯域ＯＦＤＭシステムのための低ランクチャネル推定』（ＩＥＥＥ車両技術会議議事録（ＶＴＣ２００１年秋）、アトランチックシティー、米国、７２２〜７７６頁、２００１年）において、以下のようになると示されている。
2006515481.tif 000004
【００１１】  しかしながら、これは送信アンテナ数ＮＴが増加するのに伴い、チャネル推定に必要とされるパイロット数も増加することを意味する。
【発明の開示】
【発明が解決しようとする課題】
【００１２】  上記に鑑みて、本発明の目的は２次元チャネル推定の概念をＭＩＭＯシステムに拡大することにある。
【課題を解決するための手段】
【００１３】  本発明によると上記目的は、複数の送信アンテナから作用し、パイロットシンボルが挿入された２次元データ系列を搬送する、マルチキャリア変調送信信号を使用する多重入力多重出力送信システムのための、２次元チャネル推定方法により達成される。第１ステップにおいては、複数の送信アンテナが、分割送信アンテナのサブセットに分割される。第２ステップにおいては、作用するパイロット系列が、第１段階チャネル推定を行うことにより送信アンテナサブセットに関連して分割され、送信の第１次元におけるチャネル応答の仮推定を得る。第３ステップにおいては、作用するパイロット系列が、各送信アンテナサブセットの各アンテナに対して第２段階チャネル推定を行うことにより送信アンテナサブセットのアンテナに関連して分割され、チャネル応答の推定を得る。【００１４】  本発明の重要な利点は、チャネル推定における柔軟性の向上である。なぜならば、時間および周波数方向にける重畳送信信号の分割タスクを分けることにより、パイロットシンボルのより効率的な使用が可能となるからである。よって、必要とされるパイロットシンボル数を減少できる、または能力を向上することが可能である。【００１５】  上記に鑑みて、本発明の別の重要な利点は、２段階チャネル推定手法の適用により一定数のパイロットシンボルにより推定できる送信アンテナ数が増加することである。【００１６】  本発明のさらに別の重要な利点は、２段階チャネル推定手法により、ドップラー周波数が高い場合においてもチャネル変動のトラッキングが可能となることである。これはモバイルユーザの高速度を支える必要条件であり、よって、移動に真に適した多重入力多重出力ＭＩＭＯ通信システムが実現される。【００１７】  本発明の好適な実施形態によると、第１段階チャネル推定はパイロットシンボルの２次元グリッドとして配置されたパイロット系列を使用して行われ、第１段階チャネル推定に使用されるパイロットシンボルは送信の第１次元にのみ依存し、第２段階チャネル推定に使用されるパイロットシンボルは送信の第２次元にのみ依存する。送信の第１次元および送信の第２次元においてパイロット系列を分割するために、パイロット系列が積の形で表されることが好ましい。【００１８】  本発明のこの好適な実施形態の利点は、分散パイロットグリッドの利用により、パイロットシンボルが効率よく使用できることにある。さらに、時間および周波数におけるパイロット間隔をチャネル特性が最悪の場合に合わせることにより、従来の１次元スキームに対して、より高い移動速度をサポートすることが可能である。【００１９】  本発明の別の好適な実施形態によると、１の値を有するパイロット間隔に対し、推定の次元におけるパイロットシンボルグリッド位置に関連して仮推定を得ることにより、第１段階チャネル推定および／または第２段階チャネル推定が補間なしで達成される。また、１よりも大きな値を有するパイロット間隔に対し、推定の次元においてすべてのデータ系列グリッド位置に対する仮推定を得ることにより、第１段階チャネル推定および／または第２段階チャネル推定が補間ありで達成される。【００２０】  この好適な実施形態の重要な利点は、異なる２��元のパイロットグリッドを柔軟にサポートできることにある。言い換えると、本発明は、周波数および時間方向の両方におけるいずれのタイプのパイロット間隔、および適切な補間技術を適用することにより、柔軟に適用することが可能である。【００２１】  本発明のさらなる好適な実施形態は、周波数方向または時間方向などの第１チャネル推定段階および第２チャネル推定段階に対して送信の第１次元を選択することに関し、さらに周波数領域チャネル推定または時間領域チャネル推定などのチャネル推定領域の選択に関する。ここで、本発明によると、送信次元およびチャネル推定領域のどのような組み合わせもサポートされる。【００２２】  送信次元およびチャネル推定領域が自由に選択可能であることは、本発明に係るチャネル推定手法の柔軟性を説明するもう１つの理由である。これにより、マルチキャリアに関する送信パラメータの最適な考察、選択されたパイロットグリッド構造、さらに計算上最適なチャネル推定技術の適用が可能となる。【００２３】  本発明の別の好適な実施形態によると、基地局および関連するアンテナアレイが複数の送信アンテナを形成するように、また、送信アンテナサブセットおよび関連する送信アンテナがこの複数の送信アンテナに関連して定義されるように、チャネル推定手法が周波数再利用係数が１であるセルラー通信システムに適用される。【００２４】  本発明のこの好適な実施形態の重要な利点は、上記に概説した２段階チャネル推定技術が分散したアンテナに適用できることにある。特に、これによって、モバイルユーザがセル境界においてローミングしている場合に対処可能となる。データを含むシンボルは、チャネルコードまたは拡散を用いて干渉から保護可能であるが、パイロットシンボルに対してはこれが不可能である。本発明によると、セルラー通信システムのセルに関連してサブセットを適切に定義することにより、これが可能となる。【００２５】  本発明のさらに別の好適な実施形態によると、２次元において多重入力多重出力送信チャネルを推定するための、チャネル推定装置の内部メモリに直接ロード可能なコンピュータプログラム製品であって、チャネル推定装置のプロセッサ上で製品が作動する時に本発明に係る２次元チャネル推定方法のステップを実行するソフトウエアコード部を含むコンピュータプログラム製品が提供される。【００２６】  よって、本発明はコンピュータまたはプロセッサシステムにおける進歩性を有する方法ステップの実施を実現するためにも提供される。結論として、このような実施により、２次元において多重入力多重出力送信チャネルを推定するチャネル推定装置などに構成されるコンピュータシステム、またはより特定的には、プロセッサと併用されるコンピュータプログラム製品が提供される。【００２７】  本発明の機能を定義するプログラムは、コンピュータ／プロセッサに対し、プロセッサまたはコンピュータＩ／Ｏ装置により読み出し可能なＲＯＭまたはＣＤ−ＲＯＭディスクなどの読み取り専用メモリ装置といった書き込み不可能な記録媒体に永久的に記録された情報、フロッピーディスクおよびハードライブなどの書き込み可能な記録媒体に記録された情報、または構内通信網、電話網、および／またはインターネットやその他のインターフェース装置などの通信媒体を介してコンピュータ／プロセッサに伝達される情報を含む数多くの形体で提供されるが、これらに限定されるものではない。このような媒体は、進歩性を有する概念を実施するプロセッサにより読み取り可能な指示を伝達する場合、本発叨の別の実施形態を表していると理解される。
【発明を実施するための最良の形態】
【００２８】  以下、本発明の最良の形態および好適な実施形態について図面を参照しながら説明する。まず、本発明のよりよい理解のために、差分マルチ長送信ダイバーシチおよび関連するダイバーシチ受信の基礎をなす基本的考察を説明する。【００２９】基本的考察【００３０】  本発明は、パイロットグリッドの挿入に基づくマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムのためのパイロットシンボル援用チャネル推定ＰＡＣＥを取り上げる。典型的に、マルチキャリア式多重入力多重出力ＭＩＭＯ通信システムに対しては、パイロットシンボル援用チャネル推定ＰＡＣＥ（以下、ＰＡＣＥともいう）を周波数方向におけるサブキャリアに、また、ＯＦＤＭシンボルなどのいくつかのマルチキャリア送信シンボルを介して２次元に適用され、結果として２次元の２ＤＰＡＣＥが得られる。【００３１】  以下に説明される本発明は、特定の種類のマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムに限定されるものではなく、直交周波数分割多重ＯＦＤＭ、離散マルチトーン送信ＤＭＴ、フィルタリングマルチトーン送信ＦＭＴ、または双直交周波数分割多重方式ＢＦＤＭなどに適用してもよい。【００３２】  別のマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムとして、マルチキャリア変調に加えて、周波数および／または時間方向への拡散が導入されているマルチキャリア符号分割多重アクセスＭＣ−ＣＤＭＡがある。さらに別のマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムとして、可変拡散率を用いたマルチキャリア符号分割多重アクセスＭＣ−ＣＤＭＡ、すなわち可変拡散率直交周波数・符号分割多重アクセスＶＳＦ−ＯＦＣＤＭがある。【００３３】  以上のことから、２次元において相関がとられる送信信号により作動するすべてのマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムに、パイロットシンボル援用チャネル推定ＰＡＣＥを適用してもよいことが理解される。よって、これらすべてのマルチキャリア式多重入力多重出力ＭＩＭＯ通信システムは、以下に説明されるように本発明の異なる実施形態を採用してもよい。【００３４】システムモデル【００３５】  図１は、本発明の基礎をなすシステムモデルを説明するためのＯＦＤＭに基づく多重入力多重出力ＭＩＭＯ通信システム（以下、ＯＦＤＭシステムともいう）の概略図を示す。さらに、図２は、図１に示したＯＦＤＭに基づく多重入力多重出力ＭＩＭＯ通信システムに対応する、ＯＦＤＭ変調および復調を説明する概略図を示す。【００３６】  図１に示すように、考察されるＯＦＤＭに基づくＭＩＭＯシステムにおいては、各送信アンテナに１つのＯＦＤＭ変調器が使用されているとともに、ＯＦＤＭ復調は各受信アンテナ毎に独立して行われる。典型的には、いずれのマルチキャリア変調スキームにおいても、信号ストリームはＮｃ個の並列したサブストリームに分割される。以下においてＯＦＤＭシンボルと記載される、ｌ番目のシンボルブロックのｉ番目のサブストリーム（以下、サブキャリアともいう）はＸｌ，ｉと表記する。ＮＦＦＴポイント逆ＤＦＴが各ブロックに対して行われ、続いてＮＧＩ個のサンプルを有するガードインターバルを挿入してｘｌ，ｎを得る。Ｄ／Ａ変換の後���信号ｘ（ｔ）は、応答ｈ（ｔ，τ）の移動無線チャネルを介して送信される。【００３７】  図１に示すように、受信アンテナνにおける受信信号は、ＮＴ個の送信アンテナからの重畳信号で構成される。同調が完全であると仮定すると、受信アンテナｖにおいてサンプリング時点ｔ＝［ｎ＋ｌＮｓｙｍ］Ｔｓｐｌで作用する等価ベースバンドシステムの受信信号は、以下の形式で表現される。
2006515481.tif 000005
ここで、ｘ（μ）（ｔ）はＯＦＤＭ変調後の送信アンテナμの送信信号を表し、ｎ（ｔ）は加法的白色ガウス雑音を表し、また、Ｎｓｙｍ＝ＮＦＦＴ＋ＮＧＩはＯＦＤＭシンボルあたりのサンプル数を表す。【００３８】  図２に示すように、受信側においてガードインターバルが除去され、信号サンプルの受信ブロックにＤＦＴを行うことにより情報が復元され、ＯＦＤＭ復調の出力Ｙｌ，ｉが得られる。受信アンテナνにおいて受信したＯＦＤＭ復調後の信号は、以下のように表現される。
2006515481.tif 000006
送信された情報信号と送信アンテナμのチャネル伝達関数（ＣＴＦ）とをそれぞれ表す。Ｎｌ，ｉという語は、平均ゼロで分散Ｎ０の加法的白色ガウス雑音ＡＷＧＮを表す。送信された信号は、それぞれＮｃ個のサブキャリアを有するＬＯＦＤＭシンボルで構成されるものと仮定する。【００３９】チャネル特性【００４０】  以下に、チャネル特性のモデリング方法を説明するが、ここで、本発明はこのモデルなどのいずれの方法にも限定されるものではない。反対に、本発明は実際に存在する移動無線チャネルに適用可能である。本発明の適用に際して、移動無線チャネルの帯域が制限されていればよい。それらが時間においても制限されていることが好ましく、移動無線チャネルの最大遅延を考慮すると実際上そうなる。さらに、異なる送信アンテナおよび受信アンテナの相関がとられていないことが好ましい。【００４１】  本発明では、時変周波数選択性フェージングチャネルが考察される。典型的に、非ゼロ
2006515481.tif 000007
ナμおよび受信アンテナνの間のチャネルインパルス応答ＣＩＲは、以下のように定義される。
2006515481.tif 000008
表す。【００４２】  本発明によると、Ｑ０個のチャネルタップとすべてのアンテナとは相互に相関がとられ
2006515481.tif 000009
あり、その帯域は最大ドップラー周波数νｍａｘによって制限されている。【００４３】  さらに、１つのＯＦＤＭシンボル中チャネルインパルス応答ＣＩＲはほぼ一定であると一般的に仮定されており、よって１つのＯＦＤＭシンボルにおけるＣＩＲの時間依存性を
2006515481.tif 000010
立つ。厳密には、これは非時変チャネルにのみあてはまるものであるが、この仮定は実際には多くの場合に正当化されると考えられ、良好なシステムデザインにおいてはＯＦＤＭシンボル継続時間を必ず十分に短いものとしなければならない。【００４４】  さらに、式（２）のチャネル伝達関数ＣＴＦはＣＩＲｈ（μ，ν）（ｔ，τ）のフーリエ変換を示している。時間ｔ＝ｌＴｓｙｍおよび周波数ｆ＝ｉ／Ｔにおける結果をサンプリングすると、ＯＦＤＭシンボルｌのサブキャリアｉにおけるＣＴＦは、以下のようになる。
2006515481.tif 000011
【００４５】  ここで、Ｔｓｙｍ＝（ＮＦＦＴ＋ＮＧＩ）ＴｓｐｌとＴ＝ＮＦＦＴＴｓｐｌとは、それぞれ、ガードインターバルがある場合とない場合とにおけるＯＦＤＭシンボル継続時間を表す。式（４）の行列形式は、以下のように与えられる。
2006515481.tif 000012
いて、以下のように定義される。
2006515481.tif 000013
【００４６】
2006515481.tif 000014
直交性は維持され、式（２）の受信信号が得られる。【００４７】  受信アンテナにおけるフェージングが相互に相関どりされていないと仮定すると、本発明に係るチャネル推定は各アンテナに対して独立して行われる。よって、以下の記載において、受信アンテナνへのマーカは省略される。【００４８】  さらに、本発明において、チャネルは、タップ遅延τｑがサンプリング時点Ｔｓｐｌの倍数である場合においてサンプルによって間隔が定められていると定義され、例えば、以下が成り立つ。
2006515481.tif 000015
【００４９】
2006515481.tif 000016
ここで、Ｆは次元ＮＧＩ×ＮｃのＤＦＴ−行列を表し、以下のように定義される。
2006515481.tif 000017
【００５０】  実際のチャネルに対しては、タップ遅延τｑサンプリング継続時間Ｔｓｐｌの倍数とはならないが、多くの適用において、サンプリングによって間隔が定められたチャネルモデルは、十分に実際のチャネルに近似したものとなる。【００５１】ＯＦＤＭに基づくＭＩＭＯシステムのためのパイロットシンボル援用チャネル推定の原理【００５２】  図３は、本発明に係る２次元チャネル推定に適した分散パイロットグリッドを示す。【００５３】  図３に示すように、パイロット援用チャネル推定ＰＡＣＥはデータ系列に、公知のシンボル、いわゆるパイロットシンボルを周期的に挿入することに基づく。パイロットシンボルの間隔がサンプリング定理を満たすのに十分に狭い場合、データ系列全体に渡るチャネル推定および補間が可能となる。【００５４】  本発明でいうところの２次元パイロット援用チャネル推定においては、ＯＦＤＭの場合、フェージング変動は時間および周波数の２次元であるということを考慮に入れなくてはならない。従って、２次元サンプリング定理を満たすために、パイロットシンボルは時間−周波数グリッドに渡って分散され、よって２次元パイロットグリッドが得られる。分散パイロットグリッドを選択する別の理由は、スペクトル効果を最大化するためである。【００５５】２次元におけるパイロットグリッドの説明【００５６】  本発明に係る２Ｄ面における通常のグリッドを形式面から説明すると、１つのフレーム
2006515481.tif 000018
とは、それぞれ、周波数方向と時間方向とにおけるパイロット数である）と定義される。ここで、以下の表記法を使用する。２Ｄ構造Ｘを表す行列を想定して、周波数および時間方向に対応する次元を表すサブセットを、それぞれ、Ｘ′およびＸ″と表記する。【００５７】  サブキャリアｉおよびＯＦＤＭシンボルｌのパイロットシンボルをｐ＝［ｉ，ｌ］Ｔと表記することによって、本発明の枠組内において使用されるいずれの通常２Ｄグリッドも、以下のように表される。
2006515481.tif 000019
ロットのインデックスを表し、ｇ０はパイロットグリッドオフセットを定義している。【００５８】  多重入力多重出力ＭＩＭＯ通信システム（以下、ＭＩＭＯシステムともいう）の場合、送信アンテナにおけるすべての送信信号は独自のパイロットグリッドを使用してもよい。このことによって、異なる送信アンテナからの重畳送信信号を受信側で分割することが可能となる。【００５９】  本発明でいうところのパイロットシンボル援用チャネル推定を説明すると、それぞれ、
2006515481.tif 000020
【００６０】
2006515481.tif 000021
【００６１】  図３に示す特定の例において、パイロットグリッド構造は以下のように表されるＧにより定義される。
2006515481.tif 000022
ここでＤｆとＤｔとは、それぞれ、周波数と時間とにおけるいわゆるパイロット間隔を表す。図３に示すパイロットグリッドの場合に特定の値はＤｆ＝５およびＤｔ＝５である。【００６２】  なお、送信前に、すべての送信アンテナに対して同一である外部パイロット系列
2006515481.tif 000023
列が得られる。
2006515481.tif 000024
【００６３】
2006515481.tif 000025
いて低いピークが得られる、かつ／または同期化において良好な相関特性が得られる。【００６４】  図２に示すように、受信側において、巡回プレフィックスが除去され、ＦＦＴが行われ、ＯＦＤＭ復調後の受信信号が得られる。同期が完全であると仮定すると、式（２）の受信信号Ｙｌ，ｉが得られる。チャネル推定のために、パイロット位置の受信信号をデータストリ
2006515481.tif 000026
とで、受信パイロットは以下のように得られる。
2006515481.tif 000027
【００６５】多重入力多重出力通信システムのための２×１Ｄ−パイロット援用チャネル推定の基本的考察【００６６】  図４は、本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。【００６７】  図４に示すように、チャネル推定装置１０は、第１推定装置段階１２と第２推定装置段階１４とを含む。さらに、このチャネル推定装置は、送信アンテナサブセットメモリ１６を含む。【００６８】  図５は、図４に示すチャネル推定装置の操作のフローチャートを示す。【００６９】  図５に示すように、本発明に係るチャネル推定装置の動作は、第１ステップＳ１０において、複数の送信アンテナを分割送信アンテナサブセットに分割する方法によっている。次に、ステップＳ１２において、送信アンテナサブセットに関連して作用するパイロット系列を第１段階チャネル推定を行うことにより分割し、送信の第１次元におけるチャネル応答の仮推定を得る。次に、ステップＳ１４において、送信アンテナサブセットのアンテナに関連して作用するパイロット系列を、各送信アンテナサブセットの各アンテナに対して第２段階チャネル推定を行うことにより分割し、チャネル応答の推定を得る。【００７０】  図６は、本発明に係るチャネル推定の原理を説明する概略図を示す。【００７１】  図６および上記の概説に示すように、ステップＳ１２において、チャネル推定は１次元において行われ、ＯＦＤＭシンボルのすべてのサブキャリアに対する仮推定が得られる。これらの仮推定は、フレーム全体のチャネルを推定するために、ステップＳ１４において新たなパイロットとして使用される。【００７２】  図６に示すように、第２段階はパイロットを有するＯＦＤＭシンボル間に補間するだけのものではなく、仮推定の精度を向上させるものである。【００７３】  さらに、本発明によれば、周波数方向または時間方向のチャネル推定のいずれを最初に行ってもよい。最初に周波数方向におけるチャネル推定を行う場合については、以下に２×１Ｄ−ＰＡＣＥタイプＩとして参照される。さらに、最初に時間方向におけるチャネル推定を行う場合については、以下に２×１Ｄ−ＰＡＣＥタイプＩＩとして参照される。【００７４】  問題を形式面から説明するために、ＯＦＤＭシンボルｌＤｔの受信パイロットは、（ｉＤｆ）番目のサブキャリアにおいて、以下のように考えられる。
2006515481.tif 000028
関数（ＣＴＦ）とを表す。【００７５】  さらに、ＣＴＦは、時間および周波数などの変数ｌおよびｉにおいて変動すると仮定さ
2006515481.tif 000029
ＯＦＤＭシンボル数を表し、ＮｃはＯＦＤＭシンボルあたりのサブキャリア数を表し、ＤｆおよびＤｔは周波数および時間におけるパイロット間隔を表し、ＮＴは送信アンテナ数を表す。【００７６】
2006515481.tif 000030
推定することである。【００７７】  本目的を達成するために、一般的に、ＭＩＭＯチャネル推定においては、チャネル推定および補間に加えてＮＴ個の送信アンテナから作用する信号を分割することが必要となる。ＮＴ個の送信アンテナに対応する信号をＮＴ１個のサブセットに分割することで、それぞれＮＴ２個の信号を有するＮＴ１個のグループが得られ、よってＮＴ＝ＮＴ１ＮＴ２となる。言い換えると、ＮＴ個の送信アンテナを有するＭＩＭＯシステムに対して、本発明によると、ＮＴ２個の送信アンテナに対応する信号が１つのセットにグループ化され、ＮＴ１個のサブセットを形成する。【００７８】  言い換えると、本発明の基礎をなす基本概念は、ステップＳ１２でＯＦＤＭシンボル
2006515481.tif 000031
ャリアに対して仮推定を得るという方法により、分割タスクを２つの段階に分けることにある。第２ステップＳ１４では、フレーム全体のチャネルを推定するために、これら仮推定を新たなパイロットとして使用する。【００７９】  図６に示すように、第２段階推定はパイロットを有するＯＦＤＭシンボル間に補間するだけのものではなく、仮推定の精度を向上させるものである。従って、以下に詳細に説明する本発明の異なる実施形態においては、パフォーマンスはほとんど低下しない一方で、複雑性が大幅に減少されている。【００８０】  図６によると、最初に時間方向におけるチャネル推定を行うという順序とは反対に、最初に周波数方向におけるチャネル推定を行うことは簡単である。最初に周波数方向におけるチャネル推定を行う場合については、以下に２×１Ｄ−ＰＡＣＥタイプＩとして参照される。反対に、最初に時間方向におけるチャネル推定を行う場合については、２×１Ｄ−ＰＡＣＥタイプＩＩとして参照される。【００８１】
2006515481.tif 000032
ングが必要となる。【００８２】
2006515481.tif 000033
う別の方法もある。この場合、時間方向におけるチャネル推定に対しては、先行するパイロットのみが使用される。【００８３】  予測フィルタリングがバッファリングを必要としないことは明らかであるが、平滑化に関するパフォーマンスは低下する。一方で、周波数方向におけるチャネル推定では、１つのＯＦＤＭシンボルのすべてのパイロットが同時に受信されるので、バッファリングを必要としない。しかしながら、典型的にはチャネル推定の精度が帯域端近くまで低下する。【００８４】２×１Ｄ−ＰＡＣＥタイプＩ【００８５】  図７は、チャネル推定が最初に周波数方向において行われる、本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。【００８６】  本発明によると、２×１Ｄ−ＰＡＣＥをＯＦＤＭに基づくＭＩＭＯチャネル推定に拡大することが提案されている。ＭＩＭＯチャネル推定においては、ＮＴ個の送信アンテナから作用する信号を分割することが必要となる。ＮＴ個の送信アンテナを有するＭＩＭＯシステムをセットとしてＡと表す。さらに、本発明によると、ＮＴ個の送信アンテナのセッ
2006515481.tif 000034
トにはＮＴ２個の送信アンテナが含まれ、よってＮＴ＝ＮＴ１ＮＴ２となる。【００８７】  言い換えると、本発明によるとＮＴ個の送信アンテナを有するＭＩＭＯシステムにおいては、ＮＴ２個の送信アンテナに対応する信号を１つのセットとしてグループ化し、ＡのＮＴ１個のサブセットを形成することが提案される。保護範囲を限定することなく、すべて
2006515481.tif 000035
アンテナが１つのセットとしてのみ扱われると仮定してもよい。【００８８】  本発明の基礎をなす概念は、第１推定段階において最初にＮＴ１＜ＮＴ信号のサブセットを分割し、チャネル推定を併せて行う方法により、分割タスクを２つの段階に分けること
2006515481.tif 000036
【００８９】  第２推定段階において、残りのＮＴ２個の重畳信号が第１推定段階のＮＴ１個の信号のそれぞれに対して分割され、第２次元におけるチャネル推定を併せて行うことにより、周波
2006515481.tif 000037
【００９０】  図７は、提案されたスキームのタイプＩに関するする基本的な考えを説明するものである。【００９１】  なお、図７に示すバッファは時間方向において平滑化フィルタリングを適用するために使用される。しかしながら、上記に概説したように、本発明は、平滑化および予測フィルタリングの両方に適用可能であるため、図７に示すバッファは任意に用いられ、使用される特定のチャネル推定アルゴリズムに依存する。【００９２】  図８および９は、本発明に係るＯＦＤＭに基づくＭＩＭＯチャネル推定のための２段階手法を支えるために使用されうるパイロット系列デザインを示す。ここで、図９に示すパイロットグリッドは、ＥＴＳＩ  ＥＮ３００  ７４４、Ｖ１．４．１（２００１−０１）に準拠したＤＶＢ−Ｔパイロットグリッドに対応する。保護範囲を限定するものとして考えられるべきではないが、さらなる規格として、ＩＥＥＥ  ８０２．１１ａ、またはＥＴＳＩ  ＴＳ  １０１  ４７５  ＨＩＰＥＲＬＡＮ／２が挙げられる。【００９３】  より一般的な意味において、多重送信アンテナを用いたＯＦＤＭに基づくシステムに対しパイロット系列デザインを拡大するためには、本発明の好適な実施形態によると、送信
2006515481.tif 000038
【００９４】  以下のような積の形で表されるパイロット系列を選択することが提案される。
2006515481.tif 000039
ンボルである。【００９５】
2006515481.tif 000040
存する。【００９６】
2006515481.tif 000041
あることを意味する。それに対応して、１つのサブキャリアのパイロットを表す系列
2006515481.tif 000042
【００９７】
2006515481.tif 000043
直交デザインから選択されることが好ましい。【００９８】  （１２）の提示された２Ｄパイロット系列を、（１１）の受信パイロットに代入すると、以下が得られる。
2006515481.tif 000044
2006515481.tif 000045
【００９９】
2006515481.tif 000046
【０１００】
2006515481.tif 000047
ループを分割し、さらに、周波数方向においてチャネルを推定および補間することにある。【０１０１】
2006515481.tif 000048
ことを意味する。【０１０２】  さらに、第１推定段階の出力は後に、７に示す第２推定段階のための入力として用いられる。【０１０３】  図７に示すように、第２推定段階において、チャネルは時間方向において推定され、サ
2006515481.tif 000049
【０１０４】  本発明では、この枠組みを使用して、ＯＦＤＭに基づくＭＩＭＯシステムのための、利用可能な１次元スキームを適用することでチャネル推定を行うことが提案される。【０１０５】２×１Ｄ−ＰＡＣＥタイプＩＩ【０１０６】  図１０は、チャネル推定が最初に時間方向において行われる、本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。【０１０７】  図１０に示すように、２×１Ｄ−ＰＡＣＥタイプＩＩの２×１Ｄ−ＰＡＣＥタイプＩに対する主たる相違点は、第１推定段階におけるＮＴ１個のサブセットの分割が、時間方向におけるチャネル推定とともに行われることにある。このことにより、２×１Ｄ−ＰＡＣＥタイプＩＩのパイロットシンボルが得られる。
2006515481.tif 000050
【０１０８】  上記から、パイロット系列が式（１２）における２×１Ｄ−ＰＡＣＥタイプＩと非常に
2006515481.tif 000051
2006515481.tif 000052
された２Ｄパイロット系列を、式（１１）の受信パイロットに代入すると、受信パイロット系列は以下のように表現することができる。
2006515481.tif 000053
【０１０９】
2006515481.tif 000054
信パイロット系列は、以下のように定義される。
2006515481.tif 000055
【０１１０】
2006515481.tif 000056
こと、すなわち、ＮＴ１個のグループを分割し、その後、時間方向においてチャネルを、ｌ変数などを介して推定および補間することにある。【０１１１】
2006515481.tif 000057
にΔｌＤｔＯＦＤＭシンボルはバッファリングを必要とする。第２推定段階に対しては、チャネルを周波数方向に推定してサブセットあたりの残りのＮＴ２個の信号を分割することが提案される。【０１１２】２×１Ｄ−ＰＡＣＥタイプＩとタイプＩＩとの比較【０１１３】  タイプＩおよびタイプＩＩの両スキームのパフォーマンスは、ほとんどの実施において同一であることが分かる。しかし、チャネル推定に必要とされる乗算の数において、計算複雑性が異なりうる。実際の計算複雑性はその多くを、システムパラメータ、パイロットグリッド構造、および使用されるチャネル推定アルゴリズムに依存している。さらに、ＤＳＰまたはハードウェアアーキテクチャによっては一方のスキームが好ましい場合もある。【０１１４】  しかしながら、より重要なことは、あるパイロットグリッドを選択することで、特定のスキームの実施が不可能となる点にある。図８および図９に示すパイロットグリッドにはタイプＩおよびタイプＩＩの両スキームが適用可能である。しかしながら、これは一般的ではない。この問題を掘り下げるために、図８に示す以下のグリッドを考察する。
2006515481.tif 000058
このようなグリッド構造に対しては、周波数方向のパイロットが一列に配置されるために２×１Ｄ−ＰＡＣＥタイプＩを採用することがより適切である。一方で、時間方向に連続するパイロットは１つのサブキャリア分離れて変位するため、２×１Ｄ−ＰＡＣＥタイプＩＩの採用は困難となる。【０１１５】  従って、図９に示す以下のグリッドを選択すると、２×１Ｄ−ＰＡＣＥタイプＩＩの採用がより適切となる。
2006515481.tif 000059
図９のこのグリッドはＤＶＢ−Ｔ規格に対して選択されてきたものであり、よって、ＤＶＢ−Ｔに対しては、２×１Ｄ−ＰＡＣＥが使用される場合、時間方向における本発明の好適な実施形態に係るチャネル推定を最初に行うべきである。【０１１６】２×１Ｄ−ＰＡＣＥタイプＩの適用【０１１７】  本発明に係るパイロットシンボル援用チャネル推定を説明すると、周波数方向における
2006515481.tif 000060
いてそれぞれサンプリングされた、パイロットのみを含む受信信号系列のサブセット、す
2006515481.tif 000061
【０１１８】  以下においては、２つのチャネル推定技術が２×１Ｄ−ＰＡＣＥタイプＩに対して説明される。１つ目は、ウィーナフィルタ補間による周波数領域におけるチャネル推定であり、２つ目の手法は、受信信号の時間領域への伝達である。【０１１９】周波数領域チャネル推定【０１２０】第１段階ウィーナフィルタリング【０１２１】  ２×１Ｄ−ＰＡＣＥタイプＩの場合、第１のステップはチャネルを周波数方向に推定することである。ベクトル表記法においては、ＯＦＤＭシンボルｌの受信パイロット系列は、以下のようになる。
2006515481.tif 000062
信パイロット系列を表す。【０１２２】  さらに、アンテナμｌから送信されるＯＦＤＭシンボルｌの送信パイロット系列、サブ
2006515481.tif 000063
2006515481.tif 000064
【０１２３】  さらに、本発明によると、周波数領域におけるチャネル推定はＦＩＲ補間フィルタを用いて行われることが好ましく、これは周波数方向における第１段階に対して、以下のように表現することができる。
2006515481.tif 000065
【０１２４】
2006515481.tif 000066
ボルの位置に依存している。これは、各送信アンテナだけではなく各サブキャリアに対しても、別のフィルタが必要となることを意味する。【０１２５】
2006515481.tif 000067
とである。ウィーナフィルタにより、パイロット系列および所望の応答の間における平均二乗誤差ＭＳＥが最小となる。これは最小ＭＳＥまたは同等にＭＭＳＥ推定装置としても公知である。【０１２６】  さらに、ウィーナフィルタを生成するためには、共分散行列で表されるチャネル統計の
2006515481.tif 000068
より定義される。共分散行列のｍ番目の列およびｎ番目の行における値は、以下のように与えられる。
2006515481.tif 000069
【０１２７】
2006515481.tif 000070
のように表現することができる。
2006515481.tif 000071
【０１２８】  式（２０）および式（２１）による数値はウィーナ補間フィルタを評価するために必要である。ＭＭＳＥに関して最適の解決法を、以下によるウィーナ・ホッフを使用して定義してもよい。
2006515481.tif 000072
【０１２９】第２段階ウィーナフィルタリング【０１３０】  上記に鑑みて、ＮＴ１のサブセットおよび各サブキャリアに対する式（１９）の計算は、第２段階においてさらに処理される。ここで、時間方向におけるフィルタリングおよび補間によって、以下の形式で得られる周波数応答推定が得られる。
2006515481.tif 000073
階のＦＩＲ補間フィルタを表す。【０１３１】
2006515481.tif 000074
低い推定になる代わりに、チャネル推定による時間遅延を引き起こすことなく線形予測受信側を特定する。【０１３２】  なお、一般的に、推定されるべきシンボルがそのシンボルの推定に使用されるパイロッ
2006515481.tif 000075
シンボルが関与する。【０１３３】  さらに、行列表記法によると、式（２３）は以下のようになる。
2006515481.tif 000076
【０１３４】
2006515481.tif 000077
補間フィルタを使用することである。時間方向における第１段階出力の共分散行列は、
2006515481.tif 000078
における値は、以下のように与えられる。
2006515481.tif 000079
【０１３５】
2006515481.tif 000080
の値は、以下のように表現することができる。
2006515481.tif 000081
【０１３６】  式（２０）および式（２１）による数値はウィーナ補間フィルタを評価するために必要である。ＭＭＳＥに関して最適の解決法は、以下によるウィーナ・ホッフ式を使用して得られる。
2006515481.tif 000082
【０１３７】時間領域チャネル推定【０１３８】  図１１は、本発明に係る時間領域のチャネル推定を実現するように構成された推定装置段階の概略図を示す。【０１３９】  図１１に示すように、本発明に係る周波数応答Ｈｌ，ｉを定義するためのもうひとつの手
2006515481.tif 000083
ある。μ番目の送信アンテナから作用するＯＦＤＭシンボルｌのＣＩＲベクトルに関し
2006515481.tif 000084
のようになる。
2006515481.tif 000085
【０１４０】
2006515481.tif 000086
で説明したサンプリングによって間隔が定められたチャネルモデルにのみ当てはまる。Ｈｌ，ｉのＤＦＴが時間制限されない、サンプリングによって間隔が定められていないチャネルモデルの場合は、エイリアシングを避けるためにオーバーサンプリングが必要となる。この場合、Ｑが有効なタップ数を表す。【０１４１】
2006515481.tif 000087
【０１４２】  これは、ＯＦＤＭに基づくＭＩＭＯシステムの２段階チャネル推定が、時間および周波数領域において同様に採用可能であるということを意味する。【０１４３】  ＯＦＤＭ復調後の受信パイロット系列は、以下のように与えられる。
2006515481.tif 000088
ここで、以下が成り立つ。
2006515481.tif 000089
【０１４４】
2006515481.tif 000090
域に変換する。【０１４５】
2006515481.tif 000091
わち、以下のようになる。
2006515481.tif 000092
【０１４６】
2006515481.tif 000093
わち、最小二乗ＬＳ推定装置および最小平均二乗誤差ＭＭＳＥ推定装置について説明する。【０１４７】最小二乗ＬＳ推定装置【０１４８】
2006515481.tif 000094
定義されてもよい。
2006515481.tif 000095
【０１４９】  推定装置は送信信号に依存するために、パイロット系列は適切に選択しなくてはならな
2006515481.tif 000096
ての場合に当てはまらない。ＬＳ推定装置が存在するために必要な条件として、以下のものがある。
2006515481.tif 000097
【０１５０】  本発明の好適な実施形態によると、２倍のオーバーサンプリングにより、パイロットによるシステムオーバーヘッドを最小化すること、およびパフォーマンスを最適化すること
2006515481.tif 000098
れる。【０１５１】  なお、複数の送信アンテナに対するＬＳ推定装置は時間領域のみに存在する。【０１５２】最小平均二乗誤差ＭＭＳＥ推定装置【０１５３】  ＭＭＳＥ推定装置は、時間領域チャネル推定のためのＦＩＲフィルタにより与えられる。
2006515481.tif 000099
【０１５４】  ウィーナフィルタはウィーナ・ホッフ式により定義される。
2006515481.tif 000100
【０１５５】
2006515481.tif 000101
【０１５６】
2006515481.tif 000102
【０１５７】  ここで、異なる送信アンテナ間のフェージングは相関がとられていないと仮定される。
2006515481.tif 000103
【０１５８】
2006515481.tif 000104
ーナフィルタは、以下のように簡略化可能となる。
2006515481.tif 000105
【０１５９】  なお、本発明によると、ＬＳ推定装置により行われるＮＴ１個の信号の分割は、フィルタリングタスクより分けることが可能である。【０１６０】  上記に概説したように、一般的にＭＭＳＥ推定装置はパイロットシンボルの選択に依存
2006515481.tif 000106
【０１６１】  従って、式（３３）および式（４０）において必要とされる行列反転が簡単になるので、式（４０）のＭＭＳＥ推定装置と同様に、式（３３）の上記ＬＳ推定装置も著しく簡略化可能となる。【０１６２】  さらに、式（３３）および式（４０）から、推定装置は選択したパイロット系列から独立し、よってフィルタ生成が大幅に簡単となることが分かる。【０１６３】  図１１は、直交パイロット系列を使用した時間領域におけるチャネル推定および補間のブロック図を示す。最初に、受信パイロット系列をＮＴ１個のブランチに分割し、各ブラン
2006515481.tif 000107
する。任意のフィルタリングおよび／またはウィンドウィングにより、チャネルは時間領域において推定される。第１段階推定のゼロパディングにより、Ｎｃ個のサンプルにまで拡大される。ＯＦＤＭシンボル全体（パイロットおよびデータ）のＣＴＦ推定は、ＣＩＲ推定のＮｃポイントＦＦＴにより得られる。
2006515481.tif 000108
2006515481.tif 000109
ＮＴ１Ｎｃ×ＮＴ１Ｎｃブロック対角行列を表す。【０１６４】
2006515481.tif 000110
ることが可能である。【０１６５】
2006515481.tif 000111
ＮＴ１回のＦＦＴにより周波数領域に変換される。さらに、ドップラー領域における第２段階のチャネル推定も可能であり、すなわち、時間領域と同等のアルゴリズムを使用して
2006515481.tif 000112
【０１６６】  図１２は、本発明に係る時間領域のチャネル推定を実現するように構成された推定装置段階の概略図を示す。【０１６７】  図１２より、第２段階のチャネル推定も、上記した第１推定段階に対応するＤＦＴ補間を用いて行ってもよいことが分かる。しかしながら、ゼロパディングおよびＮｃポイントＦＦＴ前などに、第２推定段階を時間領域においても行うことが計算上より効率的である。【０１６８】  図１２に示すように、送信の第１および第２次元など周波数および時間において相関関数が相互に独立しているため、第２段階のフィルタリングそれ自体は影響を受けないままとなる。よって、ウィーナフィルタリングが第２段階に選択された場合、上記の式（２
2006515481.tif 000113
プルインデックスを表す。第２段階フィルタリングの後、ＮＴ個の出力全体、例えばサブセットあたりＮＴ２個の出力はそれぞれ周波数領域に戻るように変換される。【０１６９】周波数再利用係数１を用いたセルラーシステムへの適用【０１７０】  図１３は、本発明に係る２段階チャネル推定手法の、周波数再利用係数が１であるセルラー通信システムへの適用を示す。【０１７１】  図１３に示すように、ＮＴ個のアンテナ素子を含むアンテナアレイを有する代わりに、提案したスキームは分散したアンテナにも適用可能である。例えば、２×１Ｄ−ＰＡＣＥが周波数再利用係数が１であるセルラーシステムに適用される。モバイルユーザがセル境界にいる場合、ユーザは１つの基地局から所望の信号を受信し、さらに他の基地局から１つ以上の干渉信号を受信する。【０１７２】  さらに、各基地局はＮＴ２個のアンテナ素子を有すると仮定してもよい。データを含むシンボルは、チャネルコードまたは拡散を用いて干渉から保護可能であるが、パイロットシンボルはこのようには保護されない。しかしながら、高精度のチャネル推定はシステムが効率よく作動するために最も重要である。解決法の１つはパイロットを増加することであるが、これは他の基地局のサービスを受けるユーザに対する干渉を増やすこととなり、よってシステムの容量を制限してしまう。【０１７３】  本発明によると、２×１Ｄ−ＰＡＣＥはこのような場合に以下のように適用可能である。基地局においてＮＴ１個のサブセットを形成し、各サブセットはＮＴ２個のアンテナ素子を含むアンテナアレイを有し、結果としてＮＴ＝ＮＴ１ＮＴ２個の素子を有するアレイを形成する。これには、セル間の同期が必要となる。【０１７４】略語【０１７５】ＡＷＧＮ      加法的白色ガウス雑音
ＣＩＲ        チャネルインパルス応答
ＣＴＦ        チャネル伝達関数
ＤＦＴ        離散フーリエ変換
ＦＦＴ        高速フーリエ変換
ＧＩ          ガードインターバル
ＩＤＦＴ      逆離散フーリエ変換
ＩＦＦＴ      逆高速フーリエ変換
ＬＳ          最小二乗
ＭＩＭＯ      多重入力多重出力、一般的には、いくつかの送信および受信アンテナを有するシステム
ＭＭＳＥ      最小平均二乗誤差
ＭＳＥ        平均二乗誤差
ＯＦＤＭ      直交周波数分割多重方式
ＰＡＣＥ      パイロットシンボル援用チャネル推定【０１７６】共通して使用されるシステムパラメータリスト【０１７７】ＮＦＦＴ      ＦＦＴ長
Ｎｃ          サブキャリア数
ＮＧＩ        ガードインターバルサンプル数
Ｌ            フレームあたりＯＦＤＭシンボル数
Ｔ            ＯＦＤＭシンボル継続時間
Ｔｓｐｌ      Ｔｓｐｌ＝Ｔ／ＮＦＦＴにより与えられるサンプルインターバル
Ｔｓｙｍ      ガードインターバルＴｓｐｌ＝Ｔ＋ＮＧＩＴｓｐｌを含むＯＦＤＭシンボル合計継続時間
Ｑ０          非ゼロチャネルタップ数
Ｑ            チャネルタップ合計数
ＮＲ          受信アンテナ数
ＮＴ          送信アンテナ数
ＮＴ１        サブセット数
ＮＴ２        １つのサブセット内の送信アンテナ数
2006515481.tif 000114
Ｄｆ          周波数におけるパイロット間隔
Ｄｔ          時間におけるパイロット間隔【０１７８】共通して使用される変数リスト【０１７９】
2006515481.tif 000115
Ｎｌ，ｉ      サブキャリアｉにおけるＯＦＤＭシンボルｌの平均ゼロで分散Ｎ０のＡＷＧＮサンプル
ｎ（ｔ）      ＯＦＤＭ復調前の時間ｔにおけるＡＷＧＮ過程の実現
【図面の簡単な説明】
【０１８０】
【図１】  本発明の基礎をなすシステムモデルを説明するためのＯＦＤＭに基づく多重入力多重出力ＭＩＭＯ通信システムの概略図を示す。
【図２】  それぞれ、ＯＦＤＭ変調および復調を説明する概略図を示す。
【図３】  本発明に係る２次元チャネル推定に適した分散パイロットグリッドを示す。
【図４】  本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。
【図５】  図５に示すチャネル推定装置の動作のフローチャートを示す。
【図６】  本発明の基礎をなす２×１Ｄチャネル推定の原理を説明する概略図を示す。
【図７】  チャネル推定が最初に周波数方向において行われる、本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。
【図８】  本発明に係る２次元チャネル推定に適したさらなる分散パイロットグリッドを示す。
【図９】  デジタルビデオブロードキャストＤＶＢ−Ｔアプリケーションに対応し、本発明に係る２次元チャネル推定に適したさらなる分散パイロットグリッドを示す。
【図１０】  チャネル推定が最初に時間方向において行われる、本発明に係るマルチキャリア通信システムの多重入力多重出力送信チャネルを推定するためのチャネル推定装置の概略図を示す。
【図１１】  本発明に係る時間領域のチャネル推定を実現するように構成された推定装置段階の概略図を示す。
【図１２】  本発明に係る時間領域のチャネル推定を実現するように構成されたさらなる推定装置段階の概略図を示す。
【図１３】  本発明に係る２段階チャネル推定手法の、周波数再利用係数が１であるセルラー通信システムへの適用を示す。
【図１】
2006515481.tif 000163
【図２】
2006515481.tif 000164
【図３】
2006515481.tif 000165
【図４】
2006515481.tif 000166
【図５】
2006515481.tif 000167
【図６】
2006515481.tif 000168
【図７】
2006515481.tif 000169
【図８】
2006515481.tif 000170
【図９】
2006515481.tif 000171
【図１０】
2006515481.tif 000172
【図１１】
2006515481.tif 000173
【図１２】
2006515481.tif 000174
【図１３】
2006515481.tif 000175
【手続補正書】
【提出日】平成17年12月22日(2005.12.22)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
複数（ＮＴ）の送信アンテナから作用し、パイロットシンボルが挿入された２次元データ系列を搬送するマルチキャリア変調送信信号を用いた多重入力多重出力送信システム用２次元チャネル推定方法であって、
2006515481.tif 000176
テップと、
2006515481.tif 000177
に、送信の第１次元におけるチャネル応答の仮推定を得るための第１段階チャネル推定を行うステップと、
2006515481.tif 000178
応答の推定を得るための第２段階チャネル推定を行うステップと
を有する方法であって、
前記第１段階チャネル推定は、パイロット系列を２次元グリッドのパイロットシンボルに配置して用いることにより実行され、前記第１段階チャネル推定に用いられるパイロットシンボルは送信の第１次元にのみ依存し、前記第２段階チャネル推定に用いられるパイロットシンボルは送信の第２次元にのみ依存する
ことを特徴とする方法。
【請求項２】
前記パイロット系列は送信の前記第１次元と送信の前記第２次元に分離可能な積形式で表現されることを特徴とする請求項１に記載の方法。
【請求項３】
前記パイロット系列は直交計画から選択されることを特徴とする請求項２に記載の方法。
【請求項４】
前記直交計画はウォルシュ系列および位相シフト系列を含むグループから選択されることを特徴とする請求項３に記載の方法。
【請求項５】
前記第１次元は周波数次元であることを特徴とする請求項１または２のいずれかに記載の方法。
【請求項６】
前記パイロット系列が以下の積形式で表現され、
2006515481.tif 000179
上記式において、
2006515481.tif 000180
示し、
2006515481.tif 000181
トシンボルを示し、
2006515481.tif 000182
トシンボルを示し、
2006515481.tif 000183
信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しい、すなわち、ＮＴ＝ＮＴ１・ＮＴ２である
ことを特徴とする請求項５に記載の方法。
【請求項７】
前記第１段階チャネル推定は、下記式に従い、パイロットシンボルと受信したパイロットシンボルの値との積形式を用いて行われ、
2006515481.tif 000184
上記式において、
2006515481.tif 000185
ことを特徴とする請求項６に記載の方法。
【請求項８】
前記第１段階チャネル推定は周波数領域において行われることを特徴とする請求項５乃至７のいずれかに記載の方法。
【請求項９】
前記第１段階チャネル推定は時間領域において行われることを特徴とする請求項５乃至７のいずれかに記載の方法。
【請求項１０】
前記第２段階チャネル推定は周波数領域において行われることを特徴とする請求項８または９に記載の方法。
【請求項１１】
前記第２段階チャネル推定は時間領域において行われることを特徴とする請求項８または９に記載の方法。
【請求項１２】
前記第１次元は時間次元であることを特徴とする請求項１または２に記載の方法。
【請求項１３】
前記パイロット系列が以下の積形式で表現され、
2006515481.tif 000186
上記式において、
2006515481.tif 000187
示し、
2006515481.tif 000188
トシンボルを示し、
2006515481.tif 000189
トシンボルを示し、
2006515481.tif 000190
信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しい、すなわち、ＮＴ＝ＮＴ１・ＮＴ２である
ことを特徴とする請求項１２に記載の方法。
【請求項１４】
前記第１段階チャネル推定は、下記式に従い、パイロットシンボルと受信したパイロットシンボルの値との積形式を用いて行われ、
2006515481.tif 000191
上記式において、
2006515481.tif 000192
ことを特徴とする請求項１３に記載の方法。
【請求項１５】
前記第１段階チャネル推定は周波数領域において行われることを特徴とする請求項１２乃至１４のいずれかに記載の方法。
【請求項１６】
前記第１段階チャネル推定は時間領域において行われることを特徴とする請求項１２乃至１４のいずれかに記載の方法。
【請求項１７】
前記第２段階チャネル推定は周波数領域で行われることを特徴とする請求項１５または１６に記載の方法。
【請求項１８】
前記第２段階チャネル推定は時間領域において行われることを特徴とする請求項１５または１６に記載の方法。
【請求項１９】
基地局とそれに対応するアンテナアレイが複数の送信アンテナ
2006515481.tif 000193
数の送信アンテナに関係付けて定義される周波数再利用ファクタを伴う携帯電話通信システムにおいて適用されることを特徴とする請求項１乃至１８のいずれかに記載の方法。
【請求項２０】
複数の送信アンテナ（ＮＴ）から作用し、パイロットシンボルが挿入された２次元データ系列を搬送するマルチキャリア変調送信信号を使用する多重入力多重出力送信システム用チャネル推定装置であって、
2006515481.tif 000194
記憶するメモリと、
2006515481.tif 000195
送信の第１次元におけるチャネル応答の仮推定を得るための第１段階チャネル推定を行う少なくとも１つの第１段階推定手段と、
2006515481.tif 000196
推定を得るための第２段階チャネル推定を行う少なくとも１つの第２段階推定手段と
を有し、
第１段階推定手段は、２次元グリッドのパイロットシンボルに配置されたパイロット系列を使用して前記第１段階チャネル推定を行い、前記第１段階チャネル推定に使用されるパイロットシンボルは送信の第１次元にのみ依存し、前記第２段階チャネル推定に使用されるパイロットシンボルは送信の第２次元にのみ依存する
ことを特徴とするチャネル推定装置。
【請求項２１】
前記第１段階推定手段は、第１次元送信と送信の第２次元に分離可能な積形式で表現されたパイロット系列を用いて前記第１段階チャネル推定を行うことを特徴とする請求項２０に記載のチャネル推定装置。
【請求項２２】
前記第１段階推定手段は、直交計画から選択されたパイロット系列を用いて前記第１段階推定を行うことを特徴とする請求項２０または２１に記載のチャネル推定装置。
【請求項２３】
前記第１段階推定手段は、ウォルシュ系列もしくは位相シフト系列を用いて前記第１段階チャネル推定を行うことを特徴とする請求項２２に記載のチャネル推定装置。
【請求項２４】
前記第１段階推定手段は、前記第１段階チャネル推定を周波数次元において行うことを特徴とする請求項２０または２１に記載のチャネル推定装置。
【請求項２５】
前記第１段階推定手段は、以下の積形式で表現されるパイロット系列を用いて前記第１段階チャネル推定を行い、
2006515481.tif 000197
上記式において、
2006515481.tif 000198
示し、
2006515481.tif 000199
ンボルを示し、
2006515481.tif 000200
ンボルを示し、
2006515481.tif 000201
信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しい、すなわち、ＮＴ＝ＮＴ１・ＮＴ２である
ことを特徴とする請求項２４に記載のチャネル推定装置。
【請求項２６】
前記第１段階推定手段は、パイロットシンボルと受信したパイロットシンボルの値との積形式を用いて下記式に従い前記第１段階チャネル推定を行い、
2006515481.tif 000202
上記式において、
2006515481.tif 000203
ことを特徴とする請求項２５に記載のチャネル推定装置。
【請求項２７】
前記第１段階推定手段は、前記第１段階チャネル推定を周波数領域において行うことを特徴とする請求項２４乃至２６のいずれかに記載のチャネル推定装置。
【請求項２８】
前記第１段階推定手段は、前記第１段階チャネル推定を時間領域において行うことを特徴とする請求項２４乃至２６のいずれかに記載のチャネル推定装置。
【請求項２９】
前記第２段階推定手段は、前記第２段階チャネル推定を周波数領域において行うことを特徴とする請求項２７または２８に記載のチャネル推定装置。
【請求項３０】
前記第２段階推定手段は、前記第２段階チャネル推定を時間領域において行うことを特徴とする請求項２７または２８に記載のチャネル推定装置。
【請求項３１】
第１段階推定手段は、前記第１段階チャネル推定を時間方向において行うことを特徴とする請求項２０乃至２１のいずれかに記載のチャネル推定装置。
【請求項３２】
前記第１段階推定手段は、以下の積形式で表現されるパイロット系列を用いて前記第１段階チャネル推定を行い、
2006515481.tif 000204
上記式において、
2006515481.tif 000205
示し、
2006515481.tif 000206
ンボルを示し、
2006515481.tif 000207
ンボルを示し、
2006515481.tif 000208
信アンテナサブセットに含まれるアンテナの数の積が複数の送信アンテナの数ＮＴと等しい、すなわち、ＮＴ＝ＮＴ１・ＮＴ２である
ことを特徴とする請求項３１に記載のチャネル推定装置。
【請求項３３】
前記第１段階推定手段は、下記式に従い、パイロットシンボルと受信したパイロットシンボルの値との積形式を用いて前記第１段階チャネル推定を行い、
2006515481.tif 000209
上記式において、
2006515481.tif 000210
2006515481.tif 000211
ことを特徴とする請求項３２に記載のチャネル推定装置。
【請求項３４】
前記第１段階推定手段は、前記第１段階チャネル推定を周波数領域において行うことを特徴とする請求項３１乃至３３のいずれかに記載のチャネル推定装置。
【請求項３５】
前記第１段階推定手段は、前記第１段階チャネル推定を時間領域において行うことを特徴とする請求項３１乃至３３のいずれかに記載のチャネル推定装置。
【請求項３６】
前記第２段階推定手段は、前記第２段階チャネル推定を周波数領域において行うことを特徴とする請求項３４または３５に記載のチャネル推定装置。
【請求項３７】
前記第２段階推定手段は、前記第２段階チャネル推定を時間領域において行うことを特徴とする請求項３４または３５に記載のチャネル推定装置。
【請求項３８】
チャネル推定装置の内部メモリに直接的にロード可能なコンピュータプログラム製品であって、多重入力多重出力送信チャネルを２次元で推定するチャネル推定装置のプロセッサにおいて実行された場合に、請求項１乃至１９のいずれかに記載の各ステップを実行させるソフトウェアコード部分を備える、多重入力多重出力送信チャネルを２次元で推定するチャネル推定装置の内部メモリに直接的にロード可能なコンピュータプログラム製品。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００３
【補正方法】変更
【補正の内容】
【０００３】  多重入力多重出力ＭＩＭＯ通信システムとして公知の多重送受信アンテナを採用した通信システムを、直交周波数分割多重ＯＦＤＭと併用して、移動無線システムの通信容量および品質を向上することが可能である。ナギブ（Ａ．Ｎａｇｕｉｂ）、シシャドリ（Ｎ．Ｓｅｓｈａｄｒｉ）およびカルダーバンク（Ａ．Ｃａｌｄｅｒｂａｎｋ）の『無線高速データ通信のための時空間符号化および信号処理』（ＩＥＥＥ信号処理誌、７６〜９２頁、２０００年５月）において記載されている時空間符号などの多重送信アンテナを用いた直交周波数分割多重ＯＦＤＭ通信システム、または特殊な多重方式においては、異なる信号が異なる送信アンテナから同時に送信される。従って、受信信号はこれらの信号が重畳されたものとなり、これがチャネル推定における課題を意味する。チャネルパラメータは、時空間符号を使用した場合にはダイバーシチの結合のために必要となり、また、空間多重方式を使用した場合には重畳信号の分割のために必要となる。
ウォン・ギ・ジェオン（Ｗｏｎ  Ｇｉ．  Ｇｅｏｎ  ｅｔ  ａｌ）には、送信ダイバーシチを有するＯＦＤＭシステムに用いる２次元ＭＭＳＥチャネル推定が記載されている。具体的には、送信ダイバーシチを有するコヒーレント直交周波数分割多重システムに用いる２次元パイロットシンボル型チャネル推定技術が記載されている。パイロットシンボル型チャネル推定装置と、この推定装置に対応して、送信ダイバーシチを有するＯＦＤＭシステムにおいて多重送信アンテナと受信アンテナ間のチャネルの変化を追跡するために用いる最小二乗誤差の各条件下で、パイロットシンボルを直交状態とすることが導き出されている。
ＥＰ１，２７６，２８８Ａ１では、マルチキャリア送信を用いたチャネル推定用の参照シンボルが記載されており、ここでは、ＯＦＤＭマルチキャリア変調とＯＦＤＭマルチキャリア復調の少なくともいずれかが実行される。さらには、パイロットパターンフレームを構成する少なくとも１つのパイロットサブキャリアパターンに基づいて、チャネル推定、チャネル等化、および／または同期化が実行される。このパイロットパターンフレームは、パイロットサブキャリアパターンとして送信され、送信される拡張ＯＦＤＭシンボルストリームの始点を誤りなく識別するのに用いられる参照シンボルを有する周期先頭フィールドを１つ、および／または、チャネル推定を実行するプリアンブルフィールドを少なくとも１つ有する。


【国際調査報告】
2006515481.tif 000212
2006515481.tif 000213
2006515481.tif 000214