(57)【要約】【課題】伝送特性を維持しつつ、チャネル推定のためのトレーニング信号の数を低減させることを可能とする。【解決手段】Ｎ行Ｎ列のユニタリ行列で表されるランダムウエイトの一部のベクトルを選択してチャネル情報の推定に用いられるトレーニングシンボルの重みづけとなるトレーニングウエイトとして生成し、トレーニングウエイトとトレーニングウエイトのベクトル数に応じた大きさのトレーニングシンボルに基づいてウエイトトレーニングシンボルを含んだ送信信号を生成して端末装置に送信し、端末装置側でウエイトトレーニングシンボルに基づいて推定されたウエイトチャネル情報を端末装置から受信し、ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトを選択し、トレーニングウエイトとウエイトチャネル情報と付加トレーニングウエイトに基づいて新たなトレーニングウエイトを生成する。【選択図】図１
2014045286.tif 000002
【特許請求の範囲】
【請求項１】
Ｎ本（Ｎ≧２）アンテナを有し、１つ以上の端末装置と通信を行う送信装置であって、
Ｎ行Ｎ列のユニタリ行列で表されるランダムウエイトを記憶し、当該ランダムウエイトの一部のベクトルを選択してチャネル情報の推定に用いられるトレーニングシンボルの重みづけとなるトレーニングウエイトとして生成するトレーニングウエイト生成回路と、
前記トレーニングウエイトと前記トレーニングウエイトのベクトル数に応じた大きさのトレーニングシンボルに基づいてウエイトトレーニングシンボルを生成し、当該ウエイトトレーニングシンボルを含んだ送信信号を生成する送信信号生成回路と、
前記ウエイトトレーニングシンボルを含んだ送信信号を前記端末装置に送信し、前記端末装置側で前記ウエイトトレーニングシンボルに基づいて推定されたウエイトチャネル情報を前記端末装置から受信する無線信号送受信回路と、を備え、
前記トレーニングウエイト生成回路は、
前記端末装置からのウエイトチャネル情報を取得すると、ウエイトチャネル情報推定に用いた前記トレーニングウエイトとウエイトチャネル情報から信号空間トレーニングウエイトを算出し、前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトを選択し、前記信号空間トレーニングウエイトと前記付加トレーニングウエイトに基づいて新たなトレーニングウエイトを生成する
ことを特徴とする送信装置。
【請求項２】
前記トレーニングウエイト生成回路は、
前記新たなトレーニングウエイトを生成する際に、前記ランダムウエイトの選択されてないベクトルの数が、前記付加トレーニングウエイトが必要とするベクトル数未満または０の場合、前記ウエイトチャネル情報と前記トレーニングウエイトとに基づいて新たにＮ行Ｎ列のユニタリ行列で表されるランダムウエイトを生成して記憶する
ことを特徴とする請求項１に記載の送信装置。
【請求項３】
前記トレーニングウエイト生成回路は、
計時手段が計時する経過時間に応じて、前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトとして選択するベクトル数を増加させる
ことを特徴とする請求項１または２に記載の送信装置。
【請求項４】
前記端末装置から受信するウエイトチャネル情報は、前記端末装置側で推定したウエイトチャネル行列のうちフィードバックを要求された所定の数のベクトルで表される情報であり、
前記トレーニングウエイト生成回路は、
前記所定の数の変更を前記端末装置に要求する
ことを特徴とする請求項１から３のいずれか１つに記載の送信装置。
【請求項５】
前記端末装置ごとの前記所定のベクトル数の総和と前記付加トレーニングウエイトのベクトル数との和を一定とし、
前記トレーニングウエイト生成回路は、
前記所定のベクトル数の変更を前記端末装置に要求した場合、前記所定のベクトル数の変更に応じて前記付加トレーニングウエイトのベクトル数を増減させて前記端末装置ごとの前記所定のベクトル数の総和と前記付加トレーニングウエイトのベクトル数の和を一定にする
ことを特徴とする請求項４に記載の送信装置。
【請求項６】
前記端末装置が複数の場合に、前記端末装置の組み合わせが複数定められ、
前記トレーニングウエイト生成回路は、
複数の前記端末装置の組み合わせごとにランダムウエイトを記憶し、前記トレーニングウエイトの生成を行う
ことを特徴とする請求項１から５のいずれか１つに記載の送信装置。
【請求項７】
前記端末装置が複数の場合に、前記端末装置の組み合わせが定められ、
前記トレーニングウエイト生成回路は、
更新前のトレーニングウエイトについて記憶し、組み合わせに含まれる前記端末装置のうち一部の前記端末装置からウエイトチャネル情報を取得できなかった場合、当該端末装置からのウエイトチャネル情報を正常に取得した際のトレーニングウエイトを記憶し、ウエイトチャネル情報を正常に取得した端末装置からのウエイトチャネル情報と用いたトレーニングウエイトに加え、ウエイトチャネル情報を取得できなかった端末装置の過去のウエイトチャネル情報と当該ウエイトチャネル情報を推定した際のトレーニングウエイトを用いてトレーニングウエイトを生成して更新する
ことを特徴とする請求項１から６のいずれか１つに記載の送信装置。
【請求項８】
Ｎ本（Ｎ≧２）アンテナを有し、１つ以上の端末装置と通信を行う送信装置におけるチャネル推定方法であって、
Ｎ行Ｎ列のユニタリ行列で表されるランダムウエイトの一部のベクトルを選択してチャネル情報の推定に用いられるトレーニングシンボルの重みづけとなるトレーニングウエイトとして生成し、
前記トレーニングウエイトと前記トレーニングウエイトのベクトル数に応じた大きさのトレーニングシンボルに基づいてウエイトトレーニングシンボルを含んだ送信信号を生成して前記端末装置に送信し、
前記端末装置側で前記ウエイトトレーニングシンボルに基づいて推定されたウエイトチャネル情報を前記端末装置から受信し、
前記端末装置からのウエイトチャネル情報を取得すると、ウエイトチャネル情報推定に用いた前記トレーニングウエイトとウエイトチャネル情報から信号空間トレーニングウエイトを算出し、
前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトを選択し、前記信号空間トレーニングウエイトと前記付加トレーニングウエイトに基づいて新たなトレーニングウエイトを生成する
ことを特徴とするチャネル推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば、複数のアンテナを用いて通信を行う送信装置、及びチャネル推定方法に関する。
【背景技術】
【０００２】
近年、２．４ＧＨｚ帯または５ＧＨｚ帯を用いた高速無線アクセスシステムとして、ＩＥＥＥ８０２．１１ｇ規格、ＩＥＥＥ８０２．１１ａ規格などの普及が目覚しい。これらのシステムでは、マルチパスフェージング環境での特性を安定化させるための技術である直交周波数分割多重（ＯＦＤＭ： Orthogonal Frequency Division Multiplexing）変調方式を用い、最大で５４Ｍｂｐｓの物理層伝送速度を実現している。
これらのシステムでの伝送速度とは物理レイヤ上での伝送速度であり、実際にはＭＡＣ（Medium Access Control）レイヤでの伝送効率が５０〜７０％程度であるため、実際のスループットの上限値は３０Ｍｂｐｓ程度であり、情報を必要とする通信相手が増えればこの特性は更に低下する。一方、有線ＬＡＮの世界ではＥｔｈｅｒｎｅｔ（登録商標）の１００Ｂａｓｅ−Ｔインタフェースをはじめ、各家庭にも光ファイバを用いたＦＴＴＨ（Fiber to the home）の普及から、１００Ｍｂｐｓの高速回線の提供が普及しており、無線ＬＡＮの世界においても更なる伝送速度の高速化が求められている。
高速化の技術として、ＩＥＥＥ８０２．１１ｎにおいて空間多重送信技術としてＭＩＭＯ（Multiple input multiple output）技術が導入され、オプションで４素子のアンテナ数までサポートされた。さらに、ＩＥＥＥ８０２．１１ａｃでは、マルチユーザＭＩＭＯ（ＭＵ(Multi User)−ＭＩＭＯ）通信方法が検討され、サポートするアンテナ素子数は８素子まで増やされている（例えば、非特許文献１参照）。
【０００３】
例えば、この技術を用いたシステムの一例として、図６に示すような無線通信システムがある。図６に示す無線通信システムは、基地局装置１ａと端末装置２ａ−１〜２ａ−ｉを備える。基地局装置１ａは、Ｎ本の送受信アンテナ１ａ−４−１〜１ａ−４−Ｎを有しており、各端末装置２ａ−ｉ（ｉ＝１〜Ｋ）は、それぞれＭｉ本の送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉを有している。ここで、Ｍｉは、ｉ番目の端末装置の送受信アンテナ数であり、１≦ｉ≦Ｋである。Ｋは、端末装置の数である。
【０００４】
基地局装置１ａから端末装置２ａ−ｉへの送信処理について説明する。まず、最初に基地局装置１aは、通信を行う対象となる各端末装置２a−ｉを送信方法選択回路１a−７により決定する。通信を行う対象は、データ選択出力回路から得られるデータを送る相手となる端末装置の情報、送信を行うデータの種類や優先度、端末装置に対する信号対雑音電力比の情報、などから決定することができる。基地局装置１ａは、データの送信を行う前に、各端末装置２ａ−ｉの送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉとの間の伝搬環境の情報を示すチャネル情報を収集する。送信信号生成回路１ａ−２は、各端末装置２ａ−ｉの送受信アンテナとの間のチャネル情報を推定するためトレーニングシンボルを生成し、生成したトレーニングシンボルにガードインターバルや制御信号の付加を行い信号を生成する。無線信号送受信回路１ａ−３は、生成された信号を搬送波周波数へアップコンバートして送受信アンテナ１ａ−４−１〜１ａ−４−Ｎを介して送信を行う。
【０００５】
端末装置２ａ−iの無線信号送受信回路２ａ−ｉ−２は、基地局装置１ａが送信した信号を送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉを介して受信する。受信信号復調回路２ａ−ｉ−３は、無線信号送受信回路２ａ−ｉ−２が受信した信号に対して同期、復調を行う。チャネル推定回路２ａ−ｉ−４は、復調された信号に対してチャネル推定を行い、チャネル推定したチャネル情報を送信信号生成回路２ａ−ｉ−５に出力する。送信信号生成回路２ａ−ｉ−５は、チャネル情報を含めた送信信号を生成し、予め定められたタイミングで、無線信号送受信回路２ａ−ｉ−２に出力する。無線信号送受信回路２ａ−ｉ−２は、送信信号生成回路２ａ−ｉ−５が出力したチャネル情報を含む送信信号を送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉの少なくとも１つを介して送信する。  なお、チャネル推定回路２ａ−ｉ−４は、送信信号生成回路２ａ−ｉ−５に対してチャネル情報を出力するのではなく、チャネル情報を圧縮するなどの処理をして生成したチャネル関連情報を出力するようにしてもよい。
【０００６】
基地局装置１ａの無線信号送受信回路１ａ−３は、送受信アンテナ１ａ−４−１〜１ａ−４−Ｎのうち少なくとも１つを介して端末装置２ａ−ｉから送信された信号を受信する。無線信号送受信回路１ａ−３は、受信した信号をダウンコンバートしデジタルデータに変換し、受信信号復調回路１ａ−５に出力する。受信信号復調回路１ａ−５は、受信信号を復調するとともに、受信信号に含まれるチャネル情報またはチャネル関連情報を読み出し、チャネル情報記憶回路１ａ−６に書き込んで記憶させる。
【０００７】
次に、基地局装置１ａのデータ選択出力回路１ａ−１は、通信相手となる端末装置２ａ−ｉのチャネル情報がチャネル情報記憶回路１ａ−６に記憶されると、データを送信するため、選択した通信相手となる端末装置２ａ−ｉの情報を送信方法選択回路１ａ−７に出力し、送信するデータを送信信号生成回路１ａ−２に出力する。送信方法選択回路１ａ−７は、チャネル情報記憶回路１ａ−６に記憶されたチャネル情報から送信ウエイトを算出し、空間ストリーム数、変調方式、符号化率からなる送信モードと算出した送信ウエイトを送信信号生成回路１ａ−２に出力する。送信信号生成回路１ａ−２は、通信を行う端末装置２ａ−ｉに送信する送信信号に対して、出力された送信モードに従って変調、及び符号化を行い、送信ウエイトを乗算し、信号検出や通信情報伝達に用いるパイロット信号を挿入して無線信号送受信回路１ａ−３に出力する。無線信号送受信回路１ａ−３は、送信信号生成回路１ａ−２が出力した信号を搬送波周波数にアップコンバートし、送受信アンテナ１ａ−４−１〜１ａ−４−Ｎを介して送信する。
【０００８】
データ選択出力回路１ａ−１において通信相手を選択する手法としては、送信を行うデータがメモリに記憶されており、その中から送信を行う準備ができている端末装置２ａ−ｉを選択する手法や、記憶されている複数のデータのうち最も古いデータに対応する端末装置２ａ−ｉから選択する手法がある。また、端末装置２ａ−ｉを利用するユーザに割り当てられているＱｏＳ（Quality of service）に基づいて選択したり、予めグループＩＤで決められたユーザに対応する端末装置２ａ−ｉの組み合わせを選択したり、チャネル情報の相関が低い端末装置２−ｉの組み合わせを選択したりする手法もある。
【０００９】
上記の通信手順は、基地局装置１ａと複数の端末装置２ａ−ｉとの間の全体的な送信処理の流れであり、以下、空間多重方式を用いて通信を行う一例として、ブロック直交化（ＢＤ：Block Diagonalization）指向性制御法を用いた例について説明する（例えば、非特許文献２、３参照）。
基地局装置１ａは、端末装置２ａ−１〜２ａ−Ｋの送受信アンテナに対するチャネル情報を取得するため、送信信号生成回路１ａ−２が、チャネル推定用のトレーニングシンボルを生成する。図７は、トレーニングシンボルの構成例を示した図である。図７の四角いブロックＬＴ（１，１）〜ＬＴ（Ｎ，Ｎ）は、チャネル推定のためのＯＦＤＭシンボルを表しており、ガードインターバル（ＧＩ：Guard Interval）が付与されたものとなっている。ＬＴ（ｊ，ｋ）は、トレーニングシンボルブロックのうち、ｋ番目のタイミングにおいてｊ番目の送受信アンテナ１ａ−４−ｊから送信されるＯＦＤＭシンボルに対応する。ＬＴ（ｊ，ｋ）の生成は、周波数チャネルにおける既知信号ｓ１〜ｓＦを逆フーリエ変換してＧＩを付与することによって行われる。ＬＴ（１，１）〜ＬＴ（Ｎ，Ｎ）のｎ番目の周波数チャネルの送信信号ｘｊ，ｋ，ｎからなる送信信号行列Ｘｎは、以下の式（１）として表される。
【００１０】
【数１】
2014045286.tif 000003
【００１１】
式（１）において、ｘｊ，ｋ，ｎは、ＬＴ（ｊ，ｋ）のｎ番目の周波数チャネルに対応する送信信号を表し、行列Ａは、トレーニングシンボル用の符号行列であり送信側と受信側の両側で既知の行列が用いられる。行列Ａとして、単位行列Ｉを用いれば、各送受信アンテナからそれぞれ異なるタイミングで、ｓｎを送信することができる。また、行列ＡとしてＡＨＡ＝Ｉとなるような直交行列を用いることもできる。
端末装置２ａ−ｉが、基地局装置１ａからの信号を受信すると、受信信号復調回路２ａ−ｉ−３は、受信信号に対し同期を行う。チャネル推定回路２ａ−ｉ−４は、ＬＴ（１，１）〜ＬＴ（Ｎ，Ｎ）に対応する受信信号から、チャネル情報である基地局装置１ａと端末装置２ａ−ｉの送受信アンテナ間の伝搬係数からなるチャネル行列を推定して送信信号生成回路２ａ−ｉ−５に出力する。ｊ番目のトレーニングシンボルに対応する受信信号において、送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉで受信されたｎ番目の周波数チャネルに対応する受信信号ｙ１，ｊ，ｉ，ｎ〜ｙＭｉ，ｊ，ｉ，ｎ からなる受信信号ベクトルｙｊ，ｉ，ｎは、以下の式（２）として表される。
【００１２】
【数２】
2014045286.tif 000004
【００１３】
式（２）において、行列Ｈｉ，ｎは、端末装置２ａ−ｉのｎ番目の周波数チャネルに対するチャネル情報を表すチャネル行列（Ｍｉ×Ｎ行列）である。チャネル行列Ｈｉ，ｎのｐ列ｑ行目の要素は、基地局装置１ａのｐ番目の送信アンテナ（送受信アンテナ１ａ−４−ｐ）から、端末装置２ａ−ｉのｑ番目の受信アンテナ（送受信アンテナ２ａ−ｉ―１−ｑ）間の伝搬係数を表す。ベクトルｎｊ，ｉ，ｎは、ｊ番目のトレーニングシンボルの受信タイミングで、送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉが受信した端末装置２ａ−ｉのｎ番目の周波数チャネルにおける熱雑音ベクトルを表す。したがって、端末装置２ａ−ｉの第ｎの周波数チャネルにおける１〜Ｎ番目のトレーニングシンボルに対応する受信信号行列Ｙｉ，ｎは、以下の式（３）として表すことができる。
【００１４】
【数３】
2014045286.tif 000005
【００１５】
式（３）において、行列Ｎｉ，ｎは、端末装置２ａ−ｉのｎ番目の周波数チャネルおいて、送受信アンテナ２ａ−ｉ−１−１〜２ａ−ｉ−１−Ｍｉが１〜Ｎ番目の受信タイミングで受信した信号に対する熱雑音行列を表す。端末装置２ａ−ｉにおいて、行列Ａとｓｎは既知であるため、Ｘｎは既知である。したがって、チャネル行列Ｈ’ｉ，ｎは、受信信号行列Ｙｉ，ｎから、以下の式（４）によって求められる。
【００１６】
【数４】
2014045286.tif 000006
【００１７】
端末装置２ａ−ｉは、式（４）によって算出されたチャネル行列Ｈ’ｉ，ｎをそのまま基地局装置１ａへのフィードバック情報とすることもできるし、行列Ｈ’ｉ，ｎの中からＬｉ個の行ベクトルを選択してフィードバックすることもできる。また、行列Ｈ’ｉ，ｎに対して特異値分解を行って式（５）のように表し、信号空間に対応する右特異行列Ｖ（ｓ）ｉ，ｎ（Ｎ×Ｍｉ行列）のうちＬｉ個の列ベクトルを選択してフィードバックすることもできる。
【００１８】
【数５】
2014045286.tif 000007
【００１９】
式（５）において、行列Ｕｉ，ｎは、左特異行列である。行列Σｉ，ｎは、特異値を対角要素とし、非対角項が０のＭｉ×Ｎの対角行列である。行列Ｖ（ｎ）ｉ，ｎ（Ｎ×（Ｎ−Ｍｉ）行列）は、右特異行列のうち特異値に対応しないベクトルの集まりを表す。上記の説明では、基地局装置１ａの送受信アンテナ数Ｎが、端末装置２ａ−ｉの送受信アンテナ数Ｍｉより多いことを前提としていた。これに対して、端末装置２ａ−ｉの送受信アンテナ数Ｍｉの方がＮよりも多い場合には、右特異行列Ｖ（ｓ）ｉ，ｎ（Ｎ×Ｎ行列）の列ベクトルの全てが特異値と対応するため右特異行列Ｖ（ｓ）ｉ，ｎのうちＬｉ個の列ベクトルを選択してフィードバックすることができる。また、Ｈ’ｉ，ｎのエルミート行列の行ベクトルにグラムシュミットの直交化法を用いて得られる基底ベクトルのうち、Ｌｉ個の列ベクトルを選択してフィードバックしてもよい。また、直交化法による基底ベクトルは、式（６）のようにＱＲ分解によっても得ることができる。
【００２０】
【数６】
2014045286.tif 000008
【００２１】
式（６）により得られた直交行列（Ｑ（ｓ）ｉ，ｎ  Ｑ（ｎ）ｉ，ｎ）のうちＭｉ×Ｍｉの三角行列Ｒｉ，ｎに対応するＱ（ｓ）ｉ，ｎのうちＬｉ個の列ベクトルをフィードバックするようにしてもよい。
【００２２】
上記において、Ｌｉは、基地局装置１ａが、端末装置２ａ−ｉにフィードバックを要求するチャネルの次元であり、実際に基地局装置１ａが端末装置２ａ−ｉに送信を行う空間ストリーム数Ｌ’ｉはＬｉ以下に設定される。行列Ｑ（ｓ）ｉ，ｎや行列Ｖ（ｓ）ｉ，ｎは、ユニタリ行列の一部となっており、以下、チャネル情報から得られるこれらの行列を信号空間行列Ｇｉ，ｎ、または信号空間情報と定義する。
【００２３】
端末装置２ａ−１〜２ａ−Ｋの送信信号生成回路２ａ−１−５〜２ａ−Ｋ−５は、信号空間行列Ｇｉ，ｎの情報を含んだ送信信号を生成する。無線信号送受信回路２ａ−１−２〜２ａ−Ｋ−２が、生成された送信信号を無線信号送受信回路２ａ−１−２〜２ａ−Ｋ−２を介して基地局装置１ａに送信する。このとき空間相関や周波数相関やユニタリ行列の特徴等を用いて再生可能な形にＬｉ個のベクトルを圧縮することもできる。
基地局装置１ａの無線信号送受信回路１ａ−３は、送受信アンテナ１ａ−４−１〜１ａ−４−Ｎを介して端末装置２ａ−１〜２ａ−Ｋからの信号を受信し、受信信号復調回路１ａ−５は、受信した信号からそれぞれチャネル情報、もしくはチャネル情報から算出された信号空間情報を読み出してチャネル情報記憶回路１ａ−６に書き込んで記憶させる。
【００２４】
基地局装置１ａにおいて、端末装置２ａ−１〜２ａ−Ｋに対して送信を行う際には、送信方法選択回路１ａ−７が、フィードバックされたＬｉ個のベクトルの情報をチャネル情報記憶回路１ａ−６から読み出す。端末装置２ａ−ｉに対して、第ｎ番目の周波数チャネルに対応して得られたＬｉ個のベクトルからなるＮ×Ｌｉの行列をＧｉ，ｎと定義する。行列Ｇｉ，ｎは、推定したチャネル行列のうちのＬｉ個の行ベクトルの複素共役転置行列でもよい。また、行列Ｇｉ，ｎは、Ｌｉ個のデータストリームに対する受信ウエイト行列αｉ，ｎを端末装置２ａ−ｉにおいて予め定義して式（７）としてもよい。
【００２５】
【数７】
2014045286.tif 000009
【００２６】
また、行列Ｇｉ，ｎは、特異値分解により得られた右特異行列Ｖ（ｓ）ｉ，ｎやＱＲ分解により得られたＱ（ｓ）により、式（８）、式（９）とすることもできる。
【００２７】
【数８】
2014045286.tif 000010
【００２８】
【数９】
2014045286.tif 000011
【００２９】
式（８）、（９）において、［Ａ］Ｌは、Ｎ×Ｍ行列ＡからＬ個の列ベクトルを選択してＮ×Ｌの行列を得る関数である。端末装置２ａ−ｉにおける受信ウエイト行列αｉ，ｎとしては、Ｈ’ｉ，ｎをＱＲ分解して得られるＭｉ×Ｍｉのユニタリ行列のうち、対応する三角行列の対角成分が大きいものからＬｉ個のベクトルをＭｉ×Ｍｉのユニタリ行列の中から選択して得られるＭｉ×Ｌｉの行列にエルミート共役をとることで得られるＬｉ×Ｍｉの行列を用いたり、Ｈ’ｉ，ｎに対し特異値分解で得られる左特異行列のうち、大きい特異値に対応するＬｉ個のベクトルを選択し得られるＭｉ×Ｌｉ行列に対し、エルミート共役をとることで得られるＬｉ×Ｍｉの行列などを用いることができる。
【００３０】
次に、マルチユーザＭＩＭＯの通信方法の例として、ＢＤ法による送信ウエイトを算出する手法について説明する。ここで、Ｋユーザ（端末装置２ａ−１〜２ａ−Ｋ）に対し、通信を行うことを考える。ｉ番目の端末装置２ａ−ｉに対する送信ウエイトの算出方法を示す。まず、端末装置２ａ−ｉ以外の端末装置に対応する集合信号空間行列Ｇ＋ｉ，ｎを、式（１０）として定義する。
【００３１】
【数１０】
2014045286.tif 000012
【００３２】
このＧ＋ｉ，ｊに対し、特異値分解を行うと、式（１１）として表される。
【００３３】
【数１１】
2014045286.tif 000013
【００３４】
式（１１）において、行列Ｖ（ｓ）＋ｉ，ｎは、固有値Σ＋ｉ，ｎに対応する信号空間ベクトルであり、行列Ｖ（ｎ）＋ｉ，ｎは、固有値がない、または固有値０に対応するヌル空間ベクトルである。ここで、行列Ｖ（ｎ）＋ｉ，ｎで表せるヌル空間に対して送信を行うと、端末装置２ａ−ｉ以外の端末装置の受信ウエイトに対して干渉を生じない。したがって、複数の端末装置に空間多重方式を用いて通信するには、ｎ番目の周波数チャネルに用いる送信ウエイトとして、ここで得られたＶ（ｎ）＋ｉ，ｎに線形演算を行って得られる値を用いればよいことになる。
例えば、端末装置２ａ−ｉに対応する信号空間行列Ｇｉ，ｎのエルミート行列Ｇｉ，ｎＨに、行列Ｖ（ｎ）＋ｉ，ｎを乗算してＧｉ，ｎＨＶ（ｎ）＋ｉ，ｎを得る。得られたＧｉ，ｎＨＶ（ｎ）＋ｉ，ｎの行ベクトルに対して直交化法を用いて得られる基底ベクトルのエルミート行列をＶ（ｎ）＋ｉ，ｎに乗算することにより得られる値を送信ウエイトとすることができる。また、得られたＧｉ，ｎＨＶ（ｎ）＋ｉ，ｎに特異値分解を行って得られる右特異ベクトルをＶ（ｎ）＋ｉ，ｎに乗算することにより得られる値を送信ウエイトとすることもできる。Ｇｉ，ｎＨＶ（ｎ）＋ｉ，ｎから得られた行列をＤｉ，ｎとすると、送信ウエイトＷｉ，ｎは、［Ｖ（ｎ）＋ｉ，ｎＤｉ，ｎ］L’iとして表される。
【００３５】
以上、ＢＤ法に基づくＭＵ−ＭＩＭＯ送信法について説明したが、Ｚｅｒｏ Ｆｏｒｃｉｎｇ法や、ＭＭＳＥ（Minimum Mean Square Error）法、Ｓｕｃｃｅｓｓｉｖｅ ｏｐｔｉｍｉｚａｔｉｏｎ法、Ｔｏｍｌｉｎｓｏｎ Ｈａｒａｓｈｉｍａ Ｐｒｅｃｏｄｉｎｇ、Ｄｉｒｔｙ ｐａｐｅｒ ｃｏｄｉｎｇなどを用いて送信ウエイトＷｉ，ｎを演算することもできる。
また、ユーザ数が１（Ｋ＝１）の場合には、信号空間行列Ｇｉ，ｎを送信ウエイトとして用いることができ、また信号空間行列Ｇｉ，ｎの右特異行列を送信ウエイトとして用いることもできる。
【００３６】
このように各端末装置２ａ−ｉに対してそれぞれ送信ウエイトを算出することができ、得られた端末装置２ａ−１〜２ａ−Ｋに対するｎ番目の周波数チャネルに対する送信ウエイトＷｎは、式（１２）として表される。
【００３７】
【数１２】
2014045286.tif 000014
【先行技術文献】
【非特許文献】
【００３８】
【非特許文献１】IEEE, “Proposed specification framework for TGac,” doc.: IEEE 802.11-09/0992r21, Jan. 2011.
【非特許文献２】T. Yoo, and A. Goldsmith, “Optimality of zero-forcing beamforming with multiuser diversity,” in Proc. ICC 2005, vol. 1, pp. 542-546, May 2005.
【非特許文献３】Q. H. Spencer, A. L. Swindlehurst, and M. Haardt, “Zero-forcing methods for downlink spatial multiplexingin multiuser MIMO channels,”IEEE Trans. Sig. Processing, vol. 52, issue 2, Feb. 2004.
【発明の概要】
【発明が解決しようとする課題】
【００３９】
上述したＢＤ法、Ｚｅｒｏ Ｆｏｒｃｉｎｇ法、ＭＭＳＥ法、Ｓｕｃｃｅｓｓｉｖｅ ｏｐｔｉｍｉｚａｔｉｏｎ法、Ｔｏｍｌｉｎｓｏｎ Ｈａｒａｓｈｉｍａ Ｐｒｅｃｏｄｉｎｇ、Ｄｉｒｔｙ ｐａｐｅｒ ｃｏｄｉｎｇなど、いずれのＭＵ−ＭＩＭＯの送信ウエイトを算出する手法においても、送信ウエイトを用いて送信を行い受信側で復号演算を行うことで、ユーザ間干渉をキャンセルさせるか、または、ある期待値まで低減させることができる。
【００４０】
しかしながら、いずれの手法でも、図７に示すように送信アンテナ数分のトレーニング信号が必要であり、アンテナ数が増えるほど、チャネル推定のためのパイロット信号の数が増大し、また、端末装置におけるチャネル推定の演算負荷及び端末装置からのフィードバック量も増大し、チャネル情報推定のために多くの時間を必要とし、スループットを低下させる問題があった。
【００４１】
本発明は、このような事情に鑑みてなされたもので、その目的は、伝送特性を維持しつつ、チャネル推定のためのトレーニング信号の数を低減させることを可能とする送信装置、及びチャネル推定方法を提供することにある。
【課題を解決するための手段】
【００４２】
上記問題を解決するために、本発明の一態様は、Ｎ本（Ｎ≧２）アンテナを有し、１つ以上の端末装置と通信を行う送信装置であって、Ｎ行Ｎ列のユニタリ行列で表されるランダムウエイトを記憶し、当該ランダムウエイトの一部のベクトルを選択してチャネル情報の推定に用いられるトレーニングシンボルの重みづけとなるトレーニングウエイトとして生成するトレーニングウエイト生成回路と、前記トレーニングウエイトと前記トレーニングウエイトのベクトル数に応じた大きさのトレーニングシンボルに基づいてウエイトトレーニングシンボルを生成し、当該ウエイトトレーニングシンボルを含んだ送信信号を生成する送信信号生成回路と、前記ウエイトトレーニングシンボルを含んだ送信信号を前記端末装置に送信し、前記端末装置側で前記ウエイトトレーニングシンボルに基づいて推定されたウエイトチャネル情報を前記端末装置から受信する無線信号送受信回路と、を備え、前記トレーニングウエイト生成回路は、前記端末装置からのウエイトチャネル情報を取得すると、ウエイトチャネル情報推定に用いた前記トレーニングウエイトとウエイトチャネル情報から信号空間トレーニングウエイトを算出し、前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトを選択し、前記信号空間トレーニングウエイトと前記付加トレーニングウエイトに基づいて新たなトレーニングウエイトを生成することを特徴とする送信装置である。
【００４３】
また、本発明の一態様は、上記に記載の発明において、前記トレーニングウエイト生成回路は、前記新たなトレーニングウエイトを生成する際に、前記ランダムウエイトの選択されてないベクトルの数が、前記付加トレーニングウエイトが必要とするベクトル数未満または０の場合、前記ウエイトチャネル情報と前記トレーニングウエイトとに基づいて新たにＮ行Ｎ列のユニタリ行列で表されるランダムウエイトを生成して記憶することを特徴とする。
【００４４】
また、本発明の一態様は、上記に記載の発明において、前記トレーニングウエイト生成回路は、計時手段が計時する経過時間に応じて、前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトとして選択するベクトル数を増加させることを特徴とする。
【００４５】
また、本発明の一態様は、上記に記載の発明において、前記端末装置から受信するウエイトチャネル情報は、前記端末装置側で推定したウエイトチャネル行列のうちフィードバックを要求された所定の数のベクトルで表される情報であり、前記トレーニングウエイト生成回路は、前記所定の数の変更を前記端末装置に要求することを特徴とする。
【００４６】
また、本発明の一態様は、上記に記載の発明において、前記端末装置ごとの前記所定のベクトル数の総和と前記付加トレーニングウエイトのベクトル数との和を一定とし、前記トレーニングウエイト生成回路は、前記所定のベクトル数の変更を前記端末装置に要求した場合、前記所定のベクトル数の変更に応じて前記付加トレーニングウエイトのベクトル数を増減させて前記端末装置ごとの前記所定のベクトル数の総和と前記付加トレーニングウエイトのベクトル数の和を一定にすることを特徴とする。
【００４７】
また、本発明の一態様は、上記に記載の発明において、前記端末装置が複数の場合に、前記端末装置の組み合わせが複数定められ、前記トレーニングウエイト生成回路は、複数の前記端末装置の組み合わせごとにランダムウエイトを記憶し、前記トレーニングウエイトの生成を行うことを特徴とする。
【００４８】
また、本発明の一態様は、上記に記載の発明において、前記端末装置が複数の場合に、前記端末装置の組み合わせが定められ、前記トレーニングウエイト生成回路は、更新前のトレーニングウエイトについて記憶し、組み合わせに含まれる前記端末装置のうち一部の前記端末装置からウエイトチャネル情報を取得できなかった場合、当該端末装置からのウエイトチャネル情報を正常に取得した際のトレーニングウエイトを記憶し、ウエイトチャネル情報を正常に取得した端末装置からのウエイトチャネル情報と用いたトレーニングウエイトに加え、ウエイトチャネル情報を取得できなかった端末装置の過去のウエイトチャネル情報と当該ウエイトチャネル情報を推定した際のトレーニングウエイトを用いてトレーニングウエイトを生成して更新することを特徴とする。
【００４９】
また、本発明の一態様は、Ｎ本（Ｎ≧２）アンテナを有し、１つ以上の端末装置と通信を行う送信装置におけるチャネル推定方法であって、Ｎ行Ｎ列のユニタリ行列で表されるランダムウエイトの一部のベクトルを選択してチャネル情報の推定に用いられるトレーニングシンボルの重みづけとなるトレーニングウエイトとして生成し、前記トレーニングウエイトと前記トレーニングウエイトのベクトル数に応じた大きさのトレーニングシンボルに基づいてウエイトトレーニングシンボルを含んだ送信信号を生成して前記端末装置に送信し、前記端末装置側で前記ウエイトトレーニングシンボルに基づいて推定されたウエイトチャネル情報を前記端末装置から受信し、前記端末装置からのウエイトチャネル情報を取得すると、ウエイトチャネル情報推定に用いた前記トレーニングウエイトとウエイトチャネル情報から信号空間トレーニングウエイトを算出し、前記ランダムウエイトの選択されていないベクトルから付加トレーニングウエイトを選択し、前記信号空間トレーニングウエイトと前記付加トレーニングウエイトに基づいて新たなトレーニングウエイトを生成することを特徴とするチャネル推定方法である。
【発明の効果】
【００５０】
この発明によれば、伝送特性を維持しつつ、チャネル推定のためのトレーニング信号の数を低減させることが可能となる。
【図面の簡単な説明】
【００５１】
【図１】本発明の一実施形態による無線通信システムを示す概略ブロック図である。
【図２】同実施形態によるトレーニングウエイトの生成及び更新の処理を示すフローチャートである。
【図３】同実施形態によるウエイトトレーニングシンボルの構成を示す図である。
【図４】同実施形態による伝送可能ビット数を示した図である。
【図５】同実施形態による伝送可能ビット数の累積確率分布を示した図である。
【図６】ＭＵ−ＭＩＭＯ方式による無線通信システムを示す概略ブロック図である。
【図７】トレーニングシンボルの構成を示す図である。
【発明を実施するための形態】
【００５２】
以下、本発明の一実施形態による無線通信システムおよび基地局装置１について図面を参照して説明する。
図１に本発明の一実施形態による無線通信システムを表す概略ブロック図を示す。図１において、無線通信システムは、指向性制御によるビームを前提としたシステムであり、基地局装置１と、複数の端末装置２−１〜２−ｉ〜２−Ｋを備える。ここで、Ｋは、端末装置の数であり、１≦ｉ≦Ｋである。基地局装置１は、Ｎ本の送受信アンテナ１−４−１〜１−４−Ｎ、データ選択出力回路１−１、送信信号生成回路１−２、無線信号送受信回路１−３、受信信号復調回路１−５、ウエイトチャネル情報記憶回路１−６、送信方法選択回路１−７、及びトレーニングウエイト生成回路１−８を備える。各端末装置２−ｉは、それぞれＭｉ本の送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉ、無線信号送受信回路２−ｉ−２、受信信号復調回路２−ｉ−３、チャネル推定回路２−ｉ−４、及び送信信号生成回路２−ｉ−５を備える。ここで、Ｍｉはｉ番目の端末装置の送受信アンテナ数である。なお、端末装置２−ｉと記載した場合は、いずれか１つの端末装置を代表して表す場合と、各々の端末装置を表す場合があり、端末装置が備える各機能ブロックについても同様とする。
【００５３】
基地局装置１から端末装置２−ｉへの送信処理について説明する。基地局装置１はデータの送信を行う前に、各端末装置２−ｉのアンテナとの間の伝搬環境の情報を示すチャネル情報を収集する。本実施形態の基地局装置１では、チャネル情報を収集するためのトレーニングシンボルに重みづけをするトレーニングウエイトを用いて、新たなトレーニングシンボルを生成する。トレーニングウエイト生成回路１−８は、ランダムウエイトから選択してトレーニングウエイトを生成する。また、トレーニングウエイト生成回路１−８は、各端末装置２−ｉのウエイトチャネル情報と、ランダムウエイトと、既に生成したトレーニングウエイトに基づいて新たなトレーニングウエイトを生成してトレーニングウエイトの更新を行う。ランダムウエイト及び生成されたトレーニングウエイトは、トレーニングウエイト生成回路１−８が内部に備える記憶領域に記憶される。ウエイトチャネル情報記憶回路１−６は、各端末装置２−ｉのウエイトチャネル情報を記憶する。
【００５４】
送信信号生成回路１−２は、トレーニングウエイト生成回路１−８から供給されたトレーニングウエイトを用いて、各端末装置２−ｉのアンテナとの間のチャネル情報を推定するための、トレーニングウエイトを用いたトレーニングシンボル（以下、ウエイトトレーニングシンボルという）を生成する。送信信号生成回路１−２は、ウエイトトレーニングシンボルに対し、ガードインターバルや、制御信号の付加を行い信号を生成する。無線信号送受信回路１−３は、送信信号生成回路１−２が生成した信号を搬送波周波数へアップコンバートし、送受信アンテナ１−４−１〜１−４−Ｎを介して送信を行う。
【００５５】
各端末装置２−iの無線信号送受信回路２−ｉ−２は、基地局装置１が送信した信号を送受信アンテナ２−ｉ−１−Ｍｉを介して受信する。受信信号復調回路２−ｉ−３は、無線信号送受信回路２−ｉ−２が受信した信号に対して同期、復調を行う。チャネル推定回路２−ｉ−４は、復調された信号に対してチャネル推定を行い、トレーニングウエイト込みの信号に対するチャネル情報（以下、ウエイトチャネル情報という）、を、送信信号生成回路２−ｉ−５へ出力する。送信信号生成回路２−ｉ−５は、ウエイトチャネル情報を含めた送信信号を生成し、予め定められたタイミングで、無線信号送受信回路２−ｉ−２へ出力する。無線信号送受信回路は、送信信号生成回路２−ｉ−５が出力したウエイトチャネル情報を含む信号を送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉの少なくとも１つを介して送信する。なお、チャネル推定回路２−ｉ−４は、送信信号生成回路２−ｉ−５に対して推定したウエイトチャネル情報をそのまま出力するのではなく、圧縮するなどの処理をして生成した情報をウエイトチャネル情報として出力するようにしてもよい。
【００５６】
基地局装置１の無線信号送受信回路１−３は、送受信アンテナ１−４−１〜１−４−Ｎのうち少なくとも一つを介して端末装置２−ｉから送信された信号を受信する。無線信号送受信回路１−３は、受信した信号をダウンコンバートしデジタルデータに変換し、受信信号復調回路１−５へ出力する。受信信号復調回路１−５は、受信信号を復調するとともに、受信信号に含まれるウエイトチャネル情報を抽出し、ウエイトチャネル情報記憶回路１−６に書き込んで記憶させる。
【００５７】
基地局装置１のデータ選択出力回路１−１は、通信相手となる端末装置２−ｉのウエイトチャネル情報がウエイトチャネル情報記憶回路１−６に記憶されると、選択した通信相手となる端末装置２−ｉの情報を送信方法選択回路１−７に出力し、送信するデータを送信信号生成回路１−２に出力する。送信方法選択回路１−７は、ウエイトチャネル情報記憶回路１−６に記憶されたウエイトチャネル情報と、トレーニングウエイト生成回路１−８の内部の記憶領域に記憶されているトレーニングウエイトから、送信ウエイトを算出し、空間ストリーム数、変調方式、符号化率からなる送信モードと算出した送信ウエイトを送信信号生成回路１−２に出力する。送信信号生成回路１−２は、通信を行う端末装置２−ｉに対し、送信信号に対して送信方法選択回路１−７から出力された送信モードに従い変調・符号化を行い、送信ウエイトを乗算し、信号検出や通信情報伝達に用いるパイロット信号を挿入し、無線信号送受信回路１−３へ出力する。無線信号送受信回路１−３は、送信信号生成回路１−２が出力した信号を搬送波周波数にアップコンバートし、送受信アンテナ１−４−１〜１−４−Ｎを介して送信する。
【００５８】
データ選択出力回路１−１において通信相手となる端末装置２−ｉを選択する手法としては、送信を行うデータがメモリに記憶されており、その中から送信を行う準備ができている端末装置２−ｉを選択する手法や、記憶されている複数のデータのうち最も古いデータに対応する端末装置２−ｉから選択する手法がある。また、端末装置２−ｉを利用するユーザに割り当てられているＱｏＳ（Quality of service）に基づいて選択したり、予めグループＩＤで決められたユーザに対応する端末装置２−ｉの組み合わせを選択したり、チャネル情報の相関が低い端末装置２−ｉの組み合わせを選択したりする手法もある。
【００５９】
次に、図２を参照しつつトレーニングウエイト生成回路１−８におけるトレーニングウエイトの生成処理の流れについて説明する。まず、最初に基地局装置１は、通信を行う対象となる各端末装置２−ｉを送信方法選択回路１−７により決定する。通信を行う対象は、データ選択出力回路から得られるデータを送る相手となる端末装置の情報、送信を行うデータの種類や優先度、端末装置に対する信号対雑音電力比の情報、などから決定することができる。基地局装置１は、各端末装置２−ｉのアンテナとの間の伝搬環境の情報を示すチャネル情報を収集する。通信相手となる端末装置が決定されると、送信信号生成回路１−２は、トレーニングウエイト生成回路１−８にトレーニングウエイトを要求する（ステップＳ１０１）。トレーニングウエイト生成回路１−８は、送信信号生成回路１−２からの要求を受けると、ウエイトチャネル情報記憶回路１−６に端末装置２−ｉのウエイトチャネル情報が記憶されていない場合、またはウエイトチャネル情報が、予め定めた時間Ｔｆよりも古い情報である場合、内部の記憶領域に記憶しているランダムウエイトＷＲ，ｎのみを用いてトレーニングウエイトとし、更新回数ｖに１（ｖ＝１）を設定する（ステップＳ１０２）。
【００６０】
ランダムウエイトＷＲ，ｎは、ｎ番目の周波数チャネルに対するランダムウエイトであり、初期値としては、全て同じＷＲとして設定しておくこともできる（ＷＲ，１＝ＷＲ，２＝・・・＝ＷＲ，Ｆ＝ＷＲ）。また、各周波数チャネルにおけるランダムウエイトは後述するように通信により更新されるため、全周波数チャネルに対し共通に用いることのできる共通ランダムウエイトＷＲを各周波数チャネルに記憶するランダムウエイトＷＲ，ｎとは別に記憶しておくこともできる。各周波数チャネルのランダムウエイトＷＲ，ｎとしては、通信相手となる当該端末装置の組み合わせが新たに決まった時点で、ＷＲに初期化してもよいし、または、前回の通信終了時に記憶していた各周波数チャネルに対するＷＲ，ｎをそのまま用いることもできる。ランダムウエイトＷＲ，ｎは、Ｎ×Ｎのユニタリ行列であり、対角要素１、非対角要素０とした行列、Ｎ×Ｎのウォルシュの直交符号行列を用いることができる。また、Ｎ×Ｎのフーリエ変換、逆フーリエ変換行列を用いてもよく、またランダム行列に対してＱＲ分解や特異値分解を適用することにより得られたユニタリ行列を用いるようにしてもよい。端末装置２−ｉのウエイトチャネル情報がない場合には、ｎ番目の周波数チャネルのトレーニングウエイトＷ’ｎは、式（１３）として表される。
【００６１】
【数１３】
2014045286.tif 000015
【００６２】
式（１３）において、関数［Ａ］ａ，ｂ は、行列Ａのa番目の列行列からｂ番目の列行列まで（ｂ−ａ）＋１の列ベクトルを選択する関数である。ＬＡは、トレーニング信号の数であり、式（１４）として表すことができる。
【００６３】
【数１４】
2014045286.tif 000016
【００６４】
式（１４）において、Ｌｉはi番目の端末装置２−ｉに対してチャネル情報を推定する空間多重数であり、後述するように、端末装置２−ｉが基地局装置１にフィードバックするチャネルの次元である。端末装置２−ｉに対するフィードバックの要求は、例えば、基地局装置１の送信方法選択回路１−７が、端末装置との通信環境（信号対雑音費、送信ビット数、伝送品質、スループット、データの種類、優先度、バッファに残るデータビット数）からＬｉを決定し、各端末装置２−ｉに対して送信することになる。Ｌｉは端末装置２−ｉの受信アンテナ数Ｍｉとすることで最大のチャネル情報量を得る。通信環境に応じ、ＬｉをＭｉより少ない数とすることもできる。また、端末装置側において、基地局装置１に指定されたＬｉ個のベクトルをフィードバックできず、（Ｌｉ−γ）個のベクトルをフィードバックした場合には、以降の通信においてＬｉ＝Ｌｉ−γとして、チャネル推定を行うこともできる。このようにして、端末装置２−ｉがＬｉを間接的に決定することもできる。Ｌ０は、付加トレーニングウエイトに含まれる列ベクトルの数である。Ｌｉは各端末装置に対してデータを送信する際の空間多重数と同数、またはそれより多い数に設定できる。式（１４）において、ＬＡを小さく設定するほど、ＭＡＣサブレイヤにおける伝送効率の低下を小さくすることができ、またフィードバック量を小さくすることができる一方、スループットの収束時間は長くなる。
トレーニングウエイトは、逐次的に更新されるため、ｖ回目の更新で得られるトレーニングウエイトをＷ’ｎ（ｖ）と定義する。したがって、式（１３）によって生成されたトレーニングウエイトはｖ＝１ととらえることができ、式（１５）のように表すことができる。
【００６５】
【数１５】
2014045286.tif 000017
【００６６】
このようにして、トレーニングウエイト生成回路１−８は、更新回数ｖ＝１においてランダムウエイトからＬＡ個を選択してトレーニングウエイトとし、送信信号生成回路１−２に出力する（ステップＳ１０３）。送信信号生成回路１−２は、トレーニングウエイト生成回路１−８が出力したトレーニングウエイトを用いてウエイトチャネル推定用のウエイトトレーニングシンボルを生成して各端末装置２−ｉに送信する（ステップＳ１０４）。
【００６７】
図３にトレーニングシンボルの構成例を示す。図３の四角いブロックＷＬＴ（１，１）〜ＷＬＴ（Ｎ，ＬＡ）は、チャネル推定のためのＯＦＤＭシンボルを表しており、ガードインターバル（ＧＩ）が付与されたものとなっている。ＷＬＴ（ｊ，ｋ）は、ウエイトトレーニングシンボルブロックのうち、ｋ番目のタイミングにおいてｊ番目のアンテナから送信されるＯＦＤＭシンボルに対応する。ＷＬＴ（ｊ，ｋ）の生成は、周波数チャネルにおける既知信号ｓ１〜ｓＦとトレーニングウエイトＷ’１ 〜 Ｗ’Ｆに基づく信号に逆フーリエ変換を行い、ＧＩを付与することによって行われる。ＷＬＴ（１，１）〜ＷＬＴ（Ｎ，ＬＡ）のｎ番目の周波数チャネルの送信信号ｘｊ，ｋ，ｎからなる送信信号行列Ｘｎは、式（１６）として表される。
【００６８】
【数１６】
2014045286.tif 000018
【００６９】
式（１６）において、ｘｊ，ｋ，ｎは、ＷＬＴ（ｊ，ｋ）のｎ番目の周波数チャネルに対応する送信信号を表す。行列Ｗ’ｎ（ｖ）は、ｎ番目の周波数チャネルに対応するトレーニングウエイトである。行列Ａは、トレーニングシンボル用の符号行列であり、送信側と受信側の両側で既知の行列が用いられる。行列Ａは、本実施形態では、ＬＡ×ＬＡの行列として定義される。行列Ａとして、単位行列Ｉを用いれば、それぞれ異なるタイミングでトレーニングウエイトにしたがってｓｎを送信することができる。また、行列Ａとして任意の行列やＡＨＡ＝Ｉとなるような直交行列を用いることもできる。
【００７０】
次に、端末装置２−ｉの無線信号送受信回路２−ｉ−２が、送受信アンテナ２−ｉ―１−１〜２−ｉ−１−Ｍｉを介して基地局装置１からの信号を受信すると、受信信号復調回路２−ｉ−３は、受信信号に対し同期を行う。チャネル推定回路２−ｉ−４が、ＷＬＴ（１，１）〜ＷＬＴ（Ｎ，ＬＡ）に対応する受信信号から、チャネル情報、すなわち基地局装置１と端末装置２−ｉのアンテナ間の伝搬係数からなるトレーニングウエイト込みのチャネル行列（以下、ウエイトチャネル行列という）を推定する。ｊ番目のトレーニングシンボルに対応する受信信号において、送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉで受信されたｎ番目の周波数チャネルに対応する受信信号ｙ１，ｊ，ｉ，ｎ〜ｙＭｉ，ｊ，ｉ，ｎ からなる受信信号ベクトルｙｊ，ｉ，ｎは、以下の式（１７）として表される。
【００７１】
【数１７】
2014045286.tif 000019
【００７２】
式（１７）において、行列Ｈｉ，ｎは、端末装置２−ｉのｎ番目の周波数チャネルに対するチャネル情報を表すチャネル行列（Ｍｉ×Ｎ行列）である。チャネル行列Ｈｉ，ｎのｐ列ｑ行目の要素は、基地局装置１のｐ番目の送信アンテナ（送受信アンテナ１−４−ｐ）から、端末装置２−ｉのｑ番目の受信アンテナ（送受信アンテナ２−ｉ―１−ｑ）間の伝搬係数を表す。ベクトルｎｊ，ｉ，ｎは、ｊ番目のトレーニングシンボルの受信タイミングで、送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉが受信した端末装置２−ｉのｎ番目の周波数チャネルの熱雑音ベクトルを表す。したがって、端末装置２−ｉにおける、第ｎの周波数チャネルにおける、１〜ＬＡ番目のトレーニングシンボルに対応する受信信号行列Ｙｉ，ｎは、式（１８）として表すことができる。
【００７３】
【数１８】
2014045286.tif 000020
【００７４】
式（１８）において、行列Ｎｉ，ｎは、端末装置２−ｉのｎ番目の周波数チャネルにおいて、送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉが１〜ＬＡ番目の受信タイミングで受信した信号に対する熱雑音行列を表す。端末装置２−ｉにおいて、行列Ａとｓｎは既知であるため、（Ａｓｎ）−１を受信信号行列Ｙｉ，ｎに乗算することで、ウエイトチャネル行列を式（１９）のように求めることができる。
【００７５】
【数１９】
2014045286.tif 000021
【００７６】
式（１９）において、行列Ｈ’’ｉ，ｎが、推定されたウエイトチャネル行列であり、Ｍｉ×ＬＡの行列である。端末装置２−ｉは、得られたＨ’’ｉ，ｎをそのままフィードバック情報とすることもできるし、Ｈ’’ｉ，ｎの中からＬｉ個の行ベクトルを選択してフィードバックすることもできる。また、Ｈ’’ｉ，ｎに特異値分解を行って式（２０）のように表し、信号空間に対応する右特異行列Ｖ（ｗｓ）ｉ，ｎ（ＬＡ×Ｍｉ行列）のうちＬｉ個の列ベクトルを、フィードバックすることもできる。
【００７７】
【数２０】
2014045286.tif 000022
【００７８】
式（２０）において、Ｕ（ｗ）ｉ，ｎは、左特異行列である。行列Σ（ｗ）ｉ，ｎは、特異値を対角要素とし、非対角項が０のＭｉ×ＬＡの対角行列である。行列Ｖ（ｗｎ）ｉ，ｎ（ＬＡ×（ＬＡ−Ｍｉ）行列）は、右特異行列のうち、特異値に対応しないベクトルの集まりを表す。上記の説明では、基地局装置１のアンテナ数ＮがＬＡより多いことを前提としていた。これに対して、端末装置２−ｉの送受信アンテナ数Ｍｉの方がＬＡより多い場合には、右特異行列の列ベクトルの全てが特異値と対応するため右特異行列Ｖ（ｗｓ）ｉ，ｎ（ＬＡ×ＬＡ行列）のうちＬｉ個の列ベクトルを選択してフィードバックすることができる。また、Ｈ’’ｉ，ｎのエルミート転置行列の列ベクトルにグラムシュミットの直交化法を用いて得られる基底ベクトルのうち、Ｌｉ個の列ベクトルを選択して、フィードバックすることができる。直交化法による基底ベクトルは、式（２１）のようにＱＲ分解によっても得ることができる。
【００７９】
【数２１】
2014045286.tif 000023
【００８０】
式（２１）により得られた直交行列（Ｑ（ｗｓ）ｉ，ｎ  Ｑ（ｗｎ）ｉ，ｎ）のうちＭｉ×Ｍｉの三角行列Ｒｉ，ｎに対応するＱ（ｗｓ）ｉ，ｎのうち、Ｌｉ個の列ベクトルをフィードバックするようにしてもよい。
【００８１】
上述したように、Ｌｉは、基地局装置が端末装置２−ｉにフィードバックを要求したチャネルの次元であり、実際に基地局装置が端末装置２−ｉに送信を行う空間ストリーム数Ｌ’ｉはＬｉ以下に設定される。Ｑ（ｗｓ）ｉ，ｎやＶ（ｗｓ）ｉ，ｎは、ユニタリ行列の一部となっている。チャネル情報から得られるこれらの行列、すなわちウエイトチャネル行列Ｈ’’ｉ，ｎ、またはこの行列から選択したＬｉ個のベクトルや、この行列に上述した特異値分解や直交化法を行うことにより得られた行列から選択したＬｉ個のベクトルが上述したウエイトチャネル情報であり、ウエイト信号空間行列Ｇ’ｉ，ｎ、と定義する。
【００８２】
端末装置２−ｉの送信信号生成回路２−ｉ−５は、ウエイト信号空間行列Ｇ’ｉ，ｎを再生可能な形に圧縮して送信信号を生成する。例えば、ユニタリ行列としての特徴を用いて角度情報に変換したり、周波数領域の信号を時間領域の信号に変換したり、空間相関を利用してＣｏｍｐｒｅｓｓｉｖｅ  ｓｅｎｓｉｎｇなどを用いることで情報を圧縮できる。そして、無線信号送受信回路２−ｉ−２が、生成された送信信号を基地局装置１に送受信アンテナ２−ｉ−１−１〜２−ｉ−１−Ｍｉを介して送信する。基地局装置１の無線信号送受信回路１−３は、送受信アンテナ１−４−１〜１−４−Ｎを介して端末装置２−１〜２−Ｋからの信号を受信し、受信信号復調回路１−５は、それぞれの信号からウエイトチャネル情報を読み出してウエイトチャネル情報記憶回路１−６に書き込んで記憶させる。
【００８３】
上記のようにしてウエイトチャネル情報を取得することで、基地局装置１は、チャネル推定用に用いるＷＬＴのタイムスロットの数をＬＡ個に削減でき、各端末装置２−ｉがフィードバックするウエイト信号空間行列Ｇ’ｉ，ｎのサイズも、ＬＡ×Ｌｉに削減できる。また、信号処理に用いるウエイト信号空間行列Ｇ’ｉ，ｎのサイズが小さくなっているため送信ウエイトの演算量も低減することができる。
これに対して、ウエイトチャネル行列としては基地局装置１の送受信アンテナ１−４−１〜１−４−Ｎが有するＮ次元の信号空間のうち、ＬＡ次元の部分空間しか推定していないため、Ｎ個のトレーニング信号をそのまま用いた方式よりも特性が劣化する。この特性劣化を防ぐために、本実施形態の構成ではトレーニングウエイトを逐次的に更新する。
【００８４】
基地局装置１が、端末装置２−ｉと一度以上通信を行うと、通信相手の信号空間に対するウエイト信号空間行列Ｇ’ｉ，ｎがウエイトチャネル情報記憶回路１−６に記憶される。記憶されているウエイト信号空間行列が更新されると、更新回数ｖをｖ＝ｖ＋１として、トレーニングウエイトの更新回数を増加させ、トレーニングウエイト生成回路１−８は、ウエイト信号空間行列Ｇ’ｉ，ｎをウエイトチャネル情報記憶回路１−６から読み出す（ステップＳ１０５）。更新回数が、２回目以降（ｖ≧２）の場合、新たにトレーニングウエイトを生成するにはウエイト信号空間行列Ｇ’ｉ，ｎを用いる。まず、Ｋ個の端末装置２−１〜２−Ｋの組み合わせに対し式（２２）のようにしてＱＲ分解を行う。
【００８５】
【数２２】
2014045286.tif 000024
【００８６】
式（２２）によって得られた直交行列（Ｑ（ｐｓ）ｎ  Ｑ（ｐｎ）ｉ，ｎ）のうち、ΣＬｉ×ΣＬｉの三角行列Ｒ（ｐｓ）ｎに対応するＬＡ×ΣＬｉ行列Ｑ（ｐｓ）ｉ，ｎとトレーニングウエイトＷ（ｖ）ｎを新たな（ｖ＋１）回目のトレーニングウエイトの生成に用いる。なお、ΣＬｉは、Ｌｉのｉ＝１〜Ｋまでの総和を表す。ｖ＝ｖ＋１とし、新たなトレーニングウエイトＷ’ｎ（ｖ）は、式（２３）として算出することができる（ｖ≧２）。
【００８７】
【数２３】
2014045286.tif 000025
【００８８】
式（２３）は、ΣＬｉ個のＮ×１の列ベクトルＷ（ｖ−１）ｎＱ（ｐｓ）ｎに、ランダムウエイトから新たな列ベクトルをＬ０個加えた構成となっている。ここで、Ｗ（ｖ−１）ｎＱ（ｐｓ）ｎを信号空間トレーニングウエイトと定義する。信号空間トレーニングウエイトはすでに用いられたランダムウエイトの部分空間となっているため、新たに加えたＬ０個の列ベクトルは、信号空間トレーニングウエイトＷ（ｖ−１）ｎＱ（ｐｓ）ｎと直交しており、内積を計算すると０になる。ｖ＝１の段階でのトレーニングウエイト込みのチャネル情報推定では、１〜ＬＡ個のランダムウエイトに対応するチャネル情報を算出するのに対して、毎回ＬＡ個のトレーニングシンボルを用いながら、ｖ＝２ではＬＡ＋Ｌ０、ｖ＝３ではＬＡ＋２Ｌ０と、付加トレーニングウエイトにより徐々に推定する信号空間の次元を増やすことができるのが本実施形態の処理の特徴である。
【００８９】
ここで、ランダムウエイトは、Ｎ個のベクトルであるため、Ｎと、ｖ、ＬＡ、Ｌｉの関係が式（２４）になると、新たにランダムウエイトからベクトルを選択することができない。
【００９０】
【数２４】
2014045286.tif 000026
【００９１】
そこで、トレーニングウエイト生成回路１−８は、ステップＳ１０５においてｖ＝ｖ＋１とした後、ｖが式（２４）の関係を満たすか否か、すなわち選択されていないＬ０個のベクトルがランダムウエイトに残っているか否かを判定する（ステップＳ１０６）。トレーニングウエイト生成回路１−８は、ｖが式（２４）の関係を満たさないと判定した場合、すなわち、選択されていないＬ０個のベクトルがランダムウエイトに残っている場合、式（２３）に従ってトレーニングウエイトを更新する。そして、ステップＳ１０４に戻り更新されたトレーニングウエイトを用いてチャネル推定を行い、選択されていないＬ０個のベクトルがランダムウエイトに残っている間、ステップＳ１０４からの処理を繰り返す（ステップＳ１０７）。一方、ｖが式（２４）の関係を満たす場合、すなわち、選択されていないＬ０個のベクトルがランダムウエイトに残っていない場合、トレーニングウエイト生成回路１−８は、ランダムウエイトを更新する（ステップＳ１０８）。
【００９２】
ランダムウエイトの更新は、ウエイト信号空間行列Ｇ’ｉ，ｎと前回のトレーニングウエイトＷ’ｎ（ｖ−１）を用いて式（２５）に示すようにＱＲ法により更新することができる。
【００９３】
【数２５】
2014045286.tif 000027
【００９４】
式（２５）において、左辺はＮ×ＬＡ行列のトレーニングウエイトとＬＡ×（Ｌ１＋Ｌ２＋・・・＋ＬＫ）行列となるＧ’ｉ，ｎの集合行列の乗算で得られるＮ×（Ｌ１＋Ｌ２＋・・・＋ＬＫ）＝Ｎ×ΣＬｉ行列である。直交行列Ｑ（ｒ）ｎは、Ｎ×Ｎの行列となり、信号空間（上三角行列Ｒ（ｒｓ）ｎ）に対応するＮ×ΣＬｉ行列のＱ（ｒｓ）ｎとＮ×（Ｎ−ΣＬｉ）行列のＱ（ｒｎ）ｎから形成される。ランダムウエイトＷＲ，ｎは新たに、式（２６）として求められる。
【００９５】
【数２６】
2014045286.tif 000028
【００９６】
トレーニングウエイト生成回路１−８は、式（２６）により求められたランダムウエイトを内部の記憶領域に書き込んで記憶させてステップＳ１０２に戻り、トレーニングウエイトの生成を繰り返す。これにより、ランダムウエイトは、初回は、チャネル情報となんら関係ない完全にランダムなものが用いられるのに対して、更新されることにより、ランダムウエイトＷＲ，ｎそのものがチャネル情報から算出されたものとなり、ランダムウエイトのみからトレーニングウエイトを選択してもチャネルに最適化されたものを利用することが可能となる。ステップＳ１０２の式（１５）は、Ｑ（ｒ）ｎを用いて表すと、式（２７）として表される。
【００９７】
【数２７】
2014045286.tif 000029
【００９８】
チャネルが時間に対し全く変動しなければ、Ｑ（ｒｓ）ｎは全端末装置２−１〜２−Ｋのチャネル行列の信号空間と一致し、Ｑ（ｒｓ）ｎはヌル信号空間に対応するため、Ｑ（ｒｓ）ｎのいかなる列ベクトルを送信ウエイトとしても、各端末装置２−ｉに信号は届かないことになる。これに対して、実際のシステムではチャネルは変動し、チャネル情報には熱雑音や量子化雑音などの誤差が含まれてしまう。このチャネル情報のずれを補正するため、本実施形態の基地局装置１は、逐次的にＱ（ｒｓ）ｎのトレーニングウエイトを加え、継続的にウエイトチャネル情報を効率よく推定し続けることができる。したがって、本実施形態の基地局装置１は、トレーニングウエイトの更新、及びランダムウエイトの更新を繰り返し行うことで、ＬＡ個のトレーニングシンボルだけを用いて、送信アンテナ数分のトレーニングシンボルを送信した場合と同等の信号空間に対するチャネル情報を取得することができる。
【００９９】
図２のフローチャートは、送信方法選択回路１−７において通信相手となる端末装置が決定された時に開始し、この通信相手の組み合わせが終了・変更された際に終了とすることができる。例えば、送信方法選択回路１−７の通信相手を示すグループＩＤの生成によりフローチャートを開始し、当該グループＩＤの割り当てを変更または取り消しした際に終了とすることができる。または、ステップＳ１０２からステップＳ１０８までのいずれかのステップにおいて、一定時間Ｔｆ経過しても同じステップまで戻ってこない場合に終了とすることができる。例えば、Ｓ１０４におけるウエイトトレーニングシンボルの送信が１秒間行われなかった場合に、フローチャートを終了することができる。
【０１００】
また、図２の処理の途中で、トレーニングウエイト生成回路１−８は、フィードバックを行うベクトル数Ｌｉ、付加トレーニングウエイトに含まれる列ベクトルの数Ｌ０、通信相手となる端末装置２−ｉの組み合わせを変更することができる。フィードバックを行うベクトル数Ｌｉの変更は、送信方法選択回路１−７が端末装置２−ｉとの通信環境に応じて行う。または、端末装置２−ｉからフィードバックされた信号空間行列のベクトル数を元にＬｉを更新することもできる。また、これらの変更は、基地局装置１の運用者の操作を受けて任意に変更させてもよいし、基地局装置１が内部に備える計時手段であるタイマが計測する時間の経過、あるいはネットワークを通じて受信する外部の時計やタイマなどの計時手段からの時間の情報に基づく時間の経過とともにトレーニングウエイト生成回路１−８が変化させてもよい（ステップＳ１０９）。これらが更新されることにより、ステップＳ１０３の式（１５）においてトレーニングウエイトをランダムウエイトから選択する数が変更され、ステップＳ１０４の式（２２）におけるウエイトチャネル行列Ｇ’ｉ，ｎのベクトル数が変更され、ステップＳ１０６において選択するＬ０の数が変更されることになる。
【０１０１】
ＬＡを一定値とした場合に、フィードバックを行うベクトル数Ｌｉが変更されると、ＬＡが一定となるように、式（１４）に従ってＬ０を変更することになる。このような条件では、Ｌ０が時間的に変化するため、新たに選択可能なランダムウエイトが存在するかどうかの判定式は式（２４）で表現することができなくなる。この場合には、ステップＳ１０３におけるランダムウエイトからトレーニングウエイトを選択する数を初期値として、ステップＳ１０７においてランダムウエイトから新たにベクトルを選択するごとに、選択したベクトル数を累積して記憶領域に記憶しておき、ステップＳ１０６においてまだ使われていないＬ０個のベクトルがランダムウエイトに残っているか判定する必要がある。この場合の判定条件は、ステップＳ１０３において初めにランダムウエイトから選択したベクトル数をＬ（１）、ステップＳ１０７において新たに追加するランダムウエイト数をＬ（ｖ）とした場合、式（２８）として表すことができる。
【０１０２】
【数２８】
2014045286.tif 000030
【０１０３】
式（２８）を満たしている場合には、トレーニングウエイト生成回路１−８は、ステップＳ１０６において、新たに選択できるランダムウエイトが残っていると判定することができ、ステップＳ１０７の処理に進むことができる。
フィードバックを行うベクトル数Ｌｉを変化させることにより、端末装置２−ｉからのフィードバック量を増加または軽減することができる。これに伴って基地局装置１の演算量も増加または軽減することになるため使用する環境に応じた負荷の調整を行うことが可能となる。また、Ｌ０を任意に変化させることで、例えば、上述したタイマが前回のチャネル推定から経過した時間を計測することにより、その計測した時間に従ってＬ０を増加させることで、後述するようにサイクルタイムを一定にすることができ、それにより、伝送特性を維持することが可能となる。
また、ステップＳ１０６において、Ｌ０個のベクトルが残っているか判定する代わりに、１つ以上のベクトルが残っているか判定することもできる。この場合、一つでもランダムウエイトの選択されていないベクトルが残っていれば、このランダムウエイトを付加トレーニングウエイトとして用いる。このときは、付加トレーニングウエイトの数を一時的にＮ−ΣＬ（ｉ）とするか、Ｌ０に足りない分は既に選択されているランダムウエイトから重複して選択することができる。このようにして、Ｎ個のランダムウエイトのベクトルを使い切り、全ての信号空間に対応する情報を推定することで、チャネル情報の推定精度を高めることもできる。
【０１０４】
上記のようにして更新を繰り返したトレーニングウエイトを用いて、基地局装置１が、端末装置２−ｉに対して送信を行う処理について説明する。送信方法選択回路１−７は、端末装置２−ｉからフィードバックされたウエイト信号空間行列Ｇ’ｉ，ｎ、をウエイトチャネル情報記憶回路１−６から読み出す。ウエイト信号空間行列Ｇ’ｉ，ｎ、すなわち端末装置２−ｉにおける第ｎ番目の周波数チャネルに対して得られるＬｉ個のベクトルからなるＬＡ×Ｌｉの行列は、推定したウエイトチャネル行列のうちのＬｉ個の行ベクトルの複素共役転置行列であってもよい。また、受信ウエイト行列βｉ，ｎを端末装置２−ｉにおいて予め定義し式（２９）のようにウエイト信号空間行列Ｇ’ｉ，ｎを表してもよい。
【０１０５】
【数２９】
2014045286.tif 000031
【０１０６】
また、前述したように、ウエイトチャネル行列Ｈ’’ｉ，ｎに基づいて特異値分解やＱＲ分解することにより得られるＶ（ｗｓ）ｉ，ｎやＱ（ｗｓ）ｉ，ｎを用いて以下の式（３０）、（３１）としても表すことができる。
【０１０７】
【数３０】
2014045286.tif 000032
【０１０８】
【数３１】
2014045286.tif 000033
【０１０９】
式（３０）、（３１）において、［Ａ］Ｌは、Ｎ×Ｍ行列ＡからＬ個の列ベクトルを選択してＮ×Ｌの行列を得る関数である。端末装置２−ｉにおける受信ウエイト行列βｉ，ｎとしては、Ｈ’ｉ，ｎをＱＲ分解して得られるＭｉ×Ｍｉのユニタリ行列のうち、対応する三角行列の対角成分が大きいものからＬｉ個のベクトルをＭｉ×Ｍｉのユニタリ行列の中から選択して得られるＭｉ×Ｌｉの行列にエルミート共役をとることで得られるＬｉ×Ｍｉの行列を用いたり、Ｈ’ｉ，ｎに対し特異値分解で得られる左特異行列のうち、大きい特異値に対応するＬｉ個のベクトルを選択し得られるＭｉ×Ｌｉ行列に対し、エルミート共役をとることで得られるＬｉ×Ｍｉの行列などを用いることができる。
【０１１０】
送信方法選択回路１−７は、これらのウエイト信号空間行列Ｇ’ｉ，ｎを従来における信号空間行列Ｇｉ，ｎとし、例えば、上述したＢＤ法による送信ウエイトを算出する手法などによりマルチユーザＭＩＭＯやシングルユーザＭＩＭＯなどの送信ウエイト演算手法を用いて送信ウエイトを算出することができる。ただし、算出した送信ウエイトは、トレーニングウエイトを含んだチャネル情報に基づく送信ウエイトとなる。そのため、ｖ回目の更新がされたトレーニングウエイトを含んだチャネル情報に基づいて算出されたｉ番目の端末装置のｎ番目の周波数チャネルにおける送信ウエイトＷ（ｖ）i，,nは、各送信アンテナ（送受信アンテナ１−４−１〜１−４−Ｎ）に対する送信ウエイトとはならない。各送信アンテナに対する送信ウエイトＷｉ，ｎは、ｖ回目で用いたトレーニングウエイトＷｎ’（ｖ）と送信ウエイトＷ（ｖ）i，,nにより式（３２）として算出されることになる。
【０１１１】
【数３２】
2014045286.tif 000034
【０１１２】
以下、図４、５を参照しつつ、上記の実施形態の構成による効果を説明する。チャネルモデルとして、ＩＥＥＥ８０２．１１ｎで提案されているＭｏｄｅｌ  Ｃのシナリオを用いた。アンテナ間隔０．５λのアレーアンテナ、ベルシェイプのドップラー分布を用い、ドップラースプレッドを０．４１Ｈｚ、コヒーレンスタイム８００ｍｓとしている。また、端末装置数を２とし、基地局装置１と端末装置２−１、２−２の間の距離を５ｍ、帯域幅２０ＭＨｚ、フーリエ変換による周波数チャネル数を６４、そのうちデータを送信する周波数チャネル数を５２（Ｆ＝５２）としている。また、受信アンテナ数２、フィードバックを行うベクトル数を２（Ｌｉ＝２）、データ送信時の空間多重数を２（Ｌｉ’＝２）としている。図４は、この条件において、通信を開始してからの時間（ｍｓ）に対する伝送可能ビット数（Ａｃｈｉｅｖａｂｌｅ  ｂｉｔ  ｒａｔｅ）について測定した結果である。伝送可能ビット数は、式（３３）として定義している。
【０１１３】
【数３３】
2014045286.tif 000035
【０１１４】
式（３３）において、σ２は熱雑音の分散値である。また、本実施形態では、ウエイトトレーニングシンボルを送信する時間間隔をＴ０と定義した。図４では、Ｔ０＝１０ [ｍｓ]としている。また、送信ウエイトＷｋ，ｎは、１ｍｓ前のチャネル情報を元に算出されていると仮定し、データ通信中のチャネルの時変動は考慮しないものとしている。
【０１１５】
図４において、点線、破線、及び実線で示されているラインは、それぞれ送信アンテナ数が８、１６、３２の場合に従来のチャネル推定手法により、送信アンテナ数だけのトレーニングシンボルを送信し、チャネル情報の推定を行った場合の伝送可能ビット数を示している。一方、＋および×で示されるポイントが、それぞれ送信アンテナ数が１６および３２の場合に本実施形態の構成によるウエイトトレーニングシンボルを用いて算出されたチャネル情報を用いて通信を行った場合の、伝送可能な伝送可能ビット数である。１０ｍｓごとに、徐々に特性が改善していき送信アンテナ数が１６素子の場合には６０ｍｓ、３２素子の場合には、１４０ｍｓ経過した後には、ほぼ従来のチャネル推定手法、すなわち送信アンテナ数と同数のトレーニングシンボルをチャネル推定に用いた場合と同等の特性が得られていることがわかる。この収束に必要な時間は、式（３４）で示されるサイクルタイムＴＣで表すことができる。
【０１１６】
【数３４】
2014045286.tif 000036
【０１１７】
従来のチャネル推定手法を用いた場合、３２または１６のトレーニングシンボルを用いているのに対し、本実施形態の手法では、Ｌ１＋Ｌ２＋Ｌ０＝６のウエイトトレーニングシンボルしか用いていないため、チャネル推定のためのオーバヘッドを大きく削減できる。
また、６個のウエイトトレーニングシンボルのみを用いるだけでも、従来のチャネル推定手法を用いた８×２のＭＩＭＯ方式の伝送特性（図４の点線）を大きく上回る伝送特性が得られているもわかる。これは、例えば、トレーニングシンボルを送信できる数を標準化規格などで変更することなく、送信アンテナ数を増やすことができる効果があることを意味している。
【０１１８】
図４の結果は収束するものの、従来の送信アンテナ数だけトレーニングシンボルを送信するチャネル推定手法に対して特性が少し劣化する。これは、チャネルの時間変動のため、完全には信号空間に対応するトレーニングウエイトにならないためである。この特性劣化は、サイクルタイムが大きくなるに従って増加する。図５は、サイクルタイムを７０ｍｓ、１４０ｍｓ、２８０ｍｓ、５６０ｍｓとした場合の伝送可能ビット数の累積分布関数（ＣＤＦ：cumulative distribution function）に基づく累積確率分布を示した図である。×、□、○で示されているラインは、それぞれ送信アンテナ数が８、１６、３２の場合の従来のチャネル推定手法による分布を示したものである。従来の３２素子のアンテナ数だけチャネルを全て推定する方式（図５の○で示されるポイント）に比べ、収束に必要な時間であるサイクルタイムが増加するにつれ、伝送特性が徐々に劣化しているのが分かる。それでもなお、１６×２のＭＩＭＯチャネル方式における１６本のアンテナ数のチャネルを全て推定する方式（図５の□で示されるポイント）に比べ、パイロットシンボルを６つ（ＬＡ＝６）しか用いていないにもかかわらず、高い伝送特性を有していることがわかる。
【０１１９】
上記のように本実施形態の構成では、全ての自由度（送信アンテナ数）のチャネル情報を推定し終えるまでにかかるサイクルタイムに比例して特性が若干劣化する。したがって、式（３４）に基づいて付加トレーニングウエイトに含まれる列ベクトルの数Ｌ０を設定することで、サイクルタイムＴＣをあらかじめ定めた値以下と伝送特性が劣化しすぎないようにすることができる。図４と図５のシミュレーションでは、チャネル推定を一定間隔Ｔ０で行うことを仮定しているが、実際の通信において、このように定期的なチャネル推定を行えるとは限らない。そこで、前回チャネル推定を行ってから、経過した時間に応じてＬ０を設定する。例えば、図５の結果から、ＴＣを１４０ｍｓ以下になるように設定することを考える。この場合、１０ｍｓ間隔に、Ｌ０＝２でチャネル推定が行われていれば、サイクルタイム１４０ｍｓで通信できる。チャネル推定が、２０ｍｓ間隔になってしまう場合には、Ｌ０＝４とすることで、サイクルタイムを一定にできる。すなわち式（３５）に基づいてＬ０を決めることができる。
【０１２０】
【数３５】
2014045286.tif 000037
【０１２１】
式（３５）において、Ｔは前回のチャネル推定からの経過時間であり、式（３５）において、以下の式（３６）の演算子はＡ以上で、最も小さい整数を表す関数である。
【０１２２】
【数３６】
2014045286.tif 000038
【０１２３】
なお、本実施形態の構成は、マルチユーザＭＩＭＯ通信のため選択された複数の端末装置組み合わせに対してそれぞれ独立に用いることができる。すなわち、送信方法選択回路が決定した端末装置２−ｉの組み合わせごとにランダムウエイトＷＲ，ｎをトレーニングウエイト生成回路が記憶する。トレーニングウエイト生成回路１−８が、トレーニングウエイトの生成を行う場合には、送信方法選択回路１−７が指定する通信相手の端末装置の組み合わせに対応するランダムウエイトＷＲ，ｎを用いトレーニングウエイトの生成を行う。これにより、１つの基地局装置１を用いて複数の端末装置２−ｉの組み合わせを扱うことが可能である。
【０１２４】
また、マルチユーザＭＩＭＯ通信において、予め選択された端末装置の組み合わせの中で、一部の端末装置からのみ、チャネル情報のフィードバックを要求する場合や、パケット衝突などの問題により、一部の端末装置からのチャネル情報のフィードバックが取得できない場合もある。このような場合、式（２２）のＧ’ｉ，ｎは、全ての端末装置から取得することができないため、ｚ番目の端末装置からのチャネル情報が取得できなかった場合には、Ｇ’ｚ，ｎを用いずに式（２２）を計算するか、または前回の通信でｚ番目の端末装置からフィードバックされた過去の情報であるＧ’ｚ，ｎを代わりに用いることもできる。ただし、過去のＧ’ｚ，ｎを用いる場合には、対応するトレーニングウエイトが異なるため、信号空間トレーニングウエイトの演算に工夫が必要となる。ここで端末装置２−１〜２−３との通信を考え、チャネル情報のフィードバックを求めたにもかかわらず、端末装置２−３から正常に信号空間行列が取得できなかった場合を考える（ｚ＝３）。端末装置２−１と２−２からは信号空間行列が得られているため、これらの端末装置に対しては式（２２）と同様に、式（３７）として、信号空間行列が得られた端末に対するユニタリ行列をＱ’（ｐｓ）ｎとして得ることができる。
【０１２５】
【数３７】
2014045286.tif 000039
【０１２６】
端末装置２−３からは前回信号空間行列を取得した際のトレーニングウエイトを用いる。すなわち、相手からのチャネル情報についてフィードバックがない場合に対応するため、トレーニングウエイト生成回路１−８が過去のトレーニングウエイトを記憶しておくことができる。端末装置２−３に対応する信号空間トレーニングウエイトは、信号空間行列取得時の更新回数をｑとし、Ｗ（ｑ）ｎＧ’３，ｎと表わすことができる。よって、トレーニングウエイトは、式（３８）として得られる。
【０１２７】
【数３８】
2014045286.tif 000040
【０１２８】
または、予め、信号空間トレーニングウエイトのベクトルが互いに直交するように変換しておくこともできる。ＱＲ分解を用いて、式（３９）として得られたＱ（ｐｐｓ）ｎを用いて、式（４０）とすることができる。
【０１２９】
【数３９】
2014045286.tif 000041
【０１３０】
【数４０】
2014045286.tif 000042
【０１３１】
このようにすることで、ある特定の端末装置のウエイトチャネル情報が取得できなくてもウエイトチャネル情報の更新を行うことができる。
【０１３２】
上記の実施形態の構成により、基地局装置１は、トレーニング信号の数を低減しつつ、従来の方式と同等の伝送特性のＭＵ−ＭＩＭＯ通信を行うことが可能となる。
また、上記の実施形態の構成により、トレーニングウエイトを用いたトレーニングシンボルによるチャネル推定を行うことで、システムに対するチャネル推定フレームのオーバヘッド、チャネル関連情報のフィードバック量、チャネル関連情報をフィードバックするまでの時間差、を減らすことができ、高いシステムスループットを有する無線通信システムを実現することができる。
また、上記の実施形態の構成により、基地局装置１は、チャネル推定用に用いるウエイトトレーニングシンボルのタイムスロットの数をＬＡ個に削減することができる。これにより、各端末装置２−ｉにおけるウエイトチャネル行列の推定の演算量を軽減させることができ、フィードバックするウエイト信号空間行列Ｇ’ｉ，ｎのサイズも、ＬＡ×Ｌｉに削減することができる。これにより、チャネル推定におけるトラフィック量を減らすことができるとともに、基地局装置１における送信ウエイトの演算量も低減することができる。
また、上記の実施形態の構成により、基地局装置１は、付加トレーニングウエイトを加えてトレーニングウエイトを更新し、また、付加トレーニングウエイトの選択元のランダムウエイトが不足した場合にはランダムウエイトを更新してトレーニングウエイトを逐次的に更新することが可能となる。これにより、基地局装置１の送受信アンテナ１−４−１〜１−４−Ｎが有するＮ次元の信号空間のうちＬＡ次元の部分空間しか推定していないために生じる特性劣化を防き、Ｎ個のトレーニング信号をそのまま用いた方式と同程度の伝送特性を維持することが可能となる。
なお、本発明に記載の送信装置とは、例えば、本実施形態の基地局装置１である。また、本発明に記載の所定の数のベクトルとは、例えば、本実施形態のＬｉ個のベクトルである。
【０１３３】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【０１３４】
１  基地局装置
１−１  データ選択出力回路
１−２  送信信号生成回路
１−３  無線信号送受信回路
１−４−１〜１−４−Ｎ  送受信アンテナ
１−５  受信信号復調回路
１−６  ウエイトチャネル情報記憶回路
１−７  送信方法選択回路
１−８  トレーニングウエイト生成回路
２−１〜２−ｉ  端末装置
２−１−１−１〜２−ｉ−１−Ｍｉ  送受信アンテナ
２−１−２〜２−ｉ−２  無線信号送受信回路
２−１−３〜２−ｉ−３  受信信号復調回路
２−１−４〜２−ｉ−４  チャネル推定回路
２−１−５〜２−ｉ−５  送信信号生成回路
【図１】
2014045286.tif 000043
【図２】
2014045286.tif 000044
【図３】
2014045286.tif 000045
【図４】
2014045286.tif 000046
【図５】
2014045286.tif 000047
【図６】
2014045286.tif 000048
【図７】
2014045286.tif 000049